Abstract:Decades of overconsumption of antimicrobials in the treatment and prevention of bacterial infections have resulted in the increasing emergence of drug-resistant bacteria, which poses a significant challenge to public health, driving the urgent need to find alternatives to conventional antibiotics. Bacteriophages are viruses infecting specific bacterial hosts, often destroying the infected bacterial hosts. Phages attach to and enter their potential hosts using their tail proteins, with the composition of the tail determining the range of potentially infected bacteria. To aid the exploitation of bacteriophages for therapeutic purposes, we developed the PhageTailFinder algorithm to predict tail-related proteins and identify the putative tail module in previously uncharacterized phages. The PhageTailFinder relies on a two-state hidden Markov model (HMM) to predict the probability of a given protein being tail-related. The process takes into account the natural modularity of phage tail-related proteins, rather than simply considering amino acid properties or secondary structures for each protein in isolation. The PhageTailFinder exhibited robust predictive power for phage tail proteins in novel phages due to this sequence-independent operation. The performance of the prediction model was evaluated in 13 extensively studied phages and a sample of 992 complete phages from the NCBI database. The algorithm achieved a high true-positive prediction rate (>80%) in over half (571) of the studied phages, and the ROC value was 0.877 using general models and 0.968 using corresponding morphologic models. It is notable that the median ROC value of 992 complete phages is more than 0.75 even for novel phages, indicating the high accuracy and specificity of the PhageTailFinder. When applied to a dataset containing 189,680 viral genomes derived from 11,810 bulk metagenomic human stool samples, the ROC value was 0.895. In addition, tail protein clusters could be identified for further studies by density-based spatial clustering of applications with the noise algorithm (DBSCAN). The developed PhageTailFinder tool can be accessed either as a web server (http://www.microbiome-bigdata.com/PHISDetector/index/tools/PhageTailFinder) or as a stand-alone program on a standard desktop computer (https://github.com/HIT-ImmunologyLab/PhageTailFinder).

SeqTailor: a user-friendly webserver for the extraction of DNA or protein sequences from next-generation sequencing data

Tailor: a computational framework for detecting non-templated tailing of small silencing RNAs

BioSeq-Analysis2.0: an Updated Platform for Analyzing DNA, RNA and Protein Sequences at Sequence Level and Residue Level Based on Machine Learning Approaches.

SeqCAT: Sequence Conversion and Analysis Toolbox

CoV-Seq, a New Tool for SARS-CoV-2 Genome Analysis and Visualization: Development and Usability Study

TaSer (TabAnno and SeqMiner): a toolset for annotating and querying next-generation sequence data

An Evolutionary Statistics Toolkit for Simplified Sequence Analysis on Web with Client-Side Processing

SaAlign: Multiple DNA/RNA Sequence Alignment and Phylogenetic Tree Construction Tool for Ultra-Large Datasets and Ultra-Long Sequences Based on Suffix Array

VCAT: an integrated variant function annotation tools

PolyAtailor: measuring poly(A) tail length from short-read and long-read sequencing data

BioSeq-Analysis: a platform for DNA, RNA and protein sequence analysis based on machine learning approaches

PhageTailFinder: A Tool for Phage Tail Module Detection and Annotation

DNAscan: a fast, computationally and memory efficient bioinformatics pipeline for the analysis of DNA next-generation-sequencing data

PROVEAN web server: a tool to predict the functional effect of amino acid substitutions and indels

Construction and Application of the Web Service "stdtransdb" for Database of Standard Transcript Sequences of Human and Model Species

OneStopRNAseq: A Web Application for Comprehensive and Efficient Analyses of RNA-Seq Data

Reseqtools: an Integrated Toolkit for Large-Scale Next-Generation Sequencing Based Resequencing Analysis

Seqminer2: an efficient tool to query and retrieve genotypes for statistical genetics analyses from biobank scale sequence dataset

SWAV: a web-based visualization browser for sliding window analysis

iSeq: An integrated tool to fetch public sequencing data

WebTraceMiner: a Web Service for Processing and Mining EST Sequence Trace Files.