RIP-seq identifing reverse complement reads between smallRNA and targeted mRNA
1
1
Entering edit mode
2.5 years ago

Dear Biostars, I want to find a way to search for reads that have reverse complements inside a library from a RIP-seq experiment. I have found this related question: Question: Finding complementarity between mRNA and rRNA using BLAST but nucleotide blast seems to be a bit different now. Current blastn searches for reverse complement alignments too? Is there a package that I could do it in R?

thank you in advance

RIP-seq GenomicRanges R smallRNA • 527 views
ADD COMMENT
1
Entering edit mode

I'm not sure I quite understand the question, what specifically are your inputs and your outputs? Reverse complements of what? Are you inputting a RIP-seq bam file, and an annotated gtf file and looking for reads that have the reverse complement of the gtf file? Are you looking for reads with both forward and reverse orientations in the RIP-seq experiment?

If you have a bam file and you want to find reads in the reverse complement of an annotations file, you can just use

intersectBed -a input.bam -b annotation.bed -S

The -S flag specifically searches the A file for the reverse complement of the B file.

ADD REPLY
0
Entering edit mode

I have reads from a BAM file of length <45 and reads 45< rl < 80. I have aligned them and now I'm looking to check if between these sets there is a matching.

ADD REPLY
1
Entering edit mode

You should still be able to run an intersectBed command from the bedtools suite, I just checked and both the A and B files can be bam. You'll want to be mindful of the -f flag (minimum overlap required), maybe requiring that the entire <45bp read overlaps with the 45-80bp read, but that's up to you (the default is 1bp overlap). Just be sure to raise the -S flag to require that the reads are reverse complements.

I'm also assuming that strand-specific library prep was used for these libraries, otherwise there would be no way to determine orientation.

ADD REPLY
0
Entering edit mode

Then I could use R to subset the BED and then use it with bed tools Thanks. Is it possible to post it as an answer?

ADD REPLY
3
Entering edit mode
2.5 years ago
shawn.w.foley ★ 1.3k

One option is bedtools intersect or intersectBed which can be run on two bam files.

intersectBed -a shortReads.bam -b longReads.bam -S -bed

The -S flag will force the opposite strandedness (as opposed to the -s flag that forces same strandedness). You'll want to play with the -f flag that requires the fraction of A that must be overlapped (the default is 1bp, you may want to raise -f 1.0 to require the entire short read overlaps the long read).

ADD COMMENT

Login before adding your answer.

Traffic: 2526 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