Question: How To Calculate Conservation Score Of Given Bed Regions
0
gravatar for User 6762
9.2 years ago by
User 67620
User 67620 wrote:

Is there any simple way of extracting 44way or 29 mammals conservation levels of given genomic regions in human genome (BED format) ?

conservation • 5.4k views
ADD COMMENTlink written 9.2 years ago by User 67620
4
gravatar for Pierre Lindenbaum
9.2 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum134k wrote:

44way conservation levels are available as a set of WIG/BIGWIG files at the UCSC. For example, see: ftp://hgdownload.cse.ucsc.edu/goldenPath/hg19/phastCons46way/primates/

you can query and convert those wig/bigwig files using wigtobigwig and bigWigToBedGraph in http://hgdownload.cse.ucsc.edu/admin/exe/linux.x86_64/

================================================================
========   bigWigToBedGraph   ====================================
================================================================
bigWigToBedGraph - Convert from bigWig to bedGraph format.
usage:
   bigWigToBedGraph in.bigWig out.bedGraph
options:
   -chrom=chr1 - if set restrict output to given chromosome
   -start=N - if set, restrict output to only that over start
   -end=N - if set, restict output to only that under end
   -udcDir=/dir/to/cache - place to put cache for remote bigBed/bigWigs
ADD COMMENTlink modified 9.2 years ago • written 9.2 years ago by Pierre Lindenbaum134k

29 mammals is not there. you can find it here, already in bed format. http://www.broadinstitute.org/scientific-community/science/projects/mammals-models/29-mammals-project-supplementary-info

extraction of regions can be done with bedtools i think?

ADD REPLYlink modified 17 months ago by Ram32k • written 9.2 years ago by Tafelplankje120

i removed "29 mammals"

ADD REPLYlink written 9.2 years ago by Pierre Lindenbaum134k

I used these steps to get my phastcons scores -

  1. Get the bw file from UCSC
  2. Convert bigWigToBedGraph and bedGraph to bed (awk '{print $1 "\t" $2 "\t" $3 "\t" $4}')
  3. bedtools intersect -a ChIP_peaks.bed -b phastCons.bed

Is this the best way to do this? My final goal is to plot a PhastCons score vs Distance from Binding site plot. Please advice.

P.S: I am working with an insect genome, most of the available tools seem to be for human genomes.

ADD REPLYlink modified 17 months ago by Ram32k • written 3.9 years ago by varsha61990
1
gravatar for Pascal
9.2 years ago by
Pascal160
Germany
Pascal160 wrote:

In the title you ask how to calculate the conservation score, and in the description just how to extract conservation score from multiple regions.

For extracting you might want to take a look at biopieces. It includes a program for exactly that purpose: http://code.google.com/p/biopieces/wiki/get_genome_phastcons

Of course, the genome and conservation scores have to be downloaded from UCSC and set up for biopieces first.

ADD COMMENTlink written 9.2 years ago by Pascal160
0
gravatar for User 6762
9.1 years ago by
User 67620
User 67620 wrote:

I found an easy way, using CompleteMOTIFS tool. It gives average PhastCons scores of a given BED region. http://cmotifs.tchlab.org/cgi/pipeline_v3.cgi

ADD COMMENTlink written 9.1 years ago by User 67620
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: 972 users visited in the last hour
_