Bioinformatics analysis of laryngeal squamous cell carcinoma: seeking key candidate genes and pathways.
Ontology highlight
ABSTRACT: Background
Laryngeal squamous cell carcinoma (LSCC) is the second most aggressive head and neck squamous cell carcinoma. Although much work has been done to optimize its treatment, patients with LSCC still have poor prognosis. Therefore, figuring out differentially expressed genes (DEGs) contained in the progression of LSCC and employing them as potential therapeutic targets or biomarkers for LSCC is extremely meaningful.Methods
Overlapping DEGs were screened from two standalone Gene Expression Omnibus datasets, and Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed. By applying STRING and Cytoscape, a protein-protein network was built, and module analysis was carried out. The hub genes were selected by maximal clique centrality with the CytoHubba plugin of Cytoscape. UALCAN and GEPIA data were examined to validate the gene expression findings. Moreover, the connection of the hub genes with LSCC patient overall survival was studied employing The Cancer Genome Atlas. Then, western blot, qRT-PCR, CCK-8, wound healing and transwell assays were bring to use for further verify the key genes.Results
A total of 235 DEGs were recorded, including 83 upregulated and 152 downregulated genes. A total of nine hub genes that displayed a high degree of connectivity were selected. UALCAN and GEPIA databases verified that these genes were highly expressed in LSCC tissues. High expression of the SPP1, SERPINE1 and Matrix metalloproteinases 1 (MMP1) genes was connected to worse prognosis in patients with LSCC, according to the GEPIA online tool. Western blot and qRT-PCR testify SPP1, SERPINE1 and MMP1 were upregulated in LSCC cells. Inhibition of SPP1, SERPINE1 and MMP1 suppressed cell proliferation, invasion and migration.Conclusion
The work here identified effective and reliable diagnostic and prognostic molecular biomarkers by unified bioinformatics analysis and experimental verification, indicating novel and necessary therapeutic targets for LSCC.
SUBMITTER: Ma J
PROVIDER: S-EPMC8052978 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Oral squamous cell carcinoma (OSCC) is one of the most common types of malignant head and neck tumor, which poses a serious threat to human health. In recent years, the incidence of OSCC has been increasing, while the prognosis has not significantly improved. Elucidation of the molecular mechanisms underlying the development of OSCC may provide novel therapeutic strategies. In the present study, the gene expression profiles from 4 datasets, including 244 OSCC and 95 normal oral mucosa samples, were subjected to statistical and Bioinformatics analysis. A total of 34 differentially expressed genes (DEGs) were identified, among which 14 were upregulated and 20 were downregulated in OSCC compared with normal oral mucosa tissues. Gene Ontology enrichment analysis indicated that the DEGs were mainly involved in regulation of the immune response, cell adhesion and cell proliferative processes. The Kyoto Encyclopedia of Genes and Genomes pathway analysis revealed that the DEGs were mainly associated with the phosphoinositide-3 kinase Akt and Toll-like receptor signaling pathway. The key candidate DEGs were identified from the complex protein-protein interaction network, and secreted phosphoprotein 1 (SPP1), integrin subunit α 3 and plasminogen activator, urokinase (PLAU) were confirmed to be significantly associated with the survival rate. Cell Counting Kit-8 and Transwell assays demonstrated that SPP1 and PLAU regulate cell proliferation, migration and invasion. The candidate genes/pathways identified in the present study may include promising diagnostic biomarkers or therapeutic targets for OSCC.
| S-EPMC6468404 | biostudies-literature
Project description:Ewing's sarcoma (ES) is the second most common bone tumor among children and adolescents worldwide. However, the genes and signaling pathways involved in ES tumorigenesis and progression remain unclear. The present study used two gene-expression profile datasets (GSE17674 and GSE31215) to elucidate key potential candidate genes and pathways in ES. Differentially expressed genes (DEGs) were identified and a functional enrichment analysis was performed. A protein-protein interaction (PPI) network was constructed, and the most significant module in the PPI network was selected from the Search Tool for the Retrieval of Interacting Genes/Proteins database. A total of 278 genes were identified by comparing the tumor samples with non-cancerous samples; these included 272 upregulated and 6 downregulated genes. The pathway analysis demonstrated significant enrichment in the positive regulation of transcription in the DEGs coding for RNA polymerase II promoter, plasma membrane and chromatin binding pathways in cancer in general. There were 269 nodes and 292 edges in the PPI network. Finally, MYC, IGF1, OAS1, EZH2 and ISG15 were identified as the hub genes according to the degree levels. The survival analysis revealed that EZH2 is associated with a poor prognosis in patients with ES. In conclusion, the DEGs, associated pathways and hub genes identified in the present study help elucidate the underlying molecular mechanisms of ES carcinogenesis and progression, and provide potential molecular targets and biomarkers for ES.
| S-EPMC6865160 | biostudies-literature
Project description:The aim of the present study was to identify the differentially expressed genes (DEGs) in esophageal squamous-cellcarcinoma (ESCC) and provide potential therapeutic targets. The microarray dataset GSE20347 was downloaded from the Gene Expression Omnibus (GEO) database, and included 17 tissue samples and 13 normal adjacent tissue samples from patients with ESCC. A total of 22,277 DEGs were identified. A heat map for the DEGs was constructed with the Morpheus online tool and the top 200 genes (100 upregulated and 100 downregulated) were selected for further bioinformatics analysis, including analysis of gene ontology (GO) terms, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, protein-protein interaction networks and Spearman's correlation tests. The results of the GO analysis indicated that the upregulated DEGs were most significantly enriched in membrane-bounded vesicles in the cellular component (CC) category, but were not significantly enriched in any GO terms of the categories biological process (BP) or molecular function (MF); furthermore, the downregulated DEGs were most significantly enriched in regulation of DNA metabolic processes, nucleotide binding and chromosomes in the categories BP, MF and CC, respectively. The KEGG analysis indicated that the downregulated DEGs were enriched in the regulation of cell cycle pathways. The top 10 hub proteins in the protein-protein interaction network were cyclin-dependent kinase 4, budding uninhibited by benzimidazoles 1, cyclin B2, heat shock protein 90AA1, aurora kinase A, H2A histone family member Z, replication factor C subunit 4, and minichromosome maintenance complex component 2, -4 and -7. These proteins are mainly involved in regulating tumor progression. The genes in the four top modules were mainly implicated in regulating cell cycle pathways. Secreted Ly-6/uPAR-related protein (SLURP) was the hub gene, and SLURP and its interacting genes were most enriched in the chromosomal part in the CC category, organelle organization in the BP category and protein binding in the MF category, and were involved in pathways including DNA replication, cell cycle and P53 signaling. The expression of SLURP-1 in fifteen patients with esophageal carcinoma was detected using quantitative polymerase chain reaction analysis, and the results indicated that SLURP-1 expression was significantly decreased in the tumor samples relative to that in normal adjacent tissues. These results suggest that several hub proteins and the hub gene SLURP-1 may serve as potential therapeutic targets, and that gene dysfunction may be involved in the tumorigenesis of ESCC.
| S-EPMC6090437 | biostudies-literature
Project description:BackgroundLaryngeal squamous cell carcinoma (LSCC) is the prominent cancer in head and neck, which greatly affects life quality of patients. The pathogenesis of LSCC is not clear. Presently, the LSCC treatments include chemotherapy, surgery and radiotherapy; however, these methods have poor efficacy in patients with recurrent and persistent cancer. Therefore, the study identified the hub genes accompanied with LSCC, which may be a potential therapeutic target in the future.MethodsWe extracted whole transcriptome high-throughput sequencing (HTS) LSCC data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases and calculate differentially expressed genes (DEGs) between LSCC and normal samples using statistical software RStudio. Through weighted gene co-expression network analysis (WGCNA), enrichment examination of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and Gene Ontology (GO) functions, and examination of protein-protein interaction (PPI) network, we obtained network hub genes and validated the hub genes prognostic value and expression levels of protein.ResultsThrough analysis of differential gene expression, from the GEO and TCGA databases 2,139 and 2,774 DEGs were obtained, respectively, 13 and 15 modules were screened from TCGA-LSCC and GSE127165 datasets by WGCNA, respectively. The most significant positive and negative correlation modules in the WGCNA and DEG lists were overlapped, and overall 36 co-expressed overlapping genes were retrieved. Through enrichment analysis of GO and KEGG, it was found that the gene functions were highly concentrated in cell junction assembly, basement membrane, extracellular matrix (ECM) structural constituent etc., and the pathways were mainly concentrated in ECM receptor interaction, focal adhesion, small cell lung cancer, and toxoplasmosis. Through analysis of PPI network analysis, 10 network hub genes (SNAI2, ITGA6, LAMB3, LAMC2, CAV1, COL7A1, GJA1, EHF, OAT, and GPT) were obtained. Finally, survival analysis and protein expression validation of these genes confirmed that low OAT expression and high CAV1 expression remarkably influenced the survival of patient's prognosis with LSCC.ConclusionsWe recognized the hub genes and key modules nearly associated to LSCC and these genes were validated by survival analysis and the database of Human Protein Atlas (HPA), which is of high importance for unveiling the pathogenesis of LSCC and probing for new precise biological marker and potential therapeutic targets.
| S-EPMC11319952 | biostudies-literature
Project description:The purpose of the present study was to screen the key genes associated with laryngeal carcinoma and to investigate the molecular mechanism of laryngeal carcinoma progression. The gene expression profile of GSE10935 [Gene Expression Omnibus (GEO) accession number], including 12 specimens from laryngeal papillomas and 12 specimens from normal laryngeal epithelia controls, was downloaded from the GEO database. Differentially expressed genes (DEGs) were screened in laryngeal papillomas compared with normal controls using Limma package in R language, followed by Gene Ontology (GO) enrichment analysis and pathway enrichment analysis. Furthermore, the protein-protein interaction (PPI) network of DEGs was constructed using Cytoscape software and modules were analyzed using MCODE plugin from the PPI network. Furthermore, significant biological pathway regions (sub-pathway) were identified by using iSubpathwayMiner analysis. A total of 67 DEGs were identified, including 27 up-regulated genes and 40 down-regulated genes and they were involved in different GO terms and pathways. PPI network analysis revealed that Ras association (RalGDS/AF-6) domain family member 1 (RASSF1) was a hub protein. The sub-pathway analysis identified 9 significantly enriched sub-pathways, including glycolysis/gluconeogenesis and nitrogen metabolism. Genes such as phosphoglycerate kinase 1 (PGK1), carbonic anhydrase II (CA2), and carbonic anhydrase XII (CA12) whose node degrees were >10 were identified in the disease risk sub-pathway. Genes in the sub-pathway, such as RASSF1, PGK1, CA2 and CA12 were presumed to serve critical roles in laryngeal carcinoma. The present study identified DEGs and their sub-pathways in the disease, which may serve as potential targets for treatment of laryngeal carcinoma.
| S-EPMC4950495 | biostudies-literature
Project description:Esophageal squamous cell carcinoma (ESCC) ranks as the fourth leading cause of cancer-related death in China. Although paclitaxel has been shown to be effective in treating ESCC, the prolonged use of this chemical will lead to paclitaxel resistance. In order to uncover genes and pathways driving paclitaxel resistance in the progression of ESCC, bioinformatics analyses were performed based on The Cancer Genome Atlas (TCGA) database and the Gene Expression Omnibus (GEO) database including GSE86099 and GSE161533. Differential expression analysis was performed in TCGA data and two GEO datasets to obtain differentially expressed genes (DEGs). Based on GSE161533, weighted gene co-expression network analysis (WGCNA) was conducted to identify the key modules associated with ESCC tumor status. The DEGs common to the two GEO datasets and the genes in the key modules were intersected to obtain the paclitaxel resistance-specific or non-paclitaxel resistance-specific genes, which were subjected to subsequent least absolute shrinkage and selection operator (LASSO) feature selection, whereby paclitaxel resistance-specific or non-paclitaxel resistance-specific key genes were selected. Ten machine learning models were used to validate the biological significance of these key genes; the potential therapeutic drugs for paclitaxel resistance-specific genes were also predicted. As a result, we identified 24 paclitaxel resistance-specific genes and 18 non-paclitaxel resistance-specific genes. The ESCC machine classifiers based on the key genes achieved a relatively high AUC value in the cross-validation and in an independent test set, GSE164158. A total of 207 drugs (such as bevacizumab) were predicted to be alternative therapeutics for ESCC patients with paclitaxel resistance. These results might shed light on the in-depth research of paclitaxel resistance in the context of ESCC progression.
| S-EPMC8386171 | biostudies-literature
Project description:Laryngeal squamous cell carcinoma (LSCC) is the most common malignant tumor in the head and neck, and can seriously affect the daily life of patients. To study the mechanisms of LSCC, the microarray of GSE51958 was analyzed in the present study. GSE51958 was downloaded from Gene Expression Omnibus, and included a collection of LSCC tissue samples and matched adjacent non-cancerous tissue samples from 10 patients. Differentially-expressed genes (DEGs) were identified using limma package. Next, a weighted co-expression network was constructed for the DEGs by WGCNA package in R. Modules of the weighted co-expression network were obtained through constructing a hierarchical clustering tree using the hybrid dynamic shear tree method. Using the clusterProfiler package, the potential functions of DEGs in the modules correlated with LSCC were predicted by pathway enrichment analysis. In total, 959 DEGs were screened from the LSCC samples compared with the adjacent non-cancerous samples, including 553 upregulated and 406 downregulated genes. The appointed black, brown, gray, pink and yellow modules were screened for the DEGs in the weighted co-expression network. For the DEGs in the brown and yellow modules, the enriched pathways were cytokine-cytokine receptor interaction and metabolic pathways, respectively. The DEGs in the pink module were involved in the majority of pathways. With high connectivity degrees in the pink module, TPX2, microtubule-associated (TPX2; degree, 25), minichromosome maintenance complex component 2 (MCM2; degree, 25), ubiquitin-like with PHD and ring finger domains 1 (UHRF1; degree, 22), cyclin-dependent kinase 2 (CDK2; degree, 20) and protein regulator of cytokinesis 1 (PRC1; degree, 20) may be involved in LSCC. Overall, In conclusion, from the integrated bioinformatics analysis of genes that may be associated with LSCC, 959 DEGs were obtained from LSCC samples compared with adjacent non-cancerous samples, and TPX2, MCM2, UHRF1, CDK2 and PRC1 were found to hold a possible association with the disease.
| S-EPMC4840875 | biostudies-literature
Project description:Lung squamous cell carcinoma (LUSC) is often diagnosed at the advanced stage with poor prognosis. The mechanisms of its pathogenesis and prognosis require urgent elucidation. This study was performed to screen potential biomarkers related to the occurrence, development and prognosis of LUSC to reveal unknown physiological and pathological processes. Using bioinformatics analysis, the lung squamous cell carcinoma microarray datasets from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases were analyzed to identify differentially expressed genes (DEGs). Furthermore, PPI and WGCNA network analysis were integrated to identify the key genes closely related to the process of LUSC development. In addition, survival analysis was performed to achieve a prognostic model that accomplished good prediction accuracy. Three hundred and thirty-seven up-regulated and 119 down-regulated genes were identified, in which four genes have been found to play vital roles in LUSC development, namely CCNA2, AURKA, AURKB, and FEN1. The prognostic model contained 5 genes, which were all detrimental to prognosis. The AUC of the established prognostic model for predicting the survival of patients at 1, 3, and 5 years was 0.692, 0.722, and 0.651 in the test data, respectively. In conclusion, this study identified several biomarkers of significant interest for additional investigation of the therapies and methods of prognosis of lung squamous cell carcinoma.
| S-EPMC7215920 | biostudies-literature
Project description:Hepatocellular carcinoma (HCC) accounts for approximately 85-90% of all liver cancer cases and has poor relapse-free survival. There are many gene expression studies that have been performed to elucidate the genetic landscape and driver pathways leading to HCC. However, existing studies have been limited by the sample size and thus the pathogenesis of HCC is still unclear. In this study, we performed an integrated characterization using four independent datasets including 320 HCC samples and 270 normal liver tissues to identify the candidate genes and pathways in the progression of HCC. A total of 89 consistent differentially expression genes (DEGs) were identified. Gene-set enrichment analysis revealed that these genes were significantly enriched for cellular response to zinc ion in biological process group, collagen trimer in the cellular component group, extracellular matrix (ECM) structural constituent conferring tensile strength in the molecular function group, protein digestion and absorption, mineral absorption and ECM-receptor interaction. Network system biology based on the protein-protein interaction (PPI) network was also performed to identify the most connected and important genes based on our DEGs. The top five hub genes including osteopontin (SPP1), Collagen alpha-2(I) chain (COL1A2), Insulin-like growth factor I (IGF1), lipoprotein A (LPA), and Galectin-3 (LGALS3) were identified. Western blot and immunohistochemistry analysis were employed to verify the differential protein expression of hub genes in HCC patients. More importantly, we identified that these five hub genes were significantly associated with poor disease-free survival and overall survival. In summary, we have identified a potential clinical significance of these genes as prognostic biomarkers for HCC patients who would benefit from experimental approaches to obtain optimal outcome.
| S-EPMC7396661 | biostudies-literature
Project description:Breast cancer (BC) is the leading malignancy in women worldwide, yet relatively little is known about the genes and signaling pathways involved in BC tumorigenesis and progression. The present study aimed to elucidate potential key candidate genes and pathways in BC. Five gene expression profile data sets (GSE22035, GSE3744, GSE5764, GSE21422 and GSE26910) were downloaded from the Gene Expression Omnibus (GEO) database, which included data from 113 tumorous and 38 adjacent non?tumorous tissue samples. Differentially expressed genes (DEGs) were identified using t?tests in the limma R package. These DEGs were subsequently investigated by pathway enrichment analysis and a protein?protein interaction (PPI) network was constructed. The most significant module from the PPI network was selected for pathway enrichment analysis. In total, 227 DEGs were identified, of which 82 were upregulated and 145 were downregulated. Pathway enrichment analysis results revealed that the upregulated DEGs were mainly enriched in 'cell division', the 'proteinaceous extracellular matrix (ECM)', 'ECM structural constituents' and 'ECM?receptor interaction', whereas downregulated genes were mainly enriched in 'response to drugs', 'extracellular space', 'transcriptional activator activity' and the 'peroxisome proliferator?activated receptor signaling pathway'. The PPI network contained 174 nodes and 1,257 edges. DNA topoisomerase 2?a, baculoviral inhibitor of apoptosis repeat?containing protein 5, cyclin?dependent kinase 1, G2/mitotic?specific cyclin?B1 and kinetochore protein NDC80 homolog were identified as the top 5 hub genes. Furthermore, the genes in the most significant module were predominantly involved in 'mitotic nuclear division', 'mid?body', 'protein binding' and 'cell cycle'. In conclusion, the DEGs, relative pathways and hub genes identified in the present study may aid in understanding of the molecular mechanisms underlying BC progression and provide potential molecular targets and biomarkers for BC.
| S-EPMC5983982 | biostudies-literature