extra gene ids in gene count matrix than in gtf file
Hi,
I am using gtf file from ensembl database for mapping and alignment. I used stringtie to get the gene count matrix but I see many extra gene ids in my gene count matrix than those in gtf file. Is it normal? Does anyone have any suggestion? I would appreciate any suggestion.
Thanks
1 answer
Do the ids start with "MSTRG"? You will always get additional ids while using annotation. The reason for that is that the annotation file never contains 100% of complete genes and isoforms. Therefore every transcript which is not included in annotation will be assigned with 'MSTRG' identifier. If you are interested in known transcripts only you can try Salmon or Kallisto.
Log in to answer this question.
More posts like this
-
MSTRG. gene id in differentially expressed gene list!
written by khanald2114 •
-
TPM: counts for transcripts non included in ENCODE?
written by R •
-
Compatilibity between human37 and human38 gtf
written by asumani
-
Identification of presence of genes by aligning the illumina reads to reference genome with Geneious
written by rthapa •
-
How create annotation file for WGCNA for RNA-Seq data?
written by mathavanbioinfo
-
Merged vs. Ensembl GTF for htseq-count
written by kousi31 •
-
Comparison of samples of different read lengths in differential gene expression analysis
written by rthapa •
-
Issue with StringTie to get gene_count matrix for DeSEQ
written by pixie@bioinfo
-
Obtaining a GTF with a required set of gene IDs
written by spacer •
-
Ensembl id to Gene name
written by V
Thank you for your reply. The extra gene ids in my gene count matrix are not named as 'MSTRG' identifiers. I am working on rice sample. The gtf file from ensembl has ids either starting with "Os" or "ENSRNA" but in my gene count matrix there are many genes with ids starting with "EPIOSAG". And also the total number of genes in rice is around 40,000 but in my count matrix I have around total 90,000 gene ids. Do you have any idea if it is normal or it have might have been due to mistake in alignment and mapping?
Please use
ADD COMMENTorADD REPLYto respond to existing posts to keep threads logically organized.I don't recall any annotation id starting with "EPIOSAG". Can you show command that you used while creating merged gtf file with stringtie? Maybe you used name prefix for output transcripts (the -l flag). Also keep in mind that stringtie provides two count files: gene counts and transcript counts. And yes, the increased amount of rows in your gene count matrix is expected due to the fact that stringtie provides additional unknown (considering provided annotation file) genes etc. (also noncoding transcripts). If you wish to stick with ensembl gtf only you could replace stringtie merged gtf with the ensembl gtf while estimating transcript abundances.
For the stringtie, I used the following command; stringtie /scratch/user/rthapa/RNA-seq/rnaseqS10.hisat.sorted.bam -o /scratch/user/rthapa/RNA-seq/rnaseqS10.gtf -p 8 -G /scratch/user/rthapa/RNA-seq/Oryza_sativa.IRGSP-1.0.37.gtf -A /scratch/user/rthapa/RNA-seq/gene_abund -e
And to convert the gtf file, I used the below mentioned command line; python prepDE.py -i sample3a.txt
Yes, I do have both gene count matrix and transcript count matrix. And in transcript matrix too, I can see transcripts starting with "EPIOSAT".
In fact your gtf file (Oryza_sativa.IRGSP-1.0.37.gtf) does contain "EPlOSAG" ids, (short for Ensembl Plants Oryza sativa genes), check:
However they are fine identifiers of genes/transcripts. The newest gtf file of O. sativa (1.0.42) contains unified ids (i. e. without "EPlOSAG"-like ids). Check Ensembl releases page.