SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
GSNAP gives a lot of variants!? casshyr Bioinformatics 13 08-16-2012 02:22 PM
Barcoding.... a lot of questions!!!! soleulloa Ion Torrent 1 03-30-2012 10:12 AM
Newbie. May ask a lot of stupid questions. entropy Introductions 1 03-07-2011 06:12 PM
Issues with new XLR70 lot MissDNA 454 Pyrosequencing 7 03-01-2011 07:36 AM
New York State has a lot of nerve... ECO Personalized Genomics 3 04-21-2008 10:33 AM

Reply
 
Thread Tools
Old 08-04-2014, 09:32 AM   #1
memento
Member
 
Location: UK

Join Date: Nov 2010
Posts: 49
Default A lot of Poly-A after RNAseq completed

Hi,
Is there any obvious explanation why there would be a lot (~50%) of poly-A (or Poly-T, since this was bidirectional lib prep) in the output ?
RNA was purified using Poly-A+ enrichment (purified twice). Is it possible that the library was overamplified, hence the "PCR-winners" (polyA/T have lower GC) ? Is the library usable and the only drawback would be a lot of sequencing wasted on poly-A ?

Is it possible that the RNA was already digested, hence the only fragment retained in enrichment were the poly-A them self?

Thanks!
memento is offline   Reply With Quote
Old 08-05-2014, 06:52 AM   #2
nucacidhunter
Jafar Jabbari
 
Location: Melbourne

Join Date: Jan 2013
Posts: 1,238
Default

Providing more information might make it easier to answer the questions.
For example, was the input RNA QCed by Bioanalaser or other methods? Are poly A streches the whole read or just part of the reads? How the 1st strand synthesis was primed and what kit was used for library prep? How many cycles were used for sequencing and what was average library fragment size?
nucacidhunter is offline   Reply With Quote
Old 08-05-2014, 07:12 AM   #3
memento
Member
 
Location: UK

Join Date: Nov 2010
Posts: 49
Default

1) RNA was qced using Nanodrop + 2d gel
2) Whole reads are polyA (or T), both on 50bp SE (HiSeq) and 250bp(MiSeq, did it just out of curiosity as a part of other run). mRNA polyA is between 100-500bp so that would still make sense.
3) Priming/1st strand etc was done according to Wang2011 (PLOS), 13 cycles were used to amplify.
4) Average library size (after amp) was between 250-400 (need to subtract adapters ~100bp).
memento is offline   Reply With Quote
Old 08-06-2014, 05:33 AM   #4
nucacidhunter
Jafar Jabbari
 
Location: Melbourne

Join Date: Jan 2013
Posts: 1,238
Default

Quote:
Is there any obvious explanation why there would be a lot (~50%) of poly-A (or Poly-T, since this was bidirectional lib prep) in the output ?
Quote:
Priming/1st strand etc was done according to Wang2011 (PLOS), 13 cycles were used to amplify.
I cannot find any explanation. If bidirectional lib prep here means non-stranded, it would be contrary to the reference (Wang 2011 PLOS) because the method described in that paper is stranded library prep.

Quote:
Is it possible that the RNA was already digested, hence the only fragment retained in enrichment were the poly-A them self?
Quote:
RNA was qced using Nanodrop + 2d gel
Quote:
mRNA polyA is between 100-500bp so that would still make sense.
This is less likely for two reasons: 1- you mention that you have checked RNA prior to use by gel, so you would have noticed if it was degraded, 2- I cannot think of any way that random primers used for 1st strand synthesis would prime polyA to a extent that comprises 50% of reads.

Quote:
Is it possible that the library was overamplified, hence the "PCR-winners" (polyA/T have lower GC) ?
Less likely for such an extreme bias.

Quote:
Is the library usable and the only drawback would be a lot of sequencing wasted on poly-A ?
It would be reasonable to have explanation for this before using the data.

You have not mentioned species or estimated GC content. There might be something there to explore. FastQC plots of data are also usually good start point for troubleshooting.
nucacidhunter is offline   Reply With Quote
Old 05-07-2016, 03:51 PM   #5
ipeikon
Junior Member
 
Location: New York

Join Date: Nov 2011
Posts: 4
Default Updates?

Any updates to this thread? Did you ever find out what was going wrong?
ipeikon is offline   Reply With Quote
Reply

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 07:20 PM.


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