- 1Department of Intensive Care Unit, Changxing People’s Hospital of Zhejiang, Huzhou, China
- 2Department of Anesthesiology, The First Affiliated Hospital, College of Medicine, Zhejiang University, Hangzhou, China
- 3Department of Oncology, The First Affiliated Hospital of Jiaxing University, Jiaxing, China
The role and mechanism of inflammation in breast cancer is unclear. This study aims to probe the relationship between inflammation and long non-coding RNAs (lncRNAs) and to stablish an inflammation-related competing endogenous RNA (ceRNA) network in breast cancer. Inflammation-related lncRNAs and target genes were screened based on the data from four single-cell RNA sequencing (scRNA-seq) studies and miRNAs were bioinformatically predicted according to ceRNA hypothesis. A series of in silico analyses were performed to construct an inflammation-related ceRNA network in breast cancer. Consequently, a total of seven inflammation-related lncRNAs were selected, after which LRRC75A-AS1 was identified as the most potential lncRNA in view of its expression and prognostic predictive value in breast cancer. Finally, an inflammation-related ceRNA network in breast cancer at the single cell level was established based on lncRNA LRRC75A-AS1, miR-3127-5p, miR-2114-3p, RPL36 and RPL27A mRNAs. Collectively, the lncRNA LRRC75A-AS1 and the LRRC75A-AS1-based on ceRNA network may exert crucial roles in modulating inflammation response during the initiation and progression of breast cancer.
Introduction
Breast cancer ranks the most common type of human cancer and is also the leading cause of cancer-related deaths in women all over the world (Sung et al., 2021). The prognosis of patients with breast cancer, especially triple-negative breast cancer and HER2-positive breast cancer, is still unsatisfactory, in spite of huge advancements in the aspect of its diagnosis and therapy have been achieved during the past decades (de la Mare et al., 2014). More work and efforts need to be put in the extensive and in-depth research regarding the molecular mechanism of breast carcinogenesis and progression.
When Rudolph Virchow noted leucocytes in tumor tissues in 1863 and proposed the theory of “lymphoreticular infiltrate”, the origin of cancer at sites of inflammation came into researchers’ sight (Balkwill and Mantovani, 2001). After that inflammation has been recognized to contribute to malignant transformation of several types of human malignancies, including hepatocellular carcinoma, colorectal cancer and gastric cancer (Multhoff et al., 2011). Inflammation is also a key component of breast carcinogenesis (Danforth, 2021). However, its roles and molecular mechanism in the initiation and progression of breast cancer is still unclear and need to be further elucidated.
Single-cell RNA sequencing (scRNA-seq) technology, a powerful tool to investigate cell heterogeneity, bringing new strategies to deeply understand the characteristics and behaviors of human cancer (including breast cancer) at the single-cell level (Ren et al., 2021). In this study, we first screened out several inflammation-related lncRNAs in breast cancer using four breast cancer-related scRNA-seq studies, then predicted and analyzed corresponding downstream microRNAs (miRNAs) based on the competing endogenous RNA (ceRNA) hypothesis (Salmena et al., 2011), and finally also identified inflammation-related target gens in breast cancer at the single-cell level. At the end, an inflammation-related lncRNA-miRNA-mRNA triple RNA regulatory network in breast cancer has been established.
Results
Identification of LRRC75A-AS1 as a Potential Inflammation-Related Long Non-Coding RNA in Breast Cancer
ScRNA-seq provides an unprecedented opportunity to probe the functional heterogeneity of cancer cells. It has been widely acknowledged that inflammation is closely linked to occurrence and progression of various types of human cancers, including breast cancer. To elucidate the underlying molecular mechanism of the inflammation in breast cancer, four breast cancer-associated scRNA-seq studies were collected (Balavenkatraman et al., 2011; Braune et al., 2016; Jordan et al., 2016; Chung et al., 2017). Consequently, a total of 656 lncRNAs that were significantly associated with inflammation in breast cancer were screened out (Data not shown). Among these lncRNAs, only seven members, consisting of LRRC75A-AS1, RGPD4-AS1, AC010255.1, MALAT1, AC159540.2, AC091891.2 and FAM239A, appear in more than one scRNA-seq studies. As listed in Table 1, all the seven lncRNAs were negatively correlated with inflammation in breast cancer. Next, we successively determined the expression levels of the seven lncRNAs in breast cancer by GEPIA database and starBase database. As presented in Figure 1, only LRRC75A-AS1 expression was markedly decreased in breast cancer when compared with normal controls in both two databases. These findings suggested that LRRC75A-AS1 might be the most potential inflammation-related lncRNA in breast cancer.
FIGURE 1. Identification and validation of potential inflammation-related lncRNAs in breast cancer. The expression of LRRC75A-AS1 (A), RGPD4-AS1 (B), AC010255.1 (C), MALAT1 (D), AC159540.2 (E), AC091891.2 (F) and FAM239A (G) in breast cancer compared with normal controls determined by GEPIA and starBase databases. *p < 0.05.
Survival Analysis for LRRC75A-AS1 in Various Immunocyte-Enriched/Decreased Breast Cancer
As mentioned above, LRRC75A-AS1 was identified as a potential inflammation-related lncRNA in breast cancer. As is known to all, inflammation and immunity are two ends of the balance in human body. Thus, we performed survival analysis for inflammation-related LRRC75A-AS1 in breast cancer enriched or decreased with various immunocytes. Nine different immunocytes, involving basophil, B cell, CD4+ T cell, CD8+ T cell, eosinophil, macrophage, mesenchymal stem cell, natural killer T cell and regulatory T cell, were included for this survival analysis. As shown in Figure 2, high expression of LRRC75A-AS1 indicated favorable overall survival in breast cancer enriched with basophil or natural killer T cell. Moreover, upregulation of LRRC75A-AS1 was positively correlated with good overall survival in breast cancer decreased with B cell, CD8+ T cell, eosinophil or mesenchymal stem cell (Figure 3). No significant prognostic values of LRRC75A-AS1 in breast cancer enriched or decreased with other immunocytes were observed. Taken together, some immunocytes, including basophil, natural killer T cell, B cell, CD8+ T cell, eosinophil or mesenchymal stem cell, might be involved in interaction with LRRC75A-AS1-mediated inflammation in breast cancer.
FIGURE 2. Survival analysis for LRRC75A-AS1 in breast cancer enriched with various immunocyte. (A) The prognostic value of LRRC75A-AS1 in basophil-enriched breast cancer. (B) The prognostic value of LRRC75A-AS1 in B cell-enriched breast cancer. (C) The prognostic value of LRRC75A-AS1 in CD4+ cell-enriched breast cancer. (D) The prognostic value of LRRC75A-AS1 in CD8+ cell-enriched breast cancer. (E) The prognostic value of LRRC75A-AS1 in eosinophil-enriched breast cancer. (F) The prognostic value of LRRC75A-AS1 in macrophage-enriched breast cancer. (G) The prognostic value of LRRC75A-AS1 in mesenchymal stem cell-enriched breast cancer. (H) The prognostic value of LRRC75A-AS1 in natural killer T cell-enriched breast cancer. (I) The prognostic value of LRRC75A-AS1 in regulatory T cell-enriched breast cancer.
FIGURE 3. Survival analysis for LRRC75A-AS1 in breast cancer decreased with various immunocyte. (A) The prognostic value of LRRC75A-AS1 in basophil-decreased breast cancer. (B) The prognostic value of LRRC75A-AS1 in B cell-decreased breast cancer. (C) The prognostic value of LRRC75A-AS1 in CD4+ cell-decreased breast cancer. (D) The prognostic value of LRRC75A-AS1 in CD8+ cell-decreased breast cancer. (E) The prognostic value of LRRC75A-AS1 in eosinophil-decreased breast cancer. (F) The prognostic value of LRRC75A-AS1 in macrophage-decreased breast cancer. (G) The prognostic value of LRRC75A-AS1 in mesenchymal stem cell-decreased breast cancer. (H) The prognostic value of LRRC75A-AS1 in natural killer T cell-decreased breast cancer. (I) The prognostic value of LRRC75A-AS1 in regulatory T cell-decreased breast cancer.
Prediction and Analysis of Downstream Potential miRNAs of LRRC75A-AS1 in Breast Cancer
To further explore the downstream action mechanism of inflammation-related LRRC75A-AS1 in breast cancer, miRNet was employed to predict the miRNAs that potentially bind to LRRC75A-AS1. A total of 70 miRNAs were identified. For better visualization, a LRRC75A-AS1-miRNA network was established using Cytoscape software (Figure 4). According to competing endogenous RNA (ceRNA) hypothesis, these should be negative expression correlation between lncRNA and miRNA. Therefore, expression correlation of LRRC75A-AS1 with its predicted miRNAs in breast cancer was calculated using TCGA breast cancer data. As listed in Table 2, 15 of 70 miRNAs were significantly inversely associated with LRRC75A-AS1 in breast cancer. Then, we further analyzed the expression levels of the 15 negatively correlated miRNAs in breast cancer using starBase database. As presented in Figure 5A, 6 miRNAs were significantly downregulated in breast cancer tissues compared with normal breast tissues, including miR-149-5p, miR-3127-5p, miR-330-5p, miR-2114-3p, miR-1277-3p and miR-760. The prognostic values of the overexpressed six miRNAs in breast cancer were also assessed by Kaplan-Meier plotter database. Only breast cancer patients with high expression of miR-3127-5p or miR-2114-3p had poor prognosis. By combination of correlation analysis, expression analysis and survival analysis, miR-3127-5p and miR-2114-3p might be two most potential downstream binding miRNAs of LRRC75A-AS1 in breast cancer.
TABLE 2. The expression correlation of LRRC75A-AS1 with predicted miRNAs in breast cancer determined by starBase database.
FIGURE 5. Expression analysis and survival analysis for the miRNAs negatively correlated with LRRC75A-AS1 in breast cancer. (A) The expression levels of 15 LRRC75A-AS1-related miRNAs in breast cancer detected by starBase database. Red: overexpressed in breast cancer; Green: downregulated in breast cancer; gray: no significant differences between breast cancer and normal controls. (B–G) The prognostic value of miR-149-5p (B), miR-3127-5p (C), miR-330-5p (D), miR-2114-3p (E) and miR-1277-3p (F) and miR-760 (G) in breast cancer.
Identification of Downstream Inflammation-Related Targets of LRRC75A-AS1-miRNA Pathways in Breast Cancer
Next, we intended to explore the downstream molecular mechanism of inflammation-related LRRC75A-AS1-miR-3127-5p/miR-2114-3p pathways in breast cancer. First of all, these genes that were significantly associated with inflammation in breast cancer were obtained by analyzing the four scRNA-seq studies as mentioned above. 2,452 inflammation-related genes were identified (Data not shown). Moreover, as is known to all, miRNAs usually exert their roles by suppressing target genes. Thus, the potential target genes of miR-3127-5p or miR-2114-3p were predicted by miRNet database. 64 and 75 genes were forecasted to be downstream targets of miR-3127-5p and miR-2114-3p, respectively. For better visualization, miR-3127-5p-gene and miR-2114-3p-gene networks were established as shown in Figure 6A and Figure 6B. By intersection of target genes and inflammation-related genes, 11 and 9 potential inflammation-related target genes of miR-3127-5p and miR-2114-3p were respectively identified as presented in Figures 6C,D.
FIGURE 6. Identification of downstream potential targets of miR-3127-5p and miR-2114-3p in breast cancer. (A) Establishment of a miR-3127-5p-gene network. (B) Establishment of a miR-2114-3p-gene network. (C,D) Intersection of predicted target genes of miR-3127-5p (C) or miR-2114-3p (D) and inflammation-related genes in breast cancer.
Construction of an Inflammation-Related LRRC75A-AS1-miRNA-mRNA Network in Breast Cancer
Based on the negative regulation mechanism of miRNA, expression correlation analysis was conducted. As listed in Table 3, 3 of 11 target genes (DMTN, MCL1 and RPL36) were significantly inversely correlated with miR-3127-5p expression in breast cancer. And two target genes (ULK2 and RPL36) were significantly positively correlated with LRRC75A-AS1 expression in breast cancer. Thus, RPL36 might be the most potential inflammation-related target gene of LRRC75A-AS1/miR-3127-5p pathway in breast cancer. In the same light, three target genes (FABP7, RPL27A and TMCC2) were markedly negatively associated with miR-2114-3p expression in breast cancer and two target genes (RPL27A and TMCC2) were markedly positively linked to LRRC75A-AS1 expression in breast cancer as shown in Table 4. Therefore, RPL27A and TMCC2 might be two potential inflammation-related target genes of LRRC75A-AS1/miR-2114-3p pathway in breast cancer. Finally, the prognostic values of RPL36, RPL27A and TMCC2 in breast cancer were evaluated using Kaplan-Meier plotter database. As presented in Figure 7, high expression of RPL36 and RPL27A indicated favorable prognosis in breast cancer. No statistical role of TMCC2 in predicting prognosis of breast cancer was observed. Taken together, an inflammation-related LRRC75A-AS1-mediated miRNA-mRNA pathways in breast cancer were identified as vividly depicted in Figure 8.
TABLE 3. The expression correlation of miR-3127-5p with target genes or target genes with LRRC75A-AS1 in breast cancer.
TABLE 4. The expression correlation of miR-2114-3p with target genes or target genes with LRRC75A-AS1 in breast cancer.
FIGURE 7. Survival analysis for RPL36, RPL27A and TMCC2 in breast cancer. (A) The prognostic value of RPL36 in breast cancer. (B) The plot of significance versus cutoff values between lower and upper quartiles of RPL36 expression. (C) The prognostic value of RPL27A in breast cancer. (D) The plot of significance versus cutoff values between lower and upper quartiles of RPL27A expression. (E) The prognostic value of TMCC2 in breast cancer. (F) The plot of significance versus cutoff values between lower and upper quartiles of TMCC2 expression.
FIGURE 8. The model of LRRC75A-AS1-mediated inflammation-related miRNA-mRNA pathways in breast cancer.
Discussion
Despite inflammation is present in breast cancer and may affect its outcome, the molecular mechanism of inflammation in breast carcinogenesis is still unknown. In this study, we aimed to identify a lncRNA-miRNA-mRNA regulatory network based on ceRNA hypothesis at single-cell level.
CancerSEA was employed to identify inflammation-related lncRNAs in breast cancer at single-cell resolution. By integrated analysis of lncRNAs’ expression, LRRC75A-AS1 was identified as the most potential inflammation-associated lncRNA in breast cancer. Several studies have confirmed its key roles in human cancer. For example, Pang et al. suggested that LRRC75A-AS1 suppressed multiple myeloma by targeting miR-199b-5p/PDCD4 axis (Pang et al., 2020); Chen et al. demonstrated that LRRC75A-AS1 repressed cell proliferation and migration of colorectal carcinoma (Chen et al., 2019); Li et al. confirmed that LRRC75A-AS1 facilitated cell proliferation and invasion of triple negative breast cancer by acting as a ceRNA to regulate BAALC (Li et al., 2020). Tumor infiltrating lymphocytes (TILs) are closely linked to prognosis of cancer patients (Yu S et al., 2020; Lou et al., 2021). Moreover, survival analysis also revealed that LRRC75A-AS1 possessed significant prognostic value in breast cancer enriched with basophil, natural killer T cell, or decreased with B cell, CD8+ T cell, eosinophil or mesenchymal stem cell.
ceRNA mechanism can partially account for lncRNA’s action mechanism (Lou et al., 2020). Using miRNet database, a series of miRNAs were predicted. By combination of correlation analysis, expression analysis and survival analysis, miR-3127-5p and miR-2114-3p were selected as the two most potential downstream binding miRNAs of LRRC75A-AS1 in breast cancer. Both of the two miRNAs have been reported to be involved in initiation and progression of human malignancies, including osteosarcoma (Wang et al., 2021), papillary thyroid cancer (Zhou et al., 2021), colorectal cancer (Ma et al., 2021) and ovarian cancer (Yu Z et al., 2020).
Subsequently, the downstream target genes of miR-3127-5p or miR-2114-3p were predicted by online target gene prediction tool. Moreover, the inflammation-related genes in breast cancer were also identified using four breast cancer scRNA-seq studies as mentioned above. By intersection of target gene set and inflammation-related gene sets, 11 and 9 candidate genes were screened out for miR-3127-5p and miR-2114-3p, respectively. After performing correlation analysis and survival analysis, RPL36 was identified as the most potential target for miR-3127-5p and RPL27A was for miR-2114-3p.
RPL36 was reported as a tumor suppressor to restrain KRAS-induced pancreatic cancer (Provost et al., 2014), and was a promising prognostic marker in hepatocellular carcinoma (Song et al., 2011). RPL27A was also found to be involved in development and metastasis of triple-negative breast cancer in both mouse and human (Zhao et al., 2021). Our analytic results together with these reports indicated that RPL36 and RPL27A might be downstream key targets of LRRC75A-AS1/miR-3127-5p/miR-2114-3p pathways in breast cancer.
At the end, a novel inflammation-related lncRNA-miRNA-mRNA regulatory network in breast cancer was successively established based on ceRNA hypothesis and scRNA-seq data. There were several limitations in this study. For example, the miRNAs in this network were not identified using scRNA-seq data; the results from this study were only based on in silico analysis. However, this is the first study to comprehensively explore inflammation-related lncRNA-miRNA-mRNA network in breast cancer. All these findings need to be further validated using basic experiments and clinical trials in the future.
Materials and Methods
Screening of Inflammation-Related Long Non-Coding RNAs and Genes
CancerSEA (http://biocc.hrbmu.edu.cn/CancerSEA/home.jsp) is the first dedicated database that aims to comprehensively decode distinct functional states of cancer cells at single-cell resolution (Yuan et al., 2019). In this study, CancerSEA database was employed to identify inflammation-related lncRNAs and genes in breast cancer. Only those lncRNAs and genes that were significantly associated with inflammation in breast cancer were included. p-value < 0.05 was considered as statistically significant.
Expression Analysis for Long Non-Coding RNA
Two databases, including GEPIA (http://gepia.cancer-pku.cn/) and starBase (http://starbase.sysu.edu.cn/), were used to analyze the expression of lncRNAs in breast cancer (Li et al., 2014; Tang et al., 2017). Only lncRNAs that were significantly upregulated or downregulated in both two databases were selected for subsequent analysis. p-value < 0.05 was considered as statistically significant.
Survival Analysis for Long Non-Coding RNA
The prognostic values of LRRC75A-AS1 in breast cancer were determined using Kaplan-Meier plotter (http://kmplot.com/analysis/), which is an online database capable of accessing the effects of genes or miRNAs on survival in more than 20 cancer types including breast cancer (Lou et al., 2021). Logrank p-value < 0.05 was considered as statistically significant.
MicroRNA Prediction
miRNet (http://www.mirnet.ca), a database for miRNA functional analysis and systems biology, was used to predict the binding miRNAs of LRRC75A-AS1 (Chang et al., 2020). Then, a LRRC75A-AS1-miRNA network was established using Cytoscape software.
Expression and Survival Analysis for MicroRNA
The expression levels of miRNAs of LRRC75A-AS1 in breast cancer were determined by usage of starBase (http://starbase.sysu.edu.cn/). p-value < 0.05 was considered as statistically significant. Kaplan-Meier plotter (http://kmplot.com/analysis/) was employed to assess the prognostic values of miRNAs of LRRC75A-AS1. The survival plots were downloaded from the website. Logrank p-value < 0.05 was considered as statistically significant.
Target Gene Prediction
miRNet (http://www.mirnet.ca) was introduced to predict the downstream target genes of miR-3127-5p and miR-2114-3p. Then, a miR-3127-5p-target gene network or miR-2114-3p-target gene network was constructed using Cytoscape software.
Intersection Analysis
VENNY 2.1 (https://bioinfogp.cnb.csic.es/tools/venny/index.html) was employed to perform intersection analysis for target genes of miRNAs and inflammation-related genes in breast cancer. Only the genes that commonly appeared in both two gene sets were included for subsequent analysis.
Correlation Analysis
starBase (http://starbase.sysu.edu.cn/) was used to conduct expression correlation analysis for LRRC75A-AS1-miRNA, miRNA-target gene or LRRC75A-AS1-target gene pairs in breast cancer. Only those RNA-RNA pairs with p-value less than 0.05 were considered as statistically significant.
Statistical Analysis
The statistical analyses were directly performed using online databases or tools as mentioned above. p-value < 0.05 or logrank p-value < 0.05 was considered as statistically significant.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author Contributions
JC: designed this work; JL: performed in silico analyses and wrote this manuscript. SX: performed some in silico analyses. JC: revised the manuscript. The authors have read and approved the final version of the manuscript.
Funding
This project was supported from the Scientific Research Fund of Zhejiang Provincial Education Department (Y202043573) and 2019 Jiaxing Key Discipiline of Medicine-Oncology (Supporting subject, No. 2019-ZC-11).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2022.839876/full#supplementary-material
References
Balavenkatraman, K. K., Aceto, N., Britschgi, A., Mueller, U., Bence, K. K., Neel, B. G., et al. (2011). Epithelial Protein-Tyrosine Phosphatase 1B Contributes to the Induction of Mammary Tumors by HER2/Neu but Is Not Essential for Tumor Maintenance. Mol. Cancer Res. 9, 1377–1384. doi:10.1158/1541-7786.mcr-11-0198
Balkwill, F., and Mantovani, A. (2001). Inflammation and Cancer: Back to Virchow? Lancet 357, 539–545. doi:10.1016/s0140-6736(00)04046-0
Braune, E.-B., Tsoi, Y. L., Phoon, Y. P., Landor, S., Silva Cascales, H., Ramsköld, D., et al. (2016). Loss of CSL Unlocks a Hypoxic Response and Enhanced Tumor Growth Potential in Breast Cancer Cells. Stem Cel Rep. 6, 643–651. doi:10.1016/j.stemcr.2016.03.004
Chang, L., Zhou, G., Soufan, O., and Xia, J. (2020). miRNet 2.0: Network-Based Visual Analytics for miRNA Functional Analysis and Systems Biology. Nucleic Acids Res. 48, W244–w251. doi:10.1093/nar/gkaa467
Chen, J., Lan, J., Ye, Z., Duan, S., Hu, Y., Zou, Y., et al. (2019). Long Noncoding RNA LRRC75A-AS1 Inhibits Cell Proliferation and Migration in Colorectal Carcinoma. Exp. Biol. Med. (Maywood) 244, 1137–1143. doi:10.1177/1535370219874339
Chung, W., Eum, H. H., Lee, H.-O., Lee, K.-M., Lee, H.-B., Kim, K.-T., et al. (2017). Single-cell RNA-Seq Enables Comprehensive Tumour and Immune Cell Profiling in Primary Breast Cancer. Nat. Commun. 8, 15081. doi:10.1038/ncomms15081
Danforth, D. N. (2021). The Role of Chronic Inflammation in the Development of Breast Cancer. Cancers (Basel) 13, 3918. doi:10.3390/cancers13153918
de la Mare, J.-A., Contu, L., Hunter, M., Moyo, B., Sterrenberg, J., Dhanani, K., et al. (2014). Breast Cancer: Current Developments in Molecular Approaches to Diagnosis and Treatment. Recent Pat Anticancer Drug Discov. 9, 153–175. doi:10.2174/15748928113086660046
Jordan, N. V., Bardia, A., Wittner, B. S., Benes, C., Ligorio, M., Zheng, Y., et al. (2016). HER2 Expression Identifies Dynamic Functional States within Circulating Breast Cancer Cells. Nature 537, 102–106. doi:10.1038/nature19328
Li, J.-H., Liu, S., Zhou, H., Qu, L.-H., and Yang, J.-H. (2014). starBase v2.0: Decoding miRNA-ceRNA, miRNA-ncRNA and Protein-RNA Interaction Networks from Large-Scale CLIP-Seq Data. Nucl. Acids Res. 42, D92–D97. doi:10.1093/nar/gkt1248
Li, S., Wu, D., Jia, H., and Zhang, Z. (2020). Long Non-coding RNA LRRC75A-AS1 Facilitates Triple Negative Breast Cancer Cell Proliferation and Invasion via Functioning as a ceRNA to Modulate BAALC. Cell Death Dis. 11, 643. doi:10.1038/s41419-020-02821-2
Lou, W., Ding, B., and Fu, P. (2020). Pseudogene-Derived lncRNAs and Their miRNA Sponging Mechanism in Human Cancer. Front. Cel Dev. Biol. 8, 85. doi:10.3389/fcell.2020.00085
Lou, W., Wang, W., Chen, J., Wang, S., and Huang, Y. (2021). ncRNAs-mediated High Expression of SEMA3F Correlates with Poor Prognosis and Tumor Immune Infiltration of Hepatocellular Carcinoma. Mol. Ther. Nucleic Acids 24, 845–855. doi:10.1016/j.omtn.2021.03.014
Ma, Y., Li, Y., Tang, Y., Tang, N., Wang, D., and Li, X. (2021). LINC00174 Facilitates Proliferation and Migration of Colorectal Cancer Cells via MiR-3127-5p/E2F7 Axis. J. Microbiol. Biotechnol. 31, 1098–1108. doi:10.4014/jmb.2103.03001
Multhoff, G., Molls, M., and Radons, J. (2011). Chronic Inflammation in Cancer Development. Front. Immunol. 2, 98. doi:10.3389/fimmu.2011.00098
Pang, Q., Wang, Y., Bi, D., and Lu, H. (2020). LRRC75A-AS1 Targets miR-199b-5p/PDCD4 axis to Repress Multiple Myeloma. Cancer Biol. Ther. 21, 1051–1059. doi:10.1080/15384047.2020.1831373
Provost, E., Bailey, J. M., Aldrugh, S., Liu, S., Iacobuzio-Donahue, C., and Leach, S. D. (2014). The Tumor Suppressor Rpl36 Restrains KRASG12V-Induced Pancreatic Cancer. Zebrafish 11, 551–559. doi:10.1089/zeb.2014.1024
Ren, L., Li, J., Wang, C., Lou, Z., Gao, S., Zhao, L., et al. (2021). Single Cell RNA Sequencing for Breast Cancer: Present and Future. Cell Death Discov. 7, 104. doi:10.1038/s41420-021-00485-1
Salmena, L., Poliseno, L., Tay, Y., Kats, L., and Pandolfi, P. P. (2011). A ceRNA Hypothesis: the Rosetta Stone of a Hidden RNA Language? Cell 146, 353–358. doi:10.1016/j.cell.2011.07.014
Song, M. J., Jung, C. K., Park, C.-H., Hur, W., Choi, J. E., Bae, S. H., et al. (2011). RPL36 as a Prognostic Marker in Hepatocellular Carcinoma. Pathol. Int. 61, 638–644. doi:10.1111/j.1440-1827.2011.02716.x
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA A. Cancer J. Clin. 71, 209–249. doi:10.3322/caac.21660
Tang, Z., Li, C., Kang, B., Gao, G., Li, C., and Zhang, Z. (2017). GEPIA: a Web Server for Cancer and normal Gene Expression Profiling and Interactive Analyses. Nucleic Acids Res. 45, W98–w102. doi:10.1093/nar/gkx247
Wang, X., Yu, X., Long, X., and Pu, Q. (2021). MIR205 Host Gene (MIR205HG) Drives Osteosarcoma Metastasis via Regulating the microRNA 2114-3p (miR-2114-3p)/twist Family bHLH Transcription Factor 2 (TWIST2) axis. Bioengineered 12, 1576–1586. doi:10.1080/21655979.2021.1920326
Yu S, S., Hu, C., Cai, L., Du, X., Lin, F., Yu, Q., et al. (2020). Seven-Gene Signature Based on Glycolysis Is Closely Related to the Prognosis and Tumor Immune Infiltration of Patients with Gastric Cancer. Front. Oncol. 10, 1778. doi:10.3389/fonc.2020.01778
Yu Z, Z., Wang, Y., Wang, B., and Zhai, J. (2020). Metformin Affects Paclitaxel Sensitivity of Ovarian Cancer Cells through Autophagy Mediated by Long Noncoding RNASNHG7/miR-3127-5p Axis. Cancer Biother. Radiopharm. doi:10.1089/cbr.2019.3390
Yuan, H., Yan, M., Zhang, G., Liu, W., Deng, C., Liao, G., et al. (2019). CancerSEA: a Cancer Single-Cell State Atlas. Nucleic Acids Res. 47, D900–d908. doi:10.1093/nar/gky939
Zhao, W., Li, X., Nian, W., Wang, J., Wang, X., Sun, L., et al. (2021). Ribosome Proteins Represented by RPL27A Mark the Development and Metastasis of Triple-Negative Breast Cancer in Mouse and Human. Front. Cel Dev. Biol. 9, 716730. doi:10.3389/fcell.2021.716730
Keywords: inflammation, breast cancer, LRRC75A-AS1, microRNA (miRNA), prognosis, single cell RNA sequencing (scRNA-seq)
Citation: Liu J, Xiao S and Chen J (2022) Development of an Inflammation-Related lncRNA-miRNA-mRNA Network Based on Competing Endogenous RNA in Breast Cancer at Single-Cell Resolution. Front. Cell Dev. Biol. 10:839876. doi: 10.3389/fcell.2022.839876
Received: 20 December 2021; Accepted: 06 January 2022;
Published: 25 January 2022.
Edited by:
Weiyang Lou, Zhejiang University, ChinaReviewed by:
Junjie Zhang, Fuyang City People’s Hospital, ChinaHuilin Zheng, Zhejiang University of Science and Technology, China
Copyright © 2022 Liu, Xiao and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jing Chen, b1818220@zju.edu.cn