Question: Adding gene ID in edgeR
0
gravatar for Sharon
17 months ago by
Sharon420
Sharon420 wrote:

I am trying to add genes ID to edgeR results . Below is part of my code. Gene Ids usually comes as NA. I would like to know what could be wrong or how can I add the genes ID. All what I have is transcripts ID, so not sure if the problem is the geneid line, and how to fix? Thanks

data <- readDGE(files)
geneid <-rownames(data)
mart <- useMart("ensembl", dataset="hsapiens_gene_ensembl", host="grch37.ensembl.org", path="/biomart/martservice")
datasets <- listDatasets(mart)
mart <- useDataset("hsapiens_gene_ensembl",mart)
genes <- getBM(attributes=c('ensembl_gene_id', 'external_gene_name'),filters='ensembl_gene_id',values=geneid,mart=mart)
rownames (genes) <- genes$ensembl_gene_id
data$genes <- genes [rownames(data), ]

print (data$genes)

group <- c(rep("WT", 2), rep("Tumor", 2))
dge = DGEList(counts=data, group=group)

dge <- calcNormFactors(dge)
dge <- estimateCommonDisp(dge)
dge <- estimateTagwiseDisp(dge)
edger rna-seq genes id • 1.2k views
ADD COMMENTlink modified 17 months ago • written 17 months ago by Sharon420

Show us a few of your input gene IDs. Specifically, we need the ones getting associated with NA.

ADD REPLYlink modified 17 months ago • written 17 months ago by Devon Ryan89k

Everything is NA. This is when I print data (genes). I think the problem is how to connect between transcripts and genes? I think one step is missing but don't know what it is.

NA.49969 NA NA

ADD REPLYlink modified 17 months ago • written 17 months ago by Sharon420

I don't care what data$genes looks like, you need to post what geneid looks like.

ADD REPLYlink written 17 months ago by Devon Ryan89k

It contains the transcripts ID, which I think is the main reason for the error. So I am missing how to go from transcripts ID to genes ID? [58349] "NR_144568.1" "NR_144569.1" "NR_144628.1" "NR_144629.1"
[58353] "NR_144632.1" "NR_144633.1" "NR_144634.1" "NR_144635.1"
[58357] "NR_144637.1" "NR_144638.1" "NR_144639.1" "NR_144640.1"
[58361] "NR_144641.1" "NR_144642.1" "NR_144643.1" "NR_144644.1"
[58365] "NR_144645.1" "NR_144646.1" "NR_144682.1" "NR_144683.1"
[58369] "NR_144684.1" "NR_144685.1" "NR_144686.1" "NR_144687.1"
[58373] "NR_144688.1" "NR_144752.1" "NR_144755.1" "NR_144756.1"
[58377] "NR_144757.1" "NR_144758.1" "NR_144759.1" "NR_144760.1"

ADD REPLYlink modified 17 months ago • written 17 months ago by Sharon420
1

You're searching against biomarts collection of ensembl-gene ids using NCBI transcript IDs as input. Hence the NAs.

ADD REPLYlink modified 17 months ago • written 17 months ago by russhh4.3k
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: 700 users visited in the last hour