get mRNAs by TxDb

Hi everyone.

I'm using transcripts() function in R environment to obtain transcript of Human from TxDb.Hsapiens.UCSC.hg38.knownGene package. I don't know that this transcripts contain intron or not. On the other word, are these transcripts mRNA or not?

tnx all,

Mansoor.

mrna

• 1,990 views • link

updated 8.9 years ago by benformatics

• written 8.9 years ago by hasani.iut6

1 answer

The transcripts() function will return coordinates for the entire range of a given transcript including introns. To get only the exonic regions use exonsBy()

## known exons for each transcript
exonsBy(TxDb.Hsapiens.UCSC.hg38.knownGene,'tx')

## known exons for each gene
exonsBy(TxDb.Hsapiens.UCSC.hg38.knownGene,'gene')

• 0 views • link

written 8.9 years ago by benformatics

More posts like this

should I remove the transcripts from annotate peak region in ATAC-seq
written 1.1 years ago by QX •

Hi all, I am working on ATAC-seq dataset and I have annotated the peak regions to the TSS with range from [-3000, 3000]. I am …
Transcript variants conversion from Ensembl ENST to NM
written 3.4 years ago by Carmen •

Hi all, I downloaded the 3'UTR folding dataset of human transcripts from Ensembl. However, the mRNA IDs do not contain the transcript variant (example: NM_001662). …
Error when trying to install a package (non-zero exit status)
written 4.3 years ago by BioC123 •

Hi, I'm trying to install TxDb.Hsapiens.UCSC.hg38.knownGene package but I fail. Unfortunately, the error message is not so helpful. Any ideas how I can install the …
What is the codification in genestrand 1 and 2?
written 4.9 years ago by Lila M

Hi there, I'm doing some peak annotation using `ChIPseeker` library(ChIPseeker) library(TxDb.Hsapiens.UCSC.hg38.knownGene) library(clusterProfiler) library(annotables) library(org.Hs.eg.db) txdb <- TxDb.Hsapiens.UCSC.hg38.knownGene peaks= readPeakFile("peaks_", header = F) peakAnno <- annotatePeak(peaks, …
Comparing different transcripts expression level of a list of genes from one condition
written 7.1 years ago by mr.monahton •

I have 20 RNA-Seq samples obtained from normal cells (without any treatment). I want to identify from the transcript quantification data whether for a certain …
How to annotate peaks using gencode27?
written 8.7 years ago by mr.monahton •

Hi guys, I'm trying to use ChIPseeker to annotate my peaks. I have used TxDb.Hsapiens.UCSC.hg38.knownGene as transcript database, but due to some issues I want …
Finding annotation of any region in human genome
written 9.0 years ago by hasani.iut6

Hi all, I had some sequences with 25 lengths. I aligned them to the human genome and found the hits. First I was interested to …
find genes name of a specific location
written 9.0 years ago by hasani.iut6

Hi. I'm trying to find the hits of an special pattern in human genome and then I want to know that the hits that has …
how to exclude set of genes from a txdb
written 9.3 years ago by tarek.mohamed

Hi All, How can I exclude single gene or a list of genes from a txdb object. I have a txdb from hg38 build and …
Gviz - creating geneModel object
written 10.1 years ago by bsmith030465

Hi, I was trying to make a geneModel object, as used in the Gviz package: > data(geneModels) head(geneModels,2) chromosome start end width strand feature gene …