Abstract
Long noncoding RNAs (lncRNA) play a role in carcinogenesis. However, the function of lncRNAs in human gastric cancer remains largely unknown. In this study, we identified a novel lncRNA, GClnc1, which was upregulated and associated with tumorigenesis, tumor size, metastasis, and poor prognosis in gastric cancer. GClnc1 affected gastric cancer cell proliferation, invasiveness, and metastasis in multiple gastric cancer models. Mechanistically, GClnc1 bound WDR5 (a key component of histone methyltransferase complex) and KAT2A histone acetyltransferase, acted as a modular scaffold of WDR5 and KAT2A complexes, coordinated their localization, specified the histone modification pattern on the target genes, including SOD2, and consequently altered gastric cancer cell biology. Thus, GClnc1 is mechanistically, functionally, and clinically oncogenic in gastric cancer. Targeting GClnc1 and its pathway may be meaningful for treating patients with gastric cancer.
Significance: This report documents a novel lncRNA, GClnc1, which may act as a scaffold to recruit the WDR5 and KAT2A complex and modify the transcription of target genes. This study reveals that GClnc1 is an oncogenic lncRNA in human gastric cancer. Cancer Discov; 6(7); 784–801. ©2016 AACR.
This article is highlighted in the In This Issue feature, p. 681
Introduction
Gastric cancer is the fourth most diagnosed type of cancer and the third most common cause of cancer-related death worldwide (1, 2). Patients with advanced gastric cancer have a poor prognosis (3, 4). Pathologic classification is used to assess prognosis and inform the treatment of gastric cancer. Massive efforts have been made to develop the noninvasive biomarkers to detect early cancer and/or reflect an individual's cancer risk, which is essential to reducing gastric cancer mortality (5). However, there has been little success in improving the disease-free survival rate of patients. Because the pathologic mechanisms of gastric cancer progression are not fully understood, more research is needed to discover and develop effective biomarkers and targets for gastric cancer diagnosis and treatment.
Long noncoding RNAs (lncRNA) are a class of noncoding transcripts greater than 200 nucleotides in length. Recent studies have revealed that lncRNAs may affect cancer progression (6, 7). For example, the lncRNAs HOTAIR, MEG3, MALAT1, H19, and GAPLINC may play a role in carcinogenesis (8–12). The role of lncRNAs and the underlying mechanisms in gastric cancer has been reported before (11–16). However, more specific mechanisms of lncRNAs in gastric cancer initiation and development need to be further teased out.
In the current study, we examined genome-wide lncRNA expression profiles in gastric tumors and paired adjacent normal tissue, from which we identified and characterized a novel lncRNA (BC041951) that was associated with prognosis in patients with gastric cancer, and therefore designated it as gastric cancer–associated lncRNA 1 (GClnc1). The biological roles of GClnc1 in gastric carcinogenesis were genetically assessed in several in vitro and in vivo models. Mechanistically, mass spectrometry combined with integrative analysis revealed that GClnc1 not only modified genomic binding of WDR5 and KAT2A, but also specified the pattern of histone modifications on target genes. Furthermore, GClnc1 upregulated the transcription of superoxide dismutase 2 mitochondrial (SOD2), one of the target genes of the WDR5/KAT2A complex, by acting as a scaffold to recruit the WDR5 and KAT2A complex to the SOD2 promoter and increasing H3K4 trimethylation and H3K9 acetylation levels in the SOD2 promoter region. Thus, our study has identified a novel lncRNA, GClnc1, with a biological, mechanistic, and clinical impact on human gastric cancer.
Results
LncRNA Candidate BC041951 Is Clinically Relevant in Gastric Cancer
The Arraystar Human LncRNA Microarray V2.0 (8 × 60K; Arraystar) was used to profile lncRNA expression in 10 paired gastric cancer tissues and paired adjacent normal tissues. In total, 33,045 lncRNAs and 30,215 coding transcripts were collected from the most authoritative databases, such as RefSeq, the UCSC Known Genes dataset (Known Genes 4), Ensembl 37.59, H-invDB 7.0, and RNAdb 2.0, as described in the manufacturer's instructions, and the raw data can be accessed via GSE50710 (12). In order to study these microarray data and rediscover potential biomarkers not yet mined completely, more stringent filtering criteria (raw signal intensity > 1,500, fold change > 2, P < 0.01, Supplementary Fig. S1A) were used in our investigation. We found that eight candidate lncRNAs were significantly increased in gastric cancer tissues compared with adjacent normal tissues (Supplementary Table S1). To validate microarray data, we analyzed the expression of the eight candidate lncRNAs in 165 cases of fresh gastric cancer and adjacent tissues (cohort 1; Supplementary Table S2). Real-time PCR revealed that four of the eight lncRNAs were significantly increased and two of the eight lncRNAs were significantly decreased in cancer versus adjacent tissues of cohort 1 (Supplementary Fig. S1B). These data indicate that a set of lncRNAs is aberrantly expressed in gastric cancer tissues.
We next analyzed the correlation between the four lncRNA candidates and the clinical outcome in cohort 1. The Kaplan–Meier analyses showed that the lncRNA candidates AK094163, ENST00000430239, and NR_024373 have no predictive value for the clinical outcome of patients with gastric cancer, whereas high expression of the lncRNA candidate BC041951 was significantly associated with a poor prognosis in these patients (Fig. 1A; Supplementary Fig. S1C). Furthermore, real-time PCR showed that only BC041951 expression gradually increased from normal gastric tissue to intestinal metaplasia (IM), to dysplasia (DYS), and to gastric cancer (Fig. 1B; Supplementary Fig. S1D). The expression of BC041951 was significantly increased in gastric cancer tissues compared with normal tissues in patients enrolled in Renji Hospital (Supplementary Fig. S1E) and in an external gastric cancer dataset from the Gene Expression Omnibus (GEO) database repository (Supplementary Fig. S1F). In laser capture microdissection (LCM) assays, BC041951 was mainly expressed and was significantly increased in the epithelial cells of gastric cancer tissues compared with that in normal gastric epithelial cells (Supplementary Fig. S1G) and tumor stromal cells (Supplementary Fig. S1H). Thus, an increase in BC041951 may be an early event in the multistep progression of gastric carcinogenesis.
To evaluate the pathologic and clinical value of the four overexpressed lncRNA candidates with the American Joint Committee on Cancer (AJCC) stage for prognosis, we computed their accuracy by prediction error curves. Predication error over time was calculated using the Brier Score. The predication error was the lowest using the combination of BC041951 and AJCC stage in cohort 1 (Fig. 1C), but not with the other three lncRNA candidates (Supplementary Fig. S1I). Further, ROC analysis showed that the AUC of the combination of BC041951-based prediction and AJCC-based model (0.732) was higher than the AJCC-based model alone (0.695; Supplementary Fig. S1J). The data indicate that the combination of BC041951 and AJCC stage is more precise in predicting clinical outcome than AJCC stage alone. Therefore, we focused our research on BC041951, and henceforth named this lncRNA candidate gastric cancer–associated lncRNA1 (GClnc1).
We next evaluated and compared GClnc1 expression with different clinicopathologic features in cohort 1. We found that GClnc1 expression positively correlated with pathologic differentiation, vascular invasion, tumor size, and AJCC stage (Fig. 1D). Univariate and multivariate regression analyses of cohort 1 demonstrated that GClnc1 expression was an independent predictor of gastric cancer aggressiveness with significant HRs for predicting clinical outcome. Its predictive value was comparable to that of the AJCC stage (Supplementary Fig. S1K and S1L; Fig. 1E).
To further validate the pathologic and clinical significance of GClnc1 expression in gastric cancer, we detected and compared GClnc1 expression by in situ hybridization (ISH) in an additional 105 paraffin-embedded gastric cancer and adjacent tissues (cohort 2; Supplementary Table S3). GClnc1 expression was higher in gastric cancer tissues than in adjacent tissues (Fig. 1F). Consistent with the results in cohort 1, high levels of GClnc1 expression were significantly associated with poor survival (Fig. 1G) in univariate and multivariate regression analyses (Supplementary Fig. S1M and S1N; Fig. 1H). Because H. pylori infections may be detected in 78% of all cases of patients with gastric cancer (17), we measured and found similar levels of GClnc1 expression in H. pylori–positive and H. pylori–negative tumor tissues (Supplementary Fig. S1O).
In order to determine whether GClnc1 is a novel lncRNA, we first analyzed its sequence. Northern blot revealed that the size of GClnc1 was approximately 2 kb in length in gastric cancer cell lines and gastric cancer tissues (Supplementary Fig. S1P and S1Q). The 5′ rapid amplification of cDNA ends (RACE) PCR (Supplementary Fig. S1R) and the chromatin immunoprecipitation (ChIP) assay with anti-RNA Pol II (Supplementary Fig. S1S) were performed to identify the 5′ ends and the transcription start site of GClnc1. Sequencing of PCR products revealed the boundary between the universal anchor primer and GClnc1 (Supplementary Fig. S1R). ChIP assays demonstrated that RNA Pol II bound to the predicted bind site of the 5′ ends of GClnc1 (Supplementary Fig. S1S). Thus, GClnc1 is an RNA Pol II transcript. Analysis of the GClnc1 sequences by ORF Finder from the National Center for Biotechnology Information failed to predict a protein (Supplementary Fig. S1T). GClnc1 did not contain a valid Kozak sequence. It was identified to be an lncRNA rather than a protein-coding transcript by CNCI software (18). We performed a codon substitution frequency analysis using PhyloCSF (19). GClnc1 had a very low codon substitution frequency score (–838.1; Supplementary Fig. S1U), indicating that it was unlikely to encode any protein product. We confirmed that GClnc1 was an ncRNA by in vitro translation analysis (Supplementary Fig. S1V). We also separated the nuclear and cytoplasm fractions of BGC823 and MKN45 gastric cancer cells and performed real-time PCR. We found that GClnc1 was mainly located in the nucleus of BGC823 and MKN45 cancer cells (Supplementary Fig. S1W). Collectively, GClnc1 is a novel lncRNA and highly expressed in gastric cancer tissues.
GClnc1 Is an Oncogenic LncRNA in Gastric Cancer
To elucidate whether GClnc1 plays a role in gastric cancer tumorigenesis, an RNA sequencing (RNA-seq) analysis was performed to compare the gene expression profiles of GClnc1 shRNA and control shRNA transfectants. A total of 2,050 downregulated genes and 1,894 upregulated genes (≥2-fold) were detected (raw data accessible via GEO number GSE63765) after knockdown of GClnc1 in gastric cancer cells (Supplementary Table S4). Gene set enrichment analysis (GSEA) revealed that the gene sets related to Benporath_Proliferation (cell proliferation), Provenzani_Metastasis_Up (metastasis), Kim_Gastric_Cancer_Chemosensitivity (upregulated genes in chemoresistant gastric cancer cells), and Vecchi_Gastric_Cancer_Early_Up (gastric cancer–specific signature) negatively correlated with GClnc1 downregulation in gastric cancer cells (Fig. 2A–D). The top-scoring genes recurring in the four gene sets included key cancer genes, c-MYC and CDK1. Real-time PCR confirmed that alteration of GClnc1 expression dramatically affected the key tumorigenesis gene signatures (Supplementary Fig. S2A). To gain further insight into the biological pathways involved in gastric cancer pathogenesis, based on the median of GClnc1 expression levels, we performed GSEA analysis in an independent public dataset from GEO (GSE27342; Supplementary Fig. S2B). GSEA showed that the gene signatures of cell proliferation, cell cycle, cancer-related pathways, metastasis, and chemotherapy resistance pathways were enriched in patients with high GClnc1 expression, but not in patients with low GClnc1 expression (Fig. 2E). These data suggest that GClnc1 may be an important modulator in gastric tumorigenesis.
To functionally validate the pathway findings, we transfected GClnc1 siRNAs into the gastric cancer cell lines BGC823 and MKN45. The two cell lines expressed higher levels of GClnc1 (Supplementary Fig. S3A) and exhibited higher proliferation (Supplementary Fig. S3B) and invasion capacities (Supplementary Fig. S3C) compared with normal GES-1 gastric cells (Supplementary Fig. S3A–S3C). The results showed that knockdown of GClnc1, but not the lncRNA ENST00000430239 (ENST430239), significantly impaired gastric cancer cell proliferation and colony formation in BGC823 cells and MKN45 cells (Fig. 3A; Supplementary Fig. S3D–S3F). Knockdown of GClnc1 dramatically reduced BGC823 and MKN45 tumor growth (Fig. 3B and C; Supplementary Fig. S3G and S3H) and tumor weight (Fig. 3D; Supplementary Fig. S3I) in xenograft mouse tumor models. In support of the protumor role of GClnc1, Ki67 staining revealed that downregulation of GClnc1 decreased tumor cell proliferation in vivo (Supplementary Fig. S3J and S3K). Furthermore, knockdown of GClnc1 sensitized BGC823 and MKN45 cells to 5′-fluorouracil and cisplatin treatment (Fig. 3E and F; Supplementary Fig. S3L and S3M). The data suggest that GClnc1 may be an oncogenic lncRNA in gastric cancer and control gastric cancer cell proliferation and chemoresistance.
Next, we examined the effects of GClnc1 on gastric cancer cell invasion and metastasis. In the invasion assay, we showed that downregulation of GClnc1, but not the lncRNA ENST430239, significantly reduced the invasion ability in gastric cancer cells (Fig. 3G; Supplementary Fig. S3N). In a gastric cancer metastatic model, mice inoculated with GClnc1 shRNA-expressing tumor cells had a longer overall survival time than mice that received control shRNA–expressing tumor cells or phosphate-buffered solutions (PBS; Fig. 3H; Supplementary Fig. S3O). There were fewer metastatic foci in the lungs of nude mice at 13 weeks after injection of GClnc1 shRNA adenovirus, when compared with control groups (Fig. 3I; Supplementary Fig. S3P). Overexpression of GClnc1 increased cell proliferation and invasion ability of GES-1 and MGC803 cells in vitro and in vivo (Supplementary Fig. S3Q–S3V). The data strongly suggest that GClnc1 may promote gastric cancer progression by regulating gastric cancer cell proliferation and metastasis, and increasing chemotherapy resistance.
GClnc1 Interacts with WDR5 and KAT2A Epigenetic Modification Complex
LncRNAs may function by physically interacting with transcriptional factors, histone regulators, and other cellular factors (20–23). To explore the mechanism of GClnc1-mediated carcinogenesis in gastric cancer, we sought to identify intracellular GClnc1-binding factors using an unbiased approach. Biotinylated GClnc1 or antisense GClnc1 RNA (negative control) was incubated with total protein extracts from BGC823 cells and pulled down with streptavidin (Fig. 4A). The associated proteins were analyzed by SDS-PAGE and silver staining. There were two specific bands in the GClnc1 pull-down samples (Fig. 4B). We excised and analyzed the two bands by mass spectrometry and identified 94 and 101 potential binding proteins in bands one and two, respectively (Supplementary Table S5). Among these binding proteins, we selected 10 proteins that may likely participate in transcription regulation (24–32) for further binding validation. The 10 proteins included eukaryotic translation initiation factor 3 subunit I (EIF3I), Zinc finger RAN-binding domain-containing protein 2 (ZRANB2), eukaryotic translation initiation factor 3 subunit G (EIF3G), WD repeat-containing protein 5 (WDR5), methyl-CpG-binding domain protein 3 (MBD3), AP-1 complex subunit gamma-1 (AP1G1), nucleolar transcription factor 1 (UBTF), histone acetyltransferase (HAT) KAT2A (KAT2A), eukaryotic translation initiation factor 3 subunit B (EIF3B), and transcription initiation factor TFIID subunit 5 (TAF5) proteins. Western blot analysis showed that only WDR5 and KAT2A bound specifically to GClnc1 (Fig. 4C; Supplementary Fig. S4A). The data suggest that the two proteins may form a complex to interact with GClnc1. WDR5 is a core subunit of the human mixed-lineage leukemia (MLL) and SET domain containing 1A (SET1, hCOMPASS) histone H3 Lys4 (H3K4) methyltransferase complexes, and an “effector” of H3K4 methylation in gene transactivation (33, 34). KAT2A is a HAT that functions primarily as a transcriptional activator, via modulating histone H3 Lys9 (H3K9) acetylation levels in the gene promoter region (35). To identify the WDR5/KAT2A-interacting region of GClnc1, we constructed and biotinylated four fragments of GClnc1 (1–500Δ1, 501–1,000Δ2, 1,001–1,500Δ3, 1,501–2,155Δ4) and used them in the pull-down assay with BGC823 cell lysates. We found that the 5′ fragment of GClnc1 mediated the interaction with WDR5 and KAT2A (Fig. 4D). To substantiate the observation, anti-WDR5 and anti-KAT2A antibodies were used to immunoprecipitate endogenous WDR5 and KAT2A from nuclear extracts of BGC823 cells. To this end, RNAs bound to WDR5 and KAT2A were extracted and analyzed. PCR data revealed that GClnc1 directly bound with WDR5 and KAT2A in gastric cancer cells (Fig. 4E). We also detected approximately 5.5-fold and 4-fold enrichments of GClnc1, but not GAPDH RNA, in the anti-WDR5 and anti-KAT2A immunoprecipitates, respectively, compared with the IgG control (Fig. 4F and G). Thus, GClnc1 may specifically bind with WDR5 and KAT2A in gastric cancer cells. Furthermore, we observed that WDR5 and KAT2A interacted with each other in gastric cancer cells (Fig. 4H). Downregulation of GClnc1 had no effect on the expression of WDR5 and KAT2A (Supplementary Fig. S4B). However, knockdown of GClnc1 or RNAase treatment abolished the interaction between WDR5 and KAT2A (Fig. 4I). Thus, GClnc1 is important for WDR5 and KAT2A complex interaction.
GClnc1 Coordinates the Localization of WDR5 and KAT2A Genome-Wide
To address whether GClnc1 modulates WDR5 and KAT2A genomic binding genome-wide, we performed chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-seq) for WDR5 and KAT2A in BGC823 cells. WDR5 and KAT2A ChIP-seq data revealed 5,549 and 6,485 called peaks, respectively, in BGC823 cells transfected with control shRNA. GClnc1 shRNA virus transfection caused reduced WDR5 and KAT2A occupancies in these two histone modulator binding DNA regions (Fig. 5A and B). Further analysis showed that WDR5 and KAT2A occupied 966 and 734 gene promoters, respectively (Fig. 5C). Nearly 20% of KAT2A-occupied genes (147 genes) were also occupied by WDR5, revealing a significant overlap (P < 0.01, hypergeometric distribution). Knockdown of GClnc1 led to a concordant decrease in WDR5 and KAT2A occupancies in 80% of the WDR5 and KAT2A co-occupied gene promoter regions (Supplementary Fig. S5A and S5B; Supplementary Table S6). The promoter regions of these genes correspondingly decreased the levels of H3K4me3 and H3K9Ac, the respective histone methylation and acetylation products of the WDR5 and KAT2A complex (Supplementary Fig. S5C). To further verify these findings, we performed ChIP for WDR5, KAT2A, H3K4me3, and H3K9Ac in BGC823 and MKN45 cells transfected with control shRNA or GClnc1 shRNA virus. We randomly selected 10% of the overlapped genes (147 genes) to experimentally validate the effect of GClnc1-mediated WDR5/KAT2A complex genomic binding. We found that knockdown of GClnc1 resulted in a decrease in the binding efficiency of WDR5 and KAT2A and the levels of H3K4me3/H3K9Ac in approximately 80% of the selected binding genes in both BGC823 and MKN45 cells (Supplementary Fig. S5D and S5E). Moreover, overexpression of GClnc1 dramatically increased the binding efficiency of WDR5 and KAT2A and the levels of H3K4me3/H3K9Ac of the same measured genes in normal GES-1 cells (Supplementary Fig. S5F).
Next, we used the RNA-seq data from BGC823 cells transfected with control or GClnc1 shRNA virus (GSE63765) to characterize the relationship between WDR5/KAT2A complex binding and GClnc1-mediated changes in gene expression. After identifying a gene signature with highly significant changes in expression (Supplementary Table S6D), we intersected this signature with the ChIP-seq data. We observed that a substantial subset of genes which exhibit decreased binding by the WDR5/KATA2 complex were dysregulated after knockdown of GClnc1 (Fig. 5D and Supplementary Table S6E). Integrative GSEA of the RNA-seq and WDR5/KAT2A ChIP-seq data demonstrated significant enrichment for genes that were repressed when GClnc1 was downregulated (Fig. 5E). Thus, knockdown of GClnc1 impairs formation of the WDR5 and KAT2A complex. This suggests that GClnc1 may act as a scaffold for this complex, thereby controlling its ability to regulate histone modification, gene expression, and biological activity.
GClnc1 Promotes Gastric Cancer Progression via SOD2
According to the ChIP-seq and RNA-seq data, SOD2 is one of the target genes of the WDR5/KAT2A complex (Supplementary Table S6E). Furthermore, we found the strongest association between SOD2 and GClnc1 expressions (Supplementary Table S7 and Fig. 6A) among 12 genes. The positive correlation between SOD2 and GClnc1 was also validated in the 165 cases of gastric cancer in our cohort 1 (Fig. 6B; Supplementary Fig. S6A).
GClnc1 is located in chromosome 6 (chr6:159679116-159681261) with a 2155 nt length. A 195 nt–length DNA fragment of GClnc1 at the 3′ terminal overlaps with the last exon of the two transcript variants of the SOD2 gene (NM_001024465, NM_001024466). The other part of GClnc1 is completely overlapped with the last intron of the SOD2 gene (Fig. 6C). We next designed the sense-specific PCR primers (Supplementary Fig. S6B) to determine whether GClnc1 and SOD2 were different transcripts. If GClnc1 and SOD2 were different transcripts, the PCR products of primer 1 + primer 2 should be approximately 150 bp (part of the GClnc1 fragment) and the PCR products of primer 4 + primer 2 should be approximately 200 bp (part of the SOD2 fragment). In addition, there should be no PCR product to be detected by primer 3 + primer 4. If GClnc1 and SOD2 were the same transcripts, in addition to the 150 bp and 200 bp PCR products, there should be an additional PCR product (around 3,260 bp) to be detected by primer 3 and primer 4. Our results showed that the 150 bp and 200 bp PCR products were detectable; however, the 3,260 bp PCR product was not detectable (Supplementary Fig. S6C). Furthermore, Northern blot analysis confirmed the sense-specific RT-PCR data (Supplementary Fig. S6D). Thus, GClnc1 and SOD2 are different transcripts. MiRNA is not found in the GClnc1 coding region. We next explored the mechanism by which GClnc1 regulates the expression of SOD2. Real-time PCR and Western blot data showed that SOD2 expression was significantly decreased in BGC823 cells transfected with two different GClnc1 siRNAs, but not with lncRNA ENST430239 (Fig. 6D). Similar results were observed in MKN45 cells (Supplementary Fig. S6E). Real-time PCR and Western blot data (Supplementary Fig. S6F) showed that SOD2 expression was significantly increased in GClnc1-overexpressing GES-1 cells. We next examined the effect of GClnc1 on the transcriptional activity of the SOD2 gene. Luciferase assays revealed that knockdown of GClnc1, not lncRNA ENST430239, impaired the transcriptional level of the SOD2 promoter in BGC823 cells (Fig. 6E) and in MKN45 cells (Supplementary Fig. S6G). The data suggest that GClnc1 may positively regulate SOD2 transcription in gastric cancer cells. It has been reported that high levels of SOD2 expression correlated with poor clinical outcome in gastric cancer (36), and SOD2 may regulate cancer cell proliferation and invasion (37, 38). In addition, gastric cancer cell proliferation and invasion ability were compared after knockdown of GClnc1 and SOD2. We found that downregulation of SOD2 and GClnc1 similarly impaired cell proliferation and invasion (Supplementary Fig. S6H).
We next hypothesized that SOD2 mediated the biological function of GClnc1 in gastric cancer. To test this hypothesis, we transfected SOD2 siRNA into gastric cancer cells and examined its effects on cancer cell biological function. Knockdown of SOD2 expression significantly reduced gastric cancer cell proliferation induced by GClnc1 in BGC823 cells (Fig. 6F) and MKN45 cells (Supplementary Fig. S6I). We next used the CRISPR/Cas9 system to knock out the SOD2 gene in BGC823 and MKN45 cells. We obtained similar results in these cells (Supplementary Fig. S6J). Furthermore, inhibition of SOD2 expression significantly blocked tumor growth in GClnc1-overexpressing xenograft mice (Supplementary Fig. S6K–S6M). Downregulation of SOD2 also markedly blocked the GClnc1-induced cell invasion in BGC823 cells (Fig. 6G; Supplementary Fig. S6N) and MKN45 cells (Supplementary Fig. S6O and S6P). Thus, SOD2 mediates the regulatory function of GClnc1 in gastric cancer cells.
GClnc1 Is a Molecular Link between WDR5/KAT2A and SOD2
H3K4 methylation and H3K9 acetylation may regulate the transcription of the SOD2 gene promoter (39). Real-time PCR and Western blot showed that knockdown of WDR5 or KAT2A significantly blocked GClnc1-induced SOD2 upregulation in BGC823 cells (Fig. 7A; Supplementary Fig. S7A). The data suggest that WDR5, KAT2A, and GClnc1 may cooperate to regulate SOD2 gene transactivation via histone modification. The ChIP assay indicated that WDR5 and KAT2A directly bound to the promoter region of SOD2 (Fig. 7B), and knockdown of WDR5 or KAT2A significantly decreased the binding efficiency of WDR5 or KAT2A, respectively, whereas knockdown of GClnc1 decreased the binding efficiencies of both WDR5 and KAT2A (Fig. 7C). These data indicate that GClnc1 may play an important role in the binding progression of WDR5 or KAT2A to the SOD2 gene promoter region. Furthermore, we observed that H3K4 trimethylation and H3K9 acetylation existed in the same promoter region of SOD2 (Fig. 7D). Knockdown of GClnc1 significantly decreased the H3K9 acetylation and H3K4 methylation levels in the SOD2 promoter region, but not the H3K27 trimethylation level (Fig. 7E). Overexpression of GClnc1 dramatically increased the H3K9 acetylation and H3K4 methylation levels at the promoter region of the SOD2 gene (Supplementary Fig. S7B). The chromatin isolation by RNA purification (ChIRP) assay revealed that GClnc1 directly bound to the SOD2 promoter region (Fig. 7F). The combined analyses strongly suggest that GClnc1 may act as a scaffold lncRNA to recruit WDR5 and KAT2A to the promoter of the SOD2 gene and transactivate the transcription of the SOD2 gene.
Discussion
Multiple oncogenic pathways may contribute to gastric cancer carcinogenesis (40, 41); however, the potential involvement of lncRNAs is poorly defined in human gastric cancer. Through a combination of genomic, biochemical, and cell-biological analyses, we have demonstrated that GClnc1 is a novel oncogenic lncRNA in gastric cancer. GSEA analyses have demonstrated that cell proliferation, metastasis, and chemoresistance pathways in cancer are significantly enriched in response to GClnc1 alteration in patients with gastric cancer. The bioinformatic analyses have been functionally validated in several in vitro and in vivo experimental models. Downregulation of GClnc1 markedly suppresses cell growth and metastasis and increases gastric cancer cell chemosensitivity in cultured gastric cancer cells and xenograft mouse models. Thus, the data consistently point to the notion that high GClnc1 expression is a decisive factor in controlling human gastric cancer aggressiveness.
LncRNAs may participate in chromatin remodeling complexes (42–46); however, the underlying molecular mechanisms remain unknown. The mass spectrometry and RNA pull-down data have demonstrated that GClnc1 directly interacts with WDR5, the key component of histone methyltransferase complex, and the histone acetyltransferase KAT2A. Furthermore, GClnc1 can link WDR5 and KAT2A by acting as a potential modular scaffold in gastric cancer cells. This notion is supported by three lines of experimental evidence: (i) GClnc1 directly binds with WDR5 and KAT2A via its 5′ domain; (ii) genetic deficiency of GClnc1 or RNase treatment of the immunoprecipitation products abrogates the interaction between WDR5 and KAT2A; and (iii) knockdown of GClnc1 significantly disrupts the binding of the WDR5/KAT2A complex to the promoter regions of WDR5/KAT2A co-occupied genes. In support of our observation, the lncRNAs HOTAIR and HOTTIP assist epigenetic complexes, including PRC2 and MLL, facilitate their genomic binding, and enhance their functions (20, 21, 46). Similarly, GClnc1 modulates the interaction of the WDR5 and KAT2A complex to further alter the pattern of histone modification and transactivate target genes in gastric cancer cells (Fig. 7G).
When we explored the mechanisms by which GClnc1 contributes to gastric carcinogenesis, we found the involvement of SOD2. It has been reported that SOD2 expression is increased in gastric cancers (36). The expression of SOD2 can be epigenetically silenced in human breast cancer cells, and H3K4 methylation and H3K9 acetylation may be the major histone regulatory pattern in SOD2 gene transcription (39). Consistent with these studies, we show that GClnc1 activates SOD2 transcription and that SOD2 is functionally responsible for GClnc1-mediated gastric carcinogenesis. Notably, as GClnc1 may be directly and indirectly linked to gene regulatory networks, in addition to SOD2, we do not rule out the involvement of other genes in GClnc1-associated biologic function.
Given the clinical, genetic, biochemical, and functional significance of GClnc1 in gastric cancer, we conclude that GClnc1 and its associated pathway is crucial for gastric carcinogenesis, and targeting this pathway may be pivotal in the prevention or treatment of gastric cancer.
Methods
Patient Specimens
Tumors and the adjacent gastric tissues were obtained from patients with gastric cancer who underwent surgery at Shanghai Renji Hospital from January 2000 to January 2005. The study protocol was approved by the ethics committee of Shanghai Jiao Tong University School of Medicine, Renji Hospital. Written informed consent was obtained from all participants in this study. All the research was carried out in accordance with the provisions of the Declaration of Helsinki of 1975. None of these patients had received radiotherapy or chemotherapy prior to surgery. The percentage of tumor cellularity in the tissue section was at least 70% as determined by a pathologist.
Human gastric mucosal tissues (normal tissues, tissues diagnosed with IM or DYS) were collected from patients by gastroscope inspection in Renji Hospital with written informed consent. None of the patients had taken nonsteroidal anti-inflammatory drugs, H2 receptor antagonists, proton pump inhibitors, antimicrobials, or bismuth compounds 4 weeks prior to the study. We used 20 samples per group to ensure adequate statistical power to distinguish phenotypes between the groups. The extent of inflammation in these tissues was examined according to the updated Sydney System (International Workshop on the Histopathology of Gastritis, Houston, 1994).
Bioinformatics Analysis
LncRNAs were identified from the gene expression data which can be accessed by GEO Series accession number GSE50710. The RNA-seq and ChIP-seq data can be accessed by GEO Series accession number GSE63765. The detailed Bioinformatics analyses are described in the Supplementary Materials and Methods.
Cell Culture and Treatment
The human gastric epithelial cell line GES-1 and the gastric cancer cell lines MKN45 and AGS were obtained from Dr. Jun Yu at the Chinese University of Hong Kong. The human gastric cancer cell lines BGC823, SGC7901, and MGC803 were obtained from Dr. Shuiping Tu of Rutgers, the State University of New Jersey. All cell lines were genotyped for identity by Beijing Microread Genetics Co., Ltd and tested routinely for Mycoplasma contamination. Cell lines were cultured in RPMI-1640 medium (Gibco) supplemented with 10% FBS at 37°C in an atmosphere of 5% CO2. Cell lines were transfected with different siRNA and GClnc1-expressing plasmids for genetic and functional assays. The detailed cell line treatments are described in the Supplementary Materials and Methods.
Sense-Specific RT-PCR
The sense-specific RT-PCR (Supplementary Fig. S6B) assay was developed after pilot and optimizing experiments. The detailed sense-specific RT-PCR is described in the Supplementary Materials and Methods.
5′ Rapid Amplification of cDNA Ends
5′-RACE was performed to determine the transcriptional initiation site of GClnc1 using a SMARTer RACE cDNA Amplification Kit (Clontech). The gene-specific primer used for the PCR of the RACE analysis was 5′-TATGTACACAGTGTTGCTCACC-3′ (5′RACE).
Cell Proliferation and Transwell Assay
Cell proliferation was assessed by the BrdUrd incorporation assay (Roche Molecular Biochemicals). Briefly, control and treated gastric cancer cells were seeded onto 96-well plates at an initial density of 5 × 103 cells per well. BrdUrd labeling solution (10 μL/well) was added to the cells at specified time points. After incubating for 2 hours, the reaction product was quantified according to the manufacturer's instructions.
Transwell assays were performed using 24-well transwells (8-μm pore size; Millipore) precoated with Matrigel (BD Biosciences). Cells at logarithmic phase were transfected with GClnc1 siRNA or control siRNA. Transfected cells were then harvested, and 1 × 105 cells were seeded in serum-free medium into the upper chamber, whereas medium supplemented with 20% FBS was applied to the lower chamber as a chemoattractant. After 48 hours of incubation, the migrated cells at the bottom surface of the filter were fixed, stained, and counted.
Hybridization In Situ and Immunohistochemical Staining
The in situ detection of GClnc1 was performed on 6-μm formalin-fixed, paraffin-embedded sections using DIG-labeled miRCURYTM Detection probe (Exiqon). Positive controls (RNU6B; Exiqon) and scrambled control RNAs were included for each hybridization procedure and analyzed using a Nikon 80i microscope with Nikon NIS-Elements F 2.3 software (Nikon).
The expression of Ki67 was examined with an anti-Ki67 antibody (Ki67; Abcam; 1:100) using the LSAB+ Kit (DakoCytomation). Slides were examined independently by two investigators. Protein expression was quantified using a visual grading system based on the extent of staining (the percentage of positive tumor cells on a scale of 0–4: 0, none; 1, 1%–25%; 2, 26%–50%; 3, 51%–75%; 4, >75%) and the intensity of staining (graded on a scale of 0–3: 0, no staining; 1, weak staining; 2, moderate staining; 3, strong staining). The product of the extent and intensity grades was used to define the Ki67 protein expression.
Luciferase Assay
Cell lines were transfected with the plasmids expressing the designated combinations of pGL3-SOD2PWT and other relevant siRNAs at 1.0 μg and 100 ng of phRL (Renilla luciferase) with Lipofectamine 2000 (Invitrogen) or DharmaFECT 1 siRNA transfection reagents (Dharmacon). Twenty-four hours after transfection, the cells were collected to detect luciferase activity using the Dual-Luciferase reporter assay system (Promega). Luciferase activity was measured by using a BD Monolight 3010 luminometer (BD Biosciences). Transfection efficiency was normalized by dividing the luciferase activity of the construct by the corresponding Renilla luciferase activity.
RNA Analysis, Extraction, and Quantitative Real-Time PCR
The RNA expression levels were measured using a real-time quantitative PCR system. Total RNA was extracted by TRIzol reagent (Invitrogen), and 1 μg of total RNA was reverse-transcribed using the PrimeScriptP RT Reagent Kit (Perfect Real Time; Takara). The amplified transcript level of each specific gene was normalized to that of 18S of GAPDH. The primers (Supplementary Table S8) were provided by Shenggong Company.
Western Blot
Western blot analysis was performed using standard techniques (47). An anti-GAPDH antibody was used as a control for whole-cell lysates. Antibodies were purchased mainly from three manufacturers: Cell Signaling Technology, Inc. (KAT2A, WDR5, GAPDH), Abcam (SOD2, ZRANB2, EIF3B, TAF5), and Sigma (EIF3I, EIF3G, MBD3, AP1G1, UBTF). The dilution of primary antibodies was 1:1,000.
Northern Blot
Northern blot analysis was performed with 300 to 600 ng of purified poly (A) mRNA. RNAs were resolved by denaturing agarose gel electrophoresis (Ambion) and transferred to Hybond-XL membranes (GE Healthcare). LncRNAs were detected using 32P-labeled DNA probes.
In Vitro Transcription and Translation
LncRNA GClnc1 and LacZ were cloned into pBluescript KSII downstream of the T7 promoter. The plasmid was linearized by single digestion with KpnI. Linear pCDNA3.1-GClnc1 (250 ng) was transcribed in vitro (Epicentre), and 700 ng of purified GClnc1 RNAs were translated in vitro using biotinylated leucine tRNA (Promega). Biotinylated proteins were detected using a BrightStar BioDetect Kit (Ambion). LacZ mRNA was used as a positive translation control, and mock-translated samples (no RNA template) were used as negative controls.
RNA Pull-Down Assay
Cell nuclear lysates were incubated with biotinylated RNA and streptavidin beads for RNA pull-down incubation. Then beads were collected by centrifugation. RNA-associated proteins were eluted and resolved by SDS-PAGE followed by silver staining (Bio-Rad). The detailed RNA pull-down assay is described in the Supplementary Materials and Methods.
LC-MS Analysis, Database Search, and Protein Identification
To identify specific GClnc1 interactors, GClnc1 and Antisense GClnc1 pull-down eluates were compared. The bands that were predominantly represented in the GClnc1 pull-down sample were chosen. The bands were excised to perform LC-MS analysis. The proteins sequenced with at least two peptides were considered a reliable identification. The detailed LC-MS, database, and protein identification are described in the Supplementary Materials and Methods.
RNA Immunoprecipitation
RNA immunoprecipitation (RIP) experiments were performed using the Megna RIP RNA-binding Protein Immunoprecipitation Kit (Millipore). The anti-WDR5 antibody and anti-KAT2A antibody used for RIP were purchased from CST (Cell Signaling Technology). The co-precipitated RNAs were detected by reverse transcription PCR. To demonstrate that the detected RNA signals specifically bind to WDR5 or KAT2A, total RNA (input controls) and normal rabbit IgG controls were simultaneously assayed. The gene-specific primers used for detecting GClnc1 are displayed in Supplementary Table S8.
Coimmunoprecipitation
Coimmunoprecipitation was performed as described previously (48). Briefly, both input and IP samples were analyzed by Western blot using various antibodies at the indicated dilutions: anti-WDR5 antibody (1:1,000; Cell Signal Technology), anti-KAT2A antibody (1:1,000; Cell Signal Technology), and normal rabbit IgG.
ChIRP
GClnc1 asDNA (antisense DNA), LacZ asDNA (antisense DNA), and GClnc1 sDNA (sense DNA) probes were designed using the online probe designer at singlemoleculefish.com. Oligonucleotides were biotinylated at the 3′ end with an 18-carbon spacer arm. BGC823 cells were collected and subjected to ChIRP as previously described (49).
In Vivo Experiments
In order to clarify the effect of GClnc1 on tumor growth and invasion in vivo, 4-week-old male BALB/c nude mice (Experimental Animal Centre of Shanghai Institutes for Biological Sciences) were used. Different numbers of cells were inoculated into the mice. Ten days after subcutaneous inoculation, mice were treated with a different virus. Tumor growth and invasion were monitored. All experimental procedures were approved by the Institutional Animal Care and Use Committee of Ren Ji Hospital, School of Medicine, Shanghai Jiao Tong University. The detailed in vivo experiments are described in the Supplementary Materials and Methods.
Adenovirus and Plasmids Construction
The control shRNA, GClnc1 shRNA1/2, SOD2 shRNA, GClnc1-overexpressing adenovirus, and all plasmids were constructed by Shanghai Obio Technology Company.
ChIP and High-Throughput Sequencing
Chromatin immunoprecipitation and high-throughput sequencing (ChIP-seq) was performed as follows. Chromatin was sonicated and immunoprecipitated with different antibodies. Real-time PCR was performed in triplicate. Each PCR assay was carried out in a 20 μL reaction volume using the eluted immunoprecipitated DNA. The amount of genomic DNA coprecipitated with the specific antibody was calculated as follows: CBTB = CBTB (genomic input) − CBTB (specific antibody), where CBTB (genomic input) and CBTB (specific antibody) were the mean threshold cycles of PCR performed in triplicate on DNA samples from the genomic input samples and the specific antibody samples, respectively. The forward and reverse primers are listed in Supplementary Table S9. High-throughput RNA and ChIP-seq was performed after knockdown of GClnc1 in BGC823 cells. The detailed ChIP and High-Throughput Sequencing are described in the Supplementary Materials and Methods.
Drug Sensitivity Assays to 5′-Fluorouracil and Cisplatin
Cell proliferation was assessed by the BrdUrd incorporation assay (Roche Molecular Biochemicals). Briefly, control and siRNA-transfected cells, which were seeded onto 96-well plates at an initial density of 5 × 103 cells per well, were treated with different doses of 5′-fluorouracil and cisplatin, and BrdUrd labeling solution (10 μL/well) was added to the cells at specified time points. After 2-hour incubation, culture medium was removed and the cells were fixed. Then DNA was denatured by adding FixDenat (200 μL/well), and anti–BrdUrd-POD working solution (100 μL/well) was added to the cells and incubated for 90 minutes. The immune complexes were detected by the subsequent substrate reaction. The reaction product was quantified by measuring the absorbance at 370 nm (reference wavelength: approximately 492 nm).
Tissues and LCM
Gastric cancer tissues and adjacent normal tissues were obtained from 20 patients enrolled in Renji Hospital. Tissues were embedded in Tissue-Tek OCT medium (Sakura) and snap-frozen on dry ice. The frozen sections were fixed in 70% ethanol for 30 seconds and stained with hematoxylin and eosin, followed by three dehydration steps of 5 seconds each in 70%, 95%, and 99.5% ethanol and a final 5-minute dehydration in xylene. Once air-dried, the stained tissues were laser-capture microdissected by a PixCell II LCM system (Arcturus Engineering). Normal gastric epithelial cells, gastric cancer cells, and stroma cells were selectively captured by microscopic visualization.
Generation of SOD2 Knockout BGC823 (BGC823KO) and (MKN45KO) Cell Lines
LentiCRISPR (Addgene plasmid #49535) is a lentiviral vector expressing Cas9 and single-guide RNA (sgRNA; ref. 50). Genome engineering experiments using CRISPR/Cas9 systems were performed as described previously (51). Briefly, three sgRNAs targeting SOD2 were cloned into the lentiCRISPR vector. Lentivirus was produced in HEK 293T cells. BGC823 and MKN45 cells were transduced with lentivirus and selected with 2 μg/mL puromycin. The SOD2 knockout cell clones were isolated by limited serial dilution. SOD2 deletion was confirmed by sequencing and immunoblotting analysis. The sgRNA sequences targeting SOD2 were as follows: 5′-TTCCAGGGCGCCGTAGTCGT-3′ (site 1); 5′-AAGCTGACGGCTGCATCTGT-3′ (site 2); and 5′-CTTGGTTTCAATAAGGAACG-3′ (site 3; ref. 50).
Statistical Analysis
The detailed statistical analysis is described in the Supplementary Materials and Methods.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Authors' Contributions
Conception and design: T.-T. Sun, J. Hong, W. Zou, J.-Y. Fang
Development of methodology: J. Hong
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): T.-T. Sun, J. He, Q. Liang, L.-L. Ren, T.-T. Yan, T.-C. Yu, J.-Y. Tang, Y.-J. Bao, Y. Hu, Y. Lin, D. Sun, Y.-X. Chen
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): J. He, Q. Liang, H. Chen
Writing, review, and/or revision of the manuscript: T.-T. Sun, J. Hong, W. Zou, J.-Y. Fang
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): Y.-X. Chen, W. Zou
Study supervision: H. Chen
Grant Support
This project was supported (in part) by grants from the National Natural Science Foundation of China (81421001, 81320108024, 31371420, 91129724, 31271366, 31371273, 81522008, 81402347, and 81530072); the Shanghai Natural Science Foundation (grant no. 13ZR14244000); the Program for Professor of Special Appointment (Eastern Scholar No. 201268 and 2015 Youth Eastern Scholar NO.QD2015003) at Shanghai Institutions of Higher Learning; Shanghai Municipal Education Commission—Gaofeng Clinical Medicine Grant Support (no. 20152512); and Chenxing Project of Shanghai Jiao-Tong University to H. Chen and J. Hong.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.