This is a test version of Biostars. For the public version, visit https://www.biostars.org.
gene set enrichment analysis

hi all

i have one doubt regarding the gene set enrichment analysis method i have a list of upregulated genes and i have not done the differential gene expression analysis but instead i did rank prod analysis from which i obtained the list of upregulated genes so now when i m searching about the ways to analyse them i m not understanding as most of them are analysing the result from deseq2 and i m looking for the way to analyse my gene list so is there any tutorial or any material which works only on gene list

thanks

gsea pathwayanalysis
• 4,259 views • link
written by nishimalhotra2612 • nishimalhotra2612

4 answers

So basically you have a list of genes and now want to know whether they're enriched for some terms or processed? I like g:profiler2 for these kinds of analysis. It has an R package if you want to automate the task.

• 0 views • link
written by ATpoint ATpoint

thanks for suggesting i tried to use that but the thing is that i have my gene list in form of csv the name is "upregulated.csv"and it has 720 genes so when i try to input my files and run the code it says no result to show. i tried to found out the format but they only mentioned the types of input not the format

so in what format i can input the data if you have any idea

thanks

• 0 views • link
written by nishimalhotra2612 • nishimalhotra2612

Assuming the formatting of the genes was correct (please show head upregulated.csv) it can simply mean that there are no significant enrichments.

• 0 views • link
written by ATpoint ATpoint

Yeah, if you indeed feed these IDs into the tool and it returns nothing then there is simply no statistical enrichment. Try relaxing the FDR threshold, maybe that helps. You can do that with the R package.

• 0 views • link
written by ATpoint ATpoint

In your case, I would suggest performing an over-representation analysis (ORA) rather than GSEA. One quick way to perform ORA is on WebGestalt.

If you want to perform a proper GSEA, you need to have a full count matrix with normalized read counts per each gene as rows and all samples as columns. Then you can run GSEA, as described here (https://www.gsea-msigdb.org/gsea/index.jsp).

• 0 views • link
written by Marco Pannone • Marco Pannone

In addition to ATpoint's suggestion of g:Profiler, you can try PANTHER, agriGO, and DAVID. Each of these tools takes a list of genes and returns gene set enrichment analysis results.

PANTHER

agriGO

DAVID

• 0 views • link
written by Jeremy • Jeremy

hello all i tried to use gprofiler2 and i got the result i m sharing my way on what i did to input multiple queries so that it can be helpful to others

so as it is being said in the gprofiler tutorial that we have to create a named list of vectors so i did that its very basic code but it can be helpful to the people who are new to bioinformatics

first input the csv file df<-read.csv("upregulated genes.CSV")

select the column consisting genes or gene ids df1<-df$x

converting the selected column to vector vectordata<-as.vector(df1)

then creating the list of that vectors gene_list<-list(vectordata)

and now input the gene_list into gprofiler

hope it helps!

• 0 views • link
written by nishimalhotra2612 • nishimalhotra2612

Log in to answer this question.

More posts like this