Question: Recommended tools for positional gene enrichment (PGE) analysis?
2
gravatar for Christian
3.8 years ago by
Christian2.7k
Cambridge, US
Christian2.7k wrote:

I am looking for recommended methods that perform positional gene enrichment (PGE) analysis in human genomes. What I mean by PGE is to find out, for a given set of genes, if there is any significant positional clustering of these genes along the chromosome or genome.

I think this problem consists of two convoluted tasks: (1) identify positional gene clusters (2) test for the statistical significance of these clusters. Ideally, known biases (e.g. non-random gene density) and multiple testing correction should be accounted for.

Here is one example of a paper that addresses this problem:

De Preter et al. (2008): Positional gene enrichment analysis of gene sets for high-resolution identification of overrepresented chromosomal regions

I am looking for more recent, stand-alone implementations of such a method (e.g. an R package). A useful software would not necessarily test for positional enrichment of only genes, but could more generally work with any type of feature that has a genomic coordinate (e.g. take a list of known enhancers and identify "super enhancer" regions).

EDIT: From reading the above paper, I realized that the problem is not so much finding significantly enriched regions (this is as simple as performing a hypergeometric test), but to eliminate overlaps. Since I now found out that the authors provide their algorithm as Perl script on their home page, I will start with this method.

enrichment • 1.7k views
ADD COMMENTlink modified 2.4 years ago by wiggs380 • written 3.8 years ago by Christian2.7k
0
gravatar for wiggs38
2.4 years ago by
wiggs380
wiggs380 wrote:

Did you manage to find a stand-alone implementation in R. I am not so familiar with Perl but will also use the authours script if required.

ADD COMMENTlink written 2.4 years ago by wiggs380

What about any of the GSEA R packages? MSigDB has a positional gene set collection (C1):

Gene sets corresponding to each human chromosome and each cytogenetic band that has at least one gene. These gene sets are helpful in identifying effects related to chromosomal deletions or amplifications, dosage compensation, epigenetic silencing, and other regional effects.

http://software.broadinstitute.org/gsea/msigdb/collections.jsp#C1

ADD REPLYlink modified 2.4 years ago • written 2.4 years ago by igor7.3k

I haven't looked for it, because the Perl script works well and I am still using it in my analysis pipelines (in addition to GSEA with MSigDB positional gene sets).

ADD REPLYlink written 2.4 years ago by Christian2.7k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1680 users visited in the last hour