Novel MicroRNA-Based Risk Score Identified by Integrated Analyses to Predict Metastasis and Poor Prognosis in Breast Cancer

Decoding the interaction between miR-19a and CBX7 focusing on the implications for tumor suppression in cancer therapy

Cancer is a complex and multifaceted disease characterized by uncontrolled cell growth, genetic alterations, and disruption of normal cellular processes, leading to the formation of malignant tumors with potentially devastating consequences for patients. Molecular research is important in the diagnosis and treatment, one of the molecular mechanisms involved in various cancers is the fluctuation of gene expression. Non-coding RNAs, especially microRNAs, are involved in different stages of cancer. MicroRNAs are small RNA molecules that are naturally produced within cells and bind to the 3'-UTR of target mRNA, repressing gene expression by regulating translation. Overexpression of miR-19a has been reported in human malignancies. Upregulation of miR-19a as a member of the miR-17-92 cluster is key to tumor formation, cell proliferation, survival, invasion, metastasis, and drug resistance. Furthermore. bioinformatics and in vitro data reveal that the miR-19a-3p isoform binds to the 3'UTR of CBX7 and was identified as the miR-19a-3p target gene. CBX7 is known as a tumor suppressor. This review initially describes the regulation of mir-19a in multiple cancers. Accordingly, the roles of miR-19 in affecting its target gene expression CBX7 in carcinoma also be discussed.

miRNA Expression Profiling in Human Breast Cancer Diagnostics and Therapy

Breast cancer is one of the most commonly diagnosed cancer types worldwide. Regarding molecular characteristics and classification, it is a heterogeneous disease, which makes it more challenging to diagnose. As is commonly known, early detection plays a pivotal role in decreasing mortality and providing a better prognosis for all patients. Different treatment strategies can be adjusted based on tumor progression and molecular characteristics, including personalized therapies. However, dealing with resistance to drugs and recurrence is a challenge. The therapeutic options are limited and can still lead to poor clinical outcomes. This review aims to shed light on the current perspective on the role of miRNAs in breast cancer diagnostics, characteristics, and prognosis. We discuss the potential role of selected non-coding RNAs most commonly associated with breast cancer. These include miR-21, miR-106a, miR-155, miR-141, let-7c, miR-335, miR-126, miR-199a, miR-101, and miR-9, which are perceived as potential biomarkers in breast cancer prognosis, diagnostics, and treatment response monitoring. As miRNAs differ in expression levels in different types of cancer, they may provide novel cancer therapy strategies. However, some limitations regarding dynamic alterations, tissue-specific profiles, and detection methods must also be raised.

Advances in the study of biomarkers related to bone metastasis in breast cancer

Breast cancer is by far the most common malignancy in females. And bone is the most common site of distant metastasis in breast cancer, accounting for about 65 to 75% of all metastatic breast cancer patients.1,2Bone metastasis is an important factor affecting the prognosis of breast cancer. When patients have early-stage breast cancer without metastasis, their 5-year survival rate is as high as 90%, and once metastasis occurs, their 5-year survival rate will drop to 10%.3 Bone radionuclide imaging (ECT), X-ray, CT scan, MRI and other imaging tests to diagnose breast cancer bone metastasis are commonly used in clinical, It is currently believed that breast cancer bone metastasis is a multistep process: first, breast cancer cells need to acquire invasive and metastatic properties; breast cancer cells enter the blood circulation and migrate from blood breast cancer cells enter the blood circulation and migrate from blood vessels to bone tissue in a targeted manner; breast cancer cells adhere and remain in bone tissue and colonise it; and finally, it leads to bone destruction.4 Several key molecules are involved in breast cancer bone metastasis, and serum biomarkers are generally able to detect pathological changes earlier Several key molecules are involved in breast cancer bone metastasis, and serum biomarkers are generally able to detect pathological changes earlier than imaging.5 This review describes the progress of serum biomarkers for breast cancer bone metastasis.

Chromatin-Accessible miRNA Regulons Driving Thyroid Tumorigenesis and Progression

BACKGROUND

Although papillary thyroid cancer can remain indolent, associated lymph node metastases and recurrence rates are approximately 50% and 20%, respectively. Omics-based medicine has led to the discovery of predictive biomarkers that can be used to predict tumor progression and clinical outcomes. We aimed to develop a noninvasive omics-driven blood test to allow accurate risk stratification and help tailor individual patient treatment plans.

STUDY DESIGN

RNA sequencing (seq) and microRNA analysis of The Cancer Genome Atlas and Gene Expression Omnibus datasets were employed to identify an epigenetic prognostic panel. Integrated bulk assay for transposase-accessible chromatin-seq and RNA-seq experiments confirmed the results. Sixty-two paired tumor and adjacent control thyroid tissues and 67 blood samples (62 papillary thyroid cancer and 5 controls) were analyzed for validation using sequencing and real-time polymerase chain reaction and correlated to clinical outcomes. A liposome-exosome fusion clustered regularly interspaced short palindromic repeats (CRISPR)-fluorescent detection system miRNA assay was developed. A predictive risk nomogram was generated and tested for performance.

RESULTS

Our miRNA panel (miR-146b-5p and miR-221-3p) from tissue and blood was associated with aggressive features and was located within accessible chromatin regions. The miRNA risk score and prognostic nomogram showed higher accuracy in predicting lymph node metastases (miR-146b: area under the curve [AUC] 0.816, sensitivity 76.9%; miR-221: AUC 0.740, sensitivity 79.5%) and recurrence (miR-146b: AUC 0.921, sensitivity 75.0%; miR-221: AUC 0.756, sensitivity 70.0%; p < 0.001) than staging and American Thyroid Association risk stratification. CRISPR-based miRNA assays showed upregulation in the blood of cancer cohorts.

CONCLUSIONS

CRISPR-based detection of miR-146b and miR-221 in the blood of thyroid cancer patients is a reliable and noninvasive tool for real-time assessment and prognostication that has great potential to provide a direct impact on the care of these patients.

The Role of Gender-Related Immune Genes in Childhood Acute Myeloid Leukemia

The study of immune genes and immune cells is highly focused in recent years. To find immunological genes with prognostic value, the current study examines childhood acute myeloid leukemia according to gender. The TARGET database was used to gather the "mRNA expression profile data" and relevant clinical data of children with AML. To normalize processing and find differentially expressed genes (DEG) between male and female subgroups, the limma software package is utilized. We identified prognostic-related genes and built models using LASSO, multivariate Cox, and univariate Cox analysis. The prognostic significance of prognostic genes was then examined through the processing of survival analysis and risk score (RS) calculation. We investigated the connections between immune cells and prognostic genes as well as the connections between prognostic genes and medications. Finally, five immune genes from the TARGET database have been identified. These immune genes are considerably correlated to the prognosis of male patients.

Low adipocyte hepatocellular carcinoma is associated with aggressive cancer biology and with worse survival

Hepatocellular carcinoma (HCC) is the second leading cause of cancer-related death worldwide, and non-alcoholic fatty liver disease is strongly associated with its development. To explore the role of adipocytes in HCC, we investigated intratumoral adipocytes, also known as cancer-associated adipocytes (CAA). Based on our prior breast cancer findings, we hypothesized that low intratumoral adipocytes would be associated with aggressive cancer biology, worse tumor microenvironment (TME), and clinical outcomes. The Cancer Genome Atlas (TCGA) was used and validated by the Gene Expression Omnibus (GEO) cohort. xCell algorithm was used to quantify intratumoral adipocytes and top 90% were defined as adipocyte high (AH) and bottom 10% as adipocyte low (AL). We found that AL-HCC was significantly associated with worse disease-free survival (DFS), disease-specific survival (DSS), and overall survival (OS). AL-HCC were higher-grade, had high MKI67 expression, enriched cell proliferation-related gene sets, and had increased altered fraction, aneuploidy, and homologous recombination defects. Also, anti-cancer immune cells, CD8, Th1, and M1 cells, as well as pro-cancer Th2 cells were increased in AL-HCC. Micro-RNAs miR-122 (associated with cholesterol metabolism) and miR-885 (associated with liver pathologies) were significantly increased in the AL TME. In conclusion, we found that AL-HCC has worse patient outcomes and is biologically more aggressive with enhanced cell proliferation. Our findings take initial steps to clarify the role of adipocytes in HCC.

Left sided breast cancer is associated with aggressive biology and worse outcomes than right sided breast cancer

Breast cancer is more common on the left side than the right side. We aim to evaluate differences in clinicopathological and genomic characteristics based on laterality. We analyzed survival outcomes and clinical characteristics of 881,320 patients recorded by the Surveillance, Epidemiology, and End Results (SEER) program. The Cancer Genome Atlas (TCGA) was used to explore genomic and clinical features from 1,062 patients. Gene expression data was used to quantitate cytolytic activity and hallmark gene-sets were used for gene set enrichment analysis. An institutional retrospective review was conducted on 155 patients treated with neoadjuvant chemotherapy (NACT). Patient characteristics were summarized by pathological complete response (pCR). Left sided tumors were found to be more prevalent than right sided tumors. No major clinicopathological differences were noted by laterality. Left sided breast cancer demonstrated poorer outcomes versus right sided tumors (HR 1.05, 95% CI 1.01–1.08; p = 0.011). Cell proliferation gene sets, including E2F Targets, G2M Checkpoint, Mitotic spindle, and MYC Targets, were enriched on the left side compared to the right. Left sided tumors had lower pCR rates versus right sided tumors (15.4% versus 29.9%, p = 0.036). Our findings suggest that left sided breast cancer is associated with aggressive biology and worse outcomes compared to right sided breast cancer.

Predictive and Prognostic Value of Non-Coding RNA in Breast Cancer

For decades since the central dogma, cancer biology research has been focusing on the involvement of genes encoding proteins. It has been not until more recent times that a new molecular class has been discovered, named non-coding RNA (ncRNA), which has been shown to play crucial roles in shaping the activity of cells. An extraordinary number of studies has shown that ncRNAs represent an extensive and prevalent group of RNAs, including both oncogenic or tumor suppressive molecules. Henceforth, various clinical trials involving ncRNAs as extraordinary biomarkers or therapies have started to emerge. In this review, we will focus on the prognostic and diagnostic role of ncRNAs for breast cancer.

PARP-1 modulates the expression of miR-223 through histone acetylation to involve in the hydroquinone-induced carcinogenesis of TK6 cells

The upstream regulators of microRNAs were rarely reported. Hydroquinone (HQ) is the main metabolite of benzene, one of the important environmental factors contributing to leukemia and lymphoma. In HQ-induced malignant transformed TK6 (TK6-HT) cells, the expression of PARP-1 and miR-223 were upregulated. When in PARP-1 silencing TK6-HT cells, miR-223 was downregulated and the apoptotic cell number correspondingly increased. In TK6 cells treated with HQ for different terms, the expression of miR-223 and PARP-1 were dynamically observed and found to be decreased and increased, respectively. Trichostatin A could increase the expression of miR-223, then the expression of HDAC1-2 and nuclear factor kappa B were found to be increased, but that of mH2A was decreased. PARP-1 silencing inhibited the protein expression of H3Ac, mH2A, and H3K27ac. By co-immunoprecipitation experiment, PARP-1 and HDAC2 were found to form a regulatory complex. In conclusion, we demonstrated that the upregulation of PARP-1 mediated activation of acetylation to promote the transcription of miR-223 possibly via coregulating with HDAC2, an epigenetic regulation mechanism involved in cell malignant transformation resulting from long-term exposure to HQ, in which course, H3K27ac might be a specific marker for the activation of histone H3, which also gives hints for benzene exposure research.

MicroRNAs: Emerging Regulators of Metastatic Bone Disease in Breast Cancer

Bone metastasis is a frequent complication in patients with advanced breast cancer. Once in the bone, cancer cells disrupt the tightly regulated cellular balance within the bone microenvironment, leading to excessive bone destruction and further tumor growth. Physiological and pathological interactions in the bone marrow are mediated by cell-cell contacts and secreted molecules that include soluble proteins as well as RNA molecules. MicroRNAs (miRNAs) are short non-coding RNAs that post-transcriptionally interfere with their target messenger RNA (mRNA) and subsequently reduce protein abundance. Since their discovery, miRNAs have been identified as critical regulators of physiological and pathological processes, including breast cancer and associated metastatic bone disease. Depending on their targets, miRNAs can exhibit pro-tumorigenic or anti-tumorigenic functions and serve as diagnostic and prognostic biomarkers. These properties have encouraged pre-clinical and clinical development programs to investigate miRNAs as biomarkers and therapeutic targets in various diseases, including metastatic cancers. In this review, we discuss the role of miRNAs in metastatic bone disease with a focus on breast cancer and the bone microenvironment and elaborate on their potential use for diagnostic and therapeutic purposes in metastatic bone disease and beyond.

Five crucial prognostic-related autophagy genes stratified female breast cancer patients aged 40-60 years

BACKGROUND

Autophagy is closely related to the progression of breast cancer. The aim at this study is to establish a prognostic-related model comprised of hub autophagy genes (AGs) to assess patient prognosis. Simultaneously, the model can guide clinicians to make up individualized strategies and stratify patients aged 40-60 years based on risk level.

METHODS

The hub AGs were identified with univariate COX regression and LASSO regression. The functions and alterations of these selected AGs were analyzed as well. Moreover, the multivariate COX regression and correlation analysis between hub AGs and clinicopathological parameters were done.

RESULTS

Totally, 33 prognostic-related AGs were obtained from the univariate COX regression (P < 0.05). SERPINA1, HSPA8, HSPB8, MAP1LC3A, and DIRAS3 were identified to constitute the prognostic model by the LASSO regression. The survival curve of patients in the high-risk and low-risk groups was statistically significant (P < 0.05). The 3-year and 5-year ROC displayed that their AUC value reached 0.762 and 0.825, respectively. Stage and risk scores were independent risk factors relevant to prognosis. RB1CC1, RPS6KB1, and BIRC6 were identified as the most predominant mutant genes. It was found that AGs were mainly involved in regulating the endopeptidases synthesis and played important roles in the ErbB signal pathway. SERPIN1, risk score was closely related to the stage (P < 0.05); HSPA8, risk score were closely related to T stag (P < 0.05); HSPB8 was closely related to N stag (P < 0.05).

CONCLUSIONS

Our prognostic model had the relatively robust predictive ability on prognosis for patients aged 40-60 years. If the stage was added into the prognostic model, the predictive ability would be more powerful.

MicroRNA-Based Risk Score for Predicting Tumor Progression Following Radioactive Iodine Ablation in Well-Differentiated Thyroid Cancer Patients: A Propensity-Score Matched Analysis

Simple Summary

The three-tiered American Thyroid Association (ATA) risk stratification helps clinicians tailor decisions regarding follow-up modalities and the need for postoperative radioactive iodine (RAI) ablation and radiotherapy. However, a significant number of well-differentiated thyroid cancers (DTC) progress after treatment. Current follow-up modalities have also been proposed to detect disease relapse and recurrence but have failed to be sufficiently sensitive or specific to detect, monitor, or determine progression. Therefore, we assessed the predictive accuracy of the microRNA-based risk score in DTC with and without postoperative RAI. We confirm the prognostic role of triad biomarkers (miR-2f04, miR-221, and miR-222) with higher sensitivity and specificity for predicting disease progression than the ATA risk score. Compared to indolent tumors, a higher risk score was found in progressive samples and was associated with shorter survival. Consequently, our prognostic microRNA signature and nomogram provide a clinically practical and reliable ancillary measure to determine the prognosis of DTC patients.

Abstract

To identify molecular markers that can accurately predict aggressive tumor behavior at the time of surgery, a propensity-matching score analysis of archived specimens yielded two similar datasets of DTC patients (with and without RAI). Bioinformatically selected microRNAs were quantified by qRT-PCR. The risk score was generated using Cox regression and assessed using ROC, C-statistic, and Brier-score. A predictive Bayesian nomogram was established. External validation was performed, and causal network analysis was generated. Within the eight-year follow-up period, progression was reported in 51.5% of cases; of these, 48.6% had the T1a/b stage. Analysis showed upregulation of miR-221-3p and miR-222-3p and downregulation of miR-204-5p in 68 paired cancer tissues (p < 0.001). These three miRNAs were not differentially expressed in RAI and non-RAI groups. The ATA risk score showed poor discriminative ability (AUC = 0.518, p = 0.80). In contrast, the microRNA-based risk score showed high accuracy in predicting tumor progression in the whole cohorts (median = 1.87 vs. 0.39, AUC = 0.944) and RAI group (2.23 vs. 0.37, AUC = 0.979) at the cutoff >0.86 (92.6% accuracy, 88.6% sensitivity, 97% specificity) in the whole cohorts (C-statistics = 0.943/Brier = 0.083) and RAI subgroup (C-statistic = 0.978/Brier = 0.049). The high-score group had a three-fold increased progression risk (hazard ratio = 2.71, 95%CI = 1.86–3.96, p < 0.001) and shorter survival times (17.3 vs. 70.79 months, p < 0.001). Our prognostic microRNA signature and nomogram showed excellent predictive accuracy for progression-free survival in DTC.

Low expression of miR-29a is associated with aggressive biology and worse survival in gastric cancer

Advanced gastric cancer (GC) is one of the most lethal cancer types, thus a better understanding of its biology in patients is urgently needed. MicroRNA (miR)-29a is a known tumor suppressive miR that is related to metastasis, but its clinical relevance in GC remains ambiguous. Here, using a large GC patient cohort we hypothesized that low expression of miR-29a in GC is associated with aggressive cancer biology and worse survival. We demonstrated that low miR-29a GC enriched cell proliferation, apoptosis, metastasis, and angiogenesis related gene sets, as well as the higher expression of related genes. Low miR-29a GC was associated with less anti-cancer immune cell infiltration as well as immune related scoring. Low miR-29a GC demonstrated a worse overall survival (OS) as well as disease specific survival (DSS) compared with high expressing miR-29a GC. Notably, low miR-29a expression was the only factor, other than residual tumor status, to be an independent prognostic biomarker of worse OS and DSS. In conclusion, low miR-29a GC was associated with aggressive cancer biology and worse OS as well as DSS. Additionally, low expression of miR-29a was an independent prognostic biomarker of OS and DSS in gastric cancer patients.

Low expression of miR-195 is associated with cell proliferation, glycolysis and poor survival in estrogen receptor (ER)-positive but not in triple negative breast cancer

MiR-195 is a tumor suppressive microRNA in breast cancer. Its clinical relevance remains debatable as it has only been studied via in vitro experiments or small cohort studies. We analyzed a total of 2,038 patients in the TCGA and METABRIC cohorts to assess whether low miR-195 expressing tumors are associated with aggressive cancer characteristics and poor prognostic outcomes. The median cutoff of miR-195 expression was used to split the groups into miR-195 high and low groups. Low miR-19 expressing tumors demonstrated high cell proliferating features by enriching the gene sets associated with cell proliferation, MKI67 expression and pathological grade. One-third of the top target miR-195 genes were related to cell proliferation. Low miR-195 expressing tumors were associated with both pro-cancerous and anti-cancerous immune cells. Low miR-195 expressing tumors were associated with enhanced glycolysis and poor survival in ER-positive tumors, but not other subtypes of breast cancer. In conclusion, low expression of miR-195 in ER-positive breast cancer was associated with enhanced cancer cell proliferation, glycolysis, and worse overall survival.

Octogenarians' Breast Cancer Is Associated with an Unfavorable Tumor Immune Microenvironment and Worse Disease-Free Survival

Elderly patients are known to have a worse prognosis for breast cancer. This is commonly blamed on their medical comorbidities and access to care. However, in addition to these social issues, we hypothesized that the extreme elderly (octogenarians-patients over 80 years old) have biologically worse cancer with unfavorable tumor immune microenvironment. The Cancer Genomic Atlas (TCGA) and the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) breast cancer cohorts were analyzed. The control (aged 40-65) and octogenarians numbered 668 and 53 in TCGA and 979 and 118 in METABRIC, respectively. Octogenarians had significantly worse breast cancer-specific survival in both cohorts (p < 0.01). Octogenarians had a higher ER-positive subtype rate than controls in both cohorts. Regarding PAM50 classification, luminal-A and -B subtypes were significantly higher in octogenarians, whereas basal and claudin-low subtypes were significantly lower (p < 0.05) in octogenarians. There was no difference in tumor mutation load, intratumor heterogeneity, or cytolytic activity by age. However, the octogenarian cohort was significantly associated with high infiltration of pro-cancer immune cells, M2 macrophage, and regulatory T cells in both cohorts (p < 0.05). Our results demonstrate that octogenarians' breast cancer is associated with worse survival and with an unfavorable tumor immune microenvironment.

Urine as a Source of Liquid Biopsy for Cancer

Simple Summary

Tissue biopsy is essential for diagnosis and characterization of a tumor. Recently circulating tumor cells and other tumor-derived nucleic acid can be detected from blood, which is called liquid biopsy. Now this concept has been expanded to many other body fluids including urine. Urine is the least invasive method to obtain a liquid biopsy and can be done anywhere, which allows longitudinal repeated sampling. Here, we review the latest update on urine liquid biopsy in urological and non-urological cancers.

Abstract

Tissue biopsy is the gold standard for diagnosis and morphological and immunohistochemical analyses to characterize cancer. However, tissue biopsy usually requires an invasive procedure, and it can be challenging depending on the condition of the patient and the location of the tumor. Even liquid biopsy analysis of body fluids such as blood, saliva, gastric juice, sweat, tears and cerebrospinal fluid may require invasive procedures to obtain samples. Liquid biopsy can be applied to circulating tumor cells (CTCs) or nucleic acids (NAs) in blood. Recently, urine has gained popularity due to its less invasive sampling, ability to easily repeat samples, and ability to follow tumor evolution in real-time, making it a powerful tool for diagnosis and treatment monitoring in cancer patients. With the development and advancements in extraction methods of urinary substances, urinary NAs have been found to be closely related to carcinogenesis, metastasis, and therapeutic response, not only in urological cancers but also in non-urological cancers. This review mainly highlights the components of urine liquid biopsy and their utility and limitations in oncology, especially in non-urological cancers.

Regulation of bone metastasis and metastasis suppressors by non-coding RNAs in breast cancer

Breast cancer (BC) is a critical health care issue that substantially affects women worldwide. Though surgery and chemotherapy can effectively control tumor growth, metastasis remains a primary concern. Metastatic BC cells predominantly colonize in bone, owing to their rigid osseous nutrient-rich nature. There are recently increasing studies investigating the context-dependent roles of non-coding RNAs (ncRNAs) in metastasis regulation. ncRNAs, including microRNAs, long non-coding RNAs, circular RNAs, and small interference RNAs, control the BC metastasis via altered mechanisms. Additionally, these ncRNAs have been reported in regulating a unique class of genes known as Metastatic suppressors. Metastasis suppressors like BRMS1, NM23, LIFR, and KAI1, etc., have been extensively studied for their role in inducing apoptosis, inhibiting metastasis, and maintaining homeostasis. In this review, we have emphasized the direct regulation of ncRNAs for effectively controlling the distant spread of BC. Furthermore, we have highlighted the ncRNA-mediated modulation of the metastatic suppressors, thereby delineating their indirect influence over metastasis.

Role of miRNA-19a in Cancer Diagnosis and Poor Prognosis

Cancer is a multifactorial disease that affects millions of people every year and is one of the most common causes of death in the world. The high mortality rate is very often linked to late diagnosis; in fact, nowadays there are a lack of efficient and specific markers for the early diagnosis and prognosis of cancer. In recent years, the discovery of new diagnostic markers, including microRNAs (miRNAs), has been an important turning point for cancer research. miRNAs are small, endogenous, non-coding RNAs that regulate gene expression. Compelling evidence has showed that many miRNAs are aberrantly expressed in human carcinomas and can act with either tumor-promoting or tumor-suppressing functions. miR-19a is one of the most investigated miRNAs, whose dysregulated expression is involved in different types of tumors and has been potentially associated with the prognosis of cancer patients. The aim of this review is to investigate the role of miR-19a in cancer, highlighting its involvement in cell proliferation, cell growth, cell death, tissue invasion and migration, as well as in angiogenesis. On these bases, miR-19a could prove to be truly useful as a potential diagnostic, prognostic, and therapeutic marker.

Application of "omics" sciences to the prediction of bone metastases from breast cancer: State of the art

Breast cancer (BC) is the most frequent malignancy and the first cause of cancer-related death in women. The majority of patients with advanced BC develop skeletal metastases which may ultimately lead to serious complications, termed skeletal-related events, that often dramatically impact on quality of life and survival. Therefore, the identification of biomarkers able to stratify BC patient risk to develop bone metastases (BM) is fundamental to define personalized diagnostic and therapeutic strategies, possibly at the earliest stages of the disease. In this regard, the advent of "omics" sciences boosted the investigation of several putative biomarkers of BC osteotropism, including deregulated genes, proteins and microRNAs. The present review revisits the current knowledge on BM development in BC and the most recent studies exploring potential BM-predicting biomarkers, based on the application of omics sciences to the study of primary breast malignancies.

Identification of a Novel Glycolysis-Related Gene Signature for Predicting Breast Cancer Survival

To identify a glycolysis-related gene signature for the evaluation of prognosis in patients with breast cancer, we analyzed the data of a training set from TCGA database and four validation cohorts from the GEO and ICGC databases which included 1,632 patients with breast cancer. We conducted GSEA, univariate Cox regression, LASSO, and multiple Cox regression analysis. Finally, an 11-gene signature related to glycolysis for predicting survival in patients with breast cancer was developed. And Kaplan–Meier analysis and ROC analyses suggested that the signature showed a good prognostic ability for BC in the TCGA, ICGC, and GEO datasets. The analyses of univariate Cox regression and multivariate Cox regression revealed that it’s an important prognostic factor independent of multiple clinical features. Moreover, a prognostic nomogram, combining the gene signature and clinical characteristics of patients, was constructed. These findings provide insights into the identification of breast cancer patients with a poor prognosis.

High G2M Pathway Score Pancreatic Cancer is Associated with Worse Survival, Particularly after Margin-Positive (R1 or R2) Resection

Pancreatic cancer is highly mortal due to uncontrolled cell proliferation. The G2M checkpoint pathway is an essential part of the cell cycle. We hypothesized that a high G2M pathway score is associated with cell proliferation and worse survival in pancreatic cancer patients. Gene set variation analysis using the Hallmark G2M checkpoint gene set was used as a score to analyze a total of 390 human pancreatic cancer patients from 3 cohorts (TCGA, GSE62452, GSE57495). High G2M score tumors enriched other cell proliferation genes sets as well as MKI67 expression, pathological grade, and proliferation score. Independent of other prognostic factors, G2M score was predictive of disease-specific survival in pancreatic cancer. High G2M tumor was associated with high mutation rate of KRAS and TP53 and significantly enriched these pathway gene sets, as well as high infiltration of Th2 cells. High G2M score consistently associated with worse overall survival in 3 cohorts, particularly in R1/2 resection, but not in R0. High G2M tumor in R1/2 highly enriched metabolic and cellular components' gene sets compared to R0. To our knowledge, this is the first study to use gene set variation analysis as a score to examine the clinical relevancy of the G2M pathway in pancreatic cancer.

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

PURPOSE

Advanced breast cancer commonly metastasises to bone; however, the molecular mechanisms underlying the affinity for breast cancer cells to bone remains unclear. Thus, we developed nomograms based on a competing endogenous RNA (ceRNA) network and analysed tumour-infiltrating immune cells to elucidate the molecular pathways that may predict prognosis in patients with breast cancer.

METHODS

We obtained the RNA expression profile of 1091 primary breast cancer samples included in The Cancer Genome Atlas database, 58 of which were from patients with bone metastasis. We analysed the differential RNA expression patterns between breast cancer with and without bone metastasis and developed a ceRNA network. Cibersort was employed to differentiate between immune cell types based on tumour transcripts. Nomograms were then established based on the ceRNA network and immune cell analysis. The value of prognostic factors was evaluated by Kaplan-Meier survival analysis and a Cox proportional risk model.

RESULTS

We found significant differences in long non-coding RNAs (lncRNAs), 18 microRNAs (miRNAs), and 20 messenger RNAs (mRNAs) between breast cancer with and without bone metastasis, which were used to construct a ceRNA network. We found that the protein-coding genes GJB3, CAMMV, PTPRZ1, and FBN3 were significantly differentially expressed by Kaplan-Meier analysis. We also observed significant differences in the abundance of plasma cell and follicular helper T cell populations between the two groups. In addition, the proportion of mast cells, gamma delta T cells, and plasma cells differed depending on disease location and stage. Our analysis showed that a high proportion of follicular helper T cells and a low proportion of eosinophils promoted survival and that DLX6-AS1, Wnt6, and GABBR2 expression may be associated with bone metastasis in breast cancer.

CONCLUSIONS

We developed a bioinformatic tool for exploring the molecular mechanisms of bone metastasis in patients with breast cancer and identified factors that may predict the occurrence of bone metastasis.

The E2F Pathway Score as a Predictive Biomarker of Response to Neoadjuvant Therapy in ER+/HER2- Breast Cancer

E2F transcription factors play critical roles in the cell cycle. Therefore, their activity is expected to reflect tumor aggressiveness and responsiveness to therapy. We scored 3905 tumors of nine breast cancer cohorts for this activity based on their gene expression for the Hallmark E2F targets gene set. As expected, tumors with a high score had an increased expression of cell proliferation-related genes. A high score was significantly associated with shorter patient survival, greater MKI67 expression, histological grade, stage, and genomic aberrations. Furthermore, metastatic tumors had higher E2F scores than the primary tumors from which they arose. Although tumors with a high score had greater infiltration by both pro- and anti-cancerous immune cells, they had an increased expression of immune checkpoint genes. Estrogen receptor (ER)-positive/human epidermal growth factor receptor 2 (HER2)-negative cancer with a high E2F score achieved a significantly higher pathological complete response (pCR) rate to neoadjuvant chemotherapy. The E2F score was significantly associated with the expression of cyclin-dependent kinase (CDK)-related genes and strongly correlated with sensitivity to CDK inhibition in cell lines. In conclusion, the E2F score is a marker of breast cancer aggressiveness and predicts the responsiveness of ER-positive/HER2-negative patients to neoadjuvant chemotherapy and possibly to CDK and immune checkpoint inhibitors.

Contribution of Immune Cells to Glucocorticoid Receptor Expression in Breast Cancer

Breast cancer (BC) patients experience increased stress with elevated cortisol levels, increasing risk of cancer recurrence. Cortisol binds to a cytoplasmic receptor, glucocorticoid receptor (GR) encoded by GR gene (NR3C1). We hypothesized that not only cancer cells, but even immune cells in the tumor microenvironment (TME) may contribute to GR expression in bulk tumor and influence prognosis. To test this, mRNA expression data was accessed from METABRIC and TCGA. "High" and "low" expression was based on highest and lowest quartiles of NR3C1 gene expression, respectively. Single-cell sequencing data were obtained from GSE75688 and GSE114725 cohorts. Computer algorithms CIBERSORT, Gene Set Enrichment Analysis and TIMER were used. GR-high BC has better median disease-free and disease-specific survival. Single cell sequencing data showed higher GR expression on immune cells compared to cancer and stromal cells. Positive correlation between GR-high BC and CD8+ T-cells was noted. In GR-high tumors, higher cytolytic activity (CYT) with decreased T-regulatory and T-follicular helper cells was observed. High GR expression was associated with lower proliferation index Ki67, enriched in IL-2_STAT5, apoptosis, KRAS, TGF-β signaling, and epithelial-to-mesenchymal transition. Immune cells significantly contribute to GR expression of bulk BC. GR-high BC has a favorable TME with higher CYT with favorable outcomes.

Current and Emerging Biomarkers Predicting Bone Metastasis Development

Bone is one of the preferential sites of distant metastases from malignant tumors, with the highest prevalence observed in breast and prostate cancers. Patients with bone metastases (BMs) may experience skeletal-related events, such as severe bone pain, pathological fractures, spinal cord compression, and hypercalcemia, with negative effects on the quality of life. In the last decades, a deeper understanding of the molecular mechanisms underlying the BM onset has been gained, leading to the development of bone-targeting agents. So far, most of the research has been focused on the pathophysiology and treatment of BM, with only relatively few studies investigating potential predictors of risk for BM development. The ability to select such "high-risk" patients could allow early identification of those most likely to benefit from interventions to prevent or delay BM. This review summarizes several evidences for the potential use of specific biomarkers able to predict early the BM development.

High Expression of microRNA-143 is Associated with Favorable Tumor Immune Microenvironment and Better Survival in Estrogen Receptor Positive Breast Cancer

microRNA-143 (miR-143) is a well-known tumor suppressive microRNA that exhibits anti-tumoral function by targeting KRAS signaling pathways in various malignancies. We hypothesized that miR-143 suppresses breast cancer progression by targeting KRAS and its effector molecules. We further hypothesized that high expression of miR-143 is associated with a favorable tumor immune microenvironment of estrogen receptor (ER)-positive breast cancer patients which result in improved survival. Two major publicly available breast cancer cohorts; The Cancer Genome Atlas (TCGA) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) were used. The miR-143 high expression group was associated with increased infiltration of anti-cancer immune cells and decreased pro-cancer immune cells, as well as enrichment of the genes relating to T helper (Th1) cells resulting in improved overall survival (OS) in ER-positive breast cancer patients. To the best of our knowledge, this is the first study to demonstrate that high expression of miR-143 in cancer cells associates with a favorable tumor immune microenvironment, upregulation of anti-cancer immune cells, and suppression of the pro-cancer immune cells, associating with better survival of the breast cancer patients.

High Expression of miR-34a Associated with Less Aggressive Cancer Biology but Not with Survival in Breast Cancer

Most breast cancer (BC) patients succumb to metastatic disease. MiR-34a is a well-known tumor suppressive microRNA which exerts its anti-cancer functions by playing a role in p53, apoptosis induction, and epithelial-mesenchymal transition (EMT) suppression. Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) and The Cancer Genome Atlas (TCGA) cohorts were used to test our hypothesis that miR-34a high BCs translate to less aggressive cancer biology and better survival in large cohorts. There was no association between miR-34a expression levels and clinicopathological features of BC patients except for HER2 positivity. MiR-34a high expressing tumors were associated with lower Nottingham pathological grades and lower MKI67 expression. In agreement, high miR-34a tumors demonstrated lower GSVA scores of cell cycle and cell proliferation-related gene sets. High miR-34a tumors enriched the p53 pathway and apoptosis gene sets. Unexpectedly, high miR-34a tumors also associated with elevated EMT pathway score and ZEB1 and two expressions. MiR-34a expression did not associate with any distant metastasis. Further, high miR-34a tumors did not associate with better survival compared with miR-34a low tumors. In conclusion, the clinical relevance of miR-34a high expressing tumors was associated with suppressed cell proliferation, enhanced p53 pathway and apoptosis, but enhanced EMT and these findings did not reflect better survival outcomes in large BC patient cohorts.

Effects of MIR143 on rat sarcoma signaling networks in solid tumors: A brief overview

Rat sarcoma (RAS) is a well-known oncogene that plays important roles in cancer proliferation, cell survival and cell invasion. RAS exists as three major isoforms, Kirsten rat sarcoma (KRAS), Harvey rat sarcoma (HRAS) and neuroblastoma rat sarcoma (NRAS). Mutations of these genes account for approximately 30% of all cancers. Among them, KRAS mutations are the most common, responsible for 85%, followed by NRAS (12%) and HRAS (3%). Although the development of RAS inhibitors has been explored for over the past decade, so far, no effective inhibitor has been found. MicroRNA (miRNA) are a class of small non-coding RNA that control the gene expression of pleural target genes at the post-transcriptional level. MiRNA play critical roles in the physiological and pathological processes at work in cancers, such as cell proliferation, cell death, cell invasion and metastasis. MicroRNA-143 (MIR143) is known to function as a tumor suppressor in a variety of cancers. One of its known mechanisms is suppression of RAS expression and its effector signaling pathways, such as PI3K/AKT and MAPK/ERK. Within the last five years, we developed a potent chemically modified MIR143-3p that enabled us to elucidate the details of the KRAS signaling networks at play in colon and other cancer cells. In this review, we will discuss the role of MIR143-3p in those RAS signaling networks that are related to various biological processes of cancer cells. In addition, we will discuss the possibility of the use of MIR143 as a therapeutic drug for targeting RAS signaling networks.

KRAS signaling enriched triple negative breast cancer is associated with favorable tumor immune microenvironment and better survival

KRAS signaling is associated with cancer progression in several cancers. Upregulation of KRAS signaling is often seen in cancers that harbor high KRAS mutation rate, such as pancreatic cancer and non-small cell lung cancer (NSCLC). Less than 2% of breast cancers have KRAS mutation, however, the alteration of the effector signaling such as PI3K/AKT and MAPK pathways are well known. Mutated KRAS is known to function as immune suppressor in other cancers, but the role of KRAS signaling on tumor immune microenvironment (TIME) in breast cancer is not known. We hypothesize that the enrichment of KRAS signaling is associated with reduced patient survival as well as TIME in triple negative breast cancer (TNBC). Patient cohorts from Molecular Taxonomy of Breast Cancer International Consortium (METABRIC; n = 1903) and The Cancer Genome Atlas (TCGA; n = 982) were used. Higher expression of KRAS in breast cancer cell-lines (MCF7, BT474, and MDA-MB231) compared to MCF10A, which is a model of benign mammary cells was found. Both MEK and PI3K inhibitors suppressed MB231 cell proliferation in dose dependent manner. Gene Set Variant Analysis (GSVA) of the patient cohorts demonstrated two peaks by KRAS_SIGNALING_UP gene sets which were divided into KRAS-high and -low groups using median cutoff. There was no difference in KRAS mutation between KRAS-high and low. Despite its cell proliferation promoting role, KRAS-high patients demonstrated significantly better Disease-Free Survival and Overall Survival in triple negative breast cancer (TNBC). KRAS-high TNBC was associated with favorable tumor immune microenvironment with elevated B cells and CD8 T cells, monocytes, or M1 macrophage. It was associated with decreased CD4 central memory T-cells, but not Regulatory T-cells, or M2 macrophage detected by xCell. To elucidate the mechanism of this association, Gene Set Enrichment Analysis was performed. Inflammatory response, IL6/JAK-STAT3 signaling, and Interferon gamma response gene sets were enriched in KRAS-high TNBC patients in both METABRIC and TCGA cohorts. In agreement, cytolytic activity score, interferon gamma response score, and lymphocyte infiltrating signature score, were all significantly elevated in KRAS-high TNBC. In conclusion, we found that patients with enrichment of KRAS signaling gene sets were associated with inflammation and favorable tumor immune microenvironment as well as improved survival in TNBC.

Identification of a novel microRNA recurrence-related signature and risk stratification system in breast cancer

Increasing evidence has revealed that microRNAs (miRNAs) play vital roles in breast cancer (BC) prognosis. Thus, we aimed to identify recurrence-related miRNAs and establish accurate risk stratification system in BC patients. A total of 381 differentially expressed miRNAs were confirmed by analyzing 1044 BC tissues and 102 adjacent normal samples from The Cancer Genome Atlas (TCGA). Then, based on the association between each miRNAs and disease-free survival (DFS), we identified miRNA recurrence-related signature to construct a novel prognostic nomogram using Cox regression model. Target genes of the four miRNAs were analyzed via Gene Ontology and KEGG pathway analyses. Time-dependent receiver operating characteristic analysis indicated that a combination of the miRNA signature and tumor-node-metastasis (TNM) stage had better predictive performance than that of TNM stage (0.710 vs 0.616, P<0.0001). Furthermore, risk stratification analysis suggested that the miRNA-based model could significantly classify patients into the high- and low-risk groups in the two cohorts (all P<0.0001), and was independent of other clinical features. Functional enrichment analysis demonstrated that the 46 target genes mainly enrichment in important cell biological processes, protein binding and cancer-related pathways. The miRNA-based prognostic model may facilitate individualized treatment decisions for BC patients.

Triple-Negative Breast Cancer with High Levels of Annexin A1 Expression Is Associated with Mast Cell Infiltration, Inflammation, and Angiogenesis

Annexin A1 (ANXA1) is a phospholipid-linked protein involved in inflammation, immune response, and mast cell reactivity. Recently, we reported that ANXA1 is associated with aggressive features of triple-negative breast cancer (TNBC); however, its clinical relevance remains controversial. We hypothesized that human TNBC with high expression of ANXA1 mRNA is associated with pro-cancerous immune cell infiltration, including mast cells, and with an aggressive phenotype. Clinical and RNA-seq data were obtained from The Cancer Genome Atlas (TCGA, n = 1079) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) (n = 1904). TNBC patients had significantly higher levels of ANXA1 expression compared to the other subtypes in both TCGA and METABRIC cohorts (p < 0.001). ANXA1 protein expression was assessed by immunohistochemistry in Japanese TNBC patient cohort (n = 48), where 17 cases (35.4%) had positive ANXA1 staining, and their overall survival was significantly shorter compared with negative staining group (p = 0.008). The CIBERSORT algorithm was used to calculate immune cell infiltrations. ANXA1 high tumors were associated with activated mast cells and M2 macrophages (p > 0.01), but did not show any association with tumor heterogeneity nor cytolytic activity. High expression of ANXA1 group enriched inflammation, epithelial-to-mesenchymal transition (EMT), and angiogenesis-related genes in a gene set enrichment assay in both cohorts. To our knowledge, this is the first study to demonstrate that ANXA1 is associated with infiltration of mast cells and inflammation that is associated with the aggressive phenotype of TNBC, such as EMT and angiogenesis.

Non-coding RNAs as potential therapeutic targets in breast cancer

Paradigm shifting studies especially involving non-coding RNAs (ncRNAs) during last few decades have significantly changed the scientific perspectives regarding the complexity of cellular signalling pathways. Several studies have shown that the non-coding RNAs, initially ignored as transcriptional noise or products of erroneous transcription; actually regulate plethora of biological phenomena ranging from developmental processes to various diseases including cancer. Current strategies that are employed for the management of various cancers including that of breast fall short when their undesired side effects like Cancer Stem Cells (CSC) enrichment, low recurrence-free survival and development of drug resistance are taken into consideration. This review aims at exploring the potential role of ncRNAs as therapeutics in breast cancer, by providing a comprehensive understanding of their mechanism of action and function and their crucial contribution in regulating various aspects of breast cancer progression such as cell proliferation, angiogenesis, EMT, CSCs, drug resistance and metastasis. In addition, we also provide information about various strategies that can be employed or are under development to explore them as potential moieties that may be used for therapeutic intervention in breast cancer.

Tumor Heterogeneity Correlates with Less Immune Response and Worse Survival in Breast Cancer Patients

BACKGROUND

Intratumor heterogeneity implies that subpopulations of cancer cells that differ in genetic, phenotypic, or behavioral characteristics coexist in a single tumor (Ma in Breast Cancer Res Treat 162(1):39-48, 2017; Martelotto in Breast Cancer Res 16(3):210, 2014). Tumor heterogeneity drives progression, metastasis and treatment resistance, but its relationship with tumor infiltrating immune cells is a matter of debate, where some argue that tumors with high heterogeneity may generate neoantigens that attract immune cells, and others claim that immune cells provide selection pressure that shapes tumor heterogeneity (McGranahan et al. in Science 351(6280):1463-1469, 2016; McGranahan and Swanton in Cell 168(4):613-628, 2017). We sought to study the association between tumor heterogeneity and immune cells in a real-world cohort utilizing The Cancer Genome Atlas.

METHODS

Mutant allele tumor heterogeneity (MATH) was calculated to estimate intratumoral heterogeneity, and immune cell compositions were estimated using CIBERSORT. Survival analyses were demonstrated using Kaplan-Meir curves.

RESULTS

Tumors with high heterogeneity (high MATH) were associated with worse overall survival (p = 0.049), as well as estrogen receptor-positive (p = 0.011) and non-triple-negative tumors (p = 0.01). High MATH tumors were also associated with less infiltration of anti-tumor CD8 (p < 0.013) and CD4 T cells (p < 0.00024), more tumor-promoting regulatory T cells (p < 4e-04), lower expression of T-cell exhaustion markers, specifically PDL-1 (p = 0.0031), IDO2 (p = 0.34), ADORA2A (p = 0.018), VISTA (p = 0.00013), and CCR4 (p < 0.00001), lower expression of cytolytic enzymes granzyme A (p = 0.0056) and perforin 1 (p = 0.053), and low cytolytic activity score (p = 0.0028).

CONCLUSIONS

High heterogeneity tumors are associated with less immune cell infiltration, less activation of the immune response, and worse survival in breast cancer. Our results support the notion that tumor heterogeneity is shaped by selection pressure of tumor-infiltrating immune cells.

A novel six-microRNA-based model to improve prognosis prediction of breast cancer

Current tumor-node-metastasis (TNM) stage is unable to accurately predict the overall survival (OS) in breast cancer (BC) patients. This study aimed to construct a microRNA (miRNA)-based model to improve survival prediction of BC. We confirmed 99 differentially expressed miRNAs (DEMs) in 1044 BC samples compared to 102 adjacent normal breast tissues from The Cancer Genome Atlas (TCGA) database. Prognostic DEMs were used to establish a miRNA-based nomogram via Cox regression model. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes analyses (KEGG) were executed to analyze target genes of miRNAs. A six-miRNA signature was screened to effectively distinguish high-risk patients in the primary and validation cohort (all P<0.001). Furthermore, we established a novel prognostic model incorporating the six-miRNA signature and clinical risk factors to predict 5-year OS of BC. Time-dependent receiver operating characteristic analysis suggested that the predictive accuracy of the six-miRNA-based nomogram was distinctly higher than that of TNM stage (0.758 vs 0.650, P<0.001). GO and KEGG pathway analyses showed that the 39 target genes mainly enrichment in protein binding, cytoplasm and MAPK signaling pathway. Our six-miRNA-based model is a reliable prognostic tool for survival prediction and provides information for individualized treatment decisions in BC patients.