New insight into long non-coding RNAs associated with bone metastasis of breast cancer based on an integrated analysis

doi:10.21203/rs.3.rs-391903/v1

Download PDF

Primary research

New insight into long non-coding RNAs associated with bone metastasis of breast cancer based on an integrated analysis

https://doi.org/10.21203/rs.3.rs-391903/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 13 Jul, 2021

Read the published version in Cancer Cell International →

You are reading this latest preprint version

Background

Bone is the most common site of metastatic breast cancer, and it is a leading cause of breast cancer-related death. This study aimed to explore bone metastasis-related long non-coding RNAs (lncRNAs) in breast cancer.

Methods

Four mRNA datasets and two lncRNA datasets of bone metastasis, lung metastasis and liver metastasis of breast cancer were downloaded from Gene Expression Omnibus (GEO) database. The differentially expressed mRNAs (DEmRNAs) and lncRNAs (DElncRNAs) in group of bone metastasis vs lung metastasis and bone metastasis vs liver metastasis, as well as the overlap of the two groups, were identified. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and protein-protein interaction (PPI) network construction of DEmRNAs were conducted. The cis nearby-targeted DEmRNAs of DElncRNAs were obtained. Quantitative real-time polymerase chain reactions (qRT-PCR) was used to detect the expression levels of selected DEmRNAs and DElncRNAs. LOC641518- lymphoid enhancer-binding factor 1 (LEF1) pair was selected to verify its role in migration and invasion capability of breast cancer cells by wounding healing assay and transwell invasion assay.

Results

A total of 237 DEmRNAs were obtained in bone metastasis compared with both lung metastasis and liver metastasis. A total of 3 DElncRNAs in bone metastasis compared with both lung metastasis and liver metastasis were obtained. A total of 7 DElncRNA-nearby-targeted DEmRNA pairs and 15 DElncRNA-nearby-targeted DEmRNA pairs in group of bone metastasis vs lung metastasis and bone metastasis vs liver metastasis, were detected, respectively. Four cis LncRNA-mRNA interaction pairs were identified, which are LOC641518-LEF1, FLJ35024-Very Low Density Lipoprotein Receptor (VLDLR), LOC285972-Retinoic Acid Receptor Responder 2 (RARRES2) and LOC254896- TNF receptor superfamily member 10c (TNFRSF10C). qRT-PCR using clinical samples from our hospital confirms the bioinformatics prediction. siRNA knocking down LOC641518 down-regulates LEF1 mRNA expression, and reduces the migration and invasion capability of breast cancer cells.

Conclusions

We concluded that four LncRNA-mRNA pairs, including LEF1-AS1-LEF1, may play a central role in breast cancer bone metastasis.

Cancer Biology

breast cancer

bone metastasis

lncRNA

GEO

interaction pairs

LOC641518

LEF1

Globally, breast cancer is the most common malignancy in women [1]. Bone is the most common site of metastatic breast cancer, accounting for nearly three-quarters of patients with metastatic breast cancer [2]. Bone metastasis is often accompanied by bone pain and skeletal-related events, leading to increased mortality and seriously affecting the quality of life of patients [3]. The clinical complications of breast cancer bone metastasis bring heavy burden to the individual and society [4]. Bone metastases secondary to breast cancer are incurable, and the molecular mechanisms are complex and involved in multiple process [4].

Over the past few decades, a great deal of information about lncRNAs has been generated, which facilitated the studies on diverse cancer etiology, including breast cancer [5]. Increasing evidence demonstrated that lncRNAs participate in many cellular process, such as cell fate decision, immune response, cancer cells proliferation and metastasis [6]. Gooding et al. indicated that the lncRNA BORG drives breast cancer metastasis and disease recurrence [7]. Metastasis associated lung adenocarcinoma transcript 1 (MALAT1), previously described as a metastasis-promoting lncRNA, was reported to suppress breast cancer metastasis [8]. HOXA11-AS was demonstrated to promote breast cancer invasion and metastasis through affecting epithelial-mesenchymal transition [9]. However, there are few studies on bone metastasis-related lncRNAs in breast cancer, so it is necessary to explore important lncRNAs associated with breast cancer bone metastasis.

In this present study, a comprehensive integrated analysis of lncRNAs and mRNAs expression profiles of bone metastasis, lung metastasis and liver metastasis of breast cancer downloaded from GEO database were performed. The differentially expressed mRNAs (DEmRNAs) and lncRNAs (DElncRNAs) were acquired. The cis nearby-targeted DEmRNAs of DElncRNAs were obtained as well. Finally, we established in vitro cell based assay to confirm the role of novel LncRNA-mRNA axis in tumor metastasis. This study aimed to provide additional knowledge on the molecular mechanisms of breast cancer bone metastasis, especially the role of pro-bone metastasis LncRNAs.

Dataset collection

To acquire the mRNA and lncRNA expression profiles of bone metastasis, lung metastasis and liver metastasis of breast cancer tissues, datasets in GEO database with the following criteria were retrieved: datasets should be whole-genome mRNA/lncRNA expression profile by array; these data were derived from bone metastasis, lung metastasis and liver metastasis of breast cancer tissues; datasets were normalized or original. Four mRNA datasets and two lncRNA datasets were enrolled in this study and shown in Table 1.

Differential expression analysis

MetaMA and limma package in R were utilized to acquire the DEmRNAs and DElncRNAs. P-values were calculated using a significance threshold of p-value < 0.05. With R package “pheatmap”, hierarchical clustering analysis of DEmRNAs and DElncRNAs were performed.

GO and pathway analysis of DEmRNAs

GeneCodis3 (http://genecodis.cnb.csic.es/analysis) was applied to perform GO and KEGG pathway enrichment analysis of DEmRNAs. The threshold was set as false discover rate (FDR) < 0.05.

Protein-protein interaction (PPI) network construction

Top 50 up- and down-regulated DEmRNAs were searched with the BioGrid (http://www.uniprot.org/database/DB-0184), and PPI network was constructed with Cytoscape software (version 3.6.1, http://www.cytoscape.org).

Cis nearby-targeted DEmRNAs of the DElncRNAs

DEmRNAs transcribed within a 100-kb window upstream or downstream of DElncRNAs were searched, which were defined as cis nearby-targeted DEmRNAs of DElncRNAs, to obtain the targeted DEmRNAs of DElncRNAs with cis-regulatory effects.

Quantitative real-time polymerase chain reaction (qRT-PCR) validation

Twenty-five patients with bone metastasis of breast cancer, 25 patients with lung metastasis of breast cancer, and 25 patients with liver metastasis of breast cancer were incorporated in our study. Meanwhile, 25 primary breast cancer patients without metastasis were enrolled, and their adjacent non-cancer breast tissues was considered as normal tissue to be used as control in the present study. Their tumor tissues from metastatic sites were used for RNA isolation. Ethical approval was obtained from the ethics committee of the First People’s Hospital of Chengdu and informed written consent was obtained from all of subjects. All the patients in the current study were same as described in the previous publication [10].

Total RNA of the tumor tissues from metastatic sites was extracted using the RNA liquid Reagent (TIANGEN Bio, Beijing, China) according to the manufacturer’s protocols. One microgram RNA was applied to synthesize cDNA by Fast Quant Reverse Transcriptase (TIANGEN Bio). Then real-time PCR was performed in an ABI 7300 Real-time PCR system with SYBR R Green PCR MasterMix (TIANGEN). All reactions were carried out in triplicate. GAPDH was used for the internal reference. Relative gene expression was analyzed by fold change method. And all the primer sequence was shown in supplemental table 1.

Cell culture and transfection with plasmid and siRNA

Human breast cancer cell line MCF-7 and MDA-MB-231, and non-carcinoma human breast epithelial cell line MCF-10A were purchased from ATCC (Manassas, VA, USA) and were cultured in RPMI 1640 medium (Invitrogen, CA, USA) supplemented with 10% fetal bovine serum (Invitrogen) and incubated at 37◦C and 5% CO₂. Commercial RNA interfering sets targeting LOC641518 were purchased from Ribobio technology (Shenzhen, Guangdong, China), and scramble siRNA was also transfected into breast cancer cells to regard it as control. Exogenous human full-length LEF1 plasmid vectors were purchased from Sino biological Inc. (Beijing, China); and empty vector was also transfected as negative control. Lipofectamine 3000 (Invitrogen) was used for all transfection studies by following the manufacturer’s protocol. Transfection efficiency was examined using qRT-PCR. Briefly, total RNA was extracted from cells using Trizol (Invitrogen, CA, USA) and further purified by RNAeasy mini kit plus DNase I treatment (Qiagen, Germany). The relative mRNA level for each gene was quantified by real-time RT-PCR with SYBR Green (Applied Biosystems, CA, USA), using GAPDH as a control.

Wound healing assay

The migration ability of cells was measured by wound-healing assay. Briefly, 2 × 10⁴ cells were inoculated in 6-cm tissue culture dishes, cultured overnight, scratch was performed when cell growth fusion reached to 90%. Migration images were captured 24 and 48 h after scratching, the migration rate (width) of cells were calculated.

Transwell migration

Corning^TM BioCoat^TM Matrigel^TM Invasion Chamber with Corning^TM Matrigel Matrix (Thermo Fisher Scientific, Waltham, MA, USA) was used for cell invasion assay. One hundred

microliters of cells with concentration of 5 × 10⁵/mL were added in the upper chamber with (8μm pore), 600 μL medium containing 10% serum was added in the lower chamber. After cultured for 6 h, medium in the chamber were removed and unmigrated cells were swabbed. Cells was fixed by 4% polyformaldehyde for 10min, stained with crystal violet for another 10min. The filter membrane was photographed under the inverted microscope (200×) after sealed with the neutral gum. Cells were counted by Image Pro Plus Version 6, 3 wells in each group and 5 vision fields of each well were randomly selected to calculate the average number of cells.

Identification of DEmRNAs in tumor tissues from bone metastatic sites compared with other metastatic sites of breast cancer

Compared with lung metastatic sites, a total of 1280 DEmRNAs (645 up- and 635 down-regulated mRNAs) were detected in tumor tissues from bone metastatic sites of breast cancer. Compared with liver metastatic sites, a total of 1482 DEmRNAs (851 up- and 631 down-regulated mRNAs) were detected in tumor tissues from bone metastatic sites of breast cancer. Hierarchical clustering analysis of top 100 up- and down-regulated DEmRNAs was exhibited in Figure 1A and Figure 1B, respectively. Top 10 up- and down-regulated DEmRNAs were displayed in Table 2. After overlapped these 1280 DEmRNAs and 1482 DEmRNAs, a total of 237 DEmRNAs (149 up- and 88 down-regulated mRNAs), which were differentially expressed in bone metastasis compared with both lung metastasis and liver metastasis, were obtained (Figure 1C-D, supplemental table 2 and table 3).

GO and KEGG pathway analysis of DEmRNAs

GO analysis indicated that DEmRNAs in group of bone metastasis vs lung metastasis were significantly enriched in negative regulation of cell proliferation (GO: 0008285; FDR = 5.43E-15), cytoplasm (GO: 0005737; FDR = 3.82E-65) and protein binding (GO: 0005515; FDR = 9.71E-76) (Figure 2A); DEmRNAs in bone metastasis vs liver metastasis group were significantly enriched in blood coagulation (GO: 0007596; FDR = 2.81E-26), extracellular region (GO: 0005576; FDR = 5.82E-104) and protein binding (GO: 0005515; FDR = 4.33E-51) (Figure 2B); the common DEmRNAs in two group mentioned above were significantly enriched in cell adhesion (GO: 0007155; FDR = 2.35E-07), extracellular region (GO: 0005576; FDR = 6.22E-22) and protein binding (GO: 0005515; FDR = 3.35E-08) (Figure 2C), respectively. All the top 10 GO list between these two groups were shown in supplemental table 4.

KEGG pathway enrichment analysis indicated that DEmRNAs in group of bone metastasis vs lung metastasis were significantly enriched in Focal adhesion (FDR = 3.08E-06), Pathways in cancer (FDR = 1.22E-05) and ECM-receptor interaction (FDR = 1.73E-05) (Figure 3A); DEmRNAs in bone metastasis vs liver metastasis group were significantly enriched in complement and coagulation cascades (FDR = 3.77E-41), Drug metabolism-cytochrome P450 (FDR = 2.63E-16) and Retinol metabolism (FDR = 1.07E-14) (Figure 3B); the common DEmRNAs in two group mentioned above were significantly enriched in Focal adhesion (FDR = 7.32E-06), ECM-receptor interaction (FDR = 1.39E-05) and Fatty acid metabolism (FDR = 3.61E-04) (Figure 3C), respectively.

PPI network construction

The PPI network constructed with DEmRNAs in group of bone metastasis vs lung metastasis included 147 nodes and 121 edges. TRIM24 (degree = 9), FGB (degree = 8) and FAM189A2 (degree = 7) were three hub proteins of the PPI network (Figure 4A). The PPI network constructed with DEmRNAs in group of bone metastasis vs liver metastasis included 231 nodes and 228 edges. TNFAIP6 (degree = 26), ACAN (degree = 21) and ADH1C (degree = 14) were three hub proteins of the PPI network (Figure 4B).

Identification of DElncRNAs in tumor tissues from bone metastatic sites compared with other metastatic sites of breast cancer

Compared with lung metastatic sites, a total of 71 DElncRNAs (38 up- and 33 down-regulated lncRNAs) were detected in tumor tissues from bone metastatic sites of breast cancer (Supplemental table 5). Compared with liver metastatic sites, a total of 91 DElncRNAs (40 up- and 51 down-regulated lncRNAs) were detected in tumor tissues from bone metastatic sites of breast cancer (Supplemental table 6). Hierarchical clustering analysis of DElncRNAs was exhibited in Figure 5A and Figure 5B, respectively. Top 10 up- and down-regulated DElncRNAs were displayed in Table 2. After overlapped these 71 DElncRNAs and 91 DElncRNAs, a total of 3 DElncRNAs (2 up- and 1 down-regulated lncRNAs), which were differentially expressed in bone metastasis compared with both lung metastasis and liver metastasis, were obtained (Figure 5C-D, table 3).

Correlation study of DElncRNAs and DEmRNAs

The correlation study between DELncRNA and DEmRNA expression was analyzed by Person’s test, the correlation coefficient and p-value were shown in table 4 (bone metastasis vs lung metastasis and bone metastasis vs liver metastasis), only pairs whose P value was less than 0.05 were shown.

Cis nearby-targeted DEmRNAs of the DElncRNAs

A total of 7 DElncRNA-nearby-targeted DEmRNA pairs in group of bone metastasis vs lung metastasis, involving in 7 DElncRNAs and 7 DEmRNAs, were detected (Table 5). A total of 15 DElncRNA-nearby-targeted DEmRNA pairs in group of bone metastasis vs liver metastasis, involving in 12 DElncRNAs and 15 DEmRNAs, were detected (Table 5). Combined with correlation study, total 8 pairs were selected with significant correlation, and four pairs had correlation coefficient higher than 0.7, which were LOC641518- LEF1, FLJ35024- VLDLR, LOC285972- RARRES2 and LOC254896- TNFRSF10C respectively.

qRT-PCR validation

These four pairs whose correlation coefficient was above 0.7 were selected for verification by qRT-PCR using clinical samples from our hospital. As shown in Fig.6A, results showed that the relative expression of LEF1, VLDLR, RARRES2 and TNFRSF10C were significantly up-regulated compared to adjacent non-tumor tissues and other metastatic sites (P <0.05). As shown in Fig.6B, RNA levels of LncRNA LOC641518 and FLJ35024 showed same trend with bioinformatics data. Although RNA levels of LOC285972 and LOC254896 in the bone metastasis were also higher than adjacent non-tumor tissues, they were not different compared with lung and liver metastasis. All the mRNA and LncRNA levels of tested genes from adjacent non-tumor breast cancer was significantly lower than metastasis tumor tissues (Fig. 6). After that, we used the RNA relative levels of each pair LncRNA-mRNA from metastatic bone tumor tissues to perform a correlation analysis, which was shown in Fig6.C. Results demonstrated all the four LncRNA pairs had significantly positive correlation, although the coefficient is not high as bioinformatics prediction. The validation test by qRCR proved that most of bioinformatics predictions were reliable.

LOC641518 Promotes breast cancer Metastasis via Activation LEF1

Due to LEF1 has been long considered as a metastasis mediator [11, 12], in the current study, the Loc641518-LEF1 axis was selected to investigate whether Loc641518 regulates LEF1 expression to promote breast cancer metastasis by in vitro study, which can help to confirm the speculation from bioinformatics data mining.

The RNA levels of LOC641518 and LEF1 in non-carcinoma human breast epithelial cell line MCF-10A and two breast cancer cell lines, MCF-7 and MDA-MB-231, were examined by qRT-PCR (Fig.7A and 7B), which could find that RNA levels of LOC641518 and LEF1 is higher in cancer cells than non-carcinoma cells, and highest in MDA-MB-231.

siRNA was used to knock down the expression of LOC641518 in both MCF7 and MDA-MB-231 cells. As shown in Fig.7C, the expression of LOC641518 was significantly reduced by siRNA compared with scramble siRNA control. Meanwhile, the LEF1 mRNA level was also reduced in siRNA-LOC641518 groups compared with scramble groups (Fig.7D). The wound healing assay and migration assay both demonstrated that knocking down LOC641518 significantly reduced capability of migration and invasion (Fig.7E and 7F), and exogenous overexpression of LEF1 could reverse these phenomena.

Bone is a primary site for breast cancer metastases, and patients often suffer from osteolytic bone metastases, which is incurable [13]. Previously, we conducted integrated bioinformatics analysis to find several possible key genes which are related with bone metastasis [10]. Currently, we used same gene dataset, expanding to find differential LncRNAs which are associated with bone metastasis, and further elucidate more possible LncRNA-mRNA interaction pairs. The purpose of this study is to introduce a broaden vision that more complicated gene networks are involved into the bone metastasis in breast cancer.

The most merit of the current study is that 71 DElncRNAs (bone metastatic vs lung metastasis) and 91 DElncRNAs (bone metastasis vs liver metastasis) have been identified. By correlation study between these DElncRNAs and DEmRNAs, combined with cis nearby-targeted analysis, total 9 LncRNA-mRNA interaction pairs were identified. By qRT-PCR study using clinical samples from our hospital, we confirm that most predictions are accurate and reliable. All these enriched data provide a landscape for researcher to clarify the role of LncRNA in breast cancer bone metastasis.

One of limitation is that we only verify the role of LOC641518-LEF1 interaction pair in migration and invasion capability of breast cancer cells, and the remaining still deserves further investigation. By knocking down the LOC641518, we observed reduced migration and invasion ability consistent with lower mRNA levels of LEF1, and over-expression reverse its effect.

LEF1, a member of the T cell Factor (TCF)/LEF1 family of high-mobility group transcription factors, is a downstream component of the Wnt signaling pathway [14]. LEF1 and HOXB9, two Wnt target genes, are identified as promoters of lung adenocarcinoma metastasis, and WNT/TCF pathway was identified as a determinant of metastasis to bone in lung adenocarcinoma [11]. In addition, it was reported that LEF1 was implicated in cell invasion and metastasis of breast cancer [15, 16]. LEF1 antisense RNA 1 (LEF1-AS1/LOC641518) is an antisense lncRNA encoded in the LEF1 locus, which was related to the metastasis in various tumors, including colorectal cancer and esophageal squamous cell carcinoma [17, 18]. Current study identified that LOC641518 regulated the expression of LEF1 in cis in group of bone metastasis vs lung metastasis, which may indicate that LOC641518 involved in breast cancer bone metastasis through regulating LEF1.

LOC285972-RARRES2, is another interaction pair identified in the current study. RARRES2, also known as chemerin, was considered as a chemoattractant and an adipokine [19]. The correlation between RARRES2 and cancers has been reported. Reduced expression of RARRES2 was observed in human hepatocellular carcinoma [20]. Lu et al. indicated RARRES2 that promotes tumorigenesis and metastasis in oral squamous cell carcinoma [21]. In addition, RARRES2 was identified as a potential biomarker for chordoma which is a rare, low-malignant bone tumor [22]. Low level of LINC00996 was report to be associated with colorectal carcinogenesis and metastasis [23]. In this study, LOC285972 was identified to down-regulated express in both two groups and regulate RARRES2 in cis in the group of bone metastasis vs lung metastasis. Hence, we hypothesized that LOC285972-RARRES2 interaction pair was associated with the development of breast cancer bone metastasis.

Very low density lipoprotein receptor (VLDLR) has been considered as a multiple function receptor due to binding numerous ligands, causing endocytosis and regulating cellular signaling [24, 25]. Previous studies have confirmed that high VLDLR expression is correlated with lymph node and distant metastasis in gastric and breast cancer patients, that VLDLR may be a clinical marker in cancers, and has a potential link with β-catenin signaling pathway [24]. In the current study, we predict that FLJ35024 regulate VLDLR, but very rare report about FLJ35024 so far. Therefore, their interaction deserves further study.

TNFRSF10C copy number variation has been reported that is associated with metastatic colorectal cancer [26], and its hypermethylation was associated with NSCLC [27]. One publication reports that TNFRSF10C plays a role in triple-negative breast cancer cell survival and metastasis regulated by SOX9 [28]. LOC254896 was expected to cis regulate its expression by bioinformatics analysis, but LOC254896 is still mysterious and lack of any detailed study.

In conclusion, four interaction pairs were identified to closely associate to breast cancer bone metastasis in this analysis. These findings may make contribution to develop new therapies that effectively target metastatic tumors and understand of the process of metastasis at biological and molecular levels.

Ethics approval and consent to participate

Ethical approval was obtained from the ethics committee of the First People’s Hospital of Chengdu and informed written consent was obtained from all of subjects.

Consent for publication

All the authors agree to publish this manuscript and all the data.

Availability of data and materials

All the data are available on requirement.

Competing interests

There is no competing interesting to be declared.

Funding

This project was supported by the Drug Innovation Major Project of China [2018ZX09711001-002-010]; Chinese Academy of Medical Sciences Innovation Fund for Medical Sciences [2016-I2M-3–011]; and Beijing Natural Science Foundation [7202138].

Authors' contributions

Yu zhang and Xiaofeng Huang are in charge data mining and manuscript draft, Jin Liu are in charge of cell culture and in vitro assay, Guo Chen and Chengjun Liu are in charge of data analysis; Sen Zhang is in charge of data mining; Jiaxin Li is in charge of project design.

Acknowledgements

Not available.

lncRNAs: long non-coding RNAs; qRT-PCR: Quantitative real-time polymerase chain reaction; GEO: Gene Expression Omnibus; DEmRNAs: differentially expressed mRNAs; DELncRNAs: differentially expressed LncRNAs; KEGG: Kyoto Encyclopedia of Genes and Genomes; GO: gene ontology; LEF1: Lymphoid enhancer-binding factor 1; VLDLR: very low density lipoprotein receptor. RARRES2: Retinoic Acid Receptor Responder 2; TNFRSF10C: TNF receptor superfamily member 10c; MALAT1: Metastasis associated lung adenocarcinoma transcript 1; PPI: Protein-protein interaction; FDR: false discover rate.

Zulauf N, Bruggmann D, Groneberg D, Oremek GM: Expressiveness of Bone Markers in Breast Cancer with Bone Metastases. Oncology 2019, 97(4):236-244.
Holen I, Lefley DV, Francis SE, Rennicks S, Bradbury S, Coleman RE, Ottewell P: IL-1 drives breast cancer growth and bone metastasis in vivo. Oncotarget 2016, 7(46):75571-75584.
Parkes A, Warneke CL, Clifton K, Al-Awadhi A, Oke O, Pestana RC, Alhalabi O, Litton JK, Hortobagyi GN: Prognostic Factors in Patients with Metastatic Breast Cancer with Bone-Only Metastases. The oncologist 2018, 23(11):1282-1288.
Brook N, Brook E, Dharmarajan A, Dass CR, Chan A: Breast cancer bone metastases: pathogenesis and therapeutic targets. The international journal of biochemistry & cell biology 2018, 96:63-78.
Wang Y, Wu N, Liu J, Wu Z, Dong D: FusionCancer: a database of cancer fusion genes derived from RNA-seq data. Diagnostic pathology 2015, 10:131.
Chen Y, Yu X, Xu Y, Shen H: Identification of dysregulated lncRNAs profiling and metastasis-associated lncRNAs in colorectal cancer by genome-wide analysis. Cancer medicine 2017, 6(10):2321-2330.
Gooding AJ, Zhang B, Jahanbani FK, Gilmore HL, Chang JC, Valadkhan S, Schiemann WP: The lncRNA BORG Drives Breast Cancer Metastasis and Disease Recurrence. Scientific reports 2017, 7(1):12698.
Kim J, Piao HL, Kim BJ, Yao F, Han Z, Wang Y, Xiao Z, Siverly AN, Lawhon SE, Ton BN et al: Long noncoding RNA MALAT1 suppresses breast cancer metastasis. Nature genetics 2018, 50(12):1705-1715.
Li W, Jia G, Qu Y, Du Q, Liu B, Liu B: Long Non-Coding RNA (LncRNA) HOXA11-AS Promotes Breast Cancer Invasion and Metastasis by Regulating Epithelial-Mesenchymal Transition. Medical science monitor : international medical journal of experimental and clinical research 2017, 23:3393-3403.
Zhang Y, He W, Zhang S: Seeking for Correlative Genes and Signaling Pathways With Bone Metastasis From Breast Cancer by Integrated Analysis. Frontiers in Oncology 2019, 9(138).
Nguyen DX, Chiang AC, Zhang XH, Kim JY, Kris MG, Ladanyi M, Gerald WL, Massague J: WNT/TCF signaling through LEF1 and HOXB9 mediates lung adenocarcinoma metastasis. Cell 2009, 138(1):51-62.
Min RQ, Ma Q: MicroRNA-381 inhibits metastasis and epithelial-mesenchymal transition of glioblastoma cells through targeting LEF1. Eur Rev Med Pharmacol Sci 2020, 24(12):6825-6833.
Taipaleenmaki H, Farina NH, van Wijnen AJ, Stein JL, Hesse E, Stein GS, Lian JB: Antagonizing miR-218-5p attenuates Wnt signaling and reduces metastatic bone disease of triple negative breast cancer cells. Oncotarget 2016, 7(48):79032-79046.
Cadigan KM, Waterman ML: TCF/LEFs and Wnt signaling in the nucleus. Cold Spring Harbor perspectives in biology 2012, 4(11).
Nguyen A, Rosner A, Milovanovic T, Hope C, Planutis K, Saha B, Chaiwun B, Lin F, Imam SA, Marsh JL et al: Wnt pathway component LEF1 mediates tumor cell invasion and is expressed in human and murine breast cancers lacking ErbB2 (her-2/neu) overexpression. International journal of oncology 2005, 27(4):949-956.
Zheng T, Wang A, Hu D, Wang Y: Molecular mechanisms of breast cancer metastasis by gene expression profile analysis. Molecular medicine reports 2017, 16(4):4671-4677.
Shi Q, He Y, Zhang X, Li J, Cui G, Zhang X, Wang X: Two Novel Long Noncoding RNAs - RP11-296E3.2 and LEF1-AS1can - Separately Serve as Diagnostic and Prognostic Bio-Markers of Metastasis in Colorectal Cancer. Medical science monitor : international medical journal of experimental and clinical research 2019, 25:7042-7051.
Zong MZ, Feng WT, Du N, Yu XJ, Yu WY: Upregulation of long noncoding RNA LEF1-AS1 predicts a poor prognosis in patients with esophageal squamous cell carcinoma. European review for medical and pharmacological sciences 2019, 23(18):7929-7934.
Ernst MC, Sinal CJ: Chemerin: at the crossroads of inflammation and obesity. Trends in endocrinology and metabolism: TEM 2010, 21(11):660-667.
Lin W, Chen YL, Jiang L, Chen JK: Reduced expression of chemerin is associated with a poor prognosis and a lowed infiltration of both dendritic cells and natural killer cells in human hepatocellular carcinoma. Clinical laboratory 2011, 57(11-12):879-885.
Lu Z, Liang J, He Q, Wan Q, Hou J, Lian K, Wang A: The serum biomarker chemerin promotes tumorigenesis and metastasis in oral squamous cell carcinoma. Clinical science (London, England : 1979) 2019, 133(5):681-695.
Scheil-Bertram S, Kappler R, von Baer A, Hartwig E, Sarkar M, Serra M, Bruderlein S, Westhoff B, Melzner I, Bassaly B et al: Molecular profiling of chordoma. International journal of oncology 2014, 44(4):1041-1055.
Ge H, Yan Y, Wu D, Huang Y, Tian F: Potential role of LINC00996 in colorectal cancer: a study based on data mining and bioinformatics. OncoTargets and therapy 2018, 11:4845-4855.
He L, Lu Y, Wang P, Zhang J, Yin C, Qu S: Up-regulated expression of type II very low density lipoprotein receptor correlates with cancer metastasis and has a potential link to beta-catenin in different cancers. BMC Cancer 2010, 10:601.
Takahashi S, Kawarabayasi Y, Nakai T, Sakai J, Yamamoto T: Rabbit very low density lipoprotein receptor: a low density lipoprotein receptor-like protein with distinct ligand specificity. Proc Natl Acad Sci U S A 1992, 89(19):9252-9256.
Tanenbaum DG, Hall WA, Colbert LE, Bastien AJ, Brat DJ, Kong J, Kim S, Dwivedi B, Kowalski J, Landry JC et al: TNFRSF10C copy number variation is associated with metastatic colorectal cancer. J Gastrointest Oncol 2016, 7(3):306-314.
Qi Y, Qi L, Qiu M, Yao C, Zhang M, Lin J, Zheng Z, Chen C, Li H, Duan S: Hypermethylation of tumor necrosis factor decoy receptor gene in non-small cell lung cancer. Oncol Lett 2020, 20(1):155-164.
Ma Y, Shepherd J, Zhao D, Bollu LR, Tahaney WM, Hill J, Zhang Y, Mazumdar A, Brown PH: SOX9 Is Essential for Triple-Negative Breast Cancer Cell Survival and Metastasis. Molecular Cancer Research 2020, 18(12):1825-1838.

Download PDF

Journal Publication

published 13 Jul, 2021

Read the published version in Cancer Cell International →

Editorial decision: Major revision
20 Jun, 2021
Review #1 received at journal
15 Jun, 2021
Reviewer #2 agreed at journal
03 Jun, 2021
Review #2 received at journal
03 Jun, 2021
Reviews received at journal
03 Jun, 2021
Reviewers invited by journal
03 Jun, 2021
Reviewer #1 agreed at journal
02 Jun, 2021
Editor assigned by journal
07 Apr, 2021
Submission checks completed at journal
07 Apr, 2021
Editor invited by journal
07 Apr, 2021
First submitted to journal
03 Apr, 2021

You are reading this latest preprint version

New insight into long non-coding RNAs associated with bone metastasis of breast cancer based on an integrated analysis

Status:

Journal Publication

Version 1

Abstract

Background

Methods

Results

Conclusions

Figures

Introduction

Materials and Methods

Results

Discussion

Conclusion

Declarations

Abbreviations

References

Supplementary Files

Status:

Journal Publication

Version 1