browserSession for Arabidopsis
1
0
Entering edit mode
mfahim • 0
@mfahim-7714
Last seen 8.7 years ago
Korea, Republic Of

Hi everyone,

I am trying to use run splieR for splicing analysis in my RNAseq data from Arabidopsis.

As usual, the packages and instructions in manuals are for anything but plants.

Here is my question

Arabidopsis genome is not listed in UCSC. What should be alternative to the following chunk in spliceR documentation. Cheers.

session <- browserSession("UCSC")
genome(session) <- "hg19"
query <- ucscTableQuery(session, "knownGene")
tableName(query) <- "knownGene"
cdsTable <- getTable(query)

 

rtracklayer Granges spliceR • 1.5k views
ADD COMMENT
0
Entering edit mode
@michael-lawrence-3846
Last seen 2.3 years ago
United States

Sorry, but there is no UCSC genome build for Arapidopsis. Unfortunately, it looks like there is no easy way with spliceR to use TxDb packages. But it just takes some manual work, using the transcripts() and exons() accessors, requesting columns=c("tx_id", "gene_id") and doing all the munging. The relevant TxDb package is "TxDb.Athaliana.BioMart.plantsmart25". You might also check out the "sgseq" package.

ADD COMMENT
0
Entering edit mode

thanks.. is that TAIR10?

and how should I do that?

ADD REPLY
0
Entering edit mode
> session <- browserSession("UCSC")
> genome(session) <- "araTha1"
Error in `genome<-`(`*tmp*`, value = "araTha1") : 
  Failed to set session genome to 'araTha1'

and the Error

ADD REPLY
0
Entering edit mode

Yea, sorry, I got that identifier from some sort of example of a custom hub. Arabidopsis is not in the actual UCSC database. I'll edit the answer.

ADD REPLY

Login before adding your answer.

Traffic: 764 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