Long noncoding RNA ZNFX1-AS1 suppresses growth of hepatocellular carcinoma cells by regulating the methylation of miR-9.

Many long noncoding RNAs have been reported to play pivotal roles in cancer biology. Among them, the long noncoding RNA ZNFX1-AS1 has been confirmed to function in breast cancer progression, but the role of ZNFX1-AS1 in hepatocellular carcinoma (HCC) growth and the related molecular mechanisms still remains unknown. In the present study, we first identified the expression of ZNFX1-AS1 in HCC patients' specimens and HCC cell lines through quantitative reverse transcription polymerase chain reaction. Next, the effects of ZNFX1-AS1 on HCC cell growth and apoptosis were analyzed. MTT assay was used to measure the cell numbers, and fluorescence-activated cell sorting analysis was performed to evaluate cell apoptosis. Finally, the relationship between ZNFX1-AS1 and miR-9 in HCC was studied. Our results suggest that ZNFX1-AS1 was markedly downregulated in HCC samples and cell lines. Overexpression of ZNFX1-AS1 inhibited the cell proliferation and colony formation in HCC cell lines and also induced HCC cell apoptosis. Additionally, miR-9 was lowly expressed in HCC tissues and positively correlated with ZNFX1-AS1 expression. Meanwhile, significant upregulation of miR-9 and downregulation of the methylation of miR-9 promoter CpG island were observed when ZNFX1-AS1 was overexpressed. In summary, our results indicate that ZNFX1-AS1 plays a vital role in HCC progression via regulating the methylation of miR-9 and may be a potential tumor suppressor.

Introduction

Hepatocellular carcinoma (HCC) has become the most common leading cause of cancer-related mortality globally.1 Despite increasing exploration for novel treatments, such as operation, liver transplantation, and chemoradiotherapy, there has still been little success in improving the 5-year overall survival (OS) rate of HCC patients.2,3 The poor understanding of the molecular mechanisms of development and progression of HCC is one of the most important reasons that result in high fatality and poor prognosis.4,5 Therefore, probing the molecular mechanisms underlying occurrence of HCC is of great clinical value for developing new preventable and therapeutic strategies.

Large-scale researches on long noncoding RNAs (lncRNAs), non-protein-coding transcripts that are more than 200 nt in length, have shown that lncRNAs are involved in a variety of biological processes, including silencing the X chromosome, transcriptional regulation, modifying chromatin, splicing, and so on.6–12 Recently, numerous specifically dysregulated lncRNAs have been found in several cancer types,13 which can be seen as a potential clinical advantage as these lncRNAs can function as diagnostic or prognostic biomarkers.14,15 LncRNAs, such as anti-NOS2A16 and MEG3,17–19 have been identified to be associated with the pathogenesis of brain tumor. In addition, H19,20 BC200,21,22 HOTAIR,23 and MALAT124 have been proven to play crucial roles in the pathological process of breast tumor. Since the lncRNA HULC was reported to be specifically upregulated in HCC,25 an increasing number of lncRNAs have been found to function in HCC progression, such as MALAT1, PCNA-AS1, and HOTTIP.26–28 ZNFX1-AS1, an lncRNA that transcripts from the antisense strand near the 5′-end of the protein-coding gene Znfx1, hosts three C/D box snoRNAs (SNORDs): Snord12, Snord12b, and Snord12c.29 ZNFX1-AS1 has been reported to be related to several cancers, including breast, gastric, and colorectal cancer.30 In addition, ZNFX1-AS1 is significantly downregulated in breast cancer, and hence can function as a potential biomarker for breast cancer.29 Nevertheless, its role in other cancers remains largely unknown, especially in HCC.

In this study, we identified the significant downregulation of ZNFX1-AS1 in HCC tissues, which was associated with the proliferation and apoptosis of HCC cells. Moreover, ZNFX1-AS1 was also found to function in the progression of HCC via regulating the methylation of miR-9.

Materials and methods

HCC specimens and cell lines

Eighty-four pairs of samples (including 84 HCC samples and normal adjacent liver tissues) from HCC patients were obtained from Eastern Hepatobiliary Surgery Hospital (Shanghai, People’s Republic of China). The study was approved by the Ethics Committee of the Eastern Hepatobiliary Surgery Hospital. Written informed consent was obtained for the use of tissue samples from all patients. None of the patients recruited for this study had anticancer treatments before the surgery. Liver cancer cell lines HepG2, Hep3B, SNU449, HuH-7, and normal human hepatocytes (HH) were purchased from American Type Culture Collection (ATCC, Manassas, VA, USA). The earlier cell lines were cultured in high-glucose Dulbecco’s Modified Eagle’s Medium (Thermo Fisher Scientific, Waltham, MA, USA), supplemented with 10% fetal bovine serum, 100 U/mL penicillin, and 100 mg/mL streptomycin in a humidified incubator with 5% CO2 at 37°C.

Quantitative reverse transcription polymerase chain reaction and methylation-specific polymerase chain reaction

Total RNA from HCC tissues was extracted using Trizol reagent (Thermo Fisher Scientific). Real-time polymerase chain reaction (PCR) was performed using the SYBR Green Real-Time PCR Master Mix (TaKaRa, Tokyo, Japan) as described.28 To examine the expression of miR-9, we performed quantitative PCR using synthetic micro-RNA sequences in a calibration curve. The following primers were used:

Forward: 5′-GGGTCTTTGGTTATCTAGC-3′; Reverse: 5′-TGCGTGTCGTGGAGTC3-3′.

Methylation-specific polymerase chain reaction (MSP) was performed as previously described.31 The expression of ZNFX1-AS1 was normalized with GAPDH. The primers of ZNFX1-AS1 and GAPDH were bought from Thermo Fisher Scientific. The primers used for miR-9 MSP in this study are listed as follows:

miR-9-msp-U forward: 5′-GTGAGATTTGTTTTATTTTGTTTGT-3′; Reverse: 5′-AAAACCCTAAATTCCAAATACCAA-3′; miR-9-msp-M forward: 5′-AGGTGAGATTCGTTTTATTTTGTTC-3′; Reverse: 5′-AAAACCCTAAATTCCAAATACCG-3′.

Vector construction and transfection

Knockdown lentivirus and the relative negative control lentivirus were purchased from GenePharma (Shanghai, People’s Republic of China). SNU449 and HepG2 cells were transfected with siRNA lentivirus according to the manufacturer’s instructions. Stable transfection lentivirus cell lines were used for subsequent assays. The cDNA encoding ZNFX1-AS1 was PCR-amplified by the PfuUltra II Fusion HS DNA Polymerase (Stratagene, Agilent Technologies, Santa Clara, CA, USA) and subcloned into a pcDNA vector (Invitrogen), named pcDNA-ZNFX1-AS1.

Cell proliferation and apoptosis assay

The cells were seeded in 96-well plates for cell proliferation assay using a 3-(4, 5-dimethylthiazol-2-yl)-2,5-diphenyltet-razolium bromide (MTT) Cell Proliferation/Viability Assay Kit (Sigma-Aldrich Chemie Gmbh, Munich, Germany) in accordance with the guidelines. The cell apoptosis was measured using an Annexin-V-Fluos and Propidium Iodide (PI) Apoptosis Detection Kit (Sigma) by a flow cytometer (Becton Dickinson, San Jose, CA, USA) according to the manufacturer’s instructions.

Soft sugar colony formation assay

The colony number of HepG2 and SNU449 cells was recorded after 2 weeks of soft sugar colony formation analysis. Briefly, cells were suspended in 1.5 mL complete medium supplemented with 0.45% agarose (Thermo Fisher Scientific), and then placed in 35 mm tissue culture plates containing 1.5 mL complete medium and agarose (0.75%) on the bottom layer. Finally, cell colonies were stained with 0.005% crystal violet and analyzed using a microscope.

Western blotting

Western blotting analysis was performed as described in the study by Zhang et al.32 Thirty-microgram total protein was extracted and separated by 10% sodium dodecyl sulfate poly-acrylamide gel electrophoresis (SDS-PAGE) gel, and then transferred to polyvinylidene difluoride membranes (Mil-lipore, Darmstadt, Germany). The membrane was blocked with specific primary antibodies and then incubated with secondary antibodies labeled with horseradish peroxidase and detected by chemiluminescence (Thermo, Rockford, IL, USA). DNMT1 and GAPDH antibodies were obtained from Santa Cruz Biotechnology, Santa Cruz, CA, USA and GAPDH was used as protein loading control.

Statistical analysis

Data were presented as mean ± standard deviation (SD) and analyzed using SPSS 17.0 software (SPSS Inc., Chicago, IL, USA). Statistical differences between two independent groups were determined using unpaired t-test. For multiple comparison, one-way analysis of variance was used. Linear regression analysis was performed to analyze the correlation between expressions of miR-9 and ZNFX1-AS1. The Kaplan–Meier method was used to estimate the OS rate. Survival differences according to hMOF expression were analyzed using the log-rank test. P<0.05 was considered statistically significant.

Results

Expression of ZNFX1-AS1 is significantly downregulated in HCC samples and cell lines

We measured ZNFX1-AS1 expression in HCC tissues and cell lines through quantitative PCR. As shown in Figure 1A, compared to normal liver tissues, the expression of ZNFX1-AS1 was prominently decreased in HCC tissues. Remarkably, consistent with the data derived from HCC samples, the expression of ZNFX1-AS1 was significantly decreased in all cell lines (HepG2, Hep3B, SNU449, HuH7) except in normal liver cell line HH (Figure 1B). In summary, ZNFX1-AS1 was downregulated in HCC tissues and cell lines. We selected HepG2 and SNU449 (the relatively higher expression of ZNFX1-AS1) as representatives for further research.

Figure 1

ZNFX1-AS1 is downregulated in human HCC tissue and cell lines.

Notes: (A) RNA was extracted from human HCC tissue and adjacent normal tissues, and ZNFX1-AS1 expression was analyzed by qRT-PCR. **P<0.01 by Student’s t-test. (B) ZNFX1-AS1 expression of HCC cells and HH was assessed by qRT-PCR. **P<0.01 vs HH by one-way ANOVA.

Abbreviations: HCC, hepatocellular carcinoma; qRT-PCR, quantitative reverse transcription polymerase chain reaction; HH, normal hepatocytes; ANOVA, analysis of variance.

Low expression of ZNFX1-AS1 is associated with poor prognosis of HCC

To analyze the relationship between ZNFX1-AS1 expression and HCC patients’ prognosis, we detected rates of OS and disease-free survival of patients using Kaplan–Maier analysis. As shown in Figure 2A, the 5-year OS rate of patients with high expression of ZNFX1-AS1 was significantly higher than that for the low ZNFX1-AS1 expression. In consistence, the 5-year disease-free survival rate in high ZNFX1-AS1 expression group was also higher than that in the low expression group (Figure 2B).

Figure 2

Low-level expression of ZNFX1-AS1 is correlated with poor prognosis.

Notes: Kaplan–Meier curve comparing the survival times of HCC cells with the low (<25th percentile) and high (>25th percentile) ZNFX1-AS1 expression was determined using the q-PCR results. (A) Overall survival and (B) disease-free survival; n=93 in the low ZNFX1-AS1 expression group and n=61 in the high ZNFX1-AS1 expression group.

Abbreviations: HCC, hepatocellular carcinoma; q-PCR, quantitative polymerase chain reaction.

ZNFX1-AS1 suppresses cell proliferation and promotes cell apoptosis in HCC cells

To identify the effects of ZNFX1-AS1 on proliferation and apoptosis of HCC cells, we first determined whether ZNFX1-AS1 siRNA inhibits the expression of ZNFX1-AS1 in HepG2 and SNU449 cells. The results revealed that cell proliferation and colony formation of HepG2 and SNU449 cells were increased when treated with si-ZNFX1-AS1 (Figures 3A–D; S1A). Next, ZNFX1-AS1 was overexpressed to further identify the relationship between ZNFX1-AS1 and cell proliferation and colony formation of HepG2 and SNU449 cells. Our results suggested that ZNFX1-AS1 overexpression inhibited cell proliferation and colony formation of HepG2 and SNU449 cells (Figures 3E–H; S1B). Furthermore, we investigated the effect of ZNFX1-AS1 on cancer cell apoptosis by performing fluorescence-activated cell sorting assay. We found that ZNFX1-AS1 overexpression induced cancer cell apoptosis (Figures 3I and J; S2A and B). Collectively, these results suggest that ZNFX1-AS1 inhibits HCC cell proliferation and promotes cell apoptosis.

Figure 3

ZNFX1-AS1 suppresses cancer cell growth in vitro.

Notes: (A–D) ZNFX1-AS1 silencing accelerates HCC cell growth and colony formation. HepG2 and SNU449 cells were infected with control siRNA or si-ZNFX1-AS1. (A and B) The relative cell numbers were measured at certain time points using MTT method. *P<0.05, **P<0.01 vs si-NC. (C and D) Soft sugar colony forming analysis was performed, and the colony numbers were recorded. (E–H) ZNFX1-AS1 overexpression inhibits HCC cell growth and colony formation. HepG2 and SNU449 cells were transfected with 1 μg pCDNA-ZNFX1-AS1 or empty vector. (E and F) The relative cell numbers were measured at certain time points using MTT method. **P<0.01 vs empty vector. (G and H) Soft sugar colony forming analysis was performed, and the colony numbers were recorded. (I and J) ZNFX1-AS1 overexpression induces HCC cell apoptosis. HepG2 and SNU449 cells were transfected with pCDNA-ZNFX1-AS1 or empty vector, and fluorescence-activated cell sorting analysis was performed on the cells with indicated markers.

Abbreviations: HCC, hepatocellular carcinoma; si-NC, si-negative control.

ZNFX1-AS1 induces miR-9 expression through inhibiting methylation of CpG island at its promoter

Furthermore, we aimed to explore the mechanisms underlying the potent effect of ZNFX1-AS1 on HCC. Since miR-9 has been reported to be involved in HCC progression,33 we hypothesized that ZNFX1-AS1 may regulate HCC growth by miR-9. Firstly, we examined the expression level of miR-9 in HCC tissues. Compared to the normal tissues, the expression of miR-9 was markedly downregulated in HCC tissues (Figure 4A). Moreover, ZNFX1-AS1 expression level was positively related to miR-9 expression level in HCC tissues (Figure 4B), which implied that ZNFX1-AS1 promoted miR-9 expression in HCC. Next, we explored the relationship between ZNFX1-AS1 and miR-9. As shown in Figure 4C and D, miR-9 was significantly upregulated in HepG2 and SNU449 cells when ZNFX1-AS1 was overexpressed. To further investigate whether ZNFX1-AS1 regulates HCC growth through miR-9, we examined the methylation of miR-9 promoter in HepG2 and SNU449 cells. DNMT1, a kind of methyltransferase, can regulate miRNA expression by controlling methylation of CpG island. Additionally, it was found that ZNFX1-AS1 overexpression inhibited DNMT1 expression (Figures 4E and F; S3). At the same time, the methylation inhibitor Aza increased the expression of miR-9 in HepG2 and SUN449 cells (Figure 4G and H). Indeed, we found that DNMT1 overexpression inhibited the expression of miR-9 in HepG2 cells (Figure S4). Finally, MSP analysis was used to analyze the methylation of miR-9 CpG sites in HepG2, SNU449, Hep3B, and Huh7 cells. Obviously, the overexpression of ZNFX1-AS1 inhibited the methylation of miR-9 CpG island (Figure 4I). Taken together, these results demonstrated that ZNFX1-AS1 may function in HCC progression by regulating the methylation of miR-9 CpG island.

Figure 4

ZNFX1-AS1 induces miR-9 expression through inhibiting methylation of CpG island in its promoter.

Notes: (A) miR-9 expression level is downregulated in HCC tissues. RNA was extracted from HCC tissue and adjacent normal tissues, and miR-9 expression level was analyzed. **P<0.01 vs normal. (B) ZNFX1-AS1 expression level is positively related to miR-9 level in HCC tissues. r2=0.5092; P<0.01. (C and D) The effect of ZNFX1-AS1 on miR-9 was analyzed by qRT-PCR in HepG2 (C) and SNU449 (D) cells. ZNFX1-AS1 transfection efficiency was measured by qRT-PCR. *P<0.05, **P<0.01 vs empty vector. (E and F) ZNFX1-AS1 expression inhibits the expression of DNMT1. HepG2 and SNU449 cells were transfected with 1 μg pCDNA-ZNFX1-AS1 or empty vector as in (C and D), and the DNMT1 expression level was measured by qRT-PCR. ZNFX1-AS1 transfection efficiency was measured by qRT-PCR. **P<0.01 vs empty vector. (G and H) The effect of Aza on miR-9 expression was examined by qRT-PCR in HepG2 (G) and SNU449 (H) cells. **P<0.01 vs control. (I) The methylation of miR-9 promoter CpG sites was analyzed by MSP analysis in HepG2, SNU449, Hep3B, and Huh7 cells. Bands in “U” and “M” lanes represent products obtained with methylation-specific or un-methylation-specific primers, respectively.

Abbreviations: HCC, hepatocellular carcinoma; qRT-PCR, quantitative reverse transcription polymerase chain reaction; MSP, methylation-specific PCR.

Discussion

A better understanding of the mechanisms underlying HCC metastasis is of importance. Presently, numerous lncRNAs are found to play crucial roles in the metastasis and invasion of HCC cells.34,35 Genetic variants at 20q13 are associated with some tumors.36,37 LncRNA ZNFX1-AS1 is located in the 20q13.13 region, which functions in breast cancer progression, implying a potential role of lncRNA ZNFX1-AS1 in HCC metastasis.30 However, the biological function or molecular mechanisms of ZNFX1-AS1 in HCC remain unclear. Therefore, understanding the effects of ZNFX1-AS1 on HCC is essential.

Here, we identified a novel ZNFX1-AS1-dependent molecular mechanism during the progression of HCC. We reported that ZNFX1-AS1 was prominently downregulated in HCC tissues and liver cancer cell lines. ZNFX1-AS1 low expression was associated with poor prognosis. To explore the functional role of ZNFX1-AS1 in HCC progression, we adopted lentivirus to knockdown or overexpress ZNFX1-AS1 in SNU449 and HepG2 cells. ZNFX1-AS1 was confirmed to function in the proliferation, colony formation, and apoptosis of HCC cells. Our results revealed that ZNFX1-AS1 may play a pivotal role in the growth of HCC.

As miR-9 has been identified to play a role in HCC growth,32 we tried to seek the correlations between ZNFX1-AS1 and miR-9 in HCC progression. In this article, we identified that the expression of miR-9 was downregulated in HCC tissues and positively correlated with the expression of ZNFX1-AS1. More importantly, miR-9 was markedly upregulated in HepG2 and SNU449 cells with ZNFX1-AS1 overexpression, implying that miR-9 may be a target of ZNFX1-AS1 in HCC cells. To further confirm whether ZNFX1-AS1 regulated HCC growth through miR-9, we examined the methylation of miR-9 in HepG2 and SNU449 cells, respectively. Obviously, the overexpression of ZNFX1-AS1 inhibited the methylation of miR-9 CpG island. Conclusively, these results indicated that ZNFX1-AS1 may function in HCC progression by regulating the methylation of miR-9 CpG island. As shown, ZNFX1-AS1 inhibits the expression of DNMT1, which reduces the expression of miR-9. However, one limitation of this mechanism is that DNMT1 is a pan-DNA methyl-transferase; it does not specifically regulate the expression of miR-9. Therefore, another potential target of ZNFX1 is needed to explore the functional role of ZNFX1-AS1 in HCC progression in a future study and we should further determine whether miR-9 contributes to this function.

Conclusion

Taken all together, we demonstrated that lncRNA ZNFX1-AS1 is significantly downregulated in HCC samples and cell lines, and low expression of ZNFX1-AS1 is associated with poor prognosis for HCC patients. Furthermore, we also showed that ZNFX1-AS1 suppresses HCC progression by regulating the methylation of miR-9 CpG island. LncRNA ZNFX1-AS1 may serve as a potential tumor suppressor gene.

Representative photographs show that ZNFX1-AS1 knockdown promotes colony formation of HepG2 cells (A) or ZNFX1-AS1 overexpression inhibits colony formation of HepG2 cells (B).

Abbreviation: si-NC, si-negative control.

Representative fluorescence-activated cell sorting analysis of Annexin-V and propidium iodide (PI) staining of HepG2 cells (A) or SUN-449 (B) with/without ZNFX1-AS1 overexpression.

Note: Q1: necrotic cell, Q2: late apoptotic cells, Q3: normal cells, Q4: early apoptotic cells.

Abbreviation: Q, quadrant.

ZNFX1-AS1 overexpression reduces the protein level of DNMT1 in HepG2 cells.

Abbreviation: DNMT1, DNA (cytosine-5)-methyltransferase 1.

DNMT1 overexpression inhibits miR-9 expression in HepG2 cells.

Abbreviation: Ctrl, control.