ceRNA network development and tumour-infiltrating immune cell analysis of metastatic breast cancer to bone

Genomic Alterations Affecting Competitive Endogenous RNAs (ceRNAs) and Regulatory Networks (ceRNETs) with Clinical Implications in Triple-Negative Breast Cancer (TNBC)

The concept of competitive endogenous RNA regulation has brought on a change in the way we think about transcriptional regulation by miRNA-mRNA interactions. Rather than the relatively simple idea of miRNAs negatively regulating mRNA transcripts, mRNAs and other non-coding RNAs can regulate miRNAs and, therefore, broad networks of gene products through competitive interactions. While this concept is not new, its significant roles in and implications on cancer have just recently come to light. The field is now ripe for the extrapolation of technologies with a substantial clinical impact on cancer. With the majority of the genome consisting of non-coding regions encoding regulatory RNAs, genomic alterations in cancer have considerable effects on these networks that have been previously unappreciated. Triple-negative breast cancer (TNBC) is characterized by high mutational burden, genomic instability and heterogeneity, making this aggressive breast cancer subtype particularly relevant to these changes. In the past few years, much has been learned about the roles of competitive endogenous RNA network regulation in tumorigenesis, disease progression and drug response in triple-negative breast cancer. In this review, we present a comprehensive view of the new knowledge and future perspectives on competitive endogenous RNA networks affected by genomic alterations in triple-negative breast cancer. An overview of the competitive endogenous RNA (ceRNA) hypothesis and its bearing on cellular function and disease is provided, followed by a thorough review of the literature surrounding key competitive endogenous RNAs in triple-negative breast cancer, the genomic alterations affecting them, key disease-relevant molecular and functional pathways regulated by them and the clinical implications and significance of their dysregulation. New knowledge of the roles of these regulatory mechanisms and the current acceleration of research in the field promises to generate insights into the diagnosis, classification and treatment of triple-negative breast cancer through the elucidation of new molecular mechanisms, therapeutic targets and biomarkers.

GABBR2 as a Downstream Effector of the Androgen Receptor Induces Cisplatin Resistance in Bladder Cancer

The precise molecular mechanisms responsible for resistance to cisplatin-based chemotherapy in patients with bladder cancer remain elusive, while we have indicated that androgen receptor (AR) activity in urothelial cancer is associated with its sensitivity. Our DNA microarray analysis in control vs. AR-knockdown bladder cancer sublines suggested that the expression of a GABA B receptor GABBR2 and AR was correlated. The present study aimed to determine the functional role of GABBR2 in modulating cisplatin sensitivity in bladder cancer. AR knockdown and dihydrotestosterone treatment considerably reduced and induced, respectively, GABBR2 expression, and the effect of dihydrotestosterone was at least partially restored by an antiandrogen hydroxyflutamide. A chromatin immunoprecipitation assay further revealed the binding of AR to the promoter region of GABBR2 in bladder cancer cells. Meanwhile, GABBR2 expression was significantly elevated in a cisplatin-resistant bladder cancer subline, compared with control cells. In AR-positive bladder cancer cells, knockdown of GABBR2 or treatment with a selective GABA B receptor antagonist, CGP46381, considerably enhanced the cytotoxic activity of cisplatin. However, no additional effect of CGP46381 on cisplatin-induced growth suppression was seen in GABBR2-knockdown cells. Moreover, in the absence of cisplatin, CGP46381 treatment and GABBR2 knockdown showed no significant changes in cell proliferation or migration. These findings suggest that GABBR2 represents a key downstream effector of AR signaling in inducing resistance to cisplatin treatment. Accordingly, inhibition of GABBR2 has the potential of being a means of chemosensitization, especially in patients with AR/GABBR2-positive bladder cancer.

The current advances of lncRNAs in breast cancer immunobiology research

Breast cancer is the most frequently diagnosed malignancy and the leading cause of cancer-related death in women worldwide. Breast cancer development and progression are mainly associated with tumor-intrinsic alterations in diverse genes and signaling pathways and with tumor-extrinsic dysregulations linked to the tumor immune microenvironment. Significantly, abnormal expression of lncRNAs affects the tumor immune microenvironment characteristics and modulates the behavior of different cancer types, including breast cancer. In this review, we provide the current advances about the role of lncRNAs as tumor-intrinsic and tumor-extrinsic modulators of the antitumoral immune response and the immune microenvironment in breast cancer, as well as lncRNAs which are potential biomarkers of tumor immune microenvironment and clinicopathological characteristics in patients, suggesting that lncRNAs are potential targets for immunotherapy in breast cancer.

Protein Tyrosine Phosphatase Receptor Zeta 1 as a Potential Target in Cancer Therapy and Diagnosis

Protein tyrosine phosphatase receptor zeta 1 (PTPRZ1) is a type V transmembrane tyrosine phosphatase that is highly expressed during embryonic development, while its expression during adulthood is limited. PTPRZ1 is highly detected in the central nervous system, affecting oligodendrocytes' survival and maturation. In gliomas, PTPRZ1 expression is significantly upregulated and is being studied as a potential cancer driver and as a target for therapy. PTPRZ1 expression is also increased in other cancer types, but there are no data on the potential functional significance of this finding. On the other hand, low PTPRZ1 expression seems to be related to a worse prognosis in some cancer types, suggesting that in some cases, it may act as a tumor-suppressor gene. These discrepancies may be due to our limited understanding of PTPRZ1 signaling and tumor microenvironments. In this review, we present evidence on the role of PTPRZ1 in angiogenesis and cancer and discuss the phenomenal differences among the different types of cancer, depending on the regulation of its tyrosine phosphatase activity or ligand binding. Clarifying the involved signaling pathways will lead to its efficient exploitation as a novel therapeutic target or as a biomarker, and the development of proper therapeutic approaches.

ncRNA-mediated ceRNA regulatory network: Transcriptomic insights into breast cancer progression and treatment strategies

With the rapid development of next-generation sequencing technology, several studies have shown that ncRNAs can act as competitive endogenous RNAs (ceRNAs) and are involved in various biological processes, such as proliferation, differentiation, apoptosis, and migration of breast cancer (BC) cells, and plays an important role in BC progression as a molecular target for its diagnosis, treatment, prognosis, and differentiation of subtypes and age groups of BC patients. Based on the description of ceRNA-related biological functions, this study screened and sorted the sequencing analysis and experimental verification conclusions of BC-related ceRNAs and found that the ncRNAs mediated ceRNA networks can promote the development of BC by promoting the expression of genes related to BC proliferation, drug resistance, and apoptosis, inducing the production of epithelial-mesenchymal transition (EMT) to promote metastasis and activating cancer-related signaling pathways.

Long non-coding RNAs in bone metastasis: progresses and perspectives as potential diagnostic and prognostic biomarkers

In a precision medicine perspective, among the biomarkers potentially useful for early diagnosis of cancers, as well as to define their prognosis and eventually to identify novel and more effective therapeutic targets, there are the long non-coding RNAs (lncRNAs). The term lncRNA identifies a class of non-coding RNA molecules involved in the regulation of gene expression that intervene at the transcriptional, post-transcriptional, and epigenetic level. Metastasis is a natural evolution of some malignant tumours, frequently encountered in patients with advanced cancers. Onset and development of metastasis represents a detrimental event that worsen the patient's prognosis by profoundly influencing the quality of life and is responsible for the ominous progression of the disease. Due to the peculiar environment and the biomechanical properties, bone is a preferential site for the secondary growth of breast, prostate and lung cancers. Unfortunately, only palliative and pain therapies are currently available for patients with bone metastases, while no effective and definitive treatments are available. The understanding of pathophysiological basis of bone metastasis formation and progression, as well as the improvement in the clinical management of the patient, are central but challenging topics in basic research and clinical practice. The identification of new molecular species that may have a role as early hallmarks of the metastatic process could open the door to the definition of new, and more effective, therapeutic and diagnostic approaches. Non-coding RNAs species and, particularly, lncRNAs are promising compounds in this setting, and their study may bring to the identification of relevant processes. In this review, we highlight the role of lncRNAs as emerging molecules in mediating the formation and development of bone metastases, as possible biomarkers for cancer diagnosis and prognosis, and as therapeutic targets to counteract cancer spread.

Overexpression of Estrogen Receptor α in Mammary Glands of Aging Mice Is Associated with a Proliferative Risk Signature and Generation of Estrogen Receptor α-Positive Mammary Adenocarcinomas

Age is a risk factor for human estrogen receptor-positive breast cancer, with highest prevalence following menopause. While transcriptome risk profiling is available for human breast cancers, it is not yet developed for prognostication for primary or secondary breast cancer development utilizing at-risk breast tissue. Both estrogen receptor α (ER) and aromatase overexpression have been linked to human breast cancer. Herein, conditional genetically engineered mouse models of estrogen receptor 1 (Esr1) and cytochrome P450 family 19 subfamily A member 1 (CYP19A1) were used to show that induction of Esr1 overexpression just before or with reproductive senescence and maintained through age 30 months resulted in significantly higher prevalence of estrogen receptor-positive adenocarcinomas than CYP19A1 overexpression. All adenocarcinomas tested showed high percentages of ER+ cells. Mammary cancer development was preceded by a persistent proliferative transcriptome risk signature initiated within 1 week of transgene induction that showed parallels to the Prosigna/Prediction Analysis of Microarray 50 human prognostic signature for early-stage human ER+ breast cancer. CYP19A1 mice also developed ER+ mammary cancers, but histology was more divided between adenocarcinoma and adenosquamous, with one ER- adenocarcinoma. Results demonstrate that, like humans, generation of ER+ adenocarcinoma in mice was facilitated by aging mice past the age of reproductive senescence. Esr1 overexpression was associated with a proliferative estrogen pathway-linked signature that preceded appearance of ER+ mammary adenocarcinomas.

Newly Released Advances in the Molecular Mechanisms of Osseous Metastasis and Potential Therapeutic Strategies

The appearance of bone metastases (BM) in individuals with advanced solid cancers (breast, prostate, lung) often worsens their quality of life and prognosis. Although none have been fully effective, several strategies have been used to combat BM. Hence, the need for new data that could be useful for treating bone metastasis. To this end, we reviewed the recent literature on the subject. About patients with prostate cancer, treatments with PIP5K1α inhibitors have been found to inhibit tumor invasion and metastasis, and G protein-coupled receptor class C group 5 member A (GPRC5A) could be a future therapeutic target. Regarding patients with breast cancer, we found the following: Asperolide A could be another curative drug; targeting transforming growth factor-beta (TGFβ) and bone morphogenetic protein (BMP) signaling pathways, along with osteoclast activity, could be a favorable therapeutic approach in the preclusion of osteolytic bone destruction; TRAF6 inhibitors such as 6877002 appear promising; aiming the BMP4-SMAD7 signaling axis is an innovative therapeutic approach; there is favorable proof for the plausible therapeutic utilization of bone aiming immunostimulatory MOF (BT-isMOF) nanoparticles, and inhibition of IL4R and macrophages could have therapeutic benefits. For lung cancer, the function of LIGHT in osteolytic osseous illness instigated by metastatic non-small cell lung cancer should be highlighted.

Correlation between Genes of the ceRNA Network and Tumor-Infiltrating Immune Cells and Their Biomarker Screening in Kidney Renal Clear Cell Carcinoma

This study aimed to using bioinformatics tools, qPCR, and the immunohistochemical analysis to find out factors related to the early diagnosis and prognosis of kidney renal clear cell carcinoma (KIRC). The expression profiles of lncRNA, miRNA, and mRNA of KIRC were downloaded from The Cancer Genome Atlas database. A ceRNA regulatory network was constructed based on the interaction between these three differentially expressed genes. The CIBERSORT deconvolution algorithm was used to analyze the differential distribution of 22 types of immune cells. The Kaplan–Meier survival and Cox analyses were used to screen genes of the ceRNA network and also immune cell subtypes related to the clinical and prognostic prediction of KIRC. Co-expression regulatory relationships were found among LINC01426, LINC00894, CCNA2, L1 cell adhesion molecule (L1CAM), and T follicular helper cells, which served as potential biomarkers. The results of quantitative reverse transcriptase-polymerase chain reaction showed that LINC01426 was upregulated while L1CAM was downregulated in KIRC, but no difference was found in the expression levels of LINC00894 and CCNA2 in cancer and adjacent samples. The immunohistochemical analysis showed that T follicular helper cells were more concentrated in core tissues and metastases of KIRC. In a word, co-expression relationships were found among LINC01426, L1CAM, and T follicular helper cells, and they may serve as biomarkers for early diagnosis and prognostic evaluation of KIRC.

Development and Validation of Prognostic Nomogram for Elderly Breast Cancer: A Large-Cohort Retrospective Study

Purpose

Our research aims to study the bone metastatic patterns and prognostic outcomes in elderly breast cancer (BC) and to develop elder-specific nomograms.

Methods

We downloaded the data of BC patients between 2010 and 2016 from the Surveillance, Epidemiology, and End Results database. The differences in clinical features and prognosis between young (age < 65) and elderly (age ≥ 65) BC patients were compared. The univariate and multivariate Cox analyses were used to determine the overall survival (OS)- and cancer-specific survival (CSS)-related variables and establish two nomograms of BC patients with bone metastasis (BCBM). The receiver operating characteristic (ROC) curve with area under the curve (AUC), calibration curve, decision curve analysis (DCA), and Kaplan–Meier survival curve were selected to evaluate nomograms.

Results

A total of 230,177 BC patients were enrolled in our research, including 142,025 young and 88,152 elderly patients. The prognosis of elderly BCBM patients was significantly worse than young patients. Age, race, breast subtype, tumor size, tumor grade, brain metastasis, liver metastasis, surgery, and chemotherapy were independent prognostic variables for elderly BCBM patients, including OS and CSS. The AUC values at 12, 18, and 24 months were 0.750, 0.751, and 0.739 for OS nomogram and 0.759, 0.762, and 0.752 for CSS nomogram in the training cohort, which were higher than the AUC values of all single independent prognostic variables. The survival curve showed a distinct prognosis between low-, median- and high-risk groups (p < 0.001). Finally, calibration curves and DCA indicated that both nomograms have favorable performance.

Conclusion

Elderly and young patients presented with different bone metastatic frequencies, clinical features, and prognostic outcomes. Two elder-specific nomograms incorporating nine clinical variables were established and validated to be a valuable predictor for elderly BCBM patients.

Identification of Novel Biomarkers for Predicting Prognosis and Immunotherapy Response in Head and Neck Squamous Cell Carcinoma Based on ceRNA Network and Immune Infiltration Analysis

Objectives

Patients with head and neck squamous cell carcinoma (HNSCC) have poor prognosis and show poor responses to immune checkpoint (IC) inhibitor (ICI) therapy. Competing endogenous RNA (ceRNA) networks, tumor-infiltrating immune cells (TIICs), and ICIs may influence tumor prognosis and response rates to ICI therapy. This study is aimed at identifying prognostic and IC-related biomarkers and key TIIC signatures to improve prognosis and ICI therapy response in HNSCC patients.

Methods and Results

Ninety-five long noncoding RNAs (lncRNAs), microRNAs (miRNAs), and 1746 mRNAs were identified using three independent methods. We constructed a ceRNA network and estimated the proportions of 22 immune cell subtypes. Ten ceRNAs were related to prognosis according to Kaplan–Meier analysis. Two risk signatures based, respectively, on nine ceRNAs (ANLN, CFL2, ITGA5, KDELC1, KIF23, NFIA, PTX3, RELT, and TMC7) and three immune cell types (naïve B cells, neutrophils, and regulatory T cells) via univariate Cox regression, least absolute shrinkage and selection operator, and multivariate Cox regression analyses could accurately and independently predict the prognosis of HNSCC patients. Key mRNAs in the ceRNA network were significantly correlated with naïve B cells and regulatory T cells and with stage, grade, and immune and molecular subtype. Eight IC genes exhibited higher expression in tumor tissues and were correlated with eight key mRNAs in the ceRNA network in HNSCC patients with different HPV statuses according to coexpression and TIMER 2.0 analyses. Most drugs were effective in association with expression of these key signatures (ANLN, CFL2, ITGA5, KIF23, NFIA, PTX3, RELT, and TMC7) based on GSCALite analysis. The prognostic value of key biomarkers and associations between key ceRNAs and IC genes were validated using online databases. Eight key ceRNAs were confirmed to predict response to ICI in other cancers based on TIDE analysis.

Conclusions

We constructed two risk signatures to accurately predict prognosis in HNSCC. Key IC-related signatures may be associated with response to ICI therapy. Combinations of ICIs with inhibitors of eight key mRNAs may improve survival outcomes of HNSCC patients.

Identification of a Competing Endogenous RNA Network Related to Immune Signature in Lung Adenocarcinoma

Background

The establishment of immunotherapy has led to a new era in oncotherapy. But the signature of immune-related genes (IRGs) in LUAD remains to be elucidated. Here we use integrated analysis to identify IRGs roles in immune signature and detect their relationship with competing endogenous RNA (ceRNA) networks in LUAD progression.

Methods

By analyzing the RNA-seq data from different platforms, we recognized the differentially expressed genes (DEGs) of each platform and screened out the top 20 hub IRGs related to immune responses. Then, we applied the CIBERSORT algorithm to explore the landscape of tumor-infiltrating immune cells (TILs) in LUAD and their connection with hub genes. Next, we predicted and validated the upstream miRNAs and lncRNAs according to their expression and prognostic roles. Finally, we constructed and validated an immune-related ceRNA network by co-expression analysis.

Results

A total of 71 IRGs were identified among 248 DEGs, which play key roles in immune responses. CIBERSORT analysis showed that six hub genes were closely related to TILs, such as SPP1 and naive B cells (R = −0.17), TEK and resting mast cells (R = 0.37). Stepwise prediction and validation from mRNA to lncRNA, including 6 hub genes, 5 miRNAs, and 9 lncRNAs, were applied to construct a ceRNA network. Ultimately, we confirmed the TMPO-AS1/miR-126-5p/SPP1 and CARD8-AS1/miR-21-5p/TEK as immune-related ceRNA networks in LUAD progression.

Conclusion

We elucidated two immune-related ceRNA networks in LUAD progression, which can be considered as immunotherapy targets for this disease.

Construction and analysis of a ceRNA network and patterns of immune infiltration in bladder cancer

Background

Bladder cancer (BC) is the ninth most common malignant tumor, accounting for an estimate of 549,000 new BC cases and 200,000 BC-related deaths worldwide in 2018. The prognosis of BC has not substantially improved despite significant advances in the diagnosis and treatment of the disease.

Methods

The RNA sequencing (RNA-seq) data and clinical information of BC patients were downloaded from The Cancer Genome Atlas (TCGA) database. The Cell-type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORT) algorithm was used to assess immune infiltration. The survival analyses were performed using the selected components of a ceRNA network and selected immune cell types by least absolute shrinkage and selection operator (LASSO) Cox regression to calculate the risk score. The accuracy of prognosis prediction was determined by receiver operating characteristic (ROC) curves, survival curves, and nomograms. Finally, the correlation analysis was performed to investigate the relationships between the signature components of the ceRNA network and the immune cell signature.

Results

Two completed survival analyses included selected components of the ceRNA network (ELN, SREBF1, DSC2, TTLL7, DIP2C, SATB1, hsa-miR-20a-5p, and hsa-miR-29c-3p) and selected immune cell types (M0 macrophages, M2 macrophages, resting mast cells, and neutrophils). ROC curves, survival curves (all P values <0.05), nomograms, and calibration curves indicated that the accuracy of the two survival analyses was acceptable. Moreover, the correlations between TTLL7 and resting mast cells (R=0.24, P<0.001), DSC2 and resting mast cells (R=−0.23, P<0.001), ELN and resting mast cells (R=0.44, P<0.001), and hsa-miR-29c-3p and M0 macrophages (R=−0.29, P<0.001) were significant, indicating that interactions of these factors may play significant roles in the prognosis of BC.

Conclusions

TTLL7, DSC2, ELN, hsa-miR-29c-3p, resting mast cells, and M0 macrophages may play an important role in the development of BC. However, additional studies are needed to confirm this hypothesis.

DLX6-AS1 is a potential biomarker and therapeutic target in cancer initiation and progression

Long noncoding RNAs (lncRNAs) are involved in multiple functions such as the regulation of cellular homeostasis. They play prominent roles in the pathogenesis of human cancer, and contribute to every hallmark of cancer. The novel cancer-related lncRNA DLX6 antisense RNA 1 (DLX6-AS1) plays an essential regulatory role in enhancing and initiating carcinogenesis and tumor progression. This progression is due to the aberrant regulation of downstream factors in vitro as well as in vivo. DLX6-AS1 is significantly dysregulated in various cancers. DLX6-AS1 functions in tumor initiation and progression are regulated at the epigenetic, transcription, and posttranscriptional regulation levels. DLX6-AS1 functions as an oncogene, binding to miRNA targeting sites competing endogenous RNAs and causing the upregulation of downstream tumor-related genes and carcinogenesis. The regulation and detailed molecular mechanisms of DLX6-AS1 and its potential role in malignancies are comprehensively described in this paper. DLX6-AS1 has the potential to become a novel biomarker and therapeutic target for cancer.

Construction of an RNA-Binding Protein-Related Prognostic Model for Pancreatic Adenocarcinoma Based on TCGA and GTEx Databases

Background: Recently, RNA-binding proteins (RBPs) were reported to interact with target mRNA to regulate gene posttranscriptional expression, and RBP-mediated RNA modification can regulate the expression and function of proto-oncogenes and tumor suppressor genes. We systematically analyzed the expression of RBPs in pancreatic adenocarcinoma (PAAD) and constructed an RBP-associated prognostic risk model. Methods: Gene expression data of normal pancreatic samples as well as PAAD samples were downloaded from TCGA-PAAD and GTEx databases. Wilcoxon test and univariate Cox analysis were, respectively, applied to screen differential expression RBPs (DE-RBPs) and prognostic-associated RBPs (pRBPs). Functional enrichment was analyzed by GO, KEGG, and GSEA. Protein-protein interaction (PPI) network was constructed by STRING online database. Modeling RBPs were selected by multivariate Cox analysis. Kaplan-Meier survival and Cox analysis were applied to evaluate the effects of risk score on the overall survival of PAAD patients. ROC curves and validation cohort were applied to verify the accuracy of the model. Nomogram was applied for predicting 1-, 3-, and 5-year overall survival (OS) of PAAD patients. At last, modeling RBPs were further analyzed to explore their differential expression, prognostic value, as well as enrichment pathways in PAAD. Results: RBPs (453) were differentially expressed in normal and tumor samples, besides, 28 of which were prognostic associated. DE-RBPs (453) are functionally associated with ribosome, ribonuclease, spliceosome, etc. Eight RBPs (PABPC1, PRPF6, OAS1, RBM5, LSM12, IPO7, FXR1, and RBM6) were identified to construct a prognostic risk model. Higher risk score not only predicted poor prognosis but also was an independent poor prognostic indicator, which was verified by ROC curves and validation cohort. Eight modeling RBPs were confirmed to be significantly differentially expressed between normal and tumor samples from RNA and protein level. Besides, all of eight RBPs were related with overall survival of PAAD patients. Conclusions: We successfully constructed an RBP-associated prognostic risk model in PAAD, which has a potential clinical application prospect.