scholarly journals Comprehensive analysis of lncRNA-associated competing endogenous RNA network in tongue squamous cell carcinoma

PeerJ ◽  
2019 ◽  
Vol 7 ◽  
pp. e6397 ◽  
Author(s):  
Shusen Zhang ◽  
Ruoyan Cao ◽  
Qiulan Li ◽  
Mianfeng Yao ◽  
Yu Chen ◽  
...  

Background Increasing evidence has demonstrated that long non-coding RNAs (lncRNAs) play an important role in the competitive endogenous RNA (ceRNA) networks in that they regulate protein-coding gene expression by sponging microRNAs (miRNAs). However, the understanding of the ceRNA network in tongue squamous cell carcinoma (TSCC) remains limited. Methods Expression profile data regarding mRNAs, miRNAs and lncRNAs as well as clinical information on 122 TSCC tissues and 15 normal controls from The Cancer Genome Atlas (TCGA) database were collected. We used the edgR package to identify differentially expressed mRNAs (DEmRNAs), lncRNAs (DElncRNAs) and miRNAs (DEmiRNAs) between TSCC samples and normal samples. In order to explore the functions of DEmRNAs, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis was performed. Subsequently, a ceRNA network was established based on the identified DElncRNAs–DEmiRNAs and DEmiRNAs–DEmRNAs interactions. The RNAs within the ceRNA network were analyzed for their correlation with overall disease survival. Finally, lncRNAs were specifically analyzed for their correlation with clinical features in the included TSCC patient samples. Results A total of 1867 mRNAs, 828 lncRNAs and 81 miRNAs were identified as differentially expressed in TSCC tissues (—log 2fold change— ≥ 2; adjusted P value <0.01). The resulting ceRNA network included 16 mRNAs, 56 lncRNAs and 6 miRNAs. Ten out of the 56 lncRNAs were found to be associated with the overall survival in TSCC patients (P < 0.05); 10 lncRNAs were correlated with TSCC progression (P < 0.05). Conclusion Our study deepens the understanding of ceRNA network regulatory mechanisms in TSCC. Furthermore, we identified ten lncRNAs (PART1, LINC00261, AL163952.1, C2orf48, FAM87A, LINC00052, LINC00472, STEAP3-AS1, TSPEAR-AS1 and ERVH48-1) as novel, potential prognostic biomarkers and therapeutic targets for TSCC.

2019 ◽  
Vol 2019 ◽  
pp. 1-14 ◽  
Author(s):  
Yan Yao ◽  
Tingting Zhang ◽  
Lingyu Qi ◽  
Ruijuan Liu ◽  
Gongxi Liu ◽  
...  

Background. Lung squamous cell carcinoma (LUSC) is a subtype of highly malignant lung cancer with poor prognosis, for which smoking is the main risk factor. However, the underlying genetic and molecular mechanisms of smoking-related LUSC remain largely unknown. Methods. We mined existing LUSC-related mRNA, miRNA, and lncRNA transcriptome data and corresponding clinical data from The Cancer Genome Atlas (TCGA) database and divided them into smoking and nonsmoking groups, followed by differential expression analysis. Functional enrichment analysis of the unique differentially expressed mRNAs of the two groups was performed using the DAVID database. Subsequently, the lncRNA-miRNA-mRNA competing endogenous RNA (ceRNA) network of LUSC in smoking and nonsmoking groups was constructed. Finally, survival analyses were performed to determine the effects of differentially expressed lncRNAs/mRNAs/miRNAs that were involved in the ceRNA network on overall survival and to discover the hub genes. Results. A total of 1696 lncRNAs, 125 miRNAs, and 3246 mRNAs and 1784 lncRNAs, 96 miRNAs, and 3229 mRNAs with differentially expressed profiles were identified in the smoking and nonsmoking groups, respectively. The ceRNA network and survival analysis revealed four lncRNAs (LINC00466, DLX6-AS1, LINC00261, and AGBL1), one miRNA (hsa-mir-210), and two mRNAs (CITED2 and ENPP4), with the potential as biomarkers for smoking-related LUSC diagnosis and prognosis. Conclusion. Taken together, our research has identified the differences in the ceRNA regulatory networks between smoking and nonsmoking LUSC, which could lay the foundation for future clinical research.


Biology ◽  
2021 ◽  
Vol 10 (5) ◽  
pp. 432
Author(s):  
Yaqin Xu ◽  
Yingying Dong ◽  
Yunhua Deng ◽  
Qianrong Qi ◽  
Mi Wu ◽  
...  

A cutaneous squamous cell carcinoma (cSCC) derived from keratinocytes is the second most common cause of non-melanoma skin cancer. The accumulation of the mutational burden of genes and cellular DNA damage caused by the risk factors (e.g., exposure to ultraviolet radiation) contribute to the aberrant proliferation of keratinocytes and the formation of a cSCC. A cSCC encompasses a spectrum of diseases that range from recursor actinic keratosis (AK) and squamous cell carcinoma (SCC) in situ (SCCIS) to invasive cSCCs and further metastatic SCCs. Emerging evidence has revealed that lncRNAs are involved in the biological process of a cSCC. According to the ceRNA regulatory theory, lncRNAs act as natural miRNA sponges and interact with miRNA response elements, thereby regulating the mRNA expression of their down-stream targets. This study was designed to search for the potential lncRNAs that may become potential therapeutic targets or biomarkers of a cSCC. Considering the spirit of the study to be adequately justified, we collected microarray-based datasets of 19 cSCC tissues and 12 normal skin samples from the GEO database (GSE42677 and GSE45164). After screening the differentially expressed genes via a limma package, we identified 24 differentially expressed lncRNAs (DElncRNAs) and 3221 differentially expressed mRNAs (DEmRNAs). The miRcode, miRTarBase, miRDB and TargetScan databases were used to predict miRNAs that could interact with DElncRNAs and DEmRNAs. A total of 137 miRNA-lncRNA and 221 miRNA-mRNA pairs were retained in the ceRNA network, consisting of 31 miRNAs, 11 DElncRNAs and 155 DEmRNAs. For the functional analysis, the top enriched biological process was enhancer sequence-specific DNA binding in Gene Ontology (GO) terms. The FoxO signaling pathway, autophagy and cellular senescence were the top enrichment terms based on a Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. The combination of a STRING tool and Cytoscape software (plug-in MCODE) identified five core mRNAs and built a core mRNA-associated ceRNA network. The expression for five identified core mRNAs and their related nine lncRNAs was validated using the external dataset GSE7553. Finally, one lncRNA HLA-F-AS1 and three mRNAs named AGO4, E2F1 and CCND1 were validated with the same expression patterns. We speculate that lncRNA HLA-F-AS1 may sponge miR-17-5p or miR-20b-5p to regulate the expression of CCND1 and E2F1 in the cSCC. The present study may provide potential diagnostic and therapeutic targets for cSCC patients.


Author(s):  
Chang Gu ◽  
Xin Shi ◽  
Wenli Qiu ◽  
Zhenyu Huang ◽  
Yan Yu ◽  
...  

BackgroundThere have been limited treatment therapies for lung squamous cell carcinoma (LUSC). M6A-related genes may be the next therapeutic targets for LUSC. In this study, we explored the prognostic role and mutational characteristics of m6A-related genes in LUSC.MethodsLUSC gene expression data, mutational data, and corresponding clinical information were extracted from The Cancer Genome Atlas database. Differentially expressed genes (DEGs) were identified, and the mutation characteristics of LUSC patients were explored. Then, m6A-related genes were extracted and the correlations among the genes were detected. Finally, the prognostic roles of the genes were investigated and the nomogram model was developed. Besides, the protein–protein interaction (PPI) network was used to explore the potential interactions among the genes.ResultsIn total, there are 551 LUSC samples enrolled in our study, containing 502 LUSC tumor samples and 49 adjacent normal LUSC samples, respectively. There were 2970 upregulated DEGs and 1806 downregulated DEGs were further explored. IGF2BP1 and RBM15 had significant co-occurrence frequency (p &lt; 0.05). Besides, METTL14 and ZC3H13 or YTHDF3 also had significant co-occurrence frequency (p &lt; 0.05). All the m6A-related genes represent the positive correlation. WTAP was identified as a prognostic gene in the TCGA database while YTHDC1 and YTHDF1 were identified as prognostic genes. In multivariate Cox analysis, YTHDF1, age, pN stage, pTNM stage, and smoking were all identified as significant prognostic factors for OS.ConclusionWe investigated the expression patterns and mutational characteristics of LUSC patients and identified three potential independent prognostic m6A-related genes (WTAP, YTHDC1, and YTHDF1) for OS in LUSC patients.


2021 ◽  
Vol 11 ◽  
Author(s):  
Heyang Cui ◽  
Yongjia Weng ◽  
Ning Ding ◽  
Chen Cheng ◽  
Longlong Wang ◽  
...  

Esophageal squamous cell carcinoma (ESCC) is one of the most aggressive malignant tumors in China, and its prognosis remains poor. Autophagy is an evolutionarily conserved catabolic process involved in the occurrence and development of ESCC. In this study, we described the expression profile of autophagy-related genes (ARGs) in ESCC and developed a prognostic prediction model for ESCC patients based on the expression pattern of ARGs. We used four ESCC cohorts, GSE53624 (119 samples) set as the discovery cohort, The Cancer Genome Atlas (TCGA) ESCC set (95 samples) as the validation cohort, 155 ESCC cohort, and Oncomine cohort were used to screen and verify differentially expressed ARGs. We identified 34 differentially expressed genes out of 222 ARGs. In the discovery cohort, we divided ESCC patients into three groups that showed significant differences in prognosis. Then, we analyzed the prognosis of 34 differentially expressed ARGs. Three genes [poly (ADP-ribose) polymerase 1 (PARP1), integrin alpha-6 (ITGA6), and Fas-associated death domain (FADD)] were ultimately obtained through random forest feature selection and were constructed as an ARG-related prognostic model. This model was further validated in TCGA ESCC set. Cox regression analysis confirmed that the three-gene signature was an independent prognostic factor for ESCC patients. This signature effectively stratified patients in both discovery and validation cohorts by overall survival (P = 5.162E-8 and P = 0.052, respectively). We also constructed a clinical nomogram with a concordance index of 0.713 to predict the survival possibility of ESCC patients by integrating clinical characteristics and the ARG signature. The calibration curves substantiated fine concordance between nomogram prediction and actual observation. In conclusion, we constructed a new ARG-related prognostic model, which shows the potential to improve the ability of individualized prognosis prediction in ESCC.


PeerJ ◽  
2019 ◽  
Vol 7 ◽  
pp. e6991 ◽  
Author(s):  
Yidan Song ◽  
Yihua Pan ◽  
Jun Liu

BackroundTongue squamous cell carcinoma (TSCC) is the most common malignant tumor in the oral cavity. An increasing number of studies have suggested that long noncoding RNA (lncRNA) plays an important role in the biological process of disease and is closely related to the occurrence and development of disease, including TSCC. Although many lncRNAs have been discovered, there remains a lack of research on the function and mechanism of lncRNAs. To better understand the clinical role and biological function of lncRNAs in TSCC, we conducted this study.MethodsIn this study, 162 tongue samples, including 147 TSCC samples and 15 normal control samples, were investigated and downloaded from The Cancer Genome Atlas (TCGA). We constructed a competitive endogenous RNA (ceRNA) regulatory network. Then, we investigated two lncRNAs as key lncRNAs using Kaplan–Meier curve analysis and constructed a key lncRNA-miRNA-mRNA subnetwork. Furthermore, gene set enrichment analysis (GSEA) was carried out on mRNAs in the subnetwork after multivariate survival analysis of the Cox proportional hazards regression model.ResultsThe ceRNA regulatory network consists of six differentially expressed miRNAs (DEmiRNAs), 29 differentially expressed lncRNAs (DElncRNAs) and six differentially expressed mRNAs (DEmRNAs). Kaplan-Meier curve analysis of lncRNAs in the TSCC ceRNA regulatory network showed that only two lncRNAs, including LINC00261 and PART1, are correlated with the total survival time of TSCC patients. After we constructed the key lncRNA-miRNA -RNA sub network, the GSEA results showed that key lncRNA are mainly related to cytokines and the immune system. High expression levels of LINC00261 indicate a poor prognosis, while a high expression level of PART1 indicates a better prognosis.


2020 ◽  
Vol 19 ◽  
pp. 153303382098578
Author(s):  
Zhibin Jing ◽  
Sitong Guo ◽  
Peng Zhang ◽  
Zheng Liang

Objective: This study aims to construct a systematic mRNA-miRNA-lncRNA network to identify novel lncRNAs and miRNAs biomarkers for laryngeal squamous cell carcinoma (LSCC). Methods: The mRNA, miRNA and lncRNA expression profiles of LSCC were obtained from Gene Expression Omnibus (GEO) database. The differentially expressed mRNAs, miRNAs and lncRNAs (DEmRNAs, DEmiRNAs and DElncRNAs) were screened between LSCC tissues and controls. Functional analysis of DEmRNAs, DEmRNAs targeted by DEmiRNAs and DEmRNAs targeted by DElncRNAs were respectively performed. The miRWalk, starbase and DIANA-LncBase were respectively used to predict DEmiRNAs-DEmRNAs, DElncRNAs-DEmRNAs and DElncRNAs-DEmiRNAs pairs. ceRNA network was built by DEmiRNAs-DEmRNAs and DElncRNAs-DEmiRNAs pairs. LncRNA subcellular localization was predicted using lncLocator. Using published The Cancer Genome Atlas (TCGA) and external datasets (GSE127165 and GSE133632), we also validated the expression of key DElncRNAs and DEmiRNAs in ceRNA network. The diagnostic and prognostic value of candidate genes was evaluated by ROC curve analysis and survival analysis, respectively. Results: There were 5 mRNA datasets, 3 miRNA datasets and 2 lncRNA datasets in this study. Totally, 2957 DEmRNAs, 61 DElncRNAs and 23 DEmiRNAs were identified. Functional analysis of DEmRNAs shows that they were significantly enriched in cancer-related pathways, such as DNA replication and extracellular matrix organization. There were 11 DEmiRNAs, 17 DElncRNAs and 967 DEmRNAs in the ceRNA network. Notably, up-regulated lncRNA DGCR5-down-regulated has-miR-338-3p/has-miR-139-5p pairs in this network were experimentally validated. Moreover, down-regulated AL121839.2, down-regulated LINC02147, up-regulated AC079328.2, up-regulated AC004943.2 and up-regulated HMGA2-AS1 were located in the cytoplasm. AL121839.2 and LINC02147 interacted with has-miR-1246. AC004943.2, AC079328.2 and HMGA2-AS1 targeted has-miR-3185, has-miR-3137 and has-miR-582-5p, respectively. Based on the TCGA and external datasets (GSE127165 and GSE133632), DGCR5 and AC004943.2 were significantly up-regulated while AL121839.2 and LINC02147, has-miR-338-3p, has-miR-139-5p and has-miR-582-5p were significantly down-regulated, which were consistent with our integration analysis. DGCR5, AL121839.2, LINC02147, AC004943.2, has-miR-338-3p, has-miR-139-5p and has-miR-582-5p could predict the occurrence of LSCC. Survival analysis suggested that only, AL121839.2 has potential prognostic value for LSCC. Conclusion: This study provided novel insights into the ceRNA network and uncovered novel lncRNAs and miRNAs with diagnostic value in LSCC.


PeerJ ◽  
2019 ◽  
Vol 7 ◽  
pp. e7727 ◽  
Author(s):  
Lingyu Qi ◽  
Tingting Zhang ◽  
Yan Yao ◽  
Jing Zhuang ◽  
Cun Liu ◽  
...  

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.


2020 ◽  
Vol 22 (1) ◽  
pp. 60
Author(s):  
Sichong Han ◽  
Zhe Wang ◽  
Jining Liu ◽  
Qipeng Yuan

Understanding the mechanism by which sulforaphene (SFE) affects esophageal squamous cell carcinoma (ESCC) contributes to the application of this isothiocyanate as a chemotherapeutic agent. Thus, we attempted to investigate SFE regulation of ESCC characteristics more deeply. We performed gene set enrichment analysis (GSEA) on microarray data of SFE-treated ESCC cells and found that differentially expressed genes are enriched in TNFα_Signaling_via_the_NFκB_Pathway. Coupled with the expression profile data from the GSE20347 and GSE75241 datasets, we narrowed the set to 8 genes, 4 of which (C-X-C motif chemokine ligand 10 (CXCL10), TNF alpha induced protein 3 (TNFAIP3), inhibin subunit beta A (INHBA), and plasminogen activator, urokinase (PLAU)) were verified as the targets of SFE. RNA-sequence (RNA-seq) data of 182 ESCC samples from The Cancer Genome Atlas (TCGA) were grouped into two phenotypes for GSEA according to the expression of CXCL10, TNFAIP3, INHBA, and PLAU. The enrichment results proved that they were all involved in the NFκB pathway. ChIP-seq analyses obtained from the Cistrome database indicated that NFκB-p65 is likely to control the transcription of CXCL10, TNFAIP3, INHBA, and PLAU, and considering TNFAIP3 and PLAU are the most significantly differentially expressed genes, we used chromatin immunoprecipitation-polymerase chain reaction (ChIP-PCR) to verify the regulation of p65 on their expression. The results demonstrated that SFE suppresses ESCC progression by down-regulating TNFAIP3 and PLAU expression in a p65-dependent manner.


Sign in / Sign up

Export Citation Format

Share Document