Which keytype from org.Mm.eg.db match the annotations in the ensembl gtf file?
0
0
Entering edit mode
6.2 years ago

Hi, I downloaded the following gtf file from ensembl: Mus_musculus.GRCm38.91.gtf

Using DESeq2 I got a res object that have ENSEMBL symbols on it. I want to use AnnotationDbi to further annotate my res object with gene symbol and gene names. However I am not sure which kind of keytypes are in my res object. I was expecting either ENSEMBL OR ENSEMBLTRANS

However when I look at the length of different object I found this:

length(keys(org.Mm.eg.db, keytype = "ENSEMBL"))
[1] 24420
length(keys(org.Mm.eg.db, keytype = "ENSEMBLTRANS"))
[1] 19245
 dim(res)
[1] 31379     6

There are more keys in my res object than in any of the ENSEMBL related keys from the org.Nm.eg.db object. How that is possible? Please find a brief summary of what I did:

I used the ensembl gtf file "Mus_musculus.GRCm38.91.gtf" it to create index and to run alignments with STAR

echo "STAR --runThreadN 12 --genomeDir /genome --sjdbGTFfile /Mus_musculus.GRCm38.91.gtf --readFilesIn /RNAseq/140311_UNC14-SN744_0402_AC3LB8ACX/${i}_R1.fastq.gz /${i}_R2.fastq.gz --outFileNamePrefix /star_mapping/aligments_forSalmon/${i} --outSAMtype BAM Unsorted  --readFilesCommand zcat" >> $i.sub

Then I created a summarize experiment object using genomicAligments

se <- summarizeOverlaps(exonsByGene, BamFileList(bamFiles), mode="Union", singleEnd=FALSE, ignore.strand=FALSE, fragments=FALSE)

Later on I run DESeq2 to perform differential expression analysis and I got the res object. Now I want to add genesymbol, For that I am extracting the keys in my result table with myKeys=rownames(res), However I'm having problem to find out which kind of keytypes those are

Hope someone can help me

Thanks

ALe

RNA-Seq DESeq2 • 4.1k views
ADD COMMENT

Login before adding your answer.

Traffic: 1756 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6