Take the 2-minute tour ×
Biology Stack Exchange is a question and answer site for biology researchers, academics, and students. It's 100% free, no registration required.

I have always faced a problem while analyzing small RNAseq data, at the step of adapter trimming.

Overview of small RNAseq (Illumina)

  1. RNA is size fractionated using columns or PAGE
  2. 3' and 5' adapter ligation
  3. cDNA synthesis
  4. PCR amplification
  5. Sequencing

The length of the reads is dependent on the machine and recent ones like HiSeq can provide ~200bp.

The problem, however, is with reads smaller than the max read length of the machine and this is common with small RNAs such as miRNAs (especially if the concatemer of smallRNA and the complete 3'adapter is bigger than the max read length).

The first step of analysis is trimming of the 3' adapter (Illumina Truseq: TCGTATGCCGTCTTCTGCTTGT).

Several algorithms are available to do this job and what they precisely do is to check for overlaps between the adapter sequence and 3'end of the reads and then clip the aligned region.

Now the problem is this

You can't be really sure of very small alignments because they may not really be originating from the adapters, which means that you should specify a lower limit of alignment for clipping. I usually set it as 5 (intuitively).

But if a small piece of sequence really came from the adapter then it will remain and there is no way to clip it without any doubt.

The real problem comes during aligning the reads to the reference sequence. Aligners like bowtie (which i prefer to use), generally have user defined argument for number of allowed mismatches. Bowtie generally doesn't perform very well if you allow a lots of mismatches.

Subsequently, you might lose a really valuable read.

Alternative

To avoid this problem I sometimes trim the reads to around 25nt (for miRNA profiling). This creates a new problem:

You can't really distinguish whether the read came from a pre-miRNA (a longer RNA) or from mature-miRNA (smaller RNA that arises by processing of pre-miRNA)

Does anyone have an experience or an idea about how to solve this problem?

share|improve this question
add comment

1 Answer

up vote 2 down vote accepted

I would first recommend you ask this question biostar as the subject matter you are inquiring about is much more relevant there.

That having been said, you have another option which is to use an aligner that soft clips 3' ends of reads specifically to account for adapter (or polyA, or whatever) contamination that might have flew under your radar.

STAR is such an aligner. If you search their mailing list archives for "clip", you should find several posts related to your question.

share|improve this answer
    
Thanks Steve.. But I wonder why this forum is less suitable.. my purpose of joining here was to have all biology in a single forum without having to have multiple accounts.. i have an account in seqanswers and i shall ask this question there. but i am slightly disappointed by this concept of separate niches (especially when it all comes under same area of study).. –  WYSIWYG Apr 11 '13 at 10:12
1  
@WYSIWYG: I suggested biostar, not seqanswers ;-) Still, this forum is less suitable simply because your question is more of a bioinformatics nature than strictly biology one -- and you will find more expertise related to bioinformatics on biostar than here. I would imagine that there are likely many people on biostar who are strong bioinformaticians that do not have accounts on this site, and you would benefit from getting their input. Consider, for example, why there are different SE sites for math and stats :-) –  Steve Lianoglou Apr 11 '13 at 10:19
    
thanks.. biostar is actually cool.. better than seqanswers :P –  WYSIWYG Apr 12 '13 at 3:13
add comment

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.