Question: How to get probe-specific intensity values including genomic positions from microarray ChIP-chip data?
0
gravatar for Marcel
8 months ago by
Marcel0
Marcel0 wrote:

Hello,

I am currently trying to get something like a .bed file from microarray ChIP-chip data, i.e. genomic locations and corresponding (normalized, log-fold change) intensities for all probes in this dataset: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE12283. They provide a .bed file with their called peaks (in mm8) and a .bar file with the intensities (in NCBIv33 which is mm5)

What I would have preferred is to convert the .bar file they provide to a bed file. For this I tried affy_bar2txt of CisGenome to convert it into what's basically a .bed file and then lifted the positions from mm5 to mm8. However, if I compare the output to their called peaks, I see that something must have gone wrong, since the peaks should be covered with probes, which they often are not.

I then decided to analyze the raw data using the R library Starr to get the intensities for each probe:

library(Starr)

#read original bpmap
bpmapFile <- "Mm_PromPR_v02-1_NCBIv33.bpmap"
bpmap <- readBpmap(bpmapFile)

#read CEL files and define experiment
celFiles <- list.files(path = "./CELs",full.names=TRUE)
celNames <- c("Trim28_1", "Trim28_2","Trim28_3", "Biotin_1","Biotin_2", "Biotin_3", "Biotin_4")
celTypes <- c("IP", "IP", "IP", "CONTROL", "CONTROL", "CONTROL", "CONTROL")
Trim28_chip <- readCelFile(bpmap, celFiles, celNames, celTypes, featureData=T, log.it=T)

ips <- Trim28_chip$type == "IP"
controls <- Trim28_chip$type == "CONTROL"
description <- c("Trim28vsBiotin")

#normalize and compute ratio
Trim28_loess <- normalize.Probes(Trim28_chip, method="loess")
Trim28_loess_ratio <- getRatio(Trim28_loess, ips, controls, description, fkt=median, featureData=F)

In theory, I should now have what I wanted. But since I'm still quite new to R, I don't know how to extract the genomic positions of the probes and their intensities from Trim28_loess_ratio which is an ExpressionSet object. I guess, I will have to somehow map the positions from the Trim28_chip object I created previously.

chip-chip microarray mapping R • 168 views
ADD COMMENTlink modified 8 months ago • written 8 months ago by Marcel0
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: 1457 users visited in the last hour
_