Question: Create Chromosome Ideogram w/ out reference genome
0
gravatar for jake.gendron
15 months ago by
jake.gendron0 wrote:

Hey guys,

I'm relatively new to R so bear with me. I'm trying to map gene locations to a chromosome, essentially looking at one chromosome at a time. I have looked at all the useful tools including Gviz, quantsmooth, and even in IdeoViz. The problem is I'm using Oryza sativa as my reference genome, and all of these packages require a separate download usually from UCSC which lacks my reference genome. Is there a way around this problem? I know essentially the length of Chr. 11 and have my data formatted as such.

Gene_ID         Chr   Start   End     Annotation
Os11t0100600-01 11  14258   16117   Alpha/beta hydrolase family protein.
Os11t0103050-00 11  127297  129067  Hypothetical protein.

Is there an easy was to get around this issue? Should I just create a box plot given the range of Chr 11 and then write a for loop to map each of the regions to their given location? I'm just struggling to find the simplest way to do this, so any help would be greatly appreciated.

Thanks

R gene • 533 views
ADD COMMENTlink modified 15 months ago by genecats.ucsc560 • written 15 months ago by jake.gendron0
0
gravatar for genecats.ucsc
15 months ago by
genecats.ucsc560
genecats.ucsc560 wrote:

Usually the chromosome ideograms at UCSC indicate results of Giemsa staining experiments, not gene locations.

If you want to get a whole chromosome view of your data you could create an assembly hub, and create a simple 3 column bigBed file out of your data file, and the just look at the whole chromosome view of your bigBed in dense mode.

Creating an assembly is described here (there is a simple cytoband creation section as well): http://genomewiki.ucsc.edu/index.php/Assembly_Hubs

If you don't have web-accessible space available you can make all the necessary files and upload them to CyVerse as explained here: http://genome.ucsc.edu/goldenPath/help/hgTrackHubHelp.html#Hosting

and then view your assembly hub that way.

Lastly, you may also want to look into our list of automatically generated NCBI assembly hubs and see if any of the Oryza Sativa assembly hub there will work for your purposes (you may need to rename your chromosomes to match whatever NCBI uses): http://genome-test.cse.ucsc.edu/gbdb/hubs/genbank/plant/plant.ncbi.html

Please note that these automatically generated hubs are a work in progress and liable to change.

If you have further questions about UCSC data or tools feel free to send your question to one of the below mailing lists:

  • General questions: genome@soe.ucsc.edu
  • Questions involving private data: genome-www@soe.ucsc.edu
  • Questions involving mirror sites: genome-mirror@ose.ucsc.edu

ChrisL from the UCSC Genome Browser

ADD COMMENTlink written 15 months ago by genecats.ucsc560
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: 1888 users visited in the last hour