This is a test version of Biostars. For the public version, visit https://www.biostars.org.
how to determine Genome and Gene ID in goseq for wheat

Hi, I am working on wheat RNA seq data.I have found differential expressed genes in wheat. Now I want to do gene ontology by using GOseq R package. This package take the genome ids as given in manual

pwf=nullp(genes,"hg19","ensGene")

hg19 is for human. but I could not found this type of genome id for Wheat. Can anyone tell me which one is used for wheat?

rna-seq assembly software error next-gen
• 2,372 views • link
updated by i.sudbery i.sudbery • written by muhammadaqeel24 • muhammadaqeel24

1 answer

You can use supportedOrganisms() to get a list of genomes that are supported and their IDs. However, I would be surprised if goseq had this available for wheat, in which case you'll have to supply the bias.data and the category mappings yourself. See section 5 of the GOSeq manual. In this case you'd supply the bias.data to nullp directly, rather than via genome and gene id type:

pwf <- nullp(genes, bais.data = my_length_vector)
• 0 views • link
written by i.sudbery i.sudbery

Where we can get category mappings as these are not available

• 0 views • link
written by muhammadaqeel24 • muhammadaqeel24

You can use bioMart. Go to plants.ensembl.org and select bioMart, Wheat as the the genome, and then under features/external you can select GO accession.

• 0 views • link
written by i.sudbery i.sudbery

Log in to answer this question.

More posts like this