Network analysis of TCGA and GTEx gene expression datasets for identification of trait-associated biomarkers in human cancer.
Ontology highlight
ABSTRACT:
SUBMITTER: Chen HM
PROVIDER: S-EPMC8841814 | biostudies-literature | 2022 Mar
REPOSITORIES: biostudies-literature
Publications
Network analysis of TCGA and GTEx gene expression datasets for identification of trait-associated biomarkers in human cancer.
Chen Huey-Miin HM MacDonald Justin A JA
STAR protocols 20220207 1
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 ...[more]
Similar Datasets
Project description:Reactivation of the multi-subunit ribonucleoprotein telomerase is the primary telomere maintenance mechanism in cancer, but it is rate-limited by the enzymatic component, telomerase reverse transcriptase (TERT). While regulatory in nature, TERT alternative splice variant/isoform regulation and functions are not fully elucidated and are further complicated by their highly diverse expression and nature. Our primary objective was to characterize TERT isoform expression across 7887 neoplastic and 2099 normal tissue samples using The Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression Project (GTEx), respectively. We confirmed the global overexpression and splicing shift towards full-length TERT in neoplastic tissue. Stratifying by tissue type we found uncharacteristic TERT expression in normal brain tissue subtypes. Stratifying by tumor-specific subtypes, we detailed TERT expression differences potentially regulated by subtype-specific molecular characteristics. Focusing on β-deletion splicing regulation, we found the NOVA1 trans-acting factor to mediate alternative splicing in a cancer-dependent manner. Of relevance to future tissue-specific studies, we clustered cancer cell lines with tumors from related origin based on TERT isoform expression patterns. Taken together, our work has reinforced the need for tissue and tumour-specific TERT investigations, provided avenues to do so, and brought to light the current technical limitations of bioinformatic analyses of TERT isoform expression.
| S-EPMC8070023 | biostudies-literature
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:Human disease often arises as a consequence of alterations in a set of associated genes rather than alterations to a set of unassociated individual genes. Most previous microarray-based meta-analyses identified disease-associated genes or biomarkers independent of genetic interactions. Therefore, in this study, we present the first meta-analysis method capable of taking gene combination effects into account to efficiently identify associated biomarkers (ABs) across different microarray platforms.We propose a new meta-analysis approach called MiningABs to mine ABs across different array-based datasets. The similarity between paired probe sequences is quantified as a bridge to connect these datasets together. The ABs can be subsequently identified from an "improved" common logit model (c-LM) by combining several sibling-like LMs in a heuristic genetic algorithm selection process. Our approach is evaluated with two sets of gene expression datasets: i) 4 esophageal squamous cell carcinoma and ii) 3 hepatocellular carcinoma datasets. Based on an unbiased reciprocal test, we demonstrate that each gene in a group of ABs is required to maintain high cancer sample classification accuracy, and we observe that ABs are not limited to genes common to all platforms. Investigating the ABs using Gene Ontology (GO) enrichment, literature survey, and network analyses indicated that our ABs are not only strongly related to cancer development but also highly connected in a diverse network of biological interactions.The proposed meta-analysis method called MiningABs is able to efficiently identify ABs from different independently performed array-based datasets, and we show its validity in cancer biology via GO enrichment, literature survey and network analyses. We postulate that the ABs may facilitate novel target and drug discovery, leading to improved clinical treatment. Java source code, tutorial, example and related materials are available at "http://sourceforge.net/projects/miningabs/".
| S-EPMC4068973 | biostudies-literature
Project description:Several large-scale human cancer genomics projects such as TCGA offered huge genomic and clinical data for researchers to obtain meaningful genomics alterations which intervene in the development and metastasis of the tumor. A web-based TCGA data analysis platform called TCGA4U was developed in this study. TCGA4U provides a visualization solution for this study to illustrate the relationship of these genomics alternations with clinical data. A whole genome screening of the survival related gene expression patterns in breast cancer was studied. The gene list that impacts the breast cancer patient survival was divided into two patterns. Gene list of each of these patterns was separately analyzed on DAVID. The result showed that mitochondrial ribosomes play a more crucial role in the cancer development. We also reported that breast cancer patients with low HSPA2 expression level had shorter overall survival time. This is widely different to findings of HSPA2 expression pattern in other cancer types. TCGA4U provided a new perspective for the TCGA datasets. We believe it can inspire more biomedical researchers to study and explain the genomic alterations in cancer development and discover more targeted therapies to help more cancer patients.
| S-EPMC4630377 | biostudies-literature
Project description:The mechanisms underlying the oncogenesis and progression of lung adenocarcinoma (LUAD) are currently unclear. The discovery of competitive endogenous RNA (ceRNA) regulatory networks has provided a new direction for the treatment and prognosis of patients with LUAD. However, the mechanism of action of ceRNA in LUAD remains elusive. In the present study, differentially expressed mRNAs, microRNAs (miRs) and long non-coding RNAs from the cancer genome atlas database were screened. CeRNAs for LUAD were then identified using online prediction software. Among the ceRNAs identified, family with sequence similarity 83 member A (FAM83A), miR-34c-5p, KCNQ1OT1 and FLJ26245 were observed to be significantly associated with the overall survival of patients with LUAD. Of note, FAM83A has potential significance in drug resistance, and may present a candidate biomarker for the prognosis and treatment of patients with LUAD.
| S-EPMC6485208 | biostudies-literature
Project description:DNA methylation (DNAm) plays an important role in the pathogenesis of psoriasis through regulating mRNA expressions. This study aimed to identify hub genes regulated by DNAm as biomarkers of psoriasis. Psoriatic skin tissues gene expression and methylation datasets were downloaded from Gene Expression Omnibus (GEO) database. Subsequently, multiple computational approaches, including immune infiltration analysis, enrichment analysis, protein-protein interaction (PPI) network establishment, and machine learning algorithm analysis (lasso, random forest, and SVM-RFE), were performed to analyze the regulatory networks, to recognize hub genes, and to clarify the pathogenesis of psoriasis. Finally, the hypermethylated genes were used to immune cell infiltration analysis, which revealed that psoriasis skin tissues were mainly composed of activated dendritic cells, resting mast cells, T follicular helper cells (cTfh), etc. Differentially expressed-methylated genes (DEMGs) were identified and partitioned into four subgroups and the 97 significantly hypermethylated and downregulated (hyper-down) genes accounted for the highest proportion (47%). Hyper-down genes were mainly enriched in glucose homeostasis, AMP-activated protein kinase (AMPK) signaling pathway, lipid storage disease, partial lipodystrophy, and insulin resistance. Furthermore, insulin receptor substrate 1 (IRS1), Rho guanine nucleotide exchange factor 10 (ARHGEF10) and retinoic acid induced 14 (RAI14) were identified as potential targets. These findings provided new ideas for future studies of psoriasis on the occurrence and the molecular mechanisms.
| S-EPMC8427602 | 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:Endocervical adenocarcinoma (EAC) is an aggressive type of endocervical cancer. At present, molecular research on EAC mainly focuses on the genome and mRNA transcriptome, the investigation of small RNAs in EAC has not been fully described. Here, we systematically explored small RNAs in 14 EAC patients with different subtypes using small RNA sequencing. MiRNAs and tRNA-derived RNAs (tDRs) accounted for the majority of mapped reads and the total number of miRNAs and tDRs maintained a relative balance. To explore the correlations between small RNAs expression and EAC with different clinical characteristics, we performed the weighted gene co-expression network analysis (WGCNA) and screened for hub small RNAs. From the key modules, we identified 9 small RNAs that were significantly related to clinical characteristics in EAC patients. Gene ontology and pathway analyses revealed that these molecules were involved in the pathogenesis of EAC. Our work provided new insights into EAC pathogenesis and successfully identified several small RNAs as candidate biomarkers for diagnosis and prognosis of EAC.
| S-EPMC8507262 | 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:BackgroundMuscular invasive bladder cancer (MIBC) is a common malignant tumor in the world. Because of their heterogeneity in prognosis and response to treatment, biomarkers that can predict survival or help make treatment decisions in patients with MIBC are essential for individualized treatment.AimWe aimed to integrate bioinformatics research methods to identify a set of effective biomarkers capable of predicting, diagnosing, and treating MIBC. To provide a new theoretical basis for the diagnosis and treatment of bladder cancer.Methods and resultsGene expression profiles and clinical data of MIBC were obtained by downloading from the Cancer Genome Atlas database. A dataset of 129 MIBC cases and controls was included. 2084 up-regulated genes and 2961 down-regulated genes were identified by differentially expressed gene (DEG) analysis. Then, gene ontology analysis was performed to explore the biological functions of DEGs, respectively. The up-regulated DEGs are mainly enriched in epidermal cell differentiation, mitotic nuclear division, and so forth. They are also involved in the cell cycle, p53 signaling pathway, PPAR signaling pathway, and so forth. The weighted gene co-expression network analysis yielded five modules related to pathological stages and grading, of which blue and turquoise were the most relevant modules for MIBC. Next, Using Kaplan-Meier survival analysis to identify further hub genes, the screening criteria at p ≤ .05, we found CNKSR1, HIP1R, CFL2, TPM1, CSRP1, SYNM, POPDC2, PJA2, and RBBP8NL genes associated with the progression and prognosis of MIBC patients. Finally, immunohistochemistry experiments further confirmed that CNKSR1 plays a vital role in the tumorigenic context of MIBC.ConclusionThe research suggests that CNKSR1, POPDC2, and PJA2 may be novel biomarkers as therapeutic targets for MIBC, especially we used immunohistochemical further to validate CNKSR1 as a therapeutic target for MIBC which may help to improve the prognosis for MIBC.
| S-EPMC9458504 | biostudies-literature