In addition to Weibo, there is also WeChat
Please pay attention
WeChat public account
Shulou
2025-04-05 Update From: SLTechnology News&Howtos shulou NAV: SLTechnology News&Howtos > Development >
Share
Shulou(Shulou.com)06/01 Report--
In this article, the editor introduces in detail "how to use enricher in clusterProfiler for enrichment analysis". The content is detailed, the steps are clear, and the details are handled properly. I hope that this article "how to use enricher in clusterProfiler for enrichment analysis" can help you solve your doubts.
1. problem
GO or kegg enrichment analysis of genes requires annotation information from the GO and KEGG databases of all genes in the genome. It is very simple for people who are model species to have ready-made annotation results, which can be used directly. For example, clusterProfiler can be directly used for gene set enrichment analysis. However, for the research objects of non-model animals and plants, there are often no ready-made annotation results, so it is impossible to carry out enrichment analysis directly.
Therefore, it is very important to learn the batch annotation of gene function; here is a method for GO and KEGG enrichment analysis of all non-model species.
two。 Batch annotation of gene
EggNOG is used to annotate the genome. in order to carry out enrichment analysis, we must first have GO annotation and KEGG annotation of the background data set, and eggNOG is selected for annotation here.
It is an online server, which can be annotated by uploading with a click of the mouse without complicated configuration. Although eggNOG is a web server, it can annotate up to 100000 sequences at a time, which should fully meet the requirements. Submit the corresponding cds sequence or protein sequence of your gene to this website for annotation.
Documents to be prepared for 3.GO enrichment analysis:
To study the corresponding GO files of all genes in the species genome:
Go2gene.tsv: collated from the eggNOG comment result file
GOGENE
CLASSGO:0000165Pg_S3686.2
Biological_processGO:0003674Pg_S3686.2
Molecular_function.
Function description file corresponding to go2name.tsv:GO term
First of all, you need to go to GO to download the obo file of GO. Here I use go-basic.obo, and then I write a script that parses the obo file into the following format:
Http://purl.obolibrary.org/obo/go/go-basic.obo
GODESCCLASSGO:0000001mitochondrion inheritancebiological_processGO:0000007low-affinity zinc ion transmembrane transporter activitymolecular_function
4. Documents needed for KEGG enrichment analysis
Ko2gene.tsv: collated from the eggNOG comment result file
KOGENEko00920Pg_S3686.2
Ko01100Pg_S33386.2
The name of the ko2name.tsv ko path
KODESCko00440Phosphonate and phosphinate metabolismko00450Selenocompound metabolismko00460Cyanoamino acid metabolismko00471D-Glutamine and D-glutamate metabolismko00472D-Arginine and D-ornithine metabolismko00473D-Alanine metabolismko00480Glutathione metabolismko00510N-Glycan biosynthesisko00513Various types of N-glycan biosynthesisko00512Mucin type O-glycan biosynthesis5. Enrichment analysis
Enrichment analysis is carried out by using the general function enricher in clusterProfiler:
Library (clusterProfiler) ko2name
Welcome to subscribe "Shulou Technology Information " to get latest news, interesting things and hot topics in the IT industry, and controls the hottest and latest Internet news, technology news and IT industry trends.
Views: 0
*The comments in the above article only represent the author's personal views and do not represent the views and positions of this website. If you have more insights, please feel free to contribute and share.
Continue with the installation of the previous hadoop.First, install zookooper1. Decompress zookoope
"Every 5-10 years, there's a rare product, a really special, very unusual product that's the most un
© 2024 shulou.com SLNews company. All rights reserved.