I have a list of about 80 HGNC symbols.
I want to find their descriptions and function(s) using biomart
So - first approach is to find everything using biomart. Here I provide an example for three genes.
gL <- c("LTC4S", "ALOX5", "NAT2") library(biomaRt) mart <- useMart(biomart = "ensembl", dataset = "hsapiens_gene_ensembl") results <- getBM(attributes = c("go_id", "hgnc_symbol"), filters = "hgnc_symbol", values = gL, mart = mart)
For 3 HGNCs we have 46 GO ids... With 80 HGNCs - well you can imagine - over 1500 GO categories.
However I'm only interested in the leaf terms, i.e. the most specific terms, for each branch of the GO with which the gene is annotated.
Is there any easy way to get this using biomaRt or other R tools like GO.db?