Question: Create Chromosome Ideogram w/ out reference genome
gravatar for jake.gendron
2.8 years 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.


R gene • 916 views
ADD COMMENTlink modified 11 months ago by bernatgel2.6k • written 2.8 years ago by jake.gendron0

you can just make the plots yourself custom if you know all the coordinates and chrom sizes you want to use; some guidelines for doing a custom ideogram/karyotype are here: Creating chromosome karyotype plot with R and ggplot2

ADD REPLYlink written 11 months ago by steve2.6k
gravatar for genecats.ucsc
2.8 years ago by
genecats.ucsc580 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):

If you don't have web-accessible space available you can make all the necessary files and upload them to CyVerse as explained here:

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):

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:
  • Questions involving private data:
  • Questions involving mirror sites:

ChrisL from the UCSC Genome Browser

ADD COMMENTlink written 2.8 years ago by genecats.ucsc580
gravatar for bernatgel
11 months ago by
Barcelona, Spain
bernatgel2.6k wrote:

If you have the length of your chromosomes or even an approximation, you should be able to plot them with karyoploteR.

You'd need to build a custom genome (a GenomicRanges object with one range per chromosome) and you should be good to go. You can find more information on how to plot a custom genome with karyoploteR and on how to plot the gene positions in the karyoploteR tutorial.

ADD COMMENTlink written 11 months ago by bernatgel2.6k
Please log in to add an answer.


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