- Open Access
A novel long noncoding RNA HOXC-AS3 mediates tumorigenesis of gastric cancer by binding to YBX1
Genome Biologyvolume 19, Article number: 154 (2018)
Recently, increasing evidence shows that long noncoding RNAs (lncRNAs) play a significant role in human tumorigenesis. However, the function of lncRNAs in human gastric cancer remains largely unknown.
By using publicly available expression profiling data from gastric cancer and integrating bioinformatics analyses, we screen and identify a novel lncRNA, HOXC-AS3. HOXC-AS3 is significantly increased in gastric cancer tissues and is correlated with clinical outcomes of gastric cancer. In addition, HOXC-AS3 regulates cell proliferation and migration both in vitro and in vivo. RNA-seq analysis reveals that HOXC-AS3 knockdown preferentially affects genes that are linked to proliferation and migration. Mechanistically, we find that HOXC-AS3 is obviously activated by gain of H3K4me3 and H3K27ac, both in cells and in tissues. RNA pull-down mass spectrometry analysis identifies that YBX1 interacts with HOXC-AS3, and RNA-seq analysis finds a marked overlap in genes differentially expressed after YBX1 knockdown and those transcriptionally regulated by HOXC-AS3, suggesting that YBX1 participates in HOXC-AS3-mediated gene transcriptional regulation in the tumorigenesis of gastric cancer.
Together, our data demonstrate that abnormal histone modification-activated HOXC-AS3 may play important roles in gastric cancer oncogenesis and may serve as a target for gastric cancer diagnosis and therapy.
Gastric cancer is one of the leading causes of cancer-related deaths worldwide and the most common gastrointestinal malignancy in East Asia [1, 2]. Gastric cancer is diagnosed at an advanced stage accompanied by malignant proliferation in most patients, and the prognosis for advanced stage patients remains very poor . Therefore, additional research is needed to discover and develop effective biomarkers and targets for gastric cancer diagnosis and treatment.
To date, gastric cancer research has mainly focused on the deregulation of protein-coding genes to identify oncogenes and tumor suppressors that could serve as diagnostic and therapeutic targets. However, protein-coding sequences occupy less than 2% of the human genome [4, 5]. LncRNAs are operationally defined as RNA transcripts that are > 200 nt with limited protein coding potential , which have been shown to play a key role in tumorigenesis, including GC [7,8,9]. Many studies found that lncRNAs could play an important role in regulating gene expression by different mechanisms, including chromatin modification, and transcriptional and posttranscriptional processing [10,11,12]. For example, HOTAIR is involved in the transcriptional repression of HOX loci and promotes breast metastasis by binding to PRC2 (Polycomb Repressive Complex) . However, the biological functions of lncRNAs in the control of GC tumorigenesis are not well characterized. Therefore, a better understanding of the role of lncRNAs underlying GC progression will enrich the understanding of the molecular mechanisms of GC carcinogenesis and provide information for improving the diagnosis and treatment of GC.
In our present study, we identified the full sequence of HOXC-AS3 and found that gain of H3K4me3 and H3K27acetylation could activate the expression of HOXC-AS3, both in cells and in tissues. HOXC-AS3 was also significantly upregulated in GC tissues compared with the corresponding nontumor tissues and may serve as an independent predictor for the overall survival in GC. In addition, HOXC-AS3 regulated cell proliferation and migration both in vitro and in vivo. RNA-seq analysis for whole transcriptome studies indicates an important role for HOXC-AS3 in the tumorigenesis of GC, and the activated function of HOXC-AS3 was mediated, in part, by interaction with YBX1. These results suggest that further studies to identify nonprotein-coding genes that contribute to oncogenesis are necessary for elucidating the complex genetic rewiring that is driven by HOXC-AS3 in GC.
Identification of HOXC-AS3 by analyzing gastric cancer RNA-expression profiling data
Raw microarray data was downloaded from GEO, including GSE50710 (n = 20) and GSE58828 (n = 6). To obtain differentially expressed lncRNAs, signal data were normalized and z-score-transformed (paired or group t-test according to the experimental design was used to validate statistical significance). The top 20 dysregulated lncRNAs in these GEO datasets were shown in Fig. 1a, b. For example, the lncRNA that showed clear upregulation was HOXA11-AS (the largest increase in GSE50710), and our previous study showed that HOXA11-AS could regulate cell proliferation and invasion of gastric cancer by scaffolding chromatin modification factors and patients with high HOXA11-AS expression had a shorter survival and poorer prognosis . The most obvious increase in GSE58828, LINC00473, could mediate tumor growth and elevated LINC00473 expression correlated with poor prognosis of lung cancer . Our previous study showed that FEZF1-AS1 could promote gastric cancer proliferation and higher expression of FEZF1-AS1 predicted poor prognosis . Qin et al. found that LINC01296 was upregulated in GC tissue and correlated with poor prognosis . CCAT1 is an oncogenic long noncoding RNA in human cancers including GC . Our previous work found that CCAT1 could regulate cell proliferation and migration in esophageal squamous cell carcinoma and higher expression of CCAT1 is correlated with poorer outcomes . In addition, CCAT1 could serve as important prognostic biomarkers in colorectal cancer . We focused on overexpressed lncRNAs because these lncRNAs may be more readily used as early diagnosis markers or therapeutic targets. Therefore, although a small portion of lncRNAs have been functionally characterized, many members in the class remain uncharacterized. In our analysis of the results, HOXC-AS3 was the lncRNA clearly upregulated in both GSE50710 and GSE58828 (the rank of HOXC-AS3 is in the exceeding front position, both in GSE50710 and GSE58828), and the role in GC tumorigenesis including other cancer types remained unclear. This prompted us to explore the role of HOXC-AS3 in human GC. To further validate this result, we analyzed that RNA-Seq data (from The Cancer Genome Atlas (TCGA)) of lncRNAs of gastric cancer were from TANRIC  (http://ibl.mdanderson.org/tanric/_design/basic/index.html). As shown in Fig. 2a, HOXC-AS3 was significantly upregulated in GC tissues from the TCGA data.
Gain of H3K4me3 and H3K27 acetylation-activated HOXC-AS3 is upregulated in human GC tissues and correlates with poor prognosis
For further study, we first performed rapid amplification of cDNA ends (RACE) to identify the full sequence of HOXC-AS3 in BGC-823 cells according the sequence archived in the RefSeq database of NCBI (440 bp, with poly (A) tail, Fig. 2b). To validate the expression results from high-throughput data, as shown in Fig. 2c, the HOXC-AS3 expression level in tumor tissues was significantly higher in 112 pairs of GC tumor tissues compared with adjacent normal tissues (P < 0.001). One tumor tissue showed an upregulation of HOXC-AS3 greater than 521-fold relative to normal tissue. Next, we explored the correlation between HOXC-AS3 expression and the clinicopathological factors of patients with GC. The result showed that the HOXC-AS3 level was associated with TNM stage. Patients with advanced TNM stage (III/IV) were associated with higher HOXC-AS3 expression, whereas patients with local TNM stage (I/II) were associated with a lower HOXC-AS3 level (27.1317 ± 80.26254 vs 4.7483 ± 3.02402, P = 0.046) (Fig. 2d). Furthermore, we divided the samples into relatively high (above the mean, n = 56) and relatively low (below the mean, n = 56) HOXC-AS3 expression groups according to the median value of HOXC-AS3 levels. A chi-square test was then performed to evaluate clinicopathological factors between the two groups. As shown in Additional file 1: Table S1, the relative HOXC-AS3 level was also correlated with histological grade (P = 0.002), tumor invasion depth (P = 0.008), lymph node metastasis (P = 0.035), and TNM stage (P = 0.002). No relationship between HOXC-AS3 expression and other clinical factors, such as sex (male, female) and patient age (≤ 60, > 60), was found in our study.
To determine the relationship between HOXC-AS3 expression and GC patient prognosis, we evaluated the correlation between HOXC-AS3 expression and clinical outcomes. Kaplan–Meier analysis and log-rank test were used to evaluate the effects of HOXC-AS3 expression and the clinicopathological characteristics on overall survival (OS). The median survival time for low HOXC-AS3 expression groups was 34 months, whereas for high HOXC-AS3 expression groups, it was only 20 ± 1.357 months. As shown in Fig. 2e, overexpression of HOXC-AS3 predicted a poor prognosis in patients with GC (P = 0.004). Similarly, the correlation between HOXC-AS3 expression levels and the survival of GC patients was also supported by Kaplan–Meier Plotter analysis (http://kmplot.com/analysis/, detailed steps were described in Additional file 8: Supplementary Methods), which indicated that higher HOXC-AS3 expression correlated with worse OS, using publicly available chip data from 631 GC patients (Fig. 2f).
Then, univariate and multivariate survival analyses (Cox proportional hazards regression model) were performed. Univariate analysis identified two prognostic factors: TNM stage and HOXC-AS3 expression. Multivariate analysis further revealed that HOXC-AS3 expression was an independent predictor for overall survival in patients with GC (P < 0.001), as well as TNM stage (P = 0.023) (Additional file 2: Table S2).
To explore the mechanism of high expression of HOXC-AS3 in GC, firstly, by using the UCSC Genome Bioinformatics Site (http://genome.ucsc.edu/), we found high enrichment and overlaps of H3K4me3 and H3K27Ac peaks at the promoter region of HOXC-AS3 (Fig. 2g, H3K4me3 and H3K27Ac, two markers of active promoters). Using ChIP assays, we found gain of H3K4me3 and H3K27Ac in cancer tissues compared with normal tissues (n = 4) at the promoter of HOXC-AS3. We also observed the gain of H3K4me3 and H3K27Ac in GC cells (BGC-823) compared with normal human esophageal epithelial cells (GES-1) at the promoter of HOXC-AS3 (Fig. 2g). Taken together, these data confirm that HOXC-AS3 is frequently increased in GC. Abnormal histone modification, gain of histone sites H3K4me3 and H3K27Ac of the promoter may partially account for the significant activation of HOXC-AS3.
HOXC-AS3 regulates GC cell proliferation and migration in vitro
To explore the role of HOXC-AS3 in GC, we examined the HOXC-AS3 expression levels in gastric cancer cell lines. As shown in Fig. 3a, gastric cancer cell lines expressed significantly higher levels of HOXC-AS3 than a normal gastric epithelium cell line (GES-1). To further confirm the function of HOXC-AS3, we used an LNA-ASO (Locked Nucleic Acid, antisense oligonucleotide) targeting HOXC-AS3 and control LNA-ASO. As shown in (Additional file 3: Figure S1A), the ASO-mediated knockdown and plasmid-mediated overexpression were used for exogenously manipulating the expression of HOXC-AS3 in BGC-823 and SGC-7901 cells. Then, MTT assays showed that knockdown of HOXC-AS3 expression significantly inhibited cell proliferation compared with control cells. By contrast, overexpressed HOXC-AS3 promoted cell proliferation (Fig. 3b). EdU assays demonstrated that HOXC-AS3 had a significant impact on GC cell proliferation (Fig. 3c). Similarly, the results of the colony-formation assay revealed that clonogenic survival was significantly decreased following knockdown of HOXC-AS3. Furthermore, overexpression of HOXC-AS3 could increase the number of clones (Fig. 3d). Next, transwell assays revealed that knockdown of HOXC-AS3 significantly repressed cell migration compared with the control in BGC-823 cells. By contrast, overexpression of HOXC-AS3 promoted cell migration (Fig. 3e). Flow cytometric analysis was performed to further examine whether HOXC-AS3 affected the proliferation of GC cells by altering cell-cycle progression. The results revealed that cell-cycle progression of knockdown-HOXC-AS3 cells was significantly stalled at the G1-G0 phase compared with controls both in BGC-823 and in SGC-7901 cells. Overexpressed HOXC-AS3 promoted S phase progression (Fig. 3f). Additionally, knockdown of HOXC-AS3 significantly induced BGC-823 and SGC-7901 cell apoptosis (Fig. 3g).
HOXC-AS3 regulates GC cell proliferation and migration in vivo
To further determine whether HOXC-AS3 affects cell proliferation of GC in vivo, BGC-823 cells and SGC-7901 stably transfected with sh-HOXC-AS3 or control vector were inoculated into nude mice. Sixteen days after the injection, the tumors formed in the sh-HOXC-AS3 group were substantially smaller than those in the control group (Fig. 4a). Moreover, the mean tumor weight at the end of the experiment was markedly lower in the sh-HOXC-AS3 group compared with the control vector group (Fig. 4a). Tumors formed from stably sh-HOXC-AS3-transfected BGC-823 and SGC-7901 cells exhibited decreased positivity for Ki-67 than tumors from the control cells (Fig. 4a). These findings indicate that knockdown of HOXC-AS3 inhibits tumor growth in vivo.
To validate the effects of HOXC-AS3 on cell metastasis in vivo, BGC-823 cells stably transfected with sh-HOXC-AS3 or control vector were injected into the tail veins of nine mice. Metastatic nodules on the surface of the lungs were counted after 7 weeks. Ectopic knockdown of HOXC-AS3 reduced the number of metastatic nodules compared with the control group (Fig. 4b). This difference was further confirmed following examination of the entire lungs and by hematoxylin and eosin (HE) staining of lung sections (Fig. 4b). Our in vivo data, therefore, complemented the results of the functional in vitro studies involving HOXC-AS3.
HOXC-AS3 interacts with YBX1
To explore the mechanism of HOXC-AS3, firstly, we examined whether HOXC-AS3 acts in cis, affecting nearby gene expression (HOXC10; HOXC-AS3 was an antisense transcript of HOXC10). We found that knockdown of HOXC-AS3 did not affect expression of HOXC10 (data not shown), indicating it may act in trans. Then, subcellular fractionation location assays showed a considerable increase in HOXC-AS3 expression in the nucleus versus the cytosol (Fig. 5a), thus suggesting that HOXC-AS3 may play a major regulatory function at the transcriptional level, indicating it might interact with nucleus molecules or proteins. Then, we performed an RNA pull-down assay followed by a proteomic analysis of the HOXC-AS3-associated protein complex in BGC-823 cells (Fig. 5b). We incubated the in vitro transcribed HOXC-AS3 bound to beads with BGC-823 nuclear extracts to purify the HOXC-AS3 RNA-protein complex, and protein identity was revealed by mass spectrometry. Among the highly enriched proteins (Additional file 4: Table S3), only YBX1 was detected by western blotting from three independent RNA pull-down assays. YBX1 attracted our attention because of its established role in tumorigenesis and this notable protein was a known transcription factor [22, 23], with 12 unique peptides detected in this MS analysis. To further validate the physical interaction between HOXC-AS3 and YBX1, we performed RNA pull-down followed by western blotting with YBX1 antibodies. Our results showed that labeled HOXC-AS3 RNA, but not empty vector or an antisense HOXC-AS3, specifically retrieved YBX1 from BGC-823 cell extracts (Fig. 5b). We also performed RNA immunoprecipitation (RIP) for the RNA-YBX1 complex using different YBX1 antibodies and measured the amount of HOXC-AS3 associated with YBX1 immunoprecipitates but not MEG3. LncRNA MEG3 was a negative control (Fig. 5c). Using a series of HOXC-AS3 deletion mapping, the YBX1-binding activity mapped to nucleotides 115 nt of HOXC-AS3 (Fig. 5d). Additionally, the qRT-PCR and western blots showed that knockdown of HOXC-AS3 could have had no effect on the expression of YBX1 (Fig. 5e). Our results proved that HOXC-AS3 could bind to YBX1 but have no effect on YBX1 expression. These results indicated that HOXC-AS3 may participate in the tumorigenesis of GC through transcriptional regulation of other genes via binding to YBX1.
A large set of genes that are linked to cell proliferation and cell migration were coregulated by interaction of HOXC-AS3 and YBX1
To ascertain the mechanism of action of HOXC-AS3/YBX1-driven oncogenesis of GC and to determine the gene expression changes downstream of HOXC-AS3, we evaluated the global effects of HOXC-AS3 knockdown compared with those of YBX1 knockdown using RNA transcriptome sequencing (Fig. 6a, Additional file 5: Table S4 and Additional file 6: Table S5). A total of 3122 genes were significantly upregulated or downregulated in BGC-823 cells as a consequence of HOXC-AS3 knockdown (fold change ≥ 1.5) (Fig. 6a). For RNA-seq results of HOXC-AS3 knockdown, gene ontology (GO) analysis showed that the most significantly overrepresented biological processes included pathways involved in cell division, cell proliferation, and cell death (Fig. 6b). Gene set enrichment analysis (GSEA) revealed that the gene sets were significantly related to cell proliferation and metastasis (Fig. 6c). We also built a coexpression network based on the differentially expressed genes that are functionally related. The co-expression network was constructed according to the gene normalized expression (RPKM), and the k-core was calculated based on Pearson’s relationship to represent the core status of each gene in different groups. The k-core factor was then applied to identify key regulatory genes, which likely play pivotal roles in gene interactions and regulation (Fig. 6d), including for key cancer genes, such as p21, FAS, CCND1, and CDK2; real-time PCR confirmed these genes after knockdown of HOXC-AS3 (Additional file 3: Figure S2B). Additionally, a total of 1480 genes were significantly upregulated or downregulated in BGC-823 cells as a consequence of YBX1 knockdown (fold change ≥ 1.5) (Fig. 6a and Additional file 3: Figure S2A). Intriguingly, the expression of 805 of the 3122 HOXC-AS3-regulated genes was also regulated by knockdown of YBX1 (Fig. 6a). The high overlap between HOXC-AS3 and YBX1-mediated regulation strongly indicates a functional interaction between HOXC-AS3 and YBX1. YBX1 functions as an RNA-binding protein and has been implicated in numerous cellular processes, including the regulation of transcription and translation . YBX1 was shown to be a protein with a nucleic acid-binding common domain in gene promoter, CCAAT-box, which is a high consensus sequence in eukaryotes . Moreover, YBX1 serves as a transcriptional activator and YBX1 activation was associated with cancer progression, including gastric cancer [23, 26, 27]. Our results showed that YBX1 was significantly upregulated in 60 pairs of GC tissues by qRT-PCR analysis. Knockdown of YBX1 could significantly inhibit cell proliferation and migration (Additional file 3: Figure S2C). Previous research supports a role for lncRNA-YBX1 interactions as mediators of gene transcriptional regulation . In particular, 476 genes were corepressed by knockdown of HOXC-AS3 and YBX1. These results indicated that a common set of target genes linked to cell proliferation and cell migration were shared by HOXC-AS3 and YBX1.
Genes were corepressed by knockdown of HOXC-AS3/YBX1. They included many well-known proliferation and migration-related genes (e.g., CDK2, HOXB13, IGFBP4, ATF5, MAPK4, MMP7, MMP24, BIRC2, WNT10B, and HDAC5). These genes were selectively confirmed by qRT-PCR after knockdown of HOXC-AS3 and YBX1 in BGC-823 and SGC-7901 cells (all these genes containing CCAAT-box, YBX1-binding site in their promoters) (Fig. 7a). Given that YBX1 regulates gene transcription by binding to the promoter regions of target genes, we examined whether HOXC-AS3 knockdown affected YBX1 occupancy of the promoter regions in these target genes. Then ChIP assays followed by qPCR demonstrated that knockdown of HOXC-AS3 decreased the binding of YBX1 levels across the promoters of the most coregulated genes (Fig. 7b; among the ten target genes tested, YBX1 occupancy were indeed affected in the promoter regions of seven genes following knockdown of HOXC-AS3). These findings demonstrate that HOXC-AS3 is important in regulating gene transcription expression, presumably in part by regulating the association between YBX1 and the promoter regions of its target genes.
HOXC-AS3 regulates the expression of HDAC5
Among the common targets of HOXC-AS3 and YBX1, HDAC5 is of particular interest because of its remarkable expression fold change upon HOXC-AS3 knockdown (Fig. 7a) and its significant contribution to tumorigenesis . Aberrant activation of HDAC5 in tumor cells leads to dysregulation of a diverse set of genes mainly involved in the regulation of proliferation, migration, and angiogenesis . Our results showed that HDAC5 was significantly upregulated in 60 pairs of GC tissues. Further analysis demonstrated that HOXC-AS3 was positively correlated with HDAC5 expression in GC tissues (Fig. 8a, b). Then, we verified the results at the protein level for HDAC5 after knockdown of HOXC-AS3 (Fig. 8c). Furthermore, the increase of HDAC5 by HOXC-AS3 was reversed with treatment by YBX1 siRNAs at both the mRNA and protein levels (Fig. 8c). Moreover, knockdown of HDAC5 induced the suppression of proliferation and migration in BGC-823 cells (Fig. 8d and Additional file 3: Figure S2D). In addition, knockdown of HDAC5 could reverse HOXC-AS3-mediated growth and migration promotion (Fig. 8d). Together, these results demonstrate that HOXC-AS3 exerts its function, at least in part by regulating HDAC5 expression.
Before the discovery of noncoding RNAs, explorations for cancer drivers focused on protein-coding genes that resided in recurrent alterations in cancer genomes. However, the newly discovered lncRNAs have emerged as important players in cellular development and human diseases, especially in cancer. In the present study, utilizing publicly available lncRNA expression profiling data of gastric cancer and integrating analyses of TCGA data, we screened and identified a novel lncRNAs HOXC-AS3. The high expression of HOXC-AS3 in GC patients was positively correlated with advanced TNM stage. Moreover, high HOXC-AS3 expression in GC tissues was associated with a poor prognosis and could be an independent prognostic indicator. In addition, GTEx data (https://www.gtexportal.org/) showed that HOXC-AS3 had lower basal expression in normal GC tissues (Additional file 3: Figure S1B). This further demonstrates an important role of HOXC-AS3 in the carcinogenesis of GC. These results and our functional evidence for HOXC-AS3 suggested that HOXC-AS3 might exhibit an important role in GC progression.
HOXC-AS3, which is located at chromosome 12q13.13, was an antisense transcript of HOXC10. HOX genes are essential for morphogenesis and development [31, 32]. The dysregulation of HOX gene expression has been shown in many diverse cancers [33, 34], and lncRNA generation in HOX genes may play important roles in tumorigenesis. For example, as a well-known lncRNA, HOTAIR, also located at chromosome 12q13.13, is an antisense transcript of HOXC11 and was an oncogenic lncRNA in many different types of cancer . Our previous study also found that a lncRNA of the HOX gene family, HOXA11-AS, could play an important role in GC tumorigenesis . In our present study, we found a novel lncRNA in the HOX genes family. Our results revealed that gain of H3K4me3 and H3K27Ac activation of the promoter also partly contributed to activation of HOXC-AS3 in GC, both in cells and tissues. Similar to protein coding transcripts, the transcription of lncRNAs is subject to typical epigenetics-mediated and transcription factor-mediated regulation. For example, the lncRNA MEG3 was lost in tumors due to an increase in CpG methylation within the promoter . Histone deacetylase3-suppressed the lncRNA LET in hepatocellular carcinoma by reducing the histone acetylation-mediated modulation of the promoter region .
In our study, we found that inhibition of HOXC-AS3 repressed GC proliferation and migration both in vitro and in vivo. RNA-seq found that knockdown of HOXC-AS3 affected key cancer-related genes, such as p21, FAS, and CCND1. Mechanistic investigations found that HOXC-AS3 could bind to YBX1, but not affect YBX1 expression. These results indicated that HOXC-AS3 may participate in the tumorigenesis of GC through the transcriptional regulation of other genes via binding to YBX1 in trans. To probe the HOXC-AS3-associated pathway on an unbiased basis in the tumorigenesis of GC, RNA-Seq assays were used after simultaneous knockdown HOXC-AS3 and YBX1. Notably, loss of HOXC-AS3 in GC cells recapitulated the phenotype observed after YBX1 knockdown. In addition, a significant fraction of the genes regulated by HOXC-AS3 loss were similarly regulated by loss of YBX1. Thus, HOXC-AS3 may act, in part, by regulating the interaction between YBX1 and the promoter of target genes, although HOXC-AS3 likely interacts with other RNA-binding proteins that will need to be identified to fully understand its molecular function. YBX1 was reported to play a role in regulating cell signaling, transcription, and tumorigenesis . Furthermore, YBX1 serves as a transcriptional activator and regulates much gene transcription . YBX1 was a protein with a nucleic acid-binding common domain in the gene promoter, CCAAT-box, which is a high consensus sequence in eukaryotes. We found that YBX1 is overactive in GC and knockdown of YBX1 inhibits the proliferation of GC cells. Our results showed that HOXC-AS3 could bind to YBX1, thus transcriptionally regulating a large set of genes that are linked to cell proliferation and cell migration in gastric cancer cells, such as MMP7, WNT10B, and HDAC5, thus promoting GC cell proliferation and migration.
The HDAC5 gene is a member of the histone deacetylase (HDAC) from a family of enzymes. Histone acetylation and deacetylation play important roles in chromatin remodeling and gene expression. An imbalance of these reactions leads to the growth, migration, and apoptosis of cancer cells. Histone deacetylase (HDAC) inhibitors were shown to have antitumor effects in clinical trials [39, 40]. Over-activation of HDAC5 was found in many different types of cancer [41,42,43]. A previous study showed that HDAC5 was induced in gastric cancer cells . Here, we also provide evidence for high expression of HDAC5 in gastric cancer, and knockdown of HDAC5 inhibits the proliferation of GC cells. We also found that the transcriptional activation of HDAC5 is partly mediated by HOXC-AS3 in the tumor progression of GC through binding to YBX1, thus facilitating GC cell proliferation and migration. In addition to HDAC5, there are many other important genes related to tumorigenesis, and they are also regulated by HOXC-AS3 in a similar manner.
In summary, abnormal histone modification-mediated activation of a novel lncRNA HOXC-AS3 promotes GC cell proliferation and migration through transcriptional activation of a large set of genes through an interaction with YBX1. Our data reveal a role for HOXC-AS3 in GC tumorigenesis and may provide a strategy for using HOXC-AS3 as a potential biomarker and a therapeutic target for patients with GC (Fig. 8e).
Tissue collection and ethics statement
A total of 112 patients in this study underwent resection of the GC at The Affiliated Jiangyin Hospital of Southeast University Medical College, affiliated Xuzhou Central Hospital of Southeast University Medical College. The study was approved by the Medical Ethical Committee of Southeast University Medical College (Nanjing, Jiangsu, PR China), and it was performed in compliance with the Helsinki Declaration. All patients have given written informed consent for publication. The clinicopathological characteristics of the GC patients are summarized in Additional file 1: Table S1.
Gastric cancer RNA-expression data retrieval and analysis
Microarray data analysis
Microarray datasets from the GEO database were used to test HOXC-AS3 differential expression. Raw microarray data was downloaded from GEO including GSE50710 and GSE58828. Then, the raw microarray data were normalized and z-score-transformed using RMAExpress (http://www.rmaexpress.bmbolstad.com/). RNA-Seq data (from TCGA) of lncRNAs of gastric cancer were from TANRIC (http://ibl.mdanderson.org/tanric/_design/basic/index.html) .
RACE (rapid amplification of cDNA ends)
5′-RACE, 3′-RACE, and full-length amplification of HOXC-AS3 were performed using the SMART RACE cDNA Amplification Kit (Cat. 634858, Clontech, Palo Alto, CA, USA) according to the manufacturer’s instructions. The gene-specific primers used for RACE analysis were presented in Additional file 7: Table S6.
Transfection of cell lines
LNA-ASO (Locked Nucleic Acid, antisense oligonucleotide) targeting HOXC-AS3 and negative control LNA-ASO were designed and synthesized by Exiqon (Exiqon, Vedbaek, Denmark). GC cells were transfected with the LNA-ASOs using Oligofectamine transfection reagent (RNAi MAX, Invitrogen) according to the manufacturer’s instructions. Cells were harvested for analyses 48 h after transfections. The sequences of ASO were listed in Additional file 7: Table S6. The sequences for siRNAs were listed in Additional file 7: Table S6. Scrambled negative control siRNA was purchased from Invitrogen (Invitrogen, CA, USA). The interference target sequence of YBX1 was acquired according to a previous study . The HDAC5 siRNA was from Santa Cruz (sc-35542). The plasmid was transfected into GC cells using the X-tremeGENE™ HP DNA Transfection Reagent (Roche) according to the manufacturer’s instructions.
Subcellular fractionation location
Separation of the nuclear and cytosolic fractions was performed using the PARIS Kit (Cat. AM1921, Invitrogen, CA, USA) according to the manufacturer’s instructions.
In vitro transcription assays and RNA pull-down mass spectrometry (LC-MS/MS) assays
In vitro translation assays were performed using mMESSAGE mMACHINE™ T7 Transcription Kit according to the manufacturer’s instructions (Cat. AM1344, Invitrogen, CA, USA). Then, HOXC-AS3 RNAs were labeled with desthiobiotinylation using the Pierce RNA 3′ End Desthiobiotinylation Kit (Cat. 20164, Magnetic RNA-Protein Pull-Down Kit, Components, Thermo). RNA pull-down assays were performed with Pierce Magnetic RNA-Protein Pull-Down Kit according to the manufacturer’s instructions (Cat. 20164, Magnetic RNA-Protein Pull-Down Kit, Thermo). After elution of lncRNA-interacting proteins, they were subjected to mass spectrometric analysis. LC-MS/MS experiments were performed with an LTQ linear ion trap mass spectrometer (Thermo Finnigan, San Jose, CA) equipped with a microspray source.
RNA immunoprecipitation (RIP) assays
RNA immunoprecipitation (RIP) experiments were performed using a Magna RIP™ RNA-Binding Protein Immunoprecipitation Kit (Cat. 17-701, Millipore, USA) according to the manufacturer’s instructions. The antibodies for RIP assays of YBX1 (Cat. ab12148, ab76149) were from Abcam.
Total RNA from BGC-823 cells with HOXC-AS3/YBX1 knockdown and control cells were isolated and quantified. The concentration of each sample was measured with a NanoDrop 2000 (Thermo Scientific, USA). The quality was assessed by an Agilent2200 (Agilent, USA). The sequencing library of each RNA sample was prepared using the Ion Proton Total RNA-Seq Kit v2 according to the protocol provided by the manufacturer (Life Technologies, USA). Data are available in Additional file 5: Table S4 and Additional file 6: Table S5.
Chromatin immunoprecipitation (ChIP) assays
ChIP assays were performed using the EZ-CHIP KIT according to the manufacturer’s instruction (Cat. 17-408, Millipore, USA). The antibodies for Histone H3, acetyl-histone H3 Lys27 (H3K27Ac, Cat. ab4729), and H3 trimethyl Lys4 (H3K4me3, Cat. ab8580) were from Abcam. The ChIP primer sequences were listed in Additional file 7: Table S6. The antibody for YBX1 (Cat. ab12148) was from Abcam. Quantification of immunoprecipitated DNA was performed using qPCR. ChIP data was calculated as a percentage relative to the input DNA from the equation 2[Input Ct − Target Ct] × 100 (%).
All statistical analyses were performed using SPSS 20.0 software (IBM, SPSS, USA). The significance of differences between groups was estimated by Student’s t test, χ2 test, or Wilcoxon test, as appropriate. OS rates were calculated by the Kaplan–Meier method with the log-rank test for comparison. Survival data were evaluated using univariate and multivariate Cox proportional hazards model. Variables with a value of P < 0.05 in univariate analysis were used in the subsequent multivariate analysis based on the Cox regression analyses. Two-sided P values were calculated, and a probability of 0.05 was selected for statistical significance.
Additional methods are described in Additional file 8: Supplementary Methods.
Siegel R, Naishadham D, Jemal A. Cancer statistics, 2013. CA Cancer J Clin. 2013;63(1):11–30.
Sugano K. Screening of gastric cancer in Asia. Best Pract Res Clin Gastroenterol. 2015;29(6):895–905.
Catalano V, Labianca R, Beretta GD, Gatta G, de Braud F, Van Cutsem E. Gastric cancer. Crit Rev Oncol Hematol. 2009;71(2):127–64.
Guttman M, Amit I, Garber M, French C, Lin MF, Feldser D, Huarte M, Zuk O, Carey BW, Cassady JP, et al. Chromatin signature reveals over a thousand highly conserved large non-coding RNAs in mammals. Nature. 2009;458(7235):223–7.
Cesana M, Cacchiarelli D, Legnini I, Santini T, Sthandier O, Chinappi M, Tramontano A, Bozzoni I. A long noncoding RNA controls muscle differentiation by functioning as a competing endogenous RNA. Cell. 2011;147(2):358–69.
Nagano T, Fraser P. No-nonsense functions for long noncoding RNAs. Cell. 2011;145(2):178–81.
Prensner JR, Chinnaiyan AM. The emergence of lncRNAs in cancer biology. Cancer Discov. 2011;1(5):391–407.
Hu Y, Wang J, Qian J, Kong X, Tang J, Wang Y, Chen H, Hong J, Zou W, Chen Y, et al. Long noncoding RNA GAPLINC regulates CD44-dependent cell invasiveness and associates with poor prognosis of gastric cancer. Cancer Res. 2014;74(23):6890–902.
Sun TT, He J, Liang Q, Ren LL, Yan TT, Yu TC, Tang JY, Bao YJ, Hu Y, Lin Y, et al. LncRNA GClnc1 promotes gastric carcinogenesis and may act as a modular scaffold of WDR5 and KAT2A complexes to specify the histone modification pattern. Cancer Discov. 2016;6(7):784–801.
Mercer TR, Dinger ME, Mattick JS. Long non-coding RNAs: insights into functions. Nat Rev Genet. 2009;10(3):155–9.
Wilusz JE, Sunwoo H, Spector DL. Long noncoding RNAs: functional surprises from the RNA world. Genes Dev. 2009;23(13):1494–504.
Hu X, Feng Y, Zhang D, Zhao SD, Hu Z, Greshock J, Zhang Y, Yang L, Zhong X, Wang LP, et al. A functional genomic approach identifies FAL1 as an oncogenic long noncoding RNA that associates with BMI1 and represses p21 expression in cancer. Cancer Cell. 2014;26(3):344–57.
Gupta RA, Shah N, Wang KC, Kim J, Horlings HM, Wong DJ, Tsai MC, Hung T, Argani P, Rinn JL, et al. Long non-coding RNA HOTAIR reprograms chromatin state to promote cancer metastasis. Nature. 2010;464(7291):1071–6.
Sun M, Nie F, Wang Y, Zhang Z, Hou J, He D, Xie M, Xu L, De W, Wang Z, et al. LncRNA HOXA11-AS promotes proliferation and invasion of gastric cancer by scaffolding the chromatin modification factors PRC2, LSD1, and DNMT1. Cancer Res. 2016;76(21):6299–310.
Chen Z, Li JL, Lin S, Cao C, Gimbrone NT, Yang R, Fu DA, Carper MB, Haura EB, Schabath MB, et al. cAMP/CREB-regulated LINC00473 marks LKB1-inactivated lung cancer and mediates tumor growth. J Clin Invest. 2016;126(6):2267–79.
Liu YW, Xia R, Lu K, Xie M, Yang F, Sun M, De W, Wang C, Ji G. LincRNAFEZF1-AS1 represses p21 expression to promote gastric cancer proliferation through LSD1-mediated H3K4me2 demethylation. Mol Cancer. 2017;16(1):39.
Qin QH, Yin ZQ, Li Y, Wang BG, Zhang MF. Long intergenic noncoding RNA 01296 aggravates gastric cancer cells progress through miR-122/MMP-9. Biomed Pharmacother. 2018;97:450–7.
Guo X, Hua Y. CCAT1: an oncogenic long noncoding RNA in human cancers. J Cancer Res Clin Oncol. 2017;143(4):555–62.
Zhang E, Han L, Yin D, He X, Hong L, Si X, Qiu M, Xu T, De W, Xu L, et al. H3K27 acetylation activated-long non-coding RNA CCAT1 affects cell proliferation and migration by regulating SPRY4 and HOXB13 expression in esophageal squamous cell carcinoma. Nucleic Acids Res. 2017;45(6):3086–3101
Ozawa T, Matsuyama T, Toiyama Y, Takahashi N, Ishikawa T, Uetake H, Yamada Y, Kusunoki M, Calin G, Goel A. CCAT1 and CCAT2 long noncoding RNAs, located within the 8q.24.21 ‘gene desert’, serve as important prognostic biomarkers in colorectal cancer. Ann Oncol. 2017;28(8):1882–8.
Li J, Han L, Roebuck P, Diao L, Liu L, Yuan Y, Weinstein JN, Liang H. TANRIC: an interactive open platform to explore the function of lncRNAs in cancer. Cancer Res. 2015;75(18):3728–37.
Yang F, Wei J, Zhang S, Zhang X. Shrimp miR-S8 suppresses the stemness of human melanoma stem-like cells by targeting the transcription factor YB-1. Cancer Res. 2017;77(20):5543–53
Goodarzi H, Liu X, Nguyen HC, Zhang S, Fish L, Tavazoie SF. Endogenous tRNA-derived fragments suppress breast cancer progression via YBX1 displacement. Cell. 2015;161(4):790–802.
Kosnopfel C, Sinnberg T, Schittek B. Y-box binding protein 1—a prognostic marker and target in tumour therapy. Eur J Cell Biol. 2014;93(1–2):61–70.
Mihailovich M, Militti C, Gabaldon T, Gebauer F. Eukaryotic cold shock domain proteins: highly versatile regulators of gene expression. Bioessays. 2010;32(2):109–18.
Jurchott K, Kuban RJ, Krech T, Bluthgen N, Stein U, Walther W, Friese C, Kielbasa SM, Ungethum U, Lund P, et al. Identification of Y-box binding protein 1 as a core regulator of MEK/ERK pathway-dependent gene signatures in colorectal cancer cells. PLoS Genet. 2010;6(12):e1001231.
Mo D, Fang H, Niu K, Liu J, Wu M, Li S, Zhu T, Aleskandarany MA, Arora A, Lobo DN, et al. Human helicase RECQL4 drives cisplatin resistance in gastric cancer by activating an AKT-YB1-MDR1 signaling pathway. Cancer Res. 2016;76(10):3057–66.
Diaz-Lagares A, Crujeiras AB, Lopez-Serra P, Soler M, Setien F, Goyal A, Sandoval J, Hashimoto Y, Martinez-Cardus A, Gomez A, et al. Epigenetic inactivation of the p53-induced long noncoding RNA TP53 target 1 in human cancer. Proc Natl Acad Sci U S A. 2016;113(47):E7535–44.
Falkenberg KJ, Johnstone RW. Histone deacetylases and their inhibitors in cancer, neurological diseases and immune disorders. Nat Rev Drug Discov. 2014;13(9):673–91.
Minucci S, Pelicci PG. Histone deacetylase inhibitors and the promise of epigenetic (and more) treatments for cancer. Nat Rev Cancer. 2006;6(1):38–51.
Pearson JC, Lemons D, McGinnis W. Modulating Hox gene functions during animal body patterning. Nat Rev Genet. 2005;6(12):893–904.
Rinn JL, Kertesz M, Wang JK, Squazzo SL, Xu X, Brugmann SA, Goodnough LH, Helms JA, Farnham PJ, Segal E, et al. Functional demarcation of active and silent chromatin domains in human HOX loci by noncoding RNAs. Cell. 2007;129(7):1311–23.
Shah N, Sukumar S. The Hox genes and their roles in oncogenesis. Nat Rev Cancer. 2010;10(5):361–71.
di Pietro M, Lao-Sirieix P, Boyle S, Cassidy A, Castillo D, Saadi A, Eskeland R, Fitzgerald RC. Evidence for a functional role of epigenetically regulated midcluster HOXB genes in the development of Barrett esophagus. Proc Natl Acad Sci U S A. 2012;109(23):9077–82.
Botti G, Marra L, Malzone MG, Anniciello A, Botti C, Franco R, Cantile M. LncRNA HOTAIR as prognostic circulating marker and potential therapeutic target in patients with tumor diseases. Curr Drug Targets. 2017;18(1):27–34.
Zhou Y, Zhang X, Klibanski A. MEG3 noncoding RNA: a tumor suppressor. J Mol Endocrinol. 2012;48(3):R45–53.
Yang F, Huo XS, Yuan SX, Zhang L, Zhou WP, Wang F, Sun SH. Repression of the long noncoding RNA-LET by histone deacetylase 3 contributes to hypoxia-mediated metastasis. Mol Cell. 2013;49(6):1083–96.
Dong J, Akcakanat A, Stivers DN, Zhang J, Kim D, Meric-Bernstam F. RNA-binding specificity of Y-box protein 1. RNA Biol. 2009;6(1):59–64.
Lee KH, Choi EY, Kim MK, Kim KO, Jang BI, Kim SW, Song SK, Kim JR. Inhibition of histone deacetylase activity down-regulates urokinase plasminogen activator and matrix metalloproteinase-9 expression in gastric cancer. Mol Cell Biochem. 2010;343(1–2):163–71.
Peixoto P, Castronovo V, Matheus N, Polese C, Peulen O, Gonzalez A, Boxus M, Verdin E, Thiry M, Dequiedt F, et al. HDAC5 is required for maintenance of pericentric heterochromatin, and controls cell-cycle progression and survival of human cancer cells. Cell Death Differ. 2012;19(7):1239–52.
Cao C, Vasilatos SN, Bhargava R, Fine JL, Oesterreich S, Davidson NE, Huang Y. Functional interaction of histone deacetylase 5 (HDAC5) and lysine-specific demethylase 1 (LSD1) promotes breast cancer progression. Oncogene. 2017;36(1):133–45.
Stypula-Cyrus Y, Damania D, Kunte DP, Cruz MD, Subramanian H, Roy HK, Backman V. HDAC up-regulation in early colon field carcinogenesis is involved in cell tumorigenicity through regulation of chromatin structure. PLoS One. 2013;8(5):e64600.
Feng GW, Dong LD, Shang WJ, Pang XL, Li JF, Liu L, Wang Y. HDAC5 promotes cell proliferation in human hepatocellular carcinoma by up-regulating Six1 expression. Eur Rev Med Pharmacol Sci. 2014;18(6):811–6.
Davies AH, Barrett I, Pambid MR, Hu K, Stratford AL, Freeman S, Berquin IM, Pelech S, Hieter P, Maxwell C, et al. YB-1 evokes susceptibility to cancer through cytokinesis failure, mitotic dysfunction and HER2 amplification. Oncogene. 2011;30(34):3649–60.
De W, Han L, Yin D. A novel long non-coding RNA HOXC-AS3 mediates tumorigenesis of gastric cancer by binding to YBX1. Gene Expression Omnibus. GSE119021. https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE119021. Accessed 27 August 2018.
Hu Y, Xu J, Fang J. LncRNA/mRNA expression profiling for 20 human samples. Gene Expression Omnibus.GSE50710.https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE50710. Accessed 19 Aug 2014
Shiming Y, Lü MH. Human gastric mucosa: 3 paired hTERT-negative para cancerous tissue (Control) vs. hTERT-positive gastric cancer tissue. Gene Expression Omnibus.GSE58828. https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE58828 Accessed 26 June 2014.
This work was supported by National Natural Science Foundation of China (81702266, 81502071, 81401873, and 81772479), China Postdoctoral Science Foundation (2017M610339 and 2017M611913) and Jiangsu Planned Projects for Postdoctoral Research Funds (1701041A). This work was also supported by Key Project and supported by Medical Science and Technology Development Foundation, Nanjing Department of Health (YKK15145), Foundation of Jiangsu Province Medical Youth Talent (QNRC2016057 and QNRC2016380), and Foundation of Xuzhou Central Hospital (XZB201616). This work was supported by the Scientific Foundation of Wuxi City of Jiangsu (Q201728).
Availability of data and materials
Our RNA-seq data used in this study (RNA-seq after knockdown HOXC-AS3 and YBX1) have been deposited in NCBI’s Gene Expression Omnibus and are accessible through GEO accession number GSE119021 . The lncRNA expression profiles data were obtained from GEO, with accession numbers GSE50710  and GSE58828 .
Ethics approval and consent to participate
The study was approved by the Medical Ethical Committee of Southeast University Medical College, and it was performed in compliance with the Helsinki Declaration. All patients have given written informed consent for publication.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1.The clinic-pathological factors of 112 GC patients. (XLS 10 kb)
Table S2. Univariate and multivariate analyses of clinicopathologic factors for overall survival in 112 patients with GC. (XLS 10 kb)
Figure S1. (A) HOXC-AS3 expression after ASO-mediated knockdown and plasmid-mediated overexpression in GC cells. (B) Expression of HOXC-AS3 across diverse normal human tissues from GTEx. Figure S2. (A) Western blots were performed to detect YBX1 expression. (B) The altered mRNA levels of genes were confirmed by qRT-PCR for knockdown HOXC-AS3 in BGC-823 and SGC-7901 cells. (C) Based on qRT-PCR assays, the level of YBX1 was upregulated in 60 pairs GC tissues. MTT assays and transwell assays were used to investigate the changes in proliferation and migratory abilities of BGC-823 cells after transfection. (D) Western blots were performed to detect HDAC5 expression after transfection in BGC-823 cells. (DOC 1991 kb)
Table S3. A list of the top ten potential HOXC-AS3-interacting protein candidates in BGC-823 cells based on RNA-protein pull-down assays and mass spectrometry analysis. (XLS 15 kb)
Table S4. The mRNA variation abundance (≥1.5-fold) for HOXC-AS3-knockdown in BGC-823 cells. (XLS 491 kb)
Table S5. The mRNA variation abundance (≥1.5-fold) for YBX1-knockdown in BGC-823 cells. (XLS 236 kb)
Table S6. The list of primers and siRNA /ASO sequence. (XLS 20 kb)
Supplementary Methods. (DOC 44 kb)