Downregulation of miR-199a-3p in Hepatocellular Carcinoma and Its Relevant Molecular Mechanism via GEO, TCGA Database and In Silico Analyses.

Existing reports have demonstrated that miR-199a-3p plays a role as a tumor suppressor in a variety of human cancers. This study aims to further validate the expression of miR-199a-3p in HCC and to explore its underlying mechanisms by using multiple data sets. Chip data or sequencing data and quantitative reverse transcription polymerase chain reaction (qRT-PCR) were integrated to assess the expression of miR-199a-3p in HCC. The potential targets and transcription factor regulatory network of miR-199a-3p in HCC were determined and possible biological mechanism of miR-199a-3p was analyzed with bioinformatics methods. In the results, miR-199a-3p expression was significantly lower in HCC tissues compared to normal tissues according to chip data or sequencing data and qRT-PCR. Moreover, 455 targets of miR-199a-3p were confirmed, and these genes were involved in the PI3K-Akt signaling pathway, pathways in cancer, and focal adhesions. LAMA4 was considered a key target of miR-199a-3p. In CMTCN, 11 co-regulatory pairs, 3 TF-FFLs, and 2 composite-FFLs were constructed. In conclusion, miR-199a-3p was down regulated in HCC and LAMA4 may be a potential target of miR-199a-3p in HCC.

Upregulated expression of SAC3D1 is associated with progression in gastric cancer.

SAC3 domain containing 1 (SAC3D1) has been reported to be involved in numerous types of cancer. However, the role of SAC3D1 in GC has not yet been elucidated. In the present study, the mRNA expression level of SAC3D1 between GC and normal tissues were assessed with a continuous variable meta­analysis based on multiple datasets from public databases. The protein expression level of SAC3D1 in GC and normal tissues was assessed by an in­house immunohistochemistry (IHC). The association between SAC3D1 expression and some clinical parameters was assessed based on the TCGA and IHC data. Survival analysis was performed to assess the association between SAC3D1 expression and the survival of GC patients. The co­expressed genes of SAC3D1 were determined by integrating three online tools, and the enrichment analyses were performed to determine SAC3D1­related pathways and hub co­expressed genes. SAC3D1 was significantly upregulated in GC tumor tissues in comparison to normal tissues with the SMD being 0.45 (0.12, 0.79). The IHC results also indicated that SAC3D1 protein expression in GC tissues was markedly higher than in normal tissues. The SMD following the addition of the IHC data was 0.59 (0.11, 1.07). The protein levels of SAC3D1 were positively associated with the histological grade, T stage and N stage of GC (P<0.001). The TCGA data also revealed that the SAC3D1 mRNA level was significantly associated with the N stage (P<0.001). Moreover, prognosis analysis indicated that SAC3D1 was closely associated with the prognosis of patients with GC. Moreover, 410 co­expressed genes of SAC3D1 were determined, and these genes were mainly enriched in the cell cycle. In total, 4 genes (CDK1, CCNB1, CCNB2 and CDC20) were considered key co­expressed genes. On the whole, these findings demonstrate that SAC3D1 is highly expressed in GC and may be associated with the progression of GC.

Upregulation of miR­132­3p in cholangiocarcinoma tissues: A study based on RT­qPCR, The Cancer Genome Atlas miRNA sequencing, Gene Expression Omnibus microarray data and bioinformatics analyses.

MicroRNAs (miRNAs/miRs) have been reported to be closely associated with numerous human diseases, including cholangiocarcinoma (CCA). However, the number of miRNAs known to be involved in CCA is limited, and the association between miR­132­3p and CCA remains unknown. In the present study, the clinical role of miR­132­3p and its potential signaling pathways were investigated by multiple approaches. Reverse transcription­quantitative PCR (RT­qPCR), CCA­associated Gene Expression Omnibus (GEO), ArrayExpress and Sequence Read Archive (SRA) miRNA­microarray or miRNA­sequencing data were screened, and meta­analyses were conducted, in order to calculate the receiver operating characteristic (ROC) curve and standardized mean difference (SMD). The predicted target genes of miR­132­3p were obtained from 12 online databases and were combined with the downregulated differentially expressed genes identified in the RNA­sequencing data of CCA. Gene Ontology annotation and pathway analysis were performed in WebGestalt. Protein­protein interaction analyses were conducted in STRING. The Cancer Genome Atlas (TCGA) mRNA expression profiles were used to validate the expression levels of hub genes at the mRNA level. The Human Protein Atlas was used to identify the protein expression levels of hub genes in CCA tissues and non­tumor biliary epithelium. The meta­analyses comprised 10 groups of RT­qPCR data, eight GEO microarray datasets and one TCGA miRNA­sequencing dataset. The SMD of miR­132­3p in CCA was 0.75 (95% CI: 0.25, 1.24), which indicated that miR­132­3p was overexpressed in CCA tissues. This finding was supported by a summary ROC value of 0.80 (95% CI: 0.76, 0.83). The pooled sensitivity and specificity were 0.81 (95% CI: 0.59, 0.93) and 0.71 (95% CI: 0.58, 0.81), respectively. The relative expression level of miR­132­3p in the early stage of CCA (stages I­II) was 6.8754±0.5279, which was markedly lower than that in the advanced stage (stages III­IVB), 7.3034±0.3267 (P=0.003). Consistently, the miR­132­3p level in low­grade CCA (grades G1­G2) was 6.7581±0.5297, whereas it was 7.1191±0.4651 in patients with high­grade CCA (grades G3­G4) (P=0.037). Furthermore, 555 potential target genes of miR­132­3p in CCA were mainly enriched in the 'Focal Adhesion­PI3K­Akt­mTOR­signaling pathway'. In conclusion, upregulation of miR­132­3p may serve a pivotal role in the tumorigenesis and progression of CCA by targeting different pathways. Further in vitro and in vivo studies are required to support the current findings.

Novel drug candidate for the treatment of several soft­tissue sarcoma histologic subtypes: A computational method using survival­associated gene signatures for drug repurposing.

Systemic treatment options for soft tissue sarcomas (STSs) have remained unchanged despite the need for novel drug candidates to improve STS outcomes. Drug repurposing involves the application of clinical drugs to different diseases, reducing development time, and cost. It has also become a fast and effective way to identify drug candidates. The present study used a computational method to screen three drug­gene interaction databases for novel drug candidates for the treatment of several common STS histologic subtypes through drug repurposing. STS survival­associated genes were generated by conducting a univariate cox regression analysis using The Cancer Genome Atlas survival data. These genes were then applied to three databases (the Connectivity Map, the Drug Gene Interaction Database and the L1000 Fireworks Display) to identify drug candidates for STS treatment. Additionally, pathway analysis and molecular docking were conducted to evaluate the molecular mechanisms of the candidate drug. Bepridil was identified as a potential candidate for several STS histologic subtype treatments by overlapping the screening results from three drug­gene interaction databases. The pathway analysis with the Kyoto Encyclopedia of Genes and Genomes predicted that Bepridil may target CRK, fibroblast growth factor receptor 4 (FGFR4), laminin subunit ß1 (LAMB1), phosphoinositide­3­kinase regulatory subunit 2 (PIK3R2), WNT5A, cluster of differentiation 47 (CD47), elastase, neutrophil expressed (ELANE), 15­hydroxyprostaglandin dehydrogenase (HPGD) and protein kinase cß (PRKCB) to suppress STS development. Further molecular docking simulation suggested a relatively stable binding selectivity between Bepridil and eight proteins (CRK, FGFR4, LAMB1, PIK3R2, CD47, ELANE, HPGD, and PRKCB). In conclusion, a computational method was used to identify Bepridil as a potential candidate for the treatment of several common STS histologic subtypes. Experimental validation of these in silico results is necessary before clinical translation can occur.

Exploration and validation of downregulated microRNA-199a-3p, downstream messenger RNA targets and transcriptional regulation in osteosarcoma.

Osteosarcoma (OS) is a primary bone tumor with a high incidence and mortality in children and adolescents. Emerging evidence shows that microRNAs (miRNAs) participate in biological tumor mechanisms by targeting downstream messenger RNAs (mRNAs). This article aimed to investigate the potential regulatory targets of microRNA-199a-3p (miR-199a-3p) in OS and to contribute to the understanding of miR-199a-3p-related OS regulatory mechanisms. MicroRNA-related Gene Expression Omnibus (GEO) chips, ArrayExpress chips and literature data were used to determine the expression of miR-199a-3p in OS and pooled to explore its potential clinical value. To investigate the target genes of miR-199a-3p further, we integrated the results from the following three-part gene study: Twelve online prediction tools were used to predict the target genes of miR-199a-3p; the GEO GSE89370 chip transfected with miRSelect pEP-miR-199a-3p was used to analyze the downregulated differentially expressed genes (DEGs) in OS cells; and highly expressed DEGs were derived from an in-house microarray generated from three pairs of clinical OS and normal tissue samples acquired through our department. Then, we analyzed the target genes using the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases and the protein-protein interaction (PPI) network to further identify the primary target genes. In addition, we constructed transcription factor (TF)-miRNA-joint gene feed-forward regulatory loops (FFLs) with Circuits DB using miR-199a-3p as the core. A comprehensive meta-analysis of a hub of miR-199a-3p targeted genes was performed to integrate expression level, summary ROC (sROC) curves and survival analysis results from the GEO data for verification and exploration. Finally, the expression levels of the hub genes were verified in OS tissues and U2OS cells by immunohistochemistry (IHC) and immunocytochemistry (ICC). Data on miR-199a-3p expression were obtained from three data sets (GSE65071, GSE69524, and PMID 21666078), which showed low miR-199a-3p expression levels in OS tissues. The combined data indicated the same tendency, with the SMD of the random effect model, as shown in forest plots, being -2.8 (95% CI: -4.49, -1.11). In addition, we determined that miR-199a-3p may serve as a molecular marker useful for distinguishing OS tissues from normal tissues with high sensitivity and specificity, with the measured outcomes being 0.94 (95% CI: 0.80, 0.99) and 0.96 (95% CI: 0.78, 1.00), respectively. In addition, 391 genes were considered targets of miR-199a-3p in OS, and the enrichment analysis indicated that these targets were mainly enriched in proteoglycans in cancer and in spliceosomes. Four genes, CDKI, CCNB1, AURKA and NEK2, were regarded as hub targets based on the PPI data. Subsequently, TF-miRNA-joint genes FFLs were constructed in Circuits DB and included 17 TFs and 82 joint targets. These joint targets were mainly enriched in spliceosomes. UBE2D1 and RBM25 were regarded as hub joint targets based on the enrichment analysis. All selected target genes were further verified to ensure that they were upregulated in OS and to determine their prognostic significance. At the experimental verification level, the CDK1 protein was confirmed to be positively expressed in the cytoplasm of OS tissues and the U2OS cell line. Our study verified that miR-199a-3p was obviously downregulated in OS. CDK1, CCNB1, NEK2, AURKA, UBE2D1 and RBM25 were identified as potential target genes of miR-199a-3p in OS.

Expression and potential molecular mechanisms of miR­204­5p in breast cancer, based on bioinformatics and a meta­analysis of 2,306 cases.

Breast cancer (BC) is the most common cancer among women worldwide. However, there is insufficient research that focuses on the expression and molecular mechanisms of microRNA (miR)­204­5p in BC. In the current study, data were downloaded from the Cancer Genome Atlas (TCGA), the Gene Expression Omnibus (GEO) and the University of California Santa Cruz (UCSC) Xena databases. They were then used to undertake a meta­analysis that leveraged the standard mean difference (SMD) and summarized receiver operating characteristic (sROC) to evaluate the expression of the precursor miR­204 and mature miR­204­5p in BC. Additionally, an intersection of predicted genes, differentially expressed genes (DEGs) from the TCGA database and the GEO database were plotted to acquire desirable putative genes. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and protein­protein interaction (PPI) network analyses were performed to assess the potential pathways and hub genes of miR­204­5p in BC. A decreased trend in precursor miR­204 expression was detected in 1,077 BC tissue samples in comparison to 104 para­carcinoma tissue samples in the TCGA database. Further, the expression of mature miR­204­5p was markedly downregulated in 756 BC tissue samples in comparison to 76 para­carcinoma tissue samples in the UCSC Xena database. The outcome of the SMD from meta­analysis also indicated that the expression of miR­204­5p was markedly reduced in 2,306 BC tissue samples in comparison to 367 para­carcinoma tissue samples. Additionally, the ROC and sROC values indicated that miR­204­5p had a great discriminatory capacity for BC. In GO analysis, 'cell development', 'cell surface activity', and 'receptor agonist activity' were the most enriched terms; in KEGG analysis, 'endocytosis' was significantly enriched. Rac GTPase activating protein 1 (RACGAP1) was considered the hub gene in the PPI network. In conclusion, miR­204­5p may serve a suppressor role in the oncogenesis and advancement of BC, and miR­204­5p may have crucial functions in BC by targeting RACGAP1.

Clinical value of miR-198-5p in lung squamous cell carcinoma assessed using microarray and RT-qPCR.

BACKGROUND: To examine the clinical value of miR-198-5p in lung squamous cell carcinoma (LUSC). METHODS: Gene Expression Omnibus (GEO) microarray datasets were used to explore the miR-198-5p expression and its diagnostic value in LUSC. Real-time reverse transcription quantitative polymerase chain reaction was used to evaluate the expression of miR-198-5p in 23 formalin-fixed, paraffin-embedded (FFPE) LUSC tissues and corresponding non-cancerous tissues. The correlation between miR-198-5p expression and clinic pathological features was assessed. Meanwhile, putative target messenger RNAs of miR-198-5p were identified based on the analysis of differentially expressed genes in the Cancer Genome Atlas (TCGA) and 12 miRNA prediction tools. Subsequently, the putative target genes were sent to Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses. RESULTS: MiR-198-5p was low expressed in LUSC tissues. The combined standard mean difference (SMD) values of miR-198-5p expression based on GEO datasets were - 0.30 (95% confidence interval (CI) - 0.54, - 0.06) and - 0.39 (95% CI - 0.83, 0.05) using fixed effect model and random effect model, respectively. The sensitivity and specificity were not sufficiently high, as the area under the curve (AUC) was 0.7749 (Q* = 0.7143) based on summarized receiver operating characteristic (SROC) curves constructed using GEO datasets. Based on the in-house RT-qPCR, miR-198-5p expression was 4.3826 ± 1.7660 in LUSC tissues and 4.4522 ± 1.8263 in adjacent normal tissues (P = 0.885). The expression of miR-198-5p was significantly higher in patients with early TNM stages (I-II) than that in cases with advanced TNM stages (III-IV) (5.4400 ± 1.5277 vs 3.5690 ± 1.5228, P = 0.008). Continuous variable-based meta-analysis of GEO and PCR data displayed the SMD values of - 0.26 (95% CI - 0.48, - 0.04) and - 0.34 (95% CI - 0.71, 0.04) based on fixed and random effect models, respectively. As for the diagnostic value of miR-198-5p, the AUC based on the SROC curve using GEO and PCR data was 0.7351 (Q* = 0.6812). In total, 542 genes were identified as the targets of miR-198-5p. The most enriched Gene Ontology terms were epidermis development among biological processes, cell junction among cellular components, and protein dimerization activity among molecule functions. The pathway of non-small cell lung cancer was the most significant pathway identified using Kyoto Encyclopedia of Genes and Genomes analysis. CONCLUSION: The expression of miR-198-5p is related to the TNM stage. Thus, miR-198-5p might play an important role via its target genes in LUSC.