This should be the simplest and most useful GSEA enrichment analysis tool in history

What is GSEA? Baidu is responsible for this. You should know what GSEA is before using it Then upload the tool screenshot: Generally, we need to prepare three documents for GSEA 1. Expression spectrum of gct format 2. Group Files 3. Rich background files, that is, gmt format

What is GSEA? Baidu is responsible for this. You should know what GSEA is before using it

Tools download link: http://gap.shengxin.ren/tool/4/

Then upload the tool screenshot:

 attachments-2017-09-dCFoPOwo59c07814542c Generally, we need to prepare three documents for GSEA

1. Expression spectrum of gct format

2. Group Files

3. Enriched background files, i.e. files in gmt format

With these three files, GSEA enrichment analysis can be done

However, this is extremely inconvenient for many cases. For example, if you want to see the difference in the pathway between a gene with high expression and a gene with low expression in an expression profile, you need to program your own manual grouping before doing GSEA. If one gene is OK, how about multiple genes

This tool focuses on solving these complex problems

It can be seen from the software interface that the input file only has the expression spectrum and the sample information table, but the sample information table is not necessary. If you want to classify the samples according to the sample information and then do the GSEA, you can import the sample information table. Note here that the software will automatically match the sample information, without one-to-one correspondence, which is very user-friendly.

Let's play with an example,

Use the GSE44076 data set. The data download and pre-processing tutorial is shown here: https://www.shengxin.ren/article/135

Now let's use the processed expression spectrum matrix and sample information, and import the expression spectrum matrix first

 attachments-2017-09-3w1vTzQk59c07a1fac27 Then let's look at the KEGG pathway enrichment differences between the top 25% and the bottom 25% samples with the highest expression of ERBB2, as shown in the figure:

 attachments-2017-09-esZOKoD059c07a9d8c35 Click Batch Run, select the directory to save the results, and wait a little while. Finally, three files are generated in the directory to save the results. gct and cls can be used to repeat the results themselves

 attachments-2017-09-03ehiIks59c07aec2340 The folder at the beginning of my_analysis. Gsea. stores the results of GSEA enrichment analysis. Take a look

 attachments-2017-09-oHVb4MZs59c07b6da233 Of course, the software provides six kinds of gmt data by default. If you need to customize it yourself, it is also very simple. You can only copy the gmt file to the gmt folder in the software directory

Then we just demonstrated how to do single gene. Now let's use sample information, for example, use gender classification to see the difference between men and women in KEGG pathway enrichment, and then import the sample information table

 attachments-2017-09-DjenXW8N59c07c3fce93 Click Batch Run, select the directory to save the results, and wait for a while

 attachments-2017-09-0EL2yaFY59c07c78cb01 Gorgeous results come out

 attachments-2017-09-dnXaFiMK59c07cdbd9d2 Another imperfection is the interpretation and beautification of the results. Let's wait for the editor to share with you~~



  • Published on 2017-09-19 10:12
  • Reading (57581)

56 comments

Please first Sign in Post comments
 Code farmers without code
Zhu Rangfei

Bioinformatics Engineer

118 articles

Writer list »

  1. Zhu Rangfei 118 Articles
  2. grapefruit 91 Articles
  3. Liu Yongxin 64 Articles
  4. admin 57 Articles
  5. Shengxin Analysis Flow 55 Articles
  6. SXR 44 Articles
  7. Helen Zhang 31 Articles
  8. Shuanger 25 Articles