This is a test version of Biostars. For the public version, visit https://www.biostars.org.
getting miRBase microRNAs gff Aspergillus fumigatus

hi,

I am performing differential expression of pathogen fungi Aspergillus fumigatus miRNAs in response to human blood infection.

first I trimmed adapters, then mapped trimmed fastq against whole genome fasta of Aspergillus fumigatus by bowtie2 now I am going to extract read count by htseq-count then i need miRBase microRNAs gff for Aspergillus fumigatus but after too much googling I did not find. do you know what can I do please? where I can find miRBase microRNAs gff or Aspergillus fumigatus or can I Aspergillus fumigatus gf3 from enseble neglecting I am working with miRNA????

thank you

rna-seq mirna
• 2,492 views • link
written by zizigolu zizigolu

2 answers

I should download Aspergillus_fumigatus.CADRE.32.gff3 from ensembl but it contains not only miRNA but also CDS and exons then how I can extract miRNA from this file?

• 1 views • link
written by zizigolu zizigolu

You want to extract just miRNA coordinates 

awk -F "\t"  '$3 == "miRNA"{print}' gff >out_miRNA

But the strange thing is biotype of these miRNA are protein coding. miRNAs are non coding right ?

• 1 views • link
written by Prasad Prasad

thank you I am performing DE for miRNAs then after providing bam file by bowtie2 I need a gff to extract read counts which should be unique miRNAs, am I wrong please???

• 1 views • link
written by zizigolu zizigolu

then no need to seperate miRNA coordinates. you can use complete gff to get the readcount using stringtie(faster than latter) or cufflinks. each instances of miRNA in the gff are unique

• 1 views • link
written by Prasad Prasad

thank you

may please tell me what should this code in R

awk -F "\t" '$3 == "miRNA"{print}' gff >out_miRNA

• 1 views • link
written by zizigolu zizigolu

I don't think Aspergillus fumigatus is in the miRBase database:

ftp://mirbase.org/pub/mirbase/CURRENT/genomes/

• 1 views • link
written by robm9119 robm9119

yes I know then what can I do??

• 1 views • link
written by zizigolu zizigolu

Log in to answer this question.

More posts like this