This is a test version of Biostars. For the public version, visit https://www.biostars.org.
GATK GRCh38 genome version for biomaRt

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 ensembl version does it correspond to? For example, for ncbi37 I would do:

 genome = useMart(biomart="ENSEMBL_MART_ENSEMBL",
                               host="grch37.ensembl.org",
                               path="/biomart/martservice",
                               dataset="hsapiens_gene_ensembl")

and for TCGA I do:

  grch38 = useMart(biomart="ENSEMBL_MART_ENSEMBL",
                   host="https://nov2020.archive.ensembl.org",
                   dataset="hsapiens_gene_ensembl")

In my case now it would then be for the GATK bundle genome version for GRCh38. Would it be as simple as:

  grch38 = useMart(biomart="ENSEMBL_MART_ENSEMBL",
                   dataset="hsapiens_gene_ensembl")

Thank you in advance.

biomart
• 1,184 views • link
updated by GenoMax GenoMax • written by ramiro.barrantes • ramiro.barrantes

1 answer

Current Ensembl release is GRCh38. See the genome versions and the relationship to Ensembl releases in this table: https://www.ensembl.org/info/website/archives/assembly.html

Difference between ENSEMBL releases is a past thread that offers details.

That said, if you need to get sequences from GATK's version of GRCh38 then you should use the sequence they provide in their resource bundle here: https://console.cloud.google.com/storage/browser/genomics-public-data/resources/broad/hg38/v0/ You can get sequences you need using samtools faidx (or similar tools) with the included indexed fasta file.

• 0 views • link
written by GenoMax GenoMax

Log in to answer this question.

More posts like this
  • Correct genome version for ncbi37
    written by ramiro.barrantes • ramiro.barrantes
    I have a paper that aligned data using NCBI build 37 (PMID 24429703), it said: > Paired-end sequencing reads were aligned to the human genome …
  • 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: 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") …
  • dbSNP build in ENSEMBL biomaRt?
    written by 09meh • 09meh
    I'm wondering how to determine what dbSNP build (i.e. 150 or 151) the variation annotations are based on. For example, I can find that for …
  • Ensembl biomaRt down
    written by sergio.artsc • sergio.artsc
    I am trying to get some queries from biomaRt package. I am defining the next connection: snp_db <- useMart("ENSEMBL_MART_SNP", dataset="hsapiens_snp", host="grch37.ensembl.org") and I get the …
  • BioMart: connection fails
    written by salamandra salamandra
    Hi, running this > ensembl = useMart(biomart > ="ENSEMBL_MART_ENSEMBL",dataset="hsapiens_gene_ensembl", host="may2017.archive.ensembl.org") or this: > ensembl = useMart(biomart > ="ENSEMBL_MART_ENSEMBL",dataset="hsapiens_gene_ensembl", host="http://dec2017.archive.ensembl.org") used to work in the past. …
  • using biomaRt in R to annotate SNPs to Genes within 1000kb
    written by Sheila Sheila
    Hi! In short, I have a list of SNPs that I would like to map to the closest gene within 1000kb. I am using the …
  • biomaRt GRCh37 down?
    written by alexander.nash • alexander.nash
    Hi, I am currently getting the following error: library(biomaRt) mart<-useMart(biomart="ENSEMBL_MART_ENSEMBL", dataset="hsapiens_gene_ensembl", host="grch37.ensembl.org") Extra content at the end of the document Error: 1: Extra content at …
  • Map ensembl gene id form hg38 in hg19
    written by Tanvir Ahamed Tanvir Ahamed
    I want to map all ensembl gene id form hg38 in hg19. Any help will be appreciated ? Thanks !! # Example : Loading library …
  • biomaRt GRCh37 down?
    written by andrew.j.skelton73 andrew.j.skelton73
    Hi, Using the connection string: ``` grch37 <- useMart(biomart = "ENSEMBL_MART_ENSEMBL", host = "grch37.ensembl.org", path = "/biomart/martservice", dataset = "hsapiens_gene_ensembl") ``` Doesn't appear to work, …