Question: Plotting and Visualizing large number of genes into chromsome
2
gravatar for rice.researcher
4.1 years ago by
Korea, Republic Of
rice.researcher130 wrote:

Hello....

I have a list of arounfd 5000 rice genes. I would like to display those genes in rice chromosome. I tried with Rcircos, tool is good but it can map upto a few genes. So I am in trouble.!! I know about circos but it is beyond my limit to execute. Do anyone of you have an idea about method/any software suit which can be used to visualize large number of genes to chromosome?

It is ok for me if I can represent genes as even dots (rather than with gene label). My aim is to make a distribution plot chromosome .

 

 

ADD COMMENTlink modified 2.3 years ago by bernatgel1.9k • written 4.1 years ago by rice.researcher130
1

Do you really want to show all 5000 at once? Realistically, that's going to be such a busy plot that it'll be quite difficult to get anything out of it. Presumably you're looking to show some kind of connection between genes, so likely using just a subset of genes will suffice.

ADD REPLYlink written 4.1 years ago by Devon Ryan90k

I thought there might be some tools which provide a scattered plot kind of mapping ( I mean many dots) onto chromosome. My  list is having 3450 genes. In this case, I am not using connector to link genes@Devon Ryan 

ADD REPLYlink written 4.1 years ago by rice.researcher130

I would plot the density of genes along the chromosome as a line. dots is still not simple enough

ADD REPLYlink written 4.0 years ago by Ido Tamir5.0k

"I know about circos but it is beyond my limit to execute" tell us what's the problem or ask the circos mailing list.

ADD REPLYlink written 4.1 years ago by Pierre Lindenbaum120k
What do you mean by plotting a gene? Do you want to plot the gene body (exon-intro structure), or do you have a value for each gene and you want to plot them as points, or anything else?
ADD REPLYlink modified 4.0 years ago • written 4.0 years ago by Irsan6.8k
1
gravatar for bernatgel
2.3 years ago by
bernatgel1.9k
Barcelona, Spain
bernatgel1.9k wrote:

In R you can use the Bioconductor package karyoploteR (at the moment only in Bioconductor devel) to create the ideogram and plot the genes. You can plot all genes individually (either as regions or a points) or compute the gene density over the genome using a non-overlapping sliding window, for example. To create this plot you'll need to have the BSgenome.Osativa.MSU.MSU7 package installed, wich would not be necessary for human, mouse and some other organisms.

As an example we can create 5000 random "genes" and plot them. Above the ideograms we'll plot the gene density and below them the actual genes using a semitransparent color to show the accumulation of genes in certain regions.

library(karyoploteR)
library(regioneR)

#Create 5000 regions as mock genes
genes <- createRandomRegions(5000, length.mean=3000, genome="MSU7")

windows <- unlist(tileGenome(seqinfo(getBSgenome("MSU7")), tilewidth=1000000))
gene.density <- unlist(lapply(windows, function(w) {numOverlaps(w, genes)}))

kp <- plotKaryotype(genome="MSU7", plot.type=2)
kpBars(kp, data=windows, y0=0, y1=gene.density, ymax = max(gene.density), col="#DDCCFF", border="#BBAADD")
kpPlotRegions(kp, data=genes, data.panel=2, col="#00000044")

enter image description here

ADD COMMENTlink written 2.3 years ago by bernatgel1.9k
0
gravatar for Selenocysteine
4.0 years ago by
Dublin, Ireland
Selenocysteine550 wrote:

Hi, perhaps you can use a Manhattan plot? :) I don't really know what you can use for the Y axis, but...
http://www.gettinggeneticsdone.com/2014/05/qqman-r-package-for-qq-and-manhattan-plots-for-gwas-results.html

ADD COMMENTlink modified 4.0 years ago • written 4.0 years ago by Selenocysteine550
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: 1992 users visited in the last hour