You have resolved my problem, thanks a lot!
• 0 views
•
link
written
by
Arturo_M •
How Do I Use Biomart To Get Upstream Flanking Sequence For A Gene?
Hi. I'm trying to get the 100bp upstream sequences of some genes from A. gambiae. I'm using biomaRt and my query looks like this:
vector<-useMart("vectorbase_mart_13",dataset="agambiae_eg_gene")
agambiaeseq<-getBM(attributes=c('start_position','end_position','chromosome_name','strand','ensembl_gene_id','gene_flank','upstream_flank'),filters='ensembl_gene_id',value='AGAP004677', mart=vector)
I know that for the attribute 'upstream_flank' I should put the value 100 but I just don't know where.
Thank you for your attention.
• 7,380 views
•
link
updated
by
Biostar
•
written
by
Arturo_M •
•
written
by
Arturo_M •
1 answer
From the biomaRt documentation for 'getBM' it says: "Sometimes attributes where a value needs to be specified, for example upstream_flank with value 20 for obtaining upstream sequence flank regions of length 20bp, are treated as filters in BioMarts. To enable such a query to work, one must specify the attribute as a filter and set checkFilters = FALSE for the query to work." Also note that for the 'values' argument, "If multiple filters are specified then the argument should be a list of vectors of which the position of each vector corresponds to the position of the filters in the filters argument."
So, does this do what you are looking for?
library('biomaRt')
mart=useMart("vectorbase_mart_13",dataset="agambiae_eg_gene")
agambiaeseq=getBM(attributes=c('gene_flank','start_position','end_position','chromosome_name','strand','ensembl_gene_id'),filters=c('ensembl_gene_id','upstream_flank'),values=list(ENSG='AGAP004677', Upstream=100), mart=mart, checkFilters=FALSE)
The output looks like:
gene_flank start_position end_position chromosome_name strand ensembl_gene_id
ATCTCAAAATGGCAACATGTCAAACGCTAAGAAGACACCTCTTCTATATTCCACCTTGATTTGAACGGTAACATTCAGTAGTCCGTGGCTTTCGGATTAT 157348 186936 2L -1 AGAP004677
It seems to correspond to what I imagine your query might look like at the VectorBase Biomart web interface.
• 55 views
•
link
written
by
Obi Griffith
Great. Glad to help. If you find the forum useful, please stick around, contribute more good questions, answers ... and vote! ;-)
• 0 views
•
link
written
by
Obi Griffith
Log in to answer this question.
More posts like this
-
gene ensemble ID- no symbol
written by Rob •
-
annotation - biomaRt - Error in UseMethod
written by ferreira.cristina •
-
GetBM gives controversial results
written by lizaveta •
-
Plotting a single transcript with ggbio
written by longoka •
-
Ensembl Biomart Plant database
written by ####
-
biomaRt getBM() returns 0 result
written by xioli2013 •
-
Empty gene names from ensembl query using biomart
written by kindlychung •
-
Biomart Annotation
written by int11ap1
-
Biomart, Sort Genes After Their Location On The Genome
written by tomislav.ilicic •
-
Retrieve All Genes Contained Within A Specific Chromosomal Region Using R And Biomart
written by secretjess
I changed the title of your question to be more specific (just "Biomart" is too generic for a forum with many BioMart questions) and formatted your question so the code appears more nicely (four spaces before each code line is all that is needed). Welcome to biostar and thanks for your question! I've made an attempt at answer below.