Question: Plotting Exons Of A Transcript Using Bioconductor Genomegraphs
gravatar for Tonig
8.6 years ago by
Tonig440 wrote:

Hi, I was looking the list but there was no response to my questions: How do I plot the exons of a specific transcript of my gene in question using Bioconductor Genographs package?


mart <- useMart(biomart = "ensembl", dataset = "hsapiens_gene_ensembl")

#getting the gene

gene <- makeGene(id = "ENSG00000171988",type = "ensembl_gene_id", biomart = mart)

I know to plot the gene, but how can i do to plot an specific transcript, let's say ENST00000399262 in order to view the exons?

Thanks in advance

ADD COMMENTlink modified 4.6 years ago by Biostar ♦♦ 20 • written 8.6 years ago by Tonig440
gravatar for Malachi Griffith
8.6 years ago by
Washington University School of Medicine, St. Louis, USA
Malachi Griffith18k wrote:

In addition to the reference manual, the user guide / vignette offers some nice tips. I was curious what the result would look like for the answer provided by @Michael. For anyone else interested, this script (GenomeGraphsExample.R) contains a complete example and creates the images shown below. I generated @Michael's two examples and an additional one with more annotation of a gene locus and its alternative transcripts.

Figure 1. All exons of the gene, followed by a specified transcript

alt text

Figure 2. All exons of the gene, followed by all transcripts

alt text

Figure 3. A different gene. This time the chromosome ideogram is included, the gene region is highlighted in red, the gene and transcript structures are displayed, exon 2 is highlighted, a main title is added, each track is labeled individually, and the transcript IDs are provided. alt text

ADD COMMENTlink written 8.6 years ago by Malachi Griffith18k
gravatar for Michael Dondrup
8.6 years ago by
Bergen, Norway
Michael Dondrup47k wrote:

Have a look at section 2.2 of the manual:

2.2 Adding alternative transcripts
To add alternative transcripts you rst have to create a Transcript object. Note that the order of the objects in the list determines the order in the plot.

> transcript <- makeTranscript(id = "ENSG00000095203", type="ensembl_gene_id", biomart = mart)
gdPlot(list(gene, transcript))

This works with the example gene ids but needs a change with your transcript id. To use your transript id you have to change type to "ensemble_transcript_id" like so:

mart <- useMart("ensembl", dataset="hsapiens_gene_ensembl")
gene <- makeGene(id = "ENSG00000171988",type = "ensembl_gene_id", biomart = mart)
transcript <- makeTranscript(id = "ENST00000399262", type="ensembl_transcript_id", biomart = mart)
gdPlot(list(gene, transcript))
ADD COMMENTlink modified 9 months ago by RamRS30k • written 8.6 years ago by Michael Dondrup47k

Can we also draw structure of novel transcripts using biomart by adding that information in it?

ADD REPLYlink written 13 months ago by waqaskhokhar99980
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: 1126 users visited in the last hour