Thomas Desvignes2, Karen EIlbeck3, Ioannis S. Vlachos4, Bastian Fromm5, Marc K. Halushka6, Michael Hackenberg7, Gianvito Urgese8, Elisa Ficarra8, Shruthi Bandyadka9, Jason Sydes2, Peter Batzel2, John H. Postlethwait2, Phillipe Loher10, Eric Londin10, Aristeidis G. Telonis10, Isidore Rigoutsos10, Lorena Pantano1
1 Harvard School TH Chan of Public Health, Boston, MA, USA. Email: firstname.lastname@example.org 2 Institute of Neuroscience, University of Oregon, Eugene, OR, USA. 3 University of Utah, Biomedical Informatics, UT, USA. 4 Brigham & Women’s Hospital, Broad Institute of MIT and Harvard, Harvard Medical School, Cambridge, MA, USA. 5 Department of Tumor Biology. Institute for Cancer Research. Norwegian Radium Hospital. Oslo University Hospital, Norwegian. 6 Department of Pathology, Johns Hopkins University School of Medicine, Baltimore, MD, USA. 7 Bioinformatics Group. University of Granada, Spain. 8 Politecnico di Torino, Italy. 9 Partners Personalized Medicine, Cambridge, MA, USA. 10 Computational Medicine Center, Thomas Jefferson University, Philadelphia, PA, USA
Project Website: http://mirtop.github.io Source Code: https://github.com/miRTop/mirtop License: MIT
MicroRNAs (miRNAs) are small RNA molecules (20-27 nt long) that are involved in eukaryotic gene regulation. They regulate targeted genes by RNA complementarity, generally between the miRNA seed region and the 3’ UTR of messenger RNAs. They have become relevant to the study of developmental stages and diseases. With the advance of sequencing technology, the cost of detecting the miRNA transcriptome has decreased substantially, and has led to the discovery of isomiRs (with slight sequence variations relative to the annotated miRNAs). As a consequence, the amount of miRNA data has increased exponentially in the last 5 years, together with the number of pipelines to analyze it. However, there is still a lack of consensus standards in storing and sharing such data, generating diverse output files; hindering comparison between tools, data sharing, and development of downstream analyses independent of the tools used for miRNA detection and quantification. Here, we present a community based, open source project to work toward the standardization of miRNA pipelines and encourage the development of downstream tools for visualization, differential expression, sample clustering, and model prediction analyses. This project is an international collaboration, with experts from different countries that have been developing miRNA pipelines and resources. We have described a standard format for the output of miRNA detection and quantification tools using small RNA-seq data. The format is based on the GFF3 standard in order to support reference coordinates and parent/child relationships between the features. The format contains information foreach sequence and its annotation to the miRNA precursor, the definition of reference or isomiR sequence, its quality, the isomiR type, and abundance in the data set. Moreover, we support a command line python tool to manage the miRNA GFF3 format (miRTop). Currently, miRTop can convert the output of commonly used small RNA-Seq pipelines, such as seqbuster, isomiR-SEA, sRNAbench, and Prost!, as well as BAM files to the miRNA GFF3 format. Importantly, the miRge pipeline has adapted the GFF3 format natively. miRTop can convert miRNA GFF3 files to a count matrix that can be easily imported to any downstream tool (i.e. for differential expression analysis). Furthermore, the tool can export the miRNA GFF file to the isomiRs package, improving the usability of this tool and rendering it independent of the upstream methods used for quantifying miRNA sequences. Together, we believe that the miRTop project will improve accessibility and uniformity of the results oft miRNA pipelines, the possibility to easily run benchmarking analyses, and promote the development of miRNA tools downstream of the detection and quantification steps.