RESUMEN
Thyroid cancer (TC) is the most common endocrine cancer, accounting for 1.7% of all cancer cases. It has been reported that the existing approach to diagnosing TC is problematic. Therefore, it is essential to develop molecular biomarkers to improve the accuracy of the diagnosis. This study aimed to screen hub lncRNAs in the ceRNA network (ceRNET) connected to TC formation and progression based on the overall survival rate. In this study, first, RNA-seq data from the GDC database were collected. A package called edgeR in R programming language was then used to obtain differentially expressed lncRNAs (DElncRNAs), miRNAs (DEmiRNAs), and mRNAs (DEmRNAs) in TC patients' samples compared to normal samples. Second, DEmRNAs were analyzed for their functional enrichment. Third, to identify RNAs associated with overall survival, the overall survival of these RNAs was analyzed using the Kaplan-Meier plotter database to create a survival associated with the ceRNA network (survival-related ceRNET). Next, the GeneMANIA plugin was used to construct a PPI network to better understand survival-related DEmRNA interactions. The survival ceRNET was then visualized with the Cytoscape software, and hub genes, including hub lncRNAs and hub mRNAs, were identified using the CytoHubba plugin. We found 45 DElncRNAs, 28 DEmiRNAs, and 723 DEmRNAs among thyroid tumor tissue and non-tumor tissue samples. According to KEGG, GO and DO analyses, 723 DEmRNAs were mainly enriched in cancer-related pathways. Importantly, the results found that ten DElncRNAs, four DEmiRNAs, and 68 DEmRNAs are associated with overall survival. In this account, the PPI network was constructed for 68 survival-related DEmRNAs, and ADAMTS9, DTX4, and CLDN10 were identified as hub genes. The ceRNET was created by combining six lncRNAs, 109 miRNAs, and 22 mRNAs related to survival using Cytoscape. in this network, ten hub RNAs were identified by the CytoHubba plugin, including mRNAs (CTXND1, XKRX, IGFBP2, ENTPD1, GALNT7, ADAMTS9) and lncRNAs (AC090673.1, AL162511.1, LINC02454, AL365259.1). This study suggests that three lncRNAs, including AL162511.1, AC090673.1, and AL365259.1, could be reliable diagnostic biomarkers for TC. The findings of this study provide a basis for future studies on the therapeutic potential of these lncRNAs.
RESUMEN
Hepatocellular carcinoma (HCC) is one of the most prevalent cancers worldwide, which has a high mortality rate and poor treatment outcomes with yet unknown molecular basis. It seems that gene expression plays a pivotal role in the pathogenesis of the disease. Circular RNAs (circRNAs) can interact with microRNAs (miRNAs) to regulate gene expression in various malignancies by acting as competitive endogenous RNAs (ceRNAs). However, the potential pathogenesis roles of the ceRNA network among circRNA/miRNA/mRNA in HCC are unclear. In this study, first, the HCC circRNA expression data were obtained from three Gene Expression Omnibus microarray datasets (GSE164803, GSE94508, GSE97332), and the differentially expressed circRNAs (DECs) were identified using R limma package. Also, the liver hepatocellular carcinoma (LIHC) miRNA and mRNA sequence data were retrieved from TCGA and differentially expressed miRNAs (DEMIs) and mRNAs (DEGs) were determined using the R DESeq2 package. Second, CSCD website was used to uncover the binding sites of miRNAs on DECs. The DECs' potential target miRNAs were revealed by conducting an intersection between predicted miRNAs from CSCD and downregulated DEMIs. Third, candidate genes were uncovered by intersecting targeted genes predicted by miRWalk and targetscan online tools with upregulated DEGs. The ceRNA network was then built using the Cytoscape software. The functional enrichment and the overall survival time of these potential targeted genes were analyzed, and a PPI network was constructed in the STRING database. Network visualization was performed by Cytoscape, and ten hub genes were detected using the CytoHubba plugin tool. Four DECs (hsa_circ_0000520, hsa_circ_0008616, hsa_circ_0070934, hsa_circ_0004315) were obtained and six miRNAs (hsa-miR-542-5p, hsa-miR-326, hsa-miR-511-5p, hsa-miR-195-5p, hsa-miR-214-3p, and hsa-miR-424-5p) which are regulated by the above DECs were identified. Then 543 overlapped genes regulated by six miRNAs mentioned above were predicted. Functional enrichment analysis showed that these genes are mostly associated with regulatory pathways in cancer. Ten hub genes (TTK, AURKB, KIF20A, KIF23, CEP55, CDC6, DTL, NCAPG, CENPF, PLK4) have been screened from the PPI network of the 204 survival-related genes. KIF20A, NCAPG, TTK, PLK4, and CDC6 were selected for the highest significance p-values. At the end, a circRNA-miRNA-mRNA regulatory axis was established for five final selected hub genes. This study implies the potential pathogenesis of the obtained network and proposes that the two DECs (has_circ_0070934 and has_circ_0004315) may be important prognostic markers for HCC.
