Skip to content

Instantly share code, notes, and snippets.

@slavailn
Created May 17, 2021 05:00
Show Gist options
  • Star 1 You must be signed in to star a gist
  • Fork 1 You must be signed in to fork a gist
  • Save slavailn/ab27b425a7a190b1e26752d9a3821b74 to your computer and use it in GitHub Desktop.
Save slavailn/ab27b425a7a190b1e26752d9a3821b74 to your computer and use it in GitHub Desktop.
Get paralogs for the list of genes using biomaRt
# Get paralogs for the list of genes
library(biomaRt)
human <- useMart("ensembl", dataset = "hsapiens_gene_ensembl")
gene_id <- c("TPM1", "BOD1", "ADAP1")
results <- getBM(attributes = c("ensembl_gene_id",
"external_gene_name",
"hsapiens_paralog_ensembl_gene",
"hsapiens_paralog_associated_gene_name"),
filters = "external_gene_name",
values = gene_id,
mart = human)
results
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment