欢迎您, 登录 | 注册

首页 | English

X
加载中

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.

X

Mungbean Yellow Mosaic India Virus (MYMIV)-infection, Small RNA Library Construction and Deep Sequencing for MicroRNA Identification in Vigna mungo
黑吉豆中的绿豆黄化花叶印度病毒(MYMIV)感染、sRNA库建立和鉴定microRNA的深度测序

植物科学 > 植物分子生物学 > RNA > RNA 测序
作者: 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, 217 views, 0 Q&A, How to cite
DOI: http://dx.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(MYMIV), Vigna mungo(黑吉豆), MicroRNA(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)

Equipment

  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)

Software

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

Procedure

  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).

Recipes

  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

Acknowledgments

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].

References

  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 this protocol: 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



可重复性反馈:

  • 添加图片
  • 添加视频

我们的目标是让重复别人的实验变得更轻松,如果您已经使用过本实验方案,欢迎您做出评价。我们鼓励上传实验图片或视频与小伙伴们(同行)分享您的实验心得和经验。(评论前请登录)

问题&解答:

  • 添加图片
  • 添加视频

(提问前,请先登陆)bio-protocol作为媒介平台,会将您的问题转发给作者,并将作者的回复发送至您的邮箱(在bio-protocol注册时所用的邮箱)。为了作者与用户间沟通流畅(作者能准确理解您所遇到的问题并给与正确的建议),我们鼓励用户用图片或者视频的形式来说明遇到的问题。由于本平台用Youtube储存、播放视频,作者需要google 账户来上传视频。


登陆 | 注册
分享
Twitter Twitter
LinkedIn LinkedIn
Google+ Google+
Facebook Facebook