https://orcid.org/0000-0003-3843-4173
Figures
Abstract
Despite the availability of various treatment protocols, response to therapy in patients with Acute Myeloid Leukemia (AML) remains largely unpredictable. Transcriptomic profiling studies have thus far revealed the presence of molecular subtypes of AML that are not accounted for by standard clinical parameters or by routinely used biomarkers. Such molecular subtypes of AML are predicted to vary in response to chemotherapy or targeted therapy. The Renin-Angiotensin System (RAS) is an important group of proteins that play a critical role in regulating blood pressure, vascular resistance and fluid/electrolyte balance. RAS pathway genes are also known to be present locally in tissues such as the bone marrow, where they play an important role in leukemic hematopoiesis. In this study, we asked if the RAS genes could be utilized to predict drug responses in patients with AML. We show that the combined in silico analysis of up to five RAS genes can reliably predict sensitivity to Doxorubicin as well as Etoposide in AML. The same genes could also predict sensitivity to Doxorubicin when tested in vitro. Additionally, gene set enrichment analysis revealed enrichment of TNF-alpha and type-I IFN response genes among sensitive, and TGF-beta and fibronectin related genes in resistant cancer cells. However, this does not seem to reflect an epithelial to mesenchymal transition per se. We also identified that RAS genes can stratify patients with AML into subtypes with distinct prognosis. Together, our results demonstrate that genes present in RAS are biomarkers for drug sensitivity and the prognostication of AML.
Citation: Turk S, Turk C, Akbar MW, Kucukkaraduman B, Isbilen M, Demirkol Canli S, et al. (2020) Renin angiotensin system genes are biomarkers for personalized treatment of acute myeloid leukemia with Doxorubicin as well as etoposide. PLoS ONE 15(11): e0242497. https://doi.org/10.1371/journal.pone.0242497
Editor: Francesco Bertolini, European Institute of Oncology, ITALY
Received: April 3, 2020; Accepted: November 4, 2020; Published: November 25, 2020
Copyright: © 2020 Turk et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the manuscript and its Supporting information files.
Funding: This study was supported by The Scientific and Technological Research Council of Turkey (116S350).
Competing interests: The authors have declared that no competing interests exist.
Introduction
Leukemia, lymphoma and multiple myeloma are the three main types of highly heterogeneous hematological malignancies that are derived from myeloid and lymphoid cell lineages [1]. Acute myeloid leukemia (AML) is characterized by abnormal expansion of immature myeloid cells and their accumulation in the bone marrow and blood, interfering with normal cellular growth [2]. AML is a highly aggressive cancer with poor prognosis. It is also the most common type of acute leukemia in adults. Treatment strategies and success rates vary depending on many factors, including the subtype of AML, prognostic factors, age and general health status of the patient [3]. Standard treatment regimens based on patient stratification include the combination of chemotherapeutics such as Cytarabine, Daunorubicin and Etoposide with or without radiotherapy. However, high heterogeneity of clinical outcomes in AML patients suggests that current classifications fail to distinguish patient subgroups sufficiently [4].
A not so well studied protein network in the context of AML is the Renin-Angiotensin System (RAS). RAS is composed of several gene products which play a critical role in regulating blood pressure, renal vascular resistance and the fluid/electrolyte balance [5, 6]. The idea of a local RAS operating independent of the circulating RAS was brought into light by demonstrating localized RAS elements in organs other than liver (angiotensinogen), kidney (renin) and lung (ACE). Localized RAS elements were found in many organs such as the brain, blood vessels and heart [7, 8]. It is predicted that locally produced angiotensins have important homeostatic functions and may contribute to local tissue dysfunction and diseases [8]. The presence of local RAS specific to the hematopoietic bone marrow microenvironment was reported for the first time in 1996 [9]. Major RAS molecules have been identified in the bone marrow microenvironment, such as renin, angiotensinogen, angiotensin receptors and angiotensin converting enzymes (ACEs) [10]. Locally active bone marrow RAS affects important stages of physiological and pathological blood cell production through autocrine, paracrine and intracrine pathways [11, 12]. Local bone marrow RAS peptides control the development of hematopoietic niche, myelopoiesis, erythropoiesis, thrombopoiesis and other cellular lineages [13–19]. Local RAS is also active in the primitive embryonic hematopoiesis phase [20–23]. The presence of renin, ACE, angiotensin II (Ang-II) and angiotensinogen in leukemic blast cells has been demonstrated, and local bone marrow RAS has been shown to play a role in the development of neoplastic malignant blood cells [24–26].
Establishing a role for genes involved in the development and biology of cancers, as prognostic and chemotherapeutic markers, is one of the most effective and successful approach used in the classification of malignancies. Thus, here we aimed to define AML subgroups based on expression of RAS genes. We also aimed to test if the resulting tumor subtypes differ in their responses to drugs and to demonstrate distinct prognostic profiles.
Materials and methods
In silico
Datasets.
Cancer Genome Project (CGP) gene expression data (E-MTAB-783) [27, 28] was downloaded from ArrayExpress website (https://www.ebi.ac.uk/arrayexpress/), and drug screening data [29] was downloaded from the CGP database. Microarray dataset GSE12417 [30], corresponding to AML patients, was downloaded from the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database. The training cohort of 163 patients in GSE12417 [30] was used for our analyses since this was the same microarray platform as the one used for CGP.
Data normalization and variance analysis.
CGP gene expression data was normalized by the RMA method using the BRB-Array Tools software [31].
In order to choose the RAS genes that would be used in real-time PCR for validation studies, we first aimed to choose the most variable genes that would likely give detectable fold differences in vitro by PCR. Analysis of variance was performed using all 39 probesets corresponding to the 25 genes in the RAS and genes with at least 0.8 of variance. Above these thresholds, the mean expression was at least 5.5 and the log fold difference between min to max was above 3 for all probesets. Thus, nine probesets corresponding to eight genes (CTSG, CPA3, AGT, ANPEP, IGF2R (two probesets), RNPEP, ATP6AP2 and CTSA) were selected to be used in further analyses (S1 Table).
IC50 calculation methods.
In order to calculate drug response parameters such as IC50, EC50, activity area and Amax, the growth rate of the cells were depicted as a function of drug concentration by being modeled with non-linear logistic regression as explained in De Lean et. al [32], which is also reported in NIH/NCGC assay guidelines [33]. While the non-linear logistic regression function used to model data is used widely for cytotoxicity calculations, here for the first time we used six different versions of this function and selected the one with the lowest standard error rate among all for the calculation of cytotoxicity values. We name this approach the 6-model (6M).
Thus, six different models were derived from the following non-linear logistic regression function:
where Y is the percent growth of the cells, X is the arithmetic drug concentration, a is the percent growth of the cells when the cells are not treated with the drug, d is the percent growth of the cells for infinite dose, i.e. a dose for which there is no additional effect when increased, c is the dose corresponding to percent growth exactly between a and d, and b is the Hill slope factor that is used to define the steepness of the curve fitted.
The following are the conditions required for the generation of 6-models:
- 3-Parameter model: Curves were fitted without using Hill slope factor b.
- 3-Parameter Top 100 model: Curves were fitted without using Hill slope factor b and with a = 100.
- 3-Parameter Bottom 0 model: Curves were fitted without using Hill slope factor b and with d = 0.
- 4-Parameter model: Formula is used as it is.
- 4-Parameter Top 100 model: Curves were fitted with a = 100.
- 4-Parameter Bottom 0 model: Curves were fitted with d = 0.
Six different drug response parameters are calculated out of the fitted curves as follows:
- IC50: Value of X when Ŷ = 50%
- IC90: Value of X when Ŷ = 90%
- IC95: Value of X when Ŷ = 95%
- EC50: Value of X when Ŷ = a+d
- Amax: a − d
- Activity Area: ΣŶX, (sum of Ŷs for each 0.01 increment of X fitted), where Ŷ is the predicted value of Y by the curve fitted.
With the 6M approach we recalculated IC50 values that were also included in the raw CGP data for the 17 AML cell lines treated with four drugs (ATRA, Cytarabine, Etoposide and Doxorubicin) using an in-house R script “SixModelIC50 V3.r” (https://github.com/muratisbilen/6-Model_IC50_CalculationV3.git).
These drugs were selected as we obtained AML chemotherapy treatment protocols from the Department of Hematology, Hacettepe University and compiled a list for all drugs in these protocols. Among these only ATRA, Cytarabine, Etoposide and Doxorubicin were present in the CGP database.
We referred to the recalculated IC50 data as 6M IC50 and performed a Pearson r correlation analysis between CGP IC50s and recalculated 6M IC50s to test the compatibility.
In addition, IC50 values were calculated using the 6M approach on the data obtained from in vitro analysis in which nine AML cell lines were treated with Doxorubicin and Etoposide.
Linear regression analyses.
We performed correlation analysis between expression values of the eight genes and drug data (both CGP IC50 data and 6M IC50 data) individually. To identify if multiple genes can be used to better identify the relationship between gene expression and drug sensitivity data, linear regression analyses were performed using the Minitab 17 software (https://www.minitab.com). Seventeen AML cell lines from the CGP database were either randomly divided into two groups, the discovery group (12 cell lines) and the test group (five cell lines), or chosen manually so that the sensitivity range of cells in both groups spanned as large variance as possible. To generate a linear regression model for each drug (ATRA, Cytarabine, Etoposide, Doxorubicin), IC50s of the discovery cell line group obtained either from CGP or recalculated as 6M IC50, and expression of the eight RAS genes which were selected from variance analysis, were used as predictors. As a measure of the response variable variation explained by each linear regression model, we used the adjusted (adj.) R2 values. To test consistency of the linear regression models generated with the eight genes, we replicated the random divison of groups ten times and reported the average of the adjusted R2.
Furthermore, to identify a minimal gene list for the prediction of chemosensitivity, the discovery group was used to fit a model explaining the drug response using “best subsets” function of the software, which runs all possible regression models with one variable, two variables and so on, based on a list of predictors, enabling the user to choose a smaller set of predictors that can explain the response. The subset with the highest R2 (adj.) was selected as the best model. Regression formula of the best models (y = ±a + [n1 × x1] ± [n2 × x2] ± [n3 × x3] ± [n4 × x4]…) were applied for the test group of each drug. In the regression formula y (predicted IC50 values) were calculated where a and n are the constant values, x: gene expression values of the 12 cell lines in the discovery group. Also, the goodness of fit measure Sy.x were computed by Graphpad. Sy.x is a standard deviation of the residuals that here has been used to describe the difference in standard deviations of CGP IC50 and 6M IC50 versus predicted IC50s. It is a goodness-of-fit measure used to show how well our predicted IC50s fit with CGP and 6M IC50 values. All the correlations were calculated with Graphpad software as Pearson’s r and p values.
Hierarchical clustering analysis.
Cluster 3.0 (http://bonsai.hgc.jp/~mdehoon/software/cluster/software.htm) [34] and Java Treeview (http://jtreeview.sourceforge.net/) [35] software were used for hierarchical clustering analysis with mean standardized gene expression values for each dataset. Hierarchical clustering was performed by clustering both genes and arrays using Euclidian distance as similarity metric and complete linkage as clustering method.
Gene sets enrichment analysis—GSEA.
Gene set enrichment analysis was performed using the GSEA guideline (https://www.gsea-msigdb.org/gsea/index.jsp) [36].
Briefly, dataset E-MTAB-783 [27, 28] has 22277 probesets IDs and these were collapsed into 13321 genes. For genes with more than one probeset, one with the highest expression was selected. “C5_all Gene ontology v6.1 database” was used for the analysis which has gene sets that contain genes annotated by the same GO term. We used default filtering criteria in GSEA for gene set sizes, which includes genesets with sizes between 15–500. After applying this filter, analysis was performed for 4081 gene sets.
Mutation analyses.
Mutation data of AML cell lines was downloaded from Genomics of Drug sensitivity in Cancer database (https://www.cancerrxgene.org/downloads/bulk_download) [37]. 14 out of the 17 AML cell lines used in our analyses were available. Seven genes which were mutated in at least three AML cell lines were analyzed further.
In vitro
Cell lines and cytotoxicity experiments.
HEL92.1.7 (2111706), and QIMR-WIL (86030601) cell lines were purchased from Sigma Aldrich (St. Louis, Mo., USA), KASUMI-3 (CRL-2725), GDM-1 (CRL-2627) and CESS (TIB-190) cell lines were purchased from ATCC (Virginia, USA) and P31/FUJ (JCRB0091), NOMO-1 (IFO50474), KASUMI-1 (JCRB1003) and SKM-1 (JCRB0118) cell lines were purchased from JCRB Cell Bank (Osaka, Japan). Cell lines were authenticated by manufacturers, all cell lines were morphologically checked by microscope and routine mycoplasma testing was performed by PCR. HEL92.1.7, GDM-1, CESS, P31/FUJ and NOMO-1 were cultured and maintained in RPMI-1640 medium (Sigma-Aldrich, R0883 (St. Louis, Mo., USA)) supplemented with 10% fetal bovine serum (FBS) (Sigma-Aldrich, F6178 (St. Louis, Mo., USA)), 1% penicillin-streptomycin (Sigma-Aldrich, 11074440001 (St. Louis, Mo., USA)), and 1% 200 mM L-glutamine (Sigma-Aldrich, G7513 (St. Louis, Mo., USA)). KASUMI-1, SKM-1 and KASUMI-3were cultured in RPMI-1640 medium but with 20% FBS. QIMR-WIL was cultured in DMEM medium (Sigma-Aldrich, D6546 (St. Louis, Mo., USA)) but with 10% FBS, 1% penicillin-streptomycin, and 1% 200 mM L-glutamine. All cell lines were cultured at 5% CO2 and 37 °C in a humidified incubator.
Doxorubicin (D1515) and Etoposide (E1383) were purchased from Sigma-Aldrich (St. Louis, Mo., USA) and were dissolved in DMSO. Cell viability was measured using CellTiter-Glo reagent (G7572, Promega, Fitchburg, Wisconsin, USA). 7000 cells/well in 90 μl medium were plated in each well of a 96-well plate. Cells were treated with six different concentrations of Doxorubicin or Etoposide separately (20, 10, 2, 1, 0.2, 0.1 μM). After 72 hours of drug treatment, cells were treated with CellTiter-Glo reagent and the luminescence signal was then recorded with a microplate luminometer (Turner Designs, CA, USA). All drug treatment experiments were repeated three times. Growth percentages were calculated for each drug and cell line, and cytotoxicity values were calculated using the 6M approach.
qRT-PCR.
AGT, ANPEP, ATP6AP2, CPA3, CTSA and IGF2R genes’ expression was quantified using SYBR ™ Green master mix (Bio-Rad, #1725150, (USA)). PCR reactions were run under cycling conditions according to manufacturer’s instructions. GAPDH was used as a reference gene in all reactions. qRT-PCR relative gene expression data was calculated using ddCT method [38].
Using qRT-PCR relative gene expression data, predicted IC50 values were calculated with the formulas generated by linear regression analyses of in silico data using qRT-PCR based expression values as predictors. Primers used in this study are shown in Table 1. GAPDH was used as endogenous control.
Clinical data validation
Log Rank with Multiple Cutoffs (LRMC) and survival analysis.
In regression analysis, four formulas were generated for Doxorubicin and Etoposide using both CGP and 6M IC50 data. IGF2R, CTSA, ATP6AP2 are common in three of the four formulas except for 6M IC50 data for Etoposide. Therefore, these genes were chosen to test relationships with clinical outcome.
Clinical data were obtained from the training cohort of the GSE12417 [30] dataset (AML Cooperating Group 1999). In the AMLCG 1999 cohort, patients were treated with TAD: Thioguanine, Cytarabine and Daunorubicin, or HAM: Cytarabine, Mitoxantrone protocols followed by the TAD protocol. We used an in-house R script (https://github.com/muratisbilen/LRMC.git) (Log Rank Multiple Cutoff, LRMC) by which log-rank test-based p-values associated with hazard ratio (HR) could be obtained using all possible cutoff values representing each sample in a given dataset and best cutoff is selected as in [39, 40]. Using this approach, we selected best cutoffs for IGF2R, ATP6AP2 and CTSA genes to be used for clinical correlation studies and Kaplan-Meier plots.
Patients with gene expression lower than cutoff, for each gene individually, were labeled as ‘Low’ (low expression) and higher than cutoff were labelled as ‘High’ (high expression). Univariate cox regression analyses were performed and Kaplan-Meier graphs were drawn using SPSS Statistics 19 (IBM, Chicago, IL, USA).
Additionally, the expression of all these three genes (IGF2R, CTSA and ATP6AP2) was evaluated together as good and bad prognostic groups. Patients were grouped as “Good” if they have high expression levels of IGF2R and CTSA and low expression levels of ATP6AP2 defined by expression value cutoffs in previous analysis. Rest of the patients were grouped as “Bad”. Then Kaplan Meier survival analysis was performed for these groups.
Results
Discovery of RAS drug sensitivity biomarker genes
The RAS consists of the 25 genes, corresponding to 39 probesets in Affymetrix HG-U133A, a microarray platforms used in the Cancer Genome Project (CGP) [27, 28]. For the 17 AML cell lines, both drug cytotoxicity and gene expression data are available in the CGP database [27, 28]. We focused only on genes which showed high variation in expression for further validation and therefore, selected nine probesets (eight genes) as described in the methods section (S1 Table). We recalculated IC50 values using the 6M approach applied to raw CGP cytotoxicity data (see Materials and methods). Using Pearson correlation we observed strong correlations between CGP IC50 and 6M IC50 for drugs widely used in AML (Cytarabine, Etoposide, Doxorubicin) but not for ATRA (S2 Table). To identify biomarkers of chemosensitivity, we calculated Pearson correlation between gene expression and IC50 values obtained from CGP and generated by 6M approach. We thus identified six gene/drug cytotoxicity correlations which were significant with CGP IC50 values, and seven significant correlations with 6M IC50 values. Four gene/drug associations were common to both analyses (S3 Table). Linear regression analysis was then performed to test whether the combined expression analyses of genes could correlate better with drug sensitivity data or not. Thus, we generated discovery and test groups. Each group include a wide range of cell line IC50 values as possible. Linear regression models for drug sensitivity prediction were generated for the discovery group (12 cell lines) using expression data of highly variant eight genes and IC50 values obtained from CGP and 6M IC50 of four drugs in Minitab 17. Then, obtained results tested with the validation group (five cell lines). The models generated with combined expression analyses of the eight genes resulted in high R2 (adj) values for Etoposide, Doxorubicin and Cytarabine but no model could be generated for ATRA (S4 Table). As independent datasets with drug sensitivity data for these compounds do not exist, we utilized a cross-validation method to test the robustness of the proposed models by generating the discovery and test groups 10 times, with 12 and five cell lines, respectively. The average of 10 R2 values generated from discovery groups was calculated for both CGP and 6M IC50s. Our results showed that the 10 random models of sensitivity to Doxorubicin had an average R2 above 85% for both CGP and 6M IC50s, but R2 decreased slightly for models of sensitivity to Etoposide while R2 values highly decreased for models of sensitivity to Cytarabine (S4 Table) when compared to those generated for cell lines that were manually selected. We therefore, focused on Doxorubicin and Etoposide for further analyses.
We then aimed to identify the minimal number of genes that needed to be included in combinations into the models that would give the highest correlation using the ‘best subsets function’ of Minitab. The software selected three genes/probesets for Doxorubicin when either CGP and 6M IC50 values were used and, four and five genes/probeset combinations for Etoposide using CGP and 6M IC50 values, respectively; all together corresponding to a total of six genes (AGT, ANPEP, ATP6AP2, CPA3, CTSA and IGF2R (two probesets)) (S5 Table), when the analysis was performed with the discovery group. Applying the resulting models to the test group showed the reliability of all models. As shown in Fig 1, the goodness of fit measures (R sq. and Sy.x) were 0.9 and 0.21 for Doxorubicin as modeled using 6M IC50 data and 0.89 and 0.34 when we used CGP IC50 values. Similarly, for Etoposide, these two measures were 0.78 and 0.34 for 6M IC50 and 0.77 and 0.57 for CGP IC50 values.
Linear regression models were generated using the discovery group and applied to the test group to predict sensitivity values. Reliability of sensitivity predictions was measured with goodness of fit test for Doxorubicin 6M IC50 (A) resulting 0.9 R sq. and 0.21 Sy.x Doxorubicin CGP IC50 resulting (B) 0.89 R sq. and 0.34 Sy.x Etoposide 6M IC50 resulting (C) 0.78 R sq. and 0.34 Sy.x Etoposide CGP IC50 resulting (D) 0.77 R sq. and 0.57 Sy.x with 90, 95, and 99% confidence intervals. Black dots represent cell lines used for discovery group, and red dots for the test/validation group.
In vitro validation of biomarker genes
We next asked if the linear regression models generated in silico could predict in vitro cytotoxicity. For this purpose, we determined gene expression values by qRT-PCR for the six RAS genes (AGT, ANPEP, ATP6AP2, CPA3, CTSA and IGF2R (two probesets)) and used these to predict in vitro IC50 values obtained for Etoposide and Doxorubicin calculated with 6M approach for nine AML cell lines (see Materials and methods section). Correlation analysis showed that in silico and in vitro gene expression data were highly concordant except for CTSA (r: >0.7 and p-value <0.05) (S6 Table). We applied normalized gene expression values obtained in vitro to the in silico generated linear regression models (using four regression formulas) (S7 Table). Thus, utilized linear regression formulas with qRT-PCR gene expression data showed a good correlation with in silico data for Doxorubicin but not Etoposide (Table 2).
Biological features of drug sensitive and resistant cells
Cell lines sensitive to Etoposide and Doxorubicin were almost identical (S1 Fig). To determine molecular mechanisms underlying differential response to Etoposide and Doxorubicin, we performed gene set enrichment analyses (GSEA) with sensitive and resistant subgroups for Gene Ontology (GO) gene sets. Several gene sets were significantly enriched among sensitive and resistant cell lines (FDR q-value<0.25). Gene sets enriched in sensitive cells with a FDR q-value of lower than 0.25 included TNF-receptor interacting, and response to type I IFN stimulus; while gene sets such as regulation of TGF-beta production and FN-binding were enriched in resistant cells suggesting a mesenchymal phenotype (S2 Fig and S8 Table).
To determine if the differentially expressed genes could be reflecting Epithelial-Mesenchymal Transition (EMT), we compared E-Cadherin (epithelial marker) and Vimentin (mesenchymal marker) expression using t-test, between sensitive and resistant cell groups. EMT is the process that epithelial cells lose the apical-basal polarity and cell adhesion, and transform to invasive mesenchymal cells [41]. It is known to play an important role in biological and pathological processes such as cancer progression, metastasis and drug resistance [42–44]. In our analysis, E-cadherin and Vimentin expression were not significantly different between sensitive and resistant groups defined in S1 Fig (p>0.1) (S3 Fig).
Then, in order to understand whether the mutational profile is involved in sensitivity to Doxorubicin, we analyzed mutational data of sensitive, intermediate and resistant groups of AML cell lines (see Materials and methods). Although we have a small sample size, especially in the resistant group (n = 2), we observed that both of the resistant cell lines are NRAS and P53 mutant, whereas all of the sensitive cell lines (n = 5) were wild type for these genes. However, these results need to be validated in larger sample sizes to be conclusive (S9 Table).
RAS genes are prognostic biomarkers for AML
We then asked if the RAS gene expression could help prognosticate AML patients. For this purpose, we utilized the training set within the GSE12417 [30] dataset, including 163 samples of bone marrow or peripheral blood mononuclear cells from adult patients with untreated acute myeloid leukemia. Patients in this cohort were also-treated with TAD protocol which contains Daunorubicin, which is also used as the starting material for semi-synthetic manufacturing of Doxorubicin. We found that high expression of genes IGF2R and CTSA were both associated with better overall survival, while the opposite was true for ATP6AP2 when patients were classified in either “High” or “Low” groups based upon LRMC cutoffs for each gene separately (see Materials and methods section) (Fig 2). We then stratified patients into “Good” and “Bad” prognosis groups using the best cutoff values obtained for these three genes as explained in the method section. As shown in Fig 3, it was revealed that there was a striking difference in overall survival in the groups that were predicted as "Good" and "Bad". The "Good" group showed better survival than the "Bad" group. Since the patients were all treated with Daunorubicin, these data suggest that the expression pattern of these genes was able to identify patients which are responders of this therapy.
(A) LRMCs of IGF2R (Probeset: 201392_s_at), CTSA (Probeset: 200661_at) and ATP6AP2 (Probeset: 201444_s_at). Graphic shows log rank based p values in the y axis for the “high” and “low” expression groups generated by all possible expression based cutoffs shown on the x axis (for details see Materials and methods). HRs above one and below one are shown with red and blue colors for specific cutoffs, respectively. Vertical dotted lines show 25th, 50th and 75th percentiles and horizontal dotted line shows significance cutoff 0.05 (-log10(p) = 1.301). From LRMC graphs, we selected cutoffs 7.077, 11.247 and 11.773 for IGF2R, CTSA and ATP6AP2 respectively which are highlighted in figure with red circle. Patients were divided into high and low groups based on these cutoffs. (B) Kaplan Meier plots for patients classified in high and low expression based on LRMC cutoffs. Patients classifed in high expression group of IGF2R and CTSA showed better overall survival when compared with low expression group and high expression group of ATP6AP2 showed worse survival when compared with low expression group. For all survival plots, overall survival time is shown in days for 163 AML patients. Table at the bottom shows number of patients in each group, median survival for each group and Log rank p value for Kaplan Meier analysis.
Patients were grouped as “Good” if they have high expression levels of IGF2R and CTSA and low expression levels of ATP6AP2 defined by expression value cutoffs in Fig 2. Rest of the patients were grouped as “Bad”. Kaplan Meier plot shows “Good” group showed better survival when compared with “Bad” group as expected. Table at the bottom shows number of patients in each group, median survival for each group and Log rank p value for Kaplan Meier analysis.
Discussion
RAS’ local presence in the marrow affects the most important stages of physiological and pathological blood cell proliferation, and also has important roles in the development of blood cancers. It has been shown that RAS plays important roles in drug resistance to chemotherapeutic agent in addition to angiogenesis, invasion and proliferation [9, 24, 45–47]. Inevitably, most of these processes are interdependent. Most of the increased metastasis and invasion occurs due to an active RAS results in angiogenesis [45, 48, 49]. AT1R upregulation in ovarian cancer and increased expression of AT1R and ACE in prostate cancer, and AGTR1 in breast cancer; localized RAS presence in gastric cancer and its correlation with tumor spread and progression; demonstrate strong associations of RAS with various cancers. Irregularity of RAS components in cancer is strongly associated with increased angiogenesis and metastasis, and these parameters are associated with poor prognosis [50–54].
Gene expression profiling has revealed various AML subtypes related to diagnosis, therapy response and prognosis [55, 56]. Although gene expression profiling has not yet been integrated into clinical practice, this is expected to happen in near future.
In our study, we focused on RAS genes and identified their association with Doxorubicin and Etoposide sensitivity. We also show that RAS genes can be used to stratify AML patients into groups with distinct prognoses. Similar to our findings, low expression of IGF2R in non-small cell liver cancer has been associated with poor prognosis and high expression in bladder cancer has been associated with good prognosis [57, 58]. Although high CTSA expression was associated with poorer outcome in breast ductal carcinoma in situ, it was also found to suppress invasion and metastasis of colorectal cancer, suggesting tissue-specific differential roles [59, 60]. Recent studies linked ATP6AP2 up-regulation to the progression of glioma and colorectal cancer, due to its roles in aberrant activation of the Wnt/beta-catenin signalling pathway [61, 62]. ATP6AP2 was also shown to be a key component of the pro-angiogenic/proliferative arm of the RAS, which plays a role in the growth and spread of endometrial cancer [63]. Compared to the presence in the lysosome, it is found more in the cell membrane. Thus, it is clear that in this way it induces TGF-beta pathway activation. IGF2R is located in the membrane of organelles and is responsible for the transport to lysosome, and its intracellular functions have not yet been clearly identified. CTSA is a protease found in the lysosome. The fact that these three genes function together in the lysosome suggests that lysosomal functions can contribute to cell sensitivity. ATP6AP2 gene was found to cause disruption of V-ATPase formation and defects in the lysosomal glycosylation and autophagy [64]. Supportively, Doxorubicin has been reported to cause autophagy induced cell death in AML cells [65, 66].
GSEA revealed that sensitive cells were correlated with TNF-receptor interacting and response to type I IFN gene sets and resistant cells were correlated with regulation of TGF-beta production and FN-binding gene sets in AML, suggesting a mesenchymal phenotype.
A good and reliable subgrouping which can predict Doxorubicin sensitivity in AML was performed with the ATP6AP2, IGF2R, and CTSA gene combination. For those analyses, we utilized a Daunorubicin treated cohort, which is used as the starting material for semi-synthetic manufacturing of Doxorubicin. Therefore, the combination of these genes which can predict the sensitivity of Doxorubicin in AML patients may, therefore, be confirmed ex vivo.
The mutational analyses performed in this study had a small sample size with only two resistant cells. Therefore more conclusive results would be reached when this type of analysis is performed with larger sample sizes, or when mutational profiling is performed in patients treated with Doxorubicin, which may shed light on future studies.
Conclusions
As a result, we identified IGF2R, CTSA and ATP6AP2 gene biomarkers, which can subgroup AML patients into distinct good and bad prognostic groups. ATP6AP2 was associated to resistance and IGF2R and CTSA were associated to sensitivity for Doxorubicin in silico and in vitro. In future studies, it is important to investigate whether these genes can be used for personalized treatment and improve the effectiveness of treatments.
Supporting information
S1 Fig. Hierarchical clustering of AML cell lines by sensitivity profiles for Doxorubicin and Etoposide.
The analysis reveals sensitive (six cell lines-green), intermediate (eight cell lines-orange) and resistant (three cell lines-red) subgroups for the 17 AML cell lines. Sensitivity to Doxorubicin and Etoposide is highly concordant in three subgroups. Green indicate low expression, orange indicate intermediate expression and red indicates high expression.
https://doi.org/10.1371/journal.pone.0242497.s001
(TIF)
S2 Fig. Comparative analysis of differentially enriched gene sets among drug sensitive and resistant cell lines.
(A) Plots showing gene sets enriched in sensitive cells, including genes interacting with TNF-receptor and genes affected in response to type I IFN stimulus. (B) Plots showing gene sets enriched in resistant cell lines, including genes having role in regulation of TGF-B production and genes interacting selectively and non-covalently with Fibronectin.
https://doi.org/10.1371/journal.pone.0242497.s002
(TIF)
S3 Fig. Expression levels of E-cadherin and Vimentin genes in AML cell lines.
RMA normalized gene expression values of CGP microarray data (y-axis) were used to determine EMT status of sensitive and resistant AML cell lines (x-axis) defined in S1 Fig. VIM: Vimentin (black bars), CDH1: E-cadherin (white bars). n.s. (not significant).
https://doi.org/10.1371/journal.pone.0242497.s003
(TIF)
S1 Table. Expression variance of RAS genes in AML cell lines.
RMA normalized gene expression values of 25 RAS genes were used to analyze variance, standard deviation (SD), mean and min-max difference among cell lines. Gene names is shown along with probe set ID.
https://doi.org/10.1371/journal.pone.0242497.s004
(XLSX)
S2 Table. Pearson correlation analysis between CGP IC50s and 6M IC50s.
IC50 values recalculated according to 6M approach using CGP raw cytotoxicity measurements were used to calculate Pearson correlation analysis with CGP IC50 values. Strong correlations are observed for all drugs except for ATRA.
https://doi.org/10.1371/journal.pone.0242497.s005
(PDF)
S3 Table. Pearson correlation analysis between CGP gene expression data and CGP / 6M IC50s for 17 AML cell lines.
Eight genes expression data (nine probesets) was used to calculate correlation (as R2 coefficient of determination) with IC50 values of four drugs (ATRA, Cytarabine, Etoposide, Doxorubicin) from CGP database along with recalculated data with 6M approach. Highlighted values with green and red indicate significant correlation in negative and positive manner respectively.
https://doi.org/10.1371/journal.pone.0242497.s006
(PDF)
S4 Table. Linear regression analysis between expression values of nine probesets and CGP / 6M IC50s for the discovery group (12 cell lines) and also for the ten times randomly divided different discovery groups (12 cell lines).
Adjusted R2 values were calculated in Minitab 17 with eight genes (nine probesets) for four drugs using CGP gene expression data and CGP / 6M IC50 values of the 12 AML cell lines. High correlations are observed with Etoposide and Doxorubicin (bold, for Etoposide R2 > 90%, for Doxorubicin R2 > 80%). Averages of adjusted R2 values of ten randomly divided groups were calculated in Minitab 17 with eight genes (nine probesets) for four drugs using CGP gene expression data and CGP / 6M IC50 values of the 12 AML cell lines. High correlation is observed with Doxorubicin and fine with Etoposide (bold, for Doxorubicin R2 > 85%, for Etoposide R2 > 60%). Asterisk represents the analysis in which Minitab could not perform linear regression analysis.
https://doi.org/10.1371/journal.pone.0242497.s007
(PDF)
S5 Table. Generation of linear regression models using CGP gene expression data and CGP / 6M IC50 data of the discovery group (12 AML cell lines) for drug sensitivity predictions.
(A) Individual genes and gene combinations were used to generate linear regression models using IC50 values of Doxorubicin and Etoposide from CGP and 6M IC50. Highest correlation is observed in IGF2R/ATP6AP2/CTSA combination with Doxorubicin CGP and 6M IC50 values. And, highest correlation is observed in IGF2R/ATP6AP2/CTSA/CPA3 combination with Etoposide CGP and in ANPEP/ATP6AP2/CTSA/CPA3/AGT combination with Etoposide 6M IC50 values. (B) Regression formulas for gene panels with highest correlations.
https://doi.org/10.1371/journal.pone.0242497.s008
(PDF)
S6 Table. Relative expression values of ATP6AP2, IGF2R (two probesets), CTSA, CPA3, AGT and ANPEP genes in nine AML cell lines and Pearson’s correlation analysis for six genes between CGP gene expression data and in vitro qRT-PCR gene expression data of nine AML cell lines.
(A) Expressions of all genes was normalized to GAPDH expression. (B) ATP6AP2, IGF2R (two probesets), CPA3, AGT, and ANPEP gene expression data obtained from CGP in silico and in vitro qRT-PCR expression data from nine cell lines show significant correlations with in vitro qRT-PCR expression data with the exception of CTSA.
https://doi.org/10.1371/journal.pone.0242497.s009
(PDF)
S7 Table. In vitro and predicted IC50s (from CGP and 6M IC50 linear regression formulas) of Doxorubicin and Etoposide for nine AML cell lines.
In vitro IC50 values were obtained from cell viability measurements of the cell lines that are treated with six different concentrations of Doxorubicin and Etoposide separately (20, 10, 2, 1, 0.2, 0.1 μM). Predicted IC50s were calculated using the four formulas generated (with CGP / 6M IC50s) in the linear regression analysis with the normalized gene expression data obtained from qRT-PCR.
https://doi.org/10.1371/journal.pone.0242497.s010
(PDF)
S8 Table. Gene sets enriched in (A) sensitive cell lines, and (B) resistant cell lines.
https://doi.org/10.1371/journal.pone.0242497.s011
(PDF)
S9 Table. Mutational data of sensitive, intermediate and resistant groups of AML cell lines.
Seven genes which are mutated in at least three AML cell lines were included to examine the relationship between mutational status and drug sensitivity. Red: mutations that cause change in aminoacid sequence, grey: unkown status of aminoacid change, change at the DNA level; blue: wild type.
https://doi.org/10.1371/journal.pone.0242497.s012
(PDF)
References
- 1.
Stein GS, Luebbers KP. Cancer: Prevention, Early Detection, Treatment and Recovery: John Wiley & Sons; 2019.
- 2. Chen SJ, Shen Y, Chen Z. A panoramic view of acute myeloid leukemia. Nat Genet. 2013;45(6):586–7. pmid:23715324
- 3. De Kouchkovsky I, Abdul-Hay M. ‘Acute myeloid leukemia: a comprehensive review and 2016 update’. Blood Cancer J. 2016;6(7):e441. pmid:27367478
- 4.
Chu E. Physicians’ cancer chemotherapy drug manual 2018. [S.l.]: Jones & Bartlett Learning; 2017.
- 5. Chiong JR, Aronow WS, Khan IA, Nair CK, Vijayaraghavan K, Dart RA, et al. Secondary hypertension: current diagnosis and treatment. Int J Cardiol. 2008;124(1):6–21. pmid:17462751
- 6. Bader M. Tissue renin-angiotensin-aldosterone systems: Targets for pharmacological therapy. Annu Rev Pharmacol Toxicol. 2010;50:439–65. pmid:20055710
- 7. Kuba K, Imai Y, Penninger JM. Angiotensin-converting enzyme 2 in lung diseases. Curr Opin Pharmacol. 2006;6(3):271–6. pmid:16581295
- 8. Campbell DJ. Clinical relevance of local Renin Angiotensin systems. Front Endocrinol (Lausanne). 2014;5:113. pmid:25071727
- 9. Haznedaroglu IC, Tuncer S, Gursoy M. A local renin-angiotensin system in the bone marrow. Med Hypotheses. 1996;46(6):507–10. pmid:8803932
- 10. Strawn WB, Richmond RS, Ann Tallant E, Gallagher PE, Ferrario CM. Renin-angiotensin system expression in rat bone marrow haematopoietic and stromal cells. Br J Haematol. 2004;126(1):120–6. pmid:15198742
- 11. Haznedaroglu IC, Ozturk MA. Towards the understanding of the local hematopoietic bone marrow renin-angiotensin system. Int J Biochem Cell Biol. 2003;35(6):867–80. pmid:12676173
- 12. Haznedaroglu IC, Beyazit Y. Local bone marrow renin-angiotensin system in primitive, definitive and neoplastic haematopoiesis. Clin Sci (Lond). 2013;124(5):307–23. pmid:23157407
- 13. Hara M O K, Wada H, Sasayama S, Matsumori A. Preformed angiotensin II is present in human mast cells. Cardiovasc Drugs Ther 2004;18. pmid:15770428
- 14. Hubert C, Savary K, Gasc JM, Corvol P. The hematopoietic system: a new niche for the renin-angiotensin system. Nat Clin Pract Cardiovasc Med. 2006;3(2):80–5. pmid:16446776
- 15. Kato H, Ishida J, Imagawa S, Saito T, Suzuki N, Matsuoka T, et al. Enhanced erythropoiesis mediated by activation of the renin-angiotensin system via angiotensin II type 1a receptor. FASEB J. 2005;19(14):2023–5. pmid:16223784
- 16. Kwiatkowski BA R R. Angiotensin II Receptor-Associated Protein (AGTRAP) Synergizes with Mpl Signaling to Promote Survival and to Increase Proliferation Rate of Hematopoietic Cells. ASH Annual Meeting Abstracts 2009;114:3606.
- 17. Lin C, Datta V, Okwan-Duodu D, Chen X, Fuchs S, Alsabeh R, et al. Angiotensin-converting enzyme is required for normal myelopoiesis. FASEB J. 2011;25(4):1145–55. pmid:21148418
- 18. Park TS, Zambidis ET. A role for the renin-angiotensin system in hematopoiesis. Haematologica. 2009;94(6):745–7. pmid:19483149
- 19. Shen XZ, Bernstein KE. The peptide network regulated by angiotensin converting enzyme (ACE) in hematopoiesis. Cell Cycle. 2011;10(9):1363–9. pmid:21441775
- 20. Jokubaitis VJ, Sinka L, Driessen R, Whitty G, Haylock DN, Bertoncello I, et al. Angiotensin-converting enzyme (CD143) marks hematopoietic stem cells in human embryonic, fetal, and adult hematopoietic tissues. Blood. 2008;111(8):4055–63. pmid:17993616
- 21. Sinka L, Biasch K, Khazaal I, Peault B, Tavian M. Angiotensin-converting enzyme (CD143) specifies emerging lympho-hematopoietic progenitors in the human embryo. Blood. 2012;119(16):3712–23. pmid:22282502
- 22. Tavian M, Biasch K, Sinka L, Vallet J, Peault B. Embryonic origin of human hematopoiesis. Int J Dev Biol. 2010;54(6–7):1061–5. pmid:20711983
- 23. Zambidis ET, Park TS, Yu W, Tam A, Levine M, Yuan X, et al. Expression of angiotensin-converting enzyme (CD143) identifies and regulates primitive hemangioblasts derived from human pluripotent stem cells. Blood. 2008;112(9):3601–14. pmid:18728246
- 24. Aksu S, Beyazit Y, Haznedaroglu IC, Canpinar H, Kekilli M, Uner A, et al. Over-expression of angiotensin-converting enzyme (CD 143) on leukemic blasts as a clue for the activated local bone marrow RAS in AML. Leuk Lymphoma. 2006;47(5):891–6. pmid:16753875
- 25. Beyazit Y, Aksu S, Haznedaroglu IC, Kekilli M, Misirlioglu M, Tuncer S, et al. Overexpression of the local bone marrow renin-angiotensin system in acute myeloid leukemia. J Natl Med Assoc. 2007;99(1):57–63. pmid:17304969
- 26. Wulf GG, Jahns-Streubel G, Nobiling R, Strutz F, Hemmerlein B, Hiddemann W, et al. Renin in acute myeloid leukaemia blasts. Br J Haematol. 1998;100(2):335–7. pmid:9488623
- 27. Garnett MJ, Edelman EJ, Heidorn SJ, Greenman CD, Dastur A, Lau KW, et al. Systematic identification of genomic markers of drug sensitivity in cancer cells. Nature. 2012;483(7391):570–5. pmid:22460902
- 28. Venkova L, Aliper A, Suntsova M, Kholodenko R, Shepelin D, Borisov N, et al. Combinatorial high-throughput experimental and bioinformatic approach identifies molecular pathways linked with the sensitivity to anticancer target drugs. Oncotarget. 2015;6(29):27227–38. pmid:26317900
- 29. Iorio F, Knijnenburg TA, Vis DJ, Bignell GR, Menden MP, Schubert M, et al. A Landscape of Pharmacogenomic Interactions in Cancer. Cell. 2016;166(3):740–54. pmid:27397505
- 30. Metzeler KH, Hummel M, Bloomfield CD, Spiekermann K, Braess J, Sauerland M-C, et al. An 86-probe-set gene-expression signature predicts survival in cytogenetically normal acute myeloid leukemia. Blood. 2008;112(10):4193–201. pmid:18716133
- 31.
Simon RM. Design and analysis of DNA microarray investigations. New York; London: Springer; 2011.
- 32. DeLean A, Munson PJ, Rodbard D. Simultaneous analysis of families of sigmoidal curves: application to bioassay, radioligand assay, and physiological dose-response curves. The American journal of physiology. 1978;235(2):E97–102. pmid:686171
- 33.
Beck B, Chen YF, Dere W, Devanarayan V, Eastwood BJ, Farmen MW, et al. Assay Operations for SAR Support. In: Sittampalam GS, Grossman A, Brimacombe K, Arkin M, Auld D, Austin CP, et al., editors. Assay Guidance Manual. Bethesda (MD): Eli Lilly & Company and the National Center for Advancing Translational Sciences; 2004.
- 34. de Hoon MJ, Imoto S, Nolan J, Miyano S. Open source clustering software. Bioinformatics (Oxford, England). 2004;20(9):1453–4. pmid:14871861
- 35. Saldanha AJ. Java Treeview—extensible visualization of microarray data. Bioinformatics (Oxford, England). 2004;20(17):3246–8. pmid:15180930
- 36. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. 2005;102(43):15545–50.
- 37. Yang W, Soares J, Greninger P, Edelman EJ, Lightfoot H, Forbes S, et al. Genomics of Drug Sensitivity in Cancer (GDSC): a resource for therapeutic biomarker discovery in cancer cells. Nucleic acids research. 2013;41(Database issue):D955–61. pmid:23180760
- 38. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25(4):402–8. pmid:11846609
- 39. Akbar MW, Isbilen M, Belder N, Canli SD, Kucukkaraduman B, Turk C, et al. A Stemness and EMT Based Gene Expression Signature Identifies Phenotypic Plasticity and is A Predictive but Not Prognostic Biomarker for Breast Cancer. Journal of Cancer. 2020;11(4):949–61. pmid:31949498
- 40. Demirkol S, Gomceli I, Isbilen M, Dayanc BE, Tez M, Bostanci EB, et al. A Combined ULBP2 and SEMA5A Expression Signature as a Prognostic and Predictive Biomarker for Colon Cancer. Journal of Cancer. 2017;8(7):1113–22. pmid:28607584
- 41. Du B, Shim JS. Targeting Epithelial-Mesenchymal Transition (EMT) to Overcome Drug Resistance in Cancer. Molecules (Basel, Switzerland). 2016;21(7). pmid:27455225
- 42. Hay ED. An overview of epithelio-mesenchymal transformation. Acta anatomica. 1995;154(1):8–20. pmid:8714286
- 43. Kalluri R, Weinberg RA. The basics of epithelial-mesenchymal transition. The Journal of clinical investigation. 2009;119(6):1420–8. pmid:19487818
- 44. Thiery JP. Epithelial-mesenchymal transitions in development and pathologies. Current opinion in cell biology. 2003;15(6):740–6. pmid:14644200
- 45. Ager EI, Neo J, Christophi C. The renin-angiotensin system and malignancy. Carcinogenesis. 2008;29(9):1675–84. pmid:18632755
- 46. Teresa Gomez Casares M, de la Iglesia S, Perera M, Lemes A, Campo C, Gonzalez San Miguel JD, et al. Renin expression in hematological malignancies and its role in the regulation of hematopoiesis. Leuk Lymphoma. 2002;43(12):2377–81. pmid:12613527
- 47. Tawinwung S, Ninsontia C, Chanvorachote P. Angiotensin II Increases Cancer Stem Cell-like Phenotype in Lung Cancer Cells. Anticancer Res. 2015;35(9):4789–97. pmid:26254369
- 48. Egami K, Murohara T, Shimada T, Sasaki K, Shintani S, Sugaya T, et al. Role of host angiotensin II type 1 receptor in tumor angiogenesis and growth. The Journal of clinical investigation. 2003;112(1):67–75. pmid:12840060
- 49. Fujita M, Hayashi I, Yamashina S, Fukamizu A, Itoman M, Majima M. Angiotensin type 1a receptor signaling-dependent induction of vascular endothelial growth factor in stroma is relevant to tumor-associated angiogenesis and tumor growth. Carcinogenesis. 2005;26(2):271–9. pmid:15637093
- 50. Ino K, Shibata K, Kajiyama H, Yamamoto E, Nagasaka T, Nawa A, et al. Angiotensin II type 1 receptor expression in ovarian cancer and its correlation with tumour angiogenesis and patient survival. Br J Cancer. 2006;94(4):552–60. pmid:16434990
- 51. Kinoshita J, Fushida S, Harada S, Yagi Y, Fujita H, Kinami S, et al. Local angiotensin II-generation in human gastric cancer: correlation with tumor progression through the activation of ERK1/2, NF-kappaB and survivin. Int J Oncol. 2009;34(6):1573–82. pmid:19424575
- 52. Rhodes DR, Ateeq B, Cao Q, Tomlins SA, Mehra R, Laxman B, et al. AGTR1 overexpression defines a subset of breast cancer and confers sensitivity to losartan, an AGTR1 antagonist. Proc Natl Acad Sci U S A. 2009;106(25):10284–9. pmid:19487683
- 53. Rocken C, Lendeckel U, Dierkes J, Westphal S, Carl-McGrath S, Peters B, et al. The number of lymph node metastases in gastric cancer correlates with the angiotensin I-converting enzyme gene insertion/deletion polymorphism. Clin Cancer Res. 2005;11(7):2526–30. pmid:15814629
- 54. Uemura H, Hasumi H, Ishiguro H, Teranishi J, Miyoshi Y, Kubota Y. Renin-angiotensin system is an important factor in hormone refractory prostate cancer. Prostate. 2006;66(8):822–30. pmid:16482568
- 55. Marcucci G, Haferlach T, Dohner H. Molecular genetics of adult acute myeloid leukemia: prognostic and therapeutic implications. J Clin Oncol. 2011;29(5):475–86. pmid:21220609
- 56. Haferlach T, Schmidts I. The power and potential of integrated diagnostics in acute myeloid leukaemia. Br J Haematol. 2020;188(1):36–48. pmid:31808952
- 57. Tian Z, Yao G, Song H, Zhou Y, Geng J. IGF2R expression is associated with the chemotherapy response and prognosis of patients with advanced NSCLC. Cell Physiol Biochem. 2014;34(5):1578–88. pmid:25402559
- 58. Wu H, Li Y, Cui H, Wu W, Yang H, Yang X, et al. IGF2/IGF2R expression in urothelial bladder cancer and its implications for tumor recurrence and prognosis. International Journal of Clinical and Experimental Medicine. 2017;10:881–8.
- 59. Ni S, Weng W, Xu M, Wang Q, Tan C, Sun H, et al. miR-106b-5p inhibits the invasion and metastasis of colorectal cancer by targeting CTSA. Onco Targets Ther. 2018;11:3835–45. pmid:30013364
- 60. Toss MS, Miligy IM, Haj-Ahmad R, Gorringe KL, AlKawaz A, Mittal K, et al. The prognostic significance of lysosomal protective protein (cathepsin A) in breast ductal carcinoma in situ. Histopathology. 2019;74(7):1025–35. pmid:30725481
- 61. Kouchi M, Shibayama Y, Ogawa D, Miyake K, Nishiyama A, Tamiya T. (Pro)renin receptor is crucial for glioma development via the Wnt/beta-catenin signaling pathway. J Neurosurg. 2017;127(4):819–28. pmid:28059652
- 62. Wang J, Shibayama Y, Zhang A, Ohsaki H, Asano E, Suzuki Y, et al. (Pro)renin receptor promotes colorectal cancer through the Wnt/beta-catenin signalling pathway despite constitutive pathway component mutations. Br J Cancer. 2019;120(2):229–37. pmid:30555158
- 63. Delforce SJ, Lumbers ER, Corbisier de Meaultsart C, Wang Y, Proietto A, Otton G, et al. Expression of renin-angiotensin system (RAS) components in endometrial cancer. Endocr Connect. 2017;6(1):9–19. pmid:27956412
- 64. Rujano MA, Cannata Serio M, Panasyuk G, Peanne R, Reunert J, Rymen D, et al. Mutations in the X-linked ATP6AP2 cause a glycosylation disorder with autophagic defects. J Exp Med. 2017;214(12):3707–29. pmid:29127204
- 65. Goussetis DJ, Altman JK, Glaser H, McNeer JL, Tallman MS, Platanias LC. Autophagy is a critical mechanism for the induction of the antileukemic effects of arsenic trioxide. J Biol Chem. 2010;285(39):29989–97. pmid:20656687
- 66. Ristic B, Bosnjak M, Arsikin K, Mircic A, Suzin-Zivkovic V, Bogdanovic A, et al. Idarubicin induces mTOR-dependent cytotoxic autophagy in leukemic cells. Exp Cell Res. 2014;326(1):90–102. pmid:24907655