Determination of long non-coding RNAs associated with EZH2 in neuroblastoma by RIP-seq, RNA-seq and ChIP-seq

Neuroblastoma (NB) is the most common type of extracranial solid tumor found in children. Despite several treatment options, patients with advanced stage disease have a poor prognosis. Previous studies have reported that enhancer of zeste homolog 2 (EZH2) and long non-coding RNAs (lncRNAs) have abnormal expression levels in NB and participate in tumorigenesis and NB development. However, the association between EZH2 and lncRNAs remain unclear. In the present study, RNA immunoprecipitation-sequencing (RIP-seq) was used to analyze the lncRNAs binding to EZH2. Following EZH2 knockdown via short hairpin RNA, RNA-seq was performed in shEZH2 and control groups in SH-SY5Y cells. Chromatin IP (ChIP)-seq was used to determine the genes that may be regulated by EZH2. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses were performed to identify the signaling pathways involved in NB. The results from RIP-seq identified 94 lncRNAs, including SNHG7, SNHG22, KTN-AS1 and Linc00843. Furthermore, results from RNA-seq demonstrated that, following EZH2 knockdown, 448 genes were up- and 571 genes were downregulated, with 32 lncRNAs up- and 35 downregulated and differentially expressed compared with control groups. Certain lncRNAs, including MALAT1, H19, Linc01021 and SNHG5, were differentially expressed in EZH2-knockdown group compared with the control group. ChIP-seq identified EZH2 located in the promoter region of 138 lncRNAs including CASC16, CASC15, LINC00694 and TBX5-AS1. In summary, the present study demonstrated that certain lncRNAs directly bound EZH2 and regulated EZH2 expression levels. A number of these lncRNAs that are associated with EZH2 may participate in NB tumorigenesis.