Question: about UCSCtrack in Gviz
0
gravatar for Bogdan
3.0 years ago by
Bogdan550
Palo Alto, CA, USA
Bogdan550 wrote:

Dear Florian, and all,

I am running the following code in Gviz in order to visualize the genes and the genetic variants :

ENSEMBL genes that are displayed with BiomartGeneRegionTrack are shown, but not the RefSeq or GENECODE genes; it would be very helpful if you could let me know how we can fix the display.

many thanks,

bogdan

 


library("Gviz")
library("biomaRt")
library("BSgenome.Hsapiens.UCSC.hg38")

transcript <- "ENST00000262189"

##

hg38 <- useMart("ensembl", "hsapiens_gene_ensembl")
transcript_ENSEMBL <- BiomartGeneRegionTrack(biomart = hg38, transcript = transcript)

### making a sequence track

sequence_track_ucsc <- SequenceTrack(Hsapiens)

 

start <- 152273760
end <- 152273780

 

ucscGene <- UcscTrack(genome="hg38", track="knownGene", trackType="GeneRegionTrack", chromosome="chr7")

refGene <- UcscTrack(genome="hg38", track="refGene", trackType="GeneRegionTrack", chromosome="chr7")

plotTracks(c(transcript_ENSEMBL, ucscGene, refGene,  sequence_track_ucsc), chromosome="chr7", from=start, to=end)

 

 

gviz • 1.2k views
ADD COMMENTlink modified 3.0 years ago by Robert Ivanek590 • written 3.0 years ago by Bogdan550
Answer: about UCSCtrack in Gviz
1
gravatar for Robert Ivanek
3.0 years ago by
Robert Ivanek590
Switzerland
Robert Ivanek590 wrote:

You need to provide additional parameters for the UcscTrack function (exon starts and ends, transcript names)
 

ucscGene <- UcscTrack(genome="hg38", track="knownGene", trackType="GeneRegionTrack", chromosome="chr7", rstarts = "exonStarts", rends = "exonEnds", gene = "name", symbol = "name", transcript = "name", strand = "strand", from=start, to=end)

refGene <- UcscTrack(genome="hg38", track="refGene", trackType="GeneRegionTrack", chromosome="chr7", rstarts = "exonStarts", rends = "exonEnds", gene = "name", symbol = "name2", transcript = "name", strand = "strand", from=start, to=end)

 

 
ADD COMMENTlink written 3.0 years ago by Robert Ivanek590

thanks Robert !

ADD REPLYlink written 3.0 years ago by Bogdan550
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 16.09
Traffic: 427 users visited in the last hour