This is a test version of Biostars. For the public version, visit https://www.biostars.org.
classification of small RNA-seq

How to categorize smallRNA reads according to different subtype such as miRNA,snoRNA,rRNA, tRNA, snRNA ?

Also same question is here Mapping microRNA reads to Genome and non-coding RNAs

Thanks

rna-seq alignment
• 3,384 views • link
written by Sam Sam

Please elaborate. If a question fits in one sentence you usually didn't explain it well enough.

• 0 views • link
written by WouterDeCoster WouterDeCoster

1 answer

Although your question seems lacking some depth and clarity, I would just point out one way (there are many ways, e.g. prediction, etc.): take the gene biotype information of your's organism Ensembl gtf file (i.e. count your reads using this file). There you can see different biotypes, including what you mentioned (although it's kind of weird if you want to get lncRNA from small-RNA reads).

• 0 views • link
written by aka001 aka001

Currently you are at which step (another way to put it: how did you process your reads)?

• 0 views • link
written by aka001 aka001

I performed adaptor removal and aligned filtered reads against to human miRNAs in miRBase using sRNAWorkbench

• 0 views • link
written by Sam Sam

.....the remaining unique sequences were mapped to the radish reference genome which consisted of radish GSS, EST and transcriptome sequences, to analyze the expression and distribution of sRNAs on genome using SOAP2 program [25,26]. Perfectly matched sequences were retained for following analysis. By querying against the NCBI Genbank (http://www.ncbi.nlm.nih.gov/genbank/) and Rfam (10.1) (http://www.sanger.ac.uk/resources/databases/rfam.html) databases, the sRNA sequences matching rRNA, tRNA, snRNA, snoRNA as well as sequences containing poly (A) tails were excluded. The remaining unique sequences were aligned against miRBase 20.0 (http://www.mirbase.org/index.shtml) to identify radish known miRNAs.....

I think you want to do something like this, correct? https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-015-1416-5

• 0 views • link
written by lessismore lessismore

yes exactly, are you agree with local NCBI blast against Rfam database and then grep according to each subtype?

• 0 views • link
written by Sam Sam

Log in to answer this question.

More posts like this