Thousands of small, novel genes predicted in global phage genomes

A catalog of small proteins from the global microbiome

Small open reading frames (smORFs) shorter than 100 codons are widespread and perform essential roles in microorganisms, where they encode proteins active in several cell functions, including signal pathways, stress response, and antibacterial activities. However, the ecology, distribution and role of small proteins in the global microbiome remain unknown. Here, we construct a global microbial smORFs catalog (GMSC) derived from 63,410 publicly available metagenomes across 75 distinct habitats and 87,920 high-quality isolate genomes. GMSC contains 965 million non-redundant smORFs with comprehensive annotations. We find that archaea harbor more smORFs proportionally than bacteria. We moreover provide a tool called GMSC-mapper to identify and annotate small proteins from microbial (meta)genomes. Overall, this publicly-available resource demonstrates the immense and underexplored diversity of small proteins.

Phage small proteins play large roles in phage-bacterial interactions

Phages have wide influence on bacterial physiology, and likewise, bacterial processes impinge on phage biology. Key to these interactions are phage small proteins (<100 aa). Long underappreciated, recent work has revealed millions of phage small proteins, and increasingly, mechanisms by which they function to dictate phage and/or bacterial behavior and evolution. Here, we describe select phage small proteins that mediate phage-bacterial interactions by modulating phage lifestyle decision-making components or by altering host gene expression.

Large-scale investigation of species-specific orphan genes in the human gut microbiome elucidates their evolutionary origins

Species-specific genes, also known as orphans, are ubiquitous across life's domains. In prokaryotes, species-specific orphan genes (SSOGs) are mostly thought to originate in external elements such as viruses followed by horizontal gene transfer, whereas the scenario of native origination, through rapid divergence or de novo, is mostly dismissed. However, quantitative evidence supporting either scenario is lacking. Here, we systematically analyzed genomes from 4644 human gut microbiome species and identified more than 600,000 unique SSOGs, representing an average of 2.6% of a given species' pangenome. These sequences are mostly rare within each species yet show signs of purifying selection. Overall, SSOGs use optimal codons less frequently, and their proteins are more disordered than those of conserved genes (i.e., non-SSOGs). Importantly, across species, the GC content of SSOGs closely matches that of conserved ones. In contrast, the ∼5% of SSOGs that share similarity to known viral sequences have distinct characteristics, including lower GC content. Thus, SSOGs with similarity to viruses differ from the remaining SSOGs, contrasting an external origination scenario for most of them. By examining the orthologous genomic region in closely related species, we show that a small subset of SSOGs likely evolved natively de novo and find that these genes also differ in their properties from the remaining SSOGs. Our results challenge the notion that external elements are the dominant source of prokaryotic genetic novelty and will enable future studies into the biological role and relevance of species-specific genes in the human gut.

An enterococcal phage protein broadly inhibits type IV restriction enzymes involved in antiphage defense

The prevalence of multidrug resistant (MDR) bacterial infections continues to rise as the development of antibiotics needed to combat these infections remains stagnant. MDR enterococci are a major contributor to this crisis. A potential therapeutic approach for combating MDR enterococci is bacteriophage (phage) therapy, which uses lytic viruses to infect and kill pathogenic bacteria. While phages that lyse some strains of MDR enterococci have been identified, other strains display high levels of resistance and the mechanisms underlying this resistance are poorly defined. Here, we use a CRISPR interference (CRISPRi) screen to identify a genetic locus found on a mobilizable plasmid from Enterococcus faecalis involved in phage resistance. This locus encodes a putative serine recombinase followed by a Type IV restriction enzyme (TIV-RE) that we show restricts the replication of phage phi47 in E. faecalis. We further find that phi47 evolves to overcome restriction by acquiring a missense mutation in a TIV-RE inhibitor protein. We show that this inhibitor, termed type IV restriction inhibiting factor A (tifA), binds and inactivates diverse TIV-REs. Overall, our findings advance our understanding of phage defense in drug-resistant E. faecalis and provide mechanistic insight into how phages evolve to overcome antiphage defense systems.

Bacteriophage-Based Bioanalysis

Bacteriophages, which are viral predators of bacteria, have evolved to efficiently recognize, bind, infect, and lyse their host, resulting in the release of tens to hundreds of propagated viruses. These abilities have attracted biosensor developers who have developed new methods to detect bacteria. Recently, several comprehensive reviews have covered many of the advances made regarding the performance of phage-based biosensors. Therefore, in this review, we first describe the landscape of phage-based biosensors and then cover advances in other aspects of phage biology and engineering that can be used to make high-impact contributions to biosensor development. Many of these advances are in fields adjacent to analytical chemistry such as synthetic biology, machine learning, and genetic engineering and will allow those looking to develop phage-based biosensors to start taking alternative approaches, such as a bottom-up design and synthesis of custom phages with the singular task of detecting their host.

A pangenome analysis of ESKAPE bacteriophages: the underrepresentation may impact machine learning models

Bacteriophages are the most prevalent biological entities in the biosphere. However, limitations in both medical relevance and sequencing technologies have led to a systematic underestimation of the genetic diversity within phages. This underrepresentation not only creates a significant gap in our understanding of phage roles across diverse biosystems but also introduces biases in computational models reliant on these data for training and testing. In this study, we focused on publicly available genomes of bacteriophages infecting high-priority ESKAPE pathogens to show the extent and impact of this underrepresentation. First, we demonstrate a stark underrepresentation of ESKAPE phage genomes within the public genome and protein databases. Next, a pangenome analysis of these ESKAPE phages reveals extensive sharing of core genes among phages infecting the same host. Furthermore, genome analyses and clustering highlight close nucleotide-level relationships among the ESKAPE phages, raising concerns about the limited diversity within current public databases. Lastly, we uncover a scarcity of unique lytic phages and phage proteins with antimicrobial activities against ESKAPE pathogens. This comprehensive analysis of the ESKAPE phages underscores the severity of underrepresentation and its potential implications. This lack of diversity in phage genomes may restrict the resurgence of phage therapy and cause biased outcomes in data-driven computational models due to incomplete and unbalanced biological datasets.

Phage transmission strategies: are phages farming their host?

Extensive coevolution has led to utterly intricate interactions between phages and their bacterial hosts. While both the (short-term) intracellular molecular host-subversion mechanisms during a phage infection cycle and the (long-term) mutational arms race between phages and host cells have traditionally received a lot of attention, there has been an underestimating neglect of (mid-term) transmission strategies by which phages manage to cautiously spread throughout their host population. However, recent findings underscore that phages encode mechanisms to avoid host cell scarcity and promote coexistence with the host, giving the impression that some phages manage to 'farm' their host population to ensure access to host cells for lytic consumption. Given the tremendous impact of phages on bacterial ecology, charting and understanding the complexity of such transmission strategies is of key importance.

Discovery and description of novel phage genomes from urban microbiomes sampled by the MetaSUB consortium

Bacteriophages are recognized as the most abundant members of microbiomes and have therefore a profound impact on microbial communities through the interactions with their bacterial hosts. The International Metagenomics and Metadesign of Subways and Urban Biomes Consortium (MetaSUB) has sampled mass-transit systems in 60 cities over 3 years using metagenomics, throwing light into these hitherto largely unexplored urban environments. MetaSUB focused primarily on the bacterial community. In this work, we explored MetaSUB metagenomic data in order to recover and analyze bacteriophage genomes. We recovered and analyzed 1714 phage genomes with size at least 40 kbp, from the class Caudoviricetes, the vast majority of which (80%) are novel. The recovered genomes were predicted to belong to temperate (69%) and lytic (31%) phages. Thirty-three of these genomes have more than 200 kbp, and one of them reaches 572 kbp, placing it among the largest phage genomes ever found. In general, the phages tended to be site-specific or nearly so, but 194 genomes could be identified in every city from which phage genomes were retrieved. We predicted hosts for 48% of the phages and observed general agreement between phage abundance and the respective bacterial host abundance, which include the most common nosocomial multidrug-resistant pathogens. A small fraction of the phage genomes are carriers of antibiotic resistance genes, and such genomes tended to be particularly abundant in the sites where they were found. We also detected CRISPR-Cas systems in five phage genomes. This study expands the previously reported MetaSUB results and is a contribution to the knowledge about phage diversity, global distribution, and phage genome content.

Preclinical characterization and in silico safety assessment of three virulent bacteriophages targeting carbapenem-resistant uropathogenic Escherichia coli

Phage therapy has recently been revitalized in the West with many successful applications against multi-drug-resistant bacterial infections. However, the lack of geographically diverse bacteriophage (phage) genomes has constrained our understanding of phage diversity and its genetics underpinning host specificity, lytic capability, and phage-bacteria co-evolution. This study aims to locally isolate virulent phages against uropathogenic Escherichia coli (E. coli) and study its phenotypic and genomic features. Three obligately virulent Escherichia phages (øEc_Makalu_001, øEc_Makalu_002, and øEc_Makalu_003) that could infect uropathogenic E. coli were isolated and characterized. All three phages belonged to Krischvirus genus. One-step growth curve showed that the latent period of the phages ranged from 15 to 20 min, the outbreak period ~ 50 min, and the burst size ranged between 74 and 127 PFU/bacterium. Moreover, the phages could tolerate a pH range of 6 to 9 and a temperature range of 25-37 °C for up to 180 min without significant loss of phage viability. All phages showed a broad host spectrum and could lyse up to 30% of the 35 tested E. coli isolates. Genomes of all phages were approximately ~ 163 kb with a gene density of 1.73 gene/kbp and an average gene length of ~ 951 bp. The coding density in all phages was approximately 95%. Putative lysin, holin, endolysin, and spanin genes were found in the genomes of all three phages. All phages were strictly virulent with functional lysis modules and lacked any known virulence or toxin genes and antimicrobial resistance genes. Pre-clinical experimental and genomic analysis suggest these phages may be suitable candidates for therapeutic applications.

MetaCerberus: distributed highly parallelized HMM-based processing for robust functional annotation across the tree of life

MOTIVATION

MetaCerberus is a massively parallel, fast, low memory, scalable annotation tool for inference gene function across genomes to metacommunities. MetaCerberus provides an elusive HMM/HMMER-based tool at a rapid scale with low memory. It offers scalable gene elucidation to major public databases, including KEGG (KO), COGs, CAZy, FOAM, and specific databases for viruses, including VOGs and PHROGs, from single genomes to metacommunities.

RESULTS

MetaCerberus is 1.3× as fast on a single node than eggNOG-mapper v2 on 5× less memory using an exclusively HMM/HMMER mode. In a direct comparison, MetaCerberus provides better annotation of viruses, phages, and archaeal viruses than DRAM, Prokka, or InterProScan. MetaCerberus annotates more KOs across domains when compared to DRAM, with a 186× smaller database, and with 63× less memory. MetaCerberus is fully integrated for automatic analysis of statistics and pathways using differential statistic tools (i.e. DESeq2 and edgeR), pathway enrichment (GAGE R), and pathview R. MetaCerberus provides a novel tool for unlocking the biosphere across the tree of life at scale.

AVAILABILITY AND IMPLEMENTATION

MetaCerberus is written in Python and distributed under a BSD-3 license. The source code of MetaCerberus is freely available at https://github.com/raw-lab/metacerberus compatible with Python 3 and works on both Mac OS X and Linux. MetaCerberus can also be easily installed using bioconda: mamba create -n metacerberus -c bioconda -c conda-forge metacerberus.

The Cryptic Bacterial Microproteome

Microproteins encoded by small open reading frames (smORFs) comprise the "dark matter" of proteomes. Although functional microproteins were identified in diverse organisms from all three domains of life, bacterial smORFs remain poorly characterized. In this comprehensive study of intergenic smORFs (ismORFs, 15-70 codons) in 5,668 bacterial genomes of the family Enterobacteriaceae, we identified 67,297 clusters of ismORFs subject to purifying selection. The ismORFs mainly code for hydrophobic, potentially transmembrane, unstructured, or minimally structured microproteins. Using AlphaFold Multimer, we predicted interactions of some of the predicted microproteins encoded by transcribed ismORFs with proteins encoded by neighboring genes, revealing the potential of microproteins to regulate the activity of various proteins, particularly, under stress. We compiled a catalog of predicted microprotein families with different levels of evidence from synteny analysis, structure prediction, and transcription and translation data. This study offers a resource for investigation of biological functions of microproteins.

Obtaining Detailed Phage Transcriptomes Using ONT-Cappable-Seq

Detailed transcription maps of bacteriophages are not usually explored, limiting our understanding of molecular phage biology and restricting their exploitation and engineering. The ONT-cappable-seq method described here brings phage transcriptomics to the accessible nanopore sequencing platform and provides an affordable and more detailed overview of transcriptional features compared to traditional RNA-seq experiments. With ONT-cappable-seq, primary transcripts are specifically capped, enriched, and prepared for long-read sequencing on the nanopore sequencing platform. This enables end-to-end sequencing of unprocessed transcripts covering both phage and host genome, thus providing insight on their operons. The subsequent analysis pipeline makes it possible to rapidly identify the most important transcriptional features such as transcription start and stop sites. The obtained data can thus provide a comprehensive overview of the transcription by your phage of interest.

Small proteins in bacteria - Big challenges in prediction and identification

Proteins with up to 100 amino acids have been largely overlooked due to the challenges associated with predicting and identifying them using traditional methods. Recent advances in bioinformatics and machine learning, DNA sequencing, RNA and Ribo-seq technologies, and mass spectrometry (MS) have greatly facilitated the detection and characterisation of these elusive proteins in recent years. This has revealed their crucial role in various cellular processes including regulation, signalling and transport, as toxins and as folding helpers for protein complexes. Consequently, the systematic identification and characterisation of these proteins in bacteria have emerged as a prominent field of interest within the microbial research community. This review provides an overview of different strategies for predicting and identifying these proteins on a large scale, leveraging the power of these advanced technologies. Furthermore, the review offers insights into the future developments that may be expected in this field.

Systematic and scalable genome-wide essentiality mapping to identify nonessential genes in phages

Phages are one of the key ecological drivers of microbial community dynamics, function, and evolution. Despite their importance in bacterial ecology and evolutionary processes, phage genes are poorly characterized, hampering their usage in a variety of biotechnological applications. Methods to characterize such genes, even those critical to the phage life cycle, are labor intensive and are generally phage specific. Here, we develop a systematic gene essentiality mapping method scalable to new phage-host combinations that facilitate the identification of nonessential genes. As a proof of concept, we use an arrayed genome-wide CRISPR interference (CRISPRi) assay to map gene essentiality landscape in the canonical coliphages λ and P1. Results from a single panel of CRISPRi probes largely recapitulate the essential gene roster determined from decades of genetic analysis for lambda and provide new insights into essential and nonessential loci in P1. We present evidence of how CRISPRi polarity can lead to false positive gene essentiality assignments and recommend caution towards interpreting CRISPRi data on gene essentiality when applied to less studied phages. Finally, we show that we can engineer phages by inserting DNA barcodes into newly identified inessential regions, which will empower processes of identification, quantification, and tracking of phages in diverse applications.

Genome sequence of bacteriophage Aoka, a cluster FO phage isolated using Arthrobacter globiformis

We report the discovery and genome sequence of bacteriophage Aoka, an actinobacteriophage isolated from a soil sample in Pueblo, Colorado using Arthrobacter globiformis, B2880-SEA. Its genome length is 36,744 base pairs with 54 protein-coding genes. Based on gene content similarity to other actinobacteriophages, Aoka is assigned to cluster FO.

Microviruses: A World Beyond phiX174

Two decades of metagenomic analyses have revealed that in many environments, small (∼5 kb), single-stranded DNA phages of the family Microviridae dominate the virome. Although the emblematic microvirus phiX174 is ubiquitous in the laboratory, most other microviruses, particularly those of the gokushovirus and amoyvirus lineages, have proven to be much more elusive. This puzzling lack of representative isolates has hindered insights into microviral biology. Furthermore, the idiosyncratic size and nature of their genomes have resulted in considerable misjudgments of their actual abundance in nature. Fortunately, recent successes in microvirus isolation and improved metagenomic methodologies can now provide us with more accurate appraisals of their abundance, their hosts, and their interactions. The emerging picture is that phiX174 and its relatives are rather rare and atypical microviruses, and that a tremendous diversity of other microviruses is ready for exploration.

An extended catalog of integrated prophages in the infant and adult fecal microbiome shows high prevalence of lysogeny

Background and aims

The acquisition and gradual maturation of gut microbial communities during early childhood is central to an individual's healthy development. Bacteriophages have the potential to shape the gut bacterial communities. However, the complex ecological interactions between phages and their bacterial host are still poorly characterized. In this study, we investigated the abundance and diversity of integrated prophages in infant and adult gut bacteria by detecting integrated prophages in metagenome assembled genomes (MAGs) of commensal bacteria.

Methods

Our study included 88 infants sampled at 3 weeks, 3 months, 6 months, and 12 months (n = 323 total samples), and their parents around delivery time (n = 138 total samples). Fecal DNA was extracted and characterized by using shotgun metagenomic sequencing, and a collection of prokaryotic MAGs was generated. The MAG collection was screened for the presence of integrated bacteriophage sequences, allowing their taxonomic and functional characterization.

Results

A large collection of 6,186 MAGs from infant and adult gut microbiota was obtained and screened for integrated prophages, allowing the identification of 7,165 prophage sequences longer than 10 kb. Strikingly, more than 70% of the near-complete MAGs were identified as lysogens. The prevalence of prophages in MAGs varied across bacterial families, with a lower prevalence observed among Coriobacteriaceae, Eggerthellaceae, Veillonellaceae and Burkholderiaceae, while a very high prevalence of lysogen MAGs were observed in Oscillospiraceae, Enterococcaceae, and Enterobacteriaceae. Interestingly for several bacterial families such as Bifidobacteriaceae and Bacteroidaceae, the prevalence of prophages in MAGs was higher in early infant time point (3 weeks and 3 months) than in later sampling points (6 and 12 months) and in adults. The prophage sequences were clustered into 5,616 species-like vOTUs, 77% of which were novel. Finally, we explored the functional repertoire of the potential auxiliary metabolic genes carried by these prophages, encoding functions involved in carbohydrate metabolism and degradation, amino acid metabolism and carbon metabolism.

Conclusion

Our study provides an enhanced understanding of the diversity and prevalence of lysogens in infant and adult gut microbiota and suggests a complex interplay between prophages and their bacterial hosts.

Prophage identification and molecular analysis in the genomes of Pseudomonas aeruginosa strains isolated from critical care patients

Prophages are bacteriophages integrated into the bacterial host's chromosome. This research aims to analyze and characterize the existing prophages within a collection of 53 Pseudomonas aeruginosa strains from intensive care units (ICUs) in Portugal and Spain. A total of 113 prophages were localized in the collection, with 18 of them being present in more than one strain simultaneously. After annotation, five of them were discarded as incomplete, and the 13 remaining prophages were characterized. Of 13, 10 belonged to the siphovirus tail morphology group, 2 to the podovirus tail morphology group, and 1 to the myovirus tail morphology group. All prophages had a length ranging from 20,199 to 63,401 bp and a GC% between 56.2% and 63.6%. The number of open reading frames (ORFs) oscillated between 32 and 88, and in 3/13 prophages, more than 50% of the ORFs had an unknown function. With our findings, we show that prophages are present in the majority of the P. aeruginosa strains isolated from Portuguese and Spanish critically ill patients, many of them found in more than one circulating strain at the same time and following a similar clonal distribution pattern. Although a great sum of ORFs had an unknown function, number of proteins in relation to viral defense (anti-CRISPR proteins, toxin/antitoxin modules, proteins against restriction-modification systems) as well as to prophage interference into their host's quorum sensing system and regulatory cascades were found. This supports the idea that prophages have an influence in bacterial pathogenesis and anti-phage defense. IMPORTANCE Despite being known for decades, prophages remain understudied when compared to the lytic phages employed in phage therapy. This research aims to shed some light into the nature, composition, and role of prophages found within a set of circulating strains of Pseudomas aeruginosa, with special attention to high-risk clones. Given the fact that prophages can effectively influence bacterial pathogenesis, prophage basic research constitutes a topic of growing interest. Furthermore, the abundance of viral defense and regulatory proteins within prophage genomes detected in this study evidences the importance of characterizing the most frequent prophages in circulating clinical strains and in high-risk clones if phage therapy is to be used.

Functional biology and biotechnology of thermophilic viruses

Viruses have developed sophisticated biochemical and genetic mechanisms to manipulate and exploit their hosts. Enzymes derived from viruses have been essential research tools since the first days of molecular biology. However, most viral enzymes that have been commercialized are derived from a small number of cultivated viruses, which is remarkable considering the extraordinary diversity and abundance of viruses revealed by metagenomic analysis. Given the explosion of new enzymatic reagents derived from thermophilic prokaryotes over the past 40 years, those obtained from thermophilic viruses should be equally potent tools. This review discusses the still-limited state of the art regarding the functional biology and biotechnology of thermophilic viruses with a focus on DNA polymerases, ligases, endolysins, and coat proteins. Functional analysis of DNA polymerases and primase-polymerases from phages infecting Thermus, Aquificaceae, and Nitratiruptor has revealed new clades of enzymes with strong proofreading and reverse transcriptase capabilities. Thermophilic RNA ligase 1 homologs have been characterized from Rhodothermus and Thermus phages, with both commercialized for circularization of single-stranded templates. Endolysins from phages infecting Thermus, Meiothermus, and Geobacillus have shown high stability and unusually broad lytic activity against Gram-negative and Gram-positive bacteria, making them targets for commercialization as antimicrobials. Coat proteins from thermophilic viruses infecting Sulfolobales and Thermus strains have been characterized, with diverse potential applications as molecular shuttles. To gauge the scale of untapped resources for these proteins, we also document over 20,000 genes encoded by uncultivated viral genomes from high-temperature environments that encode DNA polymerase, ligase, endolysin, or coat protein domains.

Siphophage 0105phi7-2 of Bacillus thuringiensis: Novel Propagation, DNA, and Genome-Implied Assembly

Diversity of phage propagation, physical properties, and assembly promotes the use of phages in ecological studies and biomedicine. However, observed phage diversity is incomplete. Bacillus thuringiensis siphophage, 0105phi-7-2, first described here, significantly expands known phage diversity, as seen via in-plaque propagation, electron microscopy, whole genome sequencing/annotation, protein mass spectrometry, and native gel electrophoresis (AGE). Average plaque diameter vs. plaque-supporting agarose gel concentration plots reveal unusually steep conversion to large plaques as agarose concentration decreases below 0.2%. These large plaques sometimes have small satellites and are made larger by orthovanadate, an ATPase inhibitor. Phage head-host-cell binding is observed by electron microscopy. We hypothesize that this binding causes plaque size-increase via biofilm evolved, ATP stimulated ride-hitching on motile host cells by temporarily inactive phages. Phage 0105phi7-2 does not propagate in liquid culture. Genomic sequencing/annotation reveals history as temperate phage and distant similarity, in a virion-assembly gene cluster, to prototypical siphophage SPP1 of Bacillus subtilis. Phage 0105phi7-2 is distinct in (1) absence of head-assembly scaffolding via either separate protein or classically sized, head protein-embedded peptide, (2) producing partially condensed, head-expelled DNA, and (3) having a surface relatively poor in AGE-detected net negative charges, which is possibly correlated with observed low murine blood persistence.

Forest microbiome and global change

Forests influence climate and mitigate global change through the storage of carbon in soils. In turn, these complex ecosystems face important challenges, including increases in carbon dioxide, warming, drought and fire, pest outbreaks and nitrogen deposition. The response of forests to these changes is largely mediated by microorganisms, especially fungi and bacteria. The effects of global change differ among boreal, temperate and tropical forests. The future of forests depends mostly on the performance and balance of fungal symbiotic guilds, saprotrophic fungi and bacteria, and fungal plant pathogens. Drought severely weakens forest resilience, as it triggers adverse processes such as pathogen outbreaks and fires that impact the microbial and forest performance for carbon storage and nutrient turnover. Nitrogen deposition also substantially affects forest microbial processes, with a pronounced effect in the temperate zone. Considering plant-microorganism interactions would help predict the future of forests and identify management strategies to increase ecosystem stability and alleviate climate change effects. In this Review, we describe the impact of global change on the forest ecosystem and its microbiome across different climatic zones. We propose potential approaches to control the adverse effects of global change on forest stability, and present future research directions to understand the changes ahead.

Pharokka: a fast scalable bacteriophage annotation tool

SUMMARY

In recent years, there has been an increasing interest in bacteriophages, which has led to growing numbers of bacteriophage genomic sequences becoming available. Consequently, there is a need for a rapid and consistent genomic annotation tool dedicated for bacteriophages. Existing tools either are not designed specifically for bacteriophages or are web- and email-based and require significant manual curation, which makes their integration into bioinformatic pipelines challenging. Pharokka was created to provide a tool that annotates bacteriophage genomes easily, rapidly and consistently with standards compliant outputs. Moreover, Pharokka requires only two lines of code to install and use and takes under 5 min to run for an average 50-kb bacteriophage genome.

AVAILABILITY AND IMPLEMENTATION

Pharokka is implemented in Python and is available as a bioconda package using 'conda install -c bioconda pharokka'. The source code is available on GitHub (https://github.com/gbouras13/pharokka). Pharokka has been tested on Linux-64 and MacOSX machines and on Windows using a Linux Virtual Machine.