Question: tximport using Ensembdb to create Tx2gene dataframe with bovine transcripts ?
gravatar for gilles.charpigny
7 months ago by
gilles.charpigny0 wrote:

I used Salmon to  quantify my RNA-seq. As proposed in the Salmon tutorial I downloaded and indexed a bovine transcriptome from Ensembl ( After quantification, i plan to use tximport package to aggregate transcrits to the gene level for gene-level differential expression analysis in order to  continue analysis using DESeq2. In the tximport vignette it is proposed to create a tx2gene dataframe by using ensembldb package. The recommendations of the tutorial are as follows: if you are using an Ensembl transcriptome, the easiest way to create the tx2gene data.frame is to use the ensembldb packages. The annotation packages can be found by version number, and use the pattern EnsDb.Hsapiens.vXX.  The exemple is for human data.

My question is:

how to do in bovine species knowing that I indexed on a bovine cdna database from Ensembl. Where can I find any annotation package for the bovine species and a pattern like EnsDb.Btaurus.vXX?

Thank you very much for your help





ADD COMMENTlink modified 7 months ago by James W. MacDonald49k • written 7 months ago by gilles.charpigny0

Hi Gilles

you can use the following packages :

- AnnotationDbi




Here are the following code I used for using Bovine Ensembl DB for differential expression analysis

ensembl_mart_93 <- useEnsembl(biomart="ensembl")# to use the latest ENSEMBL version
mart= useMart('ensembl') # connecting to ENSEMBL
listDatasets(mart)# List of the available datasets

mart<-useMart(biomart="ensembl", dataset="btaurus_gene_ensembl", host="",
path="/biomart/martservice", verbose = FALSE) # use of Bos taurus dataset.



ADD REPLYlink written 7 months ago by cagenet3420
Answer: tximport using Ensembdb to create Tx2gene dataframe with bovine transcripts ?
gravatar for James W. MacDonald
7 months ago by
United States
James W. MacDonald49k wrote:

I have no idea what you mean by

how to do in bovine species knowing that I indexed on a bovine cdna database from Ensembl

will you clarify? As for the second question:

> library(AnnotationHub)

> hub <- AnnotationHub()
updating metadata: retrieving 1 resource
  |======================================================================| 100%

snapshotDate(): 2018-04-30
> query(hub, c("taurus","ensdb"))
AnnotationHub with 6 records
# snapshotDate(): 2018-04-30
# $dataprovider: Ensembl
# $species: Bos Taurus
# $rdataclass: EnsDb
# additional mcols(): taxonomyid, genome, description,
#   coordinate_1_based, maintainer, rdatadateadded, preparerclass, tags,
#   rdatapath, sourceurl, sourcetype
# retrieve records with, e.g., 'object[["AH53189"]]'

  AH53189 | Ensembl 87 EnsDb for Bos Taurus
  AH53693 | Ensembl 88 EnsDb for Bos Taurus
  AH56658 | Ensembl 89 EnsDb for Bos Taurus
  AH57731 | Ensembl 90 EnsDb for Bos Taurus
  AH60745 | Ensembl 91 EnsDb for Bos Taurus
  AH60948 | Ensembl 92 EnsDb for Bos Taurus
> btens <- hub[["AH60948"]]
downloading 1 resources
retrieving 1 resource
  |======================================================================| 100%

loading from cache
> btens
EnsDb for Ensembl:
|Backend: SQLite
|Db type: EnsDb
|Type of Gene ID: Ensembl Gene ID
|Supporting package: ensembldb
|Db created by: ensembldb package from Bioconductor
|script_version: 0.3.2
|Creation time: Sun Apr  8 11:33:37 2018
|ensembl_version: 92
|ensembl_host: localhost
|Organism: bos_taurus
|taxonomy_id: 9913
|genome_build: UMD3.1
| No. of genes: 24616.
| No. of transcripts: 26740.
|Protein data available.
ADD COMMENTlink written 7 months ago by James W. MacDonald49k
Please log in to add an answer.


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