What is ChIP-seq?

R ile Bioconductor kullanarak ChIP-seq

Peter Humburg

Statistician, Macquarie University

R ile Bioconductor kullanarak ChIP-seq

R ile Bioconductor kullanarak ChIP-seq

R ile Bioconductor kullanarak ChIP-seq

R ile Bioconductor kullanarak ChIP-seq

R ile Bioconductor kullanarak ChIP-seq

The Data

  • Prostate cancer patients
    • 5 primary tumor samples
    • 3 treatment resistant samples
  • ChIP targetting Androgen Receptor
  • AR is activated by binding testosterone
  • Activated AR regulates gene expression of target genes
R ile Bioconductor kullanarak ChIP-seq

Accessing ChIP-seq data in R

  • Loading BAM files
library(GenomicAlignments)
reads <- readGAlignments('file_name')
  • Obtaining read coordinates
seqnames(reads)
start(reads)
end(reads)
  • Computing coverage
coverage(reads)
R ile Bioconductor kullanarak ChIP-seq

Accessing peak calls

  • Loading BED files
library(rtracklayer)
peaks <- import.bed('file_name')
  • Obtaining peak coordinates
chrom(peaks)
ranges(peaks)
  • Extracting peak scores
score(peaks)
R ile Bioconductor kullanarak ChIP-seq

Let's practice!

R ile Bioconductor kullanarak ChIP-seq

Preparing Video For Download...