Question: counting after aligning miRNA reads to miRBase mature miRNA instead of refe

inah@bi.vt.edu · 03-26-2018, 09:09 AM

I have human miRNA-seq data. In the past I have aligned these data using Bowtie2 to the reference genome (Homo_sapiens.GRCh38.dna.toplevel.fa), and I have subsequently performed counting using featureCounts with the annotation file hsa.gff3 from miRBase. . Now I have aligned the reads to the mature miRNA from miRBase (mature.fa), but when I look at the resulting bam files, the reads have a flag of 4 (segment unmapped) and there is no position information in mature.fa.

So can I use mature.fa as the reference for alignment and if so how is counting performed?

Thanks, Ina

sbarberan · 03-26-2018, 01:39 PM

Hi Ina,

I use samtools and picard-tools to get the counts for miRNAs.

Here is what I would do for bowtie2 alignment to miRBase
bowtie2 -L8 --local -x miRBase-mature-hsa-index -U FILE.fastq -S FILE.sam
samtools view -b -S FILE.sam > FILE.bam
samtools sort FILE.bam > FILE.sorted
picard-tools BuildBamIndex I=FILE.sorted
picard-tools BamIndexStats I=FILE.sorted > FILE.txt

Makes sense?

Cheers,
sergio

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Counting reads in miRNA experiment	sara_	Epigenetics	3	01-26-2016 04:48 PM
miRNA seq analysis - large numbers of non-aligning reads	lbeltrame	Bioinformatics	10	05-03-2013 09:30 AM
miRNA aligning/counting	kmcarr	Bioinformatics	28	02-08-2013 02:35 PM
aligning miRNA-seq reads to other short sequence data	impeachme2	Bioinformatics	2	07-03-2012 08:34 AM
GATK to discover Single Nucleotide Variation in mature miRNA from miRNA-Seq	Bioinfo83	Bioinformatics	0	01-31-2012 05:11 AM

03-26-2018, 09:09 AM	#1
inah@bi.vt.edu Junior Member Location: Blacksburg VA Join Date: Mar 2018 Posts: 1	Question: counting after aligning miRNA reads to miRBase mature miRNA instead of refe I have human miRNA-seq data. In the past I have aligned these data using Bowtie2 to the reference genome (Homo_sapiens.GRCh38.dna.toplevel.fa), and I have subsequently performed counting using featureCounts with the annotation file hsa.gff3 from miRBase. . Now I have aligned the reads to the mature miRNA from miRBase (mature.fa), but when I look at the resulting bam files, the reads have a flag of 4 (segment unmapped) and there is no position information in mature.fa. So can I use mature.fa as the reference for alignment and if so how is counting performed? Thanks, Ina

03-26-2018, 01:39 PM	#2
sbarberan Member Location: Santa Cruz, CA Join Date: Feb 2017 Posts: 17	Hi Ina, I use samtools and picard-tools to get the counts for miRNAs. Here is what I would do for bowtie2 alignment to miRBase bowtie2 -L8 --local -x miRBase-mature-hsa-index -U FILE.fastq -S FILE.sam samtools view -b -S FILE.sam > FILE.bam samtools sort FILE.bam > FILE.sorted picard-tools BuildBamIndex I=FILE.sorted picard-tools BamIndexStats I=FILE.sorted > FILE.txt Makes sense? Cheers, sergio