Identification of differentially expressed long non-coding RNAs in human ovarian cancer cells with different metastatic potentials

Shi-Ping Liu, Jia-Xin Yang, Dong-Yan Cao, Keng Shen

Department of Obstetrics and Gynecology, Peking Union Medical College Hospital, Peking Union Medical College, Chinese Academy of Medical Sciences, Beijing 100730, China

Introduction

Epithelial ovarian cancer is the leading cause of death from gynecologic cancer in developed regions worldwide and second to cervical cancer in developing countries1,2.Although new therapeutic drugs have been developed to treat ovarian cancer patients, the mortality rate among women diagnosed with metastatic disease remains high, with only 27% of them surviving more than 5 years post-diagnosis3.Therefore, the identification of genes and regulatory mechanisms involved in metastasis that could enhance the understanding of cancer progression and result in the development of new therapeutics has attracted much interest.

Long non-coding RNAs (lncRNAs) are non-coding RNAs that are longer than 200 nucleotides.Recent studies have demonstrated that lncRNAs play important roles in carcinogenesis and cancer metastasis.Aberrant lncRNA pro files have been noted in multiple cancers4-6, but the lncRNAs involved in ovarian cancer metastasis have not been identified.In this study, we presented the lncRNA expression profiles in paired ovarian cancer cell lines with different metastatic potentials and selected nine differentially expressed lncRNAs for validation using RT-qPCR.

Materials and methods

Cell lines and cell culture

The human ovarian cancer cell lines SKOV3 and SKOV3.ip1 were obtained from the American Type Culture Collection and M.D.Anderson Cancer Center, respectively.The SKOV3.ip1 cell line was derived from ascites from a nude mouse injected i.p.with SKOV3 cells and demonstrated a higher metastatic potential compared with the parental cell line7.Both cell lines were maintained in Dulbecco’s modified Eagle’s medium supplemented with 10% fetal bovine serum and cultured in a humidi fied incubator (5% CO2) at 37 ℃.

Validation of invasive capability of ovarian cancer cells

In vitro invasion assay was performed using 24-well transwell units with polycarbonate filters (pore size: 8 µm) coated on the upper side with reconstituted basement membrane matrix (BD Biosciences,USA).The cells were harvested, and 4×104cells in 100 µL of serum-free medium were placed in the upper part of the transwell unit and were allowed to invade the membrane for 72 h at 37 ℃.Successfully penetrating cells were fixed, stained, and quantified at optical density of 570 nm after extraction.The results are reported as the averages of three individual experiments containing three replicates per condition.

RNA preparation

Total RNA was extracted from the cell samples using the TRIzol reagent (Invitrogen) according to the manufacturer’s instructions.The RNA integrity was analyzed using the Agilent 2100 Bioanalyzer (Agilent Technologies).

Microarray

Microarray analysis was performed by a commercial company(Oebiotech, PRC), using SurePrint G3 Human Gene Expression 8×60K v2 (Agilent Technologies) that is designed to contain approximately 12,000 lncRNAs.Briefly, samples were used to synthesize cDNA, and labeled cRNA was then synthesized and hybridized to the microarray.After hybridization and washing,processed slides were scanned with the Agilent Microarray Scanner (Agilent Technologies), and the acquired array images were analyzed using Agilent Feature Extraction Software (Agilent Technologies), which performs background subtractions.Quantile normalization and subsequent data processing were performed using the GeneSpringGX v.11.0 software package(Agilent Technologies).A threshold of fold change ＞2 was used to screen upregulated or downregulated lncRNAs.

Real-time quantitative PCR validation

To validate the microarray data, we selected nine differentially expressed lncRNAs (MALAT1, H19, XIST, UCA1,CCAT1, LOC645249, LOC100128881, LOC728228, and LOC100292680).MALAT1, H19, XIST, UCA1, and CCAT1 were selected because they had been associated with cancer.The four remaining lncRNAs were selected randomly from all aberrantly expressed lncRNAs.

Quantification was performed with a two-step reaction process, i.e., through reverse transcription (RT) and PCR.Each RT reaction consisted of 0.5 µg RNA, 2 µL of PrimerScript Buffer, 0.5 µL of oligo dT, 0.5 µL of random 6 mers, and 0.5 µL of PrimerScript RT Enzyme Mix I (TaKaRa, Japan), in a total volume of 10 µL.Reactions were performed in a GeneAmp® PCR System 9700 (Applied Biosystems, USA) for 15 min at 37 ℃,followed by heat inactivation of RT for 5 s at 85 ℃.The 10 µL RT reaction mix was then diluted ×10 in nuclease-free water and held at –20 ℃.

Real-time PCR was performed using LightCycler® 480 II Real-time PCR Instrument (Roche, Swiss) with 10 µL PCR reaction mixture that included 1 µL of cDNA, 5 µL of 2×LightCycler® 480 SYBR Green I Master (Roche, Swiss), 0.2 µL of forward primer, 0.2 µL of reverse primer, and 3.6 µL of nucleasefree water.Reactions were incubated in a 384-well optical plate(Roche, Swiss) at 95 ℃ for 10 min, followed by 40 cycles at 95 ℃ for 10 s, 60 ℃ for 30 s.At the end of the PCR cycles,melting curve analysis was performed to validate the specific generation of the expected PCR product.Three independent experiments were performed with each sample run in triplicate.The primers were synthesized by Generay Biotech (Generay,PRC) and are listed in Table 1.The expression levels of lncRNAs were normalized with respect to GAPDH and were calculated using the 2–ΔΔCtmethod.

Table 1 Primers used in RT-qPCR

Statistical analysis

Student’s t-test was performed to evaluate any significant difference.A value of P＜0.05 was considered statistically significant.

Results

Increased invasive capability of SKOV3.ip1 cells

The abilities of the two cell lines to penetrate a uniform layer of extracellular matrix-resembling material were tested using Transwell invasive assay to compare their invasive capabilities.Cells were seeded on top of a porous membrane.Fetal bovine serum was used as a chemoattractant in the lower compartment.The results showed a marked increase in the invasive ability of the SKOV3.ip1 cell line compared with that of the SKOV3 cell line (Figure 1).

Differentially expressed metastasis-related lncRNAs

A total of 4,956 lncRNAs were detected in the microarray.Data analysis showed that, compared with the parental SKOV3 cells,more than a 2-fold deregulation was found in 1,161 lncRNAs in SKOV3.ip1 cells, of which 583 were lncRNAs were upregulated and 578 downregulated.The expression levels of nine lncRNAs were validated using RT-qPCR, and seven confirmed the deregulation found by microarray analysis (Table 2).

Discussion

Cancer metastasis is the leading cause of death in cancer patients8.In the current study, we primarily identified the lncRNAs involved in the metastasis of epithelial ovarian cancer.We analyzed paired high and low metastatic ovarian cancer cells by microarray assay and selected nine lncRNAs for validation by RT-qPCR.Interestingly, we found that 1,161 of the 4,956 detected lncRNAs were differentially expressed between the two cell lines.Among these 1,161 lncRNAs, 583 and 578 were found to be upregulated and downregulated, respectively.Most of these lncRNAs have not been functionally characterized.These lncRNAs may be identified as novel biomarkers and therapeutic targets.

Figure 1 SKOV3.ip1 cells are more invasive than SKOV3 cells(＊P＜0.01).

Table 2 Deregulated genes in SKOV3.ip1 compared to SKOV3 cells

Nine lncRNAs were selected for validation on the basis of relevant literature and the microarray data.The qPCR results of seven lncRNAs (MALAT1, H19, UCA1, CCAT1, LOC645249,LOC100128881, and LOC100292680) were consistent with the deregulation found by microarray analysis, reflecting the reliability of the microarray data to some extent.Of the seven lncRNAs, H19 was markedly downregulated in the SKOV3.ip1 cells.Similarly, Zhang et al.9reported that H19 suppressed hepatocellular carcinoma metastasis and the expression of markers of epithelial-to-mesenchymal transition by altering the miR-200 pathway.However, Yang et al.10reported that H19 upregulation increased cell proliferation in gastric cancer cells.Furthermore, another study has found that H19 was remarkably upregulated in bladder cancer tissues, and upregulated H19 promoted cancer cell migration both in vitro and in vivo by associating with EZH2 and inhibiting E-cadherin expression11.Another lncRNA, MALAT1, is known to promote cancer metastasis in several cancer types, including lung, bladder, and cervical cancers12-14, and its overexpression has been reported to predict tumor recurrence of hepatocellular carcinoma after liver transplantation15.However, in the current study, it was downregulated in the cell culture with high metastatic potential cells.Thus, the function of one lncRNA may vary with different cancer types and context.H19 and MALAT1 may play a suppressive role in ovarian cancer metastasis.Published data on UCA1 and CCAT1 are limited.One study has suggested that UCA1 might contribute to bladder cancer progression16,and CCAT1 has been reported to be overexpressed in colon cancer17,18.The remaining three lncRNAs are uncharacterized,and their biological functions remain unclear.

To our knowledge, this study is the first to describe the lncRNA profiles in paired high and low metastatic ovarian cancer cells.A number of differentially expressed lncRNAs were identified, and these lncRNAS may play a partial or key role in the metastasis of epithelial ovarian cancer.Further studies are necessary to determine the roles of these lncRNAs in ovarian cancer metastasis.

Conflict of interest statement

No potential conflicts of interest are disclosed.

1.Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D.Global cancer statistics.CA Cancer J Clin 2011;61:69-90.

2.Siegel R, Naishadham D, Jemal A.Cancer statistics, 2013.CA Cancer J Clin 2013;63:11-30.

3.SEER cancer statistics review; 1975-2008.[http://seer.cancer.gov/csr/1975_2008/, based on November 2010 SEER data submission, posted to the SEER web site, 2011.]

4.Yu G, Yao W, Wang J, Ma X, Xiao W, Li H, et al.LncRNAs expression signatures of renal clear cell carcinoma revealed by microarray.PLoS One 2012;7:e42377.

5.Yang F, Zhang L, Huo XS, Yuan JH, Xu D, Yuan SX, et al.Long noncoding RNA high expression in hepatocellular carcinoma facilitates tumor growth through enhancer of zeste homolog 2 in humans.Hepatology 2011;54:1679-1689.

6.Cao WJ, Wu HL, He BS, Zhang YS, Zhang ZY.Analysis of long non-coding RNA expression pro files in gastric cancer.World J Gastroenterol 2013;19:3658-3664.

7.Yu D, Wolf JK, Scanlon M, Price JE, Hung MC.Enhanced c-erbB-2/neu expression in human ovarian cancer cells correlates with more severe malignancy that can be suppressed by E1A.Cancer Res 1993;53:891-898.

8.Gupta GP, Massagué J.Cancer metastasis: building a framework.Cell 2006;127:679-695.

9.Zhang L, Yang F, Yuan JH, Yuan SX, Zhou WP, Huo XS, et al.Epigenetic activation of the MiR-200 family contributes to H19-mediated metastasis suppression in hepatocellular carcinoma.Carcinogenesis 2013;34:577-586.

10.Yang F, Bi J, Xue X, Zheng L, Zhi K, Hua J, et al.Up-regulated long non-coding RNA H19 contributes to proliferation of gastric cancer cells.FEBS J 2012;279:3159-3165.

11.Luo M, Li Z, Wang W, Zeng Y, Liu Z, Qiu J.Long non-coding RNA H19 increases bladder cancer metastasis by associating with EZH2 and inhibiting E-cadherin expression.Cancer Lett2013;333:213-221.

12.Gutschner T, Hämmerle M, Eissmann M, Hsu J, Kim Y, Hung G, et al.The noncoding RNA MALAT1 is a critical regulator of the metastasis phenotype of lung cancer cells.Cancer Res 2013;73:1180-1189.

13.Han Y, Liu Y, Nie L, Gui Y, Cai Z.Inducing cell proliferation inhibition, apoptosis, and motility reduction by silencing long noncoding ribonucleic acid metastasis-associated lung adenocarcinoma transcript 1 in urothelial carcinoma of the bladder.Urology 2013;81:209.e1-7.

14.Guo F, Li Y, Liu Y, Wang J, Li Y, Li G.Inhibition of metastasisassociated lung adenocarcinoma transcript 1 in CaSki human cervical cancer cells suppresses cell proliferation and invasion.Acta Biochim Biophys Sin (Shanghai) 2010;42:224-229.

15.Lai MC, Yang Z, Zhou L, Zhu QQ, Xie HY, Zhang F, et al.Long non-coding RNA MALAT-1 overexpression predicts tumor recurrence of hepatocellular carcinoma after liver transplantation.Med Oncol 2012;29:1810-1816.

16.Wang F, Li X, Xie X, Zhao L, Chen W.UCA1, a non-proteincoding RNA up-regulated in bladder carcinoma and embryo,influencing cell growth and promoting invasion.FEBS Lett2008;582:1919-1927.

17.Nissan A, Stojadinovic A, Mitrani-Rosenbaum S, Halle D,Grinbaum R, Roistacher M, et al.Colon cancer associated transcript-1: a novel RNA expressed in malignant and premalignant human tissues.Int J Cancer 2012;130:1598-1606.

18.Alaiyan B, Ilyayev N, Stojadinovic A, Izadjoo M, Roistacher M,Pavlov V, et al.Differential expression of colon cancer associated transcript1 (CCAT1) along the colonic adenoma-carcinoma sequence.BMC Cancer 2013;13:196.[Epub ahead of print].