Question: downloading annotated mouse genome
0
gravatar for Affan
5.5 years ago by
Affan290
Canada
Affan290 wrote:

I want to download the mouse genome (mm9) with some basic annotations: Intergenic region, promoter region, TSS, gene start, gene end, CpG islands and so on.

Basically, I have split up the genome (mm9) into 200bp bins. Each bin has some characteristic associated with it (mainly histone modifications). I'd like to traverse the genome in 200bp bins and look at the annotation in the specific bin (mainly to relate histone modifications with certain parts of the genome - see Ernst & Kellis).

I think I can use biomart (which has a awesome package in bioconductor), but I can't see to get the ENTIRE genome. Biomart can give me promoter regions, TSS, gene start/end and so on. Can I assume that anything outside this area is intergenic region?

genome • 2.2k views
ADD COMMENTlink modified 5.5 years ago by Devon Ryan96k • written 5.5 years ago by Affan290
0
gravatar for Devon Ryan
5.5 years ago by
Devon Ryan96k
Freiburg, Germany
Devon Ryan96k wrote:

There isn't a single annotation file with everything you want, since most people don't want all of that at the same time. Since it seems you're already using bioconductor, you might find the TxDb.Mmusculus.UCSC.mm9.knownGene package useful. It will allow you to get Promoters, TSS, gene starts/ends and so on. You can also get intergenic regions, but there's not a single command to do it (see my reply here). For things like CpG islands, just download the track in BED format from UCSC and load it with the rtracklayer package. Most other things you might want are also available via UCSC.

ADD COMMENTlink written 5.5 years ago by Devon Ryan96k

Thanks, I think this will be easier to use than biomart. Although I'll probably go ahead and use biomart as a "double check"/control.

Do you know if TxDb.Mmusculus.UCSC.mm9.knownGene contains RefSeq or Ensembl genes? I also use BSgenome.Mmusculus... instead of TxDb.Musculus... Are they the same? (nvmd, BSGenome didn't install on my R. TxDb did I guess TxDb is an older version?)

 

ADD REPLYlink modified 5.5 years ago • written 5.5 years ago by Affan290
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1946 users visited in the last hour