- Open Access
A novel miR-0308-3p revealed by miRNA-seq of HBV-positive hepatocellular carcinoma suppresses cell proliferation and promotes G1/S arrest by targeting double CDK6/Cyclin D1 genes
Cell & Bioscience volume 10, Article number: 24 (2020)
Persistent infection with hepatitis B virus (HBV) accounts for the majority of hepatocellular carcinoma (HCC), but the molecular mechanisms underlying liver carcinogenesis are still not completely understood. Increasing evidence demonstrates that microRNAs (miRNAs) play significant functional roles in virus–host interactions. The aim of this study was to explore differentially expressed miRNA profiles and investigate the molecular mechanism of miR-0308-3p in HBV-positive HCC carcinogenesis.
High-throughput sequencing was used to detect novel miRNAs in three samples of HBV-positive HCC tissue compared to matched HBV-negative HCC tissue. The Cancer Genome Atlas (TCGA) database was used to mine miRNAs related to HBV-positive HCC. Bioinformatics analyses were conducted to predict the miRNAs’ possible biological and pathway regulatory functions. Quantitative polymerase chain reaction (qPCR) was then applied to evaluate the expression levels of randomly selected miRNAs. CCK-8 was used to measure cell proliferation and cell cycles were analyzed using flow cytometry. A dual luciferase reporter gene assay was used to confirm the downstream targets of miR-0308-3p.
In total, there were 34 overlapping miRNAs in both our miRNA-seq data and the TCGA database. We found two overlapping miRNAs in both the HBV-positive HCC samples and the TCGA database, and 205 novel pre-miRNA sequences were predicted. miR-522 and miR-523 were markedly overexpressed in HBV-positive HCC and were associated with a significantly poorer long-term prognosis (miR-522, HR 2.19, 95% CI 1.33–3.6, p = 0.0015; miR-523HR 1.5, 95% CI 1–2.44, p = 0.0047). Of note, we found that the novel miR-0308-3p was markedly downregulated in HBV-positive HCC samples and HCC cancer cell lines compared with HBV-negative HCC samples and adjacent normal hepatic tissue. Moreover, elevated expression of miR-0308-3p was found to inhibit proliferation of cancer cells by promoting G1/S cell cycle arrest but did not influence the apoptosis of cancer cells. A dual luciferase reporter activity assay identified that miR-0308-3p acted directly on the target sequence of the CDK6 and Cyclin D1 mRNA 3ʹUTR to suppress CDK6 and Cyclin D1 expression.
MiR-0308-3p upregulation dramatically suppressed HCC cell proliferation and induced G1/S cell cycle arrest by directly targeting CDK6/Cyclin D1. These findings reveal a novel molecular mechanism for activation of G1/S arrest in HCC and may prove clinically useful for developing new therapeutic targets.
Hepatocellular carcinoma (HCC) comprises the vast majority of primary liver cancers ; it is the sixth most prevalent cancer and third most frequent cause of cancer-associated deaths worldwide . The incidence of HCC has risen gradually, from 1.51 cases per 100,000 in 2011 to 6.20 cases per 100,000 worldwide in 2014 , and this increase is anticipated to continue into the future . According to the 2018 cancer statistics, there are 841,000 new cases of HCC worldwide each year, and about 55% of HCC patients are from China . Notably, approximately 83% of HCC cases arise in non-Western countries, with China accounting for slightly more than 50% of all new diagnoses in 2012. There are numerous established HCC risk factors, of which chronic hepatitis infection, in particular hepatitis B virus (HBV) infection, and subsequent liver fibrosis and cirrhosis are the most critical [6, 7]. It is estimated that approximately 2 billion people have been infected with HBV worldwide. Generally, HBV contributes to the development of HCC through integration of HBV DNA into the host genome, which induces both genomic instability and direct insertional mutagenesis of diverse cancer-related genes . Globally, approximately 240 million people are chronically infected with HBV, and approximately 25% of those affected eventually develop HCC. Moreover, approximately 60% of HCC cases in Africa and Asia are associated with HBV infection . Despite its prevalence, the mechanism for HBV-positive HCC is still poorly understood. A better understanding of the molecular mechanisms underlying HBV-positive HCC carcinogenesis would help in the development of more effective molecular-targeted interventions for the primary prevention and treatment of HCC.
MicroRNAs (miRNAs) are small non-coding RNAs that participate in post-transcriptional regulation of gene expression . Accumulating evidence implies that miRNAs play essential roles in the development and progression of HCC, and they are attracting increasing attention as a possible new class of key HCC biomarkers . Current research regarding the use of miRNAs as biomarkers for HCC can be categorized into three strategies. One is to use high-throughput sequencing to construct miRNA expression profiles for HCC tissue and adjacent normal tissue . The second is to mine currently available HCC-related RNA data in online RNA-seq databases, including The Cancer Genome Atlas (TCGA) and NCBI Gene Expression Omnibus (GEO) databases . The third is to use classical qRT-PCR assays to detect miRNAs in clinical HCC patient samples . These three strategies are commonly used to investigate key miRNAs, which may lead to the diagnosis and treatment of HCC at an earlier stage. miR-139 was recently investigated as a new biomarker for chronic HBV-positive HCC using microarrays . Similarly, Chen et al. detected expression of a number of miRNAs previously associated with lung cancer in HCC patients’ serum using qRT-PCR and discovered that miR-331-3p may be a new diagnostic and prognostic marker for HCC . Lin et al. conducted bioinformatics analyses using NCBI GEO and found a number of new miRNAs that were predicted to be associated with HCC development . Zhang et al. investigated miRNA profiles in 328 HCC patient samples from the TCGA database and generated an HCC-specific 7-miRNA signature that was validated as an independent prognostic biomarker . Liang L conducted bioinformatics analyses using both the GEO and TCGA databases, subsequently validating the previously reported miR-338-5p as a biomarker for diagnosing HCC, as well as discovering 423 of its potential target genes . Although rapid advances in RNA sequencing have led to an exponential increase in the prediction of miRNAs [18,19,20,21], and an enormous number of their downstream target genes, most of these have been generated from a single analysis strategy rather than integration of HCC clinical sequencing, HCC-related RNA-seq database analysis, and classical qRT-PCR. Furthermore, functional analyses regarding most newly discovered miRNAs in HBV-positive HCC samples are limited to gene ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, or survival analysis, and further investigation of their molecular mechanisms through cell experimentation is seldom conducted. Li et al. integrated microarray and TCGA database mining to investigate the role of miR-200b in renal cell carcinoma carcinogenesis and demonstrated its previously undescribed role as a suppressor of tumor metastasis by directly destabilizing Laminin subunit alpha 4 (LAMA4 mRNA) . Hence, we hypothesized that integration of high-throughput sequencing of clinical HBV-positive HCC samples and RNA-seq database mining would lead to the discovery of more biomarkers and a better understanding of the molecular-mechanisms in HBV-positive HCC carcinogenesis.
Methods and materials
Patient information, TCGA introduction, cell lines and plasmids construction and tranfection
In total, 15 pairs of fresh specimens containing human HBV-positive and HBV-negative HCC tissue were collected from the First Affiliated Hospital of the University of South China. The baseline information of the enrolled patients is shown in Table 1. All the enrolled HCC patients’ diagnoses were based on case histories, clinical symptoms, predictors, and/or tissue biopsy according to the National Diagnosis Criteria. After the samples were collected, they were sent to RiboBio Co., Ltd. (Guangzhou, China) for immediate high-throughput sequencing analysis. This study was performed in accordance with the Declaration of Helsinki and was approved by the institutional review board (IRB) of the First Affiliated Hospital of the University of South China (IRB approval no. NHUH-2017-38). All enrolled patients signed informed-consent forms.
The HCC-related miRNA-seq and RNA-seq data were downloaded from the TCGA database available at http://portal.gdc.cancer.gov/. The sequencing data were all publicly available, therefore no ethical issues were involved .
The human HepG2 and SMMC-7721 cell lines were originally purchased from ATCC. Both cell lines were maintained in Dulbecco’s modified Eagle medium (DMEM; GIBCO BRL, Grand Island, NY) supplemented with 10% FBS and 1% PS (100 units/ml penicillin, 100 μg/ml streptomycin) and incubated at 37 °C with 5% CO2.
The CDK6 and Cyclin D1 without 3ʹ-UTR was synthesized by Beijing Genomics Institute (Beijing, China) and subcloned into the pEX-3 vector (GenePharma, Shanghai, China) to generate the miR-0308-3p-resistant pEX-3-CDK6 vector and Cyclin D1 following the manufacturer’s instructions. An empty pEX-3 vector was used as a control. miR-0308-3p-resistant pEX-3-CDK6 vector, miR-0308-3p-resistant pEX-3- Cyclin D1 vector and control plasmid were transfected using Lipofectamine 3000 Reagent based on the agent-box instruction.
High-throughput sequencing was conducted using an Illumina HiSeq 2500 (Illumina, Inc., San Diego, CA, USA) followed by DNA library construction. For construction of the DNA library, first total RNA was isolated from HCC samples with or without HBV infection using an miRNA isolation kit (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions. Then, RNA molecules 18–30 nt in size were separated by gel electrophoresis and 3ʹ and 5ʹ adapters were ligated to the RNAs. The ligation products were constructed by reverse transcription polymerase chain reaction (RT-PCR), and PCR products 140–160 bp in size were enriched to generate a cDNA library . After sequencing, the reads were filtered to achieve clean tags. All the clean tags were aligned with small RNAs using miRBase version 21 (http://www.mirbase.org), Rfam12.1 (http://rfam.xfam.org), and piRNAbank (http://pirnabank.ibab.ac.in) to identify known miRNAs. For novel miRNA candidate prediction, the following selection criteria were applied using Mireap v0.2 software: removal of isomers, hairpin structure, stable secondary structure with binding free energy less than − 20 kcal/mol, and located at an intergenic spacer.
Quantitative reverse-transcription polymerase chain reaction (qRT-PCR)
For analysis of miRNA expression, qRT-PCR was used. Briefly, total RNA was isolated from HCC samples using an miRNA isolation kit (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s instructions. Then, 1 µg of total RNA was reverse-transcribed. SYBR Premix Ex Taq II (Takara, Tokyo, Japan) was used for RT-qPCR analysis, according to the manufacturer’s instructions. The relative expression of mRNA or miRNA was calculated using the 2−ΔΔCt method and normalized to the expression of GAPDH or U6, respectively. All PCRs were performed in triplicate. The miRNA primer sets were purchased from GenePharma (Suzhou, China). The primer pairs of the relevant genes are shown in Additional file 1: Table S1.
Western blotting analysis, flow cytometry analysis, and apoptosis detection
Western blotting was used to evaluate levels of protein expression according to our previously published studies [25,26,27]. The cell-cycle distribution of HepG2 was estimated using flow cytometry analysis. HepG2 and SMMC-7721 cells were placed in 12-well plates at a density of 2 × 105 cells/well and transfected with or without an miR-0308-3p mimic followed by 24 h incubation. Dimethyl sulfoxide (DMSO) was used as a control. For estimation of DNA content, the cells were washed in PBS and fixed in ethanol at – 20 °C. This was followed by re-suspension in PBS with 40 μg/ml PI and RNase A (0.1 mg/ml) and Triton X-100 (0.1%) for 30 min in a dark room at 37 °C .
For the detection of apoptotic cell populations, HepG2 and SMMC-7721 cells at a density of 2 × 105 cells/well were seeded in six-well plates and treated with or without an miR-0308-3p mimic for 24 h. The cells were then stained with 4ʹ-6-diamidino-2-phenylindole (DAPI) and imaged under a fluorescence microscope. Annexin V/IP was used to estimate apoptotic cell populations.
For the HepG2 cell-proliferation analysis, HepG2 cells were collected at passage 3–4 and inoculated in six-well plates at a density of 4 × 106 cells/well. The effects of miR-0308 on cell viability were detected using a standard cell-counting kit-8 (CCK-8) according to the manufacturer’s instructions. The optical density (OD) of the cells in each group was tested by measuring the absorbance at 570 nm using a microplate reader .
miRNA target prediction, Gene Ontology (GO) analysis, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis
miRNA target prediction was conducted using miRDB (http://mirdb.org/). For target prediction of known miRNAs, the miRNA name was entered directly as the Target search, while for target prediction of novel miRNAs, the RNA sequence was entered in the Custom Prediction search.
The Gene Ontology (GO) database is a consortium-based dataset that offers information on gene products regarding cellular components, molecular function, and biological processes using ontologies. There are approximately 500,000 annotations about Homo sapiens in the database. Predicted targets for a miRNA are evaluated as a group using GO terms. Significant GO categories are identified using statistics for gene-function enrichment .
KEGG was developed by Kanehisa Laboratories in 1995 and is now a widely available bioinformatics database of the high-level functions and utilities of biological systems from the perspective of genomes and molecules (https://www.genome.jp/kegg/). The current version consists of 18 databases, including KEGG Pathway, KEGG Genome, and KEGG Drug. In this study, we performed GO and KEGG analyses on the significantly differentially expressed miRNAs, novel miRNAs, and miRNAs with overlapping expression patterns.
Dual luciferase reporter assay
The wild-type CDK6 mRNA 3ʹ-untranslated region (3ʹUTR) (Genbank: NM001145306.2), Cyclin D1 3'UTR (Genbank:BC023620.2) and the mutant sequences at the predicted target sites for miR-0308-3p in the CDK6 mRNA 3ʹUTR were cloned into the pmirGLO vector to generate the pmirGLO-CDK6_3ʹUTR_wt, pmirGLO-Cyclin D1_3'UTR_wt, pmirGLO-CDK6_3ʹUTR-Mut and pmirGLO-Cyclin D1_3'UTR_Mut constructs, respectively. HepG2 cells were seeded onto 24-well plates (6.0 × 104 cells/well) for 24 h before co-transfection with 1 μg of the reporter plasmids and 1 μg of the pmirGLO internal control plasmids. After 8 h, the medium was replaced and the cells were transfected with 100 nM miR-0308 mimic or control (DMSO). After incubation for 48 h, the transfected cells were lysed and their luciferase activity was detected using the Dual-Luciferase Reporter Assay System (Cat. No. E1910; Promega, Madison, WI, USA). Firefly luciferase activity was normalized to that of Renilla luciferase and each group was detected in triplicate according to previous reports [29, 31].
A Kaplan–Meier survival analysis was performed on the available Kaplan–Meier plotter website (http://www.kmplot.com/analysis/), which is capable of assessing the effects of 54,000 genes on survival in 21 types of cancer, including HCC . The following five miRNAs with overlapping expression patterns were used in the survival analysis: miR-523, miR-517a, miR-372, miR-522, and miR-524-5p.
Statistical analyses were performed using SPSS software (ver. 19.0; IBM Corp., Armonk, NY, USA). Quantitative data were expressed as mean ± standard deviation. A one-way analysis of variance (ANOVA) followed by the least significant difference post-hoc test and t-test were used for the comparison of means. We considered p-values < 0.05 to be statistically significant.
HCC patients with HBV-positive infection showed global changes in miRNA expression
Figure 1 indicates the study design. To compare miRNA expression profiles between HCC patients with or without HBV infection and to mine possible biomarkers for HCC patients with HBV-positive infection, five HBV-positive HCC samples, five HBV-negative HCC samples, and the adjacent non-tumor tissue were used. In total, 2588 miRNAs were identified, of which 775 were differentially expressed with a fold-change of more than 1. Figure 2a shows the 148 upregulated and 627 downregulated miRNAs in the HBV-positive HCC samples, and Table 2 shows the top 20 upregulated and downregulated miRNAs from Fig. 2. The most highly upregulated miRNA was miR-4686 with an approximately sixfold change (p < 0.01). Three miRNAs—miR-190a-3p, miR-4662b, and miR-2681-3phad a fold-change above four. Among the downregulated miRNAs, miR-520d-3p and miR-372-3p exhibited the highest fold-changes in the HBV-positive HCC samples (p < 0.01).
To predict the biological functions of these differentially expressed miRNAs, a GO analysis was conducted (Fig. 2b). Differentially expressed miRNAs were predicted to have the following main functions: cell-junction, cell-part, and extracellular-matrix functions in the cellular-component category; and catalytic-activity, channel-regulator-activity, and metallochaperone-activity functions in the molecular-function category. The heatmap and volcano plot illustrated that compared with HBV-negative HCC patients, the numbers of downregulated miRNAs in HBV-positive HCC patients were higher than the numbers of upregulated miRNAs (Fig. 2c, d). We then performed a KEGG analysis (Fig. 2e), which showed that the main target pathways for the differentially expressed miRNAs were purine metabolism, oxytocin signaling, miRNAs in cancer pathways, HIF-1 signaling, and focal adhesion. Taken together, the GO and KEGG data imply that HBV-infection changes miRNA expression profiles in HCC tissue. Furthermore, it may be possible to use this profile change to predict specific inflammation and tumorigenesis changes observed in HCC tissue.
Novel miRNA prediction in HBV-positive HCC samples
To screen for novel miRNAs, high-throughput sequencing was again used. Based on the results, most of the putative pre-miRNAs had a typical hairpin structure, with a high read depth in the mature arm and a much lower read depth in the other arm. In total, 205 novel pre-miRNA sequences were predicted, each derived from a unique genome location (Additional file 1: Table S1). Table 3 lists 25 novel miRNAs, together with the details of each sequence and read numbers. We classified the novel miRNAs based on their different chromosomal locations; novel miRNAs were found on all chromosomes except the Y chromosome (Fig. 3a). The high-throughput sequencing data showed that the length of the novel miRNAs ranged from 60 to 101 nt, and the majority were between 81 and 90 nt (Fig. 3a). The precursors of these potential novel miRNAs formed proper secondary hairpin structures with free energies ranging from − 20.60 to − 70.20 kcal/mol (Additional file 1: Table S1). Additional file 1: Table S1 also presents the probability that the predicted novel miRNAs were derived from the genomic spacer of a gene.
To predict the biological functions of the novel miRNAs, the same strategy as for the differentially expressed miRNAs was employed. The GO analysis predicted the novel miRNAs to have the following main functions: biological-adhesion, cell-killing, and cellular process functions in the biological-process category; cellular-structure, extracellular-matrix, and cell-membrane functions in the cellular-component category; and antioxidant-activity, catalytic-activity, and enzyme-regulator-activity functions in the molecular-function category (Fig. 3b). Next, the KEGG analysis showed that the main target pathways for the novel miRNAs were Wnt signaling, oxytocin signaling, MAPK signaling, miRNAs in cancer, HIF-1 signaling, and calcium signaling (Fig. 3c).
Compared to the differentially expressed miRNAs, these novel miRNAs appeared to have exceptional functions such as cell killing and pathway-specific signaling. The data also indicated that the HBV-positive HCC samples contained more novel miRNAs compared to the HBV-negative HCC samples. The presence of these novel miRNAs may predict their involvement in the physiological regulation of HBV-positive HCC.
Identification of miRNAs with overlapping expression patterns in both the miRNA-seq data and the TCGA database data
To explore miRNAs with overlapping expression patterns, we used both our miRNA-seq data and RNAseq data from the TCGA database related to HCC with or without HBV infection. All data were obtained in accordance with relevant laws, and any necessary approvals and informed-consent documents were obtained . In total, 34 miRNAs with overlapping expression patterns were identified in both the miRNA-seq data and TCGA database data (Fig. 4a). Table 4 lists the 20 miRNAs with the greatest overlap. Most of the miRNAs with overlapping expression patterns were downregulated in HBV-positive HCC samples, except for miR-520d-5p and miR-518e-5p, which were upregulated. Downregulated miR-520d-3p had the highest fold change, of 11.03, followed by downregulated miR-520d-5p and miR-520a-3p with fold changes of 8.99 and 8.42, respectively. The heatmap further demonstrated that most of these miRNAs were typically downregulated in the HBV-positive HCC samples (Fig. 4b). The miRNAs’ main target pathways were transcriptional misregulation in cancer, purine metabolism, endocytosis calcium signaling pathway, and hematopoietic cell lineages (Fig. 4c). The overlapping miRNAs were predicted to play specific roles in inflammation and carcinogenesis in HBV-positive HCC.
Validation of novel miRNAs and miRNAs with overlapping expression patterns using qRT-PCR in 15 pairs of HCC samples with or without HBV infection
To evaluate the miRNAseq and TCGA mining results, we conducted qRT-PCR to verify the expression patterns of several miRNAs selected at random. For validation of the novel miRNAs, we selected the following: miR-0050-3p, miR-0142-3p, miR-0157-3p, miR-0308-3p and miR-0233-5p. Compared to the HBV-negative HCC samples, miRN-0050-3p, miR-0233-5p and miR-0308-3p were significantly downregulated, whereas miR-0142-3p and miR-0157-3p were downregulated but this did not reach significance (Fig. 5a).
For validation of the miRNAs with overlapping expression patterns, we selected miR-517a-3p, miR-372-5p, miR-524-5p, miR-523, and miR-522-5p. Compared to the HBV-negative HCC samples, miRN-517a-3p, miR-523 and miR-522-5p were significantly upregulated, whereas miR-372-5p and miR-524-5p were significantly downregulated (Fig. 5b). In Table 4, the miR-517a-3p expression was downregualted in HBV-negative HCC samples whereas upregualted in the validation test. The possible explain may be associated with heterogeneity of the HCC sample which needs further investigation.
To evaluate the clinical significance of these five miRNAs with overlapping expression patterns further, we performed a Kaplan–Meier analysis of overall survival (OS) rates based on the TCGA database data. The results showed that HBV-positive HCC patients with high miR-523, miR-517a, and miR-524 expression had a poorer long-term prognosis, although this finding was significant only for miR-522 (HR 2.19, 95% CI 1.33–3.6, p = 0.0015) and miR-523 (HR 1.5, 95% CI 1–2.44, p = 0.0047) (Fig. 5c).
miR-0308-3p expression may play a role in regulation of the HCC cell cycle
To explore the role of novel miRNAs in HBV-positive HCC further, we performed miRNA-expression, cell-proliferation, cell-cycle, and apoptosis assays on HepG2 and SMMC-7721 cell lines. We then selected miR-0308-3p to investigate its effect on HepG2 cell proliferation using the CCK-8 assay. Following transfection with the miR-0308 mimic, HepG2 cell proliferation decreased significantly compared to cells transfected with NC (p < 0.05, Fig. 6a). We also investigated the effect of miR-0308-3p on the cell cycles of HepG2 and SMMC-7721 cell lines. Following transfection with the miR-0308 mimic, G1 cell numbers significantly increased in both cell lines compared to cells transfected with NC (p < 0.05, Fig. 6b–e).
Next, the effect of miR-0308-3p on cell apoptosis was investigated, but no significant effect was found in either HepG2 or SMMC-7721 (p > 0.05, Fig. 6f–i). Together, these results imply that increased miR-0308-3p expression contributes to the activation of G1/S arrest in HCC cancer cells.
miR-0308-3p downregulates double CDK6/Cyclin D1 expression by directly targeting the CDK6/Cyclin D1 3′-UTR
To uncover the underlying mechanism of miR-0308-3p in activation of HCC G1/S cell arrest, we performed qRT-PCR and western blot analysis. The results showed that the levels of mRNA and protein expression related to double CDK6/Cyclin D1 decreased gradually as the concentration of the miR-0308 mimic increased from 25 to 100 nM (Fig. 7a, b).
To explore the direct action of miR-0308-3p on the target CDK6 mRNA 3ʹUTR, pmirGLO-CDK6_3ʹUTR_wt and pmirGLO-CDK6_3ʹUTR-Mut reporter vectors were constructed. The sequences of the 3ʹUTR region of CDK6 mRNA containing the miR-0308-3p targeting sequence and its mutant were cloned separately, downstream of the pmirGLO plasmid luciferase gene, to construct the reporter vectors. The predicted miR-0308 targeting site and the sequence of the CDK6 mRNA 3ʹUTR are shown in Fig. 7c. The cells were then co-transfected with these reporter vectors and either miR-0308-3p or miR-NC. miR-0308-3p significantly decreased the luciferase activity of the pmirGLO-CDK6_3ʹUTR reporter, but had no effect on the pmirGLO-CDK6_3ʹUTR mutant (Fig. 7d). This result implies that miR-0308-3p acts directly on the target sequence of the CDK6 mRNA 3ʹUTR to suppress CDK6 expression.
The same method was employed to explore the direct action of miR-0308-3p on the target Cyclin D1 mRNA 3ʹUTR. First pmirGLO-Cyclin D1_3ʹUTR_wt and pmirGLO-Cyclin D1_3ʹUTR-Mut reporter vectors were constructed in the same manner as those for CDK6 (Fig. 7e). Cells were co-transfected with these reporter vectors and either miR-0308-3p or miR-NC, and miR-0308-3p significantly decreased the luciferase activity of the pmirGLO-Cyclin D1_3ʹUTR reporter (Fig. 7f), but had no effect on the pmirGLO-Cyclin D1_3ʹUTR mutant. This result implies that miR-0308-3p acts directly on the target sequence of the Cyclin D1 mRNA 3ʹUTR to suppress Cyclin D1 expression.
To further role out that downregualtion of CDK1 and CYClin D1 by miR-0308-3p could play a important role to mediate miR-0308-3p-induced G1/S cell cycle arrest, the expression vectors for miR-0308-3p-resistance CDK6 and Cyclin D1 vectors were constructed and trnasfected in the HepG2 cells, as Additional file 2 :Figure S1 shown, compared to the control group, G1 cell numbers in cells treated with miR-0308-3p-resistant pEX-3- Cyclin D1 vector or CDK1 vector were significantly rescued (p < 0.05).
Microarray analysis and high-throughput analysis have enabled tens of thousands of miRNAs to be identified over the past few decades . More recently, studies to identify cancer-related miRNAs and predict or evaluate their functions have been increasing rapidly . Several miRNAs have also been found to play essential roles in the carcinogenesis of HCC; [36,37,38,39] however the number of miRNAs characterized in HBV-positive HCC is still limited and worth further investigation. In light of this, we used high-throughput sequencing and re-analyzed the RNA-seq data from the TCGA database to predict new and potential miRNAs associated with the carcinogenesis of HBV-related HCC. Notably, we found two miRNAs, miR-522 and miR-523, with overlapping expression patterns in both our HBV-positive HCC samples and the TCGA database. These miRNAs were markedly overexpressed in HBV-positive HCC and were associated with a significantly poorer long-term prognosis (miR-522, HR 2.19, 95% CI 1.33–3.6, p = 0.0015; miR-523, HR 1.5, 95% CI 1–2.44, p = 0.0047).
Of note, we found a novel miR-0308-3p that was markedly downregulated in HBV-positive HCC samples compared to HBV-negative HCC samples and adjacent normal hepatic tissue. Moreover, we found that miR-0308-3p was downregulated in HepG2 cells and 7721 cancer cells, and that elevated expression of miR-0308-3p inhibited proliferation and promoted G1/S arrest but did not influence the apoptosis of cancer cells. This implies that miR-0308-3p may affect the carcinogenesis of HBV-related HCC through regulating G1/S-related genes. Based on this inference, we used bioinformatics to predict that two G1/S-related genes, CDK6 and Cyclin D1, were downstream targets of miR-0308-3p. Furthermore, a dual luciferase reporter activity assay identified that miR-0308-3p acted directly on the target sequence of the CDK6/Cyclin D1 mRNA 3ʹUTR to suppress CDK6/Cyclin D1 expression.
Many reports have revealed that miRNAs are involved in human carcinogenesis. HBV-positive HCC has been found to be associated with altered expression of miRNAs. Wang et al. investigated the miRNA profiles of 12 pairs of chronic HBV-associated HCC samples using microarray analysis  and found seven miRNAs (miR-150, miR-342-3p, miR-663, miR-20b, miR-92a-3p, miR-376c-3p and miR-92b) that were specifically altered in HBV-related HCC. Nielsen et al. screened for alterations in miRNA expression in liver-cancer cell lines post-HBV infection and found that miR-192-5p, miR-194-5p, and miR-215-5p were upregulated, providing further evidence that these miRNAs are associated with the development of HBV-related liver disease . Lou et al. mined the NCBI GEO database for HBV-related HCC RNAseq data and found seven upregulated and nine downregulated differentially expressed miRNAs . Our study integrated high-throughput sequencing and TCGA database mining and found further alterations in miRNA expression in HBV-positive versus HBV-negative HCC. Notably, we identified a series of 34 miRNAs with overlapping expression patterns. miRNAs with overlapping expression patterns such as miR-522, miR-518, miR-519, and miR-520 have previously been reported to be upregulated in HBV-positive HCC . In contrast, our results found miR-522, miR-518, miR-519, and miR-520 to be downregulated in HBV-positive HCC. This may have been due to the heterogeneity and plasticity of clinical samples, or the expression of these miRNAs may be dependent on internal or external factors.
After differentially expressed key miRNAs were identified from clinical HCC samples, we performed GO and KEGG bioinformatic analyses to predict the miRNAs’ potential biological functions and signaling pathways. miRNAs are usually used to conduct OS assays to elucidate their prognostic value in cancer patients. Wang et al. analyzed miR-522 expression in HBV-positive HCC tissue using GO and KEGG analyses and found that regulation of the actin cytoskeleton and pathway in cancer are most enrichment role in HCC carcinogenesis . Shi et al. indicated that miR-522 overexpression in HCC tissue was markedly correlated with a poor prognosis and decreased OS in HCC patients . miR-523 has also been reported to be upregulated in acute myeloid leukemia plasma , pituitary adenomas , colon cancer , and uveal melanom , but has not previously been reported to be associated with HCC. In our study, we found miR-523to be upregulated in HBV-positive HCC carcinogenesis. Since miR-523 has not been investigated in HBV-positive HCC tissue previously, our results imply that overexpression of miR-523 may serve as a novel biomarker for an unfavorable prognosis in HBV-positive HCC patients.
Exploring the mechanism underlying the carcinogenesis of HBV is critical in prevention of, and therapeutics for, HBV-positive HCC. In recent years, studies focusing on the roles of miRNAs in HBV-positive HCC have increased. A study by Liu et al. demonstrated that HBV preS2 promotes the transcriptional co-activator with PDZ binding motif (TAZ) through miR-338-3p leading to enhanced HCC proliferation and migration . Shi et al. indicated that miR-200a-3p was downregulated by HBV X protein and promoted liver cancer cell proliferation and invasion . Another study found that miR-384 was downregulated in HBV-related HCC tissue and promoted liver cancer cell proliferation and metastasis by targeting the oncogene pleiotrophin . Meanwhile, miR-34c was shown to target TGFβ-induced factor homeobox 2 and repress cell proliferation in HBV-related HCC . Xiang et al. found that upregulation of miR-499a induced HCC carcinogenesis by targeting MAPK6 . Qin et al. illustrated that miR-30b-5p repressed cell proliferation and arrested the cell cycle of a HCC cell line by targeting USP37 . Lamontagne et al. pointed out in a review that altered regulation of miRNA expression may play a critical role in HBV-related cancer development .
In general, alteration of miRNA expression leads to altered expression of its targeted genes, which may influence carcinogenesis. For instance, miR-218 can induce cell-cycle arrest at the G1 phase by targeting the 3ʹ-UTR region of CDK6/Cyclin D1 , and overexpression of miR-24 increases the percentage of cells in G1 phase in HepG2 and K562 cancer cells . Cell-cycle deregulation is a common feature of HCC, as cancer cells frequently display unscheduled proliferation resulting from disruptions in the cell cycle, such as G1/S arrest . CDK6/Cyclin D1 genes have been reported to be induced during G1/S transition, and G1/S dysregulation is usually mediated by alterations in the activity of CDK6/Cyclin D1. Although there exists limitations in this study including the evidence for that the miR-0308-3p overexpressing HCC are less aggressive and patients have better survival rate remains further investigated, as well as for such large scale screening of miRNA, the cut off as p < 0.01 may not appropriate, we found a novel miRNA, miR-0308-3p, that inhibited cancer cell proliferation and induced G1/S arrest via targeting double CDK6/Cyclin D1 genes. No miRNAs have been previously reported to suppress G1/S transition in HBV-positive HCC, and our results imply that miR-0308-3p may act as an HBV-positive HCC suppressor.
In summary, this study demonstrated the effect of miR-0308-3p on HBV-positive HCC. MiR-0308-3p upregulation dramatically suppressed HCC cell proliferation and induced G1/S arrest by directly targeting CDK6/Cyclin D1. These findings reveal a novel molecular mechanism of activation of G1/S arrest in HCC and may prove clinically useful for developing new therapeutic targets for HCC. However, the mechanism underlying HBV regulation of miR-0308-3p expression requires further study, and the effect of miR-0308-3p as well as the other newly identified miRNAs and their targets should be verified with further in vitro and in vivo experiments. Furthermore, the molecular mechanism underlying HBV-positive carcinogenesis of liver cells requires additional future studies. Answering these questions may be of significance for the prevention and treatment of HBV-positive HCC.
Availability of data and materials
Wallace MC, Preen D, Jeffrey GP, Adams LA. The evolving epidemiology of hepatocellular carcinoma: a global perspective. Expert Rev Gastroenterol Hepatol. 2015;9:765–79.
Forner A, Llovet JM, Bruix J. Hepatocellular carcinoma. Lancet. 2012;379:1245–55.
Njei B, Rotman Y, Ditah I, Lim JK. Emerging trends in hepatocellular carcinoma incidence and mortality. Hepatology. 2015;61:191–9.
Kim HS, El-Serag HB. The epidemiology of hepatocellular carcinoma in the USA. Curr Gastroenterol Rep. 2019;21:17.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2018. CA Cancer J Clin. 2018;68:7–30.
Yan SY, Fan JG, Qio L. Hepatitis B virus (HBV) infection and hepatocellular carcinoma—new insights for an old topic. Curr Cancer Drug Targets. 2017;17:505–11.
Neureiter D, Stintzing S, Kiesslich T, Ocker M. Hepatocellular carcinoma: therapeutic advances in signaling, epigenetic and immune targets. World J Gastroenterol. 2019;25:3136–50.
Levrero M, Zucman-Rossi J. Mechanisms of HBV-induced hepatocellular carcinoma. J Hepatol. 2016;64:S84–101.
Xie KL, Zhang YG, Liu J, Zeng Y, Wu H. MicroRNAs associated with HBV infection and HBV-related HCC. Theranostics. 2014;4:1176–92.
Chen EB, Zhou ZJ, Xiao K, Zhu GQ, Yang Y, Wang B, Zhou SL, et al. The miR-561-5p/CX3CL1 signaling axis regulates pulmonary metastasis in hepatocellular carcinoma involving CX3CR1(+) natural killer cells infiltration. Theranostics. 2019;9:4779–94.
Lin LJ, Lin Y, Jin Y, Zheng CQ. Investigation of key microRNAs associated with hepatocellular carcinoma using small RNA-seq data. Mol Biol Rep. 2014;41:4341–9.
Shao P, Sun D, Wang L, Fan R, Gao Z. Deep sequencing and comprehensive expression analysis identifies several molecules potentially related to human poorly differentiated hepatocellular carcinoma. FEBS Open Bio. 2017;7:1696–706.
Tsuchiya N, Sawada Y, Endo I, Saito K, Uemura Y, Nakatsura T. Biomarkers for the early diagnosis of hepatocellular carcinoma. World J Gastroenterol. 2015;21:10573–83.
Li T, Yin J, Yuan L, Wang S, Yang L, Du X, Lu J. Downregulation of microRNA-139 is associated with hepatocellular carcinoma risk and short-term survival. Oncol Rep. 2014;31:1699–706.
Chen L, Chu F, Cao Y, Shao J, Wang F. Serum miR-182 and miR-331-3p as diagnostic and prognostic markers in patients with hepatocellular carcinoma. Tumour Biol. 2015;36:7439–47.
Zhang J, Chong CC, Chen GG, Lai PB. A seven-microRNA expression signature predicts survival in hepatocellular carcinoma. PLoS ONE. 2015;10:e0128628.
Liang L, Gao L, Zou XP, Huang ML, Chen G, Li JJ, Cai XY. Diagnostic significance and potential function of miR-338-5p in hepatocellular carcinoma: a bioinformatics study with microarray and RNA sequencing data. Mol Med Rep. 2018;17:2297–312.
Li Z, Rana TM. Therapeutic targeting of microRNAs: current status and future challenges. Nat Rev Drug Discov. 2014;13:622–38.
Huang R, Xiang J, Zhou P, Vitamin D. Gut microbiota, and radiation-related resistance: a love-hate triangle. J Exp Clin Cancer Res. 2019;38:493.
Shi F, Liu Y, Li M, Wen P, Qian QQ, Fan Y, Huang R. Analysis of lncRNA and mRNA transcriptomes expression in thyroid cancer tissues among patients with exposure of medical occupational radiation. Dose Response. 2019;17:1559325819864223.
Liu X, Zhou Y, Wang S, Guan H, Hu S, Huang R, Zhou P. Impact of low-dose ionising radiation on the composition of the gut microbiota of mice. Toxicol Sci. 2019;171:258–68.
Li Y, Guan B, Liu J, Zhang Z, He S, Zhan Y, Su B, et al. MicroRNA-200b is downregulated and suppresses metastasis by targeting LAMA4 in renal cell carcinoma. EBioMedicine. 2019;44:439–51.
Zeng JH, Xiong DD, Pang YY, Zhang Y, Tang RX, Luo DZ, Chen G. Identification of molecular targets for esophageal carcinoma diagnosis using miRNA-seq and RNA-seq data from The Cancer Genome Atlas: a study of 187 cases. Oncotarget. 2017;8:35681–99.
Wu Y, Ou Y, Liao C, Liang S, Wang Y. High-throughput sequencing analysis of the expression profile of microRNAs and target genes in mechanical force-induced osteoblastic/cementoblastic differentiation of human periodontal ligament cells. Am J Transl Res. 2019;11:3398–411.
Huang R, Yu T, Li Y, Hu J. Upregulated has-miR-4516 as a potential biomarker for early diagnosis of dust-induced pulmonary fibrosis in patients with pneumoconiosis. Toxicol Res. 2018;7:415–22.
Qin J, Ning H, Zhou Y, Hu Y, Yang L, Huang R. LncRNA MIR31HG overexpression serves as poor prognostic biomarker and promotes cells proliferation in lung adenocarcinoma. Biomed Pharmacother. 2018;99:363–8.
Liu XD, Xie DF, Wang YL, Guan H, Huang RX, Zhou PK. Integrated analysis of lncRNA-mRNA co-expression networks in the alpha-particle induced carcinogenesis of human branchial epithelial cells. Int J Radiat Biol. 2019;95:144–55.
Zhang Q, Bao J, Yang J. Genistein-triggered anticancer activity against liver cancer cell line HepG2 involves ROS generation, mitochondrial apoptosis, G2/M cell cycle arrest and inhibition of cell migration. Arch Med Sci. 2019;15:1001–9.
Mo LJ, Song M, Huang QH, Guan H, Liu XD, Xie DF, Huang B, et al. Exosome-packaged miR-1246 contributes to bystander DNA damage by targeting LIG4. Br J Cancer. 2018;119:492–502.
Zhang Z, Li B, Xu P, Yang B. Integrated whole transcriptome profiling and bioinformatics analysis for revealing regulatory pathways associated with quercetin-induced apoptosis in HCT-116 cells. Front Pharmacol. 2019;10:798.
Zhang B, Shetti D, Fan C, Wei K. miR-29b-3p promotes progression of MDA-MB-231 triple-negative breast cancer cells through downregulating TRAF3. Biol Res. 2019;52:38.
Menyhart O, Nagy A, Gyorffy B. Determining consistent prognostic biomarkers of overall survival and vascular invasion in hepatocellular carcinoma. R Soc Open Sci. 2018;5:181006.
Cancer Genome Atlas Research N. Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature. 2008;455:1061–8.
Weidle UH, Epp A, Birzele F, Brinkmann U. The functional role of prostate cancer metastasis-related micro-RNAs. Cancer Genomics Proteomics. 2019;16:1–19.
Huang SW, Ali ND, Zhong L, Shi J. MicroRNAs as biomarkers for human glioblastoma: progress and potential. Acta Pharmacol Sin. 2018;39:1405–13.
Nie YH, Liu XD, Huang R, Xie DF, Yin WJ, Guan H, Yu ZJ, et al. Analysis of mRNA expression patterns in peripheral blood cells of 3 patients with cancer after the first fraction of 2 Gy irradiation: an integrated case report and systematic review. Dose Response. 2019;17:1559325819833474.
Huang R, Zhou P. Double-edged effects of noncoding RNAs in responses to environmental genotoxic insults: perspectives with regards to molecule-ecology network. Environ Pollut. 2019;247:64–71.
Zhou PK, Huang RX. Targeting of the respiratory chain by toxicants: beyond the toxicities to mitochondrial morphology. Toxicol Res. 2018;7:1008–11.
Sun Y, Ji F, Kumar MR, Zheng X, Xiao Y, Liu N, Shi J, et al. Transcriptome integration analysis in hepatocellular carcinoma reveals discordant intronic miRNA-host gene pairs in expression. Int J Biol Sci. 2017;13:1438–49.
Wang G, Dong F, Xu Z, Sharma S, Hu X, Chen D, Zhang L, et al. MicroRNA profile in HBV-induced infection and hepatocellular carcinoma. BMC Cancer. 2017;17:805.
Nielsen KO, Jacobsen KS, Mirza AH, Winther TN, Storling J, Glebe D, Pociot F, et al. Hepatitis B virus upregulates host microRNAs that target apoptosis-regulatory genes in an in vitro cell model. Exp Cell Res. 2018;371:92–103.
Lou W, Liu J, Ding B, Chen D, Xu L, Ding J, Jiang D, et al. Identification of potential miRNA-mRNA regulatory network contributing to pathogenesis of HBV-related HCC. J Transl Med. 2019;17:7.
Wang W, Zhao LJ, Tan YX, Ren H, Qi ZT. Identification of deregulated miRNAs and their targets in hepatitis B virus-associated hepatocellular carcinoma. World J Gastroenterol. 2012;18:5442–53.
Shi YH, Qi BB, Liu XB, Ding HM. Upregulation of miR-522 is associated with poor outcome of hepatocellular carcinoma. Eur Rev Med Pharmacol Sci. 2016;20:3194–8.
Fayyad-Kazan H, Bitar N, Najar M, Lewalle P, Fayyad-Kazan M, Badran R, Hamade E, et al. Circulating miR-150 and miR-342 in plasma are novel potential biomarkers for acute myeloid leukemia. J Transl Med. 2013;11:31.
Liang S, Chen L, Huang H, Zhi D. The experimental study of miRNA in pituitary adenomas. Turk Neurosurg. 2013;23:721–7.
Bobowicz M, Skrzypski M, Czapiewski P, Marczyk M, Maciejewska A, Jankowski M, Szulgo-Paczkowska A, et al. Prognostic value of 5-microRNA based signature in T2-T3N0 colon cancer. Clin Exp Metastasis. 2016;33:765–73.
Russo A, Caltabiano R, Longo A, Avitabile T, Franco LM, Bonfiglio V, Puzzo L, et al. Increased levels of miRNA-146a in serum and histologic samples of patients with uveal melanoma. Front Pharmacol. 2016;7:424.
Liu P, Zhang H, Liang X, Ma H, Luan F, Wang B, Bai F, et al. HBV preS2 promotes the expression of TAZ via miRNA-338-3p to enhance the tumorigenesis of hepatocellular carcinoma. Oncotarget. 2015;6:29048–59.
Shi T, Hua Q, Ma Z, Lv Q. Downregulation of miR-200a-3p induced by hepatitis B virus X (HBx) protein promotes cell proliferation and invasion in HBV-infection-associated hepatocarcinoma. Pathol Res Pract. 2017;213:1464–9.
Bai PS, Xia N, Sun H, Kong Y. Pleiotrophin, a target of miR-384, promotes proliferation, metastasis and lipogenesis in HBV-related hepatocellular carcinoma. J Cell Mol Med. 2017;21:3023–43.
Wang Y, Wang CM, Jiang ZZ, Yu XJ, Fan CG, Xu FF, Zhang Q, et al. MicroRNA-34c targets TGFB-induced factor homeobox 2, represses cell proliferation and induces apoptosis in hepatitis B virus-related hepatocellular carcinoma. Oncol Lett. 2015;10:3095–102.
Xiang Z, Wang S, Xiang Y. Up-regulated microRNA499a by hepatitis B virus induced hepatocellular carcinogenesis via targeting MAPK6. PLoS ONE. 2014;9:e111410.
Qin X, Chen J, Wu L, Liu Z. MiR-30b-5p acts as a tumor suppressor, repressing cell proliferation and cell cycle in human hepatocellular carcinoma. Biomed Pharmacother. 2017;89:742–50.
Lamontagne J, Steel LF, Bouchard MJ. Hepatitis B virus and microRNAs: complex interactions affecting hepatitis B virus replication and hepatitis B virus-associated diseases. World J Gastroenterol. 2015;21:7375–99.
Deng M, Zeng C, Lu X, He X, Zhang R, Qiu Q, Zheng G, et al. miR-218 suppresses gastric cancer cell cycle progression through the CDK6/Cyclin D1/E2F1 axis in a feedback loop. Cancer Lett. 2017;403:175–85.
Lal A, Navarro F, Maher CA, Maliszewski LE, Yan N, O’Day E, Chowdhury D, et al. miR-24 Inhibits cell proliferation by targeting E2F2, MYC, and other cell-cycle genes via binding to “seedless” 3ʹUTR microRNA recognition elements. Mol Cell. 2009;35:610–25.
Zhu X, Wu L, Yao J, Jiang H, Wang Q, Yang Z, Wu F. MicroRNA let-7c inhibits cell proliferation and induces cell cycle arrest by targeting CDC25A in human hepatocellular carcinoma. PLoS ONE. 2015;10:e0124266.
This study is supported by grants from the National Natural Science Foundations of China (Grant Nos. U1803124, 11705283, 81530085), and Natural Science Foundation of Hunan Province (Grant No. 2019JJ40396) .
Ethics approval and consent to participate
This study was performed in accordance with the Declaration of Helsinki and was approved by the First Affiliated Hospital of the University of South China (IRB approval no. NHUH-2017-38). All enrolled patients signed informed consent forms.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 2: Figure S1. Rescue analysis of the effects of miR-0308-3p-resistance CDK6 and Cyclin D1on G1/S cell cycle. (A) The effects of miR-0308-3p-resistance CDK6 on the distribution of G1 and S phases in cells using flow cytometry. A representative flow cytometry histogram of the cell-cycle progression of HepG2 cells is shown with and without the miR-0308 resistance CDK6. (B) Quantitative measurement of the G1 and S phases in HepG2 cells with or without the miR-0308 resistance CDK6. (C) The effects of miR-0308-3p-resistance Cyclin D1 on the distribution of G1 and S phases in cells using flow cytometry. A representative flow cytometry histogram of the cell-cycle progression of HepG2 cells is shown with and without the miR-0308 resistance Cyclin D1. (D) Quantitative measurement of the G1 and S phases in HepG2 cells with or without the miR-0308 resistance Cyclin D1.
About this article
Cite this article
Dai, X., Huang, R., Hu, S. et al. A novel miR-0308-3p revealed by miRNA-seq of HBV-positive hepatocellular carcinoma suppresses cell proliferation and promotes G1/S arrest by targeting double CDK6/Cyclin D1 genes. Cell Biosci 10, 24 (2020). https://doi.org/10.1186/s13578-020-00382-7