sharepathway

A Python package for KEGG pathway enrichment analysis with multiple gene lists


Keywords
detection, pathway, enrichment, share, multiple, gene, lists
License
Other
Install
pip install sharepathway==0.2.0

Documentation

SharePathway

SharePathway is a python package for KEGG pathway enrichment analysis with multiple gene lists.

There have been dozens of tools or web servers for enrichment analysis using a list of candidate genes from some kinds of high throughput experiments,such as Exome-Seq and RNA-Seq. But the reality is that we usually get multiple gene lists, each from one sample or patient. We can do enrichment analysis for each sample then check which pathway or module is enriched. This strategy is simple and commonly used in cancer study. But we may lose some important driver genes.

SharePathway is motivated at providing users a simple and easy-to-use tool for enrichment analysis on multiple lists of genes simultaneously, which may help gain insight into the underlying biological background of these lists of genes.

Installation

This version is for both python2 and python3. The first step is to install Python. Python is available from the Python project page . The next step is install sharepathway.

Install from PyPi using pip, a package manager for Python:

$ pip install sharepathway

Or, you can download the source code at Github or at PyPi for SharePathway, and then run:

$ python setup.py install

Usage

Assume you have put all the path of your gene list files in one summary file genelists.txt (one path per line) in the directory ~/data/. Go into this directory,open python and run the scripts below. The result will be saved in the result.html file:

import sharepathway as sp
filein="genelists.txt"
fileout="result"
sp.Run(fi=filein,fo=fileout,species='hsa',r=0.1)

The default value of species is 'hsa', represents human species. The ratio r is the min threshold. The default value is 0.01. Entrez Gene ID is supported. The result will be output to a html file.

Output Description

  • Summary

    This part summarizes the input data.

  • Details

    This part list the ranked pathways and related information as shown below.

    Column

    Description

    Pathway

    Pathway name and hyperlink to modified KEGG map

    Genes

    KEGG ID of the genes in the pathway

    pCount

    Total number of genes in the pathway

    Count

    The number of recognized genes from user input

    Ratio

    The percentage of lists that containing genes in the pathway

    Pvalue

    The combined p valude from Fisher's Method

    EASE

    EASE score defined by DAVID, from merged gene list

    FET

    P value of the Fisher Exact Test, from merged gene list

    Samples

    The number of genes in the pathway in each list

Test data

See the gene list files and genelists.txt file in data/. This is just toy data.

Contact

Author: Guipeng Li

Email: guipeng.lee(AT)gmail.com