Next Article in Journal
Aryl-4,5-dihydro-1H-pyrazole-1-carboxamide Derivatives Bearing a Sulfonamide Moiety Show Single-digit Nanomolar-to-Subnanomolar Inhibition Constants against the Tumor-associated Human Carbonic Anhydrases IX and XII
Next Article in Special Issue
Skeletal Dysplasias Caused by Sulfation Defects
Previous Article in Journal
The Potential for Natural Antioxidant Supplementation in the Early Stages of Neurodegenerative Disorders
Previous Article in Special Issue
Versican G1 Fragment Establishes a Strongly Stabilized Interaction with Hyaluronan-Rich Expanding Matrix during Oocyte Maturation
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 21
Issue 7
10.3390/ijms21072619
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Significance of MicroRNAs Expression in Regulation of Extracellular Matrix and Other Drug Resistant Genes in Drug Resistant Ovarian Cancer Cell Lines

by
Dominika Kazmierczak
,
Karol Jopek
,
Karolina Sterzynska
,
Barbara Ginter-Matuszewska
,
Michal Nowicki
,
Marcin Rucinski
and
Radoslaw Januchowski
*
Department of Histology and Embryology, Poznan University of Medical Sciences, PL-61-781 Poznan, Poland
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(7), 2619; https://doi.org/10.3390/ijms21072619
Submission received: 3 March 2020 / Revised: 2 April 2020 / Accepted: 7 April 2020 / Published: 9 April 2020
(This article belongs to the Special Issue Extracellular Matrix in Development and Disease 2.0)
Browse Figures
Versions Notes

Abstract

:
Ovarian cancer rates the highest mortality among all gynecological malignancies. The main reason for high mortality is the development of drug resistance. It can be related to increased expression of drug transporters and increased expression of extracellular matrix (ECM) proteins. Our foremost aim was to exhibit alterations in the miRNA expression levels in cisplatin (CIS), paclitaxel (PAC), doxorubicin (DOX), and topotecan (TOP)—resistant variants of the W1 sensitive ovarian cancer cell line—using miRNA microarray. The second goal was to identify miRNAs responsible for the regulation of drug-resistant genes. According to our observation, alterations in the expression of 40 miRNAs were present. We could observe that, in at least one drug-resistant cell line, the expression of 21 miRNAs was upregulated and that of 19 miRNAs was downregulated. We identified target genes for 22 miRNAs. Target analysis showed that miRNA regulates key genes responsible for drug resistance. Among others, we observed regulation of the ATP-binding cassette subfamily B member 1 gene (ABCB1) in the paclitaxel-resistant cell line by miR-363 and regulation of the collagen type III alpha 1 chain gene (COL3A1) in the topotekan-resistant cell line by miR-29a.

1. Introduction

Ovarian cancer poses a growing threat to women’s lives and health and occupies the first place in mortality among gynecological cancers [1]. Most cases of ovarian cancer are diagnosed at a very advanced stage (III or IV according to the International Federation of Gynaecology and Obstetrics (FIGO)) [2,3]. The molecular basis of ovarian cancer is very complicated, as ovarian cancer is very heterogeneous with different histological types. Even in the same histological type, a different profile of cancer-related genes among individual cases can be observed [4]. Despite the dynamic development of gynecologic oncology, the main method of ovarian cancer treatment is the surgical treatment followed by chemotherapy [5]. The first line treatment is based on combination platinum derivatives (carboplatin or cisplatin) and taxanes (paclitaxel—PAC) [2]. The basis for choosing the second line chemotherapy regimen is sensitivity to treatment with platinum derivatives, which determines the prognosis [2]. Low chemotherapy effectiveness result from primary or developed during treatment drug resistance [6]. In over 80% of patients, it comes to recurrence cancer [2]. There are two main types of mechanisms of drug resistance: cellular specific and tissue specific [7]. A cellular mechanism is based on cytostatics active removal from cancer cells. Proteins from the ATP-binding cassette (ABC) family contribute to this phenomenon [8]. Cellular mechanisms also include repairing damaged DNA, developing point mutations in the genes that encode proteins that bind cytostatic drugs, and increasing the activity of anti-apoptotic or pro-survival pathways as well as disrupting apoptotic signaling pathways [5]. On the other hand, the tissue specific mechanisms are associated with tumor vascularization, cell density in the tumor, and expression of extracellular matrix proteins (ECM) [9]. All of these change drug distribution in the tumor tissue, decreasing their availability to tumor cells. Additionally, ECM proteins can induce cell adhesion mediated drug resistance (CAM-DR), leading to increased resistance to apoptosis [10]. Although many genes associated with the development of resistance to chemotherapy are known, the mechanisms of their regulation are still poorly understood. One of the ways to regulate gene expression is regulation at the mRNA level by small noncoding RNA particles, designated as a micro RNA (miRNA). Short, non-coding RNAs were discovered and described for the first time in 1993 by Victor Ambros from Harvard University [11]. miRNAs are short (19–29 nucleotides), single-stranded, non-coding RNA molecules with a phosphate residue (5′ end) and a hydroxyl group (3′ end) [12]. They play an important regulatory role in the expression of genes in both animals and plants [12]. In the case of miRNA, the regulation of gene expression takes place at the post-transcriptional level and involves interaction with mRNA [13]. The seed sequence (6–8 nucleotides) located in the 5′ region of the miRNA plays a key role in this process. miRNA seed recognizes the sequence of the target mRNA complementary within the 3′ UTR (untranslated region). Depending on the degree of complementarity of the mRNA to the mature miRNA, either the transcript degradation (complete or almost complete complementarity) occurs or the translation is inhibited (incomplete complementarity) [12]. The effect of miRNA expression can be pleiotropic. Modulation of one miRNA expression can affect many different genes responsible for the various mechanisms of cancer [14]. The sequences encoding miRNAs in the human genome can be located both in introns of genes and in non-endogenous areas. Some of the miRNA encoding genes may also be located within exons of other genes. miRNAs are classified in groups called miRNA families, and most of them are transcribed by RNA polymerase II [15]. Family members are characterized by common origin, evolutionary conservation in sequence, similar mature miRNA/miRNA* (shared functional characteristics or biological function), and conserved mature miRNA–seed–target relationships [16,17]. It is very interesting that miRNA genes in the same miRNA family are non-randomly localized around genes connected with cancer, infectious, immune system, sensory system, neurodegenerative diseases, and development [16]. A set of two or more [18] miRNAs transcribed from physically adjacent miRNA genes is called a cluster of miRNA. miRNAs that create the same cluster are transcribed in the same orientation without separation by any transcription unit. Clusters of miRNAs take part in controlling various important cellular processes such as tumor formation and development of lungs, heart, and immune system [18]. Changes in miRNA expression have been reported in many diseases, including different cancers. An miRNA profiling study showed that the pattern of miRNA expression in cancer tissues compared to normal tissues was very different [19]. miRNA can perform suppressor or oncogenic functions in cancers. As a tumor suppressor, it may inhibit cancer growth by negative oncogenes regulation; generally, suppressor miRNA expression is decreased in cancer cells [20]. On the contrary, the oncogene miRNAs expression (called “oncomirs”) is increased in tumors. Oncomirs promote tumors development and progression by the mechanism of downregulating a tumor suppressor or other genes involved in, for example, cell proliferation or cell differentiation in cancer [19]. In ovarian cancer, expression of different suppressors (miR-26b [21], miR-29b [22], miR-30a-5p [23], miR-93-5p [24], miR-106b-5p [25], miR-125b [26]) and oncogenes (miR-183 [27], miR-376a [28], miR-383-5p [29], miR-551b-3p [30], miR-572 [31]) was described. Ovarian cancer tissue studies revealed the negative correlation between members of the miR-23 family. High expression level of miR-23a and low expression level of miR-23b were connected with the stage and the degree of cancer, the differentiation level, and the lymph nodes metastases [32]. Much attention is currently being paid to miRNA in the context of drug resistance in ovarian cancer [13,19,20]. Changes in miRNA expression are associated with the development of resistance to chemotherapy in patients [33] as well as in vitro in the sensitive-resistant cell line model [33]. In PAC-resistant patients, downregulation of miR-663 and miR-622 was correlated with better prognosis [33]. Overexpression of miR-647 in the PAC-sensitive patients was correlated with good prognosis, which suggested suppressor function of miR-647 [33]. Changes in miRNA expression were also observed in the cell line study. Downregulation of miR-31 expression correlated with taxane resistance in ovarian cancer cell lines [34]. In contrast, upregulation of miR-98-5p was observed in cisplatin (CIS)-resistant cell lines [35]. Changes in miRNAs expression were also noted in another cancer. miR-195 expression was downregulated in temozolomid-resistant glioma cells [36], and miR-203 was downregulated in prostate cancer cells resistant to doxorubicin (DOX) [37]. The use of miRNA microarrays to analyze changes in miRNAs expression is an effective molecular tool for discovering new miRNAs involved in drug resistance processes. The present study shows alterations in the miRNA expression levels in CIS (W1CR), PAC (W1PR1 and W1PR2), DOX (W1DR), and topotecan (TOP) (W1TR) resistant variants of the W1 sensitive ovarian cancer cell line. We also identified a set of target genes that can be responsible for drug resistance in these cell lines.

2. Results

2.1. Gene Chip Quality Assessment

In the present study, we used standard factors such as signal-to-noise ratio internal hybridization and controlled spike-in-controls to preliminarily determine the quality of analyzed samples. Controlled spike-in-controls were: spike_in-control-2, spike_in-control-23, spike_in-control-29, spike_in-control-31, spike_in-control-36. Oligos 2, 23, and 29 are RNA, confirming the poly(A) tailing and the ligation. Oligo 31 (poly(A) RNA) confirmed ligation. Oligo 36 is poly(day) DNA and confirmed ligation and the lack of RNases in the RNA sample.

2.2. Gene Expression Evaluation and Gene Expression Lists

We evaluated changes in the transcription level of miRNAs. Analysis of miRNA expression in five drug resistant ovarian cancer cell lines provided new information about the significance of changes in miRNAs expression in drug resistance development in ovarian cancer. Table 1 summarizes the changes in the expression of miRNAs in drug resistant sublines with respect to the drug sensitive W1 cell line. Statistically significant level changes higher than 5-fold and less than 0.2-fold (up-/down-regulation of more than/less than 5 and −5, respectively) in the drug-resistant cells relative to those in their drug-sensitive counterparts were evaluated. miRNAs with expression levels of between 5- and 0.2-fold of those of the controls were considered “not significant (NS)” when the gene lists were constructed.

2.3. miRNAs Expression in Drug Resistant Cell Lines

The general profile of miRNA expression in the appropriate experimental group is presented in Figure 1.
We observed alterations in the expression of 40 miRNAs (Table 1, Figure 2). The level of expression of 21 miRNAs was upregulated in at least one drug resistant cell line. Expression of 19 miRNAs was downregulated in at least one drug resistant cell line.
In relation to the determined cut-off criteria (fold ±5, p < 0.05), the most alternations were observed for W1PR1 (PAC resistant cell line), which showed 25 miRNAs with altered expression (14 miRNAs upregulated and 11 miRNAs downregulated). Changes in expression of 12 miRNAs were observed for W1PR2 cell line with five genes upregulated and seven genes downregulated. Much fewer changes were observed for other drug resistant cell lines. In the CIS resistant cell line, we observed changes in nine miRNAs (five miRNAs upregulated and four miRNAs downregulated). In the TOP resistant cell line, changes in eight miRNAs expressions were observed (three miRNAs upregulated and five miRNAs downregulated), and finally, in the DOX-resistant cell line, we observed changes in the expression of seven miRNAs (one miRNA upregulated and six miRNAs downregulated).
Based on the following standards, we selected 27 of 40 miRNAs for further analysis: (1) miRNAs with changes in expression at least 5-fold (up or down) in at least three drug resistant cell lines; (2) at least 10-fold changes in miRNAs expression in one drug resistant cell line; (3) 5-fold changes in expression of at least two miRNAs from the same family. Expression of miR-31 and miR-449c changed in four drug resistant cell lines (W1PR1, W1PR2, W1CR, and W1TR). Expression of miR-105, miR-1271, miR-152, miR-767-5p, and miR-875-3p changed in three drug resistant cell lines. Expression of 12 miRNAs was changed at least 10-fold in one tested line, and between them, expression of miR-10a, miR-146b-5p, miR-18b, miR-214, miR-29a, and miR-363 was downregulated, and expression of miR-145, miR-195, miR-203, miR-30a, miR-335, and miR-99a was upregulated.
Expression of miRNA from families 181 (miRNA-181a and miRNA-181b), 193 (miRNA-193a-5p and miRNA-193b), 199 (miRNA-199a-3p and miRNA-199b-3p), and 23 (miRNA-23b and miRNA-23c was also changed in investigated cell lines. Importantly, the expression of miRNAs belonging to the same family was altered in the same way. In particular, members of 193 and 23 families were upregulated in the W1PR1 cell line and members of 181 and 199 families were downregulated in the W1TR and the W1DR cell lines, respectively.
A high similarity between cell lines resistant to PAC was observed. In both cell lines, the same six miRNAs (miR-105, miR-10a, miR-31, miR-449c, miR-767-5p, and miR-875-3p) were downregulated, and miR-152 was upregulated. Some similarity was also observed between both PAC and CIS resistant cell lines—the drug used in the first line chemotherapy of ovarian cancer.
Of the 27 miRNAs selected for analysis, nine showed a very significant level of expression change with >20-fold. miR-31 was downregulated in four cell lines (W1PR1, W1PR2, W1CR, and W1TR) with very strong downregulation of more than 217-fold in the W1PR1 cell line and strong downregulation of 52-fold in the W1PR2 cell line and 36-fold in the W1TR cell line. Among other significantly downregulated miRNAs, we observed downregulation of miRNA-10a (34-fold), miRNA-146b-5p (39-fold), and miRNA-767-5p (20-fold) in the W1PR1 cell line. We also observed a significant downregulation of miR-105 in the W1PR2 (22-fold) cell line and miRNA-29a (22-fold) in the W1TR cell line as well as downregulation of miRNA-199a (23-fold), 199b (28-fold), and miRNA-214 (26-fold) in the W1DR cell line. Significant upregulation of miRNA-203 (22-fold) was observed in the W1CR cell line and of miRNA-335 (21-fold) in the W1TR cell line.

2.4. Analysis of Target Genes Expression

In the second part of our investigation, we were interested in whether described miRNAs are involved in the regulation of genes responsible for drug resistance development. By the assumption that increasing miRNA expression leads to decreasing target gene expression and vice versa, for further analysis, we selected those target genes in which the fold change was inversely correlated with the fold change of miRNA and changes at least five times up/down, respectively, with adj.p.val. <0.05. Targets with expression levels between 5- and 0.2-fold (up/down regulation between 5 and −5) compared to those of the controls were considered “not significant (NS)” when the target gene lists were constructed. For analysis of target genes expression, we used our microarray data that were published previously [38,39,40]. From target genes analysis, we excluded miRNAs expressed only in the W1PR2 cell line, because we do not possess mRNA microarray data for this cell line. Therefore, the following miRNAs were excluded from analysis: miR-1285, miR-30a, miR-4721, miR-501-5p, and miR-572. Using our criteria (changes in expression 5-fold up or down), we did not find any targets for the following miRNA: miR-1271, miR-152, miR-203, miR-335, and miR-214.
For further analysis, we selected genes involved in drug resistance, extracellular matrix, and cancer stem cell biology using the following key words in the Gene Ontology (GO) Database: response to drug, drug transport, extracellular space, extracellular matrix, collagen containing extracellular matrix, stem cell; previously, we described changes in expression of genes from these groups in investigated drug resistance cell lines [37,38,39] Using these criteria, we did not find targets for the following miRNA: miR-100, miR-10a, miR-143, miR-193a-5p, miR-199a-3p, miR-199b-3p, miR-4269, and miR-99a.
In the CIS resistant cell line, we identified targets for miR-767-5p, miR-195, and miR-708. Among them, miR-195 overexpression correlated with SLC2A14 (solute carrier) transporter downregulation (Figure 3).
In the DOX resistant cell line, we found targets for miR-146b-5p, miR-205, and miR-875-3p, and, among them, we could distinguish Semaphorin 6A (SEMA6A) regulated by miR-875-3p (Figure 4).
Many more targets were identified in the PAC and the TOP resistant cell lines and among them were important genes related to drug resistance or other features of cancer cells. In the W1PR1 cell line, we observed inverse correlation between expression of miR-363 and key drug resistant gene ABCB1 (ATP binding cassette subfamily B member 1). Expression of important collagens—COL3A1 (collagen type III alpha 1 chain) and COL5A2 (collagen type V alpha 2 chain)—was correlated with miR-767-5p. Overexpression of receptor tyrosine kinases EPHA7 (ephrin type-A receptor 7) was associated with downregulation of miR-18 and miR-20b. miR-20b and miR-146b-5p downregulation was correlated with upregulation of MAP3K8 (mitogen-activated protein kinase kinase kinase 8). SEM3A (semaphorin 3A) and MYC (myelocytomatosis viral oncogene homolog) downregulation was associated with overexpression of miR-145. MYC expression was also regulated by let-7c. Upregulation of miR-497 and miR-195 was associated with PCDH9 (protocadherin 9) downregulation (Figure 5).
Changes in different miRNAs were also associated with changes in expression of important drug resistant genes in the W1TR cell line. miR-181a and miR-181b downregulation was associated with increased expression of SPP1/OPN (secreted phosphoprotein 1/osteopontin), TGFBI (transforming growth factor beta induced), LOX (lysyl oxidase), EPHA7, DCN (decorin), and MAP3K8 (only miR-181b). Expression of SPP1 was also associated with miR-449c downregulation, and miR-449c downregulation correlated with SEMA3D (semaphorin 3D) overexpression. Similarly, in the W1PR1 cell line, MYC downregulation was associated with overexpression of miR-145. Very high upregulation of COL3A1 was associated with downregulation of miR-29a in this cell line (Figure 6).
A negative correlation between miRNAs and their mRNA targets was verified using the data from the The Cancer Genome Atlas (TCGA) database, where we selected mRNA and miRNA expression data for 139 patients treated with cisplatin. It should be noted that this group was not homogeneous because the majority of patients underwent chemotherapy with other cytostatic agents. From the whole group, we received data for 96 patients with person neoplasm cancer status = “with tumor”, considered as chemotherapy resistant patients. A total of 43 patients with person neoplasm cancer status = “tumor free” were considered as chemotherapy sensitive cases. Negative correlations with appropriate R and p < 0.05 are presented in Figure 7.

3. Discussion

The present study demonstrates the correlation between resistance to cytotoxic drugs and expression of genes encoding microRNAs.
It was reported that changes in miRNAs expression can be related to the development of chemotherapy resistance in solid tumors as well as the development of resistance of cancer cell lines to cytotoxic drugs in vitro [41]. Thus, in the present study, we analyzed changes in the expression of several miRNAs in ovarian cancer cell lines resistant to drugs from the first line (CIS and PAC) as well as from the second line of chemotherapy (DOX and TOP). We were especially interested in changes of miRNAs expression in cancer cells exposed to drugs with different mechanisms of action. In contrast to studies in which only one pair of sensitive/resistant cell lines is tested, in this study, changes in miRNA expression in five drug resistant cell lines with PAC-resistant twin cell lines were examined. It is a unique model for such research.
Among the 27 miRNAs analyzed in this experiment, we could observe significant downregulation of miR-31 and miR-449c in four out of five drug resistant cell lines. Such a significant decrease in expression in cell lines resistant to CIS, TOP, and—especially noteworthy—in both PAC resistant cell lines may indicate a role in the process of resistance to these drugs. This is supported by other studies in which a downregulation of miRNA-31 induced taxane resistance in ovarian cancer cells. miR-31 was downregulated in a PAC resistant cell line, and reintroduction of miR-31 again sensitized them to PAC both in vitro and in vivo [34]. The miR-31 suppressor effect in the ovarian cancer cell line and tissue was confirmed by Hassan M.H. et al. [42]. They observed a high level of miR-31 in the chemosensitive serous ovarian cancer cell line and a decrease in the taxane-resistant line [42]. All these results suggest that downregulation of miR-31 expression may be associated with PAC resistance in ovarian cancer. The potential role of reduced miR-31 expression in CIS resistance can also be assumed. This is demonstrated by the results of our study, in which significant downregulation in the CIS-resistant cell line was noted. This is in line with the results of other studies that showed decreased miR-31 levels in CIS-resistant tumors and the cell line of gallbladder cancer (GBC). In contrast, the ectopic expression of miR-31 in a CIS-resistant cell line causes an increased CIS sensitivity [43]. Because, in our study, for the first time, we observed downregulation of miR-31 expression in a TOP-resistant cell line, the potential contribution of miR-31 with the development of TOP resistance requires further investigation.
Unlike miR-31, expression of miR-449c has not been described in the literature in the context of drug resistance in cancer. Since we demonstrated a downregulation of miR-449c in cell lines resistant to different cytotoxic agents, this may suggest its non-specific role in drug resistance. Previously, we observed ALDH+ stem cell-like population in W1PR1, W1TR, and W1PR2 cell lines [44,45]. Here, we observed downregulation of miR-31 and miR-449c in all these cell lines, and, importantly, the level of these miRNAs downregulation was inversely correlated with the number of ALDH+ cell populations. This suggests the role of miR-31 and miR-449c downregulation in the development of the stem cell population. Lv C. et al. showed that miR-31 plays a critical role in mammary stem cell self-renewal and breast tumorigenesis by regulating the Wnt pathway [46].
Increased expression of miR-1271 in our PAC-, CIS-, and TOP-resistant cell lines contradicts the results of other investigators. Expression of miR-1271 was downregulated in gastric cancer and cell lines in comparison to healthy tissues and cells and was further downregulated in the CIS-resistant cell lines [47]. Glioblastoma patients with a low level of miR-1271 have shorter overall survival than patients with a higher level [48].
We demonstrated significant downregulation of three subsequent miRNAs (miR-105 and miR-767-5p in CIS and PAC-resistant and miR-875-3p in PAC- and TOP- resistant cell lines). Lu G. et al. showed that, in NSCLC (Non-Small Cell Lung Cancer), low miR-105 expression was associated with poor patient survival, suggesting its hypothetical role in chemotherapy resistance [49]. In contrast, in triple negative breast cancer, it was shown that miR-105 was upregulated and correlated with poor patient survival and promoted chemoresistance [50]. However, clinical study does not reflect a simple drug sensitive/resistant cell line model. In contrast, we did not find any data about the role of miR-767-5p and miR-875-3p in drug resistance. Thus, the role of miR-767-5p and miR-875-3p in drug resistance can be considered possible, but drawing specific conclusions requires more thorough research.
Upregulation of two other miRNAs (miR-152 and miR-195) in cell lines resistant to CIS and PAC could indicate their involvement in the regulation of genes associated with resistance to the first line chemotherapy of ovarian cancer. However, in contrast to our study, downregulation of miR-152 and miR-195-5p was observed in CIS-resistant ovarian cancer cell lines [51,52]. On the other hand, the upregulation of miR-195 was observed in the PAC-resistant laryngeal cancer cell line [53] and in docetaxel resistant head and neck squamous cell carcinoma [54]. Thus, the role of these miRNAs in drug resistance could be cell and drug type dependent.
Based on the second research hypothesis assumed in advance, a group of miRNAs whose expression level was changed at least 10-fold in one drug resistant cell line was selected. In this group, as many as nine different miRNAs (miR-10a, miR-146b-5p, miR-203, miR-30a, miR-363, miR-18b, miR-145, miR-195, and miR-99a) showed changes in the expression level in PAC resistant cell lines. However, we could also observe high changes in the expression of single miRNAs in CIS- (miR-195, miR-203), TOP- (miR-29a, miR-145, miR-335), and DOX-resistant (miR-146b-5p, miR-214) cell lines. We are aware that the obtained results cannot clearly indicate the role of these molecules in the development of drug resistance. Some of them were described here for the first time (miR-99a, miR-145, and miR-18b). For others, we searched any literature data describing their involvement in the context of drug resistance or, if not found, at least in the pathogenesis of various cancers.
As the result of literature search, we found that other research teams described the changes in expression levels of the same miRNAs as ours. In contrast to our result, where a decreased level of miR-10a in the PAC-resistant cell line was noted, Sun W. et al. observed an increased expression level in the CIS-resistant A549 lung cancer cell line. Furthermore, miR-10a silencing increases sensitivity of resistance in the cell line to CIS [55]. These differences can result from different drugs and different cancer types used in both studies. Further studies report that low level of miR-146b-5p expression correlates with ovarian cancer progression. In advanced cancer stages (III/IV), the level of miR-146b-5p expression was lower compared to earlier stages (I/II) [56]. In the cell line study, overexpression of miR-146b-5p enhanced ovarian cancer cell sensitivity to PAC and CIS [56]. Thus, miR-146b-5p downregulation seems to play a role in the resistance of ovarian cancer to cytotoxic drugs, especially to PAC.
In opposition to our results (upregulation in PAC and CIS resistant cell lines), Cheng R. et al. showed that, in CIS-resistant NSCLC tumors, miR-203 was downregulated compared with CIS-sensitive tumors [57]. Ectopic expression of miR-203 in the NSCLC cell line increased CIS sensitivity and apoptosis [57]. Overexpression of miR-203 also reversed PAC-resistance in colorectal cancer [58]. Regarding miR-30a and the increase in expression we observed in the PAC-resistant cell line, Sestito et al. showed miR-30a was downregulated in CIS- and PAC-resistant ovarian cancer cell lines [23], and overexpression of miR-30a caused increased sensitivity of ovarian cancer cells to CIS-induced apoptosis [23].
In contrast to our report, Mohamed Z. et al. [59] observed upregulation of miR-363 in the PAC resistance ovarian cancer cell line. Furthermore, overexpression of miR-363 confers PAC resistance, and inhibition of miR-363 restores the response to PAC [59]. In contrast, downregulation of miR-363 was observed in the CIS-resistant hepatocellular carcinoma cell line [60]. A significantly lower expression level of miR-363 was also observed by Cao L. et al. in the CIS-resistant ovarian cancer cell line [61]. Thus, the significance of miR-363 in drug resistance can be cell line dependent.
The next three miRNAs described in this report hypothesized their involvement in resistance to DOX, CIS, and TOP. In the ovarian cancer cell line study, downregulation of miR-29a correlated with increased CIS resistance, while overexpression of miR-29a correlated with sensitized ovarian cancer cell line to CIS [62]. Similar results were obtained by Yang L. et al. in glioblastoma stem cells in which overexpression of miR-29a in CD133+ GBM (glioblastoma multiforme) stem cells effectively reversed the resistance to CIS [63]. Downregulation of miR-29a was also observed in the methotrexate resistant osteosarcoma cell line [64]. Here, we observed downregulation of miR-29a in the TOP-resistant cell line. Thus, it appears that a downregulation of miR-29a could be involved in cytotoxic drug resistance.
In our study, miR-335 was significantly upregulated in the TOP-resistant cell line. In contrast, downregulation of miR-335 was observed in PAC- and CIS- resistant A2780 ovarian cancer cell lines [65]. In the DOX-resistant cell line, we could observe very significant downregulation of miR-214. Downregulation of this miR was also observed in the DOX-resistant urothelial bladder cancer tissues and cells lines [66]. Thus, the miR-214 can be a potential marker of DOX resistance. We, for the first time, observed upregulation of miR-99a and miR-145 and downregulation of miR-18b in the PAC-resistant cell line. Upregulation of miR-99a was observed in CIS-resistant gastric cancer cell lines [67] and vincristine-resistant acute lymphoblastic leukemia (ALL) [68]. In contrast to our study, downregulation of miR-145 was observed in PAC-resistant variants of A2780 and SKOV-3 ovarian cancer cell lines [69]. We did not find any information about the significance of miR-18b in drug resistance in general. Thus, its role in drug resistance requires further investigation.
It is assumed that miRNAs from the same family can play a very similar role in regulation of genes expression. Therefore, as the last step in our research, we analyzed changes in miRNAs expression in the same families. In our study, miR-193a-5p and miR-193b were upregulated in the PAC-resistant cell line, which suggests they can regulate the expression of genes involved in resistance to this drug. As demonstrated by Khordadmehr et al., the miR-193 family (miR-193a-3p, miR-193a-5p, miR-193b3p, and miR-193b-5p) plays an important role in the pathogenesis of ovarian cancer [70]. A similar result was obtained for miR-23 family members, miR-23b and miR-23c. As there are no literature data describing the contribution of these miRNAs as a cluster to cytostatic resistance, at this stage, one can only assume that they may affect the expression of genes associated with PAC resistance. However, the upregulation of miR-23b correlated with better survival of lung cancer patients [71] and, on the other hand, with DOX-resistance in thyroid cancer cells [72]. In another study, downregulation of miR-23b and miR-27b was observed in multidrug resistance protein 1 (MDR) Ehrlich ascites tumor cells [73]. Thus, the role of these miRs in drug resistance requires further investigation. On the other hand, the downregulation of miR-181a and miR-181b in the TOP-resistant ovarian cell line might be an indicator of their involvement in the development of TOP resistance. However, their role in TOP-resistance has not been investigated by others. Decreased expression of miR-181a was observed in the CIS-resistant cervical cancer cell line [74] and miR-181b in the NSCLC cell CIS-resistant cell line [75].
Members of the miR-199 family were described in ovarian cancer. Cheng W. et al. showed that miR-199a-3p and miR-199a-5p overexpression significantly decreased chemoresistance of cancer-initiating cells (CICs) to CIS, PAC, and DOX and reduced mRNA expression of the multidrug resistance gene ABCG2 in vivo [76]. miR-199b-5p silencing is significantly associated with acquired chemoresistance in ovarian cancer cell lines and cancer tissues [77]. Here, we observed significant downregulation of miR-199a-3p and miR-199b-3p in a DOX-resistant cell line. However, as previously reported, ABCG2 gene expression did not change in this cell line [38]. Determining the involvement of miR-199 family members in the regulation of drug resistant genes requires further investigation.
As presented above, both our results and studies of other researchers indicate some involvement of specific miRNAs in the process of drug resistance. However, the interpretation of miRNA expression results is much more complicated compared to gene encoding proteins. As an example, the expression of MDR1 gene encoding P-gp protein is always related to a drug resistant phenotype. In the case of miRNA genes, this phenomenon is more comprehensive. On the one hand, one miRNA can regulate the expression of different mRNAs, and, on the other hand, the same mRNA can be regulated by different miRNAs. This generates a complex network of dependencies.
Considering the high complexity of the process of acquiring resistance by cancer cells and summarizing the results obtained in this study, we believe that the contribution of miR-31, miR-146b, and miR-195 in PAC- and of miR-214, miR-199a-3p, and miR-199a-5p in DOX resistance is highly possible. However, the role of the described miRNAs requires further molecular, cellular, and bioinformatics analysis, which we intend to proceed with in the next step of our study.
To shed more light on the role of investigated miRs in drug resistance, in the second part of our study, we analyzed the correlation between the expression of miRs and their targets. We were mainly interested to see if genes described previously in the context of drug resistance could be regulated by miRNAs, although some other genes were also discussed.
In our previous microarray analysis [38,39,40,78,79], we could observe many more genes with changed expression in PAC- and TOP-resistant cell lines than in CIS- and DOX-resistant ones. This was also reflected by the number of targets to miRNAs identified in the Gene Expression Omnibus (GEO) database in these cell lines. Thus, according to the number of identified targets, we could divide our cell lines into those with a low number of identified targets, represented by CIS- and DOX-resistant cell lines, and those with a high number of identified targets, represented by PAC- and TOP-resistant cell lines.
In the W1CR1 cell line, the only gene that was described in the context of CIS-resistance was SLC2A14 regulated by miR-195 [38]. The only gene described in the context of drug resistance in the DOX-resistant cell line was SEMA6A regulated by miR-875-3p. However, in contrast to our results, downregulation of SEMA6A was observed in the DOX-resistant A2780 ovarian cancer cell line [80]. Thus, the role of this miR and its target in drug resistance requires further investigation.
As we described previously [81], the most important gene responsible for PAC resistance in the W1PR1 cell line was ABCB1 (MDR1) encoding the glycoprotein P (P-gp). Here, we observed that upregulation of this gene is related to downregulation of miR-363 in this cell line, which is the first such observation. However, downregulation of miR-363 was associated with ABCC1 overexpression in oxoplatin hepatocellular carcinoma cell lines [82].
Previously, we reported overexpression of several collagen genes in different drug resistant ovarian cancer cell lines [83]. In this experiment, we could observe that increased expression of COL3A1 and COL5A2 inversely correlates with miR-767-5p level. Thus, expression of these two genes seems to be regulated by this miRNA. Similarly, we observed downregulation of miR-767-5p in the W1PR2 cell line accompanied by very strong COL3A1 upregulation [83]. This further supports the possible role of miR-767-5p in COL3A1 gene regulation.
Drug resistant cells are characterized by strong signal transduction. Previously, we reported increased total pTYR (phosphotyrosine) level in drug resistant cell lines in comparison to parental cell lines W1 [84]. Here, we observed that upregulation of receptor tyrosine kinase—EPHA7 and MAP3K8 kinase—is inversely correlated with miR-18b (EPHA7) and miR-18b and miR-20b together (EPHA7 and MAP3K8). Regulation of EPHA7 expression by miR-18b was also described by others [85]. It was reported that signal transduction by the EPHA7 receptor may activate the extracellular-signal-regulated kinase (ERK) signaling pathway that also involves MAP3K8 [86]. Upregulation of EPHA7 was observed in many cancers, including hepatocellular carcinoma, where it was associated with tumor progression, invasiveness, and metastasis [87]. MAP3K8 was described as a mediator of vemurafenib resistance of thyroid cancer stem cells [88].
SEMA3A is a member of the semaphorin family of proteins that play a role in many developmental processes as well as in cancer progression [89]. It is a tumor suppressor gene, and its downregulation correlates with progression of gastric [90] and ovarian cancer [91]. Recently, we described SEMA3A downregulation in three PAC-resistant ovarian cancer cell lines [92]. Here we showed that expression of SEMA3A can be influenced by miR-145 upregulation in the PAC resistant cell line. Regulation of SEMA3A by miR-145 is not surprising because SEMA3A as a direct target of miR-145 was described recently [93]. The upregulation of miR-145 was also discovered to be connected with strong downregulation of the Myc proto-oncogene protein in W1PR1 cell lines. miR-145 was described as a suppressor miRNA that directly targets the MYC transcript in esophageal squamous cell carcinoma [94] and ovarian cancer [95]. Moreover, as our analysis shows, the MYC downregulation can be due to let-7c upregulation in this cell line as well, which was described previously in prostate cancer cells [96].
Another gene found to be downregulated in the PAC resistant cell line was PCDH9, the expression of which correlated with miR-195 and miR-497 upregulation. PCDH9 (protocadherin 9), a member of the protocadherin protein family [97], is a calium-dependent adhesion protein playing a role in neural cell interaction [98]. It is a tumor suppressor gene with decreased expression in ovarian [99] and gastric cancer [100], and its downregulation usually correlates with disease progression and shorter survival. Previously, we observed the PCDH9 downregulation in three PAC resistant cell lines [92], suggesting its significance in resistance to this cytotoxic agent. However, there are no literature data available describing the regulation of the PCDH9 gene by miR-195 or miR-497.
The very important role in the TOP-resistant cell line was downregulation of miR-181a and miR-181b, since their downregulation was inversely correlated with upregulation of key genes involved in drug resistance. As in the W1PR1 cell line, in the W1TR cell line, we observed an increased total pTYR level [84]. Additionally, increased EPHA7 and MAP3K8 levels were noted. It can be associated with stronger signal transduction, increased pTYR level, and eventual drug resistance. As opposed to the W1PR1 cell line, both genes were regulated here by miR-181a and miR-181b, which is the first such observation described thus far.
We could also observe that both miR-181a and miR-181b regulated DCN and LOX genes’ expression. DCN (decorin) is a small leucine reach proteoglycan, a part of ECM that binds to collagens and plays an important role in cancer development and metastasis [101]. Regulation of DCN by miR-181 was reported in skin and wound healing [102]. LOX (lysyl oxidase) is a secretory protein involved in cross-linking of collagens and elastin in the ECM that results in increased ECM stabilization [103]. Its expression seems to be an important metastatic factor in breast [104] and ovarian [105] cancers, among others. Previously, we described increased expression of LOX in PAC and TOP resistant ovarian cancer cell lines, and its expression was upregulated in ALDH1A1 positive cells [45]. However, the regulation of the LOX gene by miR-181a or miR-181b has not been described in literature thus far.
Two other genes important in drug resistance and cancer progression were also regulated by miR from the 181 family as well as miR-449c. TGFBI (transforming growth factor-beta-induced protein) is a collagen I, II, and IV binding protein present in ECM. It has been described as the metastasis promoting protein in ovarian cancer [106], and its expression was correlated with shorter survival of patients with serous ovarian cancer [107]. Recently, we also described TGFBI as a gene related to TOP resistance in three different ovarian cancer cell lines [38]. It was proven that expression of TGFBI is regulated by miR-181a in osteoblasts during differentiation [108].
SPP1 (secreted phosphoprotein 1), also known as osteopontin (OPN), was originally described in bone tissue [109]; however, its expression was also described in different cancers, including ovarian [110], as a protein involved in tumor progression, metastasis, and drug resistance [111]. It was reported that PAC-, DOX-, and CIS-resistance were related to the induction of drug transporters protein expression by SPP1 [112] or blocking caspase [113]. Recently, we also described expression of the SPP1 protein in three TOP resistant cell lines [38]. Similar to our results, the downregulation of miR-181a and the upregulation of OPN mRNA were observed in CIS-resistant cervical cancer cells [74], and overexpression of miR-181a led to reduced expression of OPN and higher sensitivity to CIS. Regulation of OPN by miR-181b was described also in eosinophilia in asthma [114]. In contrast, regulation of SPP1 by miR-449c has not been described thus far.
As we described previously, the W1TR cell line was characterized by abundant COL3A1 expression, which was secreted from the cells [83]. COL3A1 was also described as a protein involved in CIS resistance in ovarian cancer [115]. Here, we observed that upregulation of COL3A1 inversely correlated with downregulation of miR-29a. In a similar way, upregulation of COL3A1 was observed in two MTX resistant osteosarcoma cell lines and correlated with miR-29a downregulation. Furthermore, it was proven that miR-29a regulates COL3A1 expression [64].

4. Materials and Methods

4.1. Reagents

Cisplatin, Doxorubicyn, Topotekan, and Paclitaxel were obtained from Sigma (St. Louis, MO, USA). RPMI-1640 medium, fetal bovine serum, penicillin, streptomycin, amphotericin B (25 µg/mL), and L-glutamine were also purchased from Sigma. QIazol Lysys Reagent, miRNeasy Mini Kit, and RNeasy MinElute Cleanup Kit were obtained from Qiagen (Hilden, Germany). GeneChip™ miRNA 3.1 Array Strip, FlashTag™ Biotin HSR RNA Labeling Kits, GeneAtlas™ Hybridization, Wash, and Stain Kit for miRNA Arrays were obtained from Affymetrix (Santa Clara, CA, USA).

4.2. Cell Lines and Cell Culture

The human primary ovarian cancer cell line W1 was established from the tumor tissue of an untreated 54-year-old Caucasian female patient diagnosed for serous ovarian adenocarcinoma (G3, FIGO IIIc), approval of Bioethical Committee at the Medical University of Poznan of 14 June, 2012, number 699/12. Cells grow as a monolayer, present epithelial morphology, and adherent growth model. Sublines resistant to cisplatin (W1CR), doxorubicin (W1DR), topotecan (W1TR), and paclitaxel (W1PR1 and W1PR2) were obtained by exposure of the W1 cell line to stepwise increasing drug concentrations. Final concentrations of each drug were twofold greater than the concentration in the plasma 2 h after intravenous administration. The cells were 8-, 10-, 20-, 641-, and 967-folds resistant to their selective drugs, respectively, as determined by Cell Proliferation Kit I (MTT) [80]. All cell lines were maintained as a monolayer in complete medium (RPMI-1640 medium supplemented with 10% (v/v) fetal bovine serum, 2 pMl-glutamine, penicillin (100 U/mL), streptomycin (100 U/mL), and amphotericin B (25 µg/mL)) at 37 °C in a 5% CO2 atmosphere.

4.3. miRNA Isolation

miRNA was isolated using a reagents kit from Qiagen according to the manufacturer’s protocol. The RNA was quantified using spectrophotometry by measuring the absorbance values at 260 nm and 280 nm, and the 260/280 nm ratio was used to estimate the level of protein contamination. The 260/280 nm ratios of the samples ranged from 1.8 to 2.0. Four RNA samples from each experimental group were subjected to miRNA expression profiling using microarray method. (N/group = 4).

4.4. Microarray Preparation, Hybridization, and Scanning

miRNA expression profiling was performed using the microarray approach with Applied BiosystemsTM miRNA 3.1 Array Strip (ThermoFisher Scientific, Waltham, MA, USA). The detailed technical procedure was described earlier [116,117]. Each microarray was designed in accordance with the miRBase Release 17 database, including complementary probes for: 1733 human mature miRNA, 2216 human snoRNA, CDBox RNA, H/ACA Box RNA, and scaRNA, 1658 human pre-miRNA. The full procedure for preparing miRNA for hybridization was performed using the FlashTagTM Biotin HSR RNA Labeling Kit (ThermoFisher Scientific, Waltham, MA, USA). Then, 150 ng of previously isolated miRNA was subjected to the poly(A) tailing and biotin ligation procedure, according to the manufacturer’s protocol. Biotin-labeled miRNA were hybridized to Applied BiosystemsTM miRNA 3.1 Array Strip (20 h, 48 °C). Afterwards, the microarrays were rinsed and stained according to the technical protocol using the Affymetrix GeneAtlas Fluidics Station (Affymetrix, Santa Clara, CA, USA). The array strips were scanned using an Imaging Station of GeneAtlas System (Thermo Fisher Scientific, MA, USA).

4.5. Microarray Analysis and miRNA Gene Screening

The initial analysis of the scanned array strips was performed with Affymetrix GeneAtlas Operating Software (Affymetrix, Santa Clara, CA, USA). The quality of the miRNA gene expression data was verified using quality control criteria set by the software. The obtained CEL files from the scanned microarrays were imported for further data analysis using BioConductor libraries from the “R” statistical programming language. The Robust Multiarray Average (RMA) normalization algorithm (implemented in the Affy library) was used to normalize, correct the background, and calculate the expression values of all the tested miRNAs [118]. The biological annotation was obtained from the pd.mirna 3.1 library, where the annotated data frame object was merged with a normalized data set, leading to a complete miRNA data table. Differential expression and statistical evaluation were determined using a linear model for microarray data implemented in the "limma" library [119]. P-values were obtained by the empirical Bayes moderated t-test with false discovery rate (FDR) correction for multiple testing. The algorithm of such a correction was implemented in the R "limma" library. The selection criteria for significantly changed miRNA gene expression were based on the difference between folds greater than absolute five and FDR adjusted p value (adj.p.val) < 0.05. The results of such a selection were presented as a scatter plot diagram showing the total number of miRNAs up- and downregulated. Differentially expressed miRNA were also visualized as heat maps and tables. Raw data files were deposited in the Gene Expression Omnibus (GEO) repository at the National Center for Biotechnology Information (http://www.ncbi.nlm.nih.gov/geo/) under the GEO accession number GEO: GSE148251.

4.6. miRNA-Target Gene Prediction

SpidermiR package was applied to identify potential target genes for differently expressed miRNA. Differentially expressed miRNAs were used as a query for searching target genes in the following databases: for predicted targets—DIANA, Miranda, PicTar, TargetScan, and for experimentally validated targets—miRTAR, miRwalk [120]. To determine the actual expression value of target genes, mRNA transcriptomic data from our published experiment were used [37,38,39]. Obtained fold change values for mRNA were assigned to the target genes data table. For further analysis, we selected only those target genes for which fold change was inversely correlated with the fold change value of appropriate miRNA (cut-off criteria: fold ±5, adjusted p value (adj.p.val.) < 0.05). From the whole set of miRNA-targets pairs, we selected only those that were involved in drug resistance, extracellular matrix, and cancer stem cell biology using the following key words in GO terms: “collagen-containing extracellular matrix”, “extracellular matrix”, “extracellular space”, “response to drug”, and “stem cell”. Interactions between miRNA and target genes in the form of selected GO terms were visualized using Cytoscape 3.7.1 (National Institute of General Medical Sciences, Bethesda, MD, USA) [121].

4.7. Analysis of The Cancer Genome Atlas (TCGA) Dataset

The detailed procedure of analysis was described earlier [122]. Briefly, clinical description file, RNAseq, and miRNAseq data (n = 307 and 461, respectively) from ovarian serous cystadenocarcinoma patients were downloaded from the public TCGA database [123] using the FireBrowse server (http://gdac.broadinstitute.org/) [124]. Then, the voom algorithm from the “Limma” package was used for data normalization [119]. Using the discriminator "patient.drugs.drug. drug_name = cisplatin", we selected patients treated with cisplatin. Because there was no information on resistance to a given cytostatic in the clinical description file, we used "person_neoplasm_cancer_status" to select patients resistant and sensitive to treatment. We assumed that "person_neoplasm_cancer_status = tumor free" refers to chemotherapy sensitive, whilst "person_neoplasm_cancer_status = with tumor" refers to chemotherapy resistant patients. Negative correlation between appropriate miRNA and its target gene was performed using the Pearson’s product-moment correlation coefficient (PPMCC) test. Paired-log transformed normalized expression values were presented as scatter plots with a line of best fit. The 95% confidence interval level for predictions based on linear models was also shown. Results were presented only for the comparisons with p < 0.05

5. Conclusions

In summary, we identified a set of potential or already described targets for miRNAs with altered expression in drug resistant cell lines. Some of these targets were previously described as important factors in development of drug resistance and/or cancer progression. The way of regulation of some of these target genes by miRNAs was previously described. Regulation of others requires a detailed study at the molecular level, and significance of this regulation needs to be confirmed in the context of drug resistance.

Author Contributions

Conceptualization, methodology, investigation, and writing—original draft preparation, D.K.; methodology, investigation, and validation, K.J.; validation and supervision M.R.; review and editing, K.S.; methodology and investigation, B.G.-M.; validation and supervision, M.N.; supervision, conceptualization, review and original draft editing, R.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by grants from the National Science Centre, Poland no 2016/22/E/NZ5/00381.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

CISCisplatin
PACPaclitaxel
DOXDoxorubicin
TOPTopotecan
FIGOfr. Fédération internationale de gynécologie et d’obstétrique
ABCB1 (MDR1)ATP-binding cassette subfamily B member 1
(multidrug resistance protein 1)
W1PR1, W1PR2, W1CR, W1TR, W1DR Resistant Variants of W1 Sensitive Ovarian Cancer Cell Line
ECMExtracellular Matrix
FDRFalse Discovery Rate
GEOGene Expression Omnibus
GOGene Ontology
COL3A1 Collagen Type III Alpha 1 Chain
COL5A2Collagen Type V Alpha 2 Chain
MYC Myelocytomatosis Viral Oncogene Homolog
EPHA7 Ephrin Type-A Receptor 7
MAP3K8Mitogen-Activated Protein Kinase Kinase Kinase 8
LOXLysyl Oxidase
EPHA7Ephrin type-A receptor 7
DCNDecorin
MAP3K8Mitogen-Activated Protein Kinase 8
SEM3ASemaphorin 3A
PCDH9Protocadherin 9
SPP1/OPN Secreted Phosphoprotein 1/Osteopontin
TGFBITransforming Growth Factor Beta Induced

References

  1. Roett, M.A.; Evans, P. Ovarian cancer: An overview. Am. Fam. Physician 2009, 80, 609–616. [Google Scholar] [PubMed]
  2. Hennessy, B.T.; Coleman, R.L.; Markman, M. Ovarian cancer. Lancet (Lond. Engl.) 2009, 374, 1371–1382. [Google Scholar] [CrossRef]
  3. Webb, P.M.; Jordan, S.J. Epidemiology of epithelial ovarian cancer. Best Pract. Res. Clin. Obstet. Gynaecol. 2017, 41, 3–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Wei, W.; Giulia, F.; Luffer, S.; Kumar, R.; Wu, B.; Tavallai, M.; Bekele, R.T.; Birrer, M.J. How can molecular abnormalities influence our clinical approach. Ann. Oncol. 2017, 28, viii16–viii24. [Google Scholar] [CrossRef]
  5. Gottesman, M.M. Mechanisms of cancer drug resistance. Ann. Rev. Med. 2002, 53, 615–627. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Greenlee, R.T.; Murray, T.; Bolden, S.; Wingo, P.A. Cancer statistics, 2000. CA Cancer J. Clin. 2000, 50, 7–33. [Google Scholar] [CrossRef]
  7. Gottesman, M.M.; Fojo, T.; Bates, S.E. Multidrug resistance in cancer: Role of ATP-dependent transporters. Nat. Rev. Cancer 2002, 2, 48–58. [Google Scholar] [CrossRef] [Green Version]
  8. Leonard, G.D.; Fojo, T.; Bates, S.E. The role of ABC transporters in clinical practice. Oncologist 2003, 8, 411–424. [Google Scholar] [CrossRef]
  9. Mansoori, B.; Mohammadi, A.; Davudian, S.; Shirjang, S.; Baradaran, B. The Different Mechanisms of Cancer Drug Resistance: A Brief Review. Adv. Pharm. Bull. 2017, 7, 339–348. [Google Scholar] [CrossRef]
  10. Correia, A.L.; Bissell, M.J. The tumor microenvironment is a dominant force in multidrug resistance. Drug Resist. Updates Rev. Comment. Antimicrob. Anticancer Chemother. 2012, 15, 39–49. [Google Scholar] [CrossRef] [Green Version]
  11. Lee, R.C.; Feinbaum, R.L.; Ambros, V. The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell 1993, 75, 843–854. [Google Scholar] [CrossRef]
  12. Bartel, D.P. MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell 2004, 116, 281–297. [Google Scholar] [CrossRef] [Green Version]
  13. Mihanfar, A.; Fattahi, A.; Nejabati, H.R. MicroRNA-mediated drug resistance in ovarian cancer. J. Cell. Physiol. 2019, 234, 3180–3191. [Google Scholar] [CrossRef]
  14. Baranwal, S.; Alahari, S.K. miRNA control of tumor cell invasion and metastasis. Int. J. Cancer 2010, 126, 1283–1290. [Google Scholar] [CrossRef]
  15. Lee, Y.; Kim, M.; Han, J.; Yeom, K.H.; Lee, S.; Baek, S.H.; Kim, V.N. MicroRNA genes are transcribed by RNA polymerase II. EMBO J. 2004, 23, 4051–4060. [Google Scholar] [CrossRef] [PubMed]
  16. Kamanu, T.K.; Radovanovic, A.; Archer, J.A.; Bajic, V.B. Exploration of miRNA families for hypotheses generation. Sci. Rep. 2013, 3, 2940. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Zou, Q.; Mao, Y.; Hu, L.; Wu, Y.; Ji, Z. miRClassify: An advanced web server for miRNA family classification and annotation. Comput. Biol. Med. 2014, 45, 157–160. [Google Scholar] [CrossRef]
  18. Zhang, W.; Dahlberg, J.E.; Tam, W. MicroRNAs in tumorigenesis: A primer. Am. J. Pathol. 2007, 171, 728–738. [Google Scholar] [CrossRef] [Green Version]
  19. Lu, J.; Getz, G.; Miska, E.A.; Alvarez-Saavedra, E.; Lamb, J.; Peck, D.; Sweet-Cordero, A.; Ebert, B.L.; Mak, R.H.; Ferrando, A.A.; et al. MicroRNA expression profiles classify human cancers. Nature 2005, 435, 834–838. [Google Scholar] [CrossRef] [PubMed]
  20. Shenouda, S.K.; Alahari, S.K. MicroRNA function in cancer: Oncogene or a tumor suppressor? Cancer Metastasis Rev. 2009, 28, 369–378. [Google Scholar] [CrossRef] [PubMed]
  21. Lin, J.; Zhang, L.; Huang, H.; Huang, Y.; Huang, L.; Wang, J.; Huang, S.; He, L.; Zhou, Y.; Jia, W.; et al. MiR-26b/KPNA2 axis inhibits epithelial ovarian carcinoma proliferation and metastasis through downregulating OCT4. Oncotarget 2015, 6, 23793–23806. [Google Scholar] [CrossRef] [PubMed]
  22. Teng, Y.; Zhang, Y.; Qu, K.; Yang, X.; Fu, J.; Chen, W.; Li, X. MicroRNA-29B (mir-29b) regulates the Warburg effect in ovarian cancer by targeting AKT2 and AKT3. Oncotarget 2015, 6, 40799–40814. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Sestito, R.; Cianfrocca, R.; Rosano, L.; Tocci, P.; Semprucci, E.; Di Castro, V.; Caprara, V.; Ferrandina, G.; Sacconi, A.; Blandino, G.; et al. miR-30a inhibits endothelin A receptor and chemoresistance in ovarian carcinoma. Oncotarget 2016, 7, 4009–4023. [Google Scholar] [CrossRef]
  24. Chen, X.; Chen, S.; Xiu, Y.L.; Sun, K.X.; Zong, Z.H.; Zhao, Y. RhoC is a major target of microRNA-93-5P in epithelial ovarian carcinoma tumorigenesis and progression. Mol. Cancer 2015, 14, 31. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Chen, S.; Chen, X.; Xiu, Y.L.; Sun, K.X.; Zhao, Y. Inhibition of Ovarian Epithelial Carcinoma Tumorigenesis and Progression by microRNA 106b Mediated through the RhoC Pathway. PLoS ONE 2015, 10, e0125714. [Google Scholar] [CrossRef] [Green Version]
  26. Ying, X.; Wei, K.; Lin, Z.; Cui, Y.; Ding, J.; Chen, Y.; Xu, B. MicroRNA-125b Suppresses Ovarian Cancer Progression via Suppression of the Epithelial-Mesenchymal Transition Pathway by Targeting the SET Protein. Cell. Physiol. Biochem. Int. J. Exp. Cell. Physiol. Biochem. Pharmacol. 2016, 39, 501–510. [Google Scholar] [CrossRef]
  27. Chen, H.; Zhang, L.; Zhang, L.; Du, J.; Wang, H.; Wang, B. MicroRNA-183 correlates cancer prognosis, regulates cancer proliferation and bufalin sensitivity in epithelial ovarian caner. Am. J. Transl. Res. 2016, 8, 1748–1755. [Google Scholar]
  28. Yang, L.; Wei, Q.M.; Zhang, X.W.; Sheng, Q.; Yan, X.T. MiR-376a promotion of proliferation and metastases in ovarian cancer: Potential role as a biomarker. Life Sci. 2017, 173, 62–67. [Google Scholar] [CrossRef]
  29. Liu, J.; Dou, Y.; Sheng, M. Inhibition of microRNA-383 has tumor suppressive effect in human epithelial ovarian cancer through the action on caspase-2 gene. Biomed. Pharmacother. 2016, 83, 1286–1294. [Google Scholar] [CrossRef]
  30. Chaluvally-Raghavan, P.; Jeong, K.J.; Pradeep, S.; Silva, A.M.; Yu, S.; Liu, W.; Moss, T.; Rodriguez-Aguayo, C.; Zhang, D.; Ram, P.; et al. Direct Upregulation of STAT3 by MicroRNA-551b-3p Deregulates Growth and Metastasis of Ovarian Cancer. Cell Rep. 2016, 15, 1493–1504. [Google Scholar] [CrossRef] [Green Version]
  31. Zhang, X.; Liu, J.; Zang, D.; Wu, S.; Liu, A.; Zhu, J.; Wu, G.; Li, J.; Jiang, L. Upregulation of miR-572 transcriptionally suppresses SOCS1 and p21 and contributes to human ovarian cancer progression. Oncotarget 2015, 6, 15180–15193. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Su, L.; Liu, M. Correlation analysis on the expression levels of microRNA-23a and microRNA-23b and the incidence and prognosis of ovarian cancer. Oncol. Lett. 2018, 16, 262–266. [Google Scholar] [CrossRef] [PubMed]
  33. Kim, Y.W.; Kim, E.Y.; Jeon, D.; Liu, J.L.; Kim, H.S.; Choi, J.W.; Ahn, W.S. Differential microRNA expression signatures and cell type-specific association with Taxol resistance in ovarian cancer cells. Drug Des. Dev. Ther. 2014, 8, 293–314. [Google Scholar] [CrossRef] [Green Version]
  34. Mitamura, T.; Watari, H.; Wang, L.; Kanno, H.; Hassan, M.K.; Miyazaki, M.; Katoh, Y.; Kimura, T.; Tanino, M.; Nishihara, H.; et al. Downregulation of miRNA-31 induces taxane resistance in ovarian cancer cells through increase of receptor tyrosine kinase MET. Oncogenesis 2013, 2, e40. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Wang, Y.; Bao, W.; Liu, Y.; Wang, S.; Xu, S.; Li, X.; Li, Y.; Wu, S. miR-98-5p contributes to cisplatin resistance in epithelial ovarian cancer by suppressing miR-152 biogenesis via targeting Dicer1. Cell Death Dis. 2018, 9, 447. [Google Scholar] [CrossRef]
  36. Wang, H.; Ren, S.; Xu, Y.; Miao, W.; Huang, X.; Qu, Z.; Li, J.; Liu, X.; Kong, P. MicroRNA-195 reverses the resistance to temozolomide through targeting cyclin E1 in glioma cells. Anti-Cancer Drugs 2019, 30, 81–88. [Google Scholar] [CrossRef]
  37. Chen, L.Z.; Ding, Z.; Zhang, Y.; He, S.T.; Wang, X.H. MiR-203 over-expression promotes prostate cancer cell apoptosis and reduces ADM resistance. Eur. Rev. Med. Pharmacol. Sci. 2018, 22, 3734–3741. [Google Scholar] [CrossRef]
  38. Januchowski, R.; Sterzynska, K.; Zawierucha, P.; Rucinski, M.; Swierczewska, M.; Partyka, M.; Bednarek-Rajewska, K.; Brazert, M.; Nowicki, M.; Zabel, M.; et al. Microarray-based detection and expression analysis of new genes associated with drug resistance in ovarian cancer cell lines. Oncotarget 2017, 8, 49944–49958. [Google Scholar] [CrossRef] [Green Version]
  39. Januchowski, R.; Zawierucha, P.; Andrzejewska, M.; Rucinski, M.; Zabel, M. Microarray-based detection and expression analysis of ABC and SLC transporters in drug-resistant ovarian cancer cell lines. Biomed. Pharmacother. 2013, 67, 240–245. [Google Scholar] [CrossRef]
  40. Januchowski, R.; Zawierucha, P.; Rucinski, M.; Zabel, M. Microarray-based detection and expression analysis of extracellular matrix proteins in drug resistant ovarian cancer cell lines. Oncol. Rep. 2014, 32, 1981–1990. [Google Scholar] [CrossRef] [Green Version]
  41. Ma, J.; Dong, C.; Ji, C. MicroRNA and drug resistance. Cancer Gene Ther. 2010, 17, 523–531. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Hassan, M.K.; Watari, H.; Mitamura, T.; Mohamed, Z.; El-Khamisy, S.F.; Ohba, Y.; Sakuragi, N. P18/Stathmin1 is regulated by miR-31 in ovarian cancer in response to taxane. Oncoscience 2015, 2, 294–308. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Li, M.; Chen, W.; Zhang, H.; Zhang, Y.; Ke, F.; Wu, X.; Zhang, Y.; Weng, M.; Liu, Y.; Gong, W. MiR-31 regulates the cisplatin resistance by targeting Src in gallbladder cancer. Oncotarget 2016, 7, 83060–83070. [Google Scholar] [CrossRef] [PubMed]
  44. Januchowski, R.; Wojtowicz, K.; Sterzyska, K.; Sosiska, P.; Andrzejewska, M.; Zawierucha, P.; Nowicki, M.; Zabel, M. Inhibition of ALDH1A1 activity decreases expression of drug transporters and reduces chemotherapy resistance in ovarian cancer cell lines. Int. J. Biochem. Cell Biol. 2016, 78, 248–259. [Google Scholar] [CrossRef]
  45. Sterzynska, K.; Klejewski, A.; Wojtowicz, K.; Swierczewska, M.; Nowacka, M.; Kazmierczak, D.; Andrzejewska, M.; Rusek, D.; Brazert, M.; Brazert, J.; et al. Mutual Expression of ALDH1A1, LOX, and Collagens in Ovarian Cancer Cell Lines as Combined CSCs- and ECM-Related Models of Drug Resistance Development. Int. J. Mol. Sci. 2018, 20, 54. [Google Scholar] [CrossRef] [Green Version]
  46. Lv, C.; Li, F.; Li, X.; Tian, Y.; Zhang, Y.; Sheng, X.; Song, Y.; Meng, Q.; Yuan, S.; Luan, L.; et al. MiR-31 promotes mammary stem cell expansion and breast tumorigenesis by suppressing Wnt signaling antagonists. Nat. Commun. 2017, 8, 1036. [Google Scholar] [CrossRef] [Green Version]
  47. Yang, M.; Shan, X.; Zhou, X.; Qiu, T.; Zhu, W.; Ding, Y.; Shu, Y.; Liu, P. miR-1271 regulates cisplatin resistance of human gastric cancer cell lines by targeting IGF1R, IRS1, mTOR, and BCL2. Anti-Cancer Agents Med. Chem. 2014, 14, 884–891. [Google Scholar] [CrossRef] [Green Version]
  48. Yang, L.; Wang, Y.; Li, Y.J.; Zeng, C.C. Chemo-resistance of A172 glioblastoma cells is controlled by miR-1271-regulated Bcl-2. Biomed. Pharmacother. 2018, 108, 734–740. [Google Scholar] [CrossRef]
  49. Lu, G.; Fu, D.; Jia, C.; Chai, L.; Han, Y.; Liu, J.; Wu, T.; Xie, R.; Chang, Z.; Yang, H.; et al. Reduced miR-105-1 levels are associated with poor survival of patients with non-small cell lung cancer. Oncol. Lett. 2017, 14, 7842–7848. [Google Scholar] [CrossRef] [Green Version]
  50. Li, H.Y.; Liang, J.L.; Kuo, Y.L.; Lee, H.H.; Calkins, M.J.; Chang, H.T.; Lin, F.C.; Chen, Y.C.; Hsu, T.I.; Hsiao, M.; et al. miR-105/93-3p promotes chemoresistance and circulating miR-105/93-3p acts as a diagnostic biomarker for triple negative breast cancer. Breast Cancer Res. BCR 2017, 19, 133. [Google Scholar] [CrossRef]
  51. Xiang, Y.; Ma, N.; Wang, D.; Zhang, Y.; Zhou, J.; Wu, G.; Zhao, R.; Huang, H.; Wang, X.; Qiao, Y.; et al. MiR-152 and miR-185 co-contribute to ovarian cancer cells cisplatin sensitivity by targeting DNMT1 directly: A novel epigenetic therapy independent of decitabine. Oncogene 2014, 33, 378–386. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Dai, J.; Wei, R.; Zhang, P.; Kong, B. Overexpression of microRNA-195-5p reduces cisplatin resistance and angiogenesis in ovarian cancer by inhibiting the PSAT1-dependent GSK3beta/beta-catenin signaling pathway. J. Transl. Med. 2019, 17, 190. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Xu, C.Z.; Xie, J.; Jin, B.; Chen, X.W.; Sun, Z.F.; Wang, B.X.; Dong, P. Gene and microRNA expression reveals sensitivity to paclitaxel in laryngeal cancer cell line. Int. J. Clin. Exp. Pathol. 2013, 6, 1351–1361. [Google Scholar] [PubMed]
  54. Dai, Y.; Xie, C.H.; Neis, J.P.; Fan, C.Y.; Vural, E.; Spring, P.M. MicroRNA expression profiles of head and neck squamous cell carcinoma with docetaxel-induced multidrug resistance. Head Neck 2011, 33, 786–791. [Google Scholar] [CrossRef]
  55. Sun, W.; Ma, Y.; Chen, P.; Wang, D. MicroRNA-10a silencing reverses cisplatin resistance in the A549/cisplatin human lung cancer cell line via the transforming growth factor-beta/Smad2/STAT3/STAT5 pathway. Mol. Med. Rep. 2015, 11, 3854–3859. [Google Scholar] [CrossRef] [Green Version]
  56. Yan, M.; Yang, X.; Shen, R.; Wu, C.; Wang, H.; Ye, Q.; Yang, P.; Zhang, L.; Chen, M.; Wan, B.; et al. miR-146b promotes cell proliferation and increases chemosensitivity, but attenuates cell migration and invasion via FBXL10 in ovarian cancer. Cell Death Dis. 2018, 9, 1123. [Google Scholar] [CrossRef]
  57. Cheng, R.; Lu, C.; Zhang, G.; Zhang, G.; Zhao, G. Overexpression of miR-203 increases the sensitivity of NSCLC A549/H460 cell lines to cisplatin by targeting Dickkopf-1. Oncol. Rep. 2017, 37, 2129–2136. [Google Scholar] [CrossRef]
  58. Liu, Y.; Gao, S.; Chen, X.; Liu, M.; Mao, C.; Fang, X. Overexpression of miR-203 sensitizes paclitaxel (Taxol)-resistant colorectal cancer cells through targeting the salt-inducible kinase 2 (SIK2). Tumour Biol. J. Int. Soc. Oncodev. Biol. Med. 2016, 37, 12231–12239. [Google Scholar] [CrossRef]
  59. Mohamed, Z.; Hassan, M.K.; Okasha, S.; Mitamura, T.; Keshk, S.; Konno, Y.; Kato, T.; El-Khamisy, S.F.; Ohba, Y.; Watari, H. miR-363 confers taxane resistance in ovarian cancer by targeting the Hippo pathway member, LATS2. Oncotarget 2018, 9, 30053–30065. [Google Scholar] [CrossRef] [Green Version]
  60. Ou, Y.; Zhai, D.; Wu, N.; Li, X. Downregulation of miR-363 increases drug resistance in cisplatin-treated HepG2 by dysregulating Mcl-1. Gene 2015, 572, 116–122. [Google Scholar] [CrossRef]
  61. Cao, L.; Wan, Q.; Li, F.; Tang, C.E. MiR-363 inhibits cisplatin chemoresistance of epithelial ovarian cancer by regulating snail-induced epithelial-mesenchymal transition. BMB Rep. 2018, 51, 456–461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Yu, P.N.; Yan, M.D.; Lai, H.C.; Huang, R.L.; Chou, Y.C.; Lin, W.C.; Yeh, L.T.; Lin, Y.W. Downregulation of miR-29 contributes to cisplatin resistance of ovarian cancer cells. Int. J. Cancer 2014, 134, 542–551. [Google Scholar] [CrossRef] [PubMed]
  63. Yang, L.; Li, N.; Yan, Z.; Li, C.; Zhao, Z. MiR-29a-Mediated CD133 Expression Contributes to Cisplatin Resistance in CD133(+) Glioblastoma Stem Cells. J. Mol. Neurosci. MN 2018, 66, 369–377. [Google Scholar] [CrossRef] [PubMed]
  64. Xu, W.; Li, Z.; Zhu, X.; Xu, R.; Xu, Y. miR-29 Family Inhibits Resistance to Methotrexate and Promotes Cell Apoptosis by Targeting COL3A1 and MCL1 in Osteosarcoma. Med. Sci. Monit. Int. Med. J. Exp. Clin. Res. 2018, 24, 8812–8821. [Google Scholar] [CrossRef] [PubMed]
  65. Sorrentino, A.; Liu, C.G.; Addario, A.; Peschle, C.; Scambia, G.; Ferlini, C. Role of microRNAs in drug-resistant ovarian cancer cells. Gynecol. Oncol. 2008, 111, 478–486. [Google Scholar] [CrossRef]
  66. Guo, Y.; Zhang, H.; Xie, D.; Hu, X.; Song, R.; Zhu, L. Non-coding RNA NEAT1/miR-214-3p contribute to doxorubicin resistance of urothelial bladder cancer preliminary through the Wnt/beta-catenin pathway. Cancer Manag. Res. 2018, 10, 4371–4380. [Google Scholar] [CrossRef] [Green Version]
  67. Zhang, Y.; Xu, W.; Ni, P.; Li, A.; Zhou, J.; Xu, S. MiR-99a and MiR-491 Regulate Cisplatin Resistance in Human Gastric Cancer Cells by Targeting CAPNS1. Int. J. Biol. Sci. 2016, 12, 1437–1447. [Google Scholar] [CrossRef] [Green Version]
  68. Akbari Moqadam, F.; Lange-Turenhout, E.A.; Aries, I.M.; Pieters, R.; den Boer, M.L. MiR-125b, miR-100 and miR-99a co-regulate vincristine resistance in childhood acute lymphoblastic leukemia. Leuk. Res. 2013, 37, 1315–1321. [Google Scholar] [CrossRef]
  69. Zhu, X.; Li, Y.; Xie, C.; Yin, X.; Liu, Y.; Cao, Y.; Fang, Y.; Lin, X.; Xu, Y.; Xu, W.; et al. miR-145 sensitizes ovarian cancer cells to paclitaxel by targeting Sp1 and Cdk6. Int. J. Cancer 2014, 135, 1286–1296. [Google Scholar] [CrossRef]
  70. Khordadmehr, M.; Shahbazi, R.; Sadreddini, S.; Baradaran, B. miR-193: A new weapon against cancer. J. Cell. Physiol. 2019, 234, 16861–16872. [Google Scholar] [CrossRef]
  71. Janikova, M.; Zizkova, V.; Skarda, J.; Kharaishvili, G.; Radova, L.; Kolar, Z. Prognostic significance of miR-23b in combination with P-gp, MRP and LRP/MVP expression in non-small cell lung cancer. Neoplasma 2016, 63, 576–587. [Google Scholar] [CrossRef] [PubMed]
  72. Xu, Y.; Han, Y.F.; Ye, B.; Zhang, Y.L.; Dong, J.D.; Zhu, S.J.; Chen, J. miR-27b-3p is Involved in Doxorubicin Resistance of Human Anaplastic Thyroid Cancer Cells via Targeting Peroxisome Proliferator-Activated Receptor Gamma. Basic Clin. Pharmacol. Toxicol. 2018, 123, 670–677. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Husted, S.; Sokilde, R.; Rask, L.; Cirera, S.; Busk, P.K.; Eriksen, J.; Litman, T. MicroRNA expression profiles associated with development of drug resistance in Ehrlich ascites tumor cells. Mol. Pharm. 2011, 8, 2055–2062. [Google Scholar] [CrossRef] [PubMed]
  74. Xu, X.; Jiang, X.; Chen, L.; Zhao, Y.; Huang, Z.; Zhou, H.; Shi, M. MiR-181a Promotes Apoptosis and Reduces Cisplatin Resistance by Inhibiting Osteopontin in Cervical Cancer Cells. Cancer Biother. Radiopharm. 2019, 34, 559–565. [Google Scholar] [CrossRef]
  75. Wang, X.; Chen, X.; Meng, Q.; Jing, H.; Lu, H.; Yang, Y.; Cai, L.; Zhao, Y. MiR-181b regulates cisplatin chemosensitivity and metastasis by targeting TGFbetaR1/Smad signaling pathway in NSCLC. Sci. Rep. 2015, 5, 17618. [Google Scholar] [CrossRef] [Green Version]
  76. Cheng, W.; Liu, T.; Wan, X.; Gao, Y.; Wang, H. MicroRNA-199a targets CD44 to suppress the tumorigenicity and multidrug resistance of ovarian cancer-initiating cells. FEBS J. 2012, 279, 2047–2059. [Google Scholar] [CrossRef]
  77. Liu, M.X.; Siu, M.K.; Liu, S.S.; Yam, J.W.; Ngan, H.Y.; Chan, D.W. Epigenetic silencing of microRNA-199b-5p is associated with acquired chemoresistance via activation of JAG1-Notch1 signaling in ovarian cancer. Oncotarget 2014, 5, 944–958. [Google Scholar] [CrossRef] [Green Version]
  78. Januchowski, R.; Zawierucha, P.; Rucinski, M.; Andrzejewska, M.; Wojtowicz, K.; Nowicki, M.; Zabel, M. Drug transporter expression profiling in chemoresistant variants of the A2780 ovarian cancer cell line. Biomed. Pharmacother. 2014, 68, 447–453. [Google Scholar] [CrossRef]
  79. Januchowski, R.; Zawierucha, P.; Rucinski, M.; Nowicki, M.; Zabel, M. Extracellular matrix proteins expression profiling in chemoresistant variants of the A2780 ovarian cancer cell line. BioMed Res. Int. 2014, 2014, 365867. [Google Scholar] [CrossRef] [Green Version]
  80. Prislei, S.; Mozzetti, S.; Filippetti, F.; De Donato, M.; Raspaglio, G.; Cicchillitti, L.; Scambia, G.; Ferlini, C. From plasma membrane to cytoskeleton: A novel function for semaphorin 6A. Mol. Cancer Ther. 2008, 7, 233–241. [Google Scholar] [CrossRef] [Green Version]
  81. Januchowski, R.; Wojtowicz, K.; Sujka-Kordowska, P.; Andrzejewska, M.; Zabel, M. MDR gene expression analysis of six drug-resistant ovarian cancer cell lines. BioMed Res. Int. 2013, 2013, 241763. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Huang, H.; Chen, J.; Ding, C.M.; Jin, X.; Jia, Z.M.; Peng, J. LncRNA NR2F1-AS1 regulates hepatocellular carcinoma oxaliplatin resistance by targeting ABCC1 via miR-363. J. Cell. Mol. Med. 2018, 22, 3238–3245. [Google Scholar] [CrossRef] [PubMed]
  83. Januchowski, R.; Swierczewska, M.; Sterzynska, K.; Wojtowicz, K.; Nowicki, M.; Zabel, M. Increased Expression of Several Collagen Genes is Associated with Drug Resistance in Ovarian Cancer Cell Lines. J. Cancer 2016, 7, 1295–1310. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Swierczewska, M.; Sterzynska, K.; Wojtowicz, K.; Kazmierczak, D.; Izycki, D.; Nowicki, M.; Zabel, M.; Januchowski, R. PTPRK Expression Is Downregulated in Drug Resistant Ovarian Cancer Cell Lines, and Especially in ALDH1A1 Positive CSCs-Like Populations. Int. J. Mol. Sci. 2019, 20, E2053. [Google Scholar] [CrossRef] [Green Version]
  85. Zou, J.; Yin, F.; Wang, Q.; Zhang, W.; Li, L. Analysis of microarray-identified genes and microRNAs associated with drug resistance in ovarian cancer. Int. J. Clin. Exp. Pathol. 2015, 8, 6847–6858. [Google Scholar]
  86. Zhou, R. The Eph family receptors and ligands. Pharmacol. Ther. 1998, 77, 151–181. [Google Scholar] [CrossRef]
  87. Zhang, S.J.; Zhang, G.; Zhao, Y.F.; Wu, Y.; Li, J.; Chai, Y.X. [Expression of EphA7 protein in primary hepatocellular carcinoma and its clinical significance]. Zhonghua wai ke za zhi (Chine. J. Surg.) 2010, 48, 53–56. [Google Scholar]
  88. Giani, F.; Russo, G.; Pennisi, M.; Sciacca, L.; Frasca, F.; Pappalardo, F. Computational modeling reveals MAP3K8 as mediator of resistance to vemurafenib in thyroid cancer stem cells. Bioinformatics (Oxf. Engl.) 2019, 35, 2267–2275. [Google Scholar] [CrossRef]
  89. Izycka, N.; Sterzynska, K.; Januchowski, R.; Nowak-Markwitz, E. Semaphorin 3A (SEMA3A), protocadherin 9 (PCdh9), and S100 calcium binding protein A3 (S100A3) as potential biomarkers of carcinogenesis and chemoresistance of different neoplasms, including ovarian cancer-review of literature. Ginekol. Polska 2019, 90, 223–227. [Google Scholar] [CrossRef]
  90. Tang, C.; Gao, X.; Liu, H.; Jiang, T.; Zhai, X. Decreased expression of SEMA3A is associated with poor prognosis in gastric carcinoma. Int. J. Clin. Exp. Pathol. 2014, 7, 4782–4794. [Google Scholar]
  91. Jiang, H.; Qi, L.; Wang, F.; Sun, Z.; Huang, Z.; Xi, Q. Decreased semaphorin 3A expression is associated with a poor prognosis in patients with epithelial ovarian carcinoma. Int. J. Mol. Med. 2015, 35, 1374–1380. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  92. Swierczewska, M.; Klejewski, A.; Brazert, M.; Kazmierczak, D.; Izycki, D.; Nowicki, M.; Zabel, M.; Januchowski, R. New and Old Genes Associated with Primary and Established Responses to Paclitaxel Treatment in Ovarian Cancer Cell Lines. Molecules (Basel Switz.) 2018, 23, 891. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  93. Jin, Y.; Hong, F.; Bao, Q.; Xu, Q.; Duan, R.; Zhu, Z.; Zhang, W.; Ma, C. MicroRNA-145 suppresses osteogenic differentiation of human jaw bone marrow mesenchymal stem cells partially via targeting semaphorin 3A. Connect. Tissue Res. 2019, 1–9. [Google Scholar] [CrossRef] [PubMed]
  94. Wang, F.; Xia, J.; Wang, N.; Zong, H. miR-145 inhibits proliferation and invasion of esophageal squamous cell carcinoma in part by targeting c-Myc. Onkologie 2013, 36, 754–758. [Google Scholar] [CrossRef] [PubMed]
  95. Li, J.; Li, X.; Wu, L.; Pei, M.; Li, H.; Jiang, Y. miR-145 inhibits glutamine metabolism through c-myc/GLS1 pathways in ovarian cancer cells. Cell Biol. Int. 2019, 43, 921–930. [Google Scholar] [CrossRef]
  96. Nadiminty, N.; Tummala, R.; Lou, W.; Zhu, Y.; Zhang, J.; Chen, X.; eVere White, R.W.; Kung, H.J.; Evans, C.P.; Gao, A.C. MicroRNA let-7c suppresses androgen receptor expression and activity via regulation of Myc expression in prostate cancer cells. J. Biol. Chem. 2012, 287, 1527–1537. [Google Scholar] [CrossRef] [Green Version]
  97. Frank, M.; Kemler, R. Protocadherins. Curr. Opin. Cell Biol. 2002, 14, 557–562. [Google Scholar] [CrossRef]
  98. Kim, S.Y.; Yasuda, S.; Tanaka, H.; Yamagata, K.; Kim, H. Non-clustered protocadherin. Cell Adhes. Migr. 2011, 5, 97–105. [Google Scholar] [CrossRef] [Green Version]
  99. Shi, C.; Zhang, Z. Screening of potentially crucial genes and regulatory factors involved in epithelial ovarian cancer using microarray analysis. Oncol. Lett. 2017, 14, 725–732. [Google Scholar] [CrossRef] [Green Version]
  100. Chen, Y.; Xiang, H.; Zhang, Y.; Wang, J.; Yu, G. Loss of PCDH9 is associated with the differentiation of tumor cells and metastasis and predicts poor survival in gastric cancer. Clin. Exp. Metastasis 2015, 32, 417–428. [Google Scholar] [CrossRef]
  101. Sofeu Feugaing, D.D.; Gotte, M.; Viola, M. More than matrix: The multifaceted role of decorin in cancer. Eur. J. Cell Biol. 2013, 92, 1–11. [Google Scholar] [CrossRef] [PubMed]
  102. Kwan, P.; Ding, J.; Tredget, E.E. MicroRNA 181b regulates decorin production by dermal fibroblasts and may be a potential therapy for hypertrophic scar. PLoS ONE 2015, 10, e0123054. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  103. Lucero, H.A.; Kagan, H.M. Lysyl oxidase: An oxidative enzyme and effector of cell function. Cell. Mol. Life Sci. CMLS 2006, 63, 2304–2316. [Google Scholar] [CrossRef]
  104. Erler, J.T.; Bennewith, K.L.; Nicolau, M.; Dornhofer, N.; Kong, C.; Le, Q.T.; Chi, J.T.; Jeffrey, S.S.; Giaccia, A.J. Lysyl oxidase is essential for hypoxia-induced metastasis. Nature 2006, 440, 1222–1226. [Google Scholar] [CrossRef] [PubMed]
  105. Ji, F.; Wang, Y.; Qiu, L.; Li, S.; Zhu, J.; Liang, Z.; Wan, Y.; Di, W. Hypoxia inducible factor 1alpha-mediated LOX expression correlates with migration and invasion in epithelial ovarian cancer. Int. J. Oncol. 2013, 42, 1578–1588. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  106. Ween, M.P.; Oehler, M.K.; Ricciardelli, C. Transforming growth Factor-Beta-Induced Protein (TGFBI)/(betaig-H3): A matrix protein with dual functions in ovarian cancer. Int. J. Mol. Sci. 2012, 13, 10461–10477. [Google Scholar] [CrossRef] [Green Version]
  107. Karlan, B.Y.; Dering, J.; Walsh, C.; Orsulic, S.; Lester, J.; Anderson, L.A.; Ginther, C.L.; Fejzo, M.; Slamon, D. POSTN/TGFBI-associated stromal signature predicts poor prognosis in serous epithelial ovarian cancer. Gynecol. Oncol. 2014, 132, 334–342. [Google Scholar] [CrossRef]
  108. Bhushan, R.; Grunhagen, J.; Becker, J.; Robinson, P.N.; Ott, C.E.; Knaus, P. miR-181a promotes osteoblastic differentiation through repression of TGF-beta signaling molecules. Int. J. Biochem. Cell Biol. 2013, 45, 696–705. [Google Scholar] [CrossRef]
  109. Young, M.F.; Kerr, J.M.; Termine, J.D.; Wewer, U.M.; Wang, M.G.; McBride, O.W.; Fisher, L.W. cDNA cloning, mRNA distribution and heterogeneity, chromosomal location, and RFLP analysis of human osteopontin (OPN). Genomics 1990, 7, 491–502. [Google Scholar] [CrossRef] [Green Version]
  110. Kim, J.H.; Skates, S.J.; Uede, T.; Wong, K.K.; Schorge, J.O.; Feltmate, C.M.; Berkowitz, R.S.; Cramer, D.W.; Mok, S.C. Osteopontin as a potential diagnostic biomarker for ovarian cancer. Jama 2002, 287, 1671–1679. [Google Scholar] [CrossRef] [Green Version]
  111. Bao, L.H.; Sakaguchi, H.; Fujimoto, J.; Tamaya, T. Osteopontin in metastatic lesions as a prognostic marker in ovarian cancers. J. Biomed. Sci. 2007, 14, 373–381. [Google Scholar] [CrossRef] [PubMed]
  112. Das, S.; Samant, R.S.; Shevde, L.A. Nonclassical activation of Hedgehog signaling enhances multidrug resistance and makes cancer cells refractory to Smoothened-targeting Hedgehog inhibition. J. Biol. Chem. 2013, 288, 11824–11833. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  113. Cui, R.; Takahashi, F.; Ohashi, R.; Yoshioka, M.; Gu, T.; Tajima, K.; Unnoura, T.; Iwakami, S.; Hirama, M.; Ishiwata, T.; et al. Osteopontin is involved in the formation of malignant pleural effusion in lung cancer. Lung Cancer (Amst. Neth.) 2009, 63, 368–374. [Google Scholar] [CrossRef] [PubMed]
  114. Huo, X.; Zhang, K.; Yi, L.; Mo, Y.; Liang, Y.; Zhao, J.; Zhang, Z.; Xu, Y.; Zhen, G. Decreased epithelial and plasma miR-181b-5p expression associates with airway eosinophilic inflammation in asthma. Clin. Exp. Allergy J. Br. Soc. Allergy Clin. Immunol. 2016, 46, 1281–1290. [Google Scholar] [CrossRef]
  115. Helleman, J.; Jansen, M.P.; Span, P.N.; van Staveren, I.L.; Massuger, L.F.; Meijer-van Gelder, M.E.; Sweep, F.C.; Ewing, P.C.; van der Burg, M.E.; Stoter, G.; et al. Molecular profiling of platinum resistant ovarian cancer. Int. J. Cancer 2006, 118, 1963–1971. [Google Scholar] [CrossRef]
  116. Stelcer, E.; Kulcenty, K.; Rucinski, M.; Jopek, K.; Richter, M.; Trzeciak, T.; Suchorska, W.M. The Role of MicroRNAs in Early Chondrogenesis of Human Induced Pluripotent Stem Cells (hiPSCs). Int. J. Mol. Sci. 2019, 20, 4371. [Google Scholar] [CrossRef] [Green Version]
  117. Kulcenty, K.; Wroblewska, J.P.; Rucinski, M.; Kozlowska, E.; Jopek, K.; Suchorska, W.M. MicroRNA Profiling During Neural Differentiation of Induced Pluripotent Stem Cells. Int. J. Mol. Sci. 2019, 20, 3651. [Google Scholar] [CrossRef] [Green Version]
  118. Gautier, L.; Cope, L.; Bolstad, B.M.; Irizarry, R.A. Affy—Analysis of Affymetrix GeneChip data at the probe level. Bioinformatics (Oxf. Engl.) 2004, 20, 307–315. [Google Scholar] [CrossRef]
  119. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43, e47. [Google Scholar] [CrossRef]
  120. Cava, C.; Colaprico, A.; Bertoli, G.; Graudenzi, A.; Silva, T.C.; Olsen, C.; Noushmehr, H.; Bontempi, G.; Mauri, G.; Castiglioni, I. SpidermiR: An R/Bioconductor Package for Integrative Analysis with miRNA Data. Int. J. Mol. Sci. 2017, 18, 274. [Google Scholar] [CrossRef] [Green Version]
  121. Deng, M.; Brägelmann, J.; Kryukov, I.; Saraiva-Agostinho, N.; Perner, S. FirebrowseR: An R client to the Broad Institute’s Firehose Pipeline. Database 2017, 6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  122. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef] [PubMed]
  123. Szyszka, M.; Paschke, L.; Tyczewska, M.; Jopek, K.; Celichowski, P.; Milecka, P.; Sultanova, G.; Stelcer, E.; Malinska, A.; Malendowicz, L.K.; et al. Analysis of transcriptome, selected intracellular signaling pathways, proliferation and apoptosis of LNCaP cells exposed to high leptin concentrations. Int. J. Mol. Sci. 2019, 20, e5412. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  124. Weinstein, J.N.; Collisson, E.A.; Mills, G.B.; Shaw, K.R.M.; Ozenberger, B.A.; Ellrott, K.; Sander, C.; Stuart, J.M.; Chang, K.; Creighton, C.J.; et al. The cancer genome atlas pan-cancer analysis project. Nat. Genet. 2013, 45, 1113–1120. [Google Scholar] [CrossRef]
Ijms 21 02619 g001 550
Figure 1. Scatter plots displaying the miRNA with expression levels that were upregulated (green dots) or downregulated (red dots) by 5-fold or more in drug resistant cells in relation to the W1 cells. Grey dots indicate miRNAs below cut-off criteria (|fold|> 5, p < 0.05).
Figure 1. Scatter plots displaying the miRNA with expression levels that were upregulated (green dots) or downregulated (red dots) by 5-fold or more in drug resistant cells in relation to the W1 cells. Grey dots indicate miRNAs below cut-off criteria (|fold|> 5, p < 0.05).
Ijms 21 02619 g001
Ijms 21 02619 g002 550
Figure 2. Expression ratios of miRNAs in drug-resistant sublines.
Figure 2. Expression ratios of miRNAs in drug-resistant sublines.
Ijms 21 02619 g002
Ijms 21 02619 g003 550
Figure 3. Regulation of selected target genes by miRNAs in the W1CR cell line.
Figure 3. Regulation of selected target genes by miRNAs in the W1CR cell line.
Ijms 21 02619 g003
Ijms 21 02619 g004 550
Figure 4. Regulation of selected target genes expression by miRNAs in the W1DR cell line.
Figure 4. Regulation of selected target genes expression by miRNAs in the W1DR cell line.
Ijms 21 02619 g004
Ijms 21 02619 g005 550
Figure 5. Regulation of selected target genes expression by miRNAs in the W1PR1 cell line.
Figure 5. Regulation of selected target genes expression by miRNAs in the W1PR1 cell line.
Ijms 21 02619 g005
Ijms 21 02619 g006 550
Figure 6. Regulation of selected target genes expression by miRNAs in the W1PR1 cell line.
Figure 6. Regulation of selected target genes expression by miRNAs in the W1PR1 cell line.
Ijms 21 02619 g006
Ijms 21 02619 g007 550
Figure 7. Negative correlation between miRNA and its appropriate target gene in all samples (all), chemotherapy sensitive (tumor free), and chemotherapy resistant patients (with tumor). All of the expression data were log2 transformed to achieve normal distribution.
Figure 7. Negative correlation between miRNA and its appropriate target gene in all samples (all), chemotherapy sensitive (tumor free), and chemotherapy resistant patients (with tumor). All of the expression data were log2 transformed to achieve normal distribution.
Ijms 21 02619 g007
Table
Table 1. List of the miRNA fold changes and false discovery rate (FDR) corrected p values (adj.p.val.). Each comparison was performed in relation to W1 cells.
Table 1. List of the miRNA fold changes and false discovery rate (FDR) corrected p values (adj.p.val.). Each comparison was performed in relation to W1 cells.
Gene NamemiRBase Accession Number Fold Change (adj.p.val)
W1PR1W1PR2W1CRW1TRW1DR
hsa-let-7cMIMAT00000645.40 (0.001)N.SN.SN.SN.S
hsa-miR-100MIMAT0000098−5.01 (0.0003)N.SN.SN.SN.S
hsa-miR-105MIMAT0000102−19.22 (4.6 × 10−7)−21.80 (1.5 × 10−6)−15.92 (1.5 × 10−5)N.SN.S
hsa-miR-10aMIMAT0000253−34.42 (2.3 × 10−7)−9.07 (0.0002)N.SN.SN.S
hsa-miR-1271 MIMAT00057969.08 (0.004)N.S6.24 (0.048)5.54 (0.049)N.S
hsa-miR-1285MIMAT0005876N.S5.54 (0.01)N.SN.SN.S
hsa-miR-143MIMAT00004378.81 (0.003)N.SN.SN.SN.S
hsa-miR-145MIMAT000043712.07 (4.8 × 10−5)N.SN.S7.05 (0.002)N.S
hsa-miR-146b-5pMIMAT0002809−39.13 (3.5 × 10−8)N.SN.SN.S−6.93 (0.0003)
hsa-miR-152 MIMAT000043818.19 (6.5 × 10−5)16.52 (0.0003)8.53 (0.0006)N.SN.S
hsa-miR-181aMIMAT0000256N.SN.SN.S−12.28 (0.0003)N.S
hsa-miR-181bMIMAT0000257N.SN.SN.S−9.39 (1.4 × 10−7)N.S
hsa-miR-18bMIMAT0001412−10.22 (0.0006)N.SN.SN.SN.S
hsa-miR-193a-5pMIMAT00046145.32 (0.01)N.SN.SN.SN.S
hsa-miR-193bMIMAT00028197.18 (0.01)N.SN.SN.SN.S
hsa-miR-195MIMAT000046113.94 (0.003)N.S6.67 (0.046)N.SN.S
hsa-miR-199a-3pMIMAT0000462N.SN.SN.SN.S−23.09 (0.01)
hsa-miR-199b-3pMIMAT0000464N.SN.SN.SN.S−28.18 (0.02)
hsa-miR-203 MIMAT000026414.45 (6.9 × 10−6)N.S22.33 (1.5 × 10−5)N.SN.S
hsa-miR-205MIMAT0000465N.SN.SN.SN.S6.23 (0.004)
hsa-miR-20bMIMAT0000466−8.87 (6 × 10−8)N.SN.SN.SN.S
hsa-miR-214MIMAT0000271N.SN.SN.SN.S−26.05 (0.001)
hsa-miR-23bMIMAT00004686.06 (0.0009)N.SN.SN.SN.S
hsa-miR-23cMIMAT00004697.47 (0.02)N.SN.SN.SN.S
hsa-miR-27bMIMAT00004709.21 (0.0003)N.SN.SN.SN.S
hsa-miR-29aMIMAT0000086N.SN.SN.S−22.04 (0.0009)N.S
hsa-miR-30aMIMAT0000087N.S14.82 (0.001)N.SN.SN.S
hsa-miR-31MIMAT0000089−217.59 (4.6 × 10−5)−51.72 (0.003)−13.51 (0.048)−36.90 (0.008)N.S
hsa-miR-335MIMAT0000765N.SN.SN.S21.31 (0.0009)N.S
hsa-miR-363MIMAT0000707−16.37 (1.55 × 10−7)N.SN.SN.SN.S
hsa-miR-4269MIMAT0016897N.SN.SN.SN.S−5.13 (0.04)
hsa-miR-449cMIMAT0016897−11.97 (0.006)−8.09 (0.04)−5.61 (0.049)−6.53 (0.047)N.S
hsa-miR-4721MIMAT0019835N.S7.18 (0.03)N.SN.SN.S
hsa-miR-497MIMAT00028208.22 (0.02)N.SN.SN.SN.S
hsa-miR-501-5pMIMAT0002872N.S−5.60 (0.002)N.SN.SN.S
hsa-miR-572MIMAT0003237N.S6.23 (0.04)N.SN.SN.S
hsa-miR-708MIMAT0004926N.SN.S5.52 (0.007)N.SN.S
hsa-miR-767-5pMIMAT0003882−20.61 (9 × 10−6)−19.99 (3 × 10−5)−13.95 (0.0004)N.SN.S
hsa-miR-875-3pMIMAT0004923−7.93 (0.02)−7.32 (0.048)N.SN.S−6.28 (0.049)
hsa-miR-99aMIMAT000009711.02 (0.004)N.SN.SN.SN.S

Share and Cite

MDPI and ACS Style

Kazmierczak, D.; Jopek, K.; Sterzynska, K.; Ginter-Matuszewska, B.; Nowicki, M.; Rucinski, M.; Januchowski, R. The Significance of MicroRNAs Expression in Regulation of Extracellular Matrix and Other Drug Resistant Genes in Drug Resistant Ovarian Cancer Cell Lines. Int. J. Mol. Sci. 2020, 21, 2619. https://doi.org/10.3390/ijms21072619

AMA Style

Kazmierczak D, Jopek K, Sterzynska K, Ginter-Matuszewska B, Nowicki M, Rucinski M, Januchowski R. The Significance of MicroRNAs Expression in Regulation of Extracellular Matrix and Other Drug Resistant Genes in Drug Resistant Ovarian Cancer Cell Lines. International Journal of Molecular Sciences. 2020; 21(7):2619. https://doi.org/10.3390/ijms21072619

Chicago/Turabian Style

Kazmierczak, Dominika, Karol Jopek, Karolina Sterzynska, Barbara Ginter-Matuszewska, Michal Nowicki, Marcin Rucinski, and Radoslaw Januchowski. 2020. "The Significance of MicroRNAs Expression in Regulation of Extracellular Matrix and Other Drug Resistant Genes in Drug Resistant Ovarian Cancer Cell Lines" International Journal of Molecular Sciences 21, no. 7: 2619. https://doi.org/10.3390/ijms21072619

APA Style

Kazmierczak, D., Jopek, K., Sterzynska, K., Ginter-Matuszewska, B., Nowicki, M., Rucinski, M., & Januchowski, R. (2020). The Significance of MicroRNAs Expression in Regulation of Extracellular Matrix and Other Drug Resistant Genes in Drug Resistant Ovarian Cancer Cell Lines. International Journal of Molecular Sciences, 21(7), 2619. https://doi.org/10.3390/ijms21072619

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop