tsRNAsearch: a pipeline for the identification of tRNA and ncRNA fragments from small RNA-sequencing data

Motivation: tRNAs were originally considered uni-functional RNA molecules involved in the delivery of amino acids to growing peptide chains on the ribosome. More recently, the liberation of tRNA fragments from tRNAs via specific enzyme cleavage has been characterized. Detection of tRNA fragments in sequencing data is difficult due to tRNA sequence redundancy and the short length of both tRNAs and their fragments. Results: Here we introduce tsRNAsearch, a Nextflow pipeline for the identification of differentially abundant tRNA fragments and other non-coding RNAs from small RNA-sequencing data. tsRNAsearch is intended for use when comparing two groups of datasets, such as control and treatment groups. tsRNAsearch comparatively searches for tRNAs and ncRNAs with irregular read distribution profiles (a proxy for RNA cleavage) using a combined score made up of four novel methods and a differential expression analysis, and reports the top ranked results in simple PDF and TEXT files. In this study, we used publicly available small RNA-seq data to replicate the identification of tsRNAs from chronic hepatitis-infected liver tissue data. In addition, we applied tsRNAsearch to pancreatic ductal adenocarcinoma (PDAC) and matched healthy pancreatic tissue small RNA-sequencing data. Our results support the identification of miR135b from the original study as a potential biomarker of PDAC and identify other potentially stronger miRNA biomarkers of PDAC. Availability: https://github.com/GiantSpaceRobot/tsRNAsearch. Supplementary information: Supplementary data are available at Bioinformatics online.