GO-term comparison between two genomic regions
1
0
Entering edit mode
austrius • 0
@austrius-8634
Last seen 8.6 years ago
Switzerland

I'm looking for an efficient way to compare genes between two genomic regions of an annotated genome. Basically, I'm interested in comparing GO-terms of genes located in one genomic region to another genomic region to see whether and how they differ.

It seems that the annotated genomes implemented for Bioconductor can be found here: https://www.bioconductor.org/packages/release/data/annotation/.

But does anyone know how to extract GO-terms for a specific genomic region (e.g., bp-position 250,000-550,000 on chromosome 1) and compare it to another one (e.g., bp-position 1,500,000-1,750,000 on chromosome 1) within an annotated genome (e.g. bovine)?

 

thank you!

topgo GO annotation • 816 views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 14 hours ago
United States

You will first need a TranscriptDB object. For that you want to use makeTxDbPackageFromUCSC() in the GenomicFeatures package. Or if you don't need it to be persistent, you can use makeTxDbFromUCSC(), which will make the DB, but not the package. Once you have done that, you can do something like

gns <- genes(TxDb.Btaurus.UCSC.bosTau8.knownGene) ## assuming that's what you called it

library(org.Bt.eg.db)

gos.region1 <- select(org.Bt.eg.db, names(gns[gns %over% GRanges("chr1", IRanges(250000, 550000)),]), "GO")

Do the same for the other, and then you can compare the terms. I don't know what you mean by 'compare', nor how one would do that efficiently, so I leave that up to you.

ADD COMMENT

Login before adding your answer.

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