Question: RIP-seq identifing reverse complement reads between smallRNA and targeted mRNA
gravatar for Konstantinos Yeles
9 days ago by
Konstantinos Yeles40 wrote:

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

ADD COMMENTlink modified 8 days ago by shawn.w.foley420 • written 9 days ago by Konstantinos Yeles40

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 REPLYlink written 9 days ago by shawn.w.foley420

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 REPLYlink written 9 days ago by Konstantinos Yeles40

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 REPLYlink modified 9 days ago • written 9 days ago by shawn.w.foley420

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 REPLYlink modified 8 days ago • written 9 days ago by Konstantinos Yeles40
gravatar for shawn.w.foley
8 days ago by
shawn.w.foley420 wrote:

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 COMMENTlink written 8 days ago by shawn.w.foley420
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1806 users visited in the last hour