The expression of miRNA-152-3p and miRNA-185 in tumor tissues versus margin tissues of patients with chemo-treated breast cancer

Objective Breast cancer (BC) is the most significant and lethal type of cancer in women. Although there are many newly develop chemotherapy drugs for patients with BC treating at various stages, drug resistance is the most important obstacle in their effectiveness for BC treatment. On the other hand, microRNAs are considered key regulators of genes involved in carcinogenesis and chemoresistance in cancers. The purpose of this study was to evaluate the role of miR-152-3p and miR-185 in intrinsic chemoresistance and proliferation of BC. In addition, the potential role of these miRNAs during chemoresistance was evaluated through possible signaling pathways. Results Here, miR-152-3p was significantly downregulated in tumor tissues compared to the corresponding margin tissues in patients with BC (p-value ≥ 0.04407 and fold change = − 2.0552). In contrast, no statistically significant difference was observed in the miR-185 expression between the two groups. Furthermore, no significant correlation was found between the expression of these two miRNAs and subfactors, including cancer family history, abortion, and age. Downregulation of miR-152-3p could be considered a promising regulator of BC chemoresistance. Supplementary Information The online version contains supplementary material available at 10.1186/s13104-021-05647-z.


Introduction
Breast cancer (BC) is the most common malignancy among women worldwide and the main cause of cancer-related mortality [1][2][3]. According to the world health organization (WHO) reports, the incidence rates are approximately 1.7 million new cases and more than 520,000 annual deaths [4].
During the past few decades, new treatment approaches for BC have significantly developed, which not only decreased the mortality rate but also improved patients' quality of life [5][6][7]. Recent studies revealed several intricate mechanisms that are involved in intrinsic and acquired chemoresistance [8][9][10]. There are two forms of drug resistance in cancer: intrinsic resistance and acquired resistance. Intrinsic resistance refers to resistance prior to prescribing the drug, and this is due to the innate ability of cancer cells to survive in drugrelated clinical concentrations. Poor initial response to treatment indicates that a patient has intrinsic drug resistance. According to Verheul and Pinedo, more than 50% of all cancer patients are resistant to chemotherapy before starting treatment (intrinsic resistance), and other patients also resist during treatment (acquired resistance) [11]. Intrinsic and acquired resistance to Safi et al. BMC Res Notes (2021) 14:234 chemotherapeutic agents is considered the main obstacle to successful treatment in patients with BC [12]. Thus, investigation of the molecular mechanisms of drug resistance in BC may help reach better-targeted treatment by significantly improving therapeutic outcomes and prognosis for patients with chemo-resistant BC.
Increasing evidence shows that aberrant activation of the phosphatidylinositol-3-kinase (PI3K)/AKT signaling pathway is common in various types of neoplasms. Furthermore, PTEN plays a negative regulatory role in the PI3K/AKT signaling. It is involved in many important processes such as cell growth and apoptosis and drug resistance/sensitivity. Thus, targeting PI3K/AKT and its downstream pathways may provide new hope for therapeutic interventions of cancers [13][14][15][16]. MicroR-NAs (miRNAs) are a family of noncoding-RNAs, playing a crucial role in post-transcriptionally gene expression modulation [17]. Recently, it has been determined that miRNAs can play key roles during cancer development and resistance to chemotherapy by acting as oncogenes or tumor suppressors [18][19][20].
In the current study, the expression level of miR-152-3p and miR-185 was evaluated in tumor and margin tissues of patients with BC. Moreover, in combination with bioinformatics analysis, the potential role of these miRNAs during chemoresistance was evaluated through possible signaling pathways.

Tissue specimens
Two hundred tumor tissues and corresponding normal margin of tumor tissues were collected from patients with BC from 2018 to 2019 who underwent surgery at Noor Nejat Hospital, Tabriz, Iran. The collected tissues were dissected by a pathologist for diagnosing BC foci as well as margin tissues. After chemotherapy courses, it was determined-by the oncologist-that eighty tumor tissues have shown initial resistance to chemotherapy. This study was approved by the Ethics Committee of Clinical Research of Tabriz University of Medical Sciences (Ethical code: IR.TBZMED.REC.1397.1048). Demographic and clinical characteristics were collected from the patients' records (Additional file 1: 1). The specimens were obtained from mastectomy immediately frozen and stored at − 80 °C before RNA isolation.

RNA isolation
Tissue samples were cut in 20 μm-thick sections then were homogenized in liquid nitrogen. Each of the samples was directly immersed in TRIzol reagent (GeneAll, South Korea). The purity of extracted total RNAs was evaluated by the NanoDrop spectrometer (Thermo Scientific, USA) and gel electrophoresis. Finally, isolated RNAs were eluted in 30 μL of RNase-free water and stored at − 80 °C until use.

cDNA synthesis and Real-time PCR
Complementary DNA (cDNA) was synthesized using reverse-transcriptase enzyme (Thermo Fisher, US), buffer (Thermo Fisher, US), dNTP, and miRNA-specific stem-loop-primers (Sinacolon Co. Iran, Tehran) for miR-185-5p and miR-152-3p. The U6 snRNA was also considered the reference gene. This mix of cDNA synthesis was carried out by BIO-RAD-Gradient thermocycler PCR (Germany) at the temperature program of 30 min at 16 °C, 30 min at 42 °C and the last step to inactivate the enzyme the reaction was incubated at 5 min at 75 °C. Ampliqon SYBR ® Green master mix (Denmark) and miR-185-5p, miR-152-5p, and U6 specific forward primers and common reverse primers were used for quantitative Real-time PCR (qRT-PCR). All qPCR reactions were performed in triplicate. Real-time-PCR was done by LightCycler ® 96 Real-time PCR Cycler (Roche) system, in two steps as follows: for miR-185-5p: 10 Table 1.

Statistical analysis Data analysis
Shapiro-Wilk test was performed to assess the normality of data. Next, the gene expression differences were compared between tumor tissues and margins using paired t-test and delta-CTs. Two-sample paired t-test was performed to compare the miRNA expression levels with subsets of the study, including cancer family history, abortion history, and age. The gene expression levels in subsets of the study, including abortion history, age, and cancer family were compared. A t-test for the logarithm2 of fold changes was applied to compare gene expression levels among subsets of the study. Statistical analyses were done by R software version 0.5.1, and P-value of ≤ 0.05 was considered statistically significant. The results were presented as mean ± standard error of mean (SEM).

In silico analysis of miR-152-3p targets and related molecular signaling pathways
To validate targetome of miR-152-3p, miRTarBase 7.0 [29] database was applied. Further, in silico analysis was employed to achieve the vision of biological procedures that are controlled by miR-152. All target genes of miR-152-3p experimentally verified were extracted from the miRTarBase 7.0 database. Finally, to enrich related signaling pathways to miR-152-3p dysregulation, target genes were investigated in the WEB-based Gene Set Analysis Toolkit (webgestalt) database (http:// www. webge stalt. org/), and the Over-Representation enrichment (ORA) was selected and pathways were enriched based on Reactome database. The reference gene list was set to "genome-protein-coding," and the adjustment p-value was set as Benjamini-Hochberg. Other parameters were set as defaults.

In silico analysis
The 158 experimentally validated gene targets of miR-152-3p were obtained by the miRTarBase 7.0 database. The miRNA-gene interaction was depicted by Cytoscape V 3.7.1 (Additional file 1: 4). ORA database revealed that ten categories of the biological processes were enriched among miR-152-3p target genes. All 10 of these biological processes are shown in (Additional file 1: 5). Moreover, the ORA analysis discovered ten significant signaling pathways in which the targets of miR-152 are involved ( Table 2).

Discussion
Therapeutic resistance to chemotherapy is considered the most important global health challenge for patients with cancers, especially patients with BC [6,30,31]. Recently, numerous experimental and clinical studies have shown that miRNAs participate in several cellular events, particularly in resistance/sensitivity to drugs and cancer progression [32,33]. Therefore, identifying  drug resistance-specific miRNAs and their target genes is critical for understanding their role in BC chemotherapy. Further, they could be potential candidates for therapeutic intervention [34,35].
In the current study, based on previous studies and in silico analysis, it was predicted that miR-152-3p and miR-185 could be involved in BC chemoresistance. For this reason, we evaluate and compare the expression of the aforementioned miRNAs in tumor tissues relative to the margins from BC-resistant patients to chemotherapy. As mentioned in the results, we found that the expression level of miR-152-3p was significantly downregulated in BC tumor tissues compared to margin tissues. According to the previous studies, concerning chemoresistance   [37]. Shuke Ge et al. found that miR-152 plays a tumorsuppressive role in BC through negative regulation of PIK3CA expression and AKT and RPS6 inhibition, which leads to suppression of BC cell proliferation [38].
Moreover, Xu Chen et al. [26] found that the miR-148/152 family attenuates Adriamycin resistance of BC cells and tissues by downregulating the SPIN1, which is a protein highly expressed in human cancers, especially BC. [39,40]. Therefore, downregulation of miR-152-3p in tumors relative to adjacent tissues may result in an elevated expression level of SPIN1 and resistance to chemotherapy.
According to bioinformatics results, among important molecular pathways enriched by the webgestalt database, PI3K-Akt signaling pathway and regulation of PTEN mRNA translation are important enriched signaling pathways in which the targets of miR-152-3p are involved, and there may be a relationship between miR-152-3p expression and BC chemoresistance. The PI3K-Akt signaling is an intracellular signaling pathway involved in pivotal processes such as proliferation, cell survival, and angiogenesis [41]. PTEN acts as a tumor-suppressor gene and negatively regulates PI3K-Akt signaling [42]. In the upstream of Akt, PTEN blocks the formation of phosphatidylinositol-3, 4, 5-trisphosphate (PIP3) from phosphatidylinositol-4, 5-bisphosphate (PIP2), thereby inhibiting PI3 kinase (PI3K) activity. Several human tumor cell lines have been reported to evade apoptosis through the excessive activation of the PI3K/Akt pathway, which is the result of a mutation or downregulation of PTEN [43]. Additionally, it has been demonstrated that constitutive activation of PI3K-Akt signaling causes cell resistance to many chemotherapy agents by promoting proliferation and inhibiting apoptosis of cancer cells [44]. M Alam et al. showed the significant loss of PTEN expression in 26.4% BC cases, suggesting loss of PTEN expression could play a key role in breast carcinogenesis, due to lack of control of the signaling pathways such as likely PI3K-Akt signaling that mediates cellular processes like apoptosis and migration [45].
Finally, demographic characteristics of the current study showed that no significant difference was seen between the miR-152 and miR-185 expression levels. Some demographic factors such as cancer, family history, or abortion history need further studies to be clarified.

Conclusion
The findings of the current study suggested miR-152-3p as a potential biomarker in BC chemoresistance patients. Besides, the PI3K/Akt signaling was predicted as a possible key modulator during BC chemoresistance. Further studies are needed to elucidate these findings.

Limitations
The small sample size was considered the limitation of this study. Also, applying only one miRNA-target database could be mentioned as another limitation.