This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Pre-miRNA strcuture presentation

Hi all, Am trying to represent the secondary structure of the pre-miRNA in the manuscript. The strcuture generated using RNAfold is not very clear. So am trying to represent it in sum other way. In few articles i came across the structure well represented. Can anyone suggest which technique i can use for presenting the secondary stem-loop structure?

sequence rna-seq
• 3,465 views • link
updated by shenwei356 shenwei356 • written by itsanju87 • itsanju87

1 answer

Similar post: A: How to get RNAfold (structure) output in text format

I had a Perl Module to output ASCII text of miRNA precursor secondary structure.

RNA::HairpinFigure - Draw hairpin-like text figure from RNA sequence and its secondary structure in dot-bracket notation.

Usage: http://search.cpan.org/dist/RNA-HairpinFigure/lib/RNA/HairpinFigure.pm#synopsis

   >hsa-mir-92a-1 MI0000093 Homo sapiens miR-92a-1 stem-loop
    CUUUCUACACAGGUUGGGAUCGGUUGCAAUGCUGUGUUUCUGUAUGGUAUUGCACUUGUCCCGGCCUGUUGAGUUUGG
    ..(((...((((((((((((.(((.(((((((((((......)))))))))))))).)))))))))))).))).....
    ---CU   UAC            C   U           UU 
         UUC   ACAGGUUGGGAU GGU GCAAUGCUGUG  U
         |||   |||||||||||| ||| |||||||||||   
         GAG   UGUCCGGCCCUG UCA CGUUAUGGUAU  G
    GGUUU   --U            U   -           GU
• 0 views • link
written by shenwei356 shenwei356

Hi ran the first command, "perl makefile.pl"

C:\RNA-HairpinFigure-0.141212>perl Makefile.PL Generating a dmake-style Makefile Writing Makefile for RNA::HairpinFigure Writing MYMETA.yml and MYMETA.json

But the next command is not running, following error is popping up

C:\RNA-HairpinFigure-0.141212>make 'make' is not recognized as an internal or external command, operable program or batch file.

• 0 views • link
written by itsanju87 • itsanju87

Install it via cpan in Windows/Linux...

cpan install RNA::HairpinFigure

For Windows, I recommend to use ActivePerl (download), so you can the "Perl Package Manager" (GUI). Search and install RNA::HairpinFigure

• 0 views • link
written by shenwei356 shenwei356

If you only have few sequences, you can also get the ASCII plot via our miRNA predicting server: FOMmiR, MicroRNA Prediction Using a Fixed-Order Markov Model Based on the Secondary Structure Pattern.

• 0 views • link
written by shenwei356 shenwei356

thank you it got installed. How to run the program with my sequence? since i have large data i cant run FOMniR.

• 0 views • link
written by itsanju87 • itsanju87

The program runs good. BUt for few sequences, the structure is not coming properly. As like below, the base information is getting in two lines, one below the other also the pipe symbol for alignment. How to sort his out?

U              C   A              -      U       G       CGU        -   UC    C
      A         AAGUUUUU
 CCUUUGAUUAGAUA AAU GAAGGAUUCUGUAU AACGGU GACAUGA UUUAUCU   UUUUAAAG AUA  UUUG A
UUUCAU UGAGAUUUA        A
 |||||||||||||| ||| |||||||||||||| |||||| ||||||| |||||||   |||||||| |||  |||| |
|||||| |||||||||
 GGGGACUGAUUUAU UUA CUUCCUAAGACAUA UUGCCA CUGUACU AAAUAGA   GAAAUUUU UAU  AAAU U
AAAGUA ACUCUAAAU        G
-              A   A              A      -       G       AUU        A   UU    U
      C         AUAAUGGU
• 0 views • link
updated by GenoMax GenoMax • written by itsanju87 • itsanju87

Can you format it by

```
content to be formated
```
• 0 views • link
written by shenwei356 shenwei356

" content to be formated", what that mean sir?

• 0 views • link
written by itsanju87 • itsanju87

itsanju87 : It may be best to post an image here, if what you see does not convey the right information. I have not tried to change it except using the (101010) button to convert those lines to code.

• 0 views • link
written by GenoMax GenoMax

Sir u want me to post an image? Its the one which is already there in the question section.

• 0 views • link
written by itsanju87 • itsanju87

Note that: As said in the description, "RNA::HairpinFigure draws hairpin-like text figure from RNA sequence and its secondary structure in dot-bracket notation. If the hairpin have multi loops, they will be deleted and treated as a big loop, the longest stem will be the final stem."

So if the stem-loop structure is not typical hairpin-like, you may use bitmap picture produced by RNAfold.

• 0 views • link
written by shenwei356 shenwei356

Log in to answer this question.

More posts like this