Abstract
Background/Aim: Several works in the past decades pointed out the key role of long intergenic non-coding RNA (lincRNA) in breast cancer development. Here in we report for first time the importance of deregulation of lincRNA RP11-400K9.4 in breast cancer cells which played a role in cell survival and migration. Materials and Methods: After RP11-400K9.4 silencing by short hairpin RNAs or overexpression by GeneBlocks, real-time quantitative polymerase chain reaction (RT-PCR), microarray, migration, proliferation and viability assay were performed. Results: RP11-400K9.4 expression was mainly in the cytoplasmic fraction in 2D culture. Overexpression of RP11-400K9.4 led to a reduction of migration by MCF-7 and MDA-MB-368 cells and an increase in cellular survival after UV-C radiation. Bioinformatic analyses highlighted irradiation-induced DNA damage, DNA repair and cell-cycle pathways as the mainly affected by RP11-400K9.4. Furthermore RT-PCR assay demonstrated the overexpression of baculoviral IAP repeat containing 3 (BIRC3) a known oncogene that promotes radiotherapy resistance through the nuclear factor kappa B (NFĸB) pathway. Conclusion: RP11-400K9.4 participates in the modulation of migration and survival processes probably via the BIRC3/NFĸB pathway.
Breast cancer is a worldwide health problem, representing the most common life-threatening malignant neoplasm in women and the second cause of death associated with cancer (1). In the case of Mexico, breast cancer is also the most common type in women, and its mortality has been increasing in recent years; this tendency is expected to continue in the coming years (2).
The cancer cell phenotype is a dynamic state in which cells modify distinctive and complementary capabilities such as proliferation, evasion of cell growth signals, sensitivity to death via different external and internal stimuli, migration and invasion (3).
Resistance to death constitutes a genuine hallmark of cancer and is a major obstacle in the treatment of cancer. The sensitivity of a cell to death from external and internal stimuli depends on specific molecules and transduction pathways regulated at the translational and post-translational levels. The inhibitors of apoptosis (IAP) are a family of functionally and structurally related proteins that serve as endogenous inhibitors of programmed cell death and play a very important role in cancer cells obtaining resistance to death (4). A common feature of IAPs is the presence of a baculovirus IAP repeat (BIR), a 70 amino acid domain, in one to three copies which function mainly by mediating protein–protein interactions (5). This family comprises BIRC1 (NAIP), BIRC2 (cIAP1), BIRC3 (cIAP2), BIRC4 (xIAP), BIRC5 (survivin) and BIRC6.
BIRC3 can directly bind to caspases 3, 7, and 9 to inhibit their activation (6). This protein has also been implicated in the activation of several signal transduction pathways associated with malignancy (7). BIRC3 was found to be up-regulated in breast cancer cells stimulated with interleukin-1β, and this change is associated with resistance to doxorubicin, whereas silencing of BIRC3 reduced the viability of cells after exposure to the drug (8). BIRC3 plays also an important role in regulating the nuclear factor-ĸB (NF-ĸB) signaling pathway, which is involved in the development of cancer (9). Recently, it was shown that this protein can be associated with non-coding RNAs: binding of miR-124 to BIRC3 inhibited the proliferation and migration of in human hepatocellular carcinoma (10). Interestingly, in 2019, Jiang et al. reported that the long non-coding RNA lncRNA-HCP5 promoted triple-negative breast cancer progression, regulating BIRC3 by binding to miR-219a-5p (11).
After the discovery of lncRNA in eukaryotic cells in 1986, work has shown their importance as genetic and epigenetic regulators in tumor cells and therefore in the acquisition of hallmarks or emergent neoplastic characteristics (12, 13). Various studies have shown a role of lncRNA in the regulation of migration and death resistance, such as CCND1-upstream intergenic DNA repair 1 (CUPPID1) and 2, nuclear paraspeckle assembly transcript 1 (NEAT1), and DNA damage-sensitive RNA1 (DDSR1) (14-17). Recently, long intergenic non-coding RNA (lincRNA) RP11400k9.4 was pointed out as a possible marker for triple-negative breast cancer (18-20). In addition, we found that this lincRNA is associated with deregulation of inhibitor of nuclear factor kappa B kinase subunit epsilon (IKKε) kinase and the stem cell phenotype in breast cancer (21). In this study, we aimed to elucidate the role of this lincRNA in breast cancer cells.
Materials and Methods
Cell culture. Cells lines MCF-7, T47D, ZR751, MDA-MB-361, MDA-MB-453, MDA-MB-468, MDA-MB-231, SKBR3, HS578T and BT-20 were obtained from the American Type Culture Collection (ATCC, Manassas, VA, USA) and maintained at 37°C with 5% of CO2 in Dulbecco's modified Eagle's medium (DMEM), Eagle's minimum essential medium (EMEM) or Roswell Park Memorial Institute (RPMI) medium with fetal bovine serum (FBS) at 5-15%, according to recommendations from the ATCC.
Cellular transfection. To generate cell lines that stably overexpress lincRNA RP11-400k9.4, MCF-7 and MDA-MB-468 cells were transfected with 7.5 ml of lipofectamine and 3 μg of the pcDNA3.1 vector containing the complete cDNA of the lincRNA GB400 (GeneBlock for RP11-400k9.4) (gBlocks ®Gene fragments, IDT, Newark, NJ, USA), As transfection control, we used the same vector with a Luciferase insert. These cells were selected with 600 and 800 μg/ml of 418 Geneticin (Thermo Fisher Scientific, CDMX, México) respectively. For the loss-of-function experiments, two short hairpin RNAs (shRNAs) were designated with the software RNAi Tool (Clontech, Takara, CA, USA). The oligonucleotides with the sequence of shRNA or a Luciferase sequence were inserted in the pSiren RQ vector (Clontech) and 3 μg of vector and 7.5 ml lipofectamine were used for the transfection. Cells were selected with 400 ng/ml puromycin (GIBCO, Gaithersburg, MD, USA).
Quantitative real-time polymerase chain reaction (RT-PCR). Trizol reagent was employed for total RNA isolation of transfected cells. 1 μg of RNA previously treated with DNAse was used for retrotranscription with High Capacity cDNA Reversion Transcription kit (Thermo Fisher Scientific, Franklin Lakes, MA, USA). cDNA (50 ng) was mixed with 7 μl of Sybr Select Master Mix (Thermo Fisher Scientific) and 0.24 mM of each primer Supplementary Table I (https://www.dropbox.com/sh/z82pwriywvwjw7f/AADPvyeuQTLKPixFjuyvSG0ba/Suplementary%20table%20I.docx?dl=0), and the assays were performed using a QuantStudio 7 Flex Real Time PCR system. Data were analyzed with QuantStudio Real-Time PCR software (Thermo Fisher Scientific, Franklin Lakes, MA, USA) according to the 2−ΔΔCT method (22).
Cellular fractionation. Nuclear and cytoplasmic fractions of MCF-7 cells were obtained with PARIS kit (Thermo Fisher Scientific) according to the manufacturer instructions. Total RNA and protein from both fractions were extracted. The purity of the fractions was corroborated with RNA electrophoresis. The protein fraction was separated with sodium dodecyl sulfate-polyacrylamide gel and transferred to a polyvinylidene fluoride membrane. Finally, the membrane was incubated with antibodies toward tubulin and A/C lamin (ab15246, ab224816; Abcam, Cambridge, UK) for cytoplasmic and nuclear fractions, respectively.
Extraction of total RNA. MCF-7 and MDA-MB-468 cells were lysed with 1 ml Trizol reagent, 200 μl chloroform was added and was then vortexed, after centrifuging at 8000 × g for 12 min. The aqueous phase was transferred to a fresh tube and was mixed with 0.5 ml isopropyl alcohol and then the bottom was recuperated by centrifugation at 10,000 × g for 10 min and dissolved with 40 μl of water.
RNA microarray and pathway analysis. The quality of total RNA was evaluated using an Agilent RNA Bioanalyzer and only RNA of samples with RNA integrity number (RIN) of 9 or higher were used. Three independent samples of MCFF7 cells overexpressing lincRNA RP11-400K9.4 and three of cells transfected with control luciferase were treated with DNAse and hybridized to Affymetrix Gen 2.0 Human Gene slides (Thermo Fisher Scientific) on INMEGEN Microarray Department. Analysis was performed with the Transcriptome Analysis console (Affymetrix/Thermo Fisher Scientific, Waltham, MA, USA). Only genes with a fold change threshold of ±3 or higher and p≤0.05 were considered for analysis with Key Pathway Advisor software (KPA; Clarivate Analytics, Philadelphia, PA, USA), METACORETM (Philadelphia, PA, USA). Gene Set Enrichment Analysis (GSEA; http://www.broadinstitute.org/gsea) was performed by applying Hallmark and KEGG gene sets downloaded from MSigDB database. We used default settings and considered a false discovery rate of <0.25 as statistically significant (23). Likewise, REVIGO (24) (http://revigo.irb.hr/) and Enrichr(25) (https://amp.pharm.mssm.edu/Enrichr/) analyses were performed with default settings (medium similarity=0.7 and calculated p-value for Gene Ontology predicted pathways). This allowed us to obtain a scatterplot of affected biological processes such as mammary gland development, apoptosis, irradiation-inducedDNA damage, by REVIGO and a table of main deregulated signaling pathways such as ATR by combining the BioPlanet 2019, NCI-Nature 2016 and WikiPathways 2019 Human databases using Enrichr software
Multicellular spheroid formation. A total of 2×106 MCF-7 cells were cultivated in non-adherent plates with Leibovitz's L-15 media and incubated with 60 rpm orbital agitation at 37°C. Spheroids were counted and measured at 21 days and finally disaggregated with trypsin. Total RNA was isolated as described above.
Migration assay. MCF-7 and MDA-MB-368 cells were seeded at 3×105 on inserts of transwell plates (Corning, NY, USA) without fetal bovine serum (FBS) overnight. The medium was replaced by fresh medium after 12 hours and medium with 10% FBS was placed on the bottom wells. After 24 hours, the cells were fixed with 4% paraformaldehyde and stained with Gentian violet at 10%. Finally, the insert membranes were allowed to dry and the cells that had migrated across the membrane were counted with a stereoscopic microscope and analyzed using ImageJ software (https://imagej.nih.gov/ij/index.html). Proliferation assay. MCF-7 cells were seeded 3×103 on each well of a 96-well plate (Corning). DMEM with 5% FBS was added and the cells maintained at standard temperature and CO2 conditions. The cells were then exposed to 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) reagent (Promega, Madison, WI, USA) for 30 minutes in the dark and the dye solubilized for 1 hour. The dye concentration was measured at 0, 24, 48, 72, 96, 120 and 144 hours using a Beckman DTX 880 spectrophotometer (Beckman, Indianapolis, IN, USA) at a wavelength of 595 nm.
Viability assay. MCF-7 cells transfected with Luc, shRNA1, shRNA2 and GB400 were irradiated with UV light at 254 nm for 3 minutes before seeding at 1.5×104/well in three wells of a 96-well plate (Corning). After 24 h, the cells were stained with MTT (Promega) reagent and absorbance at 595 nm measured.
Statistical analysis. Statistical analysis was performed with Graph Pad Prism 6 software (San Diego, CA, USA). We considered a minimum p-value of 0.05 and used one- and two-tailed ANOVA test and Student t-test to determine significant differences between groups.
Results
LincRNA-RP11400k9.4 general features. lincRNA-RP11400k9.4 (ENSG00000237807(RP11-400K9.4), ENST00 000426023.1 (RP11-400K9.4-001), lncOPRK1-3:12) is a long non-coding RNA of 2387 base pairs, intergenic, with two exons (Exon1: 154435889-54436491, Exon 2: 254427731-54429514), located on chromosome 8, at 54,427,731-54,436,491 (Supplementary Figure 1A; https://www.dropbox.com/sh/z82pwriywvwjw7f/AABVYIdJa_Psqbj-SL6zAHXFa/Supplemmentary%20figure%201.pdf?dl=0). It has been shown that lincRNA-RP11400k9.4 is expressed in four normal tissues/cell lines: Heart, testes, ovary and prostate; and in five types of cancer: Invasive breast carcinoma, lung squamous cell carcinoma and low grade glioma (http://annolnc.cbi.pku.edu.cn/status.jsp?id=5D3C6210-200201084703) (Supplementary Figure 1B and C; https://www.dropbox.com/sh/z82pwriywvwjw7f/AABVYIdJa_Psqbj-SL6zAHXFa/Supplemmentary%20figure%201.pdf?dl=0).
To determine the subcellular distribution of lincRNA-RP11400k9.4, cellular cytoplasmatic and nuclear fractions of MCF-7 cells were obtained and its levels were determined by RT-PCR in each fraction. As shown in Figure 1A, lincRNA-RP11400k9.4 was found at highest levels in the cytoplasmatic fraction. This result is coincident with the data on the lncATLAS database (https://lncatlas.crg.eu/) (Figure 1B).
To determine if lincRNA-RP11400k9.4 was characteristic of one of the breast cancer subtypes, expression analysis was performed on 10 cell lines. We found that the expression was not associated with the hormonal receptor status, subtype, claudin level or doubling time (Supplementary Table II; https://www.dropbox.com/sh/z82pwriywvwjw7f/AAAxgWg8nE155B8eNo4GKWR9a/Supplementary%20Tablas%20II.docx?dl=0). Interestingly, we found an inverse correlation between invasive capacity and level of lincRNA-RP11400k9.4. The lowest levels of expression were found in the MDA-MB-361 and HS578T and MDA231 cell lines, which have a high invasive capacity, whereas the highest levels were found in the MDA-MB-453, T47D and MCF-7 cell lines, which have lower invasive capacity (Figure 1C).
Comprehensive transcriptional analysis of MCF-7 overexpressing lincRNA-RP11400k9.4. In a previous study we observed a decrease in the level of several lncRNAs, including lincRNA-RP11400k9.4, in MCF-7 cells which had acquired a stem-like phenotype after IKKε was overexpressed (21). It is well known that three-dimensional (3D) MCF-7 culture (spheroids) mimic the main features of human solid tumors including migration, invasion, proliferation, stemness features, and an increase of IKKε level (21). We evaluated the expression of lincRNA-RP11400k9.4 in 3D cultures. As shown in Figure 2A, we found a lower expression of lincRNA-RP11400k9.4 in spheroids when compared to cells cultured in monolayers.
To deepen the understanding of the role of lincRNA-RP11400k9.4 in breast cancer we performed loss- and gain-of-function experiments by over- and underexpressing lincRNA-RP11400k9.4 in MCF-7 and MDA-MB-468 cells by transfecting two shRNAs (Figure 2B and D) or a plasmid with the lincRNA-RP11400k9.4 gene (Figure 2C and E); in both cases we used transfected cells with a control vector to compare the effects.
Firstly, we performed whole-genome expression analysis using MCF-7 cells with overexpression of lincRNA-RP11400k9.4 versus luciferase expression. The validation of data obtained with the microarray analysis was performed using RT-PCR (Figure 3). As shown in Figure 3B, 710 transcripts were deregulated: Levels of 398 increased and 312 decreased, using a fold change of 2 or more and p<0.05 as threshold. Table I shows the top 10 altered RNAs. In order to delve into the expression changes mediated by lincRNA-RP11400k9.4, we carried out a network analysis using the KPA, Metacore, Enrichr, GSEA and REVIGO tools. We found that the main altered processes were related to DNA repair mechanisms, including DNA damage by UV radiation (Figure 4).
To corroborate the role of lincRNA-RP11400k9.4 in DNA repair mechanisms, and in particular DNA damage induced by UV radiation, we exposed MCF-7 and MDA-MB 468 cells with dysregulated lincRNA-RP11400k9.4 levels to UV light (254 nm) for 3 min. Twenty-four hours later, cellular viability was measured. As shown in Figure 5A, cells with lincRNA-RP11400k9.4 knockdown were more sensitive to UV radiation exposure and, as expected, cells with lincRNA-RP11400k9.4 overexpression presented higher cellular survival after UV radiation. Since UV light induces apoptosis, we analyzed the RNA data for differentially expressed genes associated with this process. Interestingly, we found that BIRC3, BCL2 apoptosis regulator-like protein 11 (BIM) and survivin (BIRC5) were dysregulated by lincRNA-RP11400k9.4 overexpression. We then measured BIRC3 and BIRC5 expression using RT-PCR and found that the expression of both genes was increased (Figure 5B). It is important to note that in the microarray data only BIRC3 expression increased whereas that of BIRC5 decreased, perhaps due to the presence of multiple functional isoform of this proteins with different regulatory patterns.
Since it has been reported that BIRC3 increases resistance to cell death but at the same time inhibits migration and proliferation, we evaluated the migratory capacity of these cells by transwell assay. As shown in Figure 6, both MCF-7 and MDA-MB-468 cells with knockdown of lincRNA-RP11400k9.4 expression presented a migration increase. As expected, overexpression of lincRNA-RP11400k9.4 produced a reduction in the migration of MCF-7 cells (Figure 6A-D). Finally, we measured the proliferative capacity using transfected MCF-7 and MDA-MB-468 cells. As shown in Figure 6E, changes in the level of lincRNA-RP11400k9.4 did not affect the proliferation of these cells.
Discussion
Cancer consists of at least 100 diseases shaped by an uncontrolled cell growth. Breast cancer in particular, represent the third cause of Mexican cancer-associated deaths for women of 30 years and older and has maintain this status over the past 8 years (1, 26).
Twenty-four years ago, with the discovery of the first lncRNA in eukaryotic cells and their important role in cellular physiology, a new paradigm was established and more recently, this paradigm has extended to carcinogenesis. Recently, Mathias et al. found that some lncRNAs are potential markers of breast cancer subtypes (18). More interesting, Su et al. performed bioinformatical analyses on 658 ductal carcinomas and found several lncRNAs to be differentially expressed between breast cancer subtypes, lincRNA-RP11400k9.4 being one of the RNAs associated with human epidermal growth factor receptor 2 (HER2)-overexpressing subtype (19). HER2-overexpressing tumors are characterized by their aggressiveness and short survival time for their carrier. This observation was similar to that of Van Grembergen et al. in a 995 patient cohort where lincRNA-RP11400k9.4 was also overexpressed in the HER2 subtype (20). We did not find a clear correlation between the level of lincRNA-RP11400k9.4 and HER2-positive status in a panel of breast cancer cell lines. A higher level of this lincRNA was found in the MDA-MB-453 cell line, which is HER2-positive, whereas a lower level was found in MDA-MB-361 cells, which are also HER2-positive. In addition, the T47D cell line, which is negative for HER2 showed high levels of lincRNA-RP11400k9.4. These results indicate that alternative pathways may be involved in regulating lincRNA-RP11400k9.4 expression.
LincRNA-RP11400k9.4 was found to be mainly located in the cytoplastic compartment and was dysregulated in 3D breast cancer cultures and in breast cancer cells that overexpress IKKε. When lincRNA-RP11400k9.4 was overexpressed in breast cancer cells, a relatively large number of RNAs were deregulated. Network analyses with data from dysregulated RNAs suggested that lincRNA-RP11400k9.4 may be associated with the ability of cells to survive when exposed to UV radiation. These results were validated with in vitro loss- and gain-of-function approach assays. As expected, overexpression of lincRNA-RP11400k9.4 resulted in a decrease in cell death caused by UV exposure and a decrease in migration but no changes in cell proliferation. BIRC3, a protein that is a member of the NFkappa-B survival pathway, was deregulated in cells that overexpressed lincRNA-RP11400k9.4.
BIRC3 participates in non-canonical nuclear factor kappa B (NFĸB) pathway activation. The activation of NFĸB is a commonly observed phenomenon in breast cancer, facilitating the development of a hormone-independent, invasive, high-grade and late-stage tumor phenotype (27). Activated NFĸB is associated with resistance to chemotherapy, radiotherapy and endocrine therapy (28).
BIRCs are able to abrogate death-inducing signals, by direct inhibition of caspases and modulation on NFĸB (16). Aberrant BIRC3 expression has been found in a several human cancer types. Dierlamm et al. found that the fusion BIRC3–mucosa-associated lymphoid tissue lymphoma protein in lymphoma conferred resistance to apoptosis (29). Yuan et al. reported that the BIRC3 enhanced resistance to apoptosis in HPV16E6/E7-immortalized human oral keratinocytes (30). There are reports that point to the role of BIRC3 in the regulation of UV-induced apoptosis: Ismail and Bateman reported that BIRC3 is regulated by T-box transcription factor 2 in SW13 cells and increased resistance to apoptotic stimuli including UV-radiation (31) and it has been reported that BIRC3 was overexpressed in UV- irradiated hepatocytes (32). BIRC3 might mediate the effect of lincRNA-RP11400k9.4,lincRN, perhaps operating downstream of NFĸB and be partly responsible for its effects. Further experiments are needed to address this.
Conclusion
LincRNA-RP11400k9.4 is an intergenic non-coding RNA, expressed in the cytoplasm of breast cancer cells, which regulates the sensitivity to death after exposure to UV radiation. This lincRNA also inhibits cell migration, although not proliferation in these cells. Its effects may be mediated by BIRC3, a member of the family of IAPs that is regulated by NFĸB.
Acknowledgements
This work was part of the PhD thesis of Miguel Ángel Fernández Rojas, doctoral student from the Programa de Doctorado en Ciencias Biomédicas, Universidad Nacional Autónoma de México (UNAM) and has received CONACyT fellowship 288752. This study was supported by grants Grant A1-S-8462 to J.M-Z. We greatly acknowledge the technical support of Dra. Karla Vázquez Santillán, Dra. Rosario Castro Oropeza, Dra. Marcela De la Fuente Hernández and Dr. Alfredo García Venzor for the ongoing support.
Footnotes
Authors' Contributions
M.F., J.M., and V.M. contributed to the design and implementation of the research, to the analysis of the results and to the writing of the article.
This article is freely accessible online.
Conflicts of Interest
The Authors declare that there are no conflicts of interest in regard to this study.
- Received July 28, 2020.
- Revision received August 29, 2020.
- Accepted September 2, 2020.
- Copyright© 2020, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved