Question: RNA-SEQ count matrix: conversion of ensemble transcript ID to gene ID
0
gravatar for Duckula
17 months ago by
Duckula20
Duckula20 wrote:

Hello all, I would like to do DEG analysis on my count matrix using DESEQ2 package in R. My reads are in Ensembl Transcript ID and I know that I need to convert them to gene ID and sum up the reads belonging to the same gene before DEG, I know the package tximport for kalisto files but I have never done such directly from count matrix. Does anyone have some experience to share that how can I proceed?

Thanks a lot! p.s. the species is mmusculus (mouse)

id conversion rna-seq deseq2 R • 1.3k views
ADD COMMENTlink modified 17 months ago by arup1.9k • written 17 months ago by Duckula20
2
gravatar for arup
17 months ago by
arup1.9k
India
arup1.9k wrote:

The following is an example code for annotations using AnnotationDbi R package. Change res to the DESeq result object and the annotation file org.Mm.eg.db to the desired organism.

library("AnnotationDbi")
#source("https://bioconductor.org/biocLite.R")
#biocLite("org.Mm.eg.db")
#Mouse annotations
library("org.Mm.eg.db")
#columns(org.Mm.eg.db)
res$symbol = mapIds(org.Mm.eg.db,
                    keys=row.names(res), 
                    column="SYMBOL",
                    keytype="ENSEMBL",
                    multiVals="first")
res$name =   mapIds(org.Mm.eg.db,
                    keys=row.names(res), 
                    column="GENENAME",
                    keytype="ENSEMBL",
                    multiVals="first")
write.table(res,file='res_with_annotation.tsv', sep ="\t",quote = F)
ADD COMMENTlink written 17 months ago by arup1.9k

Thanks a lot! I will try it out to see the results !!

ADD REPLYlink written 17 months ago by Duckula20

Hi again! thanks a lot ! but this package only converts ensemble gene IDs to gene symbol but I have a ensemle transcript IDs and I also need to sum them up if they belong to the same gene. Do you have any other suggestions?

ADD REPLYlink written 17 months ago by Duckula20

Then a better idea will be to fetch them from reference GTF file.

ADD REPLYlink written 17 months ago by arup1.9k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2388 users visited in the last hour