|Home | About | Journals | Submit | Contact Us | Français|
The deregulation of microRNA (miRNA) is frequently associated with a variety of cancers, including hepatocellular carcinoma (HCC). In this study, we identified 10 upregulated miRNAs (miR-217, miR-518b, miR-517c, miR-520g, miR-519a, miR-522, miR-518e, miR-525-3p, miR-512-3p and miR-518a-3p) and 10 downregulated miRNAs (miR-138, miR-214, miR-214#, miR-27a#, miR-199a-5p, miR-433, miR-511, miR-592, miR-483-5p and miR-483-3p) by Taqman miRNAs array and quantitative real-time PCR (qRT–PCR) confirmation. Additionally, we investigated the expression and possible role of miR-138 in HCC. qRT–PCR results showed that miR-138 was downregulated in 77.8%(14/18) of HCC tissues compared with adjacent non-tumor tissues. Overexpression of miR-138 reduced cell viability and colony formation by induction of cell arrest in HCC cell lines and inhibited tumor cell growth in xenograft nude mice. The use of miR-138 inhibitor increased cell viability and colony formation in HCC cell lines and tumor cell growth in xenograft nude mice. Using TargetScan predictions, CCND3 was defined as a potential direct target of miR-138. Furthermore, CCND3 protein expression was observed to be negatively correlated with miR-138 expression in HCC tissues. The dual-luciferase reporter gene assay results showed that CCND3 was a direct target of miR-138. The use of miR-138 mimic or inhibitor could decrease or increase CCND3 protein levels in HCC cell lines. We conclude that the frequently downregulated miR-138 can regulate CCND3 and function as a tumor suppressor in HCC. Therefore, miR-138 may serve as a useful therapeutic agent for miRNA-based HCC therapy.
Hepatocellular carcinoma (HCC) is one of the most common cancers worldwide (1). The development and progression of HCC is typical of a multistage process, which is believed to involve the deregulation of genes that are critical to cellular processes, such as cell cycle control, cell growth, apoptosis and cell migration and spreading. In the past decades, studies have focused on investigating the genes and proteins underlying the development and progression of HCC (2). Recently, an increasing number of reports have described a new class of small regulatory RNA molecules termed microRNAs (miRNAs) that are implicated in HCC progression (3).
MiRNAs are endogenous non-coding 20–22 nucleotide RNAs that have been identified as post-transcriptional regulators of gene expression (4). The miRNAs mainly bind to the 3′ untranslated regions (UTRs) of target messenger RNAs (mRNAs), resulting in mRNA degradation or the blockade of mRNA translation. Increasing evidence showed that miRNAs have significant roles in diverse biological processes (5). Meanwhile, deregulation of miRNAs has been observed in a wide range of human diseases, including cancer (6). In human cancer, miRNAs can function as oncogenes or tumor suppressor genes during tumor development and progression (7).
In this study, the expression profiles of 667 miRNAs were examined in human HCC and adjacent non-tumor tissues. A set of significantly differentially expressed miRNAs were identified in HCC tissues. Further investigation revealed that miR-138, a frequently downregulated miRNA in HCC, could induce cell cycle arrest in HCC cell lines. Moreover, cyclin D3 (CCND3) was characterized as a direct target of miR-138.
Eighteen pairs of human HCC and adjacent non-tumor tissues were obtained from surgical specimens immediately after resection from patients undergoing primary surgical treatment of HCC in the Eastern Hepatobiliary Surgery Hospital, Shanghai, China. No patient had received preoperative irradiation or chemotherapy. The samples were frozen in liquid nitrogen and stored at −80°C until use. Among these samples, three pairs were used for Taqman low density miRNA array (TLDA) analysis and all of them were used for quantitative real-time PCR (qRT–PCR) analysis. Clinical and pathological information was extracted from the patients’ medical charts and pathology reports (Table I). Written consent for tissue donation (for research purposes) was obtained from the patients before tissue collection and the protocol was approved by the Institutional Review Board of Eastern Hepatobiliary Surgery Hospital and Second Military Medical University.
Human embryonic kidney cells (HEK293T) and human HCC cell lines HepG2 and Huh7 were cultured in Dulbecco's modified Eagle's medium plus 10% fetal bovine serum (Invitrogen, Carlsbad, CA) at 37°C in a humidified atmosphere containing 5% CO2.
Total RNA was isolated using mirVana miRNA isolation kit (Ambion, Austin, TX). For miRNA complementary DNA (cDNA) synthesis, RNA was reverse transcribed using the miRNA reverse transcription kit (Applied Biosystems, Foster City, CA) in combination with the stem-loop Megaplex primer pool (Applied Biosystems). Then, TLDA v2.0 was performed as the manufacturer’s protocol. For each cDNA sample, 667 small RNAs were profiled. All PCR reactions were performed on the 7900HT real-time PCR system (Applied Biosystems). Cycling conditions were as follows: 95°C for 10 min followed by 40 cycles of 95°C for 15 s and 60°C for 1 min. Human U6 small RNA was used as an internal control to normalize RNA input. The data were analyzed by SDS v2.3 software. In details, the Ct value is defined as the fractional cycle number at which the fluorescence passes the fixed threshold. The fold change was calculated using the 2−ΔΔCt method and presented as the fold-expression change in tumors relative to their corresponding normal tissues after normalization to the endogenous control.
Synthesis of cDNA and qRT–PCR analysis of miRNA expression was carried out with TaqMan microRNA assay kits (Applied Biosystems) according to the manufacturer’s protocol. Briefly, total RNA was extracted using TRIzol Reagent (Invitrogen) from clinical samples or HCC cell lines and were used to synthesize cDNAs with gene-specific primers. Reverse transcriptase (RT) reactions contained 100 ng RNA, 50 nmol/l stem-loop RT primers, 1× RT buffer, 0.25 mmol/l each of the deoxyribonucleotide triphosphates, 3.33 U/μl MultiScribe reverse transcriptase and 0.25 U/μl ribonuclease inhibitor. The 15 μl reactions were incubated for 30 min at 16°C, 30 min at 42°C, 5 min at 85°C and then held at 4°C. The cDNA product was used for the following qRT–PCR analysis directly. The 20 μl PCR reaction included 1.33 μl RT product, 1× TaqMan universal PCR master mix and 1 μl primers and probe mix of the TaqMan microRNA assay kit. Reactions were incubated in a 96-well optical plate at 95°C for 5 min, followed by 40 cycles at 95°C for 15 s and 60°C for 1 min. PCR reactions were run on a StepOne Plus real-time PCR machine (Applied Biosystems) and the data were analyzed by SDS v2.3 software same as in TLDA.
The fragment of 3′UTR of CCND3 (1054–2061 nt, Genbank accession no. NM_001136017.2) containing the two putative miR-138-binding sequences (1279–1285 nt and 1346–1352 nt) was amplified with the primers 5′-CCCTGGAGAGGCCCTCTGGA-3′ (forward) and 5′-TTCCAAGAAGCCAAAGCCAG-3′ (reverse). The PCR product was cloned into Firefly luciferase reporter vector pGL3 (Promega Corporation, Madison, WI), termed as pGL3-CCND3-3′UTR. The plasmid which carried the mutated sequence in the two complementary sites for the seed region of miR-138 (Figure 2A) was generated based on pGL3-CCND3-3′UTR plasmid by MutanBEST Kit (Takara Bio, Shiga, JP), termed as pGL3-CCND3-3′UTR-mut. Briefly, a forward miR-138 mutagenic primer F1-GAT ACATAGTGGTCGATTCCTTTTGAACGCCCCCCCCCACCCC and a reverse miR-138 mutagenic primer R1-CAATTCTGTCCCATCAGCCTGGCCCSCCCCCSGCTSGSGTTG or a forward miR-138 mutagenic primer F2-CTGCATCTGTGGTCGAGCGCCTTTCCCAACT CTAGCTGGGG and a reverse miR-138 mutagenic primer R2-GGAGGAGGAGCTTGACTAGCCACCGAAATGCAGACATGG, were synthesized and utilized in PCR experiment as described by the manufacturer. Two point mutations along the miR-138 seed binding sites were confirmed by sequencing analyses.
The transfections were carried out using FuGene HD transfection reagent (Roche, Indianapolis, IN) following the manufacturer’s protocol. In brief, 2 × 104 HepG2 and Huh7 cells or 5 × 104 HEK293T cells in 24-well plate were transfected with indicated plasmid DNA, miRNA duplex (GenePharma, Shanghai, China), small interfering RNA (siRNA) (GenePharma) and collected 24–48 h after transfection for assay.
Twenty-four hours after transfection, 1000 transfected HepG2 and Huh7 cells were placed in a fresh 96-well plate in triplicate and maintained in Dulbecco's modified Eagle's medium containing 10% fetal bovine serum for 5 days. Cells were tested for proliferation per 24 h using Cell Titer-Blue cell viability assay (Promega Corporation) according to the manufacturer's instructions and the fluorescence value was recorded by multi-plate reader (Synergy 2; BioTek, Winooski, VT).
Forty-eight hours after transfection, 1 × 105 transfected HepG2 and Huh7cells were harvested, washed once in phosphate-buffered saline and fixed in 70% ethanol at 4°C overnight. Staining for DNA content was performed with 50 mg/ml propidium iodide and 1 mg/ml ribonuclease A at room temperature for 30 min. Populations in G0–G1, S and G2–M phase were measured by Cell Lab Quanta SC flow cytometry (Beckman Coulter, Fullerton, CA) and the data were analyzed by FlowJo v7.6 Software.
Twenty-four hours after transfection, 2000 transfected HepG2 and Huh7cells were placed in a fresh six-well plate in triplicate and maintained in Dulbecco's modified Eagle's medium containing 10% fetal bovine serum for 2 weeks. Cell colonies were fixed with 20% methanol and stained with 0.1% coomassie brilliant blue R250 at room temperature for 15 min. The colonies were counted by ELIspot Bioreader 5000 (BIO-SYS, Karben, GE).
Male BALB/c nude mice (5–6 weeks of age) were obtained from Shanghai Experimental Animal Center (Shanghai, China). Animal handling and experimental procedures were approved by the Animal Experiments Ethics Committee of Second Military Medical University. For in vivo tumorigenicity assay, all pyrimidine nucleotides in the miR-138 mimic, miR-138 inhibitor or nonrelative control (NC) duplex were substituted by their 2′-O-methyl analogues to improve RNA stability. MiR-138 mimic-or miR-138 inhibitor-transfected HepG2 cells (1 × 105) were suspended in 100 μl phosphate-buffered saline and then injected subcutaneously into left side of the posterior flank of six BALB/c nude mice, respectively. NC transfected or non-transfected HepG2 cells (1 × 105) were injected subcutaneously into right side of same 12 mice. Tumor growth was examined daily and the tumor volumes were calculated every week using the formula for hemi-ellipsoids: V =length (cm) × width (cm) × height (cm) × 0.5236. After 5 weeks, the mice were sacrificed and the tumors were dissected and photographed.
HEK293T cells seeded in 24-well plate in triplicate were cotransfected with pGL3-CCND3-3′UTR or pGL3-CCND3-3′UTR-mut and miRNA-138 mimic or non-relative control RNA duplex (NC duplex; GenePharma) by using FuGene HD transfection reagent. The pRL-TK (Promega Corporation) was also transfected as a normalization control. Cells were collected 48 h after transfection, and luciferase activity was measured using a dual-luciferase reporter assay kit (Promega Corporation) and recorded by multi-plate reader (Synergy 2; BioTek).
Protein extracts from HCC tissuses or HepG2 cells were prepared by a modified RIPA buffer with 0.5% sodium dodecyl sulfate in the presence of proteinase inhibitor cocktail (Complete mini; Roche). Twenty-five micrograms of protein of HCC tissues and their adjacent non-tumorous tissues were electrophoresed in 10% sodium dodecyl sulfate–polyacrylamide gel electrophoresis mini gels and transferred onto polyvinylidene difluoride membranes (Immobilon P−SQ; Millipore, Billerica, MA). After blocking with 5% non-fat milk, the membranes were incubated with rabbit anti-cyclin D3 antibody (1:1000 dilution; Epitomics, Burlingame, CA) or mouse anti-glyceraldehyde-3-phosphate dehydrogenase antibody (1:5000 dilution; Epitomics) at 4°C overnight, followed by incubation with horseradish peroxidase-conjugated goat anti-rabbit or goat anti-mouse antibody (1:10000 dilution; KPL, Gaithersburg, MA) for 1 h at room temperature. Finally, signals were developed with Super Signal West Pico chemoluminescent substrate (Pierce, Rockford, IL), visualized by the GeneGnome HR Image Capture System (Syngene, Frederick, MD) and analyzed by Gene tools (Syngene).
To investigate the target genes of miRNAs and the conserved sites bound by the seed region of miR-138, the TargetScan (http://www.targetscan.org/), MiRanda (http://www.microrna.org/microrna/home.do) and PicTar (http://pictar.mdc-berlin.de/) programs were used.
The top 25% miRNA targets that had been assigned the highest numbers of miRNA interaction sites were collected and subjected to Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis. GO analysis was applied to analyze the main function of the differential expression genes according to the GO, which is the key functional classification of National Center for Biotechnology Information (8).Generally, Fisher’s exact test and χ2 test were used to classify the GO category, and the false discovery rate (FDR) (9) was calculated to correct the P-value, the smaller the FDR, the smaller the error in judging the P-value. The FDR was defined as , where Nk refers to the number of Fisher’s test P-values less than χ2 test P-values. P-values were computed for the GOs of all the differential genes. Enrichment provides a measure of the significance of the function: as the enrichment increases, the corresponding function is more specific, which helps us to find those GOs with more concrete function description in the experiment. Within the significant category, the enrichment Re was given by ,where nf is the number of differential genes within the particular category, n is the total number of genes within the same category, Nf is the number of differential genes in the entire microarray and N is the total number of genes in the microarray (10). Similarly, pathway analysis was used to find out the significant pathway of the differential genes according to KEGG, Biocarta and Reatome. Still, the Fisher’s exact test and χ2 test were used to select the significant pathway, and the threshold of significance was defined by P-value and FDR. The enrichment Re was calculated like the equation above (11–13).
Data are presented as mean ± SD. Comparisons were made by using a two-tailed t-test or one-way analysis of variance for experiments with more than two subgroups. Correlation analysis was made by using Spearman correlation coefficient. P <0.01 was considered statistically significant.
HCC and adjacent non-tumor tissues were obtained from 14 male and 4 female patients and their profiles are shown in Table I. The average age of these patients was 48.5 years (ranged from 35 to 71 years). Eleven patients (61.1%) had tumors >5 cm, two (11.1%) had tumors <2 cm and five (27.8%) had 2–5 cm tumors. Seven (38.9%) cases were well differentiated (I + II) and 11(61.1%) cases were moderately differentiated (III). Four patients (22.2%) were diagnosed at invasion T1, seven (38.9%) at invasion T2, six (33.3%) at invasion T3 and one (5.6%) at invasion T4.
The TLDA were performed in three pairs of HCC and adjacent non-tumor tissues. Only miRNAs altered by at least 3-fold in all three pairs of the samples were considered significant candidates. Under these strict criteria, 11 upregulated miRNAs and 12 downregulated miRNAs were identified. To validate the miRNA array data, qRT–PCR was performed in 18 pairs of HCC tissues. Ten upregulated and 10 downregulated miRNAs showed consistent changes in >50% tumorous tissues (Table II), of which miR-217, miR-520g, miR-522 and miR-525-3p were upregulated in >70% tumorous tissues and miR-199a-5p, miR-138, miR-483-5p and miR-511 were downregulated in >70% tumorous tissues (Table II).
The targets of the above 20 deregulated miRNAs were predicted by TargetScan and the representative target genes were listed in Table II. To judge the most significant candidates and investigate the cellular function, the signaling pathway and GOs of these target genes were analyzed. The results showed that a wide variety of cellular processes were featured significantly in signaling pathways (Supplementary Table S1 is available at Carcinogenesis Online). Many of these signaling pathways, such as insulin, mitogen-activated protein kinase, transforming growth factor-beta and Wnt signaling pathway, have been shown to participate in the tumorigenesis (14–17). However, some other signaling pathways have never been reported to play a role in tumorigenesis, e.g. axon guidance. Among all these differentially regulated signaling pathways, ‘regulation of actin cytoskeleton’ and ‘pathway in cancer’ appeared to be the most enriched two in both upregulated and downregulated miRNA groups. A similar phenomenon was observed in GOs analysis. Many cellular functions were featured significantly, of which the ‘signal transduction’ appeared to be the most enriched one (Supplementary Table S2 is available at Carcinogenesis Online). Based on these results, the miRNAs which were involved in both signal transduction and the regulation of actin cytoskeleton or pathway in cancer might be the most significant candidates. MiR-520g, miR-483-5p, miR-138, miR-199a-5p, miR-217 and miR-518e were selected under the criteria for further studies. Here, the studies on miRNA-138 were presented.
The significant reduction of miR-138 expression in HCC tissues indicated possible biological significance in tumorgenesis. At first, the effect of miR-138 on cell growth was evaluated in HepG2 and Huh7 cells transfected with or not, miR-138 mimic, miR-138 inhibitor or NC duplex. The expression of miR-138 was increased 25-folds (HepG2) and 22-folds (Huh7) in cells transfected with 20 nM miR-138 mimic but deceased 30-fold (HepG2) and 43-folds (Huh7) in cells transfected with 20 nM miR-138 inhibitor (Figure 1A). From second day (HepG2) or third day (Huh7) after the transfection, the viability of cells transfected with miR-138 mimic significantly decreased compared with that of NC duplex transfected or non-transfected cells, but the viability of cells transfected with miR-138 inhibitor significantly increased (Figure 1B). These results indicate miR-138 could inhibit cell growth.
To validate the inhibitory effect of miR-138 on cell growth, the colony formation assay was performed in HepG2 and Huh7 transfected with or not, miR-138 mimic, miR-138 inhibitor or NC duplex. As showed in Figure 1C, HepG2 and Huh7 cells transfected with 20 nM miR-138 mimic displayed much fewer and smaller colonies (218 or 161 colonies) compared with NC duplex transfected (783 or 729 colonies) and non-transfected cells (756 or 692 colonies), but cells transfected with 20 nM miR-138 inhibitor displayed much more and larger colonies (1238 or 1349 colonies).
To further confirm the above findings, an in vivo mouse model was used. For the duration of the treatment with miR-138 mimic or miR-138 inhibitor for 5 weeks, tumor volume curves revealed a significant decrease in growth rates at the third, fourth and fifth week after treatment with miR-138 mimic and a significant increase in growth rates at the fourth and fifth week after treatment with miR-138 inhibitor whereas no significant differences in tumor growth rates were observed between the NC group and the control group (Figure 1D). These results indicate that introduction of miR-138 significantly inhibits tumorigenicity of HepG2 cells in xenograft nude mouse model.
To investigate the mechanism of inhibitory effect of miR-138, flow cytometry assay showed that the percentages of miR-138 mimic transfected HepG2 and Huh7 cells in the G0–G1 phase were 18% (HepG2) and 11% (Huh7) higher than that of NC duplex transfected or non-transfected cells, which paralleled with a 44% (HepG2) and 38% (Huh7) decrease in the S phase (Figure 1E). In miR-138 inhibitor-transfected cells, the percentages of cells in the G0–G1 phase were 8% (HepG2) and 6% (Huh7) less than that of NC duplex transfected or non-transfected cells, which paralleled with a 14% (HepG2) and 19% (Huh7) increase in the S phase (Figure 1E). These results indicate miR-138 could inhibit HepG2 and Huh7 proliferation by inducing cell cycle arrest at G1/S phase.
It is generally accepted that miRNAs exert their function through regulating the expression of their downstream target genes. CCND3 was predicted as a potential target of miR-138 by TargetScan and miRanda but not by PicTar. The 3′UTR of CCND3 mRNA contained a complementary site for the seed region of miR-138 (Figure 2A). CCND3 was found to be upregulated in 72.2% (13 of 18) HCC tissues compared with adjacent non-tumor liver tissues (Figure 2B). Furthermore, a negative correlation was found between the upregulated CCND3 protein and downregulated miRNA-138 (r = −0.730, P < 0.01) (Figure 2C). These results indicated that miR-138 may be associated with CCND3 and both of them may be involved in HCC tumorigenesis.
To validate whether CCND3 is a direct target of miR-138, a human CCND3 3′UTR fragment containing wild-type or mutant miR-138 binding sequence (Figure 2A) was cloned downstream of the firefly luciferase reporter gene in pGL3. In HEK293 cells cotransfected with the reporter plasmids and miR-138 mimic or NC duplex, the luciferase activity of the reporter that contained wild-type 3′UTR was significantly suppressed by miR-138 mimic, but the luciferase activity of mutant reporter was unaffected (Figure 2D), indicating that miR-138 may suppress gene expression through miR-138 binding sequence at the 3′UTR of CCND3. Furthermore, transfection of miR-138 mimic decreased CCND3 expression and transfection of miR-138 inhibitor increased CCND3 expression in HepG2 cells at protein (Figure 2E) but not mRNA level (data not shown), suggesting that CCND3 expression could be inhibited by miR-138 at post-transcriptional level. Together, the results show that miR-138 could regulate the expression of endogenous human CCND3 by directly targeting the 3′UTR of CCND3 mRNA and human CCND3 is a new target of miR-138.
To identify whether inhibition of CCND3, just like miR-138 restoration, also resulted in HCC repression, the effects of knockdown of CCND3 on cell growth were examined. First, HepG2 cells were transfected with or not, CCND3 siRNA or control siRNA. Seventy-two hours after transfection, a dose-dependent knockdown of CCND3 was observed in HepG2 cells (Figure 3A). In cell viability assay and cell cycle analysis, in vitro knockdown of CCND3 repressed cell viability (Figure 3B), induced cell cycle arrest (Figure 3C) and inhibited the colony formation (Figure 3D). The similar data were obtained in Huh7 cells transfected with CCND3 siRNA (data not shown). These results indicate that CCND3 is most likely involved in the induction of cell cycle arrest by miR-138.
MiRNAs were found to be frequently deregulated in HCC, and some specific miRNAs were found to be associated with the clinicopathological features of HCC, such as metastasis, recurrence and prognosis (18–20). Moreover, compelling evidence demonstrates that miRNAs have important roles in HCC progression and directly contribute to the cell proliferation, avoidance of apoptosis and metastasis of HCC. Identifying the miRNAs and their targets that are essential for HCC progression may provide promising therapeutic opportunities.
In this study, with TLDA and real-time reverse transcription–PCR confirmation, 10 upregulated miRNAs (miR-217, miR-518b, miR-517c, miR-520g, miR-519a, miR-522, miR-518e, miR-525-3p, miR-512-3p and miR-518a-3p) and 10 downregulated miRNAs (miR-138, miR-214, miR-214#, miR-27a#, miR-199a-5p, miR-433, miR-511, miR-592, miR-483-5p and miR-483-3p) were identified in HCC. Important of these 20 deregulated miRNAs, only miR-199a-5p was reported to contribute to the increase of cell invasion by functional deregulation of DDR1 activity in HCC (21) and regulate Brm subunit of SWI/SNF in human cancers (22), and the other 19 deregulated miRNAs were first reported to be involved in HCC tumorigenesis. MiR-27a#, miR-214#, miR-518a-3p and miR-518e have never been reported in literature. The rest 15 miRNAs were reported in various cancer but not HCC. For example, the upregulated miR-512-3p and miR-525-3p were associated with a cisplatin-resistant phenotype in human germ cell tumors (23). The upregulated miR-519a and downregulated miR-511 and miR-485-5p were associated with histological subtypes in ovarian cancers (24). MiR-517c and 520g promotes in vitro and in vivo oncogenicity, modulates cell survival and robustly enhances growth of untransformed human neural stem cells in neuroectodermal brain tumors (25). MiR-433 could regulate tumor-associated proteins GRB2 in gastric carcinoma (26). MiR-592 was reported to be associated with the stepwise progression for transformation from normal colon to carcinoma (27). Thus, the role of these deregulated miRNAs in HCC tumorigenesis is pretty worthy to be investigated.
Compared with the previous version, the more newly found miRNAs (337 updated miRNAs) were profiled in TLDA v2.0, which may be the reason so many new differentially expressed miRNAs were identified in this study. Meanwhile, very strict criteria of selection defined in TLDA (significantly different in all three pairs of HCC tissues) led to the absence of some well-known HCC-related miRNAs. Such as miR-21 (28), miR-122 (29,30), miR-16 (31) and miR-29 (32,33) were excluded for their significantly differential expression only in two pairs of HCC tissues, and miR-181 (34), miR-221 (35,36), miR-125 (37) and miR-101 (38) were excluded for their differential expression in only one pair of HCC tissues or no differential expression at all. Using such strict criteria, we might miss some important candidates but catch some unique ones. Fortunately, qRT–PCR and follow-up studies proved that our strategy helped find the above new deregulated miRNAs.
All of HCC patients included in this study were HBsAg positive. Therefore, the expression patterns of identified miRNAs may mainly represent the alterations in Hepatitis B virus-positive HCC, which may partially account for the inconsistency between the others and our results. Although the HCC tissues from four tumor grades (G1–G4) were included, we were unable to perform statistical comparison on the profiles of differentially expressed miRNAs among four tumor grades due to the quantity restriction. Future studies including larger size of four tumor grades samples are required to elucidate the miRNAs associated with the clinicopathological features of HCC, which might severed as the biomarker for HCC.
KEGG pathway and GO enrichment analysis, which was based on the reported and predicted target genes of these deregulated miRNAs, were applied to identify which particular functions and pathways were enriched among genes controlling distinctive characters between HCC and adjacent non-tumor tissues. As a result, KEGG pathway analysis showed that proliferative (cell cycle, mitogen-activated protein kinase and Wnt), adhesive (actin cytoskeleton, adherens junction and focal adhension), survival (transforming growth factor-beta and ErbB) and oncogenic (renal cell carcinoma and pancreatic cancer) signaling pathways were abundant among the significantly enriched ones. Furthermore, the GOs related to signal transduction (signal transduction, small guanosine triphosphatase-mediated signal transduction, protein–amino acid phosphorylation and nerve growth factor receptor signaling pathway) and cell growth (cell differentiation, cell division, positive regulation of cell proliferation, cell cycle and multicellular organism development) represented up to 37% of the significantly enriched GOs. As expected, various cell process and signal pathways were involved in HCC tumorigenesis. To narrow down the scope of study and judge the most significant candidates, the intersection of signal transduction and regulation of actin cytoskeleton and pathway in cancer was selected for follow-up study.
Based on TargetScan search, miR-138 has 388 conserved targets, 32% of which belong to the intersection of signal transduction and regulation of actin cytoskeleton and pathway in cancer. In addition, we showed that miR-138, frequently downregulated in human HCC tissues, could suppress colony formation and induce cell cycle arrest in hepatoma cell lines and inhibit tumor growth in xenograft nude mice model. These results suggest the important role of miR-138 in HCC tumorigenesis. Reduced expression of miR-138 has been observed in different types of cancers (39–43) but not HCC. MiR-138 plays an important role in tongue squamous cell carcinoma cell migration and invasion by concurrently targeting RhoC and ROCK2 (39). MiR-138 could inhibit the expression of HIF-1a and regulate the apoptosis and migration of clear cell renal cell carcinoma cells (40). MiR-138 enhanced cell migration and invasion by targeting enhancer of zeste homologue 2 (EZH2) in squamous cell carcinoma cell lines (41). MiR-138 may play an important role in cancer initiation and progression by regulation of FOSL1 in squamous cell carcinoma (42) or G protein alpha-inhibiting activity polypeptide 2 (GNAI2) in tongue squamous cell carcinoma (43). We identified CCND3 as a target of miR-138 in HCC, which may provide new insights into the mechanisms underlying tumorigenesis. CCND3 is expressed in nearly all proliferating cells and could promote the cell cycle start (44). Liu et al. (45) reported miR-16 family (including miR-16, miR-195 and miR-424) could induce cell cycle arrest by targeting CCND3, CCNE1 and cyclin-dependent kinase (Cdk)-6. In TLDA analysis, miR-16 and miR-195 showed no difference in all three pairs of HCC tissues and miR-424 showed significant difference only in two of three pairs of HCC tissues, which suggest miR-16 family may not play important role in HCC tumorigenesis. Thus, we didn’t verify the regulation of miR-16 family on CCND3. The cell cycle is regulated by a family of the CDKs and their activating partners (cyclins). The G1/S phase transition is regulated primarily by D-type cyclins (D1, D2 or D3) in complex with CDK4/CDK6 and E-type cyclins (E1 or E2) in complex with CDK2. These complexes cooperate in phosphorylating and preventing Rb binding to E2F, thus activating E2F-mediated transcription and driving cells from G1 into S phase (46). Whether the CDK/pRb/E2F pathway is involved in the miR-138 induction of cell cycle arrest needs further study. On the other hand, Wang et al. (47) showed that CCND3–cdk4/cdk6 specifically phosphorylate C/EBP at Ser193 and support growth-inhibitory C/EBP–cdk2 and C/EBP–Brm complexes in differentiated cells , but we showed the promotion of growth of CCND3 regulated by miR-138 in HCC cell lines. This discrepancy may be due to different cell lines and different subjects of study. CCND3 is expressed abroad, suggesting that it has many different cell functions. CCND3 had been reported to be involved in the promotion of cell growth by several groups (48–52). Under the specific study design, the researchers could investigate only one or some aspects of CCND3 but all of it, which is needed to be further investigated.
In summary, we report the expression patterns of altered miRNAs in HCC and the potential role of miR-138 in tumorigenesis. Our data suggest the frequently downregulated miR-138 can regulate CCND3 and function as a tumor suppressor in HCC. Therefore, miR-138 may serve as a useful therapeutic agent for miRNA-based HCC therapy.
This work was supported by Ministry of Science and Technology Key program (2008ZX10002-017 and 2012ZX10002009-004), Shanghai Leading Academic Discipline Project (B901) and Science Fund for Creative Research Groups, National Nature Science Foundation of China, China (30921006).
We thank research assistants at Orbital Instrument for their assistance in TDLA performance and Geminix for their assistance in GO and KEGG pathway analysis.
Conflict of Interest Statement: None declared.