Biomart select rRNA sequences
2
0
Entering edit mode
2.5 years ago
Vlad ▴ 10

I want to select all the ribossomal genes from mouse. I tried to use Rat Genome database rgd.mcw.edu) and select the ribossomal genes but some of the genes that are listed in my Sequencing results are not present there (such as Rps28 gene and some more). I am trying now to select from Biomart but I cannot find any attribute that can classify my Ensembl_id into for example mRNA, rRNA etc.

Is there a way to get all rRNA genes from biomart or is there a better library to achieve my goal?

Biomart rRNA • 1.4k views
ADD COMMENT
1
Entering edit mode
2.5 years ago
Mike Smith ★ 1.7k

You can use the 'Gene Type' filter in Ensembl BioMart to restrict your results to a certain group. With the biomaRt package you'd do that via:

library(biomaRt)
mart <- useEnsembl('ensembl', dataset = 'rnorvegicus_gene_ensembl')
biomaRt::getBM(values="rRNA", 
               filters="biotype", 
               attributes=c("ensembl_gene_id", "external_gene_name", "gene_biotype"), 
               mart = mart)

However I don't see Rps28 listed as a rat gene on Ensembl, so it may be that you are encountering annotation & naming difference between different sources.

ADD COMMENT
0
Entering edit mode

It seems Ensembl is not a complete annotation of rRNA genes. There are only 5S and 5.8S in Human Ensembl. Correct me if I am wrong. Thx

ADD REPLY
0
Entering edit mode
2.5 years ago
Vlad ▴ 10

Works perfectly , thanks.

ADD COMMENT

Login before adding your answer.

Traffic: 2444 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