PEPMatch: a tool to identify short peptide sequence matches in large sets of proteins.

BACKGROUND: Numerous tools exist for biological sequence comparisons and search. One case of particular interest for immunologists is finding matches for linear peptide T cell epitopes, typically between 8 and 15 residues in length, in a large set of protein sequences. Both to find exact matches or matches that account for residue substitutions. The utility of such tools is critical in applications ranging from identifying conservation across viral epitopes, identifying putative epitope targets for allergens, and finding matches for cancer-associated neoepitopes to examine the role of tolerance in tumor recognition. RESULTS: We defined a set of benchmarks that reflect the different practical applications of short peptide sequence matching. We evaluated a suite of existing methods for speed and recall and developed a new tool, PEPMatch. The tool uses a deterministic k-mer mapping algorithm that preprocesses proteomes before searching, achieving a 50-fold increase in speed over methods such as the Basic Local Alignment Search Tool (BLAST) without compromising recall. PEPMatch's code and benchmark datasets are publicly available. CONCLUSIONS: PEPMatch offers significant speed and recall advantages for peptide sequence matching. While it is of immediate utility for immunologists, the developed benchmarking framework also provides a standard against which future tools can be evaluated for improvements. The tool is available at https://nextgen-tools.iedb.org , and the source code can be found at https://github.com/IEDB/PEPMatch .

5' tRNA halves are highly expressed in the primate hippocampus and might sequence-specifically regulate gene expression.

Fragments of mature tRNAs have long been considered as mere degradation products without physiological function. However, recent reports show that tRNA-derived small RNAs (tsRNAs) play prominent roles in diverse cellular processes across a wide spectrum of species. Contrasting the situation in other small RNA pathways the mechanisms behind these effects appear more diverse, more complex, and are generally less well understood. In addition, surprisingly little is known about the expression profiles of tsRNAs across different tissues and species. Here, we provide an initial overview of tsRNA expression in different species and tissues, revealing very high levels of 5' tRNA halves (5' tRHs) particularly in the primate hippocampus. We further modulated the regulation capacity of selected 5' tRHs in human cells by transfecting synthetic tsRNA mimics ("overexpression") or antisense-RNAs ("inhibition") and identified differentially expressed transcripts based on RNA-seq. We then used a novel k-mer mapping approach to dissect the underlying targeting rules, suggesting that 5' tRHs silence genes in a sequence-specific manner, while the most efficient target sites align to the mid-region of the 5' tRH and are located within the CDS or 3' UTR of the target. This amends previous observations that tsRNAs guide Argonaute proteins to silence their targets via a miRNA-like 5' seed match and suggests a yet unknown mechanism of regulation. Finally, our data suggest that some 5' tRHs that are also able to sequence-specifically stabilize mRNAs as up-regulated mRNAs are also significantly enriched for 5' tRH target sites.