I have done de novo assembly of pair end raw read sequences and resulted transcripts sequences were separated based on coding potential into two categories: long non-coding RNA transcripts and coding RNA transcripts.

I want to find how many of long non coding can make hairpin loop structure and how can I annotate these miRNA which are generated from long non coding RNA in plants.

  • 1
    $\begingroup$ Welcome to Bioinformatics! Is there anything you tried before posting? Is there any reason why you tagged the question with Bioconductor (Maybe you want to restrict to use a Bioconductor package, or have used one and you can't manage to make it classify your lnRNA)? Also, I think there are two questions here: how to find if a lncRNA can make hairpins and annotate RNA to miRNA. It would be better to focus in just one question (perhaps the first one and later on the second one). $\endgroup$ – llrs Jun 29 '17 at 9:08

For the first question: "how many of long non coding can make hairpin loop structure?"

The easiest thing to would be to run your sequences through an RNA secondary structure prediction tool. There's a few tools for doing this but the most commonly used are RNAfold from the ViennaRNA package and mfold.

This will give you an output that looks like this, where the parentheses indicate which bases are paired with each other.

GGGCUAUUAGCUCAGUUGGUUAGAGCGCACCCCUGAUAAGGGUGAGGUCGCUGAUUCGAAUUCAGCAUAGCCCA (((((((..((((.........)))).(((((.......))))).....(((((.......)))))))))))).

To see if there's hairpins, you'll need to annotate the elements within it. Take a look at this question on biostars for an example of how to do this. In short, you can use the forgi library to produce an annotation like this:

(((((((..((((.........)))).(((((.......))))).....(((((.......)))))))))))). sssssssmmsssshhhhhhhhhssssmssssshhhhhhhsssssmmmmmssssshhhhhhhsssssssssssse

If there's an h in the output annotation, you have a hairpin. The remainder of your question about annotating miRNAs will be left to somebody else.

  • 2
    $\begingroup$ hairpins are quite easy to identify using the bracket notation; just look for brackets that are in a left-right orientation with nothing in between (e.g. (.), (..), (...)). $\endgroup$ – gringer Jun 29 '17 at 23:48

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.