Question: RNA-SEQ count matrix: conversion of ensemble transcript ID to gene ID
0
gravatar for Duckula
21 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.6k views
ADD COMMENTlink modified 21 months ago by Arup Ghosh2.1k • written 21 months ago by Duckula20
2
gravatar for Arup Ghosh
21 months ago by
Arup Ghosh2.1k
India
Arup Ghosh2.1k 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 21 months ago by Arup Ghosh2.1k

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

ADD REPLYlink written 21 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 21 months ago by Duckula20

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

ADD REPLYlink written 21 months ago by Arup Ghosh2.1k
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: 1321 users visited in the last hour