Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
cDNA or protein location to genomic DNA location Liam_Gallagher Bioinformatics 5 05-11-2012 03:34 AM
RNA-Seq: Synteny-based Mapping-by-Sequencing enabled by Targeted Enrichment. Newsbot! Literature Watch 0 03-14-2012 04:30 AM
RNA-Seq: ABMapper: a suffix array-based tool for multi-location searching and splice- Newsbot! Literature Watch 0 12-21-2010 03:00 AM
Pubmed: Microarray-based multicycle-enrichment read2009 Literature Watch 0 09-27-2009 01:20 PM
In Sequence: Bio-Reference Laboratories’ GeneDx to Launch Genetic Test Based on Illum Newsbot! Illumina/Solexa 0 07-22-2008 02:06 PM

Thread Tools
Old 03-04-2013, 07:19 PM   #1
Location: Melbourne, Australie

Join Date: Feb 2013
Posts: 11
Default Enrichment analysis based on genetic location

Hi everyone
I am integrating data analysis of data from a number of methods. Therefore my identified variants have different forms: CNVs, point mutations, indels, CpGs etc.

This, together with the fact that I wish to analyze variants located outside coding genes, complicated the use of available software.

Does anyone know a tool which can identify enrichment based on genetic location? E.g. identify if a number of my cases have variants located to e.g. chr1:1000-2000
The tool needs to function on different types of variants therefore it does not work for me if it is depending on SNPs or other features depending on one method.

An alternative approach is to annotate my variants and then evaluate for enrichment based on IDs however this is complicated by the fact that I want to include varianmts associated with all functional elements both coding and noncoding, TF binding sites, Enhancers etc.
Does anyone know of a tool that can handle this challenge?
LouiseT is offline   Reply With Quote
Old 03-04-2013, 07:37 PM   #2
Location: Seattle, WA

Join Date: Feb 2013
Posts: 45

If you have your annotations in a set of BED files, you could use the BEDOPS suite tool bedops to retrieve elements overlapping specified ranges, e.g.:

$ echo -e "chr1\t1000\t2000" | bedops --element-of -1 annotations.bed - > annotations_overlapping_chr1:1000-2000.bed

If you have a set of regions, say, in a file called my_regions.bed, then you can specify that directly:

$ bedops --element-of -1 annotations.bed my_regions.bed > annotations_overlapping_my_regions.bed

If you want both the regions and the annotations, then use the bedmap tool, e.g.:

$ bedmap --echo --echo-map --delim '\t' my_regions.bed annotations.bed > regions_with_overlapping_annotations.bed

The bedmap tool also does statistical analysis of mapped elements, if you have score or other data in your annotations. The documentation offers more detail.

Last edited by AlexReynolds; 03-04-2013 at 07:55 PM.
AlexReynolds is offline   Reply With Quote

association, data integration, enrichment, location

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 04:22 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO