Analysis of TERT Isoforms across TCGA, GTEx and CCLE Datasets.
Ontology highlight
ABSTRACT:
SUBMITTER: Subasri M
PROVIDER: S-EPMC8070023 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundAbnormal expression of long non-coding RNAs (lncRNA) play a significant role in the incidence and progression of high-grade serous ovarian cancer (HGSOC), which is a leading cause of mortality among gynecologic malignant tumor patients. In this study, our aim is to identify lncRNA-associated competing endogenous RNA (ceRNA ) axes that could define more reliable prognostic parameters of HGSOC, and to investigate the lncRNAs' potential mechanism of in lymphocyte infiltration.MethodsThe RNA-seq and miRNA expression profiles were downloaded from The Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression (GTEx) database; while for obtaining the differentially expressed lncRNAs (DELs), miRNAs (DEMs), and genes (DEGs), we used edgeR, limma and DESeq2. After validating the RNA, miRNA and gene expressions, using integrated three RNA expression profiles (GSE18520, GSE27651, GSE54388) and miRNA profile (GSE47841) from the Gene Expression Omnibus (GEO) database, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genome (KEGG) pathway analyses through ClusterProfiler. The prognostic value of these genes was determined with Kaplan-Meier survival analysis and Cox regression analysis. The ceRNA network was constructed using Cytoscape. The correlation between lncRNAs in ceRNA network and immune infiltrating cells was analyzed by using Tumor IMmune Estimation Resource (TIMER), and gene markers of tumor-infiltrating immune cells were identified using Spearman's correlation after removing the influence of tumor purity.ResultsA total of 33 DELs (25 upregulated and eight downregulated), 134 DEMs (76 upregulated and 58 downregulated), and 1,612 DEGs (949 upregulated and 663 downregulated) were detected that could be positively correlated with overall survival (OS) of HGSOC. With the 1,612 analyzed genes, we constructed a ceRNA network, which indicated a pre-dominant involvement of the immune-related pathways. Furthermore, our data revealed that LINC00665 influenced the infiltration level of macrophages and dendritic cells (DCs). On the other hand, FTX and LINC00665, which may play their possible roles through the ceRNA axis, demonstrated a potential to inhibit Tregs and prevent T-cell exhaustion.ConclusionWe defined several prognostic biomarkers for the incidence and progression of HGSOC and constructed a network for ceRNA axes; among which three were indicated to have a positive correlation with lymphocyte infiltration, namely: FTX-hsa-miR-150-5p-STK11, LINC00665-hsa-miR449b-5p-VAV3 and LINC00665-hsa-miR449b-5p-RRAGD.
| S-EPMC7211406 | biostudies-literature
Project description:Advances in high-throughput sequencing technologies now yield unprecedented volumes of OMICs data with opportunities to conduct systematic data analyses and derive novel biological insights. Here, we provide protocols to perform differential-expressed gene analysis of TCGA and GTEx RNA-Seq data from human cancers, complete integrative GO and network analyses with focus on clinical and survival data, and identify differential correlation of trait-associated biomarkers. For complete details on the use and execution of this protocol, please refer to Chen and MacDonald (2021).
| S-EPMC8841814 | biostudies-literature
Project description:We often encounter situations in which data from the TCGA that have been analyzed in papers we read or reviewed cannot be reproduced, even when TCGA datasets are used, especially in survival analyses. Therefore, we attempted to confirm the data source for TCGA survival analysis and found that several websites used to analyze the survival data of TCGA datasets inappropriately handle the survival data, causing differences in statistical analyses. This causes the misinterpretation of results because figures of survival analysis results in several papers are sometimes exactly as generated by these sites, and the results depend on only the tools provided by these sites. We would like to make this situation widely known and raise the problem for scientific soundness.
| S-EPMC8726696 | biostudies-literature
Project description:Gliomas are among the most common intracranial tumors which originated from neuroepithelial cells. Increasing evidence has revealed that long noncoding RNA (lncRNA)-microRNA (miRNA)-mRNA module regulation and tumor-infiltrating immune cells play important regulatory roles in the occurrence and progression of gliomas. However, the precise underlying molecular mechanisms remain largely unknown. Data on gliomas in The Cancer Genome Atlas lack normal control samples; to overcome this limitation, we combined 665 The Cancer Genome Atlas glioma RNA sequence datasets with 188 Genotype-Tissue Expression normal brain RNA sequences to construct an expression matrix profile after normalization. We systematically analyzed the expression of mRNAs, lncRNAs, and miRNAs between gliomas and normal brain tissues. Kaplan-Meier survival analyses were conducted to screen differentially expressed mRNAs, lncRNAs, and miRNAs. A prognostic miRNA-related competitive endogenous RNA network was constructed, and the core subnetworks were filtered using 6 miRNAs, 3 lncRNAs, and 11 mRNAs. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses were performed to investigate the biological functions of significantly dysregulated mRNAs. Co-expression network analysis was performed to analyze and screen the core genes. Furthermore, single-sample Gene Set Enrichment Analysis and immune checkpoint gene expression analysis were performed, as co-expression analysis indicated immune gene dysregulation in glioma. Finally, the expression of representative dysregulated genes was validated in U87 cells at the transcriptional level, establishing a foundation for further research. We identified 7017 mRNAs, 437 lncRNAs, and 9 miRNAs that were differentially expressed in gliomas. Kaplan-Meier survival analysis revealed 5684 mRNAs, 61 lncRNAs, and 7 miRNAs with potential as prognostic signatures in patients with glioma. The hub subnetwork of the competing endogenous RNA network between PART1-hsa-mir-25-SLC12A5/TACC2/BSN/TLN2/ZDHHC8 was screened out. Gene co-expression network, single-sample Gene Set Enrichment Analysis, and immune checkpoint expression analysis demonstrated that tumor-infiltrating immune cells are closely related to gliomas. We identified novel potential biomarkers to predict survival and therapeutic targets for patients with gliomas based on a large-scale sample. Importantly, we filtered pivotal genes that provide valuable information for further exploration of the molecular mechanisms underlying glioma tumorigenesis and progression.
| S-EPMC8516354 | biostudies-literature
Project description:Background: Recently, RNA-binding proteins (RBPs) were reported to interact with target mRNA to regulate gene posttranscriptional expression, and RBP-mediated RNA modification can regulate the expression and function of proto-oncogenes and tumor suppressor genes. We systematically analyzed the expression of RBPs in pancreatic adenocarcinoma (PAAD) and constructed an RBP-associated prognostic risk model. Methods: Gene expression data of normal pancreatic samples as well as PAAD samples were downloaded from TCGA-PAAD and GTEx databases. Wilcoxon test and univariate Cox analysis were, respectively, applied to screen differential expression RBPs (DE-RBPs) and prognostic-associated RBPs (pRBPs). Functional enrichment was analyzed by GO, KEGG, and GSEA. Protein-protein interaction (PPI) network was constructed by STRING online database. Modeling RBPs were selected by multivariate Cox analysis. Kaplan-Meier survival and Cox analysis were applied to evaluate the effects of risk score on the overall survival of PAAD patients. ROC curves and validation cohort were applied to verify the accuracy of the model. Nomogram was applied for predicting 1-, 3-, and 5-year overall survival (OS) of PAAD patients. At last, modeling RBPs were further analyzed to explore their differential expression, prognostic value, as well as enrichment pathways in PAAD. Results: RBPs (453) were differentially expressed in normal and tumor samples, besides, 28 of which were prognostic associated. DE-RBPs (453) are functionally associated with ribosome, ribonuclease, spliceosome, etc. Eight RBPs (PABPC1, PRPF6, OAS1, RBM5, LSM12, IPO7, FXR1, and RBM6) were identified to construct a prognostic risk model. Higher risk score not only predicted poor prognosis but also was an independent poor prognostic indicator, which was verified by ROC curves and validation cohort. Eight modeling RBPs were confirmed to be significantly differentially expressed between normal and tumor samples from RNA and protein level. Besides, all of eight RBPs were related with overall survival of PAAD patients. Conclusions: We successfully constructed an RBP-associated prognostic risk model in PAAD, which has a potential clinical application prospect.
| S-EPMC7873872 | biostudies-literature
Project description:The Cancer Genome Atlas (TCGA) projects have advanced our understanding of the driver mutations, genetic backgrounds, and key pathways activated across cancer types. Analysis of TCGA datasets have mostly focused on somatic mutations and translocations, with less emphasis placed on gene amplifications. Here we describe a bioinformatics screening strategy to identify putative cancer driver genes amplified across TCGA datasets. We carried out GISTIC2 analysis of TCGA datasets spanning 16 cancer subtypes and identified 486 genes that were amplified in two or more datasets. The list was narrowed to 75 cancer-associated genes with potential "druggable" properties. The majority of the genes were localized to 14 amplicons spread across the genome. To identify potential cancer driver genes, we analyzed gene copy number and mRNA expression data from individual patient samples and identified 42 putative cancer driver genes linked to diverse oncogenic processes. Oncogenic activity was further validated by siRNA/shRNA knockdown and by referencing the Project Achilles datasets. The amplified genes represented a number of gene families, including epigenetic regulators, cell cycle-associated genes, DNA damage response/repair genes, metabolic regulators, and genes linked to the Wnt, Notch, Hedgehog, JAK/STAT, NF-KB and MAPK signaling pathways. Among the 42 putative driver genes were known driver genes, such as EGFR, ERBB2 and PIK3CA. Wild-type KRAS was amplified in several cancer types, and KRAS-amplified cancer cell lines were most sensitive to KRAS shRNA, suggesting that KRAS amplification was an independent oncogenic event. A number of MAP kinase adapters were co-amplified with their receptor tyrosine kinases, such as the FGFR adapter FRS2 and the EGFR family adapters GRB2 and GRB7. The ubiquitin-like ligase DCUN1D1 and the histone methyltransferase NSD3 were also identified as novel putative cancer driver genes. We discuss the patient tailoring implications for existing cancer drug targets and we further discuss potential novel opportunities for drug discovery efforts.
| S-EPMC4038530 | biostudies-literature
Project description:A routine pipeline seems very common in many cancer studies that expression differentiation might be helpful in identifying prognostic molecules. There also exists a striking unanimity that molecules upregulated in cancer usually shorten survival, while downregulated ones have the opposite effect. In this study, based on the transcriptional profiles of 18 malignancies, cancer and corresponding adjacent normal tissues were used to calculate differential scores. Cox correlation coefficients of global genes were also calculated to denote survival association. The relationship between expression differentiation and survival association has been extensively studied in 18 malignancy types. Contradictory to our stereotypic research pattern, expression differentiation between cancer and adjacent normal tissues was proven irrelevant to corresponding survival correlation. Surprisingly, the more stringent cutoff we used in differentially expressed gene identification, the less prognostic information we would obtain from the collected gene groups. Moreover, the direction of dysregulated genes in cancer was irrelevant to the direction of corresponding survival correlation. Cancer-normal expression differentiation is irrelevant to genes' survival correlation in multiple cancers and, therefore, not helpful in identifying prognostic genes. For future studies, it is more sensible to look into another alternative rather than collect differentially expressed molecules in the initial step.
| S-EPMC5992444 | biostudies-literature
Project description:Colorectal cancer (CRC) is a major cause of cancer deaths worldwide. Unfortunately, many CRC patients are still being diagnosed at an advanced stage of the cancer, and the 5-year survival rate is only ~30%. Effective prognostic markers of CRC are therefore urgently needed. To address this issue, we performed a detailed bioinformatics analysis based on the Cancer Genome Atlas (TCGA), Genotype-Tissue Expression (GTEx), and Gene Expression Omnibus (GEO) databases to identify prognostic biomarkers for CRC, which in turn help in exploring potential drug-repurposing. We identified five hub genes (PGM2, PODXL, RHNO1, SCD, and SEPHS1), which had good performance in survival prediction and might be involved in CRC through three key pathways ("Cell cycle," "Purine metabolism," and "Spliceosome" KEGG pathways) identified by a KEGG pathway enrichment analysis. What is more, we performed a co-expression analysis between five hub genes and transcription factors to explore the upstream regulatory region. Furthermore, we screened the potential drug-repurposing for the five hub genes in CRC according to the Binding DB and ZINC15 databases. Taking together, we constructed a five-gene signature to predict overall survival of CRC and found the potential drug-repurposing, which may improve the outcome of CRC in the future.
| S-EPMC7848190 | biostudies-literature
Project description:AbsractBACKGROUND: Several studies have investigated the associations between the podocalyxin-like protein (PODXL) expression quantity or locations and cancers survival, but the results were far from conclusive. Therefore, we proceeded a meta-analysis on PODXL in various human cancers to find its prognostic value and followed confirmation using the TCGA datasets.MethodsWe performed a systematic search, and 18 citations, including 5705 patients were pooled in meta-analysis. The results were verified with TCGA datasets.ResultsTotal eligible studies comprised 5705 patients with 10 types of cancer. And the result indicated that PODXL high-expression or membrane-expression were significantly related to poor overall survival (OS). However, subgroup analysis showed a significant association between high expressed PODXL and poor OS in the colorectal cancer, pancreatic cancer, urothelial bladder cancer, renal cell carcinoma and glioblastoma multiforme. Then, we validated the inference using TCGA datasets, and the consistent results were demonstrated in patients with pancreatic cancer, glioblastoma multiforme, gastric cancer, esophageal cancer and lung adenocarcinoma.ConclusionThe result of meta-analysis showed that high expressed PODXL was significantly linked with poor OS in pancreatic cancer and glioblastoma multiforme, but not in gastric cancer, esophageal cancer or lung adenocarcinoma. And the membrane expression of PODXL might also associate with poor OS. PODXL may act as tumor promotor and may serve as a potential target for antitumor therapy.
| S-EPMC7331259 | biostudies-literature
Project description:Aim: Existing targeted therapies for hepatocellular carcinoma (HCC) are resistant and have limitations. It is crucial to find new HCC-related target genes. Methods: RNA-sequencing data of HCC were gathered from The Cancer Genome Atlas and Gene Expression Omnibus datasets. Initially, differentially expressed genes between normal and tumor tissues were identified from four Gene Expression Omnibus datasets, GSE36376, GSE102079, GSE54236, and GSE45267. GO terms and KEGG pathway enrichment analyses were performed to explore the potential biological functions of differentially expressed genes. A PPI network was constructed by using the STRING database, and up-regulated and down-regulated hub genes were defined through 12 topological approaches. Subsequently, the correlation bounded by up-regulated genes and down-regulated genes in the diagnosis, prognosis, and clinicopathological features of HCC was analyzed. Beyond a shadow of doubt, the key oncogene PBK and tumor suppressor gene F9 were screened out, and the specific mechanism was investigated through GSEA enrichment analysis and immune correlation analysis. The role of PBK in HCC was further verified by western blot, CCK8, transwell, and tube formation experiments. Results: CDCA5, CDC20, PBK, PRC1, TOP2A, and NCAPG are good indicators of HCC diagnosis and prognosis. The low expressions of F9, AFM, and C8B indicate malignant progression and poor prognosis of HCC. PBK was found to be closely related to VEGF, VEGFR, and PDGFR pathways. Experiments showed that PBK promotes HCC cell proliferation, migration, invasion, and tube formation in HUVEC cells. F9 was negatively correlated with the degree of immune infiltration, and low expression of F9 suggested a poor response to immunotherapy. Conclusion: The role of HCC-related oncogenes and tumor-suppressor genes in diagnosis and prognosis was identified. In addition, we have found that PBK may promote tumor proliferation through angiogenesis and F9 may be a predictor of tumor immunotherapy response.
| S-EPMC9845404 | biostudies-literature