Systematic overexpression of genes encoded by mycobacteriophage Waterfoul reveals novel inhibitors of mycobacterial growth

Mycobacteriophages: therapeutic approach for mycobacterial infections

Tuberculosis (TB) is a significant global health threat, and cases of infection with non-tuberculous mycobacteria (NTM) causing lung disease (NTM-LD) are rising. Bacteriophages and their gene products have garnered interest as potential therapeutic options for bacterial infections. Here, we have compiled information on bacteriophages and their products that can kill Mycobacterium tuberculosis or NTM. We summarize the mechanisms whereby viable phages can access macrophage-resident bacteria and not elicit immune responses, review methodologies of pharmaceutical product development containing mycobacteriophages and their gene products, mainly lysins, in the context of drug regulatory requirements and we discuss industrially relevant methods for producing pharmaceutical products comprising mycobacteriophages, emphasizing delivery of mycobacteriophages to the lungs. We conclude with an outline of some recent case studies on mycobacteriophage therapy.

A genome-wide cytotoxicity screen of cluster F1 mycobacteriophage Girr reveals novel inhibitors of Mycobacterium smegmatis growth

Over the past decade, thousands of bacteriophage genomes have been sequenced and annotated. A striking observation from this work is that known structural features and functions cannot be assigned for >65% of the encoded proteins. One approach to begin experimentally elucidating the function of these uncharacterized gene products is genome-wide screening to identify phage genes that confer phenotypes of interest like inhibition of host growth. This study describes the results of a screen evaluating the effects of overexpressing each gene encoded by the temperate Cluster F1 mycobacteriophage Girr on the growth of the host bacterium Mycobacterium smegmatis. Overexpression of 29 of the 102 Girr genes (~28% of the genome) resulted in mild to severe cytotoxicity. Of the 29 toxic genes described, 12 have no known function and are predominately small proteins of <125 amino acids. Overexpression of the majority of these 12 cytotoxic no known functions proteins resulted in moderate to severe growth reduction and represent novel antimicrobial products. The remaining 17 toxic genes have predicted functions, encoding products involved in phage structure, DNA replication/modification, DNA binding/gene regulation, or other enzymatic activity. Comparison of this dataset with prior genome-wide cytotoxicity screens of mycobacteriophages Waterfoul and Hammy reveals some common functional themes, though several of the predicted Girr functions associated with cytotoxicity in our report, including genes involved in lysogeny, have not been described previously. This study, completed as part of the HHMI-supported SEA-GENES project, highlights the power of parallel, genome-wide overexpression screens to identify novel interactions between phages and their hosts.

A genome-wide overexpression screen reveals Mycobacterium smegmatis growth inhibitors encoded by mycobacteriophage Hammy

During infection, bacteriophages produce diverse gene products to overcome bacterial antiphage defenses, to outcompete other phages, and to take over cellular processes. Even in the best-studied model phages, the roles of most phage-encoded gene products are unknown, and the phage population represents a largely untapped reservoir of novel gene functions. Considering the sheer size of this population, experimental screening methods are needed to sort through the enormous collection of available sequences and identify gene products that can modulate bacterial behavior for downstream functional characterization. Here, we describe the construction of a plasmid-based overexpression library of 94 genes encoded by Hammy, a Cluster K mycobacteriophage closely related to those infecting clinically important mycobacteria. The arrayed library was systematically screened in a plate-based cytotoxicity assay, identifying a diverse set of 24 gene products (representing ∼25% of the Hammy genome) capable of inhibiting growth of the host bacterium Mycobacterium smegmatis. Half of these are related to growth inhibitors previously identified in related phage Waterfoul, supporting their functional conservation; the other genes represent novel additions to the list of known antimycobacterial growth inhibitors. This work, conducted as part of the HHMI-supported Science Education Alliance Gene-function Exploration by a Network of Emerging Scientists (SEA-GENES) project, highlights the value of parallel, comprehensive overexpression screens in exploring genome-wide patterns of phage gene function and novel interactions between phages and their hosts.

The Cytotoxic Mycobacteriophage Protein Phaedrus gp82 Interacts with and Modulates the Activity of the Host ATPase, MoxR

Approximately 70% of bacteriophage-encoded proteins are of unknown function. Elucidating these protein functions represents opportunities to discover new phage-host interactions and mechanisms by which the phages modulate host activities. Here, we describe a pipeline for prioritizing phage-encoded proteins for structural analysis and characterize the gp82 protein encoded by mycobacteriophage Phaedrus. Structural and solution studies of gp82 show it is a trimeric protein containing two domains. Co-precipitation studies with the host Mycobacterium smegmatis identified the ATPase MoxR as an interacting partner protein. Phaedrus gp82-MoxR interaction requires the presence of a loop sequence within gp82 that is highly exposed and disordered in the crystallographic structure. We show that Phaedrus gp82 overexpression in M. smegmatis retards the growth of M. smegmatis on solid medium, resulting in a small colony phenotype. Overexpression of gp82 containing a mutant disordered loop or the overexpression of MoxR both rescue this phenotype. Lastly, we show that recombinant gp82 reduces levels of MoxR-mediated ATPase activity in vitro that is required for its chaperone function, and that the disordered loop plays an important role in this phenotype. We conclude that Phaedrus gp82 binds to and reduces mycobacterial MoxR activity, leading to reduced function of host proteins that require MoxR chaperone activity for their normal activity.

Identification of Toxic Proteins Encoded by Mycobacteriophage TM4 Using a Next-Generation Sequencing-Based Method

Mycobacteriophages are viruses that specifically infect mycobacteria and which, due to their diversity, represent a large gene pool. Characterization of the function of these genes should provide useful insights into host-phage interactions. Here, we describe a next-generation sequencing (NGS)-based, high-throughput screening approach for the identification of mycobacteriophage-encoded proteins that are toxic to mycobacteria. A plasmid-derived library representing the mycobacteriophage TM4 genome was constructed and transformed into Mycobacterium smegmatis. NGS and growth assays showed that the expression of TM4 gp43, gp77, -78, and -79, or gp85 was toxic to M. smegmatis. Although the genes associated with bacterial toxicity were expressed during phage infection, they were not required for lytic replication of mycobacteriophage TM4. In conclusion, we describe here an NGS-based approach which required significantly less time and resources than traditional methods and allowed the identification of novel mycobacteriophage gene products that are toxic to mycobacteria. IMPORTANCE The wide spread of drug-resistant Mycobacterium tuberculosis has brought an urgent need for new drug development. Mycobacteriophages are natural killers of M. tuberculosis, and their toxic gene products might provide potential anti-M. tuberculosis candidates. However, the enormous genetic diversity of mycobacteriophages poses challenges for the identification of these genes. Here, we used a simple and convenient screening method, based on next-generation sequencing, to identify mycobacteriophage genes encoding toxic products for mycobacteria. Using this approach, we screened and validated several toxic products encoded by mycobacteriophage TM4. In addition, we also found that the genes encoding these toxic products are nonessential for lytic replication of TM4. Our work describes a promising method for the identification of phage genes that encode proteins that are toxic to mycobacteria and which might facilitate the identification of novel antimicrobial molecules.

Genome Sequence of the Cluster CD Gordonia Phage Widow

Widow is a novel cluster CD bacteriophage isolated from a soil sample using the bacterial host Gordonia terrae. The Widow genome is 43,656 bp in length and encodes 64 protein-coding genes and no tRNAs. The genome shares 52 to 92% gene content with other cluster CD members.