Retrieve genes with more than N paralogs on ensembl
Is there a way on ensembl (or biomart?) to retrieve a list of gene IDs from a species that have more than N paralogs in this species?
I tried with biomart online but I didn't find any filters adapted.
• 1,703 views
•
link
updated
by
Jean-Karim Heriche
•
written
by
corend
•
written
by
corend
1 answer
You'll have to use the perl API for this. Something like this:
my @homologs = @{$Ensgene->get_all_homologous_Genes()};
my @paralogs = grep { $_->[2] eq 'homo_sapiens' && $_->[1]->description=~/paralog/ } @homologs;
• 0 views
•
link
written
by
Jean-Karim Heriche
Log in to answer this question.
More posts like this
-
Mapping SNPs to genes
written by Nigussie •
-
How to assign gene names to ensembl IDs if species not in Biomart
written by bioinfo •
-
Retrieve only protein coding esnsemble gene ids and gene symbols
written by iibrams07 •
-
Can't convert dog ensembl IDs into gene names
written by fifty_fifty •
-
getGeneLengthAndGCContent in EDASeq complains about Ensemble gene ID
written by An Ignorant Wanderer •
-
Converting Ensembl ID to NCBI ID
written by botloggy •
-
Ensembl IDs as input in ReactomePA package in R
written by sainikeshav6 •
-
Unable to query specific ensembl ID through biomaRt (R) or mygene (Python)
written by JourneyToAbyss
-
From Ensembl protein ID to sequence
written by Joseph Hughes
-
Downloading All Gene Families From Mammalian Species In Ensembl
written by Marki •
