SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   Annotation for S. pombe in GAGE. (http://seqanswers.com/forums/showthread.php?t=44116)

Parharn 06-11-2014 02:43 AM

Annotation for S. pombe in GAGE.
 
I am trying this workflow: "RNA-Seq Data Pathway and Gene-set Analysis Workflows
Weijun Luo" but on S. pombe and not human data.
In part 3.2 it starts with:

> library(TxDb.Hsapiens.UCSC.hg19.knownGene)

I need help with finding the corresponding package for S. pombe!

Please help!
/Parham

bigmw 06-11-2014 09:37 AM

You would need to construct your own TranscriptDb object using function like makeTranscriptDbFromGFF from the GenomicFeatures package:
library(GenomicFeatures)
?makeTranscriptDbFromGFF

You may download the GFF data for S. pombe genome from:
http://www.pombase.org/downloads/gen...asets#features
I haven’t go through the whole process, but it should work theoretically.

In addition, makeTranscriptDbFromUCSC and makeTranscriptDbFromBiomart in GenomicFeatures package are more staightfoward for species with genome annotation in UCSC or Biomart database. UCSC does not have annotation data for S. pombe genome. You may try makeTranscriptDbFromBiomart if you confirmed Biomart has the data.
More details on GenomicFeatures package:
http://www.bioconductor.org/packages...cFeatures.html

Parharn 06-11-2014 09:46 AM

Thanks for the tips bigmw, I will give it a try and see how it will goes.

Cheers,
/P

malcook 06-11-2014 09:33 PM

makeTranscriptDbFromBiomart for S. pombe
 
I just posted an approach to using Biomart at https://www.biostars.org/p/103251/

Here's the nub:

Code:

txdb<-makeTranscriptDbFromBiomart(
            ,biomart ="fungi_mart_22"
            ,dataset = "spombe_eg_gene"
            ,host="fungi.ensembl.org"
            )



All times are GMT -8. The time now is 03:08 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.