wiki:SOPs/geneFeatureHeatmaps

Version 3 (modified by gbell, 11 years ago) ( diff )

--

Creating genome feature heatmaps from sequencing experiments

This type of figure can be used to visualize chIP or expression expression data in selected genome regions (rows) that have been divided into windows (columns) using either similar or relative feature widths.

Sample image:

  1. Create a BED file with all the regions you'd like to display.
  2. Create a SAM/BAM file by mapping your reads with a tool like bowtie or tophat.
  3. Split the regions you'd like to display into windows,
    • labeling each window in a systematic way
      • Region 1 would include windows REGION_01_001 ... REGION_01_100
      • Region 2 would include windows REGION_02_001 ... REGION_02_100, etc.
    • The BaRC script /nfs/BaRC_Public/BaRC_code/Perl/split_genome_regions_into_windows.pl can do this.
  4. Quantify reads that map to each window of each region.
    • Method 1: Use coverageBed (ex: coverageBed -split -abam accepted_hits.bam -b region.windows.bed > out)
    • Method 2: Use cufflinks (after converting region.windows.bed into a GTF file)
    • Other methods?
  5. Convert list into matrix (with code such as /nfs/BaRC_Public/BaRC_code/Perl/convert_coverageBed_output_to_matrix.pl) like
NAME 1 2 3 4 5 6
REGION_0101113300
REGION_0214161111544
REGION_031000020
  1. Add another column at the beginning with the same region IDs and save the tab-delimited text file with a cdt extension.
  2. Open the cdt matrix file in Java TreeView, adjust the pixel settings and print the figure.

Any other approaches to do this?

Attachments (1)

Download all attachments as: .zip

Note: See TracWiki for help on using the wiki.