Improved Placement of Multi-mapping Small RNAs

Abstract

High-throughput sequencing of small RNAs (sRNA-seq) is a popular method used to discover and annotate microRNAs (miRNAs), endogenous short interfering RNAs (siRNAs), and Piwi-associated RNAs (piRNAs). One of the key steps in sRNA-seq data analysis is alignment to a reference genome. sRNA-seq libraries often have a high proportion of reads that align to multiple genomic locations, which makes determining their true origins difficult. Commonly used sRNA-seq alignment methods result in either very low precision (choosing an alignment at random), or sensitivity (ignoring multi-mapping reads). Here, we describe and test an sRNA-seq alignment strategy that uses local genomic context to guide decisions on proper placements of multi-mapped sRNA-seq reads. Tests using simulated sRNA-seq data demonstrated that this local-weighting method outperforms other alignment strategies using three different plant genomes. Experimental analyses with real sRNA-seq data also indicate superior performance of local-weighting methods for both plant miRNAs and heterochromatic siRNAs. The local-weighting methods we have developed are implemented as part of the sRNA-seq analysis program ShortStack, which is freely available under a general public license. Improved genome alignments of sRNA-seq data should increase the quality of downstream analyses and genome annotation efforts.

Shown for reading. Its publisher allows quoting it, not republishing it.

References 27

1. 1.

   Allen E. Xie Z. Gustafson A. M. Carrington J. C., 2005 microRNA-directed phasing during trans-acting siRNA biogenesis in plants. Cell 121: 207–221.

   Not in our library10.1016/j.cell.2005.04.004

2. 2.

   Aravin A. Gaidatzis D. Pfeffer S. Lagos-Quintana M. Landgraf P., 2006 A novel class of small RNAs bind to MILI protein in mouse testes. Nature 442: 203–207.

   Not in our library10.1038/nature04916

3. 3.

   Axtell M. J., 2013 ShortStack: comprehensive annotation and quantification of small RNA genes. RNA 19: 740–751.

   Not in our library10.1261/rna.035279.112

4. 4.
   Identification of Pol IV and RDR2-dependent precursors of 24 nt siRNAs guiding de novo DNA methylation in Arabidopsis · 2015

   No notice
5. 5.
   Relative quantification of mRNA: comparison of methods currently used for real-time PCR data analysis · 2007

   No notice
6. 6.

   Coruh C. Shahid S. Axtell M. J., 2014 Seeing the forest for the trees: annotating small RNA producing genes in plants. Curr. Opin. Plant Biol. 18: 87–95.

   Not in our library10.1016/j.pbi.2014.02.008

7. 7.

   Cuperus J. T. Fahlgren N. Carrington J. C., 2011 Evolution and functional diversification of MIRNA genes. Plant Cell 23: 431–442.

   Not in our library10.1105/tpc.110.082784

8. 8.

   Faulkner G. J. Forrest A. R. R. Chalk A. M. Schroder K. Hayashizaki Y., 2008 A rescue strategy for multimapping short sequence tags refines surveys of transcriptional activity by CAGE. Genomics 91: 281–288.

   Not in our library10.1016/j.ygeno.2007.11.003

9. 9.
   TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions · 2013

   No notice
10. 10.
    miRBase: annotating high confidence microRNAs using deep sequencing data · 2013

    No notice
11. 11.
    Ultrafast and memory-efficient alignment of short DNA sequences to the human genome · 2009

    No notice
12. 12.

    Li H. Durbin R., 2009 Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25: 1754–1760.

    Not in our library10.1093/bioinformatics/btp324

13. 13.
    Detection of Pol IV/RDR2-dependent transcripts at the genomic scale in Arabidopsis reveals features and regulation of siRNA biogenesis · 2014

    No notice
14. 14.

    Malone C. D. Brennecke J. Dus M. Stark A. McCombie W. R., 2009 Specialized piRNA pathways act in germline and somatic tissues of the Drosophila ovary. Cell 137: 522–535.

    Not in our library10.1016/j.cell.2009.03.040

15. 15.

    Martin M., 2011 Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.journal 17: 10–12.

    Not in our library

16. 16.

    Matzke M. A. Mosher R. A., 2014 RNA-directed DNA methylation: an epigenetic pathway of increasing complexity. Nat. Rev. Genet. 15: 394–408.

    Not in our library10.1038/nrg3683

17. 17.

    Mortazavi A. Williams B. A. McCue K. Schaeffer L. Wold B., 2008 Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat. Methods 5: 621–628.

    Not in our library10.1038/nmeth.1226

18. 18.

    Moxon S. Schwach F. Dalmay T. Maclean D. Studholme D. J., 2008 A toolkit for analysing large-scale plant small RNA datasets. Bioinformatics 24: 2252–2253.

    Not in our library10.1093/bioinformatics/btn428

19. 19.

    Okamura K. Chung W.-J. Ruby J. G. Guo H. Bartel D. P., 2008 The Drosophila hairpin RNA pathway generates endogenous short interfering RNAs. Nature 453: 803–806.

    Not in our library10.1038/nature07015

20. 20.

    Ramakers C. Ruijter J. M. Deprez R. H. L. Moorman A. F. M., 2003 Assumption-free analysis of quantitative real-time polymerase chain reaction (PCR) data. Neurosci. Lett. 339: 62–66.

    Not in our library10.1016/s0304-3940(02)01423-4

21. 21.

    Schmid M. W. Grossniklaus U., 2015 Rcount: simple and flexible RNA-Seq read counting. Bioinformatics 31: 436–437.

    Not in our library10.1093/bioinformatics/btu680

22. 22.

    Stocks M. B. Moxon S. Mapleson D. Woolfenden H. C. Mohorianu I., 2012 The UEA sRNA workbench: a suite of tools for analysing and visualizing next generation sequencing microRNA and small RNA datasets. Bioinformatics 28: 2059–2061.

    Not in our library10.1093/bioinformatics/bts311

23. 23.

    Taylor R. S. Tarver J. E. Hiscock S. J. Donoghue P. C. J., 2014 Evolutionary history of plant microRNAs. Trends Plant Sci. 19: 175–182.

    Not in our library10.1016/j.tplants.2013.11.008

24. 24.

    Yang D.-L. Zhang G. Tang K. Li J. Yang L., 2016 Dicer-independent RNA-directed DNA methylation in Arabidopsis. Cell Res. 26: 66–82.

    Not in our library10.1038/cr.2015.145

25. 25.
    Primer-BLAST: A tool to design target-specific primers for polymerase chain reaction · 2012

    No notice
26. 26.

    Ye R. Chen Z. Lian B. Rowley M. J. Xia N., 2016 A dicer-independent route for biogenesis of siRNAs that direct DNA methylation in Arabidopsis. Mol. Cell 61: 222–235.

    Not in our library10.1016/j.molcel.2015.11.015

27. 27.

    Zhai J. Bischof S. Wang H. Feng S. Lee T.-F., 2015 A one precursor one siRNA model for Pol IV-dependent siRNA biogenesis. Cell 163: 445–455.

    Not in our library10.1016/j.cell.2015.09.032

Cited by papers we hold 4

• Heat‐responsive microRNAs and phased small interfering RNAs in reproductive development of flax · 2022
  Not checked
• Establishment of H3K9me3-dependent heterochromatin during embryogenesis in Drosophila miranda · 2021
  Not checked
• First Sprayable Double-Stranded RNA-Based Biopesticide Product Targets Proteasome Subunit Beta Type-5 in Colorado Potato Beetle (Leptinotarsa decemlineata) · 2021
  No notice
• Biogenesis of phased siRNAs on membrane-bound polysomes in Arabidopsis · 2016
  No notice

Reads

A person counts once a day per paper. We keep no cookie and no address: a reader is a hash that changes every day. Any month under five reads shows no number.