Integrated microRNA, gene expression and transcription factors signature in papillary thyroid cancer with lymph node metastasis

Targeting miR-31 represses tumourigenesis and dedifferentiation of BRAFV600E-associated thyroid carcinoma

BACKGROUND

BRAFV600E is the most common genetic mutation in differentiated thyroid cancer (DTC) occurring in 60% of patients and drives malignant tumour cell phenotypes including proliferation, metastasis and immune-escape. BRAFV600E-mutated papillary thyroid cancer (PTC) also displays greatly reduced expression of thyroid differentiation markers, thus tendency to radioactive iodine (RAI) refractory and poor prognosis. Therefore, understanding the molecular mechanisms and main oncogenic events underlying BRAFV600E will guide future therapy development.

METHODS

Bioinformatics and clinical specimen analyses, genetic manipulation of BRAFV600E-induced PTC model, functional and mechanism exploration guided with transcriptomic screening, as well as systematic rescue experiments were applied to investigate miR-31 function within BRAFV600E-induced thyroid cancer development. Besides, nanoparticles carrying miR-31 antagomirs were testified to alleviate 131I iodide therapy on PTC models.

RESULTS

We identify miR-31 as a significantly increased onco-miR in BRAFV600E-associated PTC that promotes tumour progression, metastasis and RAI refractoriness via sustained Wnt/β-catenin signalling. Mechanistically, highly activated BRAF/MAPK pathway induces miR-31 expression via c-Jun-mediated transcriptional regulation across in vitro and transgenic mouse models. MiR-31 in turn facilitates β-catenin stabilisation via directly repressing tumour suppressors CEBPA and DACH1, which direct the expression of multiple essential Wnt/β-catenin pathway inhibitors. Genetic functional assays showed that thyroid-specific knockout of miR-31 inhibited BRAFV600E-induced PTC progression, and strikingly, enhanced expression of sodium-iodide symporter and other thyroid differentiation markers, thus promoted 131I uptake. Nanoparticle-mediated application of anti-miR-31 antagomirs markedly elevated radio-sensitivity of BRAFV600E-induced PTC tumours to 131I therapy, and efficiently suppressed tumour progression in the pre-clinical mouse model.

CONCLUSIONS

Our findings elucidate a novel BRAF/MAPK-miR-31-Wnt/β-catenin regulatory mechanism underlying clinically BRAFV600E-associated DTC tumourigenesis and dedifferentiation, also highlight a potential adjuvant therapeutic strategy for advanced DTC.

Unraveling the Roles of miR-204-5p and HMGA2 in Papillary Thyroid Cancer Tumorigenesis

Thyroid cancer is the most common endocrine malignant tumor with an increasing incidence rate. Although differentiated types of thyroid cancer generally present good clinical outcomes, some dedifferentiate into aggressive and lethal forms. However, the molecular mechanisms governing aggressiveness and dedifferentiation are still poorly understood. Aberrant expression of miRNAs is often correlated to tumor development, and miR-204-5p has previously been identified in papillary thyroid carcinoma as downregulated and associated with aggressiveness. This study aimed to explore its role in thyroid tumorigenesis. To address this, gain-of-function experiments were performed by transiently transfecting miR-204-5p in thyroid cancer cell lines. Then, the clinical relevance of our data was evaluated in vivo. We prove that this miRNA inhibits cell invasion by regulating several targets associated with an epithelial-mesenchymal transition, such as SNAI2, TGFBR2, SOX4 and HMGA2. HMGA2 expression is regulated by the MAPK pathway but not by the PI3K, IGF1R or TGFβ pathways, and the inhibition of cell invasion by miR-204-5p involves direct binding and repression of HMGA2. Finally, we confirmed in vivo the relationship between miR-204-5p and HMGA2 in human PTC and a corresponding mouse model. Our data suggest that HMGA2 inhibition offers promising perspectives for thyroid cancer treatment.

Noncoding RNAs in Thyroid-Follicular-Cell-Derived Carcinomas

Simple Summary

Thyroid tumors represent the most common neoplastic pathology of the endocrine system. Mutations occurring in oncogenes and tumor suppressor genes are responsible for thyroid carcinogenesis; however, the complete mutational landscape characterizing these neoplasias has not been completely unveiled. It has been established that only the 2% of the human genome codes for proteins, suggesting that the vast majority of the genome has regulatory capabilities, which, if altered, could account for the onset of cancer. Hence, many scientific efforts are currently focused on the characterization of the heterogeneous class of noncoding RNAs, which represent an abundant part of the transcribed noncoding genome. In this review, we mainly focus on the involvement of microRNAs, long noncoding RNAs, and pseudogenes in thyroid cancer. The determination of the diagnosis, prognosis, and treatment of thyroid cancers based on the evaluation of the noncoding RNA network could allow the implementation of a more personalized approach to fighting these pathologies.

Abstract

Among the thyroid neoplasias originating from follicular cells, we can include well-differentiated carcinomas, papillary (PTC) and follicular (FTC) thyroid carcinomas, and the undifferentiated anaplastic (ATC) carcinomas. Several mutations in oncogenes and tumor suppressor genes have already been observed in these malignancies; however, we are still far from the comprehension of their full regulation-altered landscape. Even if only 2% of the human genome has the ability to code for proteins, most of the noncoding genome is transcribed, constituting the heterogeneous class of noncoding RNAs (ncRNAs), whose alterations are associated with the development of several human diseases, including cancer. Hence, many scientific efforts are currently focused on the elucidation of their biological role. In this review, we analyze the scientific literature regarding the involvement of microRNAs (miRNAs), long noncoding RNAs (lncRNAs), and pseudogenes in FTC, PTC, and ATC. Recent findings emphasized the role of lncRNAs in all steps of cancer progression. In particular, lncRNAs may control progression steps by regulating the expression of genes and miRNAs involved in cell proliferation, apoptosis, epithelial–mesenchymal transition, and metastatization. In conclusion, the determination of the diagnosis, prognosis, and treatment of cancer based on the evaluation of the ncRNA network could allow the implementation of a more personalized approach to fighting thyroid tumors.

MicroRNAs in Papillary Thyroid Cancer: What Is New in Diagnosis and Treatment

Introduction

Papillary thyroid cancer (PTC) accounts for up to 80% of thyroid malignancies. New diagnostic and therapeutic options are suggested including innovative molecular methods. MicroRNAs (miRNAs) are nonprotein coding single-stranded RNAs that regulate many cell processes. The aim of the present study is to review the deregulated miRNAs associated with PTCs.

Methods

A bibliographic research was conducted, resulting in 272 articles referred to miRNAs and PTC. Regarding our exclusion criteria, 183 articles were finally included in our review.

Results

A remarkably large number of miRNAs have been found to be deregulated during PTC manifestation in the literature. The deregulated miRNAs are detected in tissue samples, serum/plasma, and FNA samples of patients with PTC. These miRNAs are related to several molecular pathways, involving genes and proteins responsible for important biological processes. MiRNA deregulation is associated with tumor aggressiveness, including larger tumor size, multifocality, extrathyroidal extension, lymphovascular invasion, lymph node and distant metastasis, and advanced tumor node metastasis stage.

Conclusion

MiRNAs are proposed as new diagnostic and therapeutic tools regarding PTC. They could be essential biomarkers for PTC diagnosis applied in serum and FNA samples, while their contribution to prognosis is of great importance.

Construction of a Tumor Immune Microenvironment-Related Prognostic Model in BRAF-Mutated Papillary Thyroid Cancer

BRAF mutation is a representative oncogenic mutation, with a frequency of 60% in papillary thyroid carcinoma (PTC), but the reasons for the poor prognosis and more aggressive course of BRAF-mutated PTC are controversial. Tumor immune microenvironment (TIME) is an essential factor permitting the development and progression of malignancy, but whether TIME participates in the prognosis of BRAF-mutated PTC has not yet been reported. The primary goal of the present study was to provide a comprehensive TIME-related prognostic model to increase the predictive accuracy of progression-free survival (PFS) in patients with BRAF-mutated PTC. In this study, we analyzed the mRNA-seq data and corresponding clinical data of PTC patients obtained from the TCGA database. By calculating the TIME scores (immune score, stromal score and ESTIMATE score), the BRAF mutation group (n=237) was dichotomized into the high- and low-score groups. By functional analysis of differentially expressed genes (DEGs) in different high/low score groups, we identified 2 key TIME-related genes, HTR3A and NIPAL4, which affected PFS in BRAF-mutated PTC. A risk scoring system was developed by multivariate Cox analysis based on the abovementioned 2 TIME-related genes. Then, the BRAF-mutated cohort was divided into the high- and low-risk groups using the median risk score as a cutoff. A high risk score correlated positively with a higher HTR3A/NIPAL4 expression level but negatively with PFS in BRAF-mutated PTC. Ultimately, a nomogram was constructed by combining risk score with clinical parameter (Tumor stage), and the areas under the ROC curve (AUCs) of the nomogram for predicting 1-, 3- and 5-year PFS were then calculated and found to be 0.694, 0.707 and 0.738, respectively, indicating the improved accuracy and clinical utility of the nomogram versus the risk score model in the BRAF-mutated PTC cohort. Moreover, we determined the associations between prognostic genes or risk score and immune cell infiltration by two-way ANOVA. In the high-risk score, high HTR3A expression, and high NIPAL4 expression groups, higher infiltration of immune cells was found. Collectively, these findings confirm that the nomogram is effective in predicting the outcome of BRAF-mutated PTC and will add a spatial dimension to the developing risk stratification system.

Identification of Transcriptional Pattern Related to Immune Cell Infiltration With Gene Co-Expression Network in Papillary Thyroid Cancer

BACKGROUND

A growing body of evidence suggests that immune cell infiltration in cancer is closely related to clinical outcomes. However, there is still a lack of research on papillary thyroid cancer (PTC).

METHODS

Based on single-sample gene set enrichment analysis (SSGSEA) algorithm and weighted gene co-expression network analysis (WGCNA) tool, the infiltration level of immune cell and key modules and genes associated with the level of immune cell infiltration were identified using PTC gene expression data from The Cancer Genome Atlas (TCGA) database. In addition, the co-expression network and protein-protein interactions network analysis were used to identify the hub genes. Moreover, the immunological and clinical characteristics of these hub genes were verified in TCGA and GSE35570 datasets and quantitative real-time polymerase chain reaction (PCR). Finally, receiver operating characteristic (ROC) curve analysis was used to evaluate the diagnostic value of hub genes.

RESULTS

Activated B cell, activated dendritic cell, CD56bright natural killer cell, CD56dim natural killer cell, Eosinophil, Gamma delta T cell, Immature dendritic cell, Macrophage, Mast cell, Monocyte, Natural killer cell, Neutrophil and Type 17 T helper cell were significantly changed between PTC and adjacent normal groups. WGCNA results showed that the black model had the highest correlation with the infiltration level of activated dendritic cells. We found 14 hub genes whose expression correlated to the infiltration level of activated dendritic cells in both TCGA and GSE35570 datasets. After validation in the TCGA dataset, the expression level of only 5 genes (C1QA, HCK, HLA-DRA, ITGB2 and TYROBP) in 14 hub genes were differentially expressed between PTC and adjacent normal groups. Meanwhile, the expression levels of these 5 hub genes were successfully validated in GSE35570 dataset. Quantitative real-time PCR results showed the expression of these 4 hub genes (except C1QA) was consistent with the results in TCGA and GSE35570 dataset. Finally, these 4 hub genes had diagnostic value to distinguish PTC and adjacent normal controls.

CONCLUSIONS

HCK, HLA-DRA, ITGB2 and TYROBP may be key diagnostic biomarkers and immunotherapy targets in PTC.

CNV Analysis of the Correlation between Preoperative Lymph Node Metastasis and Prognosis of Early Tongue Cancer

Objective: To investigate the relationship between preoperative concealed lymph node metastasis (LNM) and prognosis in patients with early tongue cancer.

Methods: According to preoperative lymph node dissection, 41 patients with tongue cancer were divided into the LNM group (n=19) and the non-LNM group (n=22). Analysis of single nucleotide variation (SNV), tumor mutation burden (TMB), mutant allele tumor heterogeneity (MATH), aneuploidy and copy number variation (CNV) was performed to identify differentially expressed genes (DEGs) related to LNM. While KEGG analysis was conducted to reveal the CNV differentially expressed genes and main enriched pathways, the correlation between pathways and genes was analyzed by hierarchical clustering. The clinical information of LNM and data on overall survival (OS) rate were obtained from The Cancer Genome Atlas (TCGA), and survival analysis was performed based on combined LNM information.

Results: We observed significant correlations for the mTOR signaling pathway, Hippo signaling pathway and Wnt signaling pathway with the LNM group, while leukocyte transendothelial migration and cytokine-cytokine receptor interaction were markedly correlated with the non-LNM group. Moreover, TNFRSF10C was identified as the main DEG in the leukocyte transendothelial migration pathway. TCGA-based validation revealed that the disease-free survival (DFS) of the non-LNM group was significantly higher than that of the LNM group (P<0.005) when TNFRSF10C CNV was set to a log2 ratio>-0.163.

Conclusion: The differential expression of TNFRSF10C in leukocyte transendothelial migration, an immune-related pathway, is associated with LNM and DFS. The TNFRSF10C CNV log2 ratio could potentially serve as an indicator of good prognosis for tongue cancer patients with non-LNM after neck dissection.

Next-Generation Biomarkers in Multiple Myeloma: Understanding the Molecular Basis for Potential Use in Diagnosis and Prognosis

Multiple myeloma (MM) is considered to be the second most common blood malignancy and it is characterized by abnormal proliferation and an accumulation of malignant plasma cells in the bone marrow. Although the currently utilized markers in the diagnosis and assessment of MM are showing promising results, the incidence and mortality rate of the disease are still high. Therefore, exploring and developing better diagnostic or prognostic biomarkers have drawn global interest. In the present review, we highlight some of the recently reported and investigated novel biomarkers that have great potentials as diagnostic and/or prognostic tools in MM. These biomarkers include angiogenic markers, miRNAs as well as proteomic and immunological biomarkers. Moreover, we present some of the advanced methodologies that could be utilized in the early and competent diagnosis of MM. The present review also focuses on understanding the molecular concepts and pathways involved in these biomarkers in order to validate and efficiently utilize them. The present review may also help in identifying areas of improvement for better diagnosis and superior outcomes of MM.

The combination of circRNA-UMAD1 and Galectin-3 in peripheral circulation is a co-biomarker for predicting lymph node metastasis of thyroid carcinoma

The diagnosis of lymph node metastasis (LNM) by liquid biopsy is a novel concept prompted by the necessity to develop a more convenient and accurate method to guide the clinical management of early LNM in papillary thyroid carcinoma (PTC). However, the sensitivity and specificity of many biomarkers are not high enough. We aimed to detect circRNAs from peripheral circulation that may be better associated with the prognosis of LNM in PTC. First, Galectin-3 (Gal3) in blood was determined to be highly expressed in LNM patients. Second, based on a bioinformatics analysis and miRNA sequencing analysis from 2 paired primary and LNM tumors, miR-873 was identified to directly target Gal3, which was significantly associated with clinical parameters including LNM. Third, from additional circRNA sequencing, circRNA-UMAD1 was selected as a specific sponge for miR-873 and was correlated with Gal3 levels in peripheral circulation. Fourth, circRNA-UMAD1 and Gal3 were identified to have stronger co-biomarker potential with relatively high expression in the serum of LNM patients compared with primary tumor patients, as demonstrated by the RNA expression levels in the serum of 50 PTC patients with or without LNM by quantitative real-time PCR. Overall, the combination of circRNA-UMAD1 and Gal3 is a useful and effective co-biomarker for the prognosis of LNM in PTC patients. This new molecular typing method for LNM in PTC is more precise.

hsa-miR-195-5p inhibits cell proliferation of human thyroid carcinoma cells via modulation of p21/cyclin D1 axis

Background

Based on existing evidence, microRNAs (miRs) are gene regulators that undertake key functions in the oncogenesis and tumor progression of every single human malignant disease, such as thyroid carcinoma (TC). Previous clinical findings showed that expression of miR-195 is down-regulated in TC, which implies that miR-195 may be practically involved in TC pathogenesis. Nevertheless, the function of hsa-miR-195-5p in TC is still largely unclear. Herein, we detected the conceivable involvement of hsa-miR-195-5p in TC cell proliferation.

Methods

Real time PCR examination was performed to assess the expression level of hsa-miR-195-5p in TC cell lines TPC-1 and B-CPAP. TPC-1 cells were transfected with either hsa-miR-195-5p mimics or hsa-miR-195-5p inhibitor. After confirmation of transfection efficiency, the effect of hsa-miR-195-5p on proliferation and cell cycle of TPC-1 cells was assessed. The expression of cyclin D1 and p21 was simultaneously detected by western blotting. Moreover, targetScan 6.2 was used to predict hsa-miR-195-5p target genes. Subsequently, luciferase reporter was performed to examine whether there is a possible binding of hsa-miR-195-5p to 3’-UTR of cyclin D1 mRNA. Furthermore, cyclin D1 mRNA and protein levels were measured to check whether hsa-miR-195-5p exerts its function at the post-transcriptional level. In addition, to explore the function of cyclin D1 in TPC-1 cells overexpressing hsa-miR-195-5p, cyclin D1 siRNA was used to silence the expression of cyclin D1 in TPC-1 cells overexpressing hsa-miR-195-5p.

Results

We quantified the expression of hsa-miR-195-5p in TC cells and normal thyroid cells and found a remarkable decrease in hsa-miR-195-5p expression in TC cells. Over-expression of hsa-miR-195-5p obviously resulted in downgraded proliferation of TC cells. Moreover, hsa-miR-195-5p caused cell arrest at the GO/G1 phase. Further in silico analyses and the dual-luciferase reporter assay confirmed that 3’-UTR of cyclin D1 is a direct target of hsa-miR-195-5p. Western blot analysis uncovered that hsa-miR-195-5p over-expression led to decreased levels of cyclin D1 and p21. In mechanistic analyses, we found that silencing of cyclin D1 reversed the inhibitory effect of hsa-miR-195-5p on the proliferation of TC cells, which indicates that hsa-miR-195-5p suppresses TC cell proliferation by adversely regulating cyclin D1.

Conclusions

We concluded that hsa-miR-195-5p is a candidate tumor-suppressor miRNA in TC and that the hsa-miR-195-5p/p21/cyclin D1 pathway could be a potential therapeutic target for TC.

Differential MicroRNA-Signatures in Thyroid Cancer Subtypes

Thyroid cancer is one of the most common endocrine cancers, with an increasing trend in the last few decades. Although papillary thyroid cancer is the most frequent subtype compared with follicular or anaplastic thyroid cancer, it can dedifferentiate to a more aggressive phenotype, and the recurrence rate is high. The cells of follicular adenomas and follicular carcinomas appear identical in cytology, making the preoperative diagnosis difficult. On the other hand, anaplastic thyroid cancer poses a significant clinical challenge due to its aggressive nature with no effective therapeutic options. In the past several years, the roles of genetic alterations of thyroid tumors have been documented, with a remarkable correlation between genotype and phenotype, indicating that distinct molecular changes are associated with a multistep tumorigenic process. Besides mRNA expression profiles, small noncoding microRNA (miRNA) expression also showed critical functions for cell differentiation, proliferation, angiogenesis, and resistance to apoptosis and finally activating invasion and metastasis in cancer. Several high-throughput sequencing studies demonstrate that miRNA expression signatures contribute clinically relevant information including types of thyroid cancer, tumor grade, and prognosis. This review summarizes recent findings on miRNA signatures in thyroid cancer subtypes.

Integrin β4 promotes invasion and anoikis resistance of papillary thyroid carcinoma and is consistently overexpressed in lymphovascular tumor thrombus

Although the majority of papillary thyroid cancers (PTC) are indolent, a subset of PTCs behaves aggressively due to extensive invasion and distant metastasis. Integrin β4, a member of the integrin family, has been shown to enhance the progression in some malignancies; however, its role in PTC remains unclear. Here, we demonstrated that β4 overexpression was associated with extrathyroid extension, lymph node metastasis, high TNM stage, and poor overall survival based on The Cancer Genome Atlas cohort. Immunohistochemistry showed that β4 expression was significantly upregulated in the tumors with infiltrating growth pattern, as well as those with positive lymphovascular invasion. Moreover, β4 was invariably overexpressed in the lymphovascular tumor thrombi, which has not been reported before. After shRNA-induced knockdown of β4 in vitro, the migration, invasion and scratch repair ability of the tumor cells were significantly reduced. Furthermore, β4 reduction decreased anchorage-independent growth and increased anoikis. The bioinformatics analysis revealed that approximately 70 pathways were significantly dysregulated in the high β4 expression group. The MAPK pathway and propanoate metabolism were located in the network center of those pathways. Taken together, our results suggest that β4 could promote the tumor's aggressiveness by enhancing invasion and antagonizing anoikis. The upregulated expression of β4 in the tumor thrombi is intrinsically linked to its role in strengthening the anoikis resistance.

Identification of lncRNAs associated with lung squamous cell carcinoma prognosis in the competitive endogenous RNA network

Background

Long noncoding RNAs (lncRNAs) play a role in the formation, development, and prognosis of various cancers. Our study aimed to identify prognostic-related lncRNAs in lung squamous cell carcinoma (LUSC), which may provide new perspectives for individualized treatment of patients.

Materials and Methods

The RNA sequencing (lncRNA, microRNA (miRNA), mRNA) data and clinical information related to LUSC were obtained from The Cancer Genome Atlas (TCGA) database. Differentially expressed RNA sequences were used to construct the competitive endogenous RNA (ceRNA) network. In present study, we mainly used two prognostic verification methods, Cox analysis and survival analysis, to identify the prognostic relevance of specific lncRNAs and construct prognostic model of lncRNA.

Results

Datasets on 551 samples of lncRNA and mRNA and 523 miRNA samples were retrieved from the TCGA database. Analysis of the normal and LUSC samples identified 170 DElncRNAs, 331 DEmiRNAs, and 417 DEmRNAs differentially expressed RNAs. The ceRNA network contained 27 lncRNAs, 43 miRNAs, and 11 mRNAs. Furthermore, we identified seven specific lncRNAs (ERVH48-1, HCG9, SEC62-AS1, AC022148.1, LINC00460, C5orf17, LINC00261) as potential prognostic factors after correlation analysis, and five of the seven lncRNAs (AC022148.1, HCG9, LINC00460, C5orf17, LINC00261) constructed a prognostic model of LUSC.

Conclusion

In present study, we identified seven lncRNAs in the ceRNA network that are associated with potential prognosis in LUSC patients, and constructed a prognostic model of LUSC which can be used to assess the prognosis risk of clinical patients. Further biological experiments are needed to elucidate the specific molecular mechanisms underlying them.

Comprehensive analysis of the aberrantly expressed lncRNA‑associated ceRNA network in breast cancer

Previous studies have suggested that long non-coding RNAs (lncRNAs) are closely associated with human diseases, particularly cancer, including cancer of the lung, breast and stomach. A variety of lncRNAs are abnormally expressed in cancer and participate in several pathways including cell proliferation and apoptosis; these elements are closely associated with the development of cancer. The Cancer Genome Atlas (TCGA) is an important cancer database. It consists of clinical data, genomic variation, mRNA, microRNA (miRNA) and lncRNAs expression, methylation and other data for various types of human cancer. In the present study, differential expression of RNA was identified using the edgeR package. A total 1,222 RNA sequencing profiles from patients with breast cancer were downloaded from TCGA. A competing endogenous RNA (ceRNA) network was constructed for breast cancer based on miRcode and miRTarBase. The top 10 lncRNAs were selected using Cox regression analysis. Survival analysis was performed using Kaplan-Meier analysis. A total of 1,028 breast cancer-associated lncRNAs and 89 miRNAs (fold change >2; P<0.05) were identified; among these, 93 lncRNAs and 19 miRNAs were included in the ceRNA network. Subsequently, 10 basic lncRNAs were selected and their associations with overall survival were identified. In addition, 5 lncRNAs (ADAM metallopeptidase with thrombospondin type 1 motif 9-antisense RNA 1, AL513123.1, chromosome 10 open reading frame 126, long intergenic non-protein coding RNA 536 and Wilms tumor 1 antisense RNA) were identified to be significantly associated with overall survival (P<0.05, log rank test). These results suggested that mRNAs, lncRNAs and miRNAs were involved in pathological mechanisms of breast cancer. The newly-identified ceRNA network included 93 breast cancer-specific lncRNAs, 19 miRNAs and 27 mRNAs. The results of the present study highlight the potential of lncRNAs in understanding the development and pathogenesis of breast cancer, and suggest novel concepts and an experimental basis for the identification of prognostic biomarkers and therapeutic targets for breast cancer.

Development of a prognostic index based on an immunogenomic landscape analysis of papillary thyroid cancer

BACKGROUND

Papillary thyroid cancer (PTC) is the most common subtype of thyroid cancer, and inflammation relates significantly to its initiation and prognosis. Systematic exploration of the immunogenomic landscape therein to assist in PTC prognosis is therefore urgent. The Cancer Genome Atlas (TCGA) project provides a large number of genetic PTC samples that enable a comprehensive and reliable immunogenomic study.

METHODS

We integrated the expression profiles of immune-related genes (IRGs) and progression-free intervals (PFIs) in survival in 493 PTC patients based on the TCGA dataset. Differentially-expressed and survival-associated IRGs in PTC patients were estimated a computational difference algorithm and COX regression analysis. The potential molecular mechanisms and properties of these PTC-specific IRGs were also explored with the help of computational biology. A new prognostic index based on immune-related genes was developed by using multivariable COX analysis.

RESULTS

A total of 46 differentially expressed immune-related genes were significantly correlated with clinical outcome of PTC patients. Functional enrichment analysis revealed that these genes were actively involved in a cytokine-cytokine receptor interaction KEGG pathway. A prognostic signature based on RGs (AGTR1, CTGF, FAM3B, IL11, IL17C, PTH2R and SPAG11A) performed moderately in prognostic predictions and correlated with age, tumor stage, metastasis, number of lesions, and tumor burden. Intriguingly, the prognostic index based on IRGs reflected infiltration by several types of immune cells.

CONCLUSIONS

Together, our results screened several IRGs of clinical significance, revealed drivers of the immune repertoire, and demonstrated the importance of a personalized, IRG-based immune signature in the recognition, surveillance, and prognosis of PTC.

MicroRNA-705 regulates the differentiation of mouse mandible bone marrow mesenchymal stem cells

The craniofacial skeleton is the foundation of most stomatological treatments, including prosthodontics and maxillofacial surgery. Although histologically similar to the appendicular skeleton, the craniofacial skeleton manifests many unique properties in response to external stimuli and signals. However, the mandibular or maxillary bone marrow mesenchyme, which is the intrinsic foundation of the functions of craniofacial skeleton, has not been well studied, and its homeostasis mechanism remains elusive. Osteoporosis is a systemic disease that affects all skeletons and is characterized by bone mass loss. Osteoporotic bone marrow mesenchymal stem cells (BMMSCs) exhibit disturbed homeostasis and distorted lineage commitment. Many reports have shown that microRNAs (miRNAs) play important roles in regulating MSCs homeostasis. Here, to obtain a better understanding of mandibular bone marrow MSCs homeostasis, we isolated and cultured mandible marrow MSCs from mouse mandibles. Using miR-705 mimics and an inhibitor, we demonstrated that miR-705 played a vital role in shifting the mandibular MSCs lineage commitment in vitro. Utilizing an osteoporosis mouse model, we demonstrated that MSCs from ovariectomized (OVX) mouse mandibular bone marrow exhibited impaired osteogenic and excessive adipogenic differentiation. miR-705 was found overexpressed in OVX mandibular MSCs. The knock down of miR-705 in vitro partially attenuated the differentiation disorder of the OVX mandibular MSCs by upregulating the expression of osteogenic marker genes but suppressing adipogenic genes. Taken together, our findings provide a better understanding of the homeostasis mechanism of mandibular BMMSCs and a novel potential therapeutic target for treating mandibular osteoporosis.

The construction and analysis of ceRNA networks in invasive breast cancer: a study based on The Cancer Genome Atlas

Background

Studies have shown that long noncoding RNAs (lncRNAs) make up the major proportion of the ceRNA network and can regulate gene expression by competitively binding to miRNAs. This reveals the existence of an RNA-miRNA regulatory pathway and is of great biological significance. CeRNAs, as competitive endogenous RNAs, have revealed a new mechanism of interaction between RNAs. Until now, the role of lncRNA-mediated ceRNAs in breast cancer and their regulatory mechanisms have been elucidated to some extent.

Purpose

In this study, comprehensive analysis of large-scale invasive breast cancer samples in TCGA were conducted to further explore the developmental mechanism of invasive breast cancer and the potential predictive markers for invasive breast cancer prognosis in the ceRNA network.

Methods

Abnormal expression profiles of invasive breast cancer associated mRNAs, lncRNAs and miRNAs were obtained from the TCGA database. Through further alignment and prediction of target genes, an abnormal lncRNA-miRNA-mRNA ceRNA network was constructed for invasive breast cancer. Through the overall survival analysis, Identification prognostic bio-markers for invasive breast cancer patients. In addition, we used Cytoscape plug-in BinGo for the different mRNA performance functional cluster analysis.

Results

Differential analysis revealed that 1059 lncRNAs, 86 miRNAs, and 2138 mRNAs were significantly different in invasive breast cancer samples versus normal samples. Then we construct an abnormal lncRNA-miRNA-mRNA ceRNA network for invasive breast cancer, consisting of 90 DElncRNAs, 18 DEmiRNAs and 26 DEmRNAs.Further, 4 out of 90 lncRNAs, 3 out of 26 mRNAs, and 2 out of 18 miRNAs were useful as prognostic biomarkers for invasive breast cancer patients (P value < 0.05). It is worth noting that based on the ceRNA network, we found that the LINC00466-Hsa-mir-204- NTRK2 LINC00466-hsa-mir-204-NTRK2 axis was present in 9 RNAs associated with the prognosis of invasive breast cancer.

Conclusion

This study provides an effective bioinformatics basis for further understanding of the molecular mechanism of invasive breast cancerand for predicting outcomes, which can guide the use of invasive breast cancerdrugs and subsequent related research.

Proteomics study of serum exosomes from papillary thyroid cancer patients

Lymph node metastasis (LNM) in papillary thyroid cancer (PTC) is related to increased risk of recurrence and poor prognosis. Tumour exosomes have been shown to be associated with metastasis of cancer cells. Therefore, we aim to identify the characteristics and biological functions of serum exosomes in lymph node metastases of PTC. We compared proteome profiles of serum-purified exosomes (SPEs) from PTC patients with LNM, PTC patients without LNM, and healthy donors, using a combination of liquid chromatography-tandem mass spectroscopy analyses and tandem mass tag label quantitation analysis. We identified 1569 proteins by two or more unique peptides. Compared with the SPEs of PTC patients without LNM, we found 697 differentially expressed proteins in the SPEs of PTC patients with LNM. Our results revealed overexpression of specific proteins with well-established links to cancer cell metastasis, such as SRC, TLN1, ITGB2 and CAPNS1. Consistent with mass spectrum results, we performed Western blot to detect the expression of these proteins in individual sample. Biological pathway analyses showed that integrin signalling was aberrantly activated in the SPEs of PTC patients with LNM compared to those without LNM. Our study reveals that SPEs of PTC patients with lymph node metastases promote BHT101 thyroid cancer cell invasiveness, but have no apparent influence on cell migration. In the serum exosomes of PTC patients with LNM, integrin-associated proteins are obviously upregulated. These proteomic findings will contribute to elucidation of the pathophysiological functions of tumour-derived exosomes.

A potential biomarker hsa-miR-200a-5p distinguishing between benign thyroid tumors with papillary hyperplasia and papillary thyroid carcinoma

Papillary thyroid carcinoma (PTC) is the most common endocrine cancer with a significantly increase of the incidence recently. Several cytokines, such as thyroid peroxidase (TPO), cluster of differentiation 56 (CD56), Galectin-3, mesothelial cell (MC), cytokeratin 19 (CK19) and BRAF (B-raf) were recommended to be tested by immunohistochemistry (IHC) for a definitive diagnosis, but were still limited in clinical use because of their relative lower sensitivity and specificity. MicroRNA (miRNA), as a new molecular biomarkers, however, has not been reported yet so far. To address this, hsa-miR-200a-5p, a miRNA, was selected and detected in PTC patients by in situ hybrization with benign thyroid tumor with papillary hyperplasia as a control, and the differential expression of hsa-miR-200a-5p between fresh PTC tissues and control was detected by qRT-PCR. Expressive levels of cytokines of TPO, CD56, Galectin-3, MC, CK19 and B-raf were also detected by immunohistochemistry. The correlation was analyzed by SPSS software using Spearman methods. As expected, the hsa-miR-200a-5p expressive level was significantly increased in PTC patients, compared to that of control, and was consistent with that of TPO, CD56, Galectin-3, MC, CK19 and B-raf. In addition, expression of hsa-miR-200a-5p showed negative correlation to that of TPO (rs = - 0.734; **: P < 0.01) and CD56 (rs = - 0.570; **: P < 0.01), but positive correlation to that of Galectin-3 (rs = 0.601; **: P < 0.01), MC (rs = 0.508; **: P < 0.01), CK19 (rs = 0.712; **: P < 0.01) and B-raf (rs = 0.378; **: P < 0.01). PTC and papillary benign thyroid papillary hyperplasia are difficult to distinguish in morphology, so requiring immunohistochemistry to further differentiate the diagnosis, however, for the existing clinical common diagnostic marker for immunohistochemistry, the sensitivity and accuracy are low, it is easy to miss diagnosis. Therefore, there is an urgent need for a rapid and sensitive molecular marker. So miR-200a-5p can be used to assist in the diagnosis of PTC at the molecular level, and as a biomarker, can be effectively used to distinguish between PTC and benign thyroid tumor with papillary hyperplasia.

MicroRNA-4728 mediated regulation of MAPK oncogenic signaling in papillary thyroid carcinoma

Papillary thyroid carcinoma (PTC) is the most common type of thyroid cancer that accounts for 85% of thyroid cancers. MicroRNAs (miRNAs) have been reported to play important roles in the biological processes in cancer. In this study, we analyzed the biological role of miR-4728 in human PTC process in human PTC cell lines in vitro. MiRNA-4728 was observed to down-regulated in human PTC tissues and PTC cell lines. Additionally, miR-4728 inhibited PTC cell proliferation. Further study demonstrated SOS1 was repressed by miR-4728 and overexpression of miR-4728 down-regulated both the mRNA and protein levels of SOS1. Moreover, miR-4728 overexpression also decreased the MAPK signaling activity. These observations suggested that miR-4728 could inhibit the process of human PTC through regulating MAPK signaling pathway. And, appropriate regulation of miR-4728 might be vital to improve human PTC treatment.

Development of a computational promoter with highly efficient expression in tumors

Background

Gene therapy is a potent method to increase the therapeutic efficacy against cancer. However, a gene that is specifically expressed in the tumor area has not been identified. In addition, nonspecific expression of therapeutic genes in normal tissues may cause side effects that can harm the patients’ health. Certain promoters have been reported to drive therapeutic gene expression specifically in cancer cells; however, low expression levels of the target gene are a problem for providing good therapeutic efficacy. Therefore, a specific and highly expressive promoter is needed for cancer gene therapy.

Methods

Bioinformatics approaches were utilized to analyze transcription factors (TFs) from high-throughput data. Reverse transcription polymerase chain reaction, western blotting and cell transfection were applied for the measurement of mRNA, protein expression and activity. C57BL/6JNarl mice were injected with pD5-hrGFP to evaluate the expression of TFs.

Results

We analyzed bioinformatics data and identified three TFs, nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB), cyclic AMP response element binding protein (CREB), and hypoxia-inducible factor-1α (HIF-1α), that are highly active in tumor cells. Here, we constructed a novel mini-promoter, D5, that is composed of the binding sites of the three TFs. The results show that the D5 promoter specifically drives therapeutic gene expression in tumor tissues and that the strength of the D5 promoter is directly proportional to tumor size.

Conclusions

Our results show that bioinformatics may be a good tool for the selection of appropriate TFs and for the design of specific mini-promoters to improve cancer gene therapy.

Electronic supplementary material

The online version of this article (10.1186/s12885-018-4421-7) contains supplementary material, which is available to authorized users.

Identification of potential prognostic microRNA biomarkers for predicting survival in patients with hepatocellular carcinoma

Background

The aim of the present study was to identify potential prognostic microRNA (miRNA) biomarkers for hepatocellular carcinoma (HCC) prognosis prediction based on a dataset from The Cancer Genome Atlas (TCGA).

Materials and methods

A miRNA sequencing dataset and corresponding clinical parameters of HCC were obtained from TCGA. Genome-wide univariate Cox regression analysis was used to screen prognostic differentially expressed miRNAs (DEMs), and multivariable Cox regression analysis was used for prognostic signature construction. Comprehensive survival analysis was performed to evaluate the prognostic value of the prognostic signature.

Results

Five miRNAs were regarded as prognostic DEMs and used for prognostic signature construction. The five-DEM prognostic signature performed well in prognosis prediction (adjusted P < 0.0001, adjusted hazard ratio = 2.249, 95% confidence interval =1.491-3.394), and time-dependent receiver-operating characteristic (ROC) analysis showed an area under the curve (AUC) of 0.765, 0.745, 0.725, and 0.687 for 1-, 2-, 3-, and 5-year HCC overall survival (OS) prediction, respectively. Comprehensive survival analysis of the prognostic signature suggests that the risk score model could serve as an independent factor of HCC and perform better in prognosis prediction than other traditional clinical indicators. Functional assessment of the target genes of hsa-mir-139 and hsa-mir-5003 indicates that they were significantly enriched in multiple biological processes and pathways, including cell proliferation and cell migration regulation, pathways in cancer, and the cyclic adenosine monophosphate (cAMP) signaling pathway.

Conclusion

Our study indicates that the novel miRNA expression signature may be a potential prognostic biomarker for HCC patients.

Primary tumor microRNA signature predicts recurrence and survival in patients with locally advanced esophageal adenocarcinoma

Esophageal adenocarcinoma (EAC) is an aggressive cancer necessitating the development of improved risk stratification tools for personalized care. Previously, microRNAs have been shown to correlate with the progression and prognosis of various cancer types; however, the value in EAC remains largely unexplored. We performed global microRNA profiling on 32 formalin-fixed, paraffin-embedded EAC specimens to identify microRNAs associated with progression. Literature search and pathway analysis further refined output to five significantly deregulated candidate biomarkers. Four of the five microRNAs (miR-652-5p, miR-7-2-3p, miR-3925-3p, and miR-219-3p) were validated by qRT-PCR. Survival outcomes were evaluated in testing set of 26 stage II/III EAC patients to determine the prognostic relevance of the selected microRNAs. In the testing set, miR-652-5p and miR-7-2-3p expressions were significantly associated with progression-free survival (p-value = .00771 and p-value = .00293). The highest area under receiver operating characteristic (ROC) curve was 0.8212 for the combination of miR-652-5p and miR-7-2-3p. Collectively, our findings demonstrated that the miR-652-5p/miR-7-2-3p signature may serve as a promising prognostic marker in patients with locally advanced EAC.

Integrated Characterization of MicroRNA and mRNA Transcriptome in Papillary Thyroid Carcinoma

The incidence rate of papillary thyroid carcinoma (PTC) has rapidly increased in the recent decades, and the microRNA (miRNA) is one of the potential biomarkers in this cancer. Despite good prognosis, certain features such as lymph node metastasis (LNM) and BRAF V600E mutation are associated with a poor outcome. More than 50% of PTC patients present with LNM and BRAF V600E is the most common mutation identified in this cancer. The molecular mechanisms underlying these features are yet to be elucidated. This study aims to elucidate miRNA-genes interaction networks in PTC with or without LNM and to determine the association of BRAF V600E mutation with miRNAs and genes expression profiles. Next generation sequencing was performed to characterize miRNA and gene expression profiles in 20 fresh frozen tumor and the normal adjacent tissues of PTC with LNM positive (PTC LNM-P) and PTC without LNM (PTC LNN). BRAF V600E was genotyped using Sanger sequencing. Bioinformatics integration and pathway analysis were performed to determine the regulatory networks involved. Based on network analysis, we then investigated the association between miRNA and gene biomarkers, and pathway enrichment analysis was performed to study the role of candidate biomarkers. We identified 138 and 43 significantly deregulated miRNAs (adjusted p value < 0.05; log2 fold change ≤ -1.0 or ≥1.0) in PTC LNM-P and PTC LNN compared to adjacent normal tissues, respectively. Ninety-six miRNAs had significant expression ratios of 3p-to-5p in PTC LNM-P as compared to PTC LNN. In addition, ribosomal RNA-reduced RNA sequencing analysis revealed 699 significantly deregulated genes in PTC LNM-P versus normal adjacent tissues, 1,362 genes in PTC LNN versus normal adjacent tissue, and 1,576 genes in PTC LNM-P versus PTC LNN. We provide the evidence of miRNA and gene interactions, which are involved in LNM of papillary thyroid cancer. These findings may lead to better understanding of carcinogenesis and metastasis processes. This study also complements the existing knowledge about deregulated miRNAs in papillary thyroid carcinoma development.

Role of microRNAs in endocrine cancer metastasis

The deregulation of transcription and processing of microRNAs (miRNAs), as well as their function, has been involved in the pathogenesis of several human diseases, including cancer. Despite advances in therapeutic approaches, cancer still represents one of the major health problems worldwide. Cancer metastasis is an aggravating factor in tumor progression, related to increased treatment complexity and a worse prognosis. After more than one decade of extensive studies of miRNAs, the fundamental role of these molecules in cancer progression and metastasis is beginning to be elucidated. Recent evidences have demonstrated a significant role of miRNAs on the metastatic cascade, acting either as pro-metastatic or anti-metastatic. They are involved in distinct steps of metastasis including epithelial-to-mesenchymal transition, migration/invasion, anoikis survival, and distant organ colonization. Studies on the roles of miRNAs in cancer have focused mainly on two fronts: the establishment of a miRNA signature for different tumors, which may aid in early diagnosis using these miRNAs as markers, and functional studies of specific miRNAs, determining their targets, function and regulation. Functional miRNA studies on endocrine cancers are still scarce and represent an important area of research, since some tumors, although not frequent, present a high mortality rate. Among the endocrine tumors, thyroid cancer is the most common and best studied. Several miRNAs show lowered expression in endocrine cancers (i.e. miR-200s, miR-126, miR-7, miR-29a, miR-30a, miR-137, miR-206, miR-101, miR-613, miR-539, miR-205, miR-9, miR-195), while others are commonly overexpressed (i.e. miR-21, miR-183, miR-31, miR-let7b, miR-584, miR-146b, miR-221, miR-222, miR-25, miR-595). Additionally, some miRNAs were found in serum exosomes (miR-151, miR-145, miR-31), potentially serving as diagnostic tools. In this review, we summarize studies concerning the discovery and functions of miRNAs and their regulatory roles in endocrine cancer metastasis, which may contribute for the finding of novel therapeutic targets. The review focus on miRNAs with at least some identified targets, with established functions and, if possible, upstream regulation.

Omics Approaches to Identify Potential Biomarkers of Inflammatory Diseases in the Focal Adhesion Complex

Inflammatory diseases such as inflammatory bowel disease (IBD) require recurrent invasive tests, including blood tests, radiology, and endoscopic evaluation both to diagnose and assess disease activity, and to determine optimal therapeutic strategies. Simple ‘bedside’ biomarkers could be used in all phases of patient management to avoid unnecessary investigation and guide further management. The focal adhesion complex (FAC) has been implicated in the pathogenesis of multiple inflammatory diseases, including IBD, rheumatoid arthritis, and multiple sclerosis. Utilizing omics technologies has proven to be an efficient approach to identify biomarkers from within the FAC in the field of cancer medicine. Predictive biomarkers are paving the way for the success of precision medicine for cancer patients, but inflammatory diseases have lagged behind in this respect. This review explores the current status of biomarker prediction for inflammatory diseases from within the FAC using omics technologies and highlights the benefits of future potential biomarker identification approaches.