I have used ERCC spike-in in my RNAseq data. I have aligned my data and now I have bam files. to count the reads per gene I used htseq-count(which needs gtf file). I also have to count ERCC (I have 98 spike-in). I have fasta file of ERCCs. do you know how I can count the ERCCs?
•
written
by
alirezamomeni707 •
1 answer
You can filter out and count the spike sequences (and rRNA, and linker) before alingment with TagDust 2. In the following document on GitHub, I used it to detect ArrayControl spikes, but it also work with ERCC ones. For maximal accuracy, make sure to use the translated sequences available from the NIST, and not the plasmid insert sequences (see Patch ERCC spike sequences to get their real 5-prime ends. for the long story). In any case, if you use the TagDust approach, make sure your sequences do not contain common parts such as polyA tails.
Log in to answer this question.
-
zero counts for all genes in RNAseq data of Ferret
written by Sara
-
Proper HTSeq usage on bacterial genome. Don't quite understand --t
written by SushiRoll •
-
ERCC normalization within DESeq2
written by Assa Yeroslaviz
-
from fasta to gtf format?
written by debitboro
-
Merged vs. Ensembl GTF for htseq-count
written by kousi31 •
-
HTSeq errors in the command line
written by Morris_Chair
-
Question: Error HISAT2 output for HTseq-count
written by 11yj3312 •
-
RNA spike-ins: Their applications and bioinformatics methods for the analysis
written by Javad
-
counting the reads mapped only to the coding sequence
written by alirezamomeni707 •
-
Reference genome for mapping RNA-seq spike-in dataset
written by sarahmanderni
Ideally you could have appended the ERCC fasta to the genome and then aligned your data. Since you have not done that you would need to create a new "genome" (and a GTF file to go with it) and then align/count.
is it not possible to align to ERCC fasta and it's GTF post alignment (using aligned bam)? (aligned bam = aligned with reference fasta other than ERCC)
Since ERCC sequences should be totally diverse it should not matter what you use.
thanks. actually I have aligned to ERCC (made index from fasta file). but I do not have GTF file for that. actually this is the main problem
You make one up yourself.
You can filter them out and quantify them with BBMap's Seal using the aligned bam.