Identification of lncRNA FAM83H-AS1 as a novel prognostic marker in luminal subtype breast cancer

OncoTargets and Therapy Identification of lncRNA FAM83H-AS1 as a novel prognostic marker in luminal subtype breast cancer 0 Department of Surgical Oncology, The First Affiliated Hospital of Wenzhou Medical University , Wenzhou, Zhejiang , People's Republic of China PowerdbyTCPDF(ww.tcpdf.org) Background: Luminal subtype breast cancer accounts for a predominant number of breast cancers. Considering the heterogeneity of the disease, it is urgent to develop novel biomarkers to improve risk stratification and optimize therapy choices. Long non-coding RNA (lncRNA) represents an emerging and understudied class of transcripts that play a significant role in cancer biology. Growing knowledge of cancer-associated lncRNAs contributes to the development of molecular markers for prognosis evaluation and gene therapy. Materials and methods: Three pairs of primary luminal subtype breast cancer tissues and adjacent non-cancerous tissues were collected and sequenced. EBseq algorithm was used to identify differentially expressed lncRNAs. RNA sequencing data from The Cancer Genome Atlas (TCGA) database were used to validate the robustness of our RNA-seq results. Kaplan-Meier and Cox regression analyses were utilized to assess the association between the lncRNAs and overall survival of patients in TCGA cohort. Results: A total of 796 lncRNAs were significantly dysregulated in luminal subtype breast cancer, including 436 upregulated and 360 downregulated lncRNAs. Among them, FAM83H antisense RNA 1 (FAM83H-AS1) was the most upregulated lncRNA, whereas GSN antisense RNA 1 (GSN-AS1) was the most downregulated lncRNA. Moreover, we proved that the high expression level of FAM83H-AS1 indicated unfavorable prognosis not only in luminal subtype breast cancer but also in all subtype breast cancers. To the best of our knowledge, this is the first report indicating that FAM83H-AS1 was involved in luminal subtype breast cancer and was an independent prognostic indicator. Conclusion: Our study provides a rich resource to the research community for further identifying lncRNAs with diagnostic and therapeutic potentials and exploring biological function of lncRNAs in luminal subtype breast cancer. O r i g i n a l; r e s e a r c h - 8 1 0 2 l u J 2 1 n o 3 7 . 7 1 1 . 7 3 . 4 5 y b / m o c . s s re . .vdoep lsyeon Long non-coding RNAs (lncRNAs) are a class of noncoding RNA molecules the transcripts of which are longer than 200 nt.5 Although they are not protein-encoding, lncRNAs play crucial regulatory roles in a diverse range of cellular processes and biological pathways, including cancer initiation, progression and metastasis.6,7 Displaying tissue and cell specific expression, lncRNAs have great potential as biomarkers.8 Shen et al9 have recently reported novel lncRNAs in triple-negative breast cancer. LncRNA expression profile of HER-2-enriched subtype breast cancer has also been analyzed in our previous study.10 However, there is little information on the aberrant expressed lncRNAs in luminal subtype breast cancer. In this study, we aimed to uncover the dysregulated lncRNAs in luminal subtype breast cancer, which might offer potential biomarkers for prognosis evaluation and gene therapy. Materials and methods w u sample collection and rna extraction /ww lan /: o Three pairs of primary breast cancer tissues and adjacent ttshp rrspe non-cancerous tissues were collected in the Department rom oF of Surgical Oncology at the First Affiliated Hospital of f edd Wenzhou Medical University. Tissue samples were snaplnoa frozen in liquid nitrogen immediately after dissection and odw then stored at −80°C before RNA extraction. ERα and PR yap status of the samples was confirmed positive by postoperative reh immunohistochemistry. Total RNA was extracted from tissue dnT samples using TRIzol reagent (Invitrogen, Carlsbad, CA, trsaeag wUaSsAa)papcrcoovreddinbgy ttohetheethmicasncuofmacmtuirteter’esopfrtohteocFoirls.tTAhfisfislitautdeyd cnoT Hospital of Wenzhou Medical University. Informed conO sent was obtained from all individual participants included in the study. RNA sequencing and identification of differentially expressed lncRNAs The RNA with integrity number .7.0 was optimum for cDNA library construction. The cDNA libraries were then processed for the Proton Sequencing process according to the commercial protocols. Next, Ion Proton (Life Technologies) was used to perform single-end and polyAselection sequencing of the three pairs of samples. The MapSplice program (v2.1.6) was used to align clean reads to the human genome (version: GRCH37). We quantified the lncRNA expression as Reads Per Kilobase per Million mapped reads (RPKM) and lncRNAs with sum read counts ,10 across all samples were abandoned. EBseq algorithm 7040 was used to identify differentially expressed lncRNAs between cancer tissues and adjacent non-cancerous tissues. Differences in RNA expression were regarded as significant if values of the false discovery rate were ,0.05 and the fold change was .2. All differentially expressed lncRNAs are listed in Table S1. The cancer genome atlas (Tcga) data set The RPKM expression value of lncRNAs in TCGA database was downloaded through The Atlas of Non-coding RNAs in Cancer (TANRIC), which contained 837 breast cancer tissues and 105 non-tumorous tissues.11 Corresponding clinical parameters and follow-up information about these 837 patients were also downloaded from TCGA database.12 According to ERα and PR status, we screened out 626 samples as luminal subtype breast cancer. Clinical characteristics of the 626 luminal subtype patients as well as lncRNAs expression values are listed in Table S2. statistical analysis Mann–Whitney U-test was used to analyze the expression difference in lncRNAs between luminal subtype breast cancer tissues and adjacent non-cancerous tissues in TCGA cohort. Patients were divided into low or high lncRNA expression groups according to the median value. Kaplan–Meier and Cox regression analyses were utilized to assess the association between the lncRNA and overall survival of patients. A P-value ,0.05 was considered statistically significant. All statistical analyses were performed using SPSS Version 22.0 (Chicago, IL, USA). Results Differentially expressed lncRNAs in luminal subtype breast cancer Analysis of global transcriptome expression detected 11,727 lncRNAs (Figure S1). Volcano plot filtering found 796 lncRNAs that were significantly changed in luminal subtype breast cancer tissues compared with their adjacent non-cancerous tissues (Figure 1A). Among them, 436 lncRNAs were upregulated while 360 were downregulated. The hierarchical clustering analysis showed that with the expression of these genes, the samples could be clearly classified into two groups (Figure 1B). The top 20 differentially expressed lncRNAs are listed in Table 1. The results showed that FAM83H antisense RNA 1 (FAM83H-AS1) was the most upregulated lncRNA, whereas GSN antisense RNA 1 (GSN-AS1) was the most downregulated lncRNA. submit your manuscript | www.dovepress.com Dovepress 7041 8 1 0 2 l u J 2 1 n o 3 7 . 7 1 1 . 7 3 . 4 5 y b / m o c . s s re . .vdoep lsyeon w u /ww lan /: o tsp rse th rp rom oF f d e d a o l n w o d y p a r e h T d n a s t e g r a T o c n O Tcga cohort validation Mining published transcriptome sequencing data was of low cost and feasible for exploring gene expression. To validate the fidelity of our RNA sequencing data, we selected the most differentially expressed lncRNAs and detected these lncRNAs expression in TCGA cohort containing 626 luminal subtype breast cancer samples and 105 non-tumorous samples. We found that six lncRNAs were annotated in the TCGA and the results showed that FAM83H-AS1, ST8SIA6 antisense RNA 1 (ST8SIA6-AS1) and HOX transcript antisense RNA (HOTAIR) were upregulated significantly in cancer tissues while TRHDE antisense RNA 1 (TRHDEAS1), ALDH1L1 antisense RNA 2 (ALDH1L1-AS2) and PGM5 antisense RNA 1 (PGM5-AS1) were downregulated significantly (P,0.001, Figure 2). The consistent variation tendency of each lncRNA showed the robustness of our RNA sequencing. Identification of novel prognostic marker To determine whether lncRNAs could serve as prognostic markers in luminal subtype breast cancer, we analyzed the association between lncRNAs and prognosis in TCGA cohort. We found that high expression level of FAM83H-AS1 was correlated with unfavorable survival. As shown in Figure 3A, the difference between the survival curves of the two groups was statistically significant (P=0.004). The cumulative 10-year overall survival rates were 66.9% and 33.48% in the low expression group and the high expression group, respectively (P=0.012). With respect to the prognostic value of FAM83H-AS1, Cox regression analysis was conducted to adjust for clinical features. It indicated that age (hazards ratio [HR] =1.872, 95% CI =1.028–3.409, P=0.040), high FAM83H-AS1 expression (HR =2.008, 95% CI =1.230–3.544, P=0.006) and lymph node metastasis (HR =1.987, 95% CI =1.145–3.449, P=0.015) were 7042 8 1 0 2 l u J 2 1 n o 3 7 . 7 1 1 . 7 3 . 4 5 y b / m o c . s s re . .vdoep lsyeon w u /ww lan /: o tsp rse th rp rom oF f d e d a o l n w o d y p a r e h T d n a s t e g r a T o c n O distinctively linked with the prognosis of luminal subtype breast cancer in univariate Cox proportional hazards regression analysis (Table 2). However, high FAM83H-AS1 expression was the only independent prognostic factor after multivariate Cox proportional hazards regression analysis (HR =2.440, 95% CI =1.238–4.807, P=0.010). Furthermore, we found that FAM83H-AS1 was also a prognostic marker for all breast cancer types (P=0.028, Figure 3B). Taken together, these results indicated that FAM83H-AS1 was an independent prognostic indicator and might act as an onco-lncRNA for breast cancer. However, other differential expressed lncRNAs failed to distinguish high-risked luminal subtype breast cancer patients (Figure S2). Discussion Luminal subtype breast cancer, which accounts for two-thirds of breast cancer, maintains a more differentiated state and confers a more favorable outcome than other subtypes.4,13 The application of adjuvant endocrine therapy has contributed to a remarkable decrease in mortality during recent decades.14 However, breast cancer is highly heterogeneous and it is urgent to distinguish high-risked patients. Cheang et al reported that luminal subtype patients with high expression of HER-2 proteins and the Ki67 index had worse recurrencefree and disease-specific survival.15 It has been declared that among luminal subtype patients who received tamoxifen as their sole adjuvant systemic therapy, the 10-year breast cancer specific survival was 79% for luminal A subtypes and 57% for luminal-HER2 subtypes.15 Therefore, highly sensitive and specific biomarkers would be of great value in the individualized treatment for luminal subtype patients. Deregulation of lncRNAs has also been shown to contribute to the initiation and progression of several human cancers including breast cancer.16,17 A previous study showed that HOTAIR could promote cancer metastasis by reprogramming chromatin state and was a powerful predictor of eventual metastasis and death.18 In contrast, NF-kappaB interacting long non-coding RNA (NKILA) suppressed Univariate HR (95% CI) 1.872 (1.028–3.409) 0.975 (0.546–1.743) 1.987 (1.145–3.449) 1.228 (0.620–2.430) 0.352 (0.048–2.578) 0.538 (0.276–1.048) 1.726 (0.825–3.763) 2.008 (1.230–3.544) P-value 0.040* 0.933 0.015* 0.556 0.304 0.068 0.144 0.006** Multivariate HR (95% CI) P-value 2.440 (1.238–4.807) breast cancer metastasis by negatively regulating the NF-kB FGF14 antisense RNA 2 (FGF14-AS2) has also been proven to be correlated with progression and poorer prognosis in breast cancer in our previous study.20 Hence, increasing knowledge of lncRNAs might help develop novel therapeutic targets and prognostic indicators. To obtain comprehensive RNA expression profile data, we adopted RNA sequencing in three pairs of luminal subtype breast cancer tissues and their adjacent non-cancerous tissues. Genome-wide analysis revealed a set of lncRNAs with differential expression in cancer tissues compared with non-cancerous tissues. To the best of our knowledge, this is the only study focused on aberrant lncRNAs expression profiling in luminal subtype breast cancer. Furthermore, RNA sequencing data in TCGA database were utilized to validate the reliability of our results. Our study provided a rich resource to the research community for further identi.vdoep lsyeon fying lncRNAs with diagnostic and therapeutic potentials and exploring biological function of lncRNAs in luminal subtype breast cancer. Cabanski et al discovered 229 lncRNAs including FAM83H-AS1 with differential expression across multiple cancers and hypothesized that these lncRNAs might have conserved oncogenic and tumor-suppressive functions.21 Consistent with the previous study, we proved that FAM83HAS1 was significantly upregulated in breast cancer and was related to patients’ outcome in a large cohort. It was first reported that FAM83H-AS1 was involved in luminal subtype breast cancer and as an independent prognostic indicator. Further investigation into the functions of FAM83H-AS1 may provide additional target and strategies for treatment. Conclusion A total of 796 significantly differentially expressed lncRNAs in luminal subtype breast cancer were identified in the present study. Moreover, FAM83H-AS1 was demonstrated as a novel prognostic marker not only in luminal subtype breast cancer but also in all subtype breast cancers. Acknowledgments This study was funded by the Key Project of Science and Technology Innovation Team of Zhejiang Province (2013TD10) and the National Natural Science Foundation of China (no 81372380). We thank dan-dan Sun and her colleagues from Novel Bioinformatics Company for technical support in the bioinformatics analysis process. 8 1 0 2 l u J 2 1 n o 3 7 . 7 1 1 . 7 3 . 4 5 y b / m o c . s s re . pathway and was associated with poor outcome.19 Besides, The authors report no conflicts of interest in this work. Disclosure submit your manuscript | www.dovepress.com Dovepress Publish your work in this journal OncoTargets and Therapy is an international, peer-reviewed, open access journal focusing on the pathological basis of all cancers, potential targets for therapy and treatment protocols employed to improve the management of cancer patients. The journal also focuses on the impact of management programs and new therapeutic agents and protocols on 1. Siegel RL , Miller KD , Jemal A . Cancer statistics, 2015 . CA Cancer J Clin. 2015 ; 65 ( 1 ): 5 - 29 . 2. Howell A , Cuzick J , Baum M , et al. Results of the ATAC (arimidex, tamoxifen, alone or in combination) trial after completion of 5 years' adjuvant treatment for breast cancer . Lancet . 2005 ; 365 ( 9453 ): 60 - 62 . 3. Mauri D , Pavlidis N , Polyzos NP , Ioannidis JP . Survival with aromatase inhibitors and inactivators versus standard hormonal therapy in advanced breast cancer: meta-analysis . J Natl Cancer Inst . 2006 ; 98 ( 18 ): 1285 - 1291 . 4. Chen X , Cong Y , Pan L , et al. Luminal (Her2 negative) prognostic index and survival of breast cancer patients . Cancer Epidemiol . 2014 ; 38 ( 3 ): 286 - 290 . 5. St Laurent G , Wahlestedt C , Kapranov P. The landscape of long noncoding RNA classification . Trends Genet . 2015 ; 31 ( 5 ): 239 - 251 . 6. Yang G , Lu X , Yuan L. LncRNA : a link between RNA and cancer . Biochim Biophys Acta . 2014 ; 1839 (11): 1097 - 1109 . 7. Gutschner T , Diederichs S. The hallmarks of cancer: a long non-coding RNA point of view . RNA Biol . 2012 ; 9 ( 6 ): 703 - 719 . 8. Gloss BS , Dinger ME . The specificity of long noncoding RNA expression . Biochim Biophys Acta . 2016 ; 1859 (1): 16 - 22 . 9. Shen X , Xie B , Ma Z , et al. Identification of novel long non-coding RNAs in triple-negative breast cancer . Oncotarget . 2015 ; 6 ( 25 ): 21730 . 10. Yang F , Lyu S , Dong S , Liu Y , Zhang X , Wang O . Expression profile analysis of long noncoding RNA in HER-2-enriched subtype breast cancer by next-generation sequencing and bioinformatics . Onco Targets Ther . 2016 ; 9 : 761 - 772 . 11. Li J , Han L , Roebuck P , et al. TANRIC: an interactive open platform to explore the function of lncRNAs in cancer . Cancer Res . 2015 ; 75 ( 18 ): 3728 - 3737 . 12. Weinstein JN , Collisson EA , Mills GB , et al. The cancer genome atlas pan-cancer analysis project . Nat Genet . 2013 ; 45 ( 10 ): 1113 - 1120 . 13. Perou CM , Sørlie T , Eisen MB , et al. Molecular portraits of human breast tumours . Nature . 2000 ; 406 ( 6797 ): 747 - 752 . 14. Berry DA , Cronin KA , Plevritis SK , et al. Effect of screening and adjuvant therapy on mortality from breast cancer . N Engl J Med . 2005 ; 353 ( 17 ): 1784 - 1792 . 15. Cheang MC , Chia SK , Voduc D , et al. Ki67 index, HER2 status, and prognosis of patients with luminal B breast cancer . J Natl Cancer Inst . 2009 ; 101 ( 10 ): 736 - 750 . 16. Prensner JR , Chinnaiyan AM . The emergence of lncRNAs in cancer biology . Cancer Discov . 2011 ; 1 ( 5 ): 391 - 407 . 17. Spizzo R , Almeida MI , Colombatti A , Calin GA . Long non-coding RNAs and cancer: a new frontier of translational research ? Oncogene. 2012 ; 31 ( 43 ): 4577 - 4587 . 18. Gupta RA , Shah N , Wang KC , et al. Long non-coding RNA HOTAIR reprograms chromatin state to promote cancer metastasis . Nature . 2010 ; 464 ( 7291 ): 1071 - 1076 . 19. Liu B , Sun L , Liu Q , et al. A cytoplasmic NF-kappaB interacting long noncoding RNA blocks IkappaB phosphorylation and suppresses breast cancer metastasis . Cancer Cell . 2015 ; 27 ( 3 ): 370 - 381 . 20. Yang F , Liu Y-H , Dong S-Y , et al. A novel long non-coding RNA FGF14-AS2 is correlated with progression and prognosis in breast cancer . Biochem Biophys Res Commun . 2016 ; 470 ( 3 ): 479 - 483 . 21. Cabanski CR , White NM , Dang HX , et al. Pan-cancer transcriptome analysis reveals long noncoding RNAs with conserved function . RNA Biol . 2015 ; 12 ( 6 ): 628 - 642 .