Ensembl phylogenetic context
1
0
Entering edit mode
2.8 years ago
gokberk ▴ 90

Hi, I was wondering if there is a way to download or generate 24 primates EPO-Extended alignments for a set of SNPs using Ensembl VEP. I can download it from the phylogenetic context section in Ensembl for a given SNP, but was wondering how to do it for a large number of SNPs from command line.

Thanks!

alignment ensembl phylogeny • 1.9k views
ADD COMMENT
3
Entering edit mode
2.8 years ago
Emily 23k

There isn't a method in the VEP or any of the VEP plugins that do this, so you would need to write something yourself. You could invoke the alignment REST API endpoint, to fetch the alignment ± whatever up/downstream distance you want.

ADD COMMENT
0
Entering edit mode

Got it, thank you!

ADD REPLY
0
Entering edit mode

Hi again Emily_Ensembl , I have quick follow-up question. All my SNPs are annotated in GRCh37 coordinates. Is there a way to get GRCh37 sequence using the code you pointed to?

I've been using the code chunk below:

server <- "http://rest.ensembl.org"
ext <- paste0("/alignment/region/homo_sapiens/", tmp_chr, ":", tmp_seq_str, "-", 
              tmp_seq_end, "?species_set_group=primates")

r <- GET(paste(server, ext, sep = ""), content_type("application/json"))

stop_for_status(r)

I tried my chance with /alignment/region/homo_sapiens/GRCh37/, but didn't manage to make it work. I was wondering if there is a simple way or if I have to CrossMap all my SNPs.

Thanks!

ADD REPLY
1
Entering edit mode

server <- "http://grch37.rest.ensembl.org"

ADD REPLY
0
Entering edit mode

Many thanks!!

ADD REPLY
0
Entering edit mode

I have just tried this, but it doesn't seem to be working for GRCh37. I received the following error: "No method_link_specices_set found for method LASTZ_NET and species_set_group primates "

ADD REPLY
1
Entering edit mode

Ah, yes, I forgot. There's no comparative genomics data available on GRCh37. It got out of date so we removed it.

ADD REPLY
0
Entering edit mode

Oh I see, I will lift over my SNPs to GRCh38 then, thanks a lot for all your help! :)

ADD REPLY

Login before adding your answer.

Traffic: 1077 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6