Question: Which read counts file from miRDeep2 should be considered for differential expression analysis?
1
gravatar for bioinforesearchquestions
2.8 years ago by
United States
bioinforesearchquestions280 wrote:

Hello everyone,

I would appreciate some help with my data. I want to do Differential Expression profiling for miRNA samples.

Using the miRDeep2 program and algorithm, I'm in trouble when dealing with results to make a proper table in order to use it as input for DE software such as EdgeR or DEseq2.

Which .csv should be used for read count?

miRNAs_expressed_all_samples_15_09_2017_t_07_26_27.csv

or

miRNA_expressed.csv

Also, there are duplicate values in the miRNA_expressed.csv. Should I consider the average or only one of them with highest value?

miRNA read_count precursor

ssc-let-7a 66271 ssc-let-7a-1

ssc-let-7a 66388 ssc-let-7a-2

ssc-let-7c 22747 ssc-let-7c

mirna rna-seq mirdeep2 mirbase • 1.8k views
ADD COMMENTlink modified 2.6 years ago by h.mon30k • written 2.8 years ago by bioinforesearchquestions280
1
gravatar for h.mon
2.6 years ago by
h.mon30k
Brazil
h.mon30k wrote:

You should use the miRNAs_expressed_all_samples_15_09_2017_t_07_26_27.csv as input for edgeR or DESeq.

miDeep2 runs bowtie with -a (report all alignments) and --best --strata, which means the same read will be mapped to multiple precursor miRNAs, so I believe you should use only one of the counts for each mature miRNA - I would use the one with higher counts.

ADD COMMENTlink modified 2.6 years ago • written 2.6 years ago by h.mon30k
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: 711 users visited in the last hour