Co-localized genomic regulation of miRNA and mRNA via DNA methylation affects survival in multiple tumor types

Pan-Cancer Analysis Reveals SH3TC2 as an Oncogene for Colorectal Cancer and Promotes Tumorigenesis via the MAPK Pathway

Simple Summary

SH3 domain and tetrapeptide repeat 2 (SH3TC2) is a protein-encoding gene and has previously been described as a critical signaling hub for neurological disorders, but no systematic analysis of SH3TC2 is available in cancer research. We analyzed SH3TC2 in various kinds of cancer to find its tumorigenic role in one or more specific cancers and further explored the mechanism of SH3TC2 in colorectal cancer (CRC). Our research revealed that higher expression of SH3TC2 indicated poor disease-free survival and promoted CRC progression and invasion via the MAPK signaling pathway.

Abstract

SH3 domain and tetrapeptide repeat 2 (SH3TC2) is a protein-encoding gene and has previously been described as a critical signaling hub for neurological disorders. Although increasing evidence supports a vital role of SH3TC2 in the tumorigenesis of various kinds of cancer, no systematic analysis of SH3TC2 is available. The function and mechanism of SH3TC2 in other cancers remain unknown. Thus, this study aimed to analyze SH3TC2 in various kinds of cancer to find its tumorigenic role in one or more specific cancers. In the current study, we analyzed the expression level and prognostic value of SH3TC2 in different tumors in the TCGA-GTEx pan-cancer dataset. Subsequently, the prognostic role and mechanism of SH3TC2 in colorectal cancer (CRC) were further explored via clinical samples and in vitro and in vivo experiments. We observed differential expression of SH3TC2 in colon adenocarcinoma (COAD), acute myeloid leukemia (LAML), READ (rectum adenocarcinoma), SKCM (skin cutaneous melanoma), and TGCT (testicular germ cell tumors). Subsequently, SH3TC2 showed a significant effect on the clinical stage and prognostic value in CRC, LAML, and SKCM. Moreover, we found in the TCGA database and seven GEO datasets that SH3TC2 was significantly highly expressed in tumor tissue. Through enrichment analysis of SH3TC2 and its co-expressed genes, we found that SH3TC2 may play a role in the MAPK signaling pathway. Correlation analysis indicated that SH3TC2 was significantly associated with multiple key factors in the MAPK signaling pathway. Additionally, higher expression of SH3TC2 was found in tumor tissue in our cohort including 40 CRC patients. Overexpression of SH3TC2 may imply poor prognosis. Knockdown of SH3TC2 significantly inhibited tumor invasion, migration, and proliferation. More importantly, knockdown of SH3TC2 inhibited tumor growth in a CRC mouse model. The study preliminarily conducted a pan-cancer study of SH3TC2 and further explored the mechanism of SH3TC2 in CRC. Our research revealed that higher expression of SH3TC2 may promote CRC progression and invasion via the MAPK signaling pathway.

MicroRNA‑218 inhibits the malignant phenotypes of glioma by modulating the TNC/AKT/AP‑1/TGFβ1 feedback signaling loop

Gliomas are the most malignant and common tumors of the human brain, and the prognosis of glioma patients is extremely poor MicroRNAs (miRNAs or miRs) play critical roles in different types of cancer by performing post-transcriptional regulation of gene expression Although miR-218 has been demonstrated to be decreased in gliomas, its role in gliomas remains largely unknown miR-218 expression was analyzed in gliomas and normal brain tissues (control subjects) using a dataset from The Cancer Genome Atlas A series of in vitro and in vivo studies were performed to determine the biological roles of miR-218 in glioma cells Potential targets of miR-218 were identified using a dual-luciferase reporter system Western blot and dual-luciferase reporter system experiments were performed to evaluate the regulatory effect of miR-218 on the tenascin C (TNC)/AKT/activator protein 1 (AP-1)/transforming growth factor β1 (TGFβ1) pathway It was demonstrated that miR-218 was significantly downregulated in gliomas compared with control subjects, and played potent tumor suppressor roles in glioma cells by inhibiting cell proliferation, colony formation, migration, invasion and tumorigenic potential in nude mice, as well as inducing cell cycle arrest and apoptosis Mechanistically, miR-218 inhibited malignant phenotypes of glioma cells by binding to the 3′-untranslated region of its target TNC and subsequently suppressing its expression As a result, miR-218 could reduce AKT phosphorylation and subsequently inhibit transcriptional activity of AP-1 by reducing JNK phosphorylation, downregulating the expression of TGFβ1, while TGFβ1 was able to, in turn, activate the TNC/AKT/AP-1 signaling axis Our data revealed a previously unknown tumor suppressor role of miR-218 by blocking the TNC/AKT/AP-1/TGFβ1-positive feedback loop in glioma

Epigenetic Regulation of Breast Cancer Stem Cells Contributing to Carcinogenesis and Therapeutic Implications

Globally, breast cancer has remained the most commonly diagnosed cancer and the leading cause of cancer death among women. Breast cancer is a highly heterogeneous and phenotypically diverse group of diseases, which require different selection of treatments. Breast cancer stem cells (BCSCs), a small subset of cancer cells with stem cell-like properties, play essential roles in breast cancer progression, recurrence, metastasis, chemoresistance and treatments. Epigenetics is defined as inheritable changes in gene expression without alteration in DNA sequence. Epigenetic regulation includes DNA methylation and demethylation, as well as histone modifications. Aberrant epigenetic regulation results in carcinogenesis. In this review, the mechanism of epigenetic regulation involved in carcinogenesis, therapeutic resistance and metastasis of BCSCs will be discussed, and finally, the therapies targeting these biomarkers will be presented.

Methylation in HOX Clusters and Its Applications in Cancer Therapy

HOX genes are commonly known for their role in embryonic development, defining the positional identity of most structures along the anterior–posterior axis. In postembryonic life, HOX gene aberrant expression can affect several processes involved in tumorigenesis such as proliferation, apoptosis, migration and invasion. Epigenetic modifications are implicated in gene expression deregulation, and it is accepted that methylation events affecting HOX gene expression play crucial roles in tumorigenesis. In fact, specific methylation profiles in the HOX gene sequence or in HOX-associated histones are recognized as potential biomarkers in several cancers, helping in the prediction of disease outcomes and adding information for decisions regarding the patient’s treatment. The methylation of some HOX genes can be associated with chemotherapy resistance, and its identification may suggest the use of other treatment options. The use of epigenetic drugs affecting generalized or specific DNA methylation profiles, an approach that now deserves much attention, seems likely to be a promising weapon in cancer therapy in the near future. In this review, we summarize these topics, focusing particularly on how the regulation of epigenetic processes may be used in cancer therapy.

High Expression of the SH3TC2-DT/SH3TC2 Gene Pair Associated With FLT3 Mutation and Poor Survival in Acute Myeloid Leukemia: An Integrated TCGA Analysis

Fms-like tyrosine kinase 3 (FLT3) mutation is one of the most common mutations in acute myeloid leukemia (AML). However, the effect of FLT3 mutation on survival is currently still controversial and the leukemogenic mechanisms are still under further investigation. The aim of our study is to identify differentially expressed genes (DEGs) in FLT3-mutant AML and to find crucial DEGs whose expression level is related to prognosis for further analysis. By mining the TCGA-LAML dataset, 619 differentially expressed lncRNAs (DElncRNAs) and 1,428 differentially expressed mRNAs (DEmRNAs) were identified between FLT3-mutant and FLT3-wildtype samples. Through weighted gene correlation network analysis (WGCNA) and the following Cox proportional hazards regression analysis, we constructed the prognostic risk models to identify the hub DElncRNAs and DEmRNAs associated with AML prognosis. The presence of both SH3TC2 divergent transcript (SH3TC2-DT) and SH3TC2 in respective prognostic risk models promotes us to further study the significance of this gene pair in AML. SH3TC2-DT and SH3TC2 were identified to be coordinately high expressed in FLT3-mutant AML samples. High expression of this gene pair was associated with poor survival. Using logistic regression analysis, we found that high SH3TC2-DT/SH3TC2 expression was associated with FLT3 mutation, high WBC count, and intermediate cytogenetic and molecular–genetic risk. AML with SH3TC2-DT/SH3TC2 high expression showed enrichment of transcripts associated with stemness, quiescence, and leukemogenesis. Our study suggests that the SH3TC2-DT/SH3TC2 gene pair may be a possible biomarker to further optimize AML prognosis and may function in stemness or quiescence of FLT3-mutant leukemic stem cells (LSCs).

The New Biomarker for Cervical Squamous Cell Carcinoma and Endocervical Adenocarcinoma (CESC) Based on Public Database Mining

To reconstruct the ceRNA biological network of cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC) and to select an appropriate mRNA as a biomarker that could be used for CESC early diagnosis and prognosis evaluation. We downloaded CESC data from the TCGA public database, and statistical analysis was conducted with the R software to find out differential expressed genes encoding for lncRNAs, miRNAs, and mRNAs. The differentially expressed mRNAs (DEmRNAs) screened in the ceRNA network were analyzed for survival to find the mRNAs with significantly linked to the survival prognosis. These mRNAs were searched in the Pathological Atlas to identify the final appropriate mRNAs. Differential expression analysis revealed 773 lncRNAs, 94 miRNAs, and 2466 mRNAs. Survival analysis of DEmRNAs in the ceRNA network indicated that ADGRF4, ANXA8L1, HCAR3, IRF6, and PDE2A (P < 0.05) were negatively correlated with survival time. Verification of these six DEmRNAs in the Pathology Atlas indicated that PDE2A was a possible biomarker for CESC patients. PDE2A might be a biomarker for early diagnosis and prognosis evaluation of CESC patients, but due to the lack of available data, further studies may be needed for confirmation.

Methylxanthines: Potential Therapeutic Agents for Glioblastoma

Glioblastoma (GBM) is the most common and aggressive primary brain tumor. Currently, treatment is ineffective and the median overall survival is 20.9 months. The poor prognosis of GBM is a consequence of several altered signaling pathways that favor the proliferation and survival of neoplastic cells. One of these pathways is the deregulation of phosphodiesterases (PDEs). These enzymes participate in the development of GBM and may have value as therapeutic targets to treat GBM. Methylxanthines (MXTs) such as caffeine, theophylline, and theobromine are PDE inhibitors and constitute a promising therapeutic anti-cancer agent against GBM. MTXs also regulate various cell processes such as proliferation, migration, cell death, and differentiation; these processes are related to cancer progression, making MXTs potential therapeutic agents in GBM.

Identification of key genes, pathways, and miRNA/mRNA regulatory networks of CUMS-induced depression in nucleus accumbens by integrated bioinformatics analysis

INTRODUCTION

Major depressive disorder (MDD) is a recurrent, devastating mental disorder, which affects >350 million people worldwide, and exerts substantial public health and financial costs to society. Thus, there is a significant need to discover innovative therapeutics to treat depression efficiently. Stress-induced dysfunction in the subtype of neuronal cells and the change of synaptic plasticity and structural plasticity of nucleus accumbens (NAc) are implicated in depression symptomology. However, the molecular and epigenetic mechanisms and stresses to the NAc pathological changes in depression remain elusive.

MATERIALS AND METHODS

In this study, treatment group mice were treated continually with the chronic unpredictable mild stress (CUMS) until expression of depression-like behaviors were found. Depression was confirmed with sucrose preference, novelty-suppressed feeding, forced swimming, and tail suspension tests. We applied high-throughput RNA sequencing to assess microRNA expression and transcriptional profiles in the NAc tissue from depression-like behaviors mice and control mice. The regulatory network of miRNAs/mRNAs was constructed based on the high-throughput RNA sequence and bioinformatics software predictions.

RESULTS

A total of 17 miRNAs and 10 mRNAs were significantly upregulated in the NAc of CUMS-induced mice with depression-like behaviors, and 12 miRNAs and 29 mRNAs were downregulated. A series of bioinformatics analyses showed that these altered miRNAs predicted target mRNA and differentially expressed mRNAs were significantly enriched in the MAPK signaling pathway, GABAergic synapse, dopaminergic synapse, cytokine-cytokine receptor interaction, axon guidance, regulation of autophagy, and so on. Furthermore, dual luciferase report assay and qRT-PCR results validated the miRNA/mRNA regulatory network.

CONCLUSION

The deteriorations of GABAergic synapses, dopaminergic synapses, neurotransmitter synthesis, as well as autophagy-associated apoptotic pathway are associated with the molecular pathological mechanism of CUMS-induced depression.

knnAUC: an open-source R package for detecting nonlinear dependence between one continuous variable and one binary variable

BACKGROUND

Testing the dependence of two variables is one of the fundamental tasks in statistics. In this work, we developed an open-source R package (knnAUC) for detecting nonlinear dependence between one continuous variable X and one binary dependent variables Y (0 or 1).

RESULTS

We addressed this problem by using knnAUC (k-nearest neighbors AUC test, the R package is available at https://sourceforge.net/projects/knnauc/ ). In the knnAUC software framework, we first resampled a dataset to get the training and testing dataset according to the sample ratio (from 0 to 1), and then constructed a k-nearest neighbors algorithm classifier to get the yhat estimator (the probability of y = 1) of testy (the true label of testing dataset). Finally, we calculated the AUC (area under the curve of receiver operating characteristic) estimator and tested whether the AUC estimator is greater than 0.5. To evaluate the advantages of knnAUC compared to seven other popular methods, we performed extensive simulations to explore the relationships between eight different methods and compared the false positive rates and statistical power using both simulated and real datasets (Chronic hepatitis B datasets and kidney cancer RNA-seq datasets).

CONCLUSIONS

We concluded that knnAUC is an efficient R package to test non-linear dependence between one continuous variable and one binary dependent variable especially in computational biology area.

Identification of key microRNAs and their targets in exosomes of pancreatic cancer using bioinformatics analysis

Pancreatic cancer (PC) is one of the most lethal tumors, due to late diagnosis and limited surgical strategies. It has been reported that serum exosomal microRNAs (S-Exo-miRNAs) play a pivotal role as signaling molecules and serve as noninvasive diagnosis methods for PC. The combination of S-Exo-miRNAs with the corresponding target also plays an important role in the tumor microenvironment.Here we investigated S-Exo-miRNAs involved in PC. The gene expression profile was downloaded from the Gene Expression Omnibus (GEO) database. The analysis was carried out using GEO2R. The targets of differentially expressed serum exosomal miRNAs (DE-S-Exo-miRNAs) were predicted by 4 bioinformatic algorithms (miRanda, miRDB, miRWalk, and Targetscan). Further analysis with gene ontology (GO) and Kyoto Encyclopedia of Genomes pathway (KEGG) enrichment analyses were performed with Cytoscape software version 3.4.0. Subsequently, the interaction regulatory network of target genes was performed with the Search Tool for the Retrieval of Interacting Genes (STRING) database (http://www.string-db.org/) and visualized using Cytoscape software.We downloaded the gene expression profile GSE50632, which was based on an Agilent microarray GPL17660 platform containing 4 eligible samples. In total 467 DE-S-Exo-miRNAs were obtained, including 7 overexpressed miRNAs (1.50%), and 460 remaining underexpressed miRNAs (98.50%). The databases miRWalk, miRDB, miRanda, and TargetScan were used to predict their potential targets, which were subsequently submitted to Cytoscape software version 3.4.0 (www.cytoscape.org). Next the functional and pathway enrichment analysis were used for the KEGG pathway and GO categories analysis. The enrichment analysis identified the genes involved in such processes as developmental and negative regulation of multicellular organismal processes, regulation of anatomical structure morphogenesis, regulation of cell death, apoptotic processes and mitogen-activated protein kinase (MAPK) signaling pathway, transforming growth factor - beta (TGF -β) signaling pathway, cyclic adenosine monophosphate (cAMP) signaling pathway, and the phosphatidylinositol-3 kinases/Akt (PI3K-Akt) signaling pathway. Subsequently according to the protein-protein interaction (PPI) network, the top 10 genes were obtained. The enrichment analyses of the genes involved in a significant module revealed that these genes were related to the TGF-β signaling pathway. After reviewing the literature, we identified the apoptosis genes, and their corresponding miRNAs that have a relationship with apoptosis of the tumor.This analysis provides a comprehensive understanding of the roles of S-Exo-miRNAs and the related targets in the development of PC. Additionally, the present study provides promising candidate targets for early diagnosis and therapeutic intervention. However, these predictions require further experimental validation in future studies.

Predicting miRNA targets for head and neck squamous cell carcinoma using an ensemble method

BACKGROUND

This study aimed to uncover potential microRNA (miRNA) targets in head and neck squamous cell carcinoma (HNSCC) using an ensemble method which combined 3 different methods: Pearson's correlation coefficient (PCC), Lasso and a causal inference method (i.e., intervention calculus when the directed acyclic graph (DAG) is absent [IDA]), based on Borda count election.

METHODS

The Borda count election method was used to integrate the top 100 predicted targets of each miRNA generated by individual methods. Afterwards, to validate the performance ability of our method, we checked the TarBase v6.0, miRecords v2013, miRWalk v2.0 and miRTarBase v4.5 databases to validate predictions for miRNAs. Pathway enrichment analysis of target genes in the top 1,000 miRNA-messenger RNA (mRNA) interactions was conducted to focus on significant KEGG pathways. Finally, we extracted target genes based on occurrence frequency ≥3.

RESULTS

Based on an absolute value of PCC >0.7, we found 33 miRNAs and 288 mRNAs for further analysis. We extracted 10 target genes with predicted frequencies not less than 3. The target gene MYO5C possessed the highest frequency, which was predicted by 7 different miRNAs. Significantly, a total of 8 pathways were identified; the pathways of cytokine-cytokine receptor interaction and chemokine signaling pathway were the most significant.

CONCLUSIONS

We successfully predicted target genes and pathways for HNSCC relying on miRNA expression data, mRNA expression profile, an ensemble method and pathway information. Our results may offer new information for the diagnosis and estimation of the prognosis of HNSCC.

Genomic, Pathway Network, and Immunologic Features Distinguishing Squamous Carcinomas

This integrated, multiplatform PanCancer Atlas study co-mapped and identified distinguishing molecular features of squamous cell carcinomas (SCCs) from five sites associated with smoking and/or human papillomavirus (HPV). SCCs harbor 3q, 5p, and other recurrent chromosomal copy-number alterations (CNAs), DNA mutations, and/or aberrant methylation of genes and microRNAs, which are correlated with the expression of multi-gene programs linked to squamous cell stemness, epithelial-to-mesenchymal differentiation, growth, genomic integrity, oxidative damage, death, and inflammation. Low-CNA SCCs tended to be HPV(+) and display hypermethylation with repression of TET1 demethylase and FANCF, previously linked to predisposition to SCC, or harbor mutations affecting CASP8, RAS-MAPK pathways, chromatin modifiers, and immunoregulatory molecules. We uncovered hypomethylation of the alternative promoter that drives expression of the ΔNp63 oncogene and embedded miR944. Co-expression of immune checkpoint, T-regulatory, and Myeloid suppressor cells signatures may explain reduced efficacy of immune therapy. These findings support possibilities for molecular classification and therapeutic approaches.

MicroRNA in Glioblastoma: An Overview

Glioblastoma is the most aggressive brain tumor and, even with the current multimodal therapy, is an invariably lethal cancer with a life expectancy that depends on the tumor subtype but, even in the most favorable cases, rarely exceeds 2 years. Epigenetic factors play an important role in gliomagenesis, are strong predictors of outcome, and are important determinants for the resistance to radio- and chemotherapy. The latest addition to the epigenetic machinery is the noncoding RNA (ncRNA), that is, RNA molecules that are not translated into a protein and that exert their function by base pairing with other nucleic acids in a reversible and nonmutational mode. MicroRNAs (miRNA) are a class of ncRNA of about 22 bp that regulate gene expression by binding to complementary sequences in the mRNA and silence its translation into proteins. MicroRNAs reversibly regulate transcription through nonmutational mechanisms; accordingly, they can be considered as epigenetic effectors. In this review, we will discuss the role of miRNA in glioma focusing on their role in drug resistance and on their potential applications in the therapy of this tumor.

MiR-93-5p promotes gastric cancer-cell progression via inactivation of the Hippo signaling pathway

MiR-93-5p has been previously found to be associated with gastric cancer (GC) tumorigenesis; however, the current understanding of its function in this context remains largely incomplete. In the present study, we showed that miR-93-5p was upregulated in GC tissues. We also demonstrated that miR-93-5p overexpression promoted the proliferation, migration, invasion, and chemoresistance of SGC-7901 cells in vitro, and conversely, that endogenously silencing miR-93-5p expression induced the opposite effects in HGC-27 cells. Overexpression of miR-93-5p was found to inactivate the Hippo pathway, and furthermore, miR-93-5p knockdown activated Hippo signaling. MiR-93-5p upregulation was also shown to inhibit the expression of two well-characterized Hippo pathway regulators, protocadherin Fat 4 (FAT4), and large tumor suppressors 2 (LATS2), at both the mRNA and protein level. Additionally, the results of bioinformatics analyses and luciferase reporter assays indicated that miR-93-5p directly targets the 3'-UTR of FAT4 and LATS2. Taken together, these results demonstrate that miR-93-5p promotes GC-cell progression via the inactivation of the Hippo signaling pathway, and thus, represents a potential therapeutic target for the treatment of GC.