scholarly journals OneStopRNAseq: A Web Application for Comprehensive and Efficient Analyses of RNA-Seq Data

Genes ◽  
2020 ◽  
Vol 11 (10) ◽  
pp. 1165 ◽  
Author(s):  
Rui Li ◽  
Kai Hu ◽  
Haibo Liu ◽  
Michael R. Green ◽  
Lihua Julie Zhu
Keyword(s):  
Gene Expression ◽  
Web Application ◽  
Enrichment Analysis ◽  
Gene Expression Omnibus ◽  
Gene Set Enrichment Analysis ◽  
Specific Gene ◽  
Differential Analysis ◽  
Rna Seq ◽  
Data Quality Control ◽  
Private And Public

Over the past decade, a large amount of RNA sequencing (RNA-seq) data were deposited in public repositories, and more are being produced at an unprecedented rate. However, there are few open source tools with point-and-click interfaces that are versatile and offer streamlined comprehensive analysis of RNA-seq datasets. To maximize the capitalization of these vast public resources and facilitate the analysis of RNA-seq data by biologists, we developed a web application called OneStopRNAseq for the one-stop analysis of RNA-seq data. OneStopRNAseq has user-friendly interfaces and offers workflows for common types of RNA-seq data analyses, such as comprehensive data-quality control, differential analysis of gene expression, exon usage, alternative splicing, transposable element expression, allele-specific gene expression quantification, and gene set enrichment analysis. Users only need to select the desired analyses and genome build, and provide a Gene Expression Omnibus (GEO) accession number or Dropbox links to sequence files, alignment files, gene-expression-count tables, or rank files with the corresponding metadata. Our pipeline facilitates the comprehensive and efficient analysis of private and public RNA-seq data.

scholarly journals Bioinformatics-based Screening of Key Genes for Saponin Metabolism in Quinoa

2021 ◽  
Author(s):  
Chengang Guo ◽  
Zhimin wei ◽  
Wei Lyu ◽  
Yanlou Geng
Keyword(s):  
Differentially Expressed Genes ◽  
Principal Component ◽  
Enrichment Analysis ◽  
Hierarchical Cluster ◽  
Gene Set Enrichment Analysis ◽  
Differentially Expressed ◽  
Differential Analysis ◽  
Rna Seq ◽  
Protein Coding ◽  
Key Genes

Abstract Quinoa saponins have complex, diverse and evident physiologic activities. However, the key regulatory genes for quinoa saponin metabolism are not yet well studied. The purpose of this study was to explore genes closely related to quinoa saponin metabolism. In this study, the significantly differentially expressed genes in yellow quinoa were firstly screened based on RNA-seq technology. Then, the key genes for saponin metabolism were selected by gene set enrichment analysis (GSEA) and principal component analysis (PCA) statistical methods. Finally, the specificity of the key genes was verified by hierarchical clustering. The results of differential analysis showed that 1654 differentially expressed genes were achieved after pseudogenes deletion. Therein, there were 142 long non-coding genes and 1512 protein-coding genes. Based on GSEA analysis, 116 key candidate genes were found to be significantly correlated with quinoa saponin metabolism. Through PCA dimension reduction analysis, 57 key genes were finally obtained. Hierarchical cluster analysis further demonstrated that these key genes can clearly separate the four groups of samples. The present results could provide references for the breeding of sweet quinoa and would be helpful for the rational utilization of quinoa saponins.

scholarly journals Bioinformatic Identification of Hub Genes and Biological Pathways for Sepsis-Associated Acute Lung Injury

2021 ◽  
Author(s):  
Chao Zhang ◽  
Feng Xu ◽  
Fang Fang
Keyword(s):  
Gene Expression ◽  
Acute Lung Injury ◽  
Lung Injury ◽  
Function Analysis ◽  
Enrichment Analysis ◽  
Gene Expression Omnibus ◽  
Biological Pathways ◽  
Gene Set Enrichment Analysis ◽  
Ppi Network ◽  
Hub Genes

Abstract Background: Sepsis-associated acute lung injury (ALI) is a potentially lethal complication associated with a poor prognosis and high mortality worldwide, especially in the outbreak of COVID-19. However, the fundamental mechanisms of this complication were still not fully elucidated. Thus, we conducted this study to identify hub genes and biological pathways of sepsis-associated ALI, mainly focus on two pathways of LPS and HMGB1. Methods: Gene expression profile GSE3037 were downloaded from Gene Expression Omnibus (GEO) database, including 8 patients with sepsis-induced acute lung injury, with 8 unstimulated blood neutrophils, 8 LPS- induced neutrophils and 8 HMGB1-induced neutrophils. Differentially expressed genes (DEGs) identifications, Gene Ontology (GO) function analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis, Gene Set Enrichment Analysis (GSEA) and protein-protein interaction (PPI) network constructions were performed to obtain hub genes and relevant biological pathways.Results: We identified 534 and 317 DEGs for LPS- and HMGB1-induced ALI, respectively. The biological pathways involved in LPS- and HMGB1-induced ALI were also identified accordingly. By PPI network analysis, we found that ten hub genes for LPS-induced ALI (CXCL8, TNF, IL6, IL1B, ICAM1, CXCL1, CXCL2, IL1A, IL1RN and CXCL3) and another ten hub genes for HMGB1-induced ALI (CCL20, CXCL2, CXCL1, CCL4, CXCL3, CXCL9, CCL21, CXCR6, KNG1 and SST). Furthermore, by combining analysis, the results revealed that genes of TNF, CCL20, IL1B, NFKBIA, CCL4, PTGS2, TNFAIP3, CXCL2, CXCL1 and CXCL3 were potential biomarkers for sepsis-associated ALI. Conclusions: Our study revealed that ten hub genes associated with sepsis-induced ALI were TNF, CCL20, IL1B, NFKBIA, CCL4, PTGS2, TNFAIP3, CXCL2, CXCL1 and CXCL3, which may serve as genetic biomarkers and be further verified in prospective experimental trials.

scholarly journals Integrating Literature-Based Knowledge Database and Expression Data to Explore Molecular Pathways Connecting PPARG and Myocardial Infarction

PPAR Research ◽  
2020 ◽  
Vol 2020 ◽  
pp. 1-6
Author(s):  
Rongyuan Cao ◽  
Yan Dong ◽  
Kamil Can Kural
Keyword(s):  
Gene Expression ◽  
Myocardial Infarction ◽  
Enrichment Analysis ◽  
Protective Role ◽  
Gene Expression Omnibus ◽  
Gene Set Enrichment Analysis ◽  
Peroxisome Proliferator ◽  
Peroxisome Proliferator Activated Receptor ◽  
Knowledge Database ◽  
Expression Activity

Peroxisome proliferator-activated receptor γ (PPARG) might play a protective role in the development of myocardial infarction (MI) with limited mechanisms identified. Genes associated with both PPARG and MI were extracted from Elsevier Pathway Studio to construct the initial network. The gene expression activity within the network was estimated through a mega-analysis with eight independent expression datasets derived from Gene Expression Omnibus (GEO) to build PPARG and MI connecting pathways. After that, gene set enrichment analysis (GSEA) was conducted to explore the functional profile of the genes involved in the PPARG-driven network. PPARG demonstrated a significantly low expression in MI patients (LFC=−0.52; p<1.84e−9). Consequently, PPARG could indicatively be promoting three MI inhibitors (e.g., SOD1, CAV1, and POU5F1) and three MI-downregulated markers (e.g., ALB, ACADM, and ADIPOR2), which were deactivated in MI cases (p<0.05), and inhibit two MI-upregulated markers (RELA and MYD88), which showed increased expression levels in MI cases (p=0.0077 and 0.047, respectively). These eight genes were mainly enriched in nutrient- and cell metabolic-related pathways and functionally linked by GSEA and PPCN. Our results suggest that PPARG could protect the heart against both the development and progress of MI through the regulation of nutrient- and metabolic-related pathways.

scholarly journals Overexpression of SP6 Correlates with Osteosarcoma Metastasis and Poor Prognosis

2020 ◽  
pp. 1-7
Author(s):  
Dongmei Guo ◽  
Chunpu Li ◽  
Sicheng Wang ◽  
Lili Zhao ◽  
Dongmei Guo ◽  
...  
Keyword(s):  
Gene Expression ◽  
Cox Regression ◽  
Excision Repair ◽  
Enrichment Analysis ◽  
Analysis Data ◽  
Gene Expression Omnibus ◽  
Gene Set Enrichment Analysis ◽  
Roc Curve Analysis ◽  
Gene Set ◽  
Base Excision

Background: SP6 (Specificity protein 6) has been explored as a prospective biomarker in several cancers. In this research, the prognostic value of SP6 expression in osteosarcoma was predicted by bioinformatics analysis. Data were obtained from the Gene Expression Omnibus (GEO) database. Methods: Gene expression data and clinical materials were downloaded from the GSE21257 dataset. The mRNA expression of SP6 was compared between metastatic and non-metastatic tissues with the Wilcoxon rank-sum test, and the relationship between SP6 and clinicopathological characters was analysed using logistic regression. In addition, the correlation between SP6 and survival rate was assessed using KaplanMeier and Cox regression. Moreover, receiver operating characteristic (ROC) curve analysis was conducted to determine the prognostic merit of SP6 for osteosarcoma. The biological functions of SP6 were annotated and evaluated through gene set enrichment analysis (GSEA) and gene set variation analysis (GSVA). Results: SP6 was significantly highly expressed in metastatic osteosarcoma tissues (p = 0.002). High SP6 expression showed a positive correlation with Huvos grade (OR = 6.60 for I vs. II, p = 0.028). The overall survival (OS) of the patients with high SP6 expression was significantly poorer than the low SP6 expression group (p = 0.027). The multivariate analysis revealed that SP6 expression (p = 0.002, HR = 15.40 (95% CI [2.84–83.44])) was independently correlated with OS. GSEA and GSVA showed that "spliceosome" and "base excision repair" were significantly upregulated in the high expression group of SP6. Conclusion: SP6 may serve an independent prognostic biomarker in osteosarcoma.

Identification of MEG8-miR378d-SOBP Axis as a Novel Regulatory Network and Associated With Immune Infiltrates in Ovarian Carcinoma by Integrated Bioinformatics Analysis

2020 ◽  
Author(s):  
Jian Lei ◽  
Zhen-Yu He ◽  
Jun Wang ◽  
Min Hu ◽  
Ping Zhou ◽  
...  
Keyword(s):  
Gene Expression ◽  
Ovarian Cancer ◽  
Bioinformatics Analysis ◽  
Enrichment Analysis ◽  
Functional Enrichment Analysis ◽  
Gene Expression Omnibus ◽  
Gene Set Enrichment Analysis ◽  
The Cancer Genome Atlas ◽  
Functional Enrichment ◽  
Human Protein Atlas

Abstract BackgroundTo investigate the potential molecular mechanism of ovarian cancer (OC) evolution and immunological correlation using the integrated bioinformatics analysis.MethodsData from the Gene Expression Omnibus (GEO) was used to gain differentially expressed genes (DEGs). Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genome (KEGG) pathway analysis were completed by utilizing the Database for Annotation, Visualization, and Integrated Discovery (DAVID). After multiple validation via The Cancer Genome Atlas (TCGA), Gene Expression Profiling Interactive Analysis 2 (GEPIA 2), the Human Protein Atlas (HPA) and Kaplan-Meier (KM) plotter, immune logical relationships of the key gene SOBP were evaluated based on Tumor Immune Estimation Resource (TIMER), and Gene Set Enrichment Analysis (GSEA) software. Finally, the lncRNAs-miRNAs-mRNAs sub-network was predicted by starBase, Targetscan, miRBD, and LncBase, individually. Correlation of expression and prognosis for mRNAs, miRNAs and lncRNAs were confirmed by TCGA, GEPIA 2, starBase, and KM.ResultsA total of 192 shared DEGs were discovered from the four data sets, including 125 upregulated and 67 downregulated genes. Functional enrichment analysis presented that they were mainly enriched in cartilage development, pathway in PI3K-Akt signaling pathway. Lower expression of SOBP was the independent prognostic factor for inferior prognosis in OC patients. Intriguingly, downregulated SOBP enhanced the infiltration levels of B cells, CD8+ T cells, Macrophage, Neutrophil and Dendritic cells. GSEA also disclosed low SOBP showed significantly association with the activation of various immune-related pathways. Finally, we firstly reported that MEG8-miR378d-SOBP axis was linked to development and prognosis of ovarian cancer through regulating cytokines pathway.Conclusions Our study establishes a novel MEG8-miR378d-SOBP axis in the development and prognosis of OC, and the triple sub-network probably affects the progression of ovarian tumor by regulating cytokines pathway.

scholarly journals Identification of Candidate Biomarkers Correlated Poorer Prognosis of Breast Cancer via Bioinformatics Method

2021 ◽  
Author(s):  
Gang Chen ◽  
Mingwei Yu ◽  
Jianqiao Cao ◽  
Huishan Zhao ◽  
Yuanping Dai ◽  
...  
Keyword(s):  
Breast Cancer ◽  
Gene Expression ◽  
Kegg Pathway ◽  
Expression Profiles ◽  
Gene Expression Profiles ◽  
Enrichment Analysis ◽  
Gene Expression Omnibus ◽  
Gene Set Enrichment Analysis ◽  
Molecular Therapy ◽  
Hub Genes

Abstract Background: Breast cancer (BC) is a malignancy with a high incidence among women in the world, and it is very urgent to identify significant biomarkers and molecular therapy methods.Methods: Total 58 normal tissues and 203 cancer tissues were collected from three Gene Expression Omnibus (GEO) gene expression profiles, and the differential expressed genes (DEGs) were identified. Subsequently, the Gene Ontology (GO) function and Kyoto Encyclopedia of Genes and Genome (KEGG) pathway were analyzed. Additionally, hub genes were screened by constructing a protein-protein interaction (PPI) network. Then, we explored the prognostic values and molecular mechanism of these hub genes Kaplan-Meier (KM) curve and Gene Set Enrichment Analysis (GSEA). Results: 42 up-regulated and 82 down-regulated DEGs were screened out from GEO datasets. GO and KEGG pathway analysis revealed that DEGs were mainly related to cell cycles and cell proliferation. Furthermore, 12 hub genes (FN1, AURKA, CCNB1, BUB1B, PRC1, TPX2, NUSAP1, TOP2A, KIF20A, KIF2C, RRM2, ASPM) with a high degree of genes were selected, among which, 11 hub gene were significantly correlated with the prognosis of patients with BC. From GSEA reviewed correlated with KEGG_CELL_CYCLE and HALLMARK_P53_PATHWAY. Conclusion: this study identified 11 key genes as BC potential prognosis biomarkers on the basis of integrated bioinformatics analysis. This finding will improve our knowledge of the BC progress and mechanisms.

scholarly journals Identification of Potential Crucial Genes in Atrial Fibrillation: A Bioinformatic Analysis

2020 ◽  
Author(s):  
Junguo Zhang ◽  
Xin Huang ◽  
Xiaojie Wang ◽  
Yanhui Gao ◽  
Li Liu ◽  
...  
Keyword(s):  
Gene Expression ◽  
Atrial Fibrillation ◽  
Growth Factor ◽  
Microarray Data ◽  
Enrichment Analysis ◽  
Bioinformatic Analysis ◽  
Gene Expression Omnibus ◽  
Gene Set Enrichment Analysis ◽  
Functional Enrichment ◽  
Insulin Like Growth Factor

Abstract Background Atrial fibrillation (AF) is at least partially heritable, affecting 2-3% of the population in Europe and the USA. However, a substantial proportion of heritability is still lacking. In the present study, we aim to identify potential crucial genes associated with AF through bioinformatic analyses of public datasets. Methods Microarray data sets of GSE115574, GSE31821, GSE79768, GSE41177 and GSE14975 from the Gene Expression Omnibus (GEO) database were retrieved. After merging all microarray data and adjusting batch effect, differentially expressed genes (DEGs) were identified. Functional enrichment analyses based on Gene Ontology (GO) resource, Kyoto Encyclopedia of Genes and Genomes (KEGG) resource, Gene Set Enrichment Analysis (GSEA), Reactome Pathway Database and Disease Ontology (DO) were carried out. Protein-protein interaction (PPI) network was constructed using the STRING database. Combined with aforementioned significant bioinformatics information, potential crucial genes were subsequently selected. The comparative toxicogenomics database (CTD) was carried out to explore the interaction between potential crucial genes and AF. Result We identified 27 of DEGs with gene expression fold change (FC) ≥ 1.5 or ≤ 2/3 (|log2 FC| ≥ 0.58) and 5 with FC ≥ 2 or ≤ 0.5 (|log2 FC| ≥ 1) of AF patients compared with sinus rhythm controls. The most significantly enriched pathway was regulation of insulin-like growth factor transport and uptake by insulin-like growth factor binding proteins. IGFBP2, C1orf105, FHL2, CHGB, ATP1B4, IGFBP3, SLC26A9, CXCR4 and HTR2B were considered the potential crucial genes. CTD showed CXCR4, IGFBP2, IGFBP3 and FHL2 had higher scores with AF. Conclusions The 9 potential crucial genes, especially CXCR4, IGFBP2, IGFBP3 and FHL2 , may be associated with risk of AF. Our study provided new insights of AF into genetics, molecular pathogenesis and new therapeutic targets.

scholarly journals Transcriptome-Wide Gene Expression in a Murine Model of Ventilator-Induced Lung Injury

2021 ◽  
Vol 2021 ◽  
pp. 1-23
Author(s):  
Zhao Li ◽  
Guoshao Zhu ◽  
Chen Zhou ◽  
Hui Wang ◽  
Le Yu ◽  
...  
Keyword(s):  
Gene Expression ◽  
Lung Injury ◽  
Enrichment Analysis ◽  
Gene Expression Omnibus ◽  
Gene Set Enrichment Analysis ◽  
Hub Genes ◽  
Qrt Pcr ◽  
Ventilator Induced Lung Injury ◽  
Significant Difference ◽  
Tf Gene

Background. Mechanical ventilation could lead to ventilator-induced lung injury (VILI), but its underlying pathogenesis remains largely unknown. In this study, we aimed to determine the genes which were highly correlated with VILI as well as their expressions and interactions by analyzing the differentially expressed genes (DEGs) between the VILI samples and controls. Methods. GSE11434 was downloaded from the gene expression omnibus (GEO) database, and DEGs were identified with GEO2R. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were conducted using DAVID. Next, we used the STRING tool to construct protein-protein interaction (PPI) network of the DEGs. Then, the hub genes and related modules were identified with the Cytoscape plugins: cytoHubba and MCODE. qRT-PCR was further used to validate the results in the GSE11434 dataset. We also applied gene set enrichment analysis (GSEA) to discern the gene sets that had a significant difference between the VILI group and the control. Hub genes were also subjected to analyses by CyTargetLinker and NetworkAnalyst to predict associated miRNAs and transcription factors (TFs). Besides, we used CIBERSORT to detect the contributions of different types of immune cells in lung tissues of mice in the VILI group. By using DrugBank, small molecular compounds that could potentially interact with hub genes were identified. Results. A total of 141 DEGs between the VILI group and the control were identified in GSE11434. Then, seven hub genes were identified and were validated by using qRT-PCR. Those seven hub genes were largely enriched in TLR and JAK-STAT signaling pathways. GSEA showed that VILI-associated genes were also enriched in NOD, antigen presentation, and chemokine pathways. We predicted the miRNAs and TFs associated with hub genes and constructed miRNA-TF-gene regulatory network. An analysis with CIBERSORT showed that the proportion of M0 macrophages and activated mast cells was higher in the VILI group than in the control. Small molecules, like nadroparin and siltuximab, could act as potential drugs for VILI. Conclusion. In sum, a number of hub genes associated with VILI were identified and could provide novel insights into the pathogenesis of VILI and potential targets for its treatment.

scholarly journals Identification and Validation of a Novel Immune-Related Four-lncRNA Signature for Lung Adenocarcinoma

2021 ◽  
Vol 12 ◽  
Author(s):  
Jixin Wang ◽  
Xiangjun Yin ◽  
Yin-Qiang Zhang ◽  
Xuming Ji
Keyword(s):  
Gene Expression ◽  
Lung Adenocarcinoma ◽  
Immune Cell ◽  
Cox Regression ◽  
Enrichment Analysis ◽  
Gene Expression Omnibus ◽  
Gene Set Enrichment Analysis ◽  
The Cancer Genome Atlas ◽  
Lasso Regression ◽  
Kaplan Meier

Lung adenocarcinoma (LUAD) is a major subtype of lung cancer, the prognosis of patients with which is associated with both lncRNAs and cancer immunity. In this study, we collected gene expression data of 585 LUAD patients from The Cancer Genome Atlas (TCGA) database and 605 subjects from the Gene Expression Omnibus (GEO) database. LUAD patients were divided into high and low immune-cell-infiltrated groups according to the single sample gene set enrichment analysis (ssGSEA) algorithm to identify differentially expressed genes (DEGs). Based on the 49 immune-related DE lncRNAs, a four-lncRNA prognostic signature was constructed by applying least absolute shrinkage and selection operator (LASSO) regression, univariate Cox regression, and stepwise multivariate Cox regression in sequence. Kaplan–Meier curve, ROC analysis, and the testing GEO datasets verified the effectiveness of the signature in predicting overall survival (OS). Univariate Cox regression and multivariate Cox regression suggested that the signature was an independent prognostic factor. The correlation analysis revealed that the infiltration immune cell subtypes were related to these lncRNAs.

scholarly journals Identification of Potential Key Genes and Pathways in Enzalutamide-Resistant Prostate Cancer Cell Lines: A Bioinformatics Analysis with Data from the Gene Expression Omnibus (GEO) Database

2020 ◽  
Vol 2020 ◽  
pp. 1-12
Author(s):  
Long Zheng ◽  
Xiaojie Dou ◽  
Xiaodong Ma ◽  
Wei Qu ◽  
Xiaoshuang Tang
Keyword(s):  
Gene Expression ◽  
Prostate Cancer ◽  
Cell Lines ◽  
Molecular Mechanisms ◽  
Enrichment Analysis ◽  
Gene Expression Omnibus ◽  
Gene Set Enrichment Analysis ◽  
Pathway Enrichment Analysis ◽  
Key Genes ◽  
Geo Database

Enzalutamide (ENZ) has been approved for the treatment of advanced prostate cancer (PCa), but some patients develop ENZ resistance initially or after long-term administration. Although a few key genes have been discovered by previous efforts, the complete mechanisms of ENZ resistance remain unsolved. To further identify more potential key genes and pathways in the development of ENZ resistance, we employed the GSE104935 dataset, including 5 ENZ-resistant (ENZ-R) and 5 ENZ-sensitive (ENZ-S) PCa cell lines, from the Gene Expression Omnibus (GEO) database. Integrated bioinformatics analyses were conducted, such as analysis of differentially expressed genes (DEGs), Gene Ontology (GO) enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, protein-protein interaction (PPI) analysis, gene set enrichment analysis (GSEA), and survival analysis. From these, we identified 201 DEGs (93 upregulated and 108 downregulated) and 12 hub genes (AR, ACKR3, GPER1, CCR7, NMU, NDRG1, FKBP5, NKX3-1, GAL, LPAR3, F2RL1, and PTGFR) that are potentially associated with ENZ resistance. One upregulated pathway (hedgehog pathway) and seven downregulated pathways (pathways related to androgen response, p53, estrogen response, TNF-α, TGF-β, complement, and pancreas β cells) were identified as potential key pathways involved in the occurrence of ENZ resistance. Our findings may contribute to further understanding the molecular mechanisms of ENZ resistance and provide some clues for the prevention and treatment of ENZ resistance.

Sign in / Sign up

Export Citation Format

Share Document