Question: Is this workflow for Mir-seq data correct ?
0
gravatar for sabaghianamir70
6 months ago by
iran
sabaghianamir7010 wrote:

Hi

I have small RNA-seq data and i have a workflow for analyzing it, any suggestion about it or its fine to work with.

Trim Galore > Collapse > RNA STAR > FeatureCount > Limma-voom

rna-seq mir-seq • 185 views
ADD COMMENTlink modified 6 months ago by Biostar ♦♦ 20 • written 6 months ago by sabaghianamir7010

I think if you collapse, you are going to have problems, because more or less every read from a miRNA will have the same sequence, then collapasing will mean you only ever get one read for each isoMir.

Also, STAR isn't really neccessary, because miRNAs arn't spliced. You should be okay using BWA or bowtie2.

Becareful when doing featureCounts. Many miRNAs multi-map because they are present in more than one genome location, but featureCounts will automatically discard these reads unless told not to.

ADD REPLYlink written 6 months ago by i.sudbery8.2k

So why in every workflow i look they use collapse, i dont understand. and instead of feacurecount, you suggest using HT-seq ?

ADD REPLYlink written 6 months ago by sabaghianamir7010

Most workflows collapse but still keep track of the number of identical sequences (so you don't end up with the problem of one read for each isoMir).

You can still use featureCounts and specify the -M flag to count multi-mapping reads. The issue of multi-mapping depends on what you're mapping against (e.g. are you aligning against the entire reference genome?).

ADD REPLYlink modified 6 months ago • written 6 months ago by dsull1.4k

yes, Human reference genome. and another thing is, in bowtie should i use mirbase gff file(hsa) instead of Ucsc gff file, which i download from UCSC website ?

ADD REPLYlink written 6 months ago by sabaghianamir7010

because bowtie isn't a spliced aligner, it doesn't use a GFF file.

ADD REPLYlink written 6 months ago by i.sudbery8.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: 1278 users visited in the last hour