Question: RIP-seq identifing reverse complement reads between smallRNA and targeted mRNA
1
gravatar for Konstantinos Yeles
10 months ago by
Italy
Konstantinos Yeles100 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 10 months ago by shawn.w.foley1.2k • written 10 months ago by Konstantinos Yeles100
1

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 10 months ago by shawn.w.foley1.2k

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 10 months ago by Konstantinos Yeles100
1

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 10 months ago • written 10 months ago by shawn.w.foley1.2k

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 10 months ago • written 10 months ago by Konstantinos Yeles100
3
gravatar for shawn.w.foley
10 months ago by
shawn.w.foley1.2k
USA
shawn.w.foley1.2k 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 10 months ago by shawn.w.foley1.2k
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: 1044 users visited in the last hour