[BioC] Use ChIPpeakAnno to find two-sided nearest genes to a peak

Zhu, Lihua (Julie) Julie.Zhu at umassmed.edu
Mon Dec 10 20:10:03 CET 2012


Holly,

Thanks for the link! The BDPs in ChIPpeakAnno is defined purely according to
the coordinates of known genes.

Best regards,

Julie


On 12/10/12 1:30 PM, "Holly" <xyang2 at uchicago.edu> wrote:

> Julie,
> 
> A basic question to verify your definition of the bi-directional promoters is,
> did you define them purely according to the coordinates of known genes, or,
> have you referred to the experimental data, e.g. EST experiments done by
> http://www.ncbi.nlm.nih.gov/pmc/articles/PMC1853124/ ?
> 
> I learned a lot from the discussion with you. Thanks again,
> 
> Holly
> 
> 
> 
> On 12/10/2012 10:46 AM, Zhu, Lihua (Julie) wrote:
>> Dear Holly,
>> 
>> I believe that you are interested in finding the peaks that reside in
>> bi-directional promoters. If so, you can use the following functions in
>> ChIPpeakAnno.
>> 
>> BDP = peaksNearBDP(peaks, AnnotationData=TSS, MaxDistance =5000)
>> c(BDP$percentPeaksWithBDP, BDP$n.peaksWithBDP, BDP$n.peaks)
>> all.genes = union(annotated.peaks$feature, BDP$peaksWithBDP$feature)
>> where annotated.peaks is generated from annotatePeakInBatch using TSS. To
>> learn more about peaksNearBDP, please type ?peaksNearBDP in R.
>> 
>> If you just want to find genes on both side of the peaks within certain
>> distance away from the peaks, you can use the following command.
>> 
>> Annotated.peaks = annotatePeakInBatch(peaks, AnnotationData = TSS,
>> output="both",select="all", maxgap=1000000)
>> Where maxgap can be adjusted according to your needs.
>> 
>> Please let me know if this suits your needs. Thanks!
>> 
>> Best regards,
>> 
>> Julie
>> 
>> 
>> 
>> On 12/10/12 11:19 AM, "Holly" <xyang2 at uchicago.edu> wrote:
>> 
>>> Dear Lihua,
>>> 
>>> I am trying to annotate peaks for not only the genes with the nearest
>>> TSS but the ones at the other side of the peaks.
>>> Do you think I can use ChIPpeakAnno to get both sided genes for a peak
>>> region? If so, what do you suggest?
>>> Thanks a lot,
>>> 
>>> Holly
> 



More information about the Bioconductor mailing list