The pluripotency factor, OCT4 gene is a stemness marker that is involved in the tumorigenicity of different cancer types and knowing about molecular mechanisms of its regulation is crucially important. To date, a few microRNAs (miRNAs) are known to be regulators of OCT4 gene expression. Looking for the novel miRNAs which are capable of regulating OCT4 gene expression, our bioinformatics analysis introduced hsa-miR-3658 (miR-3658) as a bona fide candidate. Then, RT-qPCR results indicated that miR-3658 expression is decreased in colorectal cancer (CRC) tumor tissues, compared with normal pairs. Furthermore, RT-qPCR and western blot analysis showed that the OCT4 gene has been down-regulated following the miR-3658 overexpression. Consistently, dual-luciferase assay supported the direct interaction of miR-3658 with the 3′-UTR sequence of OCT4 gene. Unlike in HCT116 cells, overexpression of miR-3658 in SW480 cells brought about growth inhibition, cell cycle arrest and reduced cell migration, detected by flow cytometry, and scratch test assay. Overall, these findings demonstrated that miR-3658 as a tumor suppressor miRNA exerts its effect against OCT4 gene expression, and it has the potential of being used as a prognostic marker and therapeutic target against colorectal cancer.
Colorectal cancer (CRC) as a major health problem in industrialized countries, is the third most commonly diagnosed malignancy worldwide and the fourth most common cause of death [1,2]. Despite advances in screening procedures, up to 50% of all patients with CRC will develop metastases . Although the genetic and molecular basis of CRC metastasis is poorly understood, some recent studies suggested that OCT4 (OCT3/4 or POU5F1) as a transcription factor that controls the properties of immortality, un-differentiation, and invasion of cancer cells, may be involved in CRC metastasis [4–6].
OCT4 gene, a well-known pluripotent stem cell marker, encodes a protein member of POU family and has been implicated in solid tumors including non-small cell lung carcinomas, esophageal squamous cell carcinoma, and glioma, as well as gastric, breast [6–10], bladder cancers , and CRC tumors . It is also suggested that OCT4 is closely related to tumor progression , self-renewal , and drug resistance [15,16]. Therefore, ﬁnding a new modulator capable of regulating this transcription factor is necessary for new diagnostic and therapeutic strategies.
Oct4, Sox2, and Nanog are considered to be essential transcription factors for the maintenance of the pluripotent embryonic stem cell phenotype. The cooperative interaction of Sox2 and Oct4 is required for transcription of pluripotency gene expression [17,18].
MicroRNAs (miRNAs) are a family of mature noncoding small RNAs that act as negative gene regulators through destabilization or translational repression in the cytoplasm [19,20]. However, recent studies indicated the positive regulatory effect of miRNAs by targeting promoter elements . Notably, many studies reported that miRNAs exhibit temporally and cell-specific regulated expression patterns  and might regulate diverse physiological and developmental processes  such as proliferation, differentiation , and stemness . An abnormal miRNA expression signature could be an indicator of various types of cancer including CRC [26–28].
Here, we have introduced miR-3658 as a novel negative regulator of OCT4 gene expression and suggested a tumor suppressor role for it in CRC cells. Functional analysis of this miRNA revealed a negative posttranscriptional regulatory effect of miR-3658 on OCT4 gene expression which was followed by inhibition of proliferation and migration of CRC-originated SW480 cell line.
Material and methods
A bioinformatics phase was conducted to predict the potential miRNAs capable of targeting OCT4 mRNA. The 3′UTR sequence of human OCT4 mRNA was retrieved using Entrez (http://www.ncbi.nlm.nih.gov/Entrez/), then Several online tools such as target scan (http://www.targetscan.org/vert_72/) and miRmap (http://mirmap.ezlab.org/app/), were used for prediction of miRNA target genes. Conservation of miRNA recognition elements (MREs) was evaluated via Blat search from the UCSC genome browser (https://genome.ucsc.edu/). The expression data of CRC tissues and cell lines were retrieved from TCGA (https://tcga-data.nci.nih.gov/) NCBI and GEO (https://www.ncbi.nlm.nih. gov/geo/) databases.
A 490 bp of human genomic DNA, representing pre-miR-3658 sequence along with the flanking region was amplified by PCR and cloned into the pEGFP-C1 vector downstream of the GFP sequence. To construct reporter plasmid, the entire 3′UTR sequence of human OCT4 and Her2 gene was PCR-amplified, using a pair of oligonucleotides and cloned into pTG19-R/T vector. Then, it was subcloned downstream of the Renilla luciferase gene in the psi-CHECK-2 vector. For the mutant construct, Using OCT4 3′-UTR wild type as a template DNA, the F1 and R2 primers were used to amplify an upstream fragment with no putative miR-3658 binding sites. Primers and oligo sequences used in this study are listed in Table 1. All recombinant vectors were SANGER sequenced to verify the fidelity of the insert sequence.
|Primer name .||Primer sequences 5′ to 3′ end .|
|OCT4-3′UTR mutant||Forward; CTCGAGCATTCAAACTGAGGTGCCTG|
|Anchored oligo(dT)||GCGTCGACTAGTACAACTCAAGGTTCTTCCAGTCACGACG (T)18V|
|GAPDH||Forward; GCCACATCGCTCAGACAC |
|Primer name .||Primer sequences 5′ to 3′ end .|
|OCT4-3′UTR mutant||Forward; CTCGAGCATTCAAACTGAGGTGCCTG|
|Anchored oligo(dT)||GCGTCGACTAGTACAACTCAAGGTTCTTCCAGTCACGACG (T)18V|
|GAPDH||Forward; GCCACATCGCTCAGACAC |
Human tissue samples and cell lines
Total RNA was extracted from 15 pairs of fresh frozen CRC tumor tissue samples and their non-tumor adjacent tissues obtained from patients undergoing surgery in Emam hospital. The study protocol was approved by the research committee at Tarbiat Modares University and all patients provided informed consent for using of specimens in research. The CRC cell lines, SW480, and HCT116 cells were obtained from Pasteur Institute/Iran and cultured in DMEM-HG (Invitrogen) and RPMI (Invitrogen), respectively. HT29 and HEK293t cell lines were purchased from the National Center for Genetic and biological reserves in Iran and cultured in DMEM-HG (Invitrogen) and DMEM-F12 (Invitrogen), respectively. These media were supplemented with 100 U/ml penicillin, 100 µg/ml streptomycin (Sigma), and 10% fetal bovine serum (Invitrogen) at 37°C with 5% CO2.
For miRNA overexpression analysis, SW480, HT29, and HCT116 cells were seeded onto 12-well plates (1.2 × 105 cells per well) and 32 h later (∼70% confluent), were transfected with two micrograms of a pEGFP-C1 expression vector containing miR-3658 precursor using TurboFect transfection reagent (Thermo Fisher). Twenty-four hours after transfection, green fluorescent protein (GFP) expression was visually examined by fluorescent microscopy (Nikon Eclipses Te2000-s) and cells were harvested 48 h after transfection for the following experiments.
To study the direct interaction of miR-3658 and the 3′UTR sequence of the OCT4 gene, Hek293T cells were seeded onto 48-well plates (3 × 104 per well). Cells (∼70% confluent) were transfected with psi-CHECK-2 luciferase reporters (200 ng) and pre-miRNAs (400 ng) or mock (400 ng) using TurboFect transfection reagent (Thermo Fisher). Forty-eight hours after transfection using the dual-luciferase reporter system, the activities of the Firefly and Renilla luciferases were measured sequentially from cell lysates, according to the manufacturer's instructions (Promega). Firefly luciferase units were normalized against Renilla luciferase units to control the transfection efficiency. All experiments were done in triplicate.
RNA extraction and quantitative reverse transcription-PCR
Total RNA was isolated from the cells and tissues using Riboex reagent (GeneAll) according to the manufacturer's protocol. To synthesize cDNA, the RNA quality and quantity were analyzed by using agarose gel electrophoresis and spectrophotometry, respectively. Then, RNA was treated with DNase I (Takara Bio), followed by cDNA synthesis using prime script II reverse transcriptase (Takara Bio) as recommended by the manufacturer. For miRNA detection, before cDNA synthesis, poly A tail was added to the 3′ end of RNAs by using poly-A polymerase (Takara Bio) according to the manufacturer's protocol. Real-time quantitative PCR was performed using brilliant Sybr green qPCR reagent (BioFact) on an ABI instrument (Applied Biosystems, StepOne). Warrington, U.K.). The relative expression level for miRNAs was normalized to U48, and for mRNAs were normalized to B2M using the 2-ΔΔCt method. Primers and oligo sequences used in this study are listed in Table 1.
Western blot analysis
SDS–PAGE was carried out on 12% polyacrylamide gels and then Solubilized proteins (30 μg/μl) extracted from SW480 and HCT116 cells were separated by electrophoresis. Next, the gels were electroblotted onto PVDF membranes (Santa Cruz) at 100 V for 1.5 h at room temperature. The blocking of the membrane was carried out with 5% BSA in Phosphate-buffered saline containing 0.1% Tween for 1 h at room temperature. Overnight incubation of Oct4 primary antibody (1 : 500, Santa Cruz) was done at 4°C and later, was washed and incubated with a goat anti-mouse secondary antibody solution (1 : 1000, Santa Cruz) for 1 h at room temperature. Using enhanced chemiluminescence (ECL) detection system (Amersham, Piscataway, NJ) the protein bands were visualized and then quantitated by GelQuantNET Software. For verification of protein loading, the membranes were stripped and re-probed with β-actin antibody (Santa Cruz).
Cell cycle and MTT assay
Cell cycle analysis was performed on SW480 and HCT116 cells 48 h after transfection with either pre-miR3658 or mock control. Cells were trypsinized and washed twice with cold PBS, later fixed in 70% cold ethanol, and incubated with propidium iodide (PI). All of the samples were examined by a FACS Calibur flow cytometer, and the cell cycle populations were determined by using Cell Quest software (BDBiosciences).
For MTT assay analysis, SW480 and HCT116 cells were seeded onto 96-well plates (8000 cells per well) in tetra plicate. Forty-eight hours after transfection, 20 µl of the MTT solution was added to each well and then the cells were continued to culture for 4 h in 37°C with 5% CO2 incubator. The culture medium was then discarded and 100 μl DMSO (Sigma) was added to each well and the ELISA reader was used to measure the absorbance at 570 nm (Biotek) after crystals were dissolved completely.
Cell migration assay
To assess the migration capabilities of SW480 cells, a scratch wound assay was done. SW480 cells were transfected with pre-miR-3658 or mock control vectors and an artificial wound was created 24 after transfection using a 200-μl pipette tip on the confluent cell monolayer. Photomicrographs were obtained at 0, 24, 48, and 72 h using a 4X objective on a Nikon Eclipse inverted microscope and recorded using cellSens software. The data were analyzed using Image J software to estimate the rate of cell migration.
The Real-time experiment data were analyzed using the CT method by Data Assist software V3.0 and normalized by endogenous control U48 small nuclear RNA and B2M genes. Other statistical analysis was performed using GraphPad Prism 5.04 (GraphPad, San Diego, CA). All experiments were repeated at least twice with duplicate samples in each experiment and all of the results are presented as the means ± SEM.
Analysis of miR-3658 expression in CRC tumor tissues and cell lines and prediction of its putative target gene
About 6000 target genes were predicted for miR-3658 of between, the OCT4 gene was prominent. Both target scan and miRmap software predicted one putative binding site for miR-3658 within the 3′UTR sequence of OCT4 mRNA (Figure 1A). Algorithm scores implied on strong base pairing between miR-3658 and the OCT4 3′UTR MRE sequence (Figure 1B) and also, recognition site was conserved among Human, Chimp, Rhesus, Rabbit, Pig, and Cow using UCSC blat tool and target scan algorithm (Figure 1C).
. Analysis of miR-3658 expression level in CRC tissues and prediction of the
OCT4 gene as its target gene.
Analysis of microarray data adopted from GEO (GSE81581) demonstrated a significant differential expression of miR-3658 where it showed a higher expression level in non-tumor colorectal tissue samples, compared with colorectal-liver metastasis (Figure 1D). Data analysis of another microarray study of CRC cell lines with different metastasis states (GSE72412) showed a significantly higher expression level of miR-3658 in SW480 versus SW620 cells (Figure 1E). To verify the suppression of miR-3658 in CRC specimens, the mature miR-3658 expression level was measured using RT-qPCR. The results indicated that miR-3658 expression was strong (>70%) down-regulated in 15 CRC specimens in comparison with matched normal tissues (Figure 1F). Moreover, In the TCGA cohort, OCT4 gene expression although not significant but seemed to be increased in the CRC tissues compared with the normal ones (Figure 1H). Eventually, consistent to the RT-qPCR results against OCT4 gene expression in CRC tissue specimens (Figure 1G), bioinformatics analysis suggested that miR-3658 could have important functions in CRC through regulating the OCT4 gene expression.
Experimental expression analysis of miR-3658 and OCT4 gene in CRC originated cell lines
The expression status of miR-3658 and the OCT4 gene was also investigated in SW480 cells (grade III CRC) and HCT116 cells (grade II CRC), using RT-qPCR. While the miR-3658 expression level in SW480 cells was lower than its level in HCT-116 cells (Figure 2), the expression level of its putative target gene (OCT4) was an inverse (Figure 2). Hence, a reverse expression pattern between miR-3658 and OCT4 was deduced in these two CRCs originated cell lines. Later, these cells were used for subsequent experiments.
The reverse expression level of miR-3658 with the
OCT4 gene in CRC originated cell lines.
Cell-specific overexpression effect of miR-3658 against OCT4 gene expression
Transfection of the pEGFP-C1 vector containing pre-miR-3658 into SW480 and HCT116 cells resulted in significant overexpression of miR-3658 in both of these cell lines, detected by RT-qPCR (Figure 3A). Then, RT-qPCR and western blot analysis indicated that ectopic overexpression of miR-3658 in SW480 cells was followed by reducing endogenous OCT4 gene expression both at the mRNA (Figure 3B) and protein levels (Figure 3D). However, such an effect was not detected in HCT116 cells (Figure 3B,D). Also, overexpression of miR-3658 in SW480 cells resulted in a reduced SOX2 expression level, detected through RT-qPCR (Figure 3C).
Experimental evidences for cell-specific suppression of
OCT4 gene expression following miR-3658 overexpression.
To investigate a direct interaction between miR-3658 and OCT4 gene transcripts, luciferase reporter in pSICHECK2 vector was fused to indicate 3′UTR sequences (OCT4 3′UTR wild type, mutated OCT4 3′UTR, and HER2 3′UTR sequences). The HER2 3′UTR sequence was used as non-target, showing no target site for miR-3658. Then, luciferase reporters in pSICHECK2 vector along with pEGFP-C1 vector overexpressing miR-3658 were co-transfected in Heck293 cells. The pre-miR-3658 overexpression resulted in a significant (∼40%) suppression of the luciferase activity while it was fused to wild type 3′UTR sequence of the OCT4 gene. However, it did not affect the luciferase activity of HER2 reporter or mutated OCT4 3′UTR, compared with the mock control (Figure 3E). Taken together, these data supported the direct interaction of miR-3658 with the 3′UTR sequence of the OCT4 transcript.
Hsa-miR-3658 overexpression effect on cell cycle status in SW480 and HCT116 cells
To investigate the miR-3658 overexpression effect on HCT116 and SW480 cells, the relevant expression vector was transfected into these cells. Then, flow cytometry results indicated that the population of the transfected cells in the sub-G1 stage has been significantly increased in SW480 cells, while the G1 stage has been decreased, compared with the mock control transfection (Figure 4A). On the other hand, a repeat of this experiment in HCT116 cells resulted in no significant change of the cell population in sub-G1 or G1 stages, rather it resulted in the increased cell population at S and G2 phases (Figure 4B). The tumor suppressor effect of miR-3658 overexpression in SW480 cells was further supported via qRT-PCR against marker genes. While, overexpression of miR-3658 resulted in the reduced expression level of PCNA, Cyclin-D1, and c-Myc genes in SW480 cells, it resulted in no significant change in HCT116 cells (Figure 4C,D).
The effects of miR-3658 overexpression on cell cycle progression and proliferation.
To further validate the anticancer role of miR-3658 in CRC, its effect on different CRC cell lines (SW480, HCT116, and HT29) viability was investigated using MTT assay in different time points. The results demonstrated that miR-3658 overexpression significantly attenuated SW480 cells viability at 24, 48, and 72 h after transfection. While, it did not affect HCT116 and HT29 cell viability, detected through MTT assay (Figure 4E).
Suppression of SW480 cell migration following miR-3658 overexpression
Considering the role of OCT4 in promoting cell migration, we investigated the effect of miR-3658 overexpression on the migration capacity of SW480 cells, using a wound-healing assay. Results indicated that the migration rate of SW480 cells has been significantly reduced following the overexpression of miR-3658, compared with mock control (Figure 5A). Once again, these data further denoted the biological importance of miR-3658 in CRC progression.
Effect of miR-3658 overexpression on SW480 cell migration.
Secondly, to assess the miR-3658 effect on epithelial-to-mesenchymal transition, Claudin1, B-catenin (epithelial marker [29,30]), ZEB1, and Snail (mesenchymal markers [30,31]) genes expression were checked in SW480 cells. An increase in the expression of the epithelial marker and a decrease in expression of the mesenchymal markers indicated that overexpression of miR-3658 reduced SW480 invasive capacity (Figure 5B).
CRC is one of the most prevalent malignant tumors in the world in terms of high cellular/molecular heterogeneity . Increasing evidence has demonstrated the importance of cancer stem cells in cancer recurrence, drug resistance and metastatic epithelial-to-mesenchymal transition [33–35]. Recent studies support the contribution of Oct4 protein as a key transcription factor in self-renewal and pluripotency of embryonic stem cells, in multiple pathological processes that contribute to tumorigenesis [36–38]. That is why understanding the molecular mechanism of OCT4 deregulation and finding the regulators of this gene could be beneficial for future cancer research. Recent studies have shown the presence of cancer stem cells in a variety of leukemia and tumor tissues. These cell types play an important role in the process of tumor cell invasion and metastasis. Tumor microenvironment, miRNAs and their controlling factors are also important regulators of cancer stem cell function [39,40].
Given that, miRNAs are involved in the regulation of multiple pathological processes through targeting numerous genes and play important role in cancer progression [40,41], we have searched for the candidate miRNAs which are potentially involved in CRC incidence through regulating OCT4 gene expression.
RT-qPCR expression analysis of miR-3658 in 15 pairs of CRC specimens consistent with RNA seq data indicated that this miRNA is down-regulated in CRC. Consistently, the expression status of OCT4 gene was investigated both in CRC specimens and RNAseq data and results indicated the up-regulation of this gene in tumors.
The oncogenic function of OCT4 has already been reported in several cancers such as breast, bladder, cervical, etc. [42,43]. Moreover, Zhou et al.,  have introduced OCT4 as a CRC marker so that the expression of it, is significantly associated with the development and prognosis of CRC. Besides, it is reported that the up-regulation of the OCT4 transcript in somatic tissues is associated with cancer stem cells .
We have used HT29, HCT116, and SW480 CRC originated cell lines in our experiments, representing different grads of CRC progression. However, the expression level difference of both miR-3658 and OCT4 expression in these cell lines could be attributed to the different environmental status of these cell lines.
The interactions between miRNAs and environmental factors are play critical roles in determining abnormal phenotypes and diseases . So, the cell lines selection is a critical first step in miRNA research. There are different CRC cell lines [3–5] which are representing different grads and stages of this cancer.
Caco2 and HT29 are the models of less metastatic and more differentiated colon cancer cells. However, HCT116 and SW480 cell lines contain CSCs . Specially, SW480 cells are frequently used as a model of more metastatic undifferentiated cancer cells. Besides, the SW480 had higher expression of OCT4 at both the gene and protein levels signifying SW480 has properties of CSCs [47,48]. Here, due to the presence of CSCs in the SW480 and HCT116 cell lines, they were selected for investigating the expression status of miR-3658 and OCT4 gene at RNA and protein levels.
Consistent with other studies , the results of the expression analysis, indicated that the OCT4 expression level in SW480 cells was relatively higher than its level in HCT116 cells. However, the miR-3658 expression level in both cell lines was reverse to the expression level of the OCT4 gene (Figure 2). Consistent with the notion that miR-3658 is a negative regulator of OCT4 gene expression, miR-3658 overexpression in SW480 cells resulted in reduced OCT4 expression levels both at the mRNA and protein levels. Interestingly, such inhibitory effect for miR-3658 against OCT4 gene expression was not detected in HCT116 cells, probably due to the much lower level of OCT4 gene expression in this cell line. Based on previous studies, miRNA target genes with low mRNA expression levels are probably under strong regulation by endogenous miRNAs. Likewise, miRNA target genes with high mRNA expression are less likely to have strong internal miRNA regulation . Therefore, since the expression of OCT4 in HCT116 cells is relatively lower and the expression of miR-3658 is higher than in SW480 cells (Figure 2), consistent to the observed results OCT4 gene expression was less affected in HCT116 cells, following the overexpression of exogenous miR-3658 (Figure 3).
SOX2 transcription factor and stemness marker is known to be regulated by Oct4 protein. In SW480 cells with up-regulated miR-3658 expression level, SOX2 gene expression was reduced, compared with cells transfected with the mock control construct. Down-regulation of SOX2 following the overexpression of miR-3658 in SW480 cells, is consistent to its effect on OCT4 gene expression.
Since the only ∼20% of changes in expression patterns of target genes could be attributed to the direct miRNA–mRNA interactions, investigation of direct interaction of miR-3658 with OCT4 3′-UTR seemed necessary. The results of the dual-luciferase assay supported the direct interaction of miR-3658 with the OCT4 gene 3′UTR transcript. Once again, these results supported the regulatory effect of miR-3658 against OCT4 gene expression through direct interaction.
Multiple studies supporting the role of OCT4 in cell cycle and apoptosis regulation by modulating cell cycle progression in G0/G1. Indeed, both the N-terminal and C-terminal transactivation domains of Oct-4 are important for the stimulation of S-phase entry . For instance, OCT4, regulates CCND1 expression, as a bottleneck gene, and subsequently promoting cell cycle progression . P21, a member of the universal cyclin-dependent kinase inhibitors (CDKIs), is another gene which its expression can be regulated by OCT4 and plays critical roles in the regulation of the G1/S transition . To further explore the putative tumor-suppressive effect of miR-3658 in human CRC cell lines, the mechanism of miR-3658 effect on the cell growth and viability in SW480 and HCT116 cell lines was examined. Results indicated that unlike in HCT116 cells in which OCT4 expression level is low, overexpression of miR-3658 in SW480 cells significantly inhibited cell proliferation as evidenced by cell cycle, cell cycle-associated genes expression, and MTT assays (Figure 4). The negative effect of miR-3658 on SW480 cell cycle and survival was confirmed through RT-qPCR against c-Myc (transcription factors associated with cell cycle control) and PCNA (a marker of replicating cells) genes, as well. Unlike in HCT116 cells, results supported that miR-3658 gain of function led to suppressing cell cycle progression in SW480 cells (Figure 4C,D). Together with the observed miR-3658 effects on cell cycle proliferation, these results suggest that miR-3658 regulates SW480 cell growth, at least in part, by targeting OCT4.
Metastasis is the neoplastic process that is responsible for the causes of cancer death. Recent studies reported that OCT4 may be involved in metastasis of CRC through the epithelial-mesenchymal transition (EMT) process [54,55]. EMT plays a critical role during the progression of cancers and it is necessary for metastasis of epithelial cancer . CRC cells acquire the invasive phenotypes through alterations in cell morphology, cell–cell adhesion, the activity of cellular signaling pathways, and the extracellular matrix . OCT4 gene has been described as one of the important factors through the regulation of epithelial-to-mesenchymal transition [58–60]. To further investigate whether the miR-3658 is associated with the progression of CRC, the wound-healing assay was performed in SW480 cells. Consistent with its suggested tumor suppressor effect, the up-regulation of miR-3658 was followed by inhibition of SW480 cell migration in a scratch-based assay (Figure 5). This migration suppressing effect of miR-3658, is probably related to the role of this miRNA in regulating OCT4 gene expression. Likewise, EMT-involved genes expression in SW480 cells was decreased upon miR-3658 overexpression, which suggests that miR-3658 may also function as a metastasis suppressor probably, by regulation OCT4 gene expression.
In summary, our ﬁndings suggested that down-regulated miR-3658 expression in CRC samples is related to CRC progression. The results of the current study show that, miR-3658 is a regulator of cell proliferation and migration in SW480 cells by targeting the OCT4 gene, suggesting that miR-3658 may serve as a promising target and diagnostic marker for treatment of CRC patients. Although these findings are promising, further data collection would be needed to determine exactly how miR-3658 is down-regulated and affects CRC.
The authors declare that there are no competing interests associated with the manuscript.
F.H. and B.M.S. conceived and designed the experiments. F.H. performed the experiments. F.H. and B.M.S. analyzed the data. B.M.S., S.H., and H.B., administrative and financial support. F.H and B.M.S. wrote the manuscript.
This study was funded by the Tarbiat Modares University and the National Institute for Medical Research Development (NIMAD) to B.M.S.
We express our appreciation to the members of the 4402 laboratory in TMU and Hassan Ansari at Royan Institute for their technical assistance.