- Open Access
Landscape of tumor suppressor long noncoding RNAs in breast cancer
Journal of Experimental & Clinical Cancer Researchvolume 38, Article number: 79 (2019)
The landscape and biological functions of tumor suppressor long noncoding RNAs in breast cancer are still unknown.
Data from whole transcriptome sequencing of 33 breast specimens in the Harbin Medical University Cancer Center cohort and The Cancer Genome Atlas was applied to identify and validate the landscape of tumor suppressor long noncoding RNAs, which was further validated by The Cancer Genome Atlas pancancer data including 33 cancer types and 12,839 patients. Next, the expression model, prognostic roles, potential biological functions and epigenetic regulation of tumor suppressor long noncoding RNAs were investigated and validated in the breast cancer and pancancer cohorts. Finally, EPB41L4A-AS2 was selected to validate our novel finding, and the tumor suppressive roles of EPB41L4A-AS2 in breast cancer were examined.
We identified and validated the landscape of tumor suppressor long noncoding RNAs in breast cancer. The expression of the identified long noncoding RNAs was downregulated in cancer tissue samples compared with normal tissue samples, and these long noncoding RNAs correlated with a favorable prognosis in breast cancer patients and the patients in the pancancer cohort. Multiple carcinogenesis-associated biological functions were predicted to be regulated negatively by these long noncoding RNAs. Moreover, these long noncoding RNAs were transcriptionally regulated by epigenetic modification, including DNA methylation and histone methylation modification. Finally, EPB41L4A-AS2 inhibited breast cancer cell proliferation, migration and invasion and induced cell apoptosis in vitro. Mechanistically, EPB41L4A-AS2, acting at least in part as a tumor suppressor, upregulated tumor suppressor gene expression. Moreover, ZNF217 recruited EZH2 to the EPB41L4A-AS2 locus and suppressed the expression of EPB41L4A-AS2 by epigenetically increasing H3K27me3 enrichment.
This work enlarges the functional landscape of known long noncoding RNAs in human cancer and provides novel insights into the suppressive roles of these long noncoding RNAs.
Breast cancer, which is the most common cause of death in women, has an incidence of 22.9% among all malignant tumors in women and causes 13.7% of cancer-related deaths among women with tumors [1,2,3]. Although different kinds of treatment, such as surgery, chemotherapy, hormonal therapy, radiotherapy and immunotherapy, have improved the prognosis of breast cancer patients to some extent, cancer-related death is still the most common cause of death [4, 5]. The main obstacles to breast cancer treatment arise from the insensitivity of cancer cells due to the overactivation of oncogenic signaling pathways, suppression of the expression of tumor suppressor genes or dormancy and subsequent reawakening of residual tumor cells, which may be a major problem with current treatment strategies [6,7,8]. Therefore, proper modulation of the expression of oncogenic or tumor suppressor genes is urgently needed to bypass the obstacles that appear in traditional breast cancer treatment.
As 75% of the genome is transcribed into noncoding RNAs, whether noncoding RNAs can serve as regulators in various biological processes has attracted increasing attention [9, 10]. Among the noncoding RNAs, long noncoding RNAs (lncRNAs) are defined as transcripts with little or no protein-coding capacity that are longer than 200 nucleotides [10, 11]. In a recent report, Li and his colleagues summarized the functional classification and experimental dissection of lncRNAs. LncRNAs participate in multiple biological processes, including histone methylation, gene transcription, protein translation, and posttranslational modification . Moreover, lncRNAs have been indicated to be associated with diverse diseases, such as Alzheimer’s disease, infectious diseases, and cancer . However, the comprehensive understanding of lncRNAs in cancer is still incomplete, and the exploration of the specific roles of lncRNAs in cancer is urgently needed.
In our previous study, the oncogenic lncRNA landscape was identified and validated in breast cancer . Moreover, we also identified an lncRNA, linc00152, which functions as a tumor propellant in breast cancer and pan-cancer . In addition, lncRNAs acting as oncogenic drivers have been reported in breast cancer and other cancers [12,13,14,15]. Lnc-BM functions as a JAK2-binding lncRNA and promotes breast cancer brain metastasis . LNMAT1 has been proven to promote lymphatic metastasis via CCL2-dependent macrophage recruitment . LncGata6 can maintain cell stemness and promote tumorigenesis . In contrast, few lncRNAs that function as tumor suppressors have been reported in breast cancer. Of these reported lncRNAs, EPB41L4A-AS2 and EGOT were first reported to suppress proliferation in breast cancer cells [16, 17]. Other lncRNAs, including GAS5 [18, 19], ZFAS1 , XIST [21, 22], and MEG3 [23,24,25], function as tumor suppressors. However, there is no comprehensive report on the landscape of tumor suppressor long noncoding RNAs (TSLNRs) in breast cancer.
In this study, we first defined the lncRNAs that are downregulated in cancer tissue, are correlated with favorable prognosis and function as tumor suppressors in breast cancer as TSLNRs. We identified and validated the comprehensive landscape of TSLNRs in breast cancer via whole transcriptome sequencing of the Harbin Medical University Cancer Center (HMUCC) cohort and The Cancer Genome Atlas (TCGA) breast cancer and pancancer data. The expression of the TSLNRs was downregulated in cancer tissue compared with normal tissue, and the expression of some TSLNRs was negatively associated with TNM stage in the pancancer data. Next, the TSLNRs were found to be correlated with a favorable prognosis in breast cancer patients and the patients in the pancancer cohort. Moreover, the TSLNRs were regulated by epigenetic modification, including DNA methylation and histone methylation modification. Finally, EPB41L4A-AS2, one of the validated TSLNRs, was selected to validate our novel finding, and its tumor suppressor roles in breast cancer were examined. As expected, EPB41L4A-AS2 inhibited breast cancer cell proliferation, migration and invasion and induced cell apoptosis partly via the upregulation of RARRES1 expression. In addition to mediating the DNA methylation regulation of EPB41L4A-AS2, ZNF217 recruited EZH2 to the EPB41L4A-AS2 locus and suppressed the expression of EPB41L4A-AS2 by epigenetically increasing H3K27me3 enrichment.
Materials and methods
Clinical sample collection and whole transcriptome sequencing
At the HMUCC, 33 samples including 15 breast cancer tissue samples, 15 adjacent normal tissue samples and 3 normal breast tissue samples were obtained from patients who had not received chemotherapy or radiotherapy. The protocol conformed to clinical research guidelines and was approved by the research ethics committee of HMUCC. Written informed consent was obtained from all of the patients who participated in this study. For whole transcriptome sequencing, ribosomal RNA was removed using Ribo-Zero™ Gold kits (Epicentre, Wisconsin, USA). Sequencing libraries were generated according to the manufacturer’s recommendations with varied index labeling using the NEBNext® Ultra™ Directional RNA Library Prep Kit for Illumina (New England Biolabs, Ipswich, USA). The libraries were then sequenced on an Illumina HiSeq 2500 platform, and 100-bp paired-end reads were generated. The raw sequencing data from this study have been deposited in the NCBI Gene Expression Omnibus (GEO) database under accession number GSE71651 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc = GSE71651s).
Access and analysis of public data
The breast cancer and pancancer genome-wide gene profiles were downloaded from TCGA (https://tcga-data.nci.nih.gov/). All transcripts were normalized by the log2 method. The correlations between genes were assessed by Spearman’s correlation coefficients according to the protocol described in a previous study . Unpaired Student’s t-tests were used to detect significant differences among tumors or between tumor and normal samples. Guilt-by-association analysis was performed to identify genes that were positively and negatively correlated with lncRNA expression [9, 27, 28]. During this analysis, the TCGA data were evaluated to generate a pairwise Spearman’s correlation between the expression of an lncRNA and that of all coding genes . Only associated coding genes with an absolute r ≥ 0.4 and a statistically significant correlation (P < 0.05) were retained. Gene ontology (GO) term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of these genes were performed using DAVID as previously described [29, 30]. H3K27me3 chromatin immunoprecipitation (ChIP)-seq data from Encyclopedia of DNA Elements (ENCODE) (www.encodeproject.org/) were analyzed. BATMAN-TCM (http://bionet.ncpsb.org/batman-tcm/) and PharmMapper (http://126.96.36.199/pharmmapper) were applied for pharmacological network analysis [31,32,33]. DNA methylation data were downloaded from the Illumina Infinium HumanMethylation450 Beadchip dataset in TCGA (https://tcga-data.nci.nih.gov/). The DNA methylation level ranged from 0 (least methylated) to 1 (most methylated). The methylation level was calculated as beta = Methylated probe intensity (M)/(Unmethylated probe intensity (U) + Methylated probe intensity (M) + 100) according to the method described in our previous study .
UACC812, BT549 and MDA-MB-453 cells were cultured in RPMI-1640 medium or DMEM (Gibco, Carlsbad, CA, USA). All media were supplemented with 10% fetal bovine serum (FBS). All cells were incubated at 37 °C in humidified air containing 5% CO2. Thereafter, all cells were cultured in regular medium for 24 h before receiving the appropriate treatment.
Plasmid and lentivirus production and infection
For the knockdown of EPB41L4A-AS2 expression, two human EPB41L4A-AS2-targeted RNAi sequences (RNAi #1225 sense: TTGATGGAGTTTCCTCTCATA; RNAi #1225 antisense: TATGAGAGGAAACTCCATCAA; RNAi #1237 sense: TTGATGGAGTTTCCTCTCATA; and RNAi# 1237 antisense: TATGAGAGGAAACTCCATCAA) were obtained from GeneChem, Co., Ltd. (Shanghai, China). A scrambled negative control siRNA and target gene siRNAs were purchased from Invitrogen (Invitrogen, CA, USA). The siRNAs were transfected into cells using Lipofectamine 3000 (Invitrogen, CA, USA) according to the manufacturer’s protocol. In addition, lentiviral particles were constructed and packaged by Shanghai GeneChem, Co., Ltd. Briefly, cells were infected with lentivirus to generate stable cell lines. After 24 h, the cells were transferred to medium containing 4 μg/ml puromycin and cultured for 3 days.
Cell viability assays
Cell viability assays were performed using the Cell Counting Kit-8 (CCK-8; Dojindo Laboratories, Kumamoto, Japan) assay according to the manufacturer’s instructions and as previously described . Briefly, cells were plated at a density of 5 × 103 cells/well in 96-well plates. CCK-8 solution was added to each well and incubated at 37 °C for 90 min. Then, the absorbance of the cell suspension was measured at a wavelength of 450 nm. Medium containing 10% CCK-8 served as a negative control.
Colony formation assays
Cells were plated into a 6-well plate and cultured in medium containing 10% FBS for two weeks. Colonies were fixed with methanol for 30 min, and 500 μl of 0.5% crystal violet was subsequently added (Sigma-Aldrich, St. Louis, MO, USA) to each well for 30 min for visualization and counting.
Wound healing assays and invasion assays
For wound healing assays, cells were seeded in plates in DMEM containing 1% FBS. Monolayers were wounded with a sterile 10-μl pipette tip, and phosphate-buffered saline (PBS) was applied to remove the suspended cells. Images were captured of each well at different time points. The migration distance was estimated based on the width of the wounds using ImageJ software. For invasion assays, cells were seeded in plates coated with Matrigel (Sigma-Aldrich, USA) according to the manufacturer’s instructions. Medium containing 10% FBS was added to the lower chamber for 24 h. Cotton swabs were used to remove the noninvading cells that remained in the top chamber, and the cells that had invaded to the underside of the membrane were fixed in 100% methanol for 30 min and stained with 0.5% crystal violet. Finally, images were obtained under a light microscope.
RNA preparation and qRT-PCR
Total RNA was extracted using TRIzol Reagent according to the manufacturer’s protocol (Invitrogen, Beijing, China), and cDNA was synthesized using the PrimeScript RT Reagent Kit with gDNA Eraser (Takara Bio, Otsu, Japan). mRNA expression was examined by real-time PCR using FastStart Universal SYBR Green Master Mix (Roche, Mannheim, Germany) with gene-specific primers and the ABI 7500 Fast Real-time PCR Detection System (Applied Biosystems, Foster City, CA, USA). The results were normalized to the expression of GAPDH. The primer sequences were as follows: EPB41L4A-AS2-F: 5′-CGGAGCAGGTGCAATCTGT-3′; EPB41L4A-AS2-R: 5′-CCCTCGTGTCTCCCCTAACTG-3′; RARRES1-F: 5′-AAACCCCTTGGAAATAGTCAGC-3′; RARRES1-R: 5′-GGAAAGCCAAATCCCAGATGAG-3′; ST18-F: 5′-CAAACCACCTAGAGTCCCAAAG-3′; ST18-R: 5′-ACACCTGTTCTCACAAGGGATA-3′; BMP4-F: 5′-CAAACCACCTAGAGTCCCAAAG-3′; BMP4-R: 5′-GACGGCACTCTTGCTAGGC-3′; FOXA2-F: 5′-GGAGCAGCTACTATGCAGAGC-3′; FOXA2-R: 5′-CGTGTTCATGCCGTTCATCC-3′; FOXL1-F: 5′-GCCTCGCCCATGCTGTATC-3′; FOXL1-R: 5′-CGTTGAGCGTGACCCTCTG-3′; LRIG1-F: 5′- GGACTTGCCGAACCTACAGG-3′; LRIG1-R: 5′-GCTGCGAATCTTGTTGTGCTG-3′; RASSF1-F: 5′-ATGTGCCTACCTGATACTTT-3′; RASSF1-R: 5′-ATGGTGAACCTGGAGAAC-3′; GAPDH-F: 5′-CATGTTCGTCATGGGTGTGAA-3′; and GAPDH-R: 5′-GGCATGGACTGTGGTCATGAG-3′.
Flow cytometry analysis
Flow cytometry analysis was performed according to the manufacturer’s protocol. Cells were seeded in 6-well plates for 24 h, harvested and washed twice with cold PBS. Then, the cells were stained with PE-conjugated Annexin V and 7-AAD (BD Biosciences, San Jose, CA, USA) for 15 min at room temperature in the dark. Finally, the cells were analyzed using a FACSCalibur flow cytometer (BD Biosciences, San Jose, CA, USA).
Co-IP and western blot assays
Coimmunoprecipitation (Co-IP) and western blot assays were carried out by using the Pierce™ Crosslink Magnetic IP/Co-IP Kit (Thermo Fisher Scientific, USA) according to the manufacturer’s protocol . Anti-ZNF217 (ab117798) and anti-EZH2 (ab186006) antibodies were obtained from Abcam. An anti-tubulin (sc-73,242, 1:1000) antibody was purchased from Santa Cruz Biotechnology. First, cells were lysed with lysis buffer and 2% SDS. Next, the protein concentrations were evaluated using a protein assay kit (Bio-Rad, Richmond, CA), and equal amounts of protein were separated by SDS-PAGE, followed by electroblotting onto a nitrocellulose membrane, which was then blocked with 5% nonfat milk in 0.1% Tween 20-PBS overnight at 4 °C. The membrane was immunoblotted with primary antibodies against Bcl-2 (ab32124), Bax (ab182733) and RARRES1 (ab76530). After washing with Tween 20/PBS, the membrane was incubated with horseradish peroxidase-conjugated secondary antibodies for 1 h. Finally, protein bands on the membrane were visualized by an enhanced chemiluminescence western blotting detection system (Western Lightning; Perkin-Elmer, Norwalk, CT).
Chromatin immunoprecipitation assays
Chromatin immunoprecipitation (ChIP) assays were performed using the ChIP Assay Kit (Beyotime, Shanghai, China) according to the manufacturer’s protocol with slight modifications. Cells were crosslinked with 1% formaldehyde, and this reaction was terminated after 10 min by the addition of glycine at a final concentration of 0.125 M. DNA was immunoprecipitated from the sonicated cell lysates using anti-ZNF217 and anti-H3K27me3 antibodies(Abcam); IgG (BD Biosciences, San Diego, CA, USA) served as the negative control. The DNA was subjected to PCR to amplify the ZNF217 and H3K27me3 binding sites. The amplified fragments were then analyzed on an agarose gel. Chromatin (10%) prior to immunoprecipitation was used as the input control.
Overall survival (OS) and disease-free survival (DFS) were calculated as the time from surgery until the occurrence of death and relapse, respectively. The expression of lncRNAs was dichotomized using the median expression as the cut-off to define high values (at or above the median) versus low values (below the median). The differences between groups in the results of the in vitro experiments were analyzed using Student’s t-test. Spearman’s correlation coefficients were calculated for correlation analyses. All of the experiments were performed independently in triplicate. R software was used to prepare all figures. All statistical tests were two-sided, and P < 0.05 indicated statistical significance.
Discovery and validation of the TSLNR landscape in breast cancer
First, 33 breast tissue samples in the HMUCC cohort including 15 paired cancer tissue and corresponding adjacent normal tissue samples and 3 completely normal breast tissue samples acquired during mammoplasty, were subjected to whole transcriptome sequencing. The aberrant differential expression of lncRNAs among these tissues was obtained (Fig. 1a and Additional file 1: Table S1). Moreover, the differential expression of lncRNAs from breast cancer patients in the TCGA portal was also investigated in detail (Fig. 1b and Additional file 1: Table S2). Next, the overlapping differentially expressed lncRNAs between the HMUCC cohort and the TCGA cohort were investigated, and 137 differentially expressed lncRNAs were identified in both cohorts (Fig. 1c and Additional file 1: Table S3). Finally, the filter criteria for differentially downregulated lncRNAs was further restricted in the HMUCC cohort, and 19 differentially downregulated lncRNAs were identified and then validated in the TCGA breast cancer cohort (Fig. 1d and Additional file 1: Table S4). Thus, we defined these 19 differentially downregulated lncRNAs as TSLNRs.
Extension validation of TSLNR expression in the pancancer cohort
Next, the expression of the 19 TSLNRs identified above were revalidated in the HMUCC and TCGA pancancer data that included 33 cancer types and 12,839 patients. In the HMUCC dataset, 11 lncRNAs (WWC2-AS2, PGM5-AS1, MEG3, LINC00667, LINC-PINT, HAND2-AS1, FTX, FAM66C, EPB41L4A-AS2, CYP1B1-AS1 and ACVR2B-AS1) exhibited downregulated expression in breast cancer tissue compared with normal tissue (Additional file 2: Figure S1). In the TCGA dataset, 16 lncRNAs (WWC2-AS2, WEE2-AS1, TRHDE-AS1, PGM5-AS1, NR2F1-AS1, MEG3, LINC00667, LINC-PINT, HCG11, HAND2-AS1, GRIK1-AS1, FTX, FAM66C, EPB41L4A-AS2, CYP1B1-AS1 and ACVR2B-AS1) exhibited significantly downregulated expression in cancer tissue compared to normal tissue, consistent with the breast cancer expression model (Fig. 2a). Genetic alteration was also examined for these lncRNAs and there was no significant amplification or deep deletion in genome in breast cancer (Additional file 1: Table S4; Additional file 2: Figure S2). Furthermore, the relationships between the expression of these lncRNAs and the TNM stages of the patients in the pancancer data from the TCGA cohort were examined, and the expression of eight lncRNAs (WWC2-AS2, WEE2-AS1, TRHDE-AS1, TPT1-AS1, PGM5-AS1, HAND2-AS1, GRIK1-AS1 and EPB41L4A-AS2) was negatively related to advanced stages of cancer, but the relationships of three lncRNAs (WWC2-AS2, TRHDE-AS1 and HAND2-AS1) were not statistically significant (Fig. 2b). Five of the lncRNAs (WEE2-AS1, TPT1-AS1, PGM5-AS1, GRIK1-AS1 and EPB41L4A-AS2) were expressed at significantly higher levels in the patients with early TNM stage (stage I & II) disease than in those with advanced stage (stage III & IV) disease in the pancancer data from TCGA (Fig. 2c).
TSLNR expression correlated with favorable prognosis in the breast cancer and pancancer datasets
Whether the TSLNRs discovered in this study are associated with the prognosis of breast cancer patients was unknown. Thus, we examined the relationships between the expression of the TSLNRs and OS or DFS in breast cancer patients in detail. The results showed that high expression of ACVR2B-AS1 or WEE2-AS1 was associated with favorable OS in breast cancer patients in the TCGA cohort (Additional file 2: Figure S3). Moreover, high expression of ACVR2B-AS1, WEE2-AS1, LINC-PINT or HAND2-AS1 was associated with favorable DFS in breast cancer patients in the TCGA cohort (Additional file 2: Figure S3). For the other lncRNAs, high expression of CYP1B1-AS1, LINC-PINT, LINC00667 or GRIK1-AS1 was associated with favorable OS, but the improvements were not statistically significant (Additional file 2: Figure S4). Finally, high expression of CYP1B1-AS1, FAM66C or GRIK1-AS1 was associated with favorable OS, but the improvements were not statistically significant (Additional file 2: Figure S3).
To validate the results obtained in the breast cancer cohorts, the TCGA pancancer data were applied to examine the relationship between TSLNR expression and OS or DFS. The results showed that high expression of EPB41L4A-AS2, WEE2-AS1, ACVR2B-AS1, PGM5-AS1, LINC00667, CYP1B1-AS1, FAM66C, FTX, SNHG5, HCG11, WWC2-AS2 or TRHDE-AS1 was associated with favorable OS in the patients in the pancancer data from TCGA (Fig. 3). Moreover, high expression of EPB41L4A-AS2, HCG11, CYP1B1-AS1, ACVR2B-AS1, SNHG5 or LINC-PINT was associated with favorable DFS in patients in this cohort (Fig. 4a-f).
Potential biological functions of TSLNRs in breast cancer
To explore the potential biological functions of TSLNRs, guilt-by-association analysis was applied to perform the following analyses (Additional file 1: Table S5). TSLNRs may negatively regulate multiple tumor biological behaviors, including cell proliferation, angiogenesis, cell migration, cell-matrix adhesion, Wnt signaling transduction, mitotic cell cycle phase transition, JAK-STAT signaling transduction, tumor necrosis factor (TNF) production, BMP signaling transduction, cell adhesion mediated by integrin, cAMP biosynthesis, phagocytosis, Rho protein signal transduction, and platelet-derived growth factor receptor signaling transduction (Fig. 4g and Additional file 1: Table S5). The pathways involving TSLNRs were further examined. The results indicated that TSLNRs may be involved in several vital oncogenic signaling pathways, including the PI3K-Akt signaling pathway, the Ras signaling pathway, proteoglycans in cancer, cytokine-cytokine receptor interactions, the Rap1 signaling pathway, the TGF-beta signaling pathway, the Hippo signaling pathway, the cGMP-PKG signaling pathway, the MAPK signaling pathway, the PPAR signaling pathway, the Hedgehog signaling pathway, the TNF signaling pathway, the NF-kappa B signaling pathway (Fig. 4h and Additional file 1: Table S5).
Epigenetic modification leads to the downregulation of TSLNR expression in breast cancer
Why is the expression of these TSLNRs downregulated in both the human breast cancer data and the pancancer data? The Illumina Infinium HumanMethylation450 Beadchip data in the TCGA portal was downloaded and investigated carefully to explore the beta value differences between cancer tissues and normal tissues for each TSLNR locus. The results showed that 12 TSLNR genome loci (those of WWC2-AS2, TRHDE-AS1, SMAD1-AS1, PGM5-AS1, NR2F1-AS1, MEG3, HCG11, HAND2-AS1, FTX, FAM66C, EPB41L4A-AS2 and CYP1B1-AS1) exhibited higher levels of DNA methylation in cancer tissues than normal tissues (Fig. 5a). Thus, the low expression of TSLNRs, at least in part, may be the result of the hypermethylation of each TSLNR genome locus in breast cancer.
Histone methylation modification was next investigated as it is usually accompanied by DNA methylation. The H3K27me3 enrichment peak for each TSLNR genome locus in MDA-MB-231 cells was investigated in the ENCODE data. As expected, 12 TSLNRs (WWC2-AS2, WEE2-AS1, PGM5-AS1, NR2F1-AS1, LINC-PINT, HCG11, FTX, FAM66C, EPB41L4A-AS2, SMAD5-AS1, TPT1-AS1 and SNHG5) showed significant H3K27me3 enrichment peaks at the corresponding TSLNR locus (Fig. 5b). Thus, the H3K27me3 histone methylation modification may also lead to the low expression of TSLNRs in breast cancer.
Next, EPB41L4A-AS2 was selected to validate the histone methylation modification model, as we first reported the potential function of EPB41L4A-AS2 in human cancer , and an obvious H3K27me3 enrichment peak at the EPB41L4A-AS2 locus could be observed in MDA-MB-231 cells (Fig. 5b). ZNF217 has been reported to be a marker of poor prognosis in breast cancer that drives epithelial-mesenchymal transition and invasion by recruiting EZH2 to its target genes, which are marked with an H3K27me3 enrichment peak [35, 36]. Thus, we hypothesized that EPB41L4A-AS2 could be regulated by this model. Initially, the expression of EPB41L4A-AS2 was upregulated in MDA-MB-231 breast cancer cells with the knockdown of ZHF217 expression (Fig. 6a and b). Moreover, EPB41L4A-AS2 expression was also found to be downregulated in MDA-MB-231 breast cancer cells overexpressing ZNF217 in the GEO dataset GSE35511 (Additional file 2: Figure S5). Next, a Co-IP assay showed that ZNF217 can directly bind to EZH2 (Fig. 6c). ChIP followed by PCR showed that EZH2 could bind to the promoter region of EPB41L4A-AS2 (Fig. 6d). Furthermore, H3K27me3 was also found to be enriched in the EPB41L4A-AS2 locus via ChIP (Fig. 6e). Finally, to extend our results, the data from ENCODE were examined again, and the results showed that a significant H3K27me3 enrichment peak could be observed at the EPB41L4A-AS2 locus in breast cancer cells (Fig. 6f). Thus, the ZNF217-EZH2-H3K27me3 axis epigenetically suppresses the expression of EPB41L4A-AS2 in breast cancer.
Potential biological functions of EPB41L4A-AS2 in breast cancer
A guilt-by-association approach was applied to analyze the potential biological functions of EPB41L4A-AS2 in breast cancer with the TCGA breast cancer data (Additional file 1: Table S6). EPB41L4A-AS2 may negatively regulate multiple biological functions involved in tumors, including cell proliferation, cell migration, TGF-beta signaling transduction, BMP signaling transduction and Notch signaling transduction (Fig. 7a). KEGG analysis showed that EPB41L4A-AS2 may participate in multiple oncogenic signaling pathways, including metabolic pathways, the cAMP signaling pathway, the MAPK signaling pathway, and the estrogen signaling pathway (Fig. 7b). Finally, to extend our results, the potential biological functions of EPB41L4A-AS2 in the pancancer dataset were investigated. The results indicated that EPB41L4A-AS2 may be involved in the key biological functions and pathways of carcinogenesis and tumor progression, and these results agreed with the results of the breast cancer analysis (Fig. 7c and d and Additional file 1: Table S7). Moreover, pharmacological network analysis indicates that EPB41L4A-AS2 may be involved in the regulation of paclitaxel activity in breast cancer (Additional file 2: Figure S6), which is in agreement with the predicted roles of EPB41L4A-AS2 in junction assembly and mitosis (Additional file 1: Table S6).
As EPB41L4A-AS2 was predicted to suppress cancer cell proliferation both in the breast cancer and pancancer TCGA data, we next wondered whether EPB41L4A-AS2 expression was negatively correlated with the expression of classic oncogenes in human cancer. As MKI67 and MYC are two famous driver oncogenes in carcinogenesis and tumor progression [37, 38], the relationships between the expression levels of EPB41L4A-AS2, MKI67 and MYC were next explored. As expected, the expression of EPB41L4A-AS2 was negatively correlated with the expression of MKI67 and MYC both in breast cancer patients and the patients in the pancancer cohort (Fig. 7e-h). Collectively, these results indicate that EPB41L4A-AS2 is negatively related to the malignant behaviors of cancer and negatively correlated with oncogene expression.
EPB41L4A-AS2 inhibited cell proliferation, migration and invasion in breast cancer
First, the expression level of EPB41L4A-AS2 was examined in eight breast cancer cell lines (Additional file 2: Figure S7a). UACC812, BT549 and MDA-MB-453 cells were ultimately selected to perform the following experiments as UACC812 and BT549 cells have low expression of EPB41L4A-AS2 and MDA-MD-453 has high expression of EPB41L4A-AS2. Next, the overexpression and knockdown efficiencies were examined in these three cell lines (Additional file 2: Figure S7b-d). Compared with Flag (control) overexpression, the overexpression of EPB41L4A-AS2 inhibited cell proliferation significantly in UACC812 and BT549 cells (Fig. 8a and b). In contrast, knocking down EPB41L4A-AS2 expression promoted cell proliferation in MDA-MB-453 cells (Fig. 8c). Furthermore, similar results that showed that the overexpression of EPB41L4A-AS2 inhibited clone formation while knocking down EPB41L4A-AS2 expression increased clone formation were observed in clone formation assays with the above three cell lines (Fig. 8d-f).
With regards to cell apoptosis, overexpressing EPB41L4A-AS2 triggered more cell apoptosis, while knocking down EPB41L4A-AS2 suppressed cell apoptosis in the three cell lines, as determined by flow cytometry analysis (Fig. 8g-i). Finally, cell apoptosis-associated markers were also investigated. Bcl-2 expression was reduced, and Bax expression was increased in the UACC812 and BT549 cells overexpressing EPB41L4A-AS2 (Fig. 8j). In MDA-MB-453 cells, knocking down EPB41L4A-AS2 expression increased Bcl-2 expression and reduced Bax expression (Fig. 8k). As EPB41L4A-AS2 was shown to inhibit the proliferation of breast cancer cells, we then investigated whether it could suppress cell migration and invasion in the abovementioned cell lines via wound healing and transwell assays. As expected, the overexpression of EPB41L4A-AS2 reduced migration (Fig. 9a) and invasion in UACC812 and BT549 cells (Fig. 9b). In addition to EPB41L4A-AS2, three other lncRNAs, MEG3, WEE2-AS1 and HAND2-AS1, were also applied to examine their potential tumor suppressor roles in breast cancer. As expected, the overexpression of each lncRNA inhibited clone formation in UACC812 cells (Additional file 2: Figure S8). Moreover, the overexpression of each lncRNA inhibited cell viability in UACC812 cells, matching the EPB41L4A-AS2 results (Additional file 2: Figure S8).
EPB41L4A-AS2 functions as a tumor suppressor via the upregulation of RARRES1 expression
Whole transcriptome sequencing of the UACC812 and BT549 breast cancer cell lines with or without EPB41L4A-AS2 overexpression was performed to identify potential target genes of EPB41L4A-AS2. The differentially expressed genes were measured and identified in different groups (Fig. 10a and Additional file 1: Table S8). Based on the KEGG analysis of these differentially expressed genes, EPB41L4A-AS2 may be involved in multiple cancer-associated pathways such as the hormone signaling pathway, FoxO signaling pathway and PKG signaling pathway in UACC812 cells (Fig. 10b). According to the KEGG analysis, EPB41L4A-AS2 may also be involved with the proteoglycans in the cancer-related Wnt signaling pathway, and in the TGF-beta signaling pathway, Notch signaling pathway, Hippo signaling pathway and AMPK signaling pathway in BT549 cells (Fig. 10c), and these results, to some extent, were consistent with those of our previous analysis (Fig. 7a-d). Next, a PCR array was applied to verify the differentially expressed genes identified by transcriptome sequencing, and we found that the expression levels of seven genes (RARRES1, ST18, BMP4, FOXA2, FOXL1, LRIG1 and RASSF1), which have all been reported to be tumor suppressors in breast cancer, were upregulated significantly with the overexpression of EPB41L4A-AS2 in the UACC812 and BT549 breast cancer cell lines (Fig. 10e). Moreover, RARRES1 protein expression was increased in the UACC812 cells overexpressing EPB41L4A-AS2 (Fig. 10e) and decreased in the MDA-MB-453 cells with knocked down EPB41L4A-AS2 expression (Fig. 10f). Thus, EPB41L4A-AS2 functions as a tumor suppressor at least in part via the upregulation of RARRES1 expression in breast cancer.
To our knowledge, no study has focused on the systematic identification of TSLNRs in human diseases, especially in human breast cancer. Moreover, few lncRNAs that function as tumor suppressors and are associated with a favorable prognosis have been identified in breast cancer. Thus, to solve this issue, TSLNRs were identified comprehensively and validated by whole transcriptome sequencing and analysis of TCGA breast cancer data via bioinformatic approaches and biological experiments in this study. We defined the identified lncRNAs as TSLNRs, enlarging the landscape of the functional lncRNA category. The levels of these TSLNRs are downregulated in cancer tissues compared with normal tissues. In addition, they were found to be correlated with a favorable prognosis in breast cancer patients and the patients in the pancancer cohort. Moreover, these TSLNRs are regulated by epigenetic modification, including DNA methylation and histone methylation modification. Finally, EPB41L4A-AS2, one of the verified TSLNRs, was selected to validate our novel findings, and it’s the tumor suppressor roles of EPB41L4A-AS2 in breast cancer were examined. In addition to mediating the DNA methylation regulation of EPB41L4A-AS2, ZNF217 recruited EZH2 to the EPB41L4A-AS2 locus and suppressed the expression of EPB41L4A-AS2 by epigenetically increasing H3K27me3 enrichment. Moreover, EPB41L4A-AS2 inhibited breast cancer cell proliferation, migration and invasion and induced cell apoptosis. To some extent, EPB41L4A-AS2 functioned as a tumor suppressor by upregulating the expression of RARRES1 in breast cancer.
Among the 19 TSLNRs in this study, only eight (EPB41L4A-AS2, MEG3, LINC-PINT, FTX, HCG11, HAND2-AS1, SNHG5 and TPT1-AS1) have been reported previously to be associated with malignancy, and the potential functions of the other 11 lncRNAs in human cancer remain a mystery. We were the first to report that EPB41L4A-AS2 inhibited tumor proliferation and was associated with favorable prognoses in breast cancer and other solid tumors . Moreover, EPB41L4A-AS2 was also reported to inhibit proliferation and invasion and promote cell apoptosis in non-small cell lung cancer , in accordance with our results. Among the other identified lncRNAs, MEG3 was reported to suppress cell proliferation, invasion, and angiogenesis through the AKT pathway or the transcriptional activity of p53 in breast cancer [23,24,25, 40, 41]. LINC-PINT, which is regulated by p53, inhibits tumor cell invasion through a highly conserved sequence element [42, 43]. FTX was reported to inhibit hepatocellular carcinoma proliferation and metastasis by binding MCM2 and miR-374a . In addition, HCG11 was found to suppress apoptosis via MAPK signaling transduction in hepatocellular carcinoma, and the downregulation of HCG11 expression predicts a poor prognosis in prostate cancer [45, 46]. A study also revealed that SNHG5 suppresses gastric cancer progression by trapping MTA2 in the cytosol . In contrast, other studies have indicated that SNHG5 functions as an oncogene in colorectal cancer, osteosarcoma, hepatocellular cancer, bladder cancer and gastric cancer [48,49,50,51,52]. TPT1-AS1 was reported to promote cell growth and metastasis in cervical cancer by acting as a sponge for miR-324-5p . However, our study suggested that TPT1-AS1 may be a tumor suppressor in breast cancer. The functional diversity of SNHG5 and TPT1-AS1 in various malignancies may be caused by tumor heterogeneity, and subsequent research is necessary to reveal the potential roles of these two lncRNAs.
This study addresses three distinct merits that are not evident in other studies. First, a panel of TSLNRs was identified via whole transcriptome sequencing and TCGA data analysis, thus providing a more efficient research methodology and opening research into this field in human cancer. Hence, the drawbacks of a single lncRNA-focused study can be avoided. Second, TSLNRs were comprehensively further validated in the TCGA pancancer data containing 33 cancer types and 12,839 patients, increasing the study reliability and strengthening our theoretical basis. Third, DNA methylation and histone methylation modification led to the downregulation of TSLNR expression in breast cancer, further suggesting that DNA methylation inhibitors or histone methylation inhibitors may be effective agents in human cancer management. To date, DNA methyltransferase inhibitors have been reported to inhibit tumor growth and promote cell apoptosis in breast cancer, colon cancer, ovarian cancer, gastric cancer and human endometrial cancer [54,55,56,57,58]. Moreover, the histone methyltransferase inhibitor A-366 has been reported to suppress cancer cell viability by targeting G9a/GLP in leukemia . In addition, 3-deazaneplanocin A, an inhibitor of the histone methyltransferase EZH2, inhibits the growth of non-small cell lung cancer cells . Thus, the anticancer roles of epigenetic therapies targeting DNA methylation or histone methylation modification could be result from direct tumor suppressor functions and the reversed expression of TSLNRs, which also function as effective antitumor agents.
We acknowledge several limitations of our study. First, not all of the gene expression profiles in TCGA can be accessed. Thus, the TSLNR landscape in human breast cancer described in this study may not be sufficient. Second, the biological functions and detailed mechanisms of other TSLNRs were not elucidated. Further in-depth efforts will be needed to investigate TSLNRs other than EPB41L4A-AS2 and focus on molecular mechanisms. Third, pharmacological network analysis indicated that EPB41L4A-AS2 may be involved in the regulation of paclitaxel activity in breast cancer; however, this result should be examined and validated in our future investigation. Lastly, as gene therapy has its own complications, whether the upregulation of TSLNR expression with or without epigenetic therapy could be applied in cancer therapy requires more clinical trials to ensure safety and specificity.
This work enlarges the functional landscape of known lncRNAs in human cancer and provides novel insights into their suppressive roles.
Encyclopedia of DNA Elements
Head and neck squamous cell carcinoma
Kyoto Encyclopedia of Genes and Genomes
Small interfering RNA
The Cancer Genome Atlas
tumor suppressor long noncoding RNA
DeSantis C, Ma J, Bryan L, Jemal A. Breast cancer statistics, 2013. CA Cancer J Clin. 2014;64(1):52–62.
Ponde NF, Zardavas D, Piccart M. Progress in adjuvant systemic therapy for breast cancer. Nat Rev Clin Oncol. 2018.
Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA Cancer J Clin. 2015;65(2):87–108.
Song M, Vogelstein B, Giovannucci EL, Willett WC, Tomasetti C. Cancer prevention: molecular and epidemiologic consensus. Science. 2018;361(6409):1317–8.
Aguirre-Ghiso JA. How dormant cancer persists and reawakens. Science. 2018;361(6409):1314–5.
Landry BD, Leete T, Richards R, Cruz-Gordillo P, Schwartz HR, Honeywell ME, et al. Tumor-stroma interactions differentially alter drug sensitivity based on the origin of stromal cells. Mol Syst Biol. 2018;14(8):e8322.
Deshmukh AA, Shirvani SM, Lal L, Swint JM, Cantor SB, Smith BD, et al. Cost-effectiveness analysis comparing conventional, Hypofractionated, and intraoperative radiotherapy for early-stage breast Cancer. J Natl Cancer Inst. 2017;109(11).
Gomis RR, Gawrzak S. Tumo r cell dormancy. Mol Oncol. 2016;11(1):62-78.
Li Y, Xu J, Shao T, Zhang Y, Chen H, Function Prediction LXRNA. Methods Mol Biol. 2017;1654:17–28.
Xu S, Kong D, Chen Q, Ping Y, Pang D. Oncogenic long noncoding RNA landscape in breast cancer. Mol Cancer. 2017;16(1):129.
Xu S, Wan L, Yin H, Xu H, Zheng W, Shen M, et al. Long noncoding RNA Linc00152 functions as a tumor propellant in pan-Cancer. Cell Physiol Biochem. 2017;44(6):2476–90.
Chen R, Liu Y, Zhuang H, Yang B, Hei K, Xiao M, et al. Quantitative proteomics reveals that long non-coding RNA MALAT1 interacts with DBC1 to regulate p53 acetylation. Nucleic Acids Res. 2017;45(17):9947–59.
Han Z, Shi L. Long non-coding RNA LUCAT1 modulates methotrexate resistance in osteosarcoma via miR-200c/ABCB1 axis. Biochem Biophys Res Commun. 2018;495(1):947–53.
Ren Y, Jia HH, Xu YQ, Zhou X, Zhao XH, Wang YF, et al. Paracrine and epigenetic control of CAF-induced metastasis: the role of HOTAIR stimulated by TGF-ss1 secretion. Mol Cancer. 2018;17(1):5.
Wu J, Zhao W, Wang Z, Xiang X, Zhang S, Liu L. Long non-coding RNA SNHG20 promotes the tumorigenesis of oral squamous cell carcinoma via targeting miR-197/LIN28 axis. J Cell Mol Med. 2019;23(1):680–8.
Xu S, Wang P, You Z, Meng H, Mu G, Bai X, et al. The long non-coding RNA EPB41L4A-AS2 inhibits tumor proliferation and is associated with favorable prognoses in breast cancer and other solid tumors. Oncotarget. 2016;7(15):20704–17.
Xu SP, Zhang JF, Sui SY, Bai NX, Gao S, Zhang GW, et al. Downregulation of the long noncoding RNA EGOT correlates with malignant status and poor prognosis in breast cancer. Tumour Biol. 2015;36(12):9807–12.
Li S, Zhou J, Wang Z, Wang P, Gao X, Wang Y. Long noncoding RNA GAS5 suppresses triple negative breast cancer progression through inhibition of proliferation and invasion by competitively binding miR-196a-5p. Biomed Pharmacother. 2018;104:451–7.
Mourtada-Maarabouni M, Pickard MR, Hedge VL, Farzaneh F, Williams GT. GAS5, a non-protein-coding RNA, controls apoptosis and is downregulated in breast cancer. Oncogene. 2009;28(2):195–208.
Fan S, Fan C, Liu N, Huang K, Fang X, Wang K. Downregulation of the long non-coding RNA ZFAS1 is associated with cell proliferation, migration and invasion in breast cancer. Mol Med Rep. 2018;17(5):6405–12.
Huang YS, Chang CC, Lee SS, Jou YS, Shih HM. Xist reduction in breast cancer upregulates AKT phosphorylation via HDAC3-mediated repression of PHLPP1 expression. Oncotarget. 2016;7(28):43256–66.
Zheng R, Lin S, Guan L, Yuan H, Liu K, Liu C, et al. Long non-coding RNA XIST inhibited breast cancer cell growth, migration, and invasion via miR-155/CDX1 axis. Biochem Biophys Res Commun. 2018;498(4):1002–8.
Sun L, Li Y, Yang B. Downregulated long non-coding RNA MEG3 in breast cancer regulates proliferation, migration and invasion by depending on p53's transcriptional activity. Biochem Biophys Res Commun. 2016;478(1):323–9.
Zhang CY, Yu MS, Li X, Zhang Z, Han CR, Yan B. Overexpression of long non-coding RNA MEG3 suppresses breast cancer cell proliferation, invasion, and angiogenesis through AKT pathway. Tumour Biol. 2017;39(6):1010428317701311.
Zhang W, Shi S, Jiang J, Li X, Lu H, Ren F. LncRNA MEG3 inhibits cell epithelial-mesenchymal transition by sponging miR-421 targeting E-cadherin in breast cancer. Biomed Pharmacother. 2017;91:312–9.
Mitra R, Chen X, Greenawalt EJ, Maulik U, Jiang W, Zhao Z, et al. Decoding critical long non-coding RNA in ovarian cancer epithelial-to-mesenchymal transition. Nat Commun. 2017;8(1):1604.
Khadirnaikar S, Narayanan SP, Shukla SK. Decoding the LncRNA transcriptome of esophageal cancer: identification of clinically relevant LncRNAs. Biomark Med. 2018;12(10):1083–93.
Lefever S, Anckaert J, Volders PJ, Luypaert M, Vandesompele J, Mestdagh P. decodeRNA- predicting non-coding RNA functions using guilt-by-association. Database (Oxford). 2017.
Huang d W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
Zhang B, Horvath S. A general framework for weighted gene co-expression network analysis. Stat Appl Genet Mol Biol. 2005;4:Article17.
Liu X, Ouyang S, Yu B, Liu Y, Huang K, Gong J, et al. PharmMapper server: a web server for potential drug target identification using pharmacophore mapping approach. Nucleic Acids Res. 2010;38(Web Server issue):W609–14.
Liu Z, Guo F, Wang Y, Li C, Zhang X, Li H, et al. BATMAN-TCM: a bioinformatics analysis tool for molecular mechANism of traditional Chinese medicine. Sci Rep. 2016;6:21146.
Cui Y, Li C, Zeng C, Li J, Zhu Z, Chen W, et al. Tongmai Yangxin pills anti-oxidative stress alleviates cisplatin-induced cardiotoxicity: network pharmacology analysis and experimental evidence. Biomed Pharmacother. 2018;108:1081–9.
Shan L, Zhou X, Liu X, Wang Y, Su D, Hou Y, et al. FOXK2 elicits massive transcription repression and suppresses the hypoxic response and breast Cancer carcinogenesis. Cancer Cell. 2016;30(5):708–22.
Littlepage LE, Adler AS, Kouros-Mehr H, Huang G, Chou J, Krig SR, et al. The transcription factor ZNF217 is a prognostic biomarker and therapeutic target during breast cancer progression. Cancer Discov. 2012;2(7):638–51.
Collins C, Volik S, Kowbel D, Ginzinger D, Ylstra B, Cloutier T, et al. Comprehensive genome sequence analysis of a breast cancer amplicon. Genome Res. 2001;11(6):1034–42.
Ren Y, Bi C, Zhao X, Lwin T, Wang C, Yuan J, et al. PLK1 stabilizes a MYC-dependent kinase network in aggressive B cell lymphomas. J Clin Invest. 2018.
Sales Gil R, Vagnarelli P. Ki-67: more hidden behind a 'Classic proliferation Marker'. Trends Biochem Sci. 2018;43(10):747–8.
Shu J, Li S, Chen YB, Zhu QF, Yu XH. Long non-coding RNA EPB41L4A-AS2 inhibited non-small cell lung cancer proliferation, invasion and promoted cell apoptosis. Neoplasma. 2018;65(5):664–72.
Cui X, Yi Q, Jing X, Huang Y, Tian J, Long C, et al. Mining prognostic significance of MEG3 in human breast Cancer using bioinformatics analysis. Cell Physiol Biochem. 2018;50(1):41–51.
Zhang JJ, Guo SH, Jia BQ. Down-regulation of long non-coding RNA MEG3 serves as an unfavorable risk factor for survival of patients with breast cancer. Eur Rev Med Pharmacol Sci. 2016;20(24):5143–7.
Garitano-Trojaola A, Jose-Eneriz ES, Ezponda T, Unfried JP, Carrasco-Leon A, Razquin N, et al. Deregulation of linc-PINT in acute lymphoblastic leukemia is implicated in abnormal proliferation of leukemic cells. Oncotarget. 2018;9(16):12842–52.
Marin-Bejar O, Mas AM, Gonzalez J, Martinez D, Athie A, Morales X, et al. The human lncRNA LINC-PINT inhibits tumor cell invasion through a highly conserved sequence element. Genome Biol. 2017;18(1):202.
Liu F, Yuan JH, Huang JF, Yang F, Wang TT, Ma JZ, et al. Long noncoding RNA FTX inhibits hepatocellular carcinoma proliferation and metastasis by binding MCM2 and miR-374a. Oncogene. 2016;35(41):5422–34.
Xu Y, Zheng Y, Liu H, Li T. Modulation of IGF2BP1 by long non-coding RNA HCG11 suppresses apoptosis of hepatocellular carcinoma cells via MAPK signaling transduction. Int J Oncol. 2017;51(3):791–800.
Zhang Y, Zhang P, Wan X, Su X, Kong Z, Zhai Q, et al. Downregulation of long non-coding RNA HCG11 predicts a poor prognosis in prostate cancer. Biomed Pharmacother. 2016;83:936–41.
Zhao L, Guo H, Zhou B, Feng J, Li Y, Han T, et al. Long non-coding RNA SNHG5 suppresses gastric cancer progression by trapping MTA2 in the cytosol. Oncogene. 2016;35(44):5770–80.
Damas ND, Marcatti M, Come C, Christensen LL, Nielsen MM, Baumgartner R, et al. SNHG5 promotes colorectal cancer cell survival by counteracting STAU1-mediated mRNA destabilization. Nat Commun. 2016;7:13875.
Ju C, Zhou R, Sun J, Zhang F, Tang X, Chen KK, et al. LncRNA SNHG5 promotes the progression of osteosarcoma by sponging the miR-212-3p/SGK3 axis. Cancer Cell Int. 2018;18:141.
Li Y, Guo D, Zhao Y, Ren M, Lu G, Wang Y, et al. Long non-coding RNA SNHG5 promotes human hepatocellular carcinoma progression by regulating miR-26a-5p/GSK3beta signal pathway. Cell Death Dis. 2018;9(9):888.
Ma Z, Xue S, Zeng B, Qiu D. lncRNA SNHG5 is associated with poor prognosis of bladder cancer and promotes bladder cancer cell proliferation through targeting p27. Oncol Lett. 2018;15(2):1924–30.
Zhao L, Han T, Li Y, Sun J, Zhang S, Liu Y, et al. The lncRNA SNHG5/miR-32 axis regulates gastric cancer cell proliferation and migration by targeting KLF4. FASEB J. 2017;31(3):893–903.
Jiang H, Huang G, Zhao N, Zhang T, Jiang M, He Y, et al. Long non-coding RNA TPT1-AS1 promotes cell growth and metastasis in cervical cancer via acting AS a sponge for miR-324-5p. J Exp Clin Cancer Res. 2018;37(1):169.
Billam M, Sobolewski MD, Davidson NE. Effects of a novel DNA methyltransferase inhibitor zebularine on human breast cancer cells. Breast Cancer Res Treat. 2010;120(3):581–92.
Cui M, Wen Z, Chen J, Yang Z, Zhang H. 5-Aza-2′-deoxycytidine is a potent inhibitor of DNA methyltransferase 3B and induces apoptosis in human endometrial cancer cell lines with the up-regulation of hMLH1. Med Oncol. 2010;27(2):278–85.
Kang KA, Kim HS, Kim DH, Hyun JW. The role of a ginseng saponin metabolite as a DNA methyltransferase inhibitor in colorectal cancer cells. Int J Oncol. 2013;43(1):228–36.
Qiu H, Yashiro M, Shinto O, Matsuzaki T, Hirakawa K. DNA methyltransferase inhibitor 5-aza-CdR enhances the radiosensitivity of gastric cancer cells. Cancer Sci. 2009;100(1):181–8.
Srivastava P, Paluch BE, Matsuzaki J, James SR, Collamat-Lai G, Taverna P, et al. Immunomodulatory action of the DNA methyltransferase inhibitor SGI-110 in epithelial ovarian cancer cells and xenografts. Epigenetics. 2015;10(3):237–46.
Pappano WN, Guo J, He Y, Ferguson D, Jagadeeswaran S, Osterling DJ, et al. The histone methyltransferase inhibitor A-366 uncovers a role for G9a/GLP in the epigenetics of leukemia. PLoS One. 2015;10(7):e0131716.
Kikuchi J, Takashina T, Kinoshita I, Kikuchi E, Shimizu Y, Sakakibara-Konishi J, et al. Epigenetic therapy with 3-deazaneplanocin a, an inhibitor of the histone methyltransferase EZH2, inhibits growth of non-small cell lung cancer cells. Lung Cancer. 2012;78(2):138–43.
The authors thank the study investigators and staff who participated in this study.
This work was supported by funding from the National Natural Science Foundation of China (Grant Number 81602323, 81872149), Wu Lien-teh Science Foundation of Harbin Medical University (Grant number WLD-QN1706), Distinguished Young Scholars of Harbin Medical University Cancer Hospital (Grant number JCQN2018–03), Yong Elite Training Foundation Grant of Harbin Medical University Cancer Hospital (Grant number JY2016–02), Natural Science Foundation of Heilongjiang Province (Grant number H201406) and Haiyan Fund Project of Harbin Medical University Cancer Hospital (Grant number JJQN 2018–10).
Availability of data and materials
The authors declare that the data supporting the findings of this study are available within the article and its supplementary information files.
Ethics approval and consent to participate
This study protocol conformed to clinical research guidelines and was approved by the research ethics committee of Harbin Medical University Cancer Hospital.
Consent for publication
The authors declare that they have no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Differential expression of lncRNAs between breast cancer tissues and normal tissues in HMUCC cohort. Table S2. Differential expression of lncRNAs between breast cancer tissues and normal tissues in TCGA cohort. Table S3. Overlapped differential lncRNAs between HMUCC cohort and TCGA cohort. Table S4. Restricted filter criterial (FC < 0.833, P < 0.001) identified 19 differential downregulated lncRNAs in HMUCC cohort and chromosome position of these lncRNAs. Table S5. Guilt-by-association followed by GO and KEGG analyses of TSLCRs in breast cancer. Table S6. Guilt-by-association followed by GO and KEGG analyses of EPB41L4A-AS2 in breast cancer. Table S7. Guilt-by-association analyses followed by GO and KEGG of EPB41L4A-AS2 in pan-cancer. Table S8. Differential gene expression between EPB41L4A-AS2 overexpression groups and control groups in UACC812 and BT549 breast cancer cell lines. (ZIP 3607 kb)
Figure S1. TSLNR expression in breast cancer samples and normal tissue samples in HMUCC. Figure S2. Genetic alteration was also examined for these lncRNAs in breast cancer data in TCGA. Figure S3. A&B Patients with high expression (N = 266) of TSLNRs (ACVR2B-AS1 and WEE2-AS1) had favorable OS than those with low expression (N = 266) in breast cancer in TCGA. C-F Patients with high expression (N = 266) of TSLNRs (ACVR2B-AS1, WEE2-AS1, LINC-PINT and HAND2-AS1) had favorable DFS than those with low expression (N = 266) in breast cancer in TCGA. Figure S4. A-D Patients with high expression (N = 266) of TSLNRs (CYP1B1-AS1, LINC-PINT, LINC00667 and GRIK1-AS1) had favorable OS than those with low expression (N = 266) in breast cancer in TCGA. E-G Patients with high expression (N = 266) of TSLNRs (CYP1B1-AS1, FAM66C and GRIK1-AS1) had favorable DFS than those with low expression (N = 266) in breast cancer in TCGA. Figure S5. EPB41L4A-AS2 was downregulated in MDA-MB-231 breast cancer cells with ZNF217 overexpression in GEO dataset GSE35511. Figure S6. A Overlapping genes of EPB41L4A-AS2 correlated genes and paclitaxel related genes in BETMAN-TCM. B KEGG pathway analysis for EPB41L4A-AS2 correlated genes in BETMAN-TCM. C GO analysis for EPB41L4A-AS2 correlated genes in BETMAN-TCM.D OMIM analysis for EPB41L4A-AS2 correlated genes in BETMAN-TCM.E Pharmacological network analysis indicates that EPB41L4A-AS2 may be involved in paclitaxel related process in breast cancer. F Pharmacological network analysis indicates that EPB41L4A-AS2 may be involved in crosstalk with paclitaxel related genes in breast cancer. Figure S7. A Expression of EPB41L4A-AS2 in breast cancer cell lines. B&C overexpression efficiency of EPB41L4A-AS2 in UACC812 and BT549 cells. D Knockdown efficiency of EPB41L4A-AS2 in MDA-MB-453 cells. Figure S8. A-C Overexpression of each lncRNA (MEG3, WEE2-AS1 and HAND2-AS1) inhibited clone formation in UACC812 cells. D-F Overexpression of each lncRNA (MEG3, WEE2-AS1 and HAND2-AS1) inhibited cell viability in UACC812 cells. (ZIP 19903 kb)