SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Mixing strand-specific and non strand-specific reads in RNA-seq alexanderxcy Bioinformatics 2 08-27-2013 01:39 AM
How to identify sense/antisense pairs with strand-specific RNA-seq data vivienne_lovely Bioinformatics 1 07-08-2013 11:08 PM
Quantifying antisense transcription from non-strand Specific RNA seq aslihan Bioinformatics 1 05-06-2013 06:39 AM
RNA-Seq: Comprehensive comparative analysis of strand-specific RNA sequencing methods Newsbot! Literature Watch 0 08-17-2010 02:00 AM

Reply
 
Thread Tools
Old 02-19-2014, 04:17 AM   #1
Pengfei Liu
Member
 
Location: Beijing

Join Date: Nov 2012
Posts: 16
Default Antisense RNA analysis from strand specific RNA-seq

Dear all

I am working on 12 samples of strand specific bacterial RNA-seq data(4 reps for 3 treatment). I have finished the DE analysis by using BWA-MEM for mapping, HTSeq for read count summarization and DESeq2 for DE test.

Now, I want to investigate the antisense RNA (find where are they located and do different expression quantification), because my libraries were constructed from strand specific protocol. And test from RseQC also confirmed that they were strand specific.

I searched the web and list several steps I need to work through:

1,separate reads/fragments mapped to + and -strand(my data is PE, so maybe Fragments) (Do I need to do?) 2,identify reads that mapped on the opposite strand of the genome annotation. (Here I mean if a gene is annotated on + strand of the genome, I need to identify reads mapped to the - strand of this location, if there were);

3, summarized those counts!

4, there might be some contamination, so a test was needed to calculation the base level and eliminate predictions below noise level.

5, use IGV to visualize each of them and check them.

6, use DESeq2 or edgeR to do differential expression test of those predicted antisense RNA

7, Moreover, do some wet expression to verify the expression those antisense RNAs.

Were these steps enough to do the job, please give some suggestions.

I get to know that SeqMonk has an antisense RNA analysis pipeline (it has the 3-UTR’ overlapping drawbacks.) And I am going to use seqmonk for this task. But for the steps list above, I still have some questions that confused me.

1, the results of HTSeq has a large proportion of reads get the marks of “ no features”.

I think antisense RNAs were in this proportion, because I had set the strand when using HTSeq for fragments count. And those fragments mapped to opposite strand of a gene would get “no features”, am I right?

2, how to separate reads mapped to the + and – strand and how to summarize antisense counts?

Could I use HTSeq to summarize them by just set the strand parameter opposite to which I used to summarize reads mapped to transcripts/genes.

Thanks!
Pengfei Liu is offline   Reply With Quote
Reply

Tags
antisense transcription, directional mrna seq

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 02:56 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO