cancers Article ERK5 Is a Major Determinant of Chemical Sarcomagenesis: Implications in Human Pathology Elena Arconada-Luque 1,† , Jaime Jiménez-Suarez 1,† , Raquel Pascual-Serra 1 , Syong Hyun Nam-Cha 2 , Teresa Moline 3 , Francisco J. Cimas 4 , Germán Fliquete 3 , Marta Ortega-Muelas 1 , Olga Roche 1,5 , Diego M. Fernández-Aroca 1 , Raúl Muñoz Velasco 6,7 , Natalia García-Flores 1 , Cristina Garnés-García 1 , Adrián Sánchez-Fdez 8,9,10 , Sofía Matilla-Almazán 8,9,10 , Víctor J. Sánchez-Arévalo Lobo 6,7 , Javier Hernández-Losa 3 , Borja Belandia 11 , Atanasio Pandiella 8,9,10 , Azucena Esparís-Ogando 8,9,10 , Santiago Ramón y Cajal 3 , Luis del Peso 12,13,14 , Ricardo Sánchez-Prieto 1,5,11,15, *,‡ and María José Ruiz-Hidalgo 1,16,‡ 1 Laboratorio de Oncología Molecular, Unidad de Medicina Molecular, Centro Regional de Investigaciones Biomédicas, Unidad Asociada de Biomedicina UCLM, Unidad Asociada al CSIC, Universidad de Castilla-La Mancha, 02008 Albacete, Spain; elena.arconada@alu.uclm.es (E.A.-L.); jaime.jimenez@uclm.es (J.J.-S.); rpascuals@sescam.jccm.es (R.P.-S.); marta.ortega.muelas@gmail.com (M.O.-M.); oroche@gmail.com (O.R.); diego.fernandez@uclm.es (D.M.F.-A.); natalia.garcia16@alu.uclm.es (N.G.-F.); cristina.garnes@alu.uclm.es (C.G.-G.); maria.rhidalgo@uclm.es (M.J.R.-H.) 2 Servicio de Anatomía Patológica, Hospital General de Albacete, 02008 Albacete, Spain; shnam@sescam.jccm.es 3 Grupo de Patología Molecular Traslacional, Vall d’Hebron Institut de Recerca, Universitat Autònoma de Barcelona Centro de Investigación Biomédica en RED de Cancer CIBERONC, 08035 Barcelona, Spain; teresa.moline@vhir.org (T.M.); gfliquete@vhebron.net (G.F.); jahernan@vhebron.net (J.H.-L.); sramon@vhebron.net (S.R.y.C.) 4 Unidad de Bioquímica y Biología Molecular, Servicio de Instrumentación Biomédica, Citation: Arconada-Luque, E.; Universidad de Castilla-La Mancha, 02008 Albacete, Spain; franciscojose.cimas@uclm.es Jiménez-Suarez, J.; Pascual-Serra, R.; 5 Departamento de Ciencias Médicas, Facultad de Medicina, Universidad de Castilla-La Mancha, Nam-Cha, S.H.; Moline, T.; Cimas, 02008 Albacete, Spain F.J.; Fliquete, G.; Ortega-Muelas, M.; 6 Grupo de Oncología Molecular, Facultad de Ciencias Experimentales, Instituto de Investigación Biosanitaria, Roche, O.; Fernández-Aroca, D.M.; Universidad Francisco de Vitoria, Pozuelo de Alarcón, 28223 Madrid, Spain; raul.munoz@ufv.es (R.M.V.); et al. ERK5 Is a Major Determinant of victor.sanchezarevalo@ufv.es (V.J.S.-A.L.) 7 Departamento de Anatomía Patológica, Instituto de Investigación Hospital 12 de Octubre, Av. Córdoba, s/n, Chemical Sarcomagenesis: Implications in Human Pathology. 28041 Madrid, Spain 8 Instituto de Biología Molecular y Celular del Cáncer-CSIC, 37007 Salamanca, Spain; adriansf@usal.es (A.S.-F.); Cancers 2022, 14, 3509. https:// sofiamatilla@gmail.com (S.M.-A.); atanasio@usal.es (A.P.); esparis@usal.es (A.E.-O.) doi.org/10.3390/cancers14143509 9 Instituto de Investigación Biomédica de Salamanca (IBSAL), Hospital Universitario de Salamanca, Academic Editor: Domenico Universidad de Salamanca, CSIC, 37007 Salamanca, Spain 10 Centro de Investigación Biomédica en RED de Cancer CIBERONC, 37007 Salamanca, Spain Andrea Campanacci 11 Departamento de Biología del Cáncer, Instituto de Investigaciones Biomédicas ‘Alberto Sols’ (CSIC-UAM), Received: 20 June 2022 Unidad Asociada de Biomedicina UCLM, Unidad Asociada al CSIC, 28029 Madrid, Spain; Accepted: 16 July 2022 bbelandia@iib.uam.es 12 Departamento de Bioquímica, Universidad Autónoma de Madrid (UAM) and Instituto de Investigaciones Published: 19 July 2022 Biomédicas ‘Alberto Sols’ (CSIC-UAM), 28029 Madrid, Spain; lpeso@uam.es Publisher’s Note: MDPI stays neutral 13 Unidad Asociada de Biomedicina UCLM, Unidad Asociada al CSIC, 28029 Madrid, Spain with regard to jurisdictional claims in 14 Centro de Investigación Biomédica en Red de Enfermedades Respiratorias CIBERES, 28029 Madrid, Spain published maps and institutional affil- 15 Instituto de Investigaciones Biomédicas ‘Alberto Sols’, Consejo Superior de Investigaciones iations. Científicas (IIBM-CSIC)-Universidad de Castilla-La Mancha, 02008 Albacete, Spain 16 Departamento de Química Inorgánica, Orgánica y Bioquímica, Área de Bioquímica y Biología Molecular, Facultad de Medicina, Universidad de Castilla-La Mancha, 02008 Albacete, Spain * Correspondence: rsprieto@iib.uam.es † These authors contributed equally to this work. Copyright: © 2022 by the authors. ‡ Both senior authors have contributed by equal to this work. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and Simple Summary: Sarcoma is a heterogeneous group of tumors poorly studied with few therapeutic conditions of the Creative Commons opportunities. Interestingly, the role of MAPKs still remains unclear in sarcomatous pathology. Here, Attribution (CC BY) license (https:// we describe for the first time the critical role of ERK5 in the biology of soft tissue sarcoma by using creativecommons.org/licenses/by/ in vitro and in vivo approaches in a murine experimental model of chemical sarcomagenesis. Indeed, 4.0/). Cancers 2022, 14, 3509. https://doi.org/10.3390/cancers14143509 https://www.mdpi.com/journal/cancers Cancers 2022, 14, 3509 2 of 19 our observations were extrapolated to a short series of human leiomyosarcoma and rhabdomyosarco- mas. Furthermore, transcriptome analysis allows us to demonstrate the critical role of KLF2 in the biological effects of ERK5. Therefore, the data presented here open new windows in the diagnosis and therapy of soft tissue sarcomas. Abstract: Sarcomas are a heterogeneous group of tumors in which the role of ERK5 is poorly studied. To clarify the role of this MAPK in sarcomatous pathology, we used a murine 3-methyl-cholanthrene (3MC)-induced sarcoma model. Our data show that 3MC induces pleomorphic sarcomas with muscle differentiation, showing an increased expression of ERK5. Indeed, this upregulation was also observed in human sarcomas of muscular origin, such as leiomyosarcoma or rhabdomyosarcoma. Moreover, in cell lines derived from these 3MC-induced tumors, abrogation of Mapk7 expression by using specific shRNAs decreased in vitro growth and colony-forming capacity and led to a marked loss of tumor growth in vivo. In fact, transcriptomic profiling in ERK5 abrogated cell lines by RNAseq showed a deregulated gene expression pattern for key biological processes such as angiogenesis, migration, motility, etc., correlating with a better prognostic in human pathology. Finally, among the various differentially expressed genes, Klf2 is a key mediator of the biological effects of ERK5 as indicated by its specific interference, demonstrating that the ERK5–KLF2 axis is an important determinant of sarcoma biology that should be further studied in human pathology. Keywords: ERK5; MAPK7; soft tissue sarcoma; leiomyosarcoma; rhabdomyosarcoma; KLF2; 3-methyl-cholanthrene 1. Introduction Sarcomas are a group of heterogeneous tumors that develop from the connective tissue, which provides a supporting matrix throughout the organism. More than 150 types and subtypes of sarcomas have been described. The classification of sarcomas divides them into soft tissue sarcomas (STS) and bone sarcomas. In addition, there is a third group, gastrointestinal stromal tumors, which are soft tissue sarcomas considered a separate group due to their diagnostic and therapeutic characteristics (for a review, see [1]). Sarcomas are rare tumors that, in adults, comprise 1% of all cancers; however, in children, they represent about 15% of cancer cases. In 2022, it is estimated that more than 13,000 new cases of STS will be diagnosed in the USA, and more than 5000 people will die of this disease (ASCO, https://www.cancer.net/, accessed on 1 April 2022, ACS www.cancer.org, accessed on 1 April 2022). Therapeutically, the most common treatments are surgery and chemo/radiotherapy. However, the use of novel anti-tumor agents such as kinase inhibitors or immunotherapy is now being considered (for a review, see [2]). Currently, there are different animal models to study STS [3–5], which have revealed the implication of key proteins in sarcomagenesis, as in the case of p53 [6] or the AKT sig- naling pathway [7,8]. Indeed, most of this experimental observation has been extrapolated to human pathologies, which will certainly have clinical implications in the coming years. A method that has been widely used for the induction of sarcomas since the first half of the last century is the use of carcinogens, being 3-methyl-cholanthrene (3MC) a great exam- ple [9]. The tumor developed after intramuscular injection of this carcinogen is qualified as “sarcoma”, and it has been used as a basis for demonstrating the oncogenic/tumor sup- pressor properties of a multitude of proteins, including cohesin SA1, p53, ATM, and PIM kinases, among others [10–12]. Furthermore, the model of 3MC-derived sarcoma has been characterized by WES, showing specific alterations in oncogenes and tumor suppressors such as K-ras or Tp53 [13], also implicated in human tumors as leiomyosarcoma or rhab- domyosarcoma [14]. However, it is notorious that a definitive pathological characterization has not yet been proposed for this experimental model. Interestingly, MAPK signaling has been related to several types of sarcomas [15,16], and recent evidence supports that modulation of MAPKs could have therapeutic implica- Cancers 2022, 14, 3509 3 of 19 tions [17–19]. However, it is remarkable that there is no work detailing the role of MAPKs in the sarcoma model of 3MC. In fact, ERK5, one of the canonical members of MAPKs superfamily described more than 25 years ago [20,21], is highly expressed in sarcomatous pathology, as indicated in The Cancer Genome Atlas (TCGA) series (see below). Inter- estingly, this particular MAPK signaling pathway has been linked to cancer [22,23], and it is expected to exploit its clinical implications in the coming years [24]. Indeed, ERK5 has been related to several types of tumors such as breast [25], multiple myeloma [26], or cholangiocarcinoma [27]. However, no deep studies evaluating the role of the ERK5 signaling pathway in STS have been performed to date. Here, we describe how the ERK5 pathway, with the direct implication of KLF2, is necessary for chemical sarcomagenesis triggered by 3MC, suggesting that this signaling axis needs to be deeply studied in order to fully exploit its potential in diagnostic and therapy in human sarcomatous pathology. 2. Materials and Methods 2.1. Cell Lines HEK293T and C3H10t1/2 cells were purchased from ATCC (LGC, Barcelona, Spain). Cells were maintained in 5% CO2 at 37 ◦ C and grown in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% fetal bovine serum (FBS), 1% glutamine plus antibiotics. All cell culture reagents were provided by Lonza (Culteck, Madrid, Spain). 2.2. Antibodies and Chemicals The antibodies used are listed in Supplementary Table S1. 3MC, etoposide, and crystal violet were purchased from Merck (Madrid, Spain). XMD8-92 was obtained from Selleckchem (Madrid, Spain). 2.3. Animal Studies All the animal experimentation was carried out according to the NIH-Intramural Animal Guide for the Care and Use of Laboratory Animals and approved by the Ethics in Animal Care Committee of the University of Castilla-La Mancha (reference PR-2019-07-18). 2.4. Induction of Murine Sarcoma and Cell Lines Derivation C57BL/6 females (n = 4) were injected in the right hind with 3MC as previously described [28]. When the tumors reached a size of 1–1.5 cm in diameter, mice were euthanized by cervical dislocation. Cell lines were derived from the tumors as previously described [28]. In this study, two cell lines out of four (referred to as 3MC-C1 and 3MC- C3) were chosen as a representative model of murine sarcoma. Cells were maintained in 5% CO2 at 37 ◦ C and grown in DMEM supplemented with 10% FBS, 1% glutamine plus antibiotics. 2.5. Xenograft Assays For xenograft assays, 5 × 105 cells from 3MC-C1, 3MC-C3, or derived cell lines were sub- cutaneously injected into the back of 5/6-week-old female mice of J:NMRIFoxn1nu/Foxn1nu strain (Janvier). Tumors were measured by caliper every two days, and tumor volume was calculated according to the formula V = (D × d2 )/2 (where D is tumor length and d tumor width). Figures show an experiment out of 2 performed with 2 different pools of infection that rendered nearly identical results. 2.6. RNA Interference Assays, Lentiviral Production, and Infections Lentiviral production and cell infection were performed as previously described [28] by using the following shRNAs from Merck (Madrid, Spain): PLKO.1-shRNA ERK5-1 (TRC-N-0000232396), PLKO.1-shRNA ERK5-2 (TRC-N-0000232397), PLKO.1-shRNA KLF2 (TRC-N-0000295770) or PLKO.1-empty vector (SHC001). Forty-eight hours post-infection, cells were selected with puromycin (2.5 µg/mL Merck) for 72 h. Each experiment was Cancers 2022, 14, 3509 4 of 19 performed with at least two different pools of infection. Infected cells were discarded 15–20 days after selection, and new pools were generated. 2.7. Western Blotting Protein quantification and western blot were performed as previously described [29]. For tumor samples, 400 µL of lysis buffer was added to 0.2 g of tumor, and it was disaggre- gated using a Polytron PT-2100 (Kinematica AG, Malters, Switzerland). Once the tissue was homogenized, the above-mentioned protocol was followed. 2.8. Clonogenic Assay Clonogenic assays were performed as previously described [28] by using 200 cells/well seeded in 6-well plates and maintained for 12–14 days. The colonies with less than 5 mm diameter were discarded. 2.9. Growth Curves Growth curves were performed as previously described [30]. Briefly, 3 × 105 cells for 3MC-C1 and 5 × 105 cells for 3MC-C3 were seeded into 100 mm plates and counted on days 3, 6, and 9 by using an automated cell counter (Bio-Rad) and replated in the same manner up to day 9. This experiment was performed with 3 different pools of infection for each cell line. 2.10. Adhesion Assay For adhesion assays, 24-well plates were coated with fibronectin (10 µg/mL in PBS, 1 h, 37 ◦ C, Merck). Then, fibronectin was removed, and plates were blocked (0.5% BSA in DMEM, 45 min) and washed (0.1% BSA in DMEM). After cold shock, 1 × 105 cells were seeded. The medium was removed at different time points, and adhered cells were washed 3× with PBS. Cells were stained with crystal violet for 15 min, washed with deionized water, and dried at RT for 24 h. Then, crystal violet was dissolved in acetic acid (10% in deionized water, 15 min on mild shaking), and absorbance was measured at 590 nm using a Biokinetics plate reader. 2.11. RNA Isolation, Reverse Transcription, and Real-Time Quantitative PCR (RT-qPCR) Total RNA from cell and mice tumor samples (after tissue homogenization with a Polytron) was obtained as previously described [28]. For samples embedded in paraffin, total RNA was isolated using the Roche High Pure FFPET RNA Isolation Kit (06650775001, (Merck Madrid, Spain) following the manufacturer’s protocol. cDNA synthesis and PCR conditions were performed as previously described [28]. Primers were designed by using the NCBI BLAST software and purchased from Sigma-Aldrich. The primers used are listed in Supplementary Table S2. 2.12. Histology and Immunohistochemistry Histology and immunohistochemistry assays were performed as previously described [15]. Human tumor samples were provided by the Tumor Banks of the Vall d’Hebron University Hospital and the General Hospital of Albacete with appropriate ethical approval (number 2019/07/071). 2.13. RNA Sequencing Cells infected with PLKO.1-empty vector or PLKO.1-shRNA ERK5-1 were grown up to 90% confluency after selection, split 1/3, and allowed to grow again. Two different pools for 3MC-C1 and three different pools for 3MC-C3 were used. For each pool of infection, new viruses were produced with at least 2 weeks of difference. Total RNA was extracted, and RNA integrity was determined by Agilent 2100 Bioanalyzer (RIN range, 9.1–9.9). cDNA libraries were generated from 600 ng of RNA from each sample using the Illumina TruSeq Stranded mRNA Sample Preparation kit (Illumina, San Diego, CA, USA). Libraries were Cancers 2022, 14, 3509 5 of 19 sequenced on an Ilumina HiSeq2000 apparatus according to the manufacturer’s protocols, and a minimum of 9 million 100-base single-reads were generated from each sample. Pseudocounts for each gene were obtained from sequencing reads with salmon soft- ware [31] using RefSeq [32] mRNA sequences for mouse genome assembly GRCm38/mm10 as reference. Differential expression in individual subsets was calculated with the R pack- age DESeq2 [33] using default settings. RNAseq data are publicly available from the NCBI’s GEO repository with accession number GSE199395. 2.14. Functional Enrichment Analysis Enrichment of Gene Ontology terms was performed with the Bioconductor’s cluster- Profiler package [34] using a q cut-off value of 0.05. 2.15. In Silico Analysis MAPKs expression in TGCA series was analyzed by using the Timer2.0 website. The top 30 up-regulated and down-regulated gene signature after Mapk7 abrogation was identified in terms of log2 of Fold Change, and prognosis correlation with expression levels was evaluated. To associate the level of gene expression with prognosis, the median expression values for different transcripts were used as a cut-off to discriminate “high” and “low” expression cohorts, which are compared using a Cox survival analysis (pro- portional hazards). Graphics were drawn using the ggplot2 package running in R Studio Version 1.2.5033. For ERK5 protein-protein interaction network analysis, the STRING website (https://www.string-db.org/cgi/input?sessionId=bfLuXEqaflvN&input_page_ show_search=off, accessed on 3 April 2022) was used. 2.16. Statistical Analysis The data are reported as the mean ± standard deviation (SD). Statistical analysis was performed using the GraphPadPrism 9 and Office Excel 2020 (Microsoft). Significance was determined using a t-test and a two-way ANOVA test. The statistical significance of differences is indicated in Figures by asterisks as follows: * p < 0.05; ** p < 0.01; and *** p < 0.001. 3. Results 3.1. ERK5 Is Upregulated in 3MC-Induced Sarcoma and in Human Leiomyosarcoma After 120 days of the 3MC injection in the left rear, mice developed a tumor mass between 1 and 1.5 cm in diameter. Then, mice were sacrificed, and tumors obtained were used for histologic assays, biochemical determinations and to derive stable cell lines. Histologically, the obtained tumors are solid, highly compacted, and infiltrating local muscle tissue. The tumor tissue is highly vascularized. The cells are spindle-shaped, with a high number of mitoses, large irregular nuclei, and almost without nucleoli (Figure 1A,B), and necrotic areas are present (Figure 1C). To properly identify the tumors generated, an immunohistochemical study was carried out, showing that these tumors are positive for vimentin, smooth muscle actin, and caldesmon and negative for desmin and S100 (Figure 1D–H). Moreover, the expression pattern of CD34 has almost a negative display with scarce positive cells (Figure 1I). All these results indicate that the type of sarcoma obtained in 3MC-induced carcinogenesis in mice with a C57BL/6J genetic background could be a pleomorphic sarcoma with muscle differentiation, which could be consistent with a human pathology as leiomyosarcoma. Cancers 2022, 14, x FOR PEER REVIEW 6 Cancers 2022, 14, 3509 6 of 19 pleomorphic sarcoma with muscle differentiation, which could be consistent with a man pathology as leiomyosarcoma. Figureand Figure 1. Histological 1. Histological and immunohistochemical immunohistochemical characterizationcharacterization of 3MC-induced of 3MC-induced tumors. Rep- tumors. R resentative hematoxylin and eosin-stained fixed sections of tumors induced by 3MC. (A) At resentative hematoxylin and eosin-stained fixed sections of tumors induced by 3MC. (A) At low magnification (10×), the tumors show high cellular density and infiltrated local muscle fibers. (B magnification (10× ), the high tumors show(40×), magnification high the cellular cellsdensity displayand infiltrated spindle local muscle or epithelial fibers. (B) morphology, At atypia, nuclear high magnification (40 × ), the cells display spindle or epithelial morphology, nuclear atypia, and frequent mitosis (indicated by yellow arrows). (C) There are large areas of necrotic cells locate frequent mitosis (indicated the center of bythe yellow arrows). tumor. Original(C) There are large magnification areas is 10×; of are insets necrotic shown cells at alocated at 40× magnification. R resentative the center of the tumor. imagesmagnification Original of immunohistochemical is 10×; insetsstaining in fixed are shown at section of tumors induced by 3 a 40× magnification. The tumor Representative images cells were positive forstaining of immunohistochemical vimentin in(D), fixedsmooth sectionmuscle actininduced of tumors (E), and caldesmon by 3MC. (F). Des (G), protein S100 (H), and CD34 (I) were negative. Pictures are shown at a 20× magnification. The tumor cells were positive for vimentin (D), smooth muscle actin (E), and caldesmon (F). Desmin (G), protein S100 (H), and CD34 (I) were negative. Pictures are shown at a 20× magnification. To characterize the MAPK signaling in 3MC-induced sarcomas, we initially analy in silico the expression To characterize the MAPK patterns signaling of differentsarcomas, in 3MC-induced MAPKs in wethe TGCAanalyzed initially series, showing h sarcoma is the pathology with the highest level of ERK5 (MAPK7) in silico the expression patterns of different MAPKs in the TGCA series, showing how expression, while differences for other MAPKs were observed (Supplementary sarcoma is the pathology with the highest level of ERK5 (MAPK7) expression, while noFigure S1). To extrapo this in silico observation to our model, we performed a western blotting assay of ER differences for other MAPKs were observed (Supplementary Figure S1). To extrapolate comparing 3MC-induced tumors with normal muscle in each case. As it is shown in this in silico observation to our model, we performed a western blotting assay of ERK5 ure 2A, a marked increase in ERK5 expression was observed in tumor samples, w comparing 3MC-induced tumors with normal muscle in each case. As it is shown in ERK2, the MAPK most closely related to ERK5, showed no differences. Furthermore, Figure 2A, a marked increase in ERK5 expression was observed in tumor samples, while observation was confirmed by RT-qPCR evaluating Mapk7 mRNA expression levels ERK2, the MAPK most closely related to ERK5, showed no differences. Furthermore, malized to Mapk1 (Figure 2B). To evaluate our observation of human pathology, a s this observation was confirmed by RT-qPCR evaluating Mapk7 mRNA expression levels normalized to Mapk1 (Figure 2B). To evaluate our observation of human pathology, a short collection of human uterine leiomyosarcomas (Supplementary Table S3) was analyzed by RT-qPCR compared to smooth muscle samples. As it is shown in Figure 2C, human leiomyosarcomas showed a significant increase in the mRNA level of MAPK7 compared Cancers 2022, 14, 3509 7 of 19 to controls. In fact, immunohistochemistry was performed, showing a clear expression of this MAPK in agreement with the RT-qPCR assay (Figure 2E). Furthermore, another sarcomatous pathology of muscular origin rhabdomyosarcoma (Supplementary Table S3), was also evaluated in terms of ERK5 expression showing an even more acute expression of ERK5 compared to control striated muscle (Figure 2D,F). Therefore, this set of experiments suggests that our observations about ERK5 in the murine experimental model of 3MC- induced sarcoma can be extrapolated to human pathologies, such as leiomyosarcoma Cancers 2022, 14, x FOR PEER REVIEW 8 of 22 or rhabdomyosarcoma. Figure 2. Figure 2. ERK5 ERK5levels areare levels increased increasedin 3MC-derived in 3MC-derived tumors and inand tumors human leiomyosarcomas in human and leiomyosarcomas and rhabdomyosarcomas. (A) Lysates from tumors (T) induced by 3MC and from control (N) normal rhabdomyosarcomas. (A) Lysates from tumors (T) induced by 3MC and from control muscles (n = 4) were collected, and protein extracts (100 µg) were blotted against ERK5 and ERK2 (N) normal muscles (n = 4) were MAPKs. Vinculin collected, was used and protein as a loading control. extracts (100analysis (B) RT-qPCR µg) were blotted of Mapk7 against mRNA ERK5 in the sameand ERK2 MAPKs. Vinculin four murine tumorswas usedby induced as3MC a loading control. and control (B) RT-qPCR normal analysis muscles. The mRNAof Mapk7 mRNA expression in the same of Mapk1 (ERK2) was used as an endogenous control. (C) RT-qPCR analysis of MAPK7 mRNA in human leiomyosarcomas and healthy smooth muscle samples. GAPDH was used as an endogenous control. (D) RT-qPCR analysis of MAPK7 mRNA in human rhabdomyosarcomas and healthy striated mus- cle samples, with GAPDH as an endogenous control. (E) Representative images of immunohisto- chemical ERK5 staining in fixed sections of one control smooth muscle and two human leiomyosar- comas. (F) Representative images of immunohistochemical ERK5 staining in fixed sections of one Cancers 2022, 14, 3509 8 of 19 four murine tumors induced by 3MC and control normal muscles. The mRNA expression of Mapk1 (ERK2) was used as an endogenous control. (C) RT-qPCR analysis of MAPK7 mRNA in human leiomyosarcomas and healthy smooth muscle samples. GAPDH was used as an endogenous control. (D) RT-qPCR analysis of MAPK7 mRNA in human rhabdomyosarcomas and healthy striated muscle samples, with GAPDH as an endogenous control. (E) Representative images of immunohistochemical ERK5 staining in fixed sections of one control smooth muscle and two human leiomyosarcomas. (F) Representative images of immunohistochemical ERK5 staining in fixed sections of one control striated muscle and two human rhabdomyosarcomas. Images are depicted at a 20× magnification. Original blots could be found in Figure S7. ** p < 0.01; and *** p < 0.001. 3.2. ERK5 Is a Major Determinant in the Biology of Sarcoma-Derived Cell Lines Induced by 3MC To study the role of ERK5 in our experimental model based on 3MC, we derived cell lines from different 3MC-induced tumors. Among the several cell lines generated, we chose two of them (referred to as 3MC-C1 and 3MC-C3) that showed a potent ability to induce in vivo tumors (Supplementary Figure S2A). Indeed, tumors induced by 3MC-C1 and 3MC-C3 cell lines showed a histological pattern almost identical to the original tumors induced by 3MC (Supplementary Figure S2B). As expected, both cell lines exhibited a lack of p53 functionality (Supplementary Figure S2C), supporting previous observations in 3MC-derived tumors and also observed in human sarcomatous pathology [13,14,35,36]. Moreover, both cell lines showed a functional ERK5 signaling pathway determined by EGF stimulation (Supplementary Figure S2D). To fully investigate the role of ERK5 in our murine sarcoma-derived cell lines, we used specific shRNA targeting the murine Mapk7 gene. After achieving an effective knock-down at the protein and RNA levels (Figure 3A), we evaluated the effects of ERK5 abrogated expression in terms of cell growth (Figure 3B), showing how the lack of ERK5 promoted a marked delay in cell growth, especially at later time points. Next, we decided to perform clonogenic assays to evaluate biological effects associated with ERK5 abrogation in terms of proliferation as single cells and its ability to grow as a colony. To this end, cells were seeded in 6-well plates at low density and allowed to grow for 12–14 days. Interestingly, cells with ERK5 abrogated expression exhibited a marked decrease in the number of colonies (Figure 3C), suggesting that ERK5 could have some effect not only in the clonogenic growth of our sarcoma-derived cells but also in the attachment ability, as it has been described in other experimental models [37,38]. Therefore, adhesion assays in fibronectin-coated wells were performed, but no significant differences were observed (Figure 3D). All these in vitro observations were further confirmed in the 3MC-C3 cell line (Figure 4A–D) as well as by using a second shRNA for Mapk7 (shERK5-2) (Supplementary Figure S3A–H), which showed nearly the same results. In sum, all the evidence suggests that ERK5 signaling is directly implicated in the growth and survival abilities of our sarcoma-derived cell lines. Next, we decided to evaluate in vivo the effect of ERK5 downmodulation. Therefore, new infection pools were obtained for both cell lines, showing an effective knock-down of ERK5 (data not shown) and injected in nude mice. As it is depicted (Figure 3E), the lack of ERK5 promoted a marked delay in the in vivo growth. In fact, when the control group was euthanized for ethical reasons on day 28 (average weight of tumors 0.32 g ± 0.9), tumors derived from the interfered lines were 5–7 times lower (Figure 3E). However, at day 39 (average weight tumors 0.31 g ± 0.13), animals inoculated with cells with abrogated ERK5 expression had developed tumors of similar size to controls, probably due to the increased proliferation rate of those cells in which the interference was less effective. In fact, analysis of tumors showed a marked recovery of ERK5 expression in tumors originating from cell lines with abrogated ERK5 expression, as well as no morphological differences between both groups (Figure 3F,G). All these observations were further confirmed in the 3MC-C3 cell line (Figure 4E–G) as well as by using a second shRNA for Mapk7 (shERK5-2) (Supplementary Figure S3I,J), which showed nearly the same results. Taking all together, these data support the critical role of ERK5 in cell growth in vivo. Cancers 2022, 14, 3509 9 of 19 Cancers 2022, 14, x FOR PEER REVIEW 10 of 22 Figure 3. ERK5 modulates in vitro and in vivo growth of 3MC-C1 cell line. (A) 3MC-C1 cells were Figure 3. lentiviruses infected with ERK5 modulates carrying PLKO.1 in vitro empty vectorand (E.V.) in vivo or the growthERK5-1 PLKO.1-shRNA of 3MC-C1 vec- cell line. (A) 3MC-C1 cells tor (shERK5-1). Interference was evaluated by RT-qPCR using β-2-microglobulin (B2m) as an en- were infected with lentiviruses carrying PLKO.1 empty vector (E.V.) dogenous control (left panel). E.V. cells were considered as 1. Right panel shows a representative or the PLKO.1-shRNA ERK5-1 image of the vector interference by western (shERK5-1). blot using vinculin Interference was asevaluated a loading control. by(B)RT-qPCR For growth curves, using β-2-microglobulin (B2m) as an endogenous control (left panel). E.V. cells were considered as 1. Right panel shows a representative image of the interference by western blot using vinculin as a loading control. (B) For growth curves, 3 × 105 E.V. or shERK51 3MC-C1 cells were seeded in 100 mm plates. Every 3 days, cells were counted and replated in the same manner up to day 9. Graphic shows the cumulative cell number from a representative experiment out of 3 with nearly identical results in different pools of infections. (C) Upper panel: Relative number of colonies obtained in clonogenic assays of E.V. and shERK5-1 3MC-C1 cells. Lower panel: Representative image of a colony formation assay from both cell lines. (D) Upper panel: Relative adhesion of E.V. or shERK5-1 3MC-C1 cells was assessed by crystal violet staining at the indicated time points. Lower panel: Representative image of adhesion assays at 7.5 and 60 min. (E) Nude mice (n = 4) were inoculated with 5 × 105 cells of each cell line derived from 3MC-C1, and tumorigenesis was analyzed at indicated times. Graphics represent the mean ± SD. The experiment was performed by using another pool of infection with nearly identical results. (F) Western blot analysis of the expression level of ERK5 in tumors recovered from 3MC-C1 E.V. and shERK5-1. Vinculin was used as a loading control. (G) Representative images of the histological study of tumors obtained from 3MC-C1 derived cell lines. Pictures are shown at a 20× magnification. The graphs/histograms represent the mean ± SD of 3 independent experiments performed in triplicate cultures with different pools of infections, if not otherwise indicated. Original blots could be found in Figure S7. * p < 0.05; ** p < 0.01; and *** p < 0.001. Cancers 2022, 14, 3509 10 of 19 Figure 4. ERK5 modulates in vitro and in vivo growth of 3MC-C3 cell line. (A) 3MC-C3 cells were infected with lentiviruses carrying PLKO.1 empty vector (E.V.) or the PLKO.1-shRNA ERK5-1 vector (shERK5-1). Interference was evaluated by RT-qPCR using β-2-microglobulin (B2m) as an endogenous control (left panel). E.V. cells were considered as 1. Right panel shows a representative image of the interference by western blot using vinculin as a loading control. (B) For growth curves, 3 × 105 E.V. or shERK5-1 3MC-C3 cells were seeded in 100 mm plates. Every 3 days, cells were counted and replated in the same manner up to day 9. Graphic shows the cumulative cell number from a representative experiment out of 3 with nearly identical results in different pools of infections. (C) Upper panel: Relative number of colonies obtained in clonogenic assays of E.V. and shERK5-1 3MC-C3 cells. Lower panel: Representative image of a colony formation assay from both cell lines. (D) Upper panel: Relative adhesion of E.V. or shERK5-1 3MC-C3 cells was assessed by crystal violet staining at the indicated time points. Lower panel: Representative image of adhesion assays at 7.5 and 60 min. (E) Nude mice (n = 4) were inoculated with 5 × 105 cells of each cell line derived from 3MC-C3, and tumorigenesis was analyzed at indicated times. Graphics represent the mean ± SD. Final tumor weight for E.V.-derived tumors was 0.27 g ± 0.17 and for shERK5-1-derived tumors was 0.26 ± 0.09). The experiment was performed by using another pool of infection with nearly identical results. (F) Western blot analysis of the expression level of ERK5 in tumors recovered from 3MC-C1 E.V. and shERK5-1. Vinculin was used as a loading control. (G) Representative images of the histological study of tumors obtained from 3MC-C3-derived cell lines. Pictures are shown at a 20× magnification. The graphs/histograms represent the mean ± SD of 3 independent experiments performed in triplicate cultures with different pools of infections, if not otherwise indicated. Original blots could be found in Figure S7. * p < 0.05; and *** p < 0.001. Cancers 2022, 14, 3509 11 of 19 3.3. ERK5 Modulates Transcriptional Landscape in Sarcoma-Derived Cell Lines Induced by 3MC In order to identify the mechanisms by which ERK5 affects tumor growth in vivo, we determined the gene expression profile of 3MC-C1 and 3MC-C3 cells with or without ERK5 expression by means of RNA sequencing. Transcriptomic analysis showed a total of 519 differentially expressed genes (DEG) (Figure 5A) (FDR < 0.05). There were 252 genes (2.1% of the total genome) whose expression was significantly increased (LFC > 0) and 267 repressed genes (LFC < 0, 2.2% of the total genome, see Gene Expression Omnibus under the accession code GSE199395). Validation of RNAseq was performed by means of RT-qPCR by testing the expression of 9 different genes previously reported to be related to ERK5 or cancer biology (Klf2, Rgs2, Vcan, Cbx6, Maff, Cepbd, Npnt, Cdk18, and Thbs1) plus Mapk7, in independent pools of infection from the 3MC-C1 and 3MC-C3 cell lines (Figure 5B). Functional enrichment analysis showed several Gene Ontology biological pro- cesses that were significantly affected by ERK5 abrogation, including adhesion, migration, muscle development, etc. (Figure 5C). Interestingly, despite not being included within the top 20 significant Gene Ontology biological processes, cell proliferation (GO:0008283, 20 out of 231 genes, p-value 1.45× 10−39 ) and regulation of cell growth (GO:0001558, 23 out 324 p-value 1.57× 10−3 ), among others, were also affected. These data reinforce our previ- ous in vitro observation of ERK5 role in cell proliferation. Next, we evaluated the potential of ERK5 expression as a biomarker for human sarcomatous pathology by performing in silico analysis in patients from the TGCA sarcoma series. As shown in Figure 5D, MAPK7 expression does not correlate significantly with prognosis. Next, the top 30 upregulated and downregulated genes associated with ERK5 abrogation were identified, and their corre- lation with prognosis was explored. Interestingly, both gene expression signatures showed a statistical correlation with prognosis, finding better prognosis in terms of overall survival in patients with high levels of the genes upregulated after ERK5 abrogation (Figure 5E) and in patients with low levels of the downregulated genes (Figure 5F). These correlations suggest a link between the experimental and the clinical setup. 3.4. KLF2 Is a Critical Mediator of the Biological Effects of ERK5 in Sarcoma-Derived Cell Lines Induced by 3MC Among the genes whose expression was significantly altered (FDR < 0.05) by ERK5 knock-down, we observed 37 mouse transcription factors [39], with five of them being strongly affected (absolute log2FC ≥ 1, Figure 6A). Next, we used a protein-protein inter- action network analysis based on the STRING website. This analysis indicates that only KLF2 interacts with ERK5 (Supplementary Figure S4), in agreement with previous reports showing that KLF2 is a known target of the ERK5 signaling pathway [40–42]. Indeed, it has been recently suggested that this transcription factor can also be a novel substrate of ERK5 kinase activity [43]. Therefore, to elucidate the role of this transcription factor in the biolog- ical effects associated with the ERK5 signaling pathway, we decided to downregulate Klf2 expression by shRNA in 3MC-C1 cells. After achieving an effective knock-down confirmed by RT-qPCR (Figure 6B), cells were challenged in terms of proliferation, reproductive capac- ity, adhesion (Figure 6C–E), and in vivo growth (Figure 6F–H), obtaining the same results that in the case of ERK5 abrogation or even more striking differences. In addition, similar results were obtained in the 3MC-C3 cell line (Supplementary Figure S5). Furthermore, to fully clarify whether the effect of KLF2 could be due to a feedback regulation between KLF2 and ERK5, the RNA and protein levels and functionality of ERK5 were analyzed in control and Klf2-interfered cells showing no differences (Supplementary Figure S6). In sum, this set of experiments reinforces the critical role of the signaling axis ERK5→KLF2 in the sarcoma biology induced by intramuscular injection of 3MC. Cancers 2022, 14, Cancers 2022, 14, 3509 x FOR PEER REVIEW 12 of 14 of 22 19 Figure Figure 5.5. ERK5 ERK5 mediates mediates transcriptional transcriptional landscape landscape in in 3MC-derived 3MC-derived cell cell lines. lines. (A) (A) Volcano Volcano plot plot showing the effect size of ERK5 attenuation on gene expression (log2 fold change in ERK5-depleted showing the effect size of ERK5 attenuation on gene expression (log2 fold change in ERK5-depleted over control cells) versus statistical significance of the difference in expression between conditions over control cells) versus statistical significance of the difference in expression between condi- (−log10 of the FDR-adjusted p-values). Each dot represents an independent gene, and DEG (FDR < 0.05) (−shown tions is log10 of in the red FDR-adjusted color. Text labels p-values). show the Each dot represents identity and locationan of independent gene, and the genes selected for DEG vali- (FDR < 0.05) is shown in red color. Text labels show the identity and location dation. (B) Effect of ERK5 attenuation on the expression of the genes selected for validation. The of the genes selected for validation. figure shows the (B)effect Effectsize of ERK5 (log2 attenuation fold changeon inthe expression ofover ERK5-depleted the genes control selected cells) asfordetermined validation. fromfigure The the RNAseq shows the (RNAseq) data effect size andfold (log2 from RT-qPCR change determinations in ERK5-depleted in cell over pools control derived cells) from two as determined independent from the RNAseq infections (3MC-C3 (RNAseq) data and and 3MC-C1). from RT-qPCR (C) Top 20 Gene Ontology determinations in cellterms pools significantly derived fromasso- two ciated with genes independent differentially infections (3MC-C3expressed and 3MC-C1). upon(C) ERK5 Top suppression. 20 Gene OntologyThe graph terms represents significantly theassoci- frac- tion ated of DEG with with genes the indicated differentially GO term. expressed upon The colorsuppression. ERK5 of the symbols The indicates the statistical graph represents signifi- the fraction of cance of the association (adjusted p-value) and their size, the absolute number of DEG genes with DEG with the indicated GO term. The color of the symbols indicates the statistical significance of the the indicated GO label. (D) Kaplan–Meyer comparing prognosis in terms of Overall Survival (OS) association (adjusted p-value) and their size, the absolute number of DEG genes with the indicated for two groups of patients, those with high (in red) and low (in black) expression levels of ERK5 for GO label.patients sarcoma (D) Kaplan–Meyer from TCGAcomparing prognosis dataset. Hazard in (HR) ratio termsandof Overall p-valueSurvival (log rank (OS) P) for two groups showed in the upper right corner of the panel. (E) Kaplan–Meyer comparing prognosis in terms of OS patients of patients, those with high (in red) and low (in black) expression levels of ERK5 for sarcoma for two from TCGA groups dataset.those of patients, Hazard withratio high(HR) andand (in red) p-value (logblack) low (in rank P) showed in expression the upper levels right corner of the upregulated genes of the associated panel. (E) with ERK5 abrogation Kaplan–Meyer comparing for sarcoma prognosis patients in termsfromof TCGA OS fordataset. two groups HR and p-value of patients, (log thoserank withP)high showed in the (in red) andupper low (in right corner black) of the panel. expression levels (F) Kaplan–Meyer of the upregulated comparing genes associatedprognosis with in terms ERK5 of OS forfor abrogation two groupspatients sarcoma of patients, from those TCGAwith highHR dataset. (in and red)p-value and low(log (inrank black)P) expression showed in levels of the the upper rightdownregulated genes (F) corner of the panel. associated with ERK5 Kaplan–Meyer abrogation comparing for sarcoma prognosis in terms patients of OS forfromtwo TCGA dataset. HR and p-value (log rank P) showed in the upper right corner of the panel. groups of patients, those with high (in red) and low (in black) expression levels of the downregulated genes associated with ERK5 abrogation for sarcoma patients from TCGA dataset. HR and p-value 3.4. KLF2 Is a Critical Mediator of the Biological Effects of ERK5 in Sarcoma-Derived Cell Lines (log rank P) showed in the upper right corner of the panel. Induced by 3MC Among the genes whose expression was significantly altered (FDR < 0.05) by ERK5 knock-down, we observed 37 mouse transcription factors [39], with five of them being strongly affected (absolute log2FC ≥ 1, Figure 6A). Next, we used a protein-protein inter- action network analysis based on the STRING website. This analysis indicates that only KLF2 interacts with ERK5 (Supplementary Figure S4), in agreement with previous reports showing that KLF2 is a known target of the ERK5 signaling pathway [40–42]. Indeed, it has been recently suggested that this transcription factor can also be a novel substrate of ERK5 kinase activity [43]. Therefore, to elucidate the role of this transcription factor in the biological effects associated with the ERK5 signaling pathway, we decided to downregu- late Klf2 expression by shRNA in 3MC-C1 cells. After achieving an effective knock-down confirmed by RT-qPCR (Figure 6B), cells were challenged in terms of proliferation, repro- ductive capacity, adhesion (Figure 6C–E), and in vivo growth (Figure 6F–H), obtaining the same results that in the case of ERK5 abrogation or even more striking differences. In addition, similar results were obtained in the 3MC-C3 cell line (Supplementary Figure S5). Furthermore, to fully clarify whether the effect of KLF2 could be due to a feedback regu- Cancers 2022, 14, 3509 lation between KLF2 and ERK5, the RNA and protein levels and functionality of ERK5 13 of 19 were analyzed in control and Klf2-interfered cells showing no differences (Supplementary Figure S6). In sum, this set of experiments reinforces the critical role of the signaling axis ERK5KLF2 in the sarcoma biology induced by intramuscular injection of 3MC. Figure 6. KLF2 is a key effector of the ERK5 biological properties in 3MC-C1 cell line. (A) Volcano Figure plot showing6. transcription KLF2 is factors a keyaffected effector by ERK5ofabrogation. the ERK5 biological Volcano properties plot showing the effect in 3MC-C1 cell line. (A) Volcano size of ERK5 attenuation on gene expression (log2 fold change in ERK5-depleted over control cells) plot versus showing transcription statistical significance of the differencefactors affected in expression by ERK5 between conditions (−log10 abrogation. of the FDR- Volcano plot showing the effect size of ERK5 attenuation on gene expression (log2 fold change in ERK5-depleted over control cells) versus statistical significance of the difference in expression between conditions (−log10 of the FDR- adjusted p-values). Each dot represents an independent gene. DEG as determined by statistical significance (FDR < 0.05), effect size (|Log2FC| > 0.99) or both are shown in blue, green, and red colors respectively. Text labels show the identity and location of the transcription factors significantly affected (FDR < 0.05 and |Log2FC| > 0.99) by ERK5. (B) RT-qPCR showing effective knock-down of Klf2 expression in 3MC-C1 cells infected with lentiviral vectors coding for empty vector (E.V.) or an shRNA specific for Klf2 (shKLF2), using β-2-microglobulin (B2m) as an endogenous control. E.V. cells were refereed as 1. Histogram shows the average of 3 independent pools of infection. (C) For growth curves, 3 × 105 E.V. or shKLF2 3MC-C1 cells were seeded in 100 mm plates. Every 3 days, cells were counted and replated in the same manner up to day 9. Graphic shows the cumulative cell number from a representative experiment out of 3 with nearly identical results in different pools of infections. (D) Upper panel: Relative number of colonies obtained in clonogenic assays of E.V. or shKLF2 3MC-C1 cells. Colonies formed by E.V. cells were considered as 1. Lower panel: Representative image of a colony formation assay at 7.5 and 60 min. (E) Upper panel: Relative adhesion of E.V. or shKLF2 3MC-C1 cells was assessed by crystal violet staining at different time points. Graphic shows mean of 3 independent experiments performed in triplicated cultures with 3 different pools of infection. Lower panel: Representative image of adhesion assays at 7.5 and 60 min. (F) Nude mice (n = 4) were inoculated with 5 × 105 E.V. and shKLF2 3MC-C1 cells, and tumor growth was analyzed at the indicated times. Final tumor weight for E.V.-derived tumors was 0.30 g ± 0.05 and for shLKF2-derived tumors was 0.29 ± 0.09. The experiment was performed by using another pool of infections with nearly identical results. (G) Representative images of the histological study of tumors obtained from E.V. and shKLF2 3MC-C1 cells. Pictures are depicted at a 20× magnification. (H) RNA from tumors induced by E.V. and shKLF2 3MC-C1 cells was extracted, and Klf2 expression levels were measured by RT-qPCR in triplicate using Mapk1 as an endogenous control. The graphs/histograms represent the mean ± SD of 3 independent experiments performed in triplicate cultures from different pools of infections, if not otherwise indicated. *** p < 0.001, **** p < 0.0001. Cancers 2022, 14, 3509 14 of 19 4. Discussion Several conclusions can be drawn from the present report. The first and most obvious is that 3MC-induced tumors are not sarcomas in a general sense as it has been previ- ously considered [13]. Our data support that the best filiation for these tumors should be pleomorphic sarcoma with muscle differentiation that could be consistent with human leiomyosarcoma, as the histologic and immunohistochemical studies support. Therefore, our work highlights this model as an appropriate tool to study the biology of muscular sarcomas as leiomyosarcomas, both showing common genetic alterations as in the case of Tp53 [13,44,45]. In addition, this model holds several advantages, such as the single intramuscular injection necessary to induce the tumor formation or the high effectiveness and reproducibility. However, it is remarkable that 3MC can promote rhabdomyosarcomas by subcutaneous injection in C3H/HeJ mice [46], suggesting that the use of 3MC could be a valuable tool for sarcoma research, with special interest for tumors with a muscular origin such as leiomyosarcoma or rhabdomyosarcoma. Regarding ERK5, it is notorious how tumors induced by 3MC showed a marked upregulation of this MAPK at the protein and RNA levels, which is also observed in a human sarcomatous pathology such as leiomyosarcoma. Our data from 3MC-derived cell lines suggest that sarcoma, or at least pathologies such as leiomyosarcoma, should be included in the list of tumors in which ERK5 plays a determinant role, as for example, hepatocellular carcinoma [47], breast cancer [38,48], lung cancer [49], or myeloma [26], among others. In fact, we were able to detect a potent expression of this MAPK in a short series of rhabdomyosarcomas, indicating that further study is necessary to assess whether our observation is restricted to sarcoma of the muscular origin or can also be extrapolated to other sarcomatous pathologies. Regarding cell proliferation, although some reports suggest a lack of implication of ERK5 [50–52], our results obtained in growth curves and foci assays reinforce the critical role of this MAPK in cell proliferation and survival, as it has been previously described [27,49,53]. However, we cannot discard that cell type/genetic background issues, as well as different methodological approaches, could be critical to ex- plaining the discrepancy between the different reports. Nonetheless, the oncogenic effect of ERK5 in vivo, regardless of their effects on proliferation and survival, seems to be universal and could be explained by different mechanisms. The fact that key biological processes such as motility, adhesion, blood vessel development, or angiogenesis could be affected by the lack of ERK5 fits perfectly with our observations, supporting a definitive role of ERK5 in tumor growth and progression [23]. However, some of these putative explanations can be discarded in our 3MC-based experimental model of sarcoma, for example, the effect on cell attachment described in other experimental models [38,51]. On the contrary, other possibilities, such as the role of ERK5 signaling pathway in angiogenesis [54–56] or even in cancer stem cell biology [57–59], could be considered. Therefore, further studies are necessary to fully clarify the biological mechanisms by which ERK5 is implicated in sarcomatous pathology, as could be the case of leiomyosarcomas or rhabdomyosarcomas. Regarding KLF2, our data show that this transcription factor is a key component in the role of ERK5 signaling pathway, being able to mimic the effects associated with ERK5 in our experimental model. It is notorious how KLF2 seems to play a dual role depending on the model under study. Thus, on the one hand, there are some studies associating KLF2 with an oncogenic role, as in the case of multiple myeloma [60] or pancreatic tumors [61], whereas in other models, it seems to have a character closer to tumor suppression, as in hepatocarcinoma [62,63] or, more recently, in renal cancer [64]. Of note, Klf2 seems to have oncogenic potential in our experimental model, being this report the first connection between this gene and sarcoma biology. Furthermore, KLF2, as an ERK5 downstream target, could be implicated in cancer through different mechanisms that remain unexplored, as could be the case of stem cell biology [65,66] or other processes such as inflammatory response [67], immune response [68], or angiogenesis [61,69]. Therefore, one possibility could be that the ERK5–KLF2 axis acts as a node in sarcoma biology by controlling several key processes. Cancers 2022, 14, 3509 15 of 19 Finally, from a clinical point of view, our murine model opens a new window of thera- peutic and diagnostic possibilities for sarcomatous pathologies based on ERK5-dependent signaling [24]. In this regard, the effectivity of ERK5 inhibitors in preclinical models has been shown [49,70,71]. However, the use of ERK5 inhibitors should be carefully considered, as it has been described that they can induce the opposite effect, meaning hyperactiva- tion of the targeted signaling pathway, as in the case of ERK1/2 [72], AKT [73], or even ERK5 [74]. From the diagnostic point of view, our in silico data suggest that ERK5 could be a weak biomarker, supporting that high levels of ERK5 may be a constitutive feature of sarcomatous pathology. However, the gene expression signature associated with ERK5 abrogation showed a strong correlation with better outcomes in the TCGA sarcoma series. This approximation, which considers the joint expression levels of 30 subrogated markers, avoids sampling effects due to the reduced size of the cohort evaluated. Therefore, the obtained results suggest that downstream effectors could be potential novel biomarkers for sarcomatous pathology. Nonetheless, the potential of ERK5 and its substrates as putative biomarkers or novel targets for therapy is a challenging issue in sarcoma biology that needs to be fully clarified. 5. Conclusions In summary, the present report demonstrates a critical role for ERK5 in murine carcinogen-induced sarcoma, which correlates with an overexpression of this MAPK, also observed in human pathologies such as leiomyosarcoma or rhabdomyosarcoma. Moreover, the use of NGS allows us to conclude that ERK5 is a modulator of the transcriptional pattern in sarcoma-derived cell lines induced by 3MC and that its biological activity is mediated through the effect exerted onto KLF2. Therefore, all the previous results reveal the necessity of a more in-depth study to exploit the diagnostic and therapeutic potential of ERK5 and its downstream targets in human sarcoma. Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/cancers14143509/s1, Figure S1: In silico analysis of RNA expression patterns of MAPK1 (ERK2), MAPK3 (ERK1), MAPK8 (JNK1), MAPK14 (p38α) and MAPK7 (ERK5) in the TGCA series, Figure S2: Characterization of cell lines from 3MC-induced tumors, Figure S3: A second shRNA for ERK5 (TRC-N-0000232397) showed identical behavior in 3MC-C1 and 3MC-C3 cell lines, Figure S4: Protein-protein interaction network, Figure S5: Abrogation of KLF2 expression in 3MC-3C cell line, Figure S6: Lack of KLF2 does not affect ERK5 expression or functionality, Figure S7: original blots; Table S1: List of antibodies used in this work, Table S2: List of RT-qPCR used in this report, Table S3: Data of patient studied, Table S4: Excel files with differential expressed genes in the RNAseq analysis. Author Contributions: R.S.-P. and M.J.R.-H.: designed research, analyzed data, obtained funding, wrote original draft, and edited manuscript. E.A.-L. and J.J.-S.: led investigation, corrected, and edited manuscript, analyzed data, designed, and developed Figures A.E.-O., A.P. and B.B.: analyzed data, corrected the manuscript, and helped in the design of research. R.P.-S., M.O.-M., D.M.F.-A., O.R., N.G.-F., C.G.-G., R.M.V., A.S.-F. and S.M.-A.: performed experiments in vitro and in vivo, analyzed data and corrected the manuscript. T.M., S.H.N.-C., G.F., J.H.-L. and S.R.y.C.: performed histological and immunohistochemical studies, corrected the manuscript, analyzed data, and helped in the design of research. L.d.P., V.J.S.-A.L. and F.J.C.: performed RNAseq analysis and in silico searches, analyzed data, corrected the manuscript, and helped in the design of research. All authors have read and agreed to the published version of the manuscript. Funding: This work has been supported with Grant RTI2018-094093-B-I00 funded by MCIN/AEI/10.1 3039/501100011033, “ERDF A way of making Europe” to RSP. Also supported with funds from Fun- dación Leticia Castillejo Castillo, Roche España and ACEPAIN to RSP and MJRH. RSP and MJRH’s Research Institute and the work carried out in their laboratory, received partial support from the Eu- ropean Community through the FEDER. JJ and EAL hold a predoctoral research contract cofounded by the European Social Fund and UCLM. OR holds a contract for accessing the Spanish System of Science, Technology, and Innovation (SECTI) funded by the University of Castilla-La Mancha (UCLM) Cancers 2022, 14, 3509 16 of 19 and received partial support from the European Social Fund (FSE) through its Operative Program for Castilla-La Mancha (2007–2013). Institutional Review Board Statement: All the animal experimentation was carried out according to the NIH-Intramural Animal Guide for the Care and Use of Laboratory Animals and approved by the Ethics in Animal Care Committee of the University of Castilla-La Mancha (reference PR-2019-07-18). Human tumor samples were provided by the Tumor Banks of the Vall d’Hebron University Hospital and the General Hospital of Albacete with appropriate ethical approval (number 2019/07/071). Informed Consent Statement: Informed consent was obtained from all subjects involved in the study. Data Availability Statement: All data and materials are available on reasonable request. RNAseq analysis data for 3MC-C1 and 3MC-C3 cells with abrogated ERK5 expression have been deposited in the Gene Expression Omnibus (GEO) under the accession code GSE199395. Conflicts of Interest: The authors declare no conflict of interest. References 1. Blay, J.-Y.; Kang, Y.-K.; Nishida, T.; von Mehren, M. Gastrointestinal stromal tumours. Nat. Rev. Dis. Primers 2021, 7, 22. [CrossRef] [PubMed] 2. Grünewald, T.G.; Alonso, M.; Avnet, S.; Banito, A.; Burdach, S.; Cidre-Aranaz, F.; Di Pompo, G.; Distel, M.; Dorado-Garcia, H.; Garcia-Castro, J.; et al. Sarcoma treatment in the era of molecular medicine. EMBO Mol. Med. 2020, 12, e11131. [CrossRef] [PubMed] 3. Dodd, R.D.; Mito, J.K.; Kirsch, D.G. Animal models of soft-tissue sarcoma. Dis. Model. Mech. 2010, 3, 557–566. [CrossRef] 4. Sampson, V.B.; Kamara, D.F.; Kolb, E.A. Xenograft and genetically engineered mouse model systems of osteosarcoma and Ewing’s sarcoma: Tumor models for cancer drug discovery. Expert Opin. Drug Discov. 2013, 8, 1181–1189. [CrossRef] [PubMed] 5. Dodd, R.D.; Añó, L.; Blum, J.M.; Li, Z.; Van Mater, D.; Kirsch, D.G. Methods to generate genetically engineered mouse models of soft tissue sarcoma. Methods Mol. Biol. 2015, 1267, 283–295. [CrossRef] 6. Donehower, L.A.; Harvey, M.; Slagle, B.L.; McArthur, M.J.; Montgomery, C.A.; Butel, J.S.; Bradley, A. Mice deficient for p53 are developmentally normal but susceptible to spontaneous tumours. Nature 1992, 356, 215–221. [CrossRef] 7. Hernando, E.; Charytonowicz, E.; Dudas, M.E.; Menendez, S.; Matushansky, I.; Mills, J.; Socci, N.D.; Behrendt, N.; Ma, L.; Maki, R.G.; et al. The AKT-mTOR pathway plays a critical role in the development of leiomyosarcomas. Nat. Med. 2007, 13, 748–753. [CrossRef] 8. Langdon, C.G.; Gadek, K.E.; Garcia, M.R.; Evans, M.K.; Reed, K.B.; Bush, M.; Hanna, J.A.; Drummond, C.J.; Maguire, M.C.; Leavey, P.J.; et al. Synthetic essentiality between PTEN and core dependency factor PAX7 dictates rhabdomyosarcoma identity. Nat. Commun. 2021, 12, 5520. [CrossRef] 9. Syverton, J.T.; Berry, G.P. Sarcomata and carcinomata induced in cottontail rabbits by methylcholanthrene. Science 1940, 91, 298. [CrossRef] 10. Efeyan, A.; Murga, M.; Martinez-Pastor, B.; Ortega-Molina, A.; Soria, R.; Collado, M.; Fernandez-Capetillo, O.; Serrano, M. Limited role of murine ATM in oncogene-induced senescence and p53-dependent tumor suppression. PLoS ONE 2009, 4, e5475. [CrossRef] 11. Narlik-Grassow, M.; Blanco-Aparicio, C.; Cecilia, Y.; Peregrina, S.; Garcia-Serelde, B.; Muñoz-Galvan, S.; Cañamero, M.; Carnero, A. The essential role of PIM kinases in sarcoma growth and bone invasion. Carcinogenesis 2012, 33, 1479–1486. [CrossRef] [PubMed] 12. Remeseiro, S.; Cuadrado, A.; Gómez-López, G.; Pisano, D.G.; Losada, A. A unique role of cohesin-SA1 in gene regulation and development. EMBO J. 2012, 31, 2090–2102. [CrossRef] [PubMed] 13. Lee, C.-L.; Mowery, Y.M.; Daniel, A.R.; Zhang, D.; Sibley, A.B.; Delaney, J.R.; Wisdom, A.J.; Qin, X.; Wang, X.; Caraballo, I.; et al. Mutational landscape in genetically engineered, carcinogen-induced, and radiation-induced mouse sarcoma. JCI Insight 2019, 4, 128698. [CrossRef] 14. Abeshouse, A.; Adebamowo, C.; Adebamowo, S.N.; Akbani, R.; Akeredolu, T.; Ally, A.; Anderson, M.L.; Anur, P.; Appelbaum, E.L.; Armenia, J.; et al. Comprehensive and Integrated Genomic Characterization of Adult Soft Tissue Sarcomas. Cell 2017, 171, 950–965.e28. [CrossRef] [PubMed] 15. Serrano, C.; Romagosa, C.; Hernández-Losa, J.; Simonetti, S.; Valverde, C.; Moliné, T.; Somoza, R.; Pérez, M.; Vélez, R.; Vergés, R.; et al. RAS/MAPK pathway hyperactivation determines poor prognosis in undifferentiated pleomorphic sarcomas. Cancer 2016, 122, 99–107. [CrossRef] 16. Czarnecka, A.M.; Synoradzki, K.; Firlej, W.; Bartnik, E.; Sobczuk, P.; Fiedorowicz, M.; Grieb, P.; Rutkowski, P. Molecular Biology of Osteosarcoma. Cancers 2020, 12, 2130. [CrossRef] 17. Lin, Y.K.; Wu, W.; Ponce, R.K.; Kim, J.W.; Okimoto, R.A. Negative MAPK-ERK regulation sustains CIC-DUX4 oncoprotein expression in undifferentiated sarcoma. Proc. Natl. Acad. Sci. USA 2020, 117, 20776–20784. [CrossRef] 18. Liu, H.; Nazmun, N.; Hassan, S.; Liu, X.; Yang, J. BRAF mutation and its inhibitors in sarcoma treatment. Cancer Med. 2020, 9, 4881–4896. [CrossRef] Cancers 2022, 14, 3509 17 of 19 19. Nagabushan, S.; Lau, L.M.S.; Barahona, P.; Wong, M.; Sherstyuk, A.; Marshall, G.M.; Tyrrell, V.; Wegner, E.A.; Ekert, P.G.; Cowley, M.J.; et al. Efficacy of MEK inhibition in a recurrent malignant peripheral nerve sheath tumor. NPJ Precis. Oncol. 2021, 5, 9. [CrossRef] 20. Zhou, G.; Bao, Z.Q.; Dixon, J.E. Components of a new human protein kinase signal transduction pathway. J. Biol. Chem. 1995, 270, 12665–12669. [CrossRef] 21. Lee, J.D.; Ulevitch, R.J.; Han, J. Primary structure of BMK1: A new mammalian map kinase. Biochem. Biophys. Res. Commun. 1995, 213, 715–724. [CrossRef] [PubMed] 22. Hoang, V.T.; Yan, T.J.; Cavanaugh, J.E.; Flaherty, P.T.; Beckman, B.S.; Burow, M.E. Oncogenic signaling of MEK5-ERK5. Cancer Lett. 2017, 392, 51–59. [CrossRef] 23. Stecca, B.; Rovida, E. Impact of ERK5 on the Hallmarks of Cancer. Int. J. Mol. Sci. 2019, 20, 1426. [CrossRef] [PubMed] 24. Monti, M.; Celli, J.; Missale, F.; Cersosimo, F.; Russo, M.; Belloni, E.; Di Matteo, A.; Lonardi, S.; Vermi, W.; Ghigna, C.; et al. Clinical Significance and Regulation of ERK5 Expression and Function in Cancer. Cancers 2022, 14, 348. [CrossRef] 25. Pavan, S.; Meyer-Schaller, N.; Diepenbruck, M.; Kalathur, R.K.R.; Saxena, M.; Christofori, G. A kinome-wide high-content siRNA screen identifies MEK5-ERK5 signaling as critical for breast cancer cell EMT and metastasis. Oncogene 2018, 37, 4197–4213. [CrossRef] [PubMed] 26. Carvajal-Vergara, X.; Tabera, S.; Montero, J.C.; Esparís-Ogando, A.; López-Pérez, R.; Mateo, G.; Gutiérrez, N.; Parmo-Cabañas, M.; Teixidó, J.; San Miguel, J.F.; et al. Multifunctional role of Erk5 in multiple myeloma. Blood 2005, 105, 4492–4499. [CrossRef] 27. Gentilini, A.; Lori, G.; Caligiuri, A.; Raggi, C.; Di Maira, G.; Pastore, M.; Piombanti, B.; Lottini, T.; Arcangeli, A.; Madiai, S.; et al. Extracellular Signal-Regulated Kinase 5 Regulates the Malignant Phenotype of Cholangiocarcinoma Cells. Hepatology 2021, 74, 2007–2020. [CrossRef] 28. Pascual-Serra, R.; Fernández-Aroca, D.M.; Sabater, S.; Roche, O.; Andrés, I.; Ortega-Muelas, M.; Arconada-Luque, E.; Garcia- Flores, N.; Bossi, G.; Belandia, B.; et al. p38β (MAPK11) mediates gemcitabine-associated radiosensitivity in sarcoma experimental models. Radiother. Oncol. 2021, 156, 136–144. [CrossRef] 29. Ortega-Muelas, M.; Roche, O.; Fernández-Aroca, D.M.; Encinar, J.A.; Albandea-Rodríguez, D.; Arconada-Luque, E.; Pascual-Serra, R.; Muñoz, I.; Sánchez-Pérez, I.; Belandia, B.; et al. ERK5 signalling pathway is a novel target of sorafenib: Implication in EGF biology. J. Cell Mol. Med. 2021, 25, 10591–10603. [CrossRef] 30. Campaner, S.; Doni, M.; Hydbring, P.; Verrecchia, A.; Bianchi, L.; Sardella, D.; Schleker, T.; Perna, D.; Tronnersjö, S.; Murga, M.; et al. Cdk2 suppresses cellular senescence induced by the c-myc oncogene. Nat. Cell Biol. 2010, 12, 54–59. [CrossRef] 31. Patro, R.; Duggal, G.; Love, M.I.; Irizarry, R.A.; Kingsford, C. Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods 2017, 14, 417–419. [CrossRef] [PubMed] 32. O’Leary, N.A.; Wright, M.W.; Brister, J.R.; Ciufo, S.; Haddad, D.; McVeigh, R.; Rajput, B.; Robbertse, B.; Smith-White, B.; Ako-Adjei, D.; et al. Reference sequence (RefSeq) database at NCBI: Current status, taxonomic expansion, and functional annotation. Nucleic Acids Res. 2016, 44, D733–D745. [CrossRef] [PubMed] 33. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [CrossRef] [PubMed] 34. Yu, G.; Wang, L.-G.; Han, Y.; He, Q.-Y. clusterProfiler: An R package for comparing biological themes among gene clusters. OMICS 2012, 16, 284–287. [CrossRef] 35. Choi, J.; Manzano, A.; Dong, W.; Bellone, S.; Bonazzoli, E.; Zammataro, L.; Yao, X.; Deshpande, A.; Zaidi, S.; Guglielmi, A.; et al. Integrated mutational landscape analysis of uterine leiomyosarcomas. Proc. Natl. Acad. Sci. USA 2021, 118, e2025182118. [CrossRef] 36. Chudasama, P.; Mughal, S.S.; Sanders, M.A.; Hübschmann, D.; Chung, I.; Deeg, K.I.; Wong, S.-H.; Rabe, S.; Hlevnjak, M.; Zapatka, M.; et al. Integrative genomic and transcriptomic analysis of leiomyosarcoma. Nat. Commun. 2018, 9, 144. [CrossRef] 37. Jiang, W.; Cai, F.; Xu, H.; Lu, Y.; Chen, J.; Liu, J.; Cao, N.; Zhang, X.; Chen, X.; Huang, Q.; et al. Extracellular signal regulated kinase 5 promotes cell migration, invasion and lung metastasis in a FAK-dependent manner. Protein Cell 2020, 11, 825–845. [CrossRef] 38. Xu, Q.; Zhang, J.; Telfer, B.A.; Zhang, H.; Ali, N.; Chen, F.; Risa, B.; Pearson, A.J.; Zhang, W.; Finegan, K.G.; et al. The extracellular- regulated protein kinase 5 (ERK5) enhances metastatic burden in triple-negative breast cancer through focal adhesion protein kinase (FAK)-mediated regulation of cell adhesion. Oncogene 2021, 40, 3929–3941. [CrossRef] 39. Zhou, Q.; Liu, M.; Xia, X.; Gong, T.; Feng, J.; Liu, W.; Liu, Y.; Zhen, B.; Wang, Y.; Ding, C.; et al. A mouse tissue transcription factor atlas. Nat. Commun. 2017, 8, 15089. [CrossRef] 40. Ananieva, O.; Macdonald, A.; Wang, X.; McCoy, C.E.; McIlrath, J.; Tournier, C.; Arthur, J.S.C. ERK5 regulation in naïve T-cell activation and survival. Eur. J. Immunol. 2008, 38, 2534–2547. [CrossRef] 41. Paez-Mayorga, J.; Chen, A.L.; Kotla, S.; Tao, Y.; Abe, R.J.; He, E.D.; Danysh, B.P.; Hofmann, M.-C.C.; Le, N.-T. Ponatinib Activates an Inflammatory Response in Endothelial Cells via ERK5 SUMOylation. Front. Cardiovasc. Med. 2018, 5, 125. [CrossRef] 42. Paolini, A.; Fontana, F.; Pham, V.-C.; Rödel, C.J.; Abdelilah-Seyfried, S. Mechanosensitive Notch-Dll4 and Klf2-Wnt9 signaling pathways intersect in guiding valvulogenesis in zebrafish. Cell Rep. 2021, 37, 109782. [CrossRef] [PubMed] 43. Brown, H.A.; Williams, C.A.C.; Zhou, H.; Rios-Szwed, D.; Fernandez-Alonso, R.; Mansoor, S.; McMulkin, L.; Toth, R.; Gourlay, R.; Peltier, J.; et al. An ERK5-KLF2 signalling module regulates early embryonic gene expression and telomere rejuvenation in stem cells. Biochem. J. 2021, 478, 4119–4136. [CrossRef] [PubMed] Cancers 2022, 14, 3509 18 of 19 44. Astolfi, A.; Nannini, M.; Indio, V.; Schipani, A.; Rizzo, A.; Perrone, A.M.; De Iaco, P.; Pirini, M.G.; De Leo, A.; Urbini, M.; et al. Genomic Database Analysis of Uterine Leiomyosarcoma Mutational Profile. Cancers 2020, 12, 2126. [CrossRef] [PubMed] 45. Cuppens, T.; Moisse, M.; Depreeuw, J.; Annibali, D.; Colas, E.; Gil-Moreno, A.; Huvila, J.; Carpén, O.; Zikán, M.; Matias-Guiu, X.; et al. Integrated genome analysis of uterine leiomyosarcoma to identify novel driver genes and targetable pathways. Int. J. Cancer 2018, 142, 1230–1243. [CrossRef] 46. Wu, H.; Inoue, M. Immunohistochemical analysis for Mdm2 and p53 proteins in methylcholanthrene-induced mouse rhab- domyosarcomas. J. Vet. Med. Sci. 2006, 68, 427–431. [CrossRef] 47. Rovida, E.; Di Maira, G.; Tusa, I.; Cannito, S.; Paternostro, C.; Navari, N.; Vivoli, E.; Deng, X.; Gray, N.S.; Esparís-Ogando, A.; et al. The mitogen-activated protein kinase ERK5 regulates the development and growth of hepatocellular carcinoma. Gut 2015, 64, 1454–1465. [CrossRef] 48. Esparís-Ogando, A.; Díaz-Rodríguez, E.; Montero, J.C.; Yuste, L.; Crespo, P.; Pandiella, A. Erk5 participates in neuregulin signal transduction and is constitutively active in breast cancer cells overexpressing ErbB2. Mol. Cell Biol. 2002, 22, 270–285. [CrossRef] 49. Sánchez-Fdez, A.; Re-Louhau, M.F.; Rodríguez-Núñez, P.; Ludeña, D.; Matilla-Almazán, S.; Pandiella, A.; Esparís-Ogando, A. Clinical, genetic and pharmacological data support targeting the MEK5/ERK5 module in lung cancer. NPJ Precis. Oncol. 2021, 5, 78. [CrossRef] 50. Ali, M.; Mutahir, Z.; Riaz, A. CRISPR/Cas9 engineering of ERK5 identifies its FAK/PYK2 dependent role in adhesion-mediated cell survival. Biochem. Biophys. Res. Commun. 2019, 513, 179–185. [CrossRef] 51. Matossian, M.D.; Hoang, V.T.; Burks, H.E.; La, J.; Elliott, S.; Brock, C.; Rusch, D.B.; Buechlein, A.; Nephew, K.P.; Bhatt, A.; et al. Constitutive activation of MEK5 promotes a mesenchymal and migratory cell phenotype in triple negative breast cancer. Oncoscience 2021, 8, 64–71. [CrossRef] [PubMed] 52. Nguyen, D.; Lemos, C.; Wortmann, L.; Eis, K.; Holton, S.J.; Boemer, U.; Moosmayer, D.; Eberspaecher, U.; Weiske, J.; Lechner, C.; et al. Discovery and Characterization of the Potent and Highly Selective (Piperidin-4-yl)pyrido[3,2-d]pyrimidine Based in Vitro Probe BAY-885 for the Kinase ERK5. J. Med. Chem. 2019, 62, 928–940. [CrossRef] [PubMed] 53. Tusa, I.; Gagliardi, S.; Tubita, A.; Pandolfi, S.; Menconi, A.; Lulli, M.; Dello Sbarba, P.; Stecca, B.; Rovida, E. The Hedgehog-GLI Pathway Regulates MEK5-ERK5 Expression and Activation in Melanoma Cells. Int. J. Mol. Sci. 2021, 22, 11259. [CrossRef] [PubMed] 54. Doebele, R.C.; Schulze-Hoepfner, F.T.; Hong, J.; Chlenski, A.; Zeitlin, B.D.; Goel, K.; Gomes, S.; Liu, Y.; Abe, M.K.; Nor, J.E.; et al. A novel interplay between Epac/Rap1 and mitogen-activated protein kinase kinase 5/extracellular signal-regulated kinase 5 (MEK5/ERK5) regulates thrombospondin to control angiogenesis. Blood 2009, 114, 4592–4600. [CrossRef] [PubMed] 55. Hultgren, N.W.; Fang, J.S.; Ziegler, M.E.; Ramirez, R.N.; Phan, D.T.T.; Hatch, M.M.S.; Welch-Reardon, K.M.; Paniagua, A.E.; Kim, L.S.; Shon, N.N.; et al. Slug regulates the Dll4-Notch-VEGFR2 axis to control endothelial cell activation and angiogenesis. Nat. Commun. 2020, 11, 5400. [CrossRef] 56. Radu, M.; Lyle, K.; Hoeflich, K.P.; Villamar-Cruz, O.; Koeppen, H.; Chernoff, J. p21-Activated Kinase 2 Regulates Endothelial Development and Function through the Bmk1/Erk5 Pathway. Mol. Cell Biol. 2015, 35, 3990–4005. [CrossRef] 57. Kedika, S.R.; Shukla, S.P.; Udugamasooriya, D.G. Design of a dual ERK5 kinase activation and autophosphorylation inhibitor to block cancer stem cell activity. Bioorg. Med. Chem. Lett. 2020, 30, 127552. [CrossRef] 58. Pereira, D.M.; Gomes, S.E.; Borralho, P.M.; Rodrigues, C.M.P. MEK5/ERK5 activation regulates colon cancer stem-like cell properties. Cell Death Discov. 2019, 5, 68. [CrossRef] 59. Wang, L.; Xiong, X.; Yao, Z.; Zhu, J.; Lin, Y.; Lin, W.; Li, K.; Xu, X.; Guo, Y.; Chen, Y.; et al. Chimeric RNA ASTN2-PAPPAas aggravates tumor progression and metastasis in human esophageal cancer. Cancer Lett. 2021, 501, 1–11. [CrossRef] 60. Salama, Y.; Heida, A.H.; Yokoyama, K.; Takahashi, S.; Hattori, K.; Heissig, B. The EGFL7-ITGB3-KLF2 axis enhances survival of multiple myeloma in preclinical models. Blood Adv. 2020, 4, 1021–1037. [CrossRef] 61. Wei, R.; Kong, L.; Xiao, Y.; Yuan, H.; Song, Y.; Wang, J.; Yu, H.; Mao, S.; Xu, W. CDK8 regulates the angiogenesis of pancreatic cancer cells in part via the CDK8-β-catenin-KLF2 signal axis. Exp. Cell Res. 2018, 369, 304–315. [CrossRef] [PubMed] 62. Huang, M.-D.; Chen, W.-M.; Qi, F.-Z.; Sun, M.; Xu, T.-P.; Ma, P.; Shu, Y.-Q. Long non-coding RNA TUG1 is up-regulated in hepatocellular carcinoma and promotes cell growth and apoptosis by epigenetically silencing of KLF2. Mol. Cancer 2015, 14, 165. [CrossRef] [PubMed] 63. Li, Y.; Tu, S.; Zeng, Y.; Zhang, C.; Deng, T.; Luo, W.; Lian, L.; Chen, L.; Xiong, X.; Yan, X. KLF2 inhibits TGF-β-mediated cancer cell motility in hepatocellular carcinoma. Acta Biochim. Biophys. Sin. 2020, 52, 485–494. [CrossRef] [PubMed] 64. Lu, Y.; Qin, H.; Jiang, B.; Lu, W.; Hao, J.; Cao, W.; Du, L.; Chen, W.; Zhao, X.; Guo, H. KLF2 inhibits cancer cell migration and invasion by regulating ferroptosis through GPX4 in clear cell renal cell carcinoma. Cancer Lett. 2021, 522, 1–13. [CrossRef] 65. Morikawa, M.; Koinuma, D.; Mizutani, A.; Kawasaki, N.; Holmborn, K.; Sundqvist, A.; Tsutsumi, S.; Watabe, T.; Aburatani, H.; Heldin, C.-H.; et al. BMP Sustains Embryonic Stem Cell Self-Renewal through Distinct Functions of Different Krüppel-like Factors. Stem Cell Rep. 2016, 6, 64–73. [CrossRef] 66. Wu, N.; Sun, H.; Zhao, X.; Zhang, Y.; Tan, J.; Qi, Y.; Wang, Q.; Ng, M.; Liu, Z.; He, L.; et al. MAP3K2-regulated intestinal stromal cells define a distinct stem cell niche. Nature 2021, 592, 606–610. [CrossRef] 67. Turpaev, K.T. Transcription Factor KLF2 and Its Role in the Regulation of Inflammatory Processes. Biochemistry 2020, 85, 54–67. [CrossRef] 68. Wittner, J.; Schuh, W. Krüppel-like Factor 2 (KLF2) in Immune Cell Migration. Vaccines 2021, 9, 1171. [CrossRef] Cancers 2022, 14, 3509 19 of 19 69. Chen, G.-H.; Li, X.-L.; Deng, Y.-Q.; Zhou, F.-M.; Zou, W.-Q.; Jiang, W.-X.; Shangguan, S.-Q.; Lu, Z.-N. The Molecular Mechanism of EPO Regulates the Angiogenesis after Cerebral Ischemia through AMPK-KLF2 Signaling Pathway. Crit. Rev. Eukaryot. Gene Expr. 2019, 29, 105–112. [CrossRef] 70. Myers, S.M.; Miller, D.C.; Molyneux, L.; Arasta, M.; Bawn, R.H.; Blackburn, T.J.; Cook, S.J.; Edwards, N.; Endicott, J.A.; Golding, B.T.; et al. Identification of a novel orally bioavailable ERK5 inhibitor with selectivity over p38α and BRD4. Eur. J. Med. Chem. 2019, 178, 530–543. [CrossRef] 71. Yang, Q.; Deng, X.; Lu, B.; Cameron, M.; Fearns, C.; Patricelli, M.P.; Yates, J.R.; Gray, N.S.; Lee, J.-D. Pharmacological inhibition of BMK1 suppresses tumor growth through promyelocytic leukemia protein. Cancer Cell 2010, 18, 258–267. [CrossRef] [PubMed] 72. Poulikakos, P.I.; Zhang, C.; Bollag, G.; Shokat, K.M.; Rosen, N. RAF inhibitors transactivate RAF dimers and ERK signalling in cells with wild-type BRAF. Nature 2010, 464, 427–430. [CrossRef] 73. O’Reilly, K.E.; Rojo, F.; She, Q.-B.; Solit, D.; Mills, G.B.; Smith, D.; Lane, H.; Hofmann, F.; Hicklin, D.J.; Ludwig, D.L.; et al. mTOR inhibition induces upstream receptor tyrosine kinase signaling and activates Akt. Cancer Res. 2006, 66, 1500–1508. [CrossRef] [PubMed] 74. Lochhead, P.A.; Tucker, J.A.; Tatum, N.J.; Wang, J.; Oxley, D.; Kidger, A.M.; Johnson, V.P.; Cassidy, M.A.; Gray, N.S.; Noble, M.E.M.; et al. Paradoxical activation of the protein kinase-transcription factor ERK5 by ERK5 kinase inhibitors. Nat. Commun. 2020, 11, 1383. [CrossRef] [PubMed]