Question: creating a miRNA gtf by mirdeep2
1
gravatar for A
3.8 years ago by
A3.8k
A3.8k wrote:

hi,

how I can create my own annotation gtf for miRNA DE analysis by mirdeep2?

thank you

mirna rna-seq mirdeep2 gtf • 1.4k views
ADD COMMENTlink modified 3.8 years ago by John12k • written 3.8 years ago by A3.8k
1

I'm not sure about mirdeep2, but there is something else you could try. If I'm not mistaken, based on your other posts, miRNA annotation is not available for your species of interest. A possibility might be to map all the reads (you probably already did that?), convert bam files to bed files (e.g. bedtools), merge all bed files (bedtools) and create your own annotation file... It could be used for differential expression analysis, but interpretation would be harder. But it would allow you to identify which coding elements are differential expressed between your two conditions. Investigating the function of differential expressed small RNAs (perhaps in other species which richer annotation) is then a next step...

ADD REPLYlink written 3.8 years ago by WouterDeCoster44k
1

thank you for paying attention. yes in mirbase there is not miRNA annotation. in ensembl gff3 contains miRNA but when I converted that to gtf, tophat gave me error

ADD REPLYlink written 3.8 years ago by A3.8k
1

I'm not sure if you posted it already elsewhere, but if you can be more specific about that error then that's perhaps something that can be fixed. Also, a gtf isn't absolutely required for mapping...

ADD REPLYlink written 3.8 years ago by WouterDeCoster44k
1

thank you. I will post errors. but I need gtf for featurecounts or cufflinks

ADD REPLYlink modified 3.8 years ago • written 3.8 years ago by A3.8k
1

Definitely. But if your annotation is poor, you can create your own gtf as described above based on all reads in all samples combined.

ADD REPLYlink written 3.8 years ago by WouterDeCoster44k
1

sorry I converted my gff3 to gtf but there is no more miRNA in my resulted gtf

head(gtf[,1:9])

I JCVI exon X216 X836 . X. ..1

1 I JCVI CDS 216 836 . + 0

2 I JCVI exon 3952 4253 . + .

3 I JCVI exon 4774 4963 . + .

4 I JCVI CDS 3952 4253 . + 0

5 I JCVI CDS 4774 4963 . + 1

6 I JCVI exon 6805 6881 . - .

ADD REPLYlink modified 3.8 years ago • written 3.8 years ago by A3.8k
1

sorry if I use gff3 contains miRNA from ensembl for featurecounts, is this something like I have used miRNA spesific annotation like which there are in mirbase or this gff3 is different and DE is not proper?

ADD REPLYlink written 3.8 years ago by A3.8k
2
gravatar for John
3.8 years ago by
John12k
Germany
John12k wrote:

Congratulation Angel for getting the rarest badge on all of Biostars! :)

http://i.imgur.com/HuL4eGU.png

ADD COMMENTlink written 3.8 years ago by John12k
2

thank you John :) :) :)

ADD REPLYlink written 3.8 years ago by A3.8k
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: 919 users visited in the last hour