Small RNA data analysis
Hi,
I am in the process of writing a script to automate small RNA-seq data analysis (18-30nt read length) to identify miRNAs and other types of small RNAs. Based on my reading this is what I came up with.
- Trim adaptors with Cutadapt.
- Alignment using Bowtie v1 (old version, not bowtie2).
- Count miRNAs using featureCounts (in the case of miRNA).
Does this look good? If not what other steps are recommended?
Thanks in advance,
KB
• 1,958 views
•
link
updated
by
Biostar
•
written
by
jkkbuddika
•
written
by
jkkbuddika
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
How does featureCounts handle multimapped reads from Bowtie2 -k 100 in default mode?
written by rachana •
-
miRNA low mapping ratings
written by Ant •
-
Is removing rRNAs, tRNAs, etc really a necessary step in the miRNA analysis?
written by tianshenbio
-
Very low alignment rate for miRNAs reads against known miRNAs from mirbase
written by LynxLynx •
-
Identification of piRNA/siRNAs from small RNA-seq data
written by jkkbuddika
-
Is there any recommended steps/pipeline for identifying differential miRNA expression profile?
written by bioinforesearchquestions
-
Non-coding RNA prediction using genomic data
written by nancydong20
-
Many 0 length reads after trimming
written by marcon •
-
miRNA seq trimming
written by bioLife •
-
Mirna Differential Expression
written by Ashutosh Pandey
Some miRNA kits require kit specific data handling (e.g. they may use a 3' adapter that is directly ligated to ends of miRNA and needs to be specifically removed, only reads containing that adapter may be real miRNA). I would definitely check on that possibility.
Thanks a lot. Will look into that. Do you think Bowtie v1 is the best choice for alignment, over let's say STAR or bowtie2?
Yes. You want to do ungapped alignments due to the small sizes of miRNA.
Awesome. Thanks a lot for all the suggestions.
@genomax, I haven't tried bowtie yet. If needed, what parameters would you recommend changing for alignment? I was thinking first running in the default mode and see.
There is a study looking at different aligners for performing miRNA alignment which could be of use for you (table 3 discusses parameters used): Evaluation of microRNA alignment techniques
That's helpful. Thanks. I ran with default settings and the alignment rate looked pretty good!