Welcome guest, Sign in



This protocol describes small RNA library preparation from Vigna mungo total RNA followed by deep sequencing and analysis for microRNA identification.​

Thanks for your further question/comment. It has been sent to the author(s) of this protocol. You will receive a notification once your question/comment is addressed again by the author(s).
Meanwhile, it would be great if you could help us to spread the word about Bio-protocol.


Mungbean Yellow Mosaic India Virus (MYMIV)-infection, Small RNA Library Construction and Deep Sequencing for MicroRNA Identification in Vigna mungo

Plant Science > Plant molecular biology > RNA > RNA sequencing
Authors: Anirban Kundu
Anirban KunduAffiliation: Division of Plant Biology, Bose Institute, Kolkata, West Bengal, India
Bio-protocol author page: a2047
Sujay Paul
Sujay PaulAffiliation: Division of Plant Biology, Bose Institute, Kolkata, West Bengal, India
Bio-protocol author page: a2046
Amita Pal
Amita PalAffiliation: Division of Plant Biology, Bose Institute, Kolkata, West Bengal, India
For correspondence: amita@jcbose.ac.in
Bio-protocol author page: a2048
 and Genotypic Technology
Genotypic TechnologyAffiliation: Genotypic Technology Private Limited, #2/13, Balaji Complex, Bangalore, India
Bio-protocol author page: a3589
Vol 6, Iss 20, 10/20/2016, 366 views, 0 Q&A
DOI: https://doi.org/10.21769/BioProtoc.1961

[Abstract] This protocol describes small RNA library preparation from Vigna mungo total RNA followed by deep sequencing and analysis for microRNA identification.​
Keywords: MYMIV, Vigna mungo, MicroRNA, Deep sequencing, Regards

Materials and Reagents

  1. Filter papers
  2. Plastic pots
  3. Young Vigna mungo plants
  4. Whiteflies (Bemisia tabaci)
  5. HgCl2
  6. TruSeq small RNA Library Preparation Kit (Illumina, catalog number: RS-200-0012 )
  7. Trizol reagent (Invitrogen, USA)
  8. Isopropanol
  9. Acrylamide/bisacrylamide (Merck Millipore, catalog number: 623100281001730 )
  10. T4 RNA ligase 2, truncated (New England Biolabs, catalog number: M0242 )
  11. Tris-HCl pH 7.5
  12. MgCl2
  13. DTT
  14. ATP
  15. PEG8000
  16. SuperScript® II reverse transcriptase (Thermo Fisher Scientific, InvitrogenTM, catalog number: 18064014 )
  17. TEMED (Merck Millipore, catalog number: 623171280051730 )
  18. Ammonium persulphate (Merck Millipore, catalog number: 623171100101730 )
  19. 10x TBE buffer (Thermo Fisher Scientific, AmbionTM, catalog number: AM9865 )
  20. UltraPureTM glycogen (Thermo Fisher Scientific, InvitrogenTM, catalog number: 10814-010 )
  21. 3 M sodium acetate (Sigma-Aldrich, catalog number: S7899-100ml )
  22. Absolute ethanol (EMD Millipore, catalog number: 100983 )
  23. Qubit® dsDNA HS Assay Kit (Thermo Fisher Scientific, InvitrogenTM, catalog number: Q32854 )
  24. High Sensitivity DNA Kit (Agilent Technologies, catalog number: 5067-4626 )
  25. T4 RNA Ligase 2 buffer (see Recipes)


  1. Electrophoresis apparatus (Bio-Rad Laboratories, model: Mini-Sub® GT Cell )
  2. Bioanalyzer (Agilent Technologies, model: Agilent 2100 )
  3. Thermocycler, DNA engine cycler (Bio-Rad Laboratories, model: PTC-0200G )
  4. Qubit 3.0 fluorometer (Life Technologies)
  5. Illumina genome analyzer IIx (Illumina)


  1. Cutadapt-0.9.3
  2. Bowtie-0.12.7
  3. SeqQCv2.1
  4. BLAST (Rfam and miRBase)


  1. Plant material and growth conditions
    1. Use young V. mungo plants for this experiment. Surface sterilize mature seeds (0.1% HgCl2 for 10 min) and rinse twice in double-distilled water and keep in moistened filter papers for germination at 28 ± 2 °C, 16/8 h light/dark photoperiod and 70% relative humidity for 1 day. Perform three replicates of this experiment.
    2. Sown pre-germinated seeds in plastic pots and grow under controlled greenhouse conditions (temperature 26 ± 2 °C, relative humidity 70% and photoperiod 16 h day and 8 h night) in Soilrite soil mix (mixture of exfoliated vermiculite, Irish peat moss and perlite in the ratio 1/3:1/3:1/3).
    3. MYMIV-stress was given to plants with similar growth when their first trifoliate leaf expanded completely (approx. 21 days after sprouting).

      Figure1. Yellow mosaic symptoms on Vigna leaves

  2. MYMIV-stress treatments
    1. Rear white fly (Bemisia tabaci) populations on susceptible black gram plants grown in insect proof cages maintained in the insectory facility at the experimental farm, Madhyamgram, Bose Institute.
    2. Capture about 25-30 adult white flies in transparent glass trappers and allowed for a 24 h acquisition access period (AAP) on symptomatic leaves of naturally infected black gram plants (Figure 1).
    3. Give a 24 h inoculation access period (IAP) in the glass trappers with viruliferous whiteflies on the leaves of healthy plants, thereby transmitting the virus particles from the viruliferous flies to the plant (Figure 2).

      Figure 2. Trapper method for MYMIV inoculation

    4. Confirm the infection by PCR amplification of the MYMIV coat protein (575 bp fragment; Gen-Bank Accession number: HQ221570) using the primers MYMIV- coat protein (CP) gene (CP-F: 5’-GAA ACCTCGGTTTTACCG ACTGTATA G-3’ and CP-R: 5’-TTGCATA CACAGGATTTG AGGCATGAG-3’).

  3. Small RNA library preparation
    1. Construct small RNA libraries for sequencing according to the Illumina TruSeq small RNA library protocol outlined in 'TruSeq small RNA sample preparation guide' (Part # 15004197; Rev. E; February 2013) (Figure 3).
    2. Manually isolate10 µg of total RNA (with a RIN value of 8 and above) using Trizol reagent and concentrate by 2 volumes chilled isopropanol, then resolve in a 15% denaturing polyacrylamide gel, and then excise small RNA fragments in the range of 18-0 nt from the gel and purify. Ligate 3’ and 5’ adaptors to the small RNAs using T4 RNA Ligase 2 in a buffer (see Recipe 1) by incubating at 25 °C for 1 h.
    3. Reverse transcribe the ligated products with Superscript II reverse transcriptase by priming with adapter specific primers in 15 µl reaction mixture containing 0.15 µl dNTP mix (100 mM), 1 µl reverse transcriptase (50 U/µl), 10x buffer (1.5 µl), RNase inhibitor (0.2 µl), water (4.2 µl), 3 µl RT primer and 5 µl RNA. Incubated the reaction mixture at 42 °C for 1 h followed by a heat inactivation of 5 min at 85 °C.
    4. Enriched the cDNA by a low cycle PCR amplification of 15 cycles (initial denaturation of 30 sec at 98 °C; 15 cycles of 98 °C for 10 sec, 60 °C for 30 sec, and 72 °C for 15 sec; followed by a final extension at 72 °C for 10 min) to generate ample template for deep sequencing and clean using polyacrylamide gel followed by overnight gel elution and salt precipitation using glycogen, 3 M sodium acetate and absolute ethanol. Re-suspended the precipitate in nuclease free water.
    5. Quantify the prepared library using Qubit fluorometer (50 ng/µl), and validate for quality by running an aliquot on high sensitivity DNA kits.

      Figure3. Flow chart for TruSeq small RNA sample preparation

    6. Analyse small RNA using softwares/programmes like cutadapt-0.9.3 (Martin, 2011), bowtie-0.12.7 (Langmead and Salzberg, 2012), SeqQCv2.1 (DeLuca et al., 2012) and BLAST (Altschul et al., 1990). Steps of analysis are as under:
      1. Sequencing samples read files
        Generate fastq files from sequencing on GAIIx and use for downstream analysis.
      2. Sequence reads quality control (QC) reports
        Perform QC for all the generated sequence reads from different samples in different lanes (Mugasimangalam et al., manuscript under preparation). Consider QC report for all the sequencing reads on the basis of percentage of high quality reads which has > 70% of bases in a read with > 20 phred score. Based on the analysis as per the QC (HQ reads > 90%) report process the sequenced reads further.
      3. Trimming reads
        Trimming is done using cutadapt-0.9.3 to remove ligated adapter sequences from the original small RNA sequences. Cutadapt tool removes user specified sequences from the fastq files (Martin, 2011).
      4. Generating unique reads and counting
        Identify unique reads by writing the custom scripts and generate the read counts. Briefly below mentioned steps were carried out:
        Remove the reads which were less than or equal to 15 bases in length. Consider reads whose lengths were in the range of 16 to 36 bases. Identify unique sequences and remove duplicate sequences. Generate read count for each unique reads by identifying how many times it is present in the sequence data.
      5. Annotations of SnoRNA, SnRNA, rRNA, tRNA
        Align the sequences against Glycine max genome using bowtie-0.12.7. Bowtie supports end-to-end alignment of fastq reads to the reference fasta file (Langmead et al., 2009). End-to-end alignment considers only those reads which align completely with upto 3 mismatches to the reference fasta file. Check the reads aligning with upto 3 mismatches by BLASTn (Altschul et al., 1990) against Rfam database (using default parameters) to identify other small RNA in the sequenced data (Meyers et al., 2008).
      6. miRNA annotation
        Check the unique sequences by BLASTn (with an E-value cutoff of 10) against miRbase database to identify known miRNAs. Blast identifies similar sequences based on the sequence homology between the reference and query (unique reads).


  1. T4 RNA Ligase 2 buffer
    31 mM Tris-HCl pH 7.5
    6.3 mM MgCl2
    0.63 mM DTT
    0.2 mM ATP
    5.2% PEG8000


The original version of this protocol was described in Paul et al. (2014). This work was supported by the Council of Scientific and Industrial Research, New-Delhi, India for the Emeritus Scientist’s Project [Sanction No. 21 (0884)/12/EMR-II].


  1. Altschul, S. F., Gish, W., Miller, W., Myers, E. W. and Lipman, D. J. (1990). Basic local alignment search tool. J Mol Biol 215(3): 403-410.
  2. DeLuca, D. S., Levin, J. Z., Sivachenko, A., Fennell, T., Nazaire, M. D., Williams, C., Reich, M., Winckler, W. and Getz, G. (2012). RNA-SeQC: RNA-seq metrics for quality control and process optimization. Bioinformatics 28(11): 1530-1532.
  3. Langmead, B., Trapnell, C., Pop, M. and Salzberg, S. L. (2009). Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol 10(3): R25.
  4. Langmead, B. and Salzberg, S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat Methods 9(4): 357-359.
  5. Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMB J 17(1): pp. 10-12.
  6. Meyers, B. C., Axtell, M. J., Bartel, B., Bartel, D. P., Baulcombe, D., Bowman, J. L., Cao, X., Carrington, J. C., Chen, X., Green, P. J., Griffiths-Jones, S., Jacobsen, S. E., Mallory, A. C., Martienssen, R. A., Poethig, R. S., Qi, Y., Vaucheret, H., Voinnet, O., Watanabe, Y., Weigel, D. and Zhu, J. K. (2008). Criteria for annotation of plant MicroRNAs. Plant Cell 20(12): 3186-3190.
  7. Mugasimangalam, R., Mohan, K., Jayakumar, V., Ashick, M. M., Kapila, G. and Niranjan, V. SeqQC: rapid quality control of next generation sequence data (manuscript under preparation).
  8. Paul, S., Kundu, A. and Pal, A. (2014). Identification and expression profiling of Vigna mungo microRNAs from leaf small RNA transcriptome by deep sequencing. J Integr Plant Biol 56(1): 15-23.
  9. TruSeqsmall RNA sample preparation guide. Illumina Inc.

How to cite: Kundu, A., Paul, S., Pal, A. and Technology, G. (2016). Mungbean Yellow Mosaic India Virus (MYMIV)-infection, Small RNA Library Construction and Deep Sequencing for MicroRNA Identification in Vigna mungo. Bio-protocol 6(20): e1961. DOI: 10.21769/BioProtoc.1961; Full Text

Share Your Feedback:

  • Add Photo
  • Add Video

Bio-protocol's major goal is to make reproducing an experiment an easier task. If you have used this protocol, it would be great if you could share your experience by leaving some comments, uploading images or even sharing some videos. Please login to post your feedback.

Q&A and Troubleshooting:

  • Add Photo
  • Add Video

Please login to post your questions/comments. Your questions will be directed to the authors of the protocol. The authors will be requested to answer your questions at their earliest convenience. Once your questions are answered, you will be informed using the email address that you register with bio-protocol.
You are highly recommended to post your data (images or even videos) for the troubleshooting. For uploading videos, you may need a Google account because Bio-protocol uses YouTube to host videos.

Login | Register
How to cite
Twitter Twitter
LinkedIn LinkedIn
Google+ Google+
Facebook Facebook