Что такое ChIP-seq?

ChIP-seq с Bioconductor в R

Peter Humburg

Statistician, Macquarie University

ChIP-seq с Bioconductor в R

ChIP-seq с Bioconductor в R

ChIP-seq с Bioconductor в R

ChIP-seq с Bioconductor в R

ChIP-seq с Bioconductor в R

Данные

  • Пациенты с раком предстательной железы
    • 5 образцов первичной опухоли
    • 3 образца с устойчивостью к лечению
  • ChIP нацелен на рецептор андрогенов
  • Рецептор андрогенов активируется при связывании тестостерона
  • Активированный рецептор регулирует экспрессию целевых генов
ChIP-seq с Bioconductor в R

Работа с данными ChIP-seq в R

  • Загрузка BAM-файлов
library(GenomicAlignments)
reads <- readGAlignments('file_name')
  • Получение координат ридов
seqnames(reads)
start(reads)
end(reads)
  • Вычисление покрытия
coverage(reads)
ChIP-seq с Bioconductor в R

Работа с вызовами пиков

  • Загрузка BED-файлов
library(rtracklayer)
peaks <- import.bed('file_name')
  • Получение координат пиков
chrom(peaks)
ranges(peaks)
  • Извлечение оценок пиков
score(peaks)
ChIP-seq с Bioconductor в R

Давайте потренируемся!

ChIP-seq с Bioconductor в R

Preparing Video For Download...