[BioC] Get gene structure information: use biomaRt or other annotation packages?

[michael watson (IAH-C)]

library(biomaRt)

mart.obj <- useMart(biomart = 'ensembl', dataset = 'hsapiens_gene_ensembl')

atb <- c('ensembl_gene_id', 'ensembl_transcript_id','5_utr_start',
		'5_utr_end','3_utr_start','3_utr_end','ensembl_exon_id','exon_chrom_start',
		'exon_chrom_end','rank','phase')


locs <- getBM(attributes=atb, mart=mart.obj)
________________________________________
From: bioconductor-bounces at stat.math.ethz.ch [bioconductor-bounces at stat.math.ethz.ch] On Behalf Of Zhu, Julie [Julie.Zhu at umassmed.edu]
Sent: 20 March 2010 19:33
To: bioconductor at stat.math.ethz.ch
Subject: [BioC] Get gene structure information: use biomaRt or other annotation packages?

Hi,

I need to get the following information in one file and thinking of using
getBM method in biomaRt to get each feature type separately and then combine
them in one file. I am wondering if there is a more elegant way to do this
in biomaRt or does there exist an annotation package that contains the
following information for different species?

EnsembleGeneID  Chromosome FeatureStart  FeatureEnd  Strand  FeatureType
FeatureID

The feature types I am interested in are 5' UTR, 3' UTR and Exon.

Thank you very much for your kind help!

Best regards,

Julie


*******************************************
Lihua Julie Zhu, Ph.D
Research Associate Professor
Program in Gene Function and Expression
University of Massachusetts Medical School
364 Plantation Street, Room 613
Worcester, MA 01605
508-856-5256
http://www.umassmed.edu/pgfe/faculty/zhu.cfm

_______________________________________________
Bioconductor mailing list
Bioconductor at stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/bioconductor
Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor