ABSTRACT
Conclusion:
miR-34a-5p and miR-145-5p may have an impact on the breast cancer cell line MDA-MB-231 proliferation and EMT mechanism. At the same time, according to our study results, it was revealed that E2F1, MEK1 and CDK4 genes, whose expression level decreased after transfection of mimic miR-34a-5p, could be targeted by miR-34a-5p in breast cancer, and that the expression level of KLF4, which decreased as a result of mimic miR-145-5p transfection, could be the target.
Results:
The proliferation of the miR-34a-5p or miR-145-5p transfected cells decreased compared to the control groups. The expression of E2F transcription factor 1 (E2F1) (p=0.009), mitogen activated protein kinase 1 (MEK1) (p=0.001) and cyclin dependent kinase 4 (CDK4) (p=0.005) genes, which were among the genes targeted by miR-34a-5p, were significantly reduced. EMT markers were significantly changed in miR-34a-5p transfected cells (E-Cad increase p=0.01; Vimentin decrease p=0.008). Kruppel-like factor 4 (KLF4) (p=0.007) targeted miR-145-5p were significantly reduced and EMT markers were significantly changed in miR-145-5p transfected cells (E-Cad increase p=0.0005; Vimentin decrease p=0.006).
Materials and Methods:
Firstly, literature search and in silico analysis were performed to detect possible target genes of miR-34a-5p and miR-145-5p, which are known to be tumor suppressors. Mimic miR-34a-5p and miR-145-5p were transfected to the breast cancer cell line MDA-MB-231. Deregulated genes were investigated by the quantitative real-time polymerase chain reaction compared to control cells. Also, the effect of these miRNAs on proliferation was determined using the Water Soluble Tetrazolium Salt-8 method. Finally, the expressions of epithelial mesenchymal transition (EMT) markers, which are known to be important in the metastatic process, are examined.
Aim:
It was aimed to investigate the effect of miR-34a-5p and miR-145-5p on breast cancer cell line MDA-MB-231 and to determine the expression of target genes of these microRNAs (miRNAs).
INTRODUCTION
Breast cancer is the most common type of cancer in women, and it occupies a large place among cancer-related deaths in women all over the world, especially among women in underdeveloped countries1. MicroRNAs (miRNAs) are single-stranded, non-protein-coding RNAs containing approximately 19-24 nucleotides and repress their target mRNAs. Depending on the functions of target genes, miRNAs are biomolecules that play a vital role in various cellular processes such as differentiation, proliferation and apoptosis, and their deregulation is effective in the development of many types of cancer2. The miR-34 family, which was characterized in 2007, has three members, miR-34a-5p, miR-34b and miR-34c. The miR-34 family of tumor suppressor miRNAs is known to target many cancer-related genes, including p53. It has been reported that these miRNA family members are deregulated in a large number of cancer types3. One of the other most important tumor suppressor miRNA is miR-145-5p. The expression level of miR-145-5p has been shown to be frequently decreased in many cancers, including colorectal cancer, hematological cancers, prostate cancer, ovarian cancer, B-cell tumor, and breast cancer4. Undoubtedly, it is extremely important in terms of the molecular mechanism of cancer to learn that this decrease in the expression levels of miR-34a-5p and miR-145-5p, whose expression levels are found to be decreased in many cancers, contribute to the cancerization process by increasing their expression of target genes4,5. However, studies on miR-34a-5p and miR-145-5p miRNAs targeting which genes in breast cancer and through which pathways are effective in the cancer process are limited.
Within the scope of the study, it was aimed to investigate the effects of related miRNAs on cell proliferation by transfection of mimic miR-145-5p and miR-34a-5p to MDA-MB-231 cells, which are triple negative breast cancer cell line.
Afterwards, the expression levels of one of the target genes of miR-34a-5p, Mouse double minute 4 (MDM4), E2F1, Insulin-like growth factor 1 receptor (IGF1R), MEK1, CDK4, cluster of differentiation 44 (CD44), eukaryotic translation initiation factor 4E (EIF4E) targeted by neurogenic locus notch homolog protein 2 (NOTCH2) and target genes of miR-145-5p, CDK4, Cyclin-dependent kinase 6 (CDK6), Epidermal growth factor receptor (EGFR), SRY-Box Transcription Factor (SOX2), KLF4, jagged canonical notch ligand 1 (JAG1) genes were investigated and it was aimed to investigate whether the expression levels of related genes are correlated with miR-34a-5p and miR-145-5p. Finally, it was aimed to investigate the relationship of miR-34a-5p and miR-145-5p with metastasis by looking at the expression of epithelial mesenchymal transition (EMT) markers that are known to contribute to the metastatic process.
MATERIALS AND METHODS
Cell Culture and miRNA Transfection
For the functional study of miR-34a-5p and miR-145-5p non-coding RNAs in vitro, MDA-MB-231 breast cancer cell line 1% antibiotic (penicillin/streptomycin) and 10% fetal bovine serum (Gibco, Thermo Fisher Scientific Inc., Waltham, MA, USA) was cultured and replicated in RPMI-1640 (Roswell Park Memorial Institute-1640) (Gibco, Thermo Fisher Scientific, Waltham, MA, USA) medium at 37 °C in a 5% CO2 incubator.
In our study, in which an in vitro functional experiment of miR-34a-5p and miR-145-5p synthetic mimic (Thermo Fisher Scientific Inc.) was planned separately in replicated cells, an oligonucleotide sequence control miRNA (non-targeting control miRNA, control nt miRNA) (Thermo Fisher Scientific Inc.) that did not target any gene was used (miR-145-5p 5’-GUCCAGUUUUCCCAGGAAUCCCU-3’), non-targeting control miRNA (5’-UCACAACCUCCUAGAAAGAGUAGA-3’), miR-34a-5p (5’-UGGCAGUGUCUUAGCUGGUUGU-3’).
24 hours after seeding the appropriate number of cells (well/300,000 cells) on 6-well plates, transfection of miRNA mimic was performed in the cells. Mimics were transported into the cell with Lipofectamine 2000 reagent (Thermo Fisher Scientific Inc.) by applying the manufacturer’s protocol in transfection. Then, the cells were incubated at 37 °C and 24 hours after transfection, cells transfected with miR-145-5p, miR-34a-5p and cells transfected with non-targeting control miRNA were taken separately in 1000 µL Trizol (MRC Inc.) in eppendorfs and RNA was then brought to -80 °C for isolation.
Isolation of RNA from Cells After Transfection
Total RNA was isolated from cells taken in trizol (MRC Inc.) in accordance with the manufacturer’s protocol, and RNA concentration and purity were measured with the help of NanoDrop ND-2000c (Thermo Fisher Scientific Inc.) spectrophotometer.
Confirmation of miRNA cDNA Synthesis and Transfection by miRNA qRT-PCR
Using TaqMan probes, quantitative real-time polymerase chain reaction (qRT-PCR) was performed to determine the expression levels of miR-34a-5p and miR-145-5p in cell lines and to check whether the transfection processes were effective. MiRNA cDNA synthesis using total 30 ng RNA was performed using TaqMan MicroRNA reverse transcriptase kit (Applied Bio., Foster City, CA, United States) and miRNA RT primers according to the manufacturer’s protocol. TaqMan assays were used for qRT-PCR, and qRT-PCR was performed on LightCycler® 480 (Roche) instrument using TaqMan Universal Master Mix (Thermo Fisher Scientific Inc.) kit with TaqMan miRNA probes (Thermo Fisher Scientific Inc.) and RNU43 (control miRNA) probes (Thermo Fisher Scientific Inc.) according to the manufacturer’s protocol. Experiments were performed in duplicate and 2-ΔΔCt method was used for relative quantitation analysis.
In Silico Prediction of Target Genes for Expression of miR-34a-5p and miR-145-5p After Transfection
Target prediction tools for prediction of related genes which are the tools such as miRBase (http://www.mirbase.org/), miRWalk (http://mirwalk.umm.uni-heidelberg.de/) and miRTargetLink (https://ccb-web.cs.uni-saarland.de/mirtargetlink/) and databases such as Genecards, NCBI, PubMed were used. That the relationship of the gene selected to be investigated selectively in the prediction phase of possible target genes and the related miRNA has not been shown in the MDA-MB-231 breast cancer cell line before, that the expression level of the selected gene in breast cancer has been shown to increase in breast cancer previously in the literature, and that the selected gene is targeted by the relevant miRNA at least in two different cancers were taken into account. The miRTargetLink database was also used to schematize genes with strong interactions of miR-34a-5p and miR-145-5p6.
cDNA Synthesis and qRT-PCR Procedures for the Target Genes of miR-34a-5p and miR-145-5p
The cDNA extraction with SCRIPT cDNA synthesis kit (Jena Bioscience) and qPCR SybrMaster (Jena Bioscience) and qRT-PCR processes were performed in accordance with the manufacturer’s protocol, using a total of 1000 ng of RNA obtained from the study group and control group cells.
While the expression levels of MDM4, E2F1, IGF1R, MEK1, CDK4, CD44 and NOTCH2 genes were determined in miR-34a-5p transfected cells, the expression levels of EIF4E, EGFR, SOX2, CDK4, CDK6, JAG1 and KLF4 genes were examined in miR-145-5p transfected cells. The primer sequences of the examined genes are shown in Table 1. qRT-PCR processes were performed on a LightCycler® 480 (Roche) device.
Investigation of the Effect of miR-34a-5p and miR-145-5p on Cell Proliferation
After cells were transfected with miR-34a-5p and miR-145-5p synthetic mimic and control miRNA (non-targeting control miRNA), whether there was a change in the proliferation of cells were examined using the Cell Proliferation Reagent Water Soluble Tetrazolium Salt-8 (WST-8) (EcoTech Biotechnology) according to the manufacturer’s protocol.
Cells seeded in triplicate at 5000 cells per well on 96-well plates were cultured at 37 °C in a 5% CO2 incubator. Proliferation changes in cells transfected after 24 hours were measured at the following 24th, 48th and 72nd hours. For measurement, 10 µL of WST-8 (EcoTech Biotechnology) solution was added to each well and cells were incubated at 37 °C for three hours in the dark. Cell viability was evaluated using a plate reader (Thermo Varioskan flash) device capable of measuring absorbance at 450 nm.
Investigation of the Effect of miR-34a-5p and miR-145-5p on Cell Migration with EMT Marker
In order to examine whether the transfected miRNAs contribute to the migration and invasion potential of cells, the expression levels of E-cadherin, N-cadherin, Vimentin genes, which are important EMT markers in epithelial mesenchymal transition, were evaluated by qRT-PCR technique. qRT-PCR processes were performed using qPCR SybrMaster (Jena Bioscience) in accordance with the manufacturer’s protocol.
Statistical Analysis
For normalization of gene expression, β-actin primers were used as internal control and experiments were performed in duplicate. The 2-ΔΔCt method was used for the relative quantitation analysis. Statistical analyzes were performed using Student’s t-test and data were presented as mean±standard deviation. Data with a p value less than 0.05 were considered statistically significant. GraphPad Prism 5 and Statistical Package for the Social Sciences 21 programs were used for graphic drawings.
RESULTS
Verification of Transfection Procedures
It was confirmed that the transfection occurred significantly for both miRNAs in miR-34a-5p and miR-145-5p transfected cells and control cells (Figure 1).
Effect of miR-34a-5p and miR-145-5p on Expression Level of Possible Target Genes in MDA-MB-231 Cells
As a result of the literature review and in silico studies, the expression of E2F1, MEK1 and CDK4 genes, which are among the genes that are predicted to be strongly associated with miR-34a-5p, were significantly suppressed in miR-34a-5p transfected cells compared to the control group (respectively p=0.009; p=0.001; p=0.005) and KLF4 gene, which is among the genes predicted to be strongly associated with miR-145-5p, was significantly suppressed in miR-145-5p transfected cells compared to the control group (p=0.007) (Figure 2, 3, 4, 5).
Effect of miR-34a-5p and miR-145-5p on Proliferation and Migration of MDA-MB-231 Cells
It was determined that proliferation in MDA-MB-231 cells transfected with both miR-34a-5p and miR-145-5p decreased in correlation with the literature (Figures 6, 7).
As a result of expression evaluation of epithelial-mesenchymal cell markers with qRT-PCR to understand whether selected miRNAs have any regulatory role in cell migration, it was determined that E-Cad expression level was increased significantly (p=0.01; p=0.0005 respectively), Vimentin expression level was decreased significantly (p=0.008;p=0.006, respectively), and N-Cad expression level was not changed significantly (p>0.05) (Figure 8, 9) in cells transfected with miR-34a-5p and miR-145-5p compared to control cells.
As a result, it was determined that miR-34a-5p and miR-145-5p contributed to the suppression of cell migration in MDA-MB-231 breast cancer cells by regulating through EMT markers.
DISCUSSION
MiRNAs, which can act as oncogene (onco-mir) or tumor suppressor (TSmir) according to the genes they target, are among the most interesting research topics in recent years19. Abnormal expression changes in miRNAs may be involved in the cancer process by changing the expression of many genes targeted by the relevant miRNAs. The fact that miR-34 was the first miRNA to be tested in phase 1 in relation to the use of miRNAs in therapy indicates that miR-34 may be one of the most promising miRNAs in the future20.
Studies with miR-34a-5p localized at 1p36.22 have shown that ectopic expression of miR-34a-5p directly or indirectly causes a decrease in the expression of many proto-oncogenes, resulting in decreased cancer cell proliferation and elimination of metastasis5. The fact that p53, which is shown as the guardian of the genome, and its homologous p63 and p73 genes are controlled by miR-145-5p, and many important genes are targeted by this miRNA, indicates that miR-145-5p can have critical importance in early diagnosis and possible miRNA-based therapies4. Determining which genes control various cellular processes by regulating miRNAs is extremely important for understanding diseases at the molecular level21. miRNAs are involved in regulating the expression of most protein-coding genes at the post-transcriptional stage. Therefore, many bioinformatics tools have been created to predict the interaction between miRNA and target mRNA. Supporting in vitro and in vivo predictions based on base pairing between miRNA and mRNA in silico may be important for the development and progression of many diseases, especially cancer, and the development of biomarkers for these diseases2. Figure 2 and Figure 4 show that miR-34a-5p and miR-145-5p can target hundreds of genes in silico. In the literature, it has been shown that the expression levels of miR-34a-5p and miR-145-5p are considerably reduced in breast cancer patient tissue samples and MDA-MB-231 cells, as in many other cancer types22,23. Therefore, according to the results of our transfection study to investigate the effects of miR-34a-5p and miR-145-5p on MDA-MB-231 cells, miR-34a-5p and miR-145-5p have been found to have a significant effect on MDA-MB-231 cell proliferation. The effects of miR-34a-5p and miR-145-5p on cell proliferation may be related to decreased expression of various genes targeted by these miRNAs. The KLF4 gene is a gene that has been found to be highly expressed in more than 70% of breast cancers and has been reported to function as an oncogene24. The E2F1 gene is also a gene with increased expression in breast cancer and associated with breast cancer metastasis25. MEK1 and CDK4 genes, whose expression levels have been found to be increased in many cancers, including breast cancer, are important oncogenes that play critical roles in apoptosis, migration, invasion and cell proliferation processes26,27. miR-34a-5p of the E2F1, MEK1 and CDK4 genes; despite the fact that KLF4 has been shown to be effective in various cancer processes by being targeted by miR-145-5p, there are no studies in the literature on this subject in breast cancer28-32. In our study to fill this gap in the literature, it was determined that E2F1, MEK1, CDK4 and KLF4 gene expression levels were significantly decreased in MDA-MB-231 cells transfected with mimic miR-34a-5p and mimic miR-145-5p.
Study Limitations
Considering our study results and the information in the literature together, it can be predicted that miR-34a-5p may contribute to the cancer process by targeting E2F1, MEK1 and CDK4 genes, and miR-145-5p targeting KLF4 in breast cancer. Among the limitations of the study; the fact that the research was carried out in cell line, it was not confirmed in the tissue, and the genes whose expression level was examined were examined only at the mRNA level. It is obvious that more meaningful results can be obtained if the findings obtained in our study at the level of RNA expression are supported by methods such as Western-Blot and luciferase activity determination in new studies on the subject.
It was seen that CD4433, NOTCH234, MDM435, IGF1R36 genes, which have been shown as targets of miR-34a-5p both in silico and in vivo, and similarly the expression levels of EIF4E37, EGFR38, SOX221, CDK437, CDK639, JAG140 genes which have been shown as targets of miR-145-5p did not have a statistically significant change in MDA-MB-231 cells in which transfection of miRNA mimics was carried out in our study.
CONCLUSION
The data obtained from our study may be evidence that the miRNA-mRNA pairing created by theoretical prediction or the miRNA-mRNA relationship shown in vivo media in different cancer types may not be seen likewise in all cancer types due to the complex nature of cancer.