Question: mirdeep2: UMI information lost when reads are collapsed
gravatar for Rituriya
6 months ago by
Rituriya20 wrote:

Hi All,

I have UMI extracted miRNA reads which I want to analyse using mirDeep2 for known and novel miRNA discovery. But when I give these as input, collapses reads drastically and brings down 7 million reads to 1200 reads with count tag on each read.

My question:

1) Is there a way to retain UMI information on the Fastq header so as to do deduplication later post alignment by mirdeep2?

2) Has anyone successfully analyzed UMI tagged miRNA data using mirDeep2?

Thank you, Pratibha.

mirna mirdeep2 umi • 303 views
ADD COMMENTlink written 6 months ago by Rituriya20

You might need to align to a genome with a regular aligner, use something like UMI_tools on that alignment, then put those reads through mirDeep2.

ADD REPLYlink written 6 months ago by swbarnes26.0k

That was exactly my thought process initially, swbarnes2. But if I want to do that, requires a reads_collapsed.fa and .arf file mandatorily. Shall I convert bowtie output to arf format using command: reads_vs_refdb.bwt > reads_vs_refdb.arf

Let me try and see if it works downstream. Do let me know if you think otherwise.

Thanks, Rituriya.

ADD REPLYlink modified 6 months ago • written 6 months ago by Rituriya20

I found to be exactly what I need but it does not retain UMI tag in header (so I am back to square one, where I will have to manually add that using programming) and also I am unable to run the perl script:

./ -i /home/xyz/aligned-bowtie1.sam -o xyz_reads_collapsed.fa -a reads_aligned-bowtie1.arf -c

It always throws an error: Sam file not found. I checked the file is very much there. Further, I learnt that all SAM files will not work from this link, but I need to know which SAM fields need to be present/absent. Any idea?

ADD REPLYlink written 6 months ago by Rituriya20
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1307 users visited in the last hour