This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Ensembl version for mapping and annotating genes

Hello,

I am mapping my data to mmGRCm38 v92. However, when I use tximport and biomart to annotate my transcripts to gene IDs I think that it annotates using the most recent ensembl version for mouse. Is that a problem?

I use the following code:

mart <- biomaRt::useMart("ensembl", "mmusculus_gene_ensembl", host = "uswest.ensembl.org")

Should I be using the code below to make sure that I am using the same version?

mart <- biomaRt::useMart("ensembl", host="http://apr2018.archive.ensembl.org", "mmusculus_gene_ensembl")

Thank you

rna-seq kallisto
• 3,399 views • link
updated by Ram Ram • written by bioinfo • bioinfo

1 answer

Ensembl Transcript IDs and GeneIDs are stable. A transcriptID mapped to a geneID remains consistently the same across all versions. In newer versions newly identified transcripts/isoforms get assigned to a gene. Older TranscriptID's and GeneID's don't change. The key point to pay attention to is that anything that is co-ordinate based will change between different versions, e.g. genomic location of a feature (Exon, Intron, UTR etc) between 2 different versions can be different. If you are using transcriptIDs and GeneIDs, without anything to do with coordinates then you are good.

• 0 views • link
written by sure • sure

Thank you for replying. I am mapping my data using Kallisto and then using tximport and biomart to change the gene IDs. Do you think that coordinates are involved in any of the steps?

• 0 views • link
written by bioinfo • bioinfo

mm39 became the default mouse genome for Ensembl starting with v.103 (Ensembl is in version 106 as of today). I think biomaRt would retrieve coordinates for the default version which would indeed be incorrect if your data was aligned to GRCm38. Verify which database you are using. To be consistent it may be safer to use the same version of genome build in biomaRt. While the ID's may not change, their versions may have changed.

• 0 views • link
written by GenoMax GenoMax

Log in to answer this question.

More posts like this
  • GATK GRCh38 genome version for biomaRt
    written by ramiro.barrantes • ramiro.barrantes
    Hello, I am using iGenomes GATK.GRCh38 with sarek. Now I need to use the R package biomaRt to fetch the sequence and other information. What …
  • Can´t acces Biomart via R anymore
    written by Beginner • Beginner
    Hi, my R code previously accessed Biomart via: library("biomaRt") mart=useMart(biomart="ENSEMBL_MART_ENSEMBL", dataset="mmusculus_gene_ensembl", host="www.ensembl.org") resMArt=getBM(attributes=c("ensembl_gene_id","description","mgi_symbol"),mart=mart) resMArt2=getBM(attributes=c("ensembl_gene_id","gene_biotype"),mart=mart) however, mart=useMart(biomart="ENSEMBL_MART_ENSEMBL", dataset="mmusculus_gene_ensembl", host="www.ensembl.org") now gives the error: Error in listMarts(host …
  • biomaRt: Timeout on getBM().
    written by jon.klonowski jon.klonowski
    My Biomart getBM() command is timing out, and I do not know why. Failed <- getBM(attributes = c("ensembl_transcript_id", "ensembl_gene_id", "transcript_tsl"), mart = ensembl) Error in …
  • Estimated gene counts with tximport
    written by bioinfo • bioinfo
    Hello, I am analyzing bulk RNA seq data and I used Kallisto to align my data to the transcriptome. Then, I used tximport to assign …
  • Biomart mmusculus_gene_ensembl mm39
    written by ced.gob.cg • ced.gob.cg
    Hi, It seems that combining Biomart and Gviz to visualize gene annotation generates strange behavior. For mouse assembly, GRCm39/mm39 is the default now. In case …
  • Biomart mapping ensembl_transcript_id_version
    written by christoph.neu • christoph.neu
    Hello, i have a list of ensembl_transcript_id_version, and I want to map the ensembl_gene_id and external_gene_name to it. However, it does not find a hit …
  • BiomaRt: archive version of ensembl not working
    written by salamandra salamandra
    Hi, I want to access an archived version of biomart and select the human dataset. mart=useMart("ensembl_mart_92",host="http://apr2018.archive.ensembl.org") works but if I do: mart=useMart("ensembl_mart_92",host="http://apr2018.archive.ensembl.org", dataset = "hsapiens_gene_ensembl") …
  • biomaRt mmusculus_gene_ensembl dataset
    written by igor igor
    For several years, I used the following code to get mouse gene annotations using the biomaRt Bioconductor package: mart = useMart(host="useast.ensembl.org", biomart="ENSEMBL_MART_ENSEMBL", dataset="mmusculus_gene_ensembl", verbose=F) Recently, …
  • how to retrieve mouse (mm10) gene information from Ensemble using Biomart inR
    written by M K M K
    I am tying to retrieve mouse mm10 gene information using biomart library in R, but I don't know how to do that (The information that …
  • Problem Querying Archived Ensembl 73 With Biomart
    written by neelroop • neelroop
    <p>I regularly use biomaRt, the R package, to query archived versions of Ensembl. I had been querying Ensembl 65 for over a year with no …