the closest gene to a breakpoint
3
0
Entering edit mode
Bogdan ▴ 670
@bogdan-2367
Last seen 13 months ago
Palo Alto, CA, USA

Dear all, please would you advise :

given a list of chromosomes coordinates (i.e. BREAKPOINTS), and the list of genes promoters( i.e. TSS), how shall I get the closest gene TSS to a BREAKPOINT ? which function shall I use ?

I am looking for a similar function as closestBed in bedtools  ...

thanks a lot,

-- bogdan

granges plyranges • 1.7k views
ADD COMMENT
2
Entering edit mode
@michael-lawrence-3846
Last seen 3.0 years ago
United States

If you know of something that works in bedtools, then just use HelloRanges to find the equivalent. It's going to be based on nearest().

ADD COMMENT
0
Entering edit mode

Another option is to use join_nearest() or pair_nearest() from plyranges

ADD REPLY
1
Entering edit mode
Julie Zhu ★ 4.3k
@julie-zhu-3596
Last seen 13 months ago
United States
Bogdan, AnnotatePeakInBatch in the ChIPpeakAnno package should do. https://bioconductor.org/packages/release/bioc/html/ChIPpeakAnno.html Best regards, Julie
ADD COMMENT
0
Entering edit mode

Thank you Julie. I believe that I shall use the function " annoPeaks " ?

ADD REPLY
0
Entering edit mode
Bogdan, The function annoPeaks works as well. Best, Julie Activity
ADD REPLY
1
Entering edit mode
@gordon-smyth
Last seen 9 hours ago
WEHI, Melbourne, Australia
> library(edgeR)
> nearestTSS(chr="chr6", locus=300000, species="Hs")
      gene_id symbol width    tss strand distance
20824   56940 DUSP22 59299 292057      +    -7943

 

ADD COMMENT
0
Entering edit mode

Thank you Gordon. Happy to see that edgeR expands its set of functions ...

ADD REPLY

Login before adding your answer.

Traffic: 846 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6