functional enichment with custom gene annotation
3
0
Entering edit mode
5.0 years ago

Hi,

I have a list of gene-annotation pairs and want to do enrichment analysis on specific sets of genes. There are many functional enrichment tools, but few allow me to input my own gene-annotation (each gene belongs to a network) info for the enrichment. What are the tools that can be used in a plug-and-play manner? Preferably, it provides multiple ways of visualisation on the enrichment output.

Thanks.

enrichment • 5.5k views
ADD COMMENT
0
Entering edit mode

You are commenting that you have a list of genes that are networks. How do you define those networks ? And what type of enrichment are you looking for?

ADD REPLY
0
Entering edit mode

They are defined according to their expression profiles measured by RNA-seq. I want to enrich memberships in some gold standard networks so as to understand how my own networks reflect reality.

ADD REPLY
0
Entering edit mode

You keep using the word network but you don't specify in which format do you have this network. How are the edges and nodes of this network stored? Or these custom network is just a list of co-expressed or a selection of genes?

ADD REPLY
0
Entering edit mode
5.0 years ago

You can use the R Bioconductor package topGO with custom annotations which I think you can use with the ViSEAGO package for visualization.

ADD COMMENT
0
Entering edit mode

On the tutorial of topGO (https://bioconductor.org/packages/release/bioc/vignettes/topGO/inst/doc/topGO.pdf). It seems this package is only applicable to GO enrichment. My annotation is about which gene belongs to which network (I did not state clearly my question at the beginning). So topGO is not applicable in this case, right?

ADD REPLY
0
Entering edit mode

Think of each network as a pathway/biological process. You simply annotate each gene with the network(s) it belongs to. If there's no relation between the networks then this doesn't have the complexity of dealing with the GO hierarchy of terms. So if your question is about enrichment in membership in some networks, you just need to build a contingency table for each network and do a chi-squared test and correct for multiple testing. If that's not it, then you need to explain more about what you're trying to do.

ADD REPLY
0
Entering edit mode

Yes, my question is about enrichment in membership in some gold standard networks. I constructed many networks, and each needs such enrichment on gold standard networks, so it should be a programmatic process. I like the powerful visualisations in some GO enrichment tools such as clusterProfiler, so I wonder if some of them can input gene-network annotations for enrichment. If not, I have to build the enrichment and visulisation pipeline by myself.

ADD REPLY
0
Entering edit mode
ADD COMMENT
0
Entering edit mode

There is instruction (https://yulab-smu.github.io/clusterProfiler-book/chapter5.html) of using clusterProfiler on custom annotations, but it only allows custom GO annotations. My custom annotation is predefined networks, e.g. each gene belongs to a network.

ADD REPLY
0
Entering edit mode
4.9 years ago
padwalmk ▴ 140

Hi, If you have got your answer than ok, otherwise you can try Cluego plugin of the Cytoscape. It requires a license but once you register you will get a license. It's pretty easy and helpful for gene-ontology analysis. It can take direct input from your network, and in addition, it has a feature for the gene-ontology term fusion. It allows not only the gene-enrichment but pathway from KEGG, Reactome, and immune system-related. Different layout for visualization, check out particularly cerebral layout which gives you a real feel for cell organization with terms and genes labeled in a hierarchical manner.

ADD COMMENT

Login before adding your answer.

Traffic: 2025 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6