Identification of prognostic biomarkers for breast cancer based on miRNA and mRNA co-expression network

miRNA deregulation and relationship with metabolic parameters after Mediterranean dietary intervention in BRCA-mutated women

Background

Breast cancer onset is determined by a genetics-environment interaction. BRCA1/2 gene alterations are often genetically shared in familial context, but also food intake and hormonal assessment seem to influence the lifetime risk of developing this neoplasia. We previously showed the relationship between a six-months Mediterranean dietary intervention and insulin, glucose and estradiol levels in BRCA1/2 carrier subjects. The aim of the present study was to evidence the eventual influence of this dietary intervention on the relationship between circulating miRNA expression and metabolic parameters in presence of BRCA1/2 loss of function variants.

Methods

Plasma samples of BRCA-women have been collected at the baseline and at the end of the dietary intervention. Moreover, subjects have been randomized in two groups: dietary intervention and placebo. miRNA profiling and subsequent ddPCR validation have been performed in all the subjects at both time points.

Results

ddPCR analysis confirmed that five (miR-185-5p, miR-498, miR-3910, miR-4423 and miR-4445) of seven miRNAs, deregulated in the training cohort, were significantly up-regulated in subjects after dietary intervention compared with the baseline measurement. Interestingly, when we focused on variation of miRNA levels in the two timepoints, it could be observed that miR-4423, miR-4445 and miR-3910 expressions are positively correlated with variation in vitaminD level; whilst miR-185-5p difference in expression is related to HDL cholesterol variation.

Conclusions

We highlighted the synergistic effect of a healthy lifestyle and epigenetic regulation in BC through the modulation of specific miRNAs. Different miRNAs have been reported involved in the tumor onset acting as tumor suppressors by targeting tumor-associated genes that are often downregulated.

A review of prognostic and predictive biomarkers in breast cancer

Breast cancer (BC) is a common cancer all over the world that affects women. BC is one of the leading causes of cancer mortality in women, which today has decreased with the advancement of technology and new diagnostic and therapeutic methods. BCs are histologically divided into in situ and invasive carcinoma, and both of them can be divided into ductal and lobular. The main function after the diagnosis of invasive breast cancer is which patient should use chemotherapy, which patient should receive adjuvant therapy, and which should not. If the decision is for adjuvant therapy, the next challenge is to identify the most appropriate treatment or combination of treatments for a particular patient. Addressing the first challenge can be helped by prognostic biomarkers, while addressing the second challenge can be done by predictive biomarkers. Among the molecular markers related to BC, ER, PR, HER2, and the Mib1/Ki-67 proliferation index are the most significant ones and are tightly confirmed in the standard care of all primary, recurrent, and metastatic BC patients. CEA and CA-15-3 antigens are the most valuable markers of serum tumors in BC patients. Determining the series of these markers helps monitor response to the treatment and early detection of recurrence or metastasis. miRNAs have been demonstrated to be intricate in mammary gland growth, proliferation, and formation of BC known to be incriminated in BC biology. By combining established prognostic factors with valid prognostic/predicted biomarkers, we can start the journey to personalized treatment for every recently diagnosed BC patient.

A Comprehensive Analysis for Expression, Diagnosis, and Prognosis of m5C Regulator in Breast Cancer and Its ncRNA–mRNA Regulatory Mechanism

Recent studies have well demonstrated that 5-methylcytosine (m⁵C) regulators play pivotal roles in pathological conditions, including cancer. This study first tried to identify potential 5-methylcytosine (m⁵C) regulators in breast cancer by combination of expression, diagnosis, and survival analyses, and then established an ncRNA–mRNA network accounting for m⁵C regulators’ roles in breast cancer. Among 13 m⁵C regulators, DNMT3B and ALYREF were significantly upregulated in breast cancer and their high expression indicated unfavorable prognosis. Both DNMT3B and ALYREF possessed the statistical abilities to distinguish breast cancer from normal breast samples. Moreover, five potential upstream miRNAs (let-7b-5p, miR-195-5p, miR-29a-3p, miR-26a-5p, and miR-26b-5p) of m⁵C regulators could not only serve as independent prognostic predictors but also together made up a promising miRNA prognostic signature in breast cancer. Next, upstream potential lncRNAs of the five miRNAs were predicted and analyzed. Pathway enrichment analysis revealed that the target genes of these miRNAs were markedly enriched in some cancer-related pathways, and further investigation indicated VEGFA and EZH2 were found to be the most potential target genes in the m⁵C regulators-related ncRNA–mRNA network in breast cancer. These findings comprehensively provided key clues for developing m⁵C regulators-related effective therapeutic targets and promising diagnostic biomarkers in breast cancer.

A Novel TCGA-Validated, MiRNA-Based Signature for Prediction of Breast Cancer Prognosis and Survival

Breast cancer (BC) is the most common cancer affecting women and the leading cause of cancer-related deaths worldwide. Compelling evidence indicates that microRNAs (miRNAs) are inextricably involved in the development of cancer. Here, we constructed a novel model, based on miRNA-seq and clinical data downloaded from The Cancer Genome Atlas (TCGA). Data from a total of 962 patients were included in this study, and the relationships among their clinicopathological features, survival, and miRNA-seq expression levels were analyzed. Hsa-miR-186 and hsa-miR-361 were identified as internal reference miRNAs and used to normalize miRNA expression data. A five-miRNA signature, constructed using univariate and multivariate Cox regression, was significantly associated with disease-specific survival (DSS) of patients with BC. Kaplan–Meier (KM) and receiver operating characteristic (ROC) analyses were conducted to confirm the clinical significance of the five-miRNA signature. Finally, a nomogram was constructed based on the five-miRNA signature to evaluate its clinical value. Cox regression analysis revealed that a five-miRNA signature was significantly associated with DSS of patients with BC. KM analysis demonstrated that the signature could efficiently distinguish high- and low-risk patients. Moreover, ROC analysis showed that the five-miRNA signature exhibited high sensitivity and specificity in predicting the prognosis of patients with BC. Patients in the high-risk subgroup who received adjuvant chemotherapy had a significantly lower incidence of mortality than those who did not. A nomogram constructed based on the five-miRNA signature was effective in predicting 5-year DSS. This study presents a novel five-miRNA signature as a reliable prognostic tool to predict DSS and provide theoretical reference significance for individualized clinical decisions for patients with BC.

Circular RNA hsa_circ_0005046 and hsa_circ_0001791 May Become Diagnostic Biomarkers for Breast Cancer Early Detection

Breast cancer (BC) is one of the most common lethal diseases in women worldwide. Recent evidence has shown that covalently closed Circular RNA (circRNA) deregulation is observed in different human malignancies and cancers. Lately, circRNAs are being considered as a new diagnostic biomarker; however, the mechanism and the correlation of action between circRNAs and BC are still unclear. In the present study, we try to investigate the expression level of hsa_circ_0005046 and hsa_circ_0001791 in BC. By using quantitative real-time polymerase chain reaction (qRT-PCR), expression profiles of candidate circRNAs were detected in 60 BC tissue and paired adjacent normal tissues. Furthermore, the clinicopathological relation and diagnostic value were estimated. Our results showed the higher expression levels of hsa_circ_0005046 and hsa_circ_0001791 in BC tissues compared to paired adjacent normal tissues with P value (P < 0.0001) for both circRNAs, and the area under the receiver operating characteristic (ROC) curve was 0.857 and 1.0, respectively; in addition, a total 10 miRNAs that can be targeted by each candidate circRNAs was predicted base on bioinformatics databases. Taken together, for the first time, the results of our study presented high expression levels of hsa_circ_0005046 and hsa_circ_00017916 in BC; although there was no direct correlation between the high expression level of both circRNAs with clinic pathological factors, except hsa_circ_0001791 association with estrogen receptors (ER), high ROC curve in expressed samples indicated that both circRNAs could be used as a new diagnostic biomarker for BC. Moreover, miRNAs selection tools predicted that miR-215 and mir-383-5p which have a tumor suppressor role in BC can be targeted by our candidate circRNAs to affect the PI3K/AKT pathway; in conclusion, further studies are required to validate the oncogene role of our candidate circRNAs through the PI3k pathway.

Regulation of breast cancer metastasis signaling by miRNAs

Despite the decline in death rate from breast cancer and recent advances in targeted therapies and combinations for the treatment of metastatic disease, metastatic breast cancer remains the second leading cause of cancer-associated death in U.S. women. The invasion-metastasis cascade involves a number of steps and multitudes of proteins and signaling molecules. The pathways include invasion, intravasation, circulation, extravasation, infiltration into a distant site to form a metastatic niche, and micrometastasis formation in a new environment. Each of these processes is regulated by changes in gene expression. Noncoding RNAs including microRNAs (miRNAs) are involved in breast cancer tumorigenesis, progression, and metastasis by post-transcriptional regulation of target gene expression. miRNAs can stimulate oncogenesis (oncomiRs), inhibit tumor growth (tumor suppressors or miRsupps), and regulate gene targets in metastasis (metastamiRs). The goal of this review is to summarize some of the key miRNAs that regulate genes and pathways involved in metastatic breast cancer with an emphasis on estrogen receptor α (ERα+) breast cancer. We reviewed the identity, regulation, human breast tumor expression, and reported prognostic significance of miRNAs that have been documented to directly target key genes in pathways, including epithelial-to-mesenchymal transition (EMT) contributing to the metastatic cascade. We critically evaluated the evidence for metastamiRs and their targets and miRNA regulation of metastasis suppressor genes in breast cancer progression and metastasis. It is clear that our understanding of miRNA regulation of targets in metastasis is incomplete.

Cancer-Associated Stemness and Epithelial-to-Mesenchymal Transition Signatures Related to Breast Invasive Carcinoma Prognostic

Simple Summary

Breast cancer is one of the most common oncological diseases in women, as its incidence is rapidly growing. In this study, we have investigated the mechanism of epithelial-to-mesenchymal transition (EMT) and cancer stem cells (CSCs), demonstrating presence of an interconnectedness between them. This interconnectedness plays important roles in patient prognostic, as well as in diagnostic and therapeutic targets. It is identified that there is a common signature between CSCs and EMT, and this is represented by ALDH1A1, SFRP1, miR-139, miR-21, and miR-200c. This finding will provide a better understanding of this mechanism, and will facilitate the development of novel treatment options.

Abstract

Breast cancer is one of the most common oncological diseases in women, as its incidence is rapidly growing, rendering it unpredictable and causing more harm than ever before on an annual basis. Alterations of coding and noncoding genes are related to tumorigenesis and breast cancer progression. In this study, several key genes associated with epithelial-to-mesenchymal transition (EMT) and cancer stem cell (CSC) features were identified. EMT and CSCs are two key mechanisms responsible for self-renewal, differentiation, and self-protection, thus contributing to drug resistance. Therefore, understanding of the relationship between these processes may identify a therapeutic vulnerability that can be further exploited in clinical practice, and evaluate its correlation with overall survival rate. To determine expression levels of altered coding and noncoding genes, The Cancer Omics Atlas (TCOA) are used, and these data are overlapped with a list of CSCs and EMT-specific genes downloaded from NCBI. As a result, it is observed that CSCs are reciprocally related to EMT, thus identifying common signatures that allow for predicting the overall survival for breast cancer genes (BRCA). In fact, common CSCs and EMT signatures, represented by ALDH1A1, SFRP1, miR-139, miR-21, and miR-200c, are deemed useful as prognostic biomarkers for BRCA. Therefore, by mapping changes in gene expression across CSCs and EMT, suggesting a cross-talk between these two processes, we have been able to identify either the most common or specific genes or miRNA markers associated with overall survival rate. Thus, a better understanding of these mechanisms will lead to more effective treatment options.

Identifying of miR-98-5p/IGF1 axis contributes breast cancer progression using comprehensive bioinformatic analyses methods and experiments validation

BACKGROUND

Breast cancer (BC) is a huge health threat for women worldwide. Although numerous microRNAs (miRNA) have been found to be aberrantly expressed in BC, the construction of a comprehensive miRNA-messenger RNA (mRNA) network is still needed.

METHODS

Limma package was used to identify differentially expressed miRNAs (DEMs) in microarray datasets downloaded from GEO database. Genes targeted by DEMs were analyzed using mirTarBase. Gene Ontology and pathway enrichment analysis for these genes were performed at DAVID. Expression correlations of DEMs and target genes were analyzed at ENCORI. Based on these results, a miRNA-mRNA regulatory network was constructed.

RESULTS

A total of 17 overlapping DEMs were identified at these two microarray datasets. Expression of DEMs in BC tissues compared with normal tissues were further validated by ENCORI. By utilizing miRTarBase, a total of 167 target genes for DEMs were obtained. 10 hub genes (AKT1, MYC, VEGFA, CCND1, PTEN, IL6, CASP3, KRAS, IGF1, ESR1) were identified. Through analyzing the effects of hub genes on overall survival of BC patients and their expression correlation with miRNAs, we found hsa-miR-98-5p/IGF1 axis may play a crucial role in BC progression. The connections of hsa-miR-98-5p and IGF1 were further validated by luciferase activity reporter assay and functional assays.

CONCLUSIONS

In this work, a miRNA-mRNA network related to BC progression was built, and identified one important miRNA-mRNA axis in BC.

Comprehensive analysis of transcriptome data for identifying biomarkers and therapeutic targets in head and neck squamous cell carcinoma

Background

Head and neck squamous cell carcinoma (HNSCC) is one of the most common malignancy worldwide. Accumulating evidences have highlighted the importance of transcriptome data during HNSCC tumorigenesis. The aim of this study was to identify significant genes as effective biomarkers for HNSCC and constructed miRNA-mRNA regulatory network for a more comprehensive understanding of the underlying molecular mechanisms.

Methods

A total of four independent microarrays conducted on HNSCC samples were downloaded from the Gene Expression Omnibus (GEO) and analyzed through R software. FunRich was applied to predict potential transcription factors and targeted genes of miRNAs. Protein-protein interaction (PPI) network and miRNA-mRNA regulatory network were constructed in Cytoscape. Additionally, the database for annotation, visualization, and integrated discovery (DAVID) was utilized to perform GO and KEGG pathway enrichment analyses. Validation of gene expression levels was conducted by online databases and qPCR experiments.

Results

A total of 35 and 193 differentially expressed miRNAs (DEMs) and mRNAs (DEGs) were screened out by the limma package in R. The interactive network of the overlapping DEGs presented three significant modules and ten hub genes (FN1, MMP3, SPP1, STAT1, LOX, CXCL5, CXCL11, ISG15, IFIT3, and RSAD2). Predicted target genes of DEMs were visualized in Cytoscape and six miRNA-mRNA regulatory pairs were identified. Further validation demonstrated the upregulation of SLC16A1 and COL4A1 in HNSCC.

Conclusions

We performed an integrated and comprehensive bioinformatics analysis of miRNAs and mRNAs in HNSCC, contributing to explore the underlying regulatory mechanisms and to identify genetic biomarkers and therapeutic targets for HNSCC.

Identifying prognostic biomarkers in endometrial carcinoma based on ceRNA network

PURPOSE

Endometrial carcinoma (EC), a common gynecological malignancy with high incidence, affects the mental and physical health of women. Mounting evidence shows that long noncoding RNAs (lncRNAs), messenger RNAs (mRNAs), and microRNAs (miRNAs) have instrumental roles in various biological processes associated with the pathogenesis of EC. In this research, we intend to further study the mechanism of EC and the potential predictive markers of EC.

METHODS

First, we obtained original data of EC RNA transcripts from The Cancer Genome Atlas database and performed differential analysis. Subsequently, according to the miRcode online software, relationship pairs of lncRNA-miRNA were constructed, and miRNA-mRNA pairs were established based on miRDB, TargetScan, and miRTarBase. Then, we constructed the competing endogenous RNA (ceRNA) network based on lncRNA-miRNA and miRNA-mRNA pairs. To further explain the function of the ceRNA network and explore the potential prognostic markers, functional enrichment analysis, and survival analysis were carried out.

RESULTS

The research showed that there were 744 differential expression lncRNAs (DElncRNAs), 164 differential expression miRNAs (DEmiRNAs), and 2447 differential expression mRNAs (DEmRNAs) between EC tissues and normal tissues. Subsequently, we built 103 DEmiRNA-DEmRNA interaction pairs and 369 DElncRNA-DEmiRNA pairs. Then, we established the ceRNA network of EC, including 62 DElncRNAs, 26 DEmiRNAs, and 70 DEmRNAs. Moreover, 10 of 62 lncRNAs, 19 of 70 mRNAs, and 4 of 26 miRNAs that closely related to the survival of EC with P < .05 were obtained. Notably, based on this network, it was found that LINC00261-hsa-mir-31 pair and LINC00261-hsa-mir-211 target pairs could be used as the potential prognostic markers of EC.

CONCLUSION

This research recommended an available basis for the molecular mechanism of EC and prognosis prediction, which could help guide the subsequent treatments and predict the prognosis for patients with EC.