Question: miRNAs detected in bulk mRNA-Seq data - how to annotate?
0
gravatar for Kristin Muench
8 months ago by
United States
Kristin Muench480 wrote:

Hello,

In our bulk mRNA-Seq data, about ~600 of our ~21,000 detected genes were miRNAs. All of these fall within the bounds of expression of non-miRNA genes, and about ~20 miRNAs fall within the upper half of gene expression in the dataset.

I was surprised by this because I thought most miRNAs would be removed via polyA selection. Also, we did NOT use a kit for small RNA capture and sequencing.

I'm wondering if these reads are aligning to a pre-miRNA that are longer than 75 basepairs. I'd like it if I could take my list of miRNAs in R and combine them with a database that has information about how long the pre- and post-processed miRNAs are, to sanity check my theory.

Does such an miRNA database exist? I'm trying to use mirbase.db, but I'm confused about its use, as well as whether it has the information I'm looking for: http://bioconductor.org/packages/release/data/annotation/html/mirbase.db.html

Thanks!

mirna rna-seq R • 261 views
ADD COMMENTlink modified 8 months ago • written 8 months ago by Kristin Muench480
1

About poly-A:

Human microRNAs are processed from capped, polyadenylated transcripts that can also function as mRNAs

ADD REPLYlink written 8 months ago by h.mon29k

This is a great point. Thank you

ADD REPLYlink written 8 months ago by Kristin Muench480
1

I suggest you try RNACentral. Filter according to your needs (organism, type of RNA).

ADD REPLYlink modified 8 months ago • written 8 months ago by genomax78k

I'm a little confused...was your prep designed to catch miRNAs or not?

ADD REPLYlink written 8 months ago by swbarnes27.5k

It was not designed to catch miRNAs, thus my confusion. It was a KAPA stranded mRNA-Seq kit: https://www.kapabiosystems.com/product-applications/products/next-generation-sequencing-2/rna-library-preparation-2/kapa-stranded-mrna-seq-kits/

ADD REPLYlink modified 8 months ago • written 8 months ago by Kristin Muench480

You wrote "bulk miRNA-Seq" but I don't think you meant that. I think you need to scrutinize the library protocol, or talk to whoever prepped it, because I strongly suspect that this library prep is supposed to filter away small fragments, so you can't treat the things you think are short fragments as legitimate.

ADD REPLYlink written 8 months ago by swbarnes27.5k

Thanks for catching that typo, fixed above. I'm wondering that too, whether pri-miRNAs are even a possibility in an mRNA-Seq prep.

ADD REPLYlink written 8 months ago by Kristin Muench480
2
gravatar for Asaf
8 months ago by
Asaf7.0k
Israel
Asaf7.0k wrote:

mirbase.org -> Download -> hairpin.fa will give you all the sequences of pre-processed miRNAs, you can then filter them according to your organism and map the reads against the resulting database.

ADD COMMENTlink written 8 months ago by Asaf7.0k
2

Be on the lookout for U in the sequence if you get data from MiRBase. You will need to change those to T before aligning.

ADD REPLYlink modified 8 months ago • written 8 months ago by genomax78k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1338 users visited in the last hour