Mycobacteriophages: genes and genomes

PRFect: a tool to predict programmed ribosomal frameshifts in prokaryotic and viral genomes

BACKGROUND

One of the stranger phenomena that can occur during gene translation is where, as a ribosome reads along the mRNA, various cellular and molecular properties contribute to stalling the ribosome on a slippery sequence and shifting the ribosome into one of the other two alternate reading frames. The alternate frame has different codons, so different amino acids are added to the peptide chain. More importantly, the original stop codon is no longer in-frame, so the ribosome can bypass the stop codon and continue to translate the codons past it. This produces a longer version of the protein, a fusion of the original in-frame amino acids, followed by all the alternate frame amino acids. There is currently no automated software to predict the occurrence of these programmed ribosomal frameshifts (PRF), and they are currently only identified by manual curation.

RESULTS

Here we present PRFect, an innovative machine-learning method for the detection and prediction of PRFs in coding genes of various types. PRFect combines advanced machine learning techniques with the integration of multiple complex cellular properties, such as secondary structure, codon usage, ribosomal binding site interference, direction, and slippery site motif. Calculating and incorporating these diverse properties posed significant challenges, but through extensive research and development, we have achieved a user-friendly approach. The code for PRFect is freely available, open-source, and can be easily installed via a single command in the terminal. Our comprehensive evaluations on diverse organisms, including bacteria, archaea, and phages, demonstrate PRFect's strong performance, achieving high sensitivity, specificity, and an accuracy exceeding 90%. The code for PRFect is freely available and installs with a single terminal command.

CONCLUSION

PRFect represents a significant advancement in the field of PRF detection and prediction, offering a powerful tool for researchers and scientists to unravel the intricacies of programmed ribosomal frameshifting in coding genes.

Diversity and conservation of the genome architecture of phages infecting the Alphaproteobacteria

IMPORTANCE

This study reports the results of the largest analysis of genome sequences from phages that infect the Alphaproteobacteria class of bacterial hosts. We analyzed over 100 whole genome sequences of phages to construct dotplots, categorize them into genetically distinct clusters, generate a bootstrapped phylogenetic tree, compute protein orthologs, and predict packaging strategies. We determined that the phage sequences primarily cluster by the bacterial host family, phage morphotype, and genome size. We expect that the findings reported in this seminal study will facilitate future analyses that will improve our knowledge of the phages that infect these hosts.

The heterogenous and diverse population of prophages in Mycobacterium genomes

IMPORTANCE

Mycobacterium species include several human pathogens and mycobacteriophages show potential for therapeutic use to control Mycobacterium infections. However, phage infection profiles vary greatly among Mycobacterium abscessus clinical isolates and phage therapies must be personalized for individual patients. Mycobacterium phage susceptibility is likely determined primarily by accessory parts of bacterial genomes, and we have identified the prophage and phage-related genomic regions across sequenced Mycobacterium strains. The prophages are numerous and diverse, especially in M. abscessus genomes, and provide a potentially rich reservoir of new viruses that can be propagated lytically and used to expand the repertoire of therapeutically useful phages.

Mycobacteriophages in diagnosis and alternative treatment of mycobacterial infections

Antimicrobial resistance is an increasing threat to human populations. The emergence of multidrug-resistant "superbugs" in mycobacterial infections has further complicated the processes of curing patients, thereby resulting in high morbidity and mortality. Early diagnosis and alternative treatment are important for improving the success and cure rates associated with mycobacterial infections and the use of mycobacteriophages is a potentially good option. Since each bacteriophage has its own host range, mycobacteriophages have the capacity to detect specific mycobacterial isolates. The bacteriolysis properties of mycobacteriophages make them more attractive when it comes to treating infectious diseases. In fact, they have been clinically applied in Eastern Europe for several decades. Therefore, mycobacteriophages can also treat mycobacteria infections. This review explores the potential clinical applications of mycobacteriophages, including phage-based diagnosis and phage therapy in mycobacterial infections. Furthermore, this review summarizes the current difficulties in phage therapy, providing insights into new treatment strategies against drug-resistant mycobacteria.

Phenotypic Characterization and Comparative Genomic Analyses of Mycobacteriophage WIVsmall as A New Member Assigned to F1 Subcluster

In this study, we conducted the morphological observation, biological and genomic characterization, evolutionary analysis, comparative genomics description, and proteome identification of a recently isolated mycobacteriophage, WIVsmall. Morphologically, WIVsmall is classified as a member of the Siphoviridae family, characterized by a flexible tail, measuring approximately 212 nm in length. The double-stranded phage genome DNA of WIVsmall spans 53,359 base pairs, and exhibits a G + C content of 61.01%. The genome of WIVsmall comprises 103 protein-coding genes, while no tRNA genes were detected. The genome annotation unveiled the presence of functional gene clusters responsible for mycobacteriophage assembly and maturation, replication, cell lysis, and functional protein synthesis. Based on the analysis of the phylogenetic tree, the genome of WIVsmall was classified as belonging to subgroup F1. A comparative genomics analysis indicated that the WIVsmall genome exhibited the highest similarity to the phage SG4, with a percentage of 64%. The single-step growth curve analysis of WIVsmall revealed a latent period of 120 min, and an outbreak period of 200 min.

PRFect: A tool to predict programmed ribosomal frameshifts in prokaryotic and viral genomes

Background

One of the stranger phenomena that can occur during gene translation is where, as a ribosome reads along the mRNA, various cellular and molecular properties contribute to stalling the ribosome on a slippery sequence, shifting the ribosome into one of the other two alternate reading frames. The alternate frame has different codons, so different amino acids are added to the peptide chain, but more importantly, the original stop codon is no longer in-frame, so the ribosome can bypass the stop codon and continue to translate the codons past it. This produces a longer version of the protein, a fusion of the original in-frame amino acids, followed by all the alternate frame amino acids. There is currently no automated software to predict the occurrence of these programmed ribosomal frameshifts (PRF), and they are currently only identified by manual curation.

Results

Here we present PRFect, an innovative machine-learning method for the detection and prediction of PRFs in coding genes of various types. PRFect combines advanced machine learning techniques with the integration of multiple complex cellular properties, such as secondary structure, codon usage, ribosomal binding site interference, direction, and slippery site motif. Calculating and incorporating these diverse properties posed significant challenges, but through extensive research and development, we have achieved a user-friendly approach. The code for PRFect is freely available, open-source, and can be easily installed via a single command in the terminal. Our comprehensive evaluations on diverse organisms, including bacteria, archaea, and phages, demonstrate PRFect's strong performance, achieving high sensitivity, specificity, and an accuracy exceeding 90%.

Conclusion

PRFect represents a significant advancement in the field of PRF detection and prediction, offering a powerful tool for researchers and scientists to unravel the intricacies of programmed ribosomal frameshifting in coding genes.

Characterization and genome analysis of G1 sub-cluster mycobacteriophage Lang

Phage therapy is revitalized as an alternative to antibiotics therapy against antimicrobials resistant pathogens. Mycobacteriophages are genetically diverse viruses that can specifically infect Mycobacterium genus including Mycobacterium tuberculosis and Mycobacterium smegmatis. Here, we isolated and annotated the genome of a mycobacteriophage Lang, a temperate mycobacteriophage isolated from the soil of Hohhot, Inner Mongolia, China, by using Mycolicibacterium smegmatis mc² 155 as the host. It belongs to the Siphoviridae family of Caudovirales as determined by transmission electron microscopy. The morphological characteristics and certain biological properties of the phage were considered in detail. Phage Lang genomes is 41,487 bp in length with 66.85% GC content and encodes 60 putative open reading frames and belongs to the G1 sub-cluster. Genome annotation indicated that genes for structure proteins, assembly proteins, replications/transcription and lysis of the host are present in function clucters. The genome sequence of phage Lang is more than 95% similar to that of mycobacteriophage Grizzly and Sweets, differing in substitutions, insertions and deletions in Lang. One-step growth curve revealed that Lang has a latent period of 30 min and a outbreak period of 90 min. The short latent period and rapid outbreak mark the unique properties of phage Lang, which can be another potential source for combating M. tuberculosis.

Genome Sequence of a Cluster CR2 Gordonia terrae Phage, StarStruck

Bacteriophage StarStruck is a lytic Siphoviridae phage that infects Gordonia terrae 3612. The 68,128-bp genome of StarStruck has a GC content of 65.4% and contains 92 protein-coding genes, including the gene for a HicA-like toxin. StarStruck was assigned to subcluster CR2 based on >35% shared gene content with other cluster CR genomes in the Actinobacteriophage Database.

Genome Sequence of a Cluster DN1 Gordonia terrae Phage, Periwinkle

Periwinkle is a temperate bacteriophage that was isolated on the host Gordonia terrae 3612. The genome has a length of 55,657 bp and a GC content of 62.9% and contains 109 protein-coding genes and no tRNA genes. An 8-kb region after the structural protein genes encodes eight membrane proteins, a tyrosine integrase, and an immunity repressor.

Genome Sequence of Singleton Gordonia terrae Bacteriophage Finkle

Bacteriophage Finkle is a temperate siphovirus isolated from soil on Gordonia terrae. The 47,895-bp genome has a GC content of 66.6% and encodes 84 protein-coding genes. The genome is not closely related to sequences in the Actinobacteriophage database, sharing less than 35% gene content, and was classified as a singleton.

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.

A monomeric mycobacteriophage immunity repressor utilizes two domains to recognize an asymmetric DNA sequence

Regulation of bacteriophage gene expression involves repressor proteins that bind and downregulate early lytic promoters. A large group of mycobacteriophages code for repressors that are unusual in also terminating transcription elongation at numerous binding sites (stoperators) distributed across the phage genome. Here we provide the X-ray crystal structure of a mycobacteriophage immunity repressor bound to DNA, which reveals the binding of a monomer to an asymmetric DNA sequence using two independent DNA binding domains. The structure is supported by small-angle X-ray scattering, DNA binding, molecular dynamics, and in vivo immunity assays. We propose a model for how dual DNA binding domains facilitate regulation of both transcription initiation and elongation, while enabling evolution of other superinfection immune specificities.

Bacteriophage repressor proteins downregulate viral lytic gene expression. Herein, the authors present the X-ray crystal structure of a monomeric repressor that binds an asymmetric DNA sequence using two independent domains.

A Review on Mycobacteriophages: From Classification to Applications

Mycobacterial infections are a group of life-threatening conditions triggered by fast- or slow-growing mycobacteria. Some mycobacteria, such as Mycobacterium tuberculosis, promote the deaths of millions of lives throughout the world annually. The control of mycobacterial infections is influenced by the challenges faced in the diagnosis of these bacteria and the capability of these pathogens to develop resistance against common antibiotics. Detection of mycobacterial infections is always demanding due to the intracellular nature of these pathogens that, along with the lipid-enriched structure of the cell wall, complicates the access to the internal contents of mycobacterial cells. Moreover, recent studies depicted that more than 20% of M. tuberculosis (Mtb) infections are multi-drug resistant (MDR), and only 50% of positive MDR-Mtb cases are responsive to standard treatments. Similarly, the susceptibility of nontuberculosis mycobacteria (NTM) to first-line tuberculosis antibiotics has also declined in recent years. Exploiting mycobacteriophages as viruses that infect mycobacteria has significantly accelerated the diagnosis and treatment of mycobacterial infections. This is because mycobacteriophages, regardless of their cycle type (temperate/lytic), can tackle barriers in the mycobacterial cell wall and make the infected bacteria replicate phage DNA along with their DNA. Although the infectivity of the majority of discovered mycobacteriophages has been evaluated in non-pathogenic M. smegmatis, more research is still ongoing to find mycobacteriophages specific to pathogenic mycobacteria, such as phage DS6A, which has been shown to be able to infect members of the M. tuberculosis complex. Accordingly, this review aimed to introduce some potential mycobacteriophages in the research, specifically those that are infective to the three troublesome mycobacteria, M. tuberculosis, M. avium subsp. paratuberculosis (MAP), and M. abscessus, highlighting their theranostic applications in medicine.

Hidden diversity of double-stranded DNA phages in symbiotic Rhizobium species

In this study, we addressed the extent of diversification of phages associated with nitrogen-fixing symbiotic Rhizobium species. Despite the ecological and economic importance of the Rhizobium genus, little is known about the diversity of the associated phages. A thorough assessment of viral diversity requires investigating both lytic phages and prophages harboured in diverse Rhizobium genomes. Protein-sharing networks identified 56 viral clusters (VCs) among a set of 425 isolated phages and predicted prophages. The VCs formed by phages had more proteins in common and a higher degree of synteny, and they group together in clades in the associated phylogenetic tree. By contrast, the VCs of prophages showed significant genetic variation and gene loss, with selective pressure on the remaining genes. Some VCs were found in various Rhizobium species and geographical locations, suggesting that they have wide host ranges. Our results indicate that the VCs represent distinct taxonomic units, probably representing taxa equivalent to genera or even species. The finding of previously undescribed phage taxa indicates the need for further exploration of the diversity of phages associated with Rhizobium species. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

Phylogenetic relationships and codon usage bias amongst cluster K mycobacteriophages

Bacteriophages infecting pathogenic hosts play an important role in medical research, not only as potential treatments for antibiotic-resistant infections but also offering novel insights into pathogen genetics and evolution. A prominent example is cluster K mycobacteriophages infecting Mycobacterium tuberculosis, a causative agent of tuberculosis in humans. However, as handling M. tuberculosis as well as other pathogens in a laboratory remains challenging, alternative nonpathogenic relatives, such as Mycobacterium smegmatis, are frequently used as surrogates to discover therapeutically relevant bacteriophages in a safer environment. Consequently, the individual host ranges of the majority of cluster K mycobacteriophages identified to date remain poorly understood. Here, we characterized the complete genome of Stinson, a temperate subcluster K1 mycobacteriophage with a siphoviral morphology. A series of comparative genomic analyses revealed strong similarities with other cluster K mycobacteriophages, including the conservation of an immunity repressor gene and a toxin/antitoxin gene pair. Patterns of codon usage bias across the cluster offered important insights into putative host ranges in nature, highlighting that although all cluster K mycobacteriophages are able to infect M. tuberculosis, they are less likely to have shared an evolutionary infection history with Mycobacterium leprae (underlying leprosy) compared to the rest of the genus’ host species. Moreover, subcluster K1 mycobacteriophages are able to integrate into the genomes of Mycobacterium abscessus and Mycobacterium marinum—two bacteria causing pulmonary and cutaneous infections which are often difficult to treat due to their drug resistance.

Wildy Prize Lecture, 2020-2021: Who wouldn't want to discover a new virus?

Innovations in science education are desperately needed to find ways to engage and interest students early in their undergraduate careers. Exposing students to authentic research experiences is highly beneficial, but finding ways to include all types of students and to do this at large scale is especially challenging. An attractive solution is the concept of an inclusive research education community (iREC) in which centralized research leadership and administration supports multiple institutions, including diverse groups of schools and universities, faculty and students. The Science Education Alliance Phage Hunters Advancing Genomics and Evolutionary Sciences (SEA-PHAGES) programme is an excellent example of an iREC, in which students explore viral diversity and evolution through discovery and genomic analysis of novel bacteriophages. The SEA-PHAGES programme has proven to be sustainable, to be implemented at large scale, and to enhance student persistence in science, as well as to produce substantial research advances. Discovering a new virus with the potential for new biological insights and clinical applications is inherently exciting. Who wouldn't want to discover a new virus?

Isolation and Characterisation of the Bundooravirus Genus and Phylogenetic Investigation of the Salasmaviridae Bacteriophages

Bacillus is a highly diverse genus containing over 200 species that can be problematic in both industrial and medical settings. This is mainly attributed to Bacillus sp. being intrinsically resistant to an array of antimicrobial compounds, hence alternative treatment options are needed. In this study, two bacteriophages, PumA1 and PumA2 were isolated and characterized. Genome nucleotide analysis identified the two phages as novel at the DNA sequence level but contained proteins similar to phi29 and other related phages. Whole genome phylogenetic investigation of 34 phi29-like phages resulted in the formation of seven clusters that aligned with recent ICTV classifications. PumA1 and PumA2 share high genetic mosaicism and form a genus with another phage named WhyPhy, more recently isolated from the United States of America. The three phages within this cluster are the only candidates to infect B. pumilus. Sequence analysis of B. pumilus phage resistant mutants revealed that PumA1 and PumA2 require polymerized and peptidoglycan bound wall teichoic acid (WTA) for their infection. Bacteriophage classification is continuously evolving with the increasing phages' sequences in public databases. Understanding phage evolution by utilizing a combination of phylogenetic approaches provides invaluable information as phages become legitimate alternatives in both human health and industrial processes.

Actinobacteriophages: Genomics, Dynamics, and Applications

Actinobacteriophages are viruses that infect bacterial hosts in the phylum Actinobacteria. More than 17,000 actinobacteriophages have been described and over 3,000 complete genome sequences reported, resulting from large-scale, high-impact, integrated research-education initiatives such as the Science Education Alliance Phage Hunters Advancing Genomics and Evolutionary Sciences (SEA-PHAGES) program. Their genomic diversity is enormous; actinobacteriophages comprise many architecturally mosaic genomes with distinct DNA sequences. Their genome diversity is driven by the highly dynamic interactions between phages and their hosts, and prophages can confer a variety of systems that defend against attack by genetically distinct phages; phages can neutralize these defense systems by coding for counter-defense proteins. These phages not only provide insights into diverse and dynamic phage populations but also have provided numerous tools for mycobacterial genetics. A case study using a three-phage cocktail to treat a patient with a drug-resistant Mycobacterium abscessus suggests that phages may have considerable potential for the therapeutic treatment of mycobacterial infections.

Phage Therapy for Mycobacterium Abscessus and Strategies to Improve Outcomes

Members of Mycobacterium abscessus complex are known for causing severe, chronic infections. Members of M. abscessus are a new "antibiotic nightmare" as one of the most resistant organisms to chemotherapeutic agents. Treatment of these infections is challenging due to the either intrinsic or acquired resistance of the M. abscessus complex to the available antibiotics. Recently, successful phage therapy with a cocktail of three phages (one natural lytic phage and two engineered phages) every 12 h for at least 32 weeks has been reported against a severe case of the disseminated M. abscessus subsp. massiliense infection, which underlines the high value of phages against drug-resistant superbugs. This report also highlighted the limitations of phage therapy, such as the absence of lytic phages with a broad host-range against all strains and subspecies of the M. abscessus complex and also the risk of phage resistant bacteria over treatment. Cutting-edge genomic technologies have facilitated the development of engineered phages for therapeutic purposes by introducing new desirable properties, changing host-range and arming the phages with additional killing genes. Here, we review the available literature and suggest new potential solutions based on the progress in phage engineering that can help to overcome the present limitations of M. abscessus treatment.

Mycobacteriophages as Potential Therapeutic Agents against Drug-Resistant Tuberculosis

The current emergence of multi-, extensively-, extremely-, and total-drug resistant strains of Mycobacterium tuberculosis poses a major health, social, and economic threat, and stresses the need to develop new therapeutic strategies. The notion of phage therapy against bacteria has been around for more than a century and, although its implementation was abandoned after the introduction of drugs, it is now making a comeback and gaining renewed interest in Western medicine as an alternative to treat drug-resistant pathogens. Mycobacteriophages are genetically diverse viruses that specifically infect mycobacterial hosts, including members of the M. tuberculosis complex. This review describes general features of mycobacteriophages and their mechanisms of killing M. tuberculosis, as well as their advantages and limitations as therapeutic and prophylactic agents against drug-resistant M. tuberculosis strains. This review also discusses the role of human lung micro-environments in shaping the availability of mycobacteriophage receptors on the M. tuberculosis cell envelope surface, the risk of potential development of bacterial resistance to mycobacteriophages, and the interactions with the mammalian host immune system. Finally, it summarizes the knowledge gaps and defines key questions to be addressed regarding the clinical application of phage therapy for the treatment of drug-resistant tuberculosis.

Mycobacteriophage SWU1-Functionalized magnetic particles for facile bioluminescent detection of Mycobacterium smegmatis

Mycobacterium tuberculosis (M. tuberculosis), the causative agent of tuberculosis, ranks one of the most dangerous pathogens for its large deaths toll. Due to its characteristic extremely slow growth, the conventional culture-based protocol cannot meet the requirement for the efficient diagnosis of M. tuberculosis-induced tuberculosis. With our previously isolated mycobacteriophage SWU1, we tried to develop a mycobacteriophage-based protocol for detecting Mycobacterium genus. In this work, Mycobacterium smegmatis (M. smegmatis) was used as a model due to its similar physiological features as pathogenic M. tuberculosis, much faster growth and nonpathogenic property. Mycobacteriophage SWU1-functionalized magnetic particles (SWU1-MPs) were used as highly efficient separation carriers for the viable host M. smegmatis. After a replication cycle of approximate 60 min, the cells of M. smegmatis were disrupted by the progeny mycobacteriophages to release intracellular adenosine triphosphate (ATP). The bioluminescent (BL) signal of released ATP was collected to quantitate the amount of M. smegmatis. For the developed protocol, the detection range is 5.0 × 10² to 5.0 × 10⁵ CFU mL^-1, and the detection limit is 3.8 × 10² CFU mL^-1 (S/N = 3). Furthermore, the protocol can exclude the potential interference of 3 non-pathogenic mycobacteria and 6 other bacterial species. It has been successfully applied to quantitate M. smegmatis in human urine, human saliva, and human serum. The results demonstrate its application potential for a simple, fast, and specific diagnosis of M. tuberculosis infection.

Genome Sequences of 12 Phages That Infect Klebsiella pneumoniae

Klebsiella pneumoniae is a pathogen responsible for significant proportions of nosocomial and health care-associated infections and is known to acquire multiple antibiotic resistance genes. Here, we announce the full genome sequences of 12 K. pneumoniae bacteriophages from samples collected in wastewater treatment facilities across the western United States.

Klebsiella pneumoniae is a pathogen responsible for significant proportions of nosocomial and health care-associated infections and is known to acquire multiple antibiotic resistance genes. Here, we announce the full genome sequences of 12 K. pneumoniae bacteriophages from samples collected in wastewater treatment facilities across the western United States.

Characterization and genome analysis of B1 sub-cluster mycobacteriophage PDRPxv

Mycobacteriophages are viruses specific to mycobacteria that have gained attention as alternative therapeutic strategies for treating antibiotic-resistant infections. Mycobacteriophages are highly diverse and have been grouped into 29 clusters, 71 sub-clusters and 10 singletons based on the genome sequence. Here, we annotate the genome of PDRPxv, a lytic mycobacteriophage isolated from New Delhi; it belongs to the Siphoviridae family as determined by transmission electron microscopy. This phage survives at higher temperatures (up to 55 °C) and in alkaline conditions (up to pH11). PDRPxv phage genome is 69,171 bp in length with 66.35 % GC content and encodes 107 putative open reading frames and belongs to the B1 sub-cluster. Genome annotation indicated that genes for DNA encapsidation, structural proteins, replication/transcription and lysis of the host are present in functional clusters. Structural proteins encoded by Gp10-Gp12, Gp18, Gp25 and Gp28-Gp33 were identified by mass spectrometry. Interestingly, no gene encoding a holin function was found. Single-step growth curve revealed that PDRPxv has an adsorption time of 45 min, a latency time of 135 min and an average burst size of 99 phage particles per infected cell. The short latency period and the large burst size mark the lytic nature of the PDRPxv phage, which could therefore be a promising therapeutic candidate against pathogenic Mycobacterium species.

Antimicrobial activity of Mycobacteriophage D29 Lysin B during Mycobacterium ulcerans infection

Buruli Ulcer (BU) is a cutaneous disease caused by Mycobacterium ulcerans. The pathogenesis of this disease is closely related to the secretion of the toxin mycolactone that induces extensive destruction of the skin and soft tissues. Currently, there are no effective measures to prevent the disease and, despite availability of antibiotherapy and surgical treatments, these therapeutic options are often associated with severe side effects. Therefore, it is important to develop alternative strategies for the treatment of BU. Endolysins (lysins) are phage encoded enzymes that degrade peptidoglycan of bacterial cell walls. Over the past years, lysins have been emerging as alternative antimicrobial agents against bacterial infections. However, mycobacteria have an unusual outer membrane composed of mycolylarabinogalactan-peptidoglycan. To overcome this complex barrier, some mycobacteriophages encode a lipolytic enzyme, Lysin B (LysB). In this study, we demonstrate for the first time that recombinant LysB displays lytic activity against M. ulcerans isolates. Moreover, using a mouse model of M. ulcerans footpad infection, we show that subcutaneous treatment with LysB prevented further bacterial proliferation, associated with IFN-γ and TNF production in the draining lymph node. These findings highlight the potential use of lysins as a novel therapeutic approach against this neglected tropical disease.

Buruli Ulcer (BU) is a necrotizing skin disease caused by Mycobacterium ulcerans. Although the current antibiotic treatment for BU is effective, daily administrations for a prolonged period of time, combined with potential risk of severe side effects, negatively impact on patient adherence. In that sense, we tested the efficacy of an alternative strategy based on Lysin B (LysB), a phage encoded lipolytic enzyme that degrades the mycolylarabinogalactan-peptidoglycan complex present in the mycobacterial cell wall. In this study, we show that LysB not only displays lytic activity against M. ulcerans isolates in vitro, but also leads to a decrease of M. ulcerans proliferation in infected mouse footpads. These findings highlight the potential use of lysins as a novel therapeutic approach against this neglected tropical disease.

Evolution of Superinfection Immunity in Cluster A Mycobacteriophages

Many aspects regarding superinfection, immunity, virulence, and the evolution of immune specificities are poorly understood due to the lack of large collections of isolated and sequenced phages with a spectrum of genetic diversity. Using a genetically diverse collection of Cluster A phages, we show that the classical and relatively straightforward patterns of homoimmunity, heteroimmunity, and virulence result from interactions between homotypic and heterotypic phages at the extreme edges of an evolutionary continuum of immune specificities. Genetic interactions between mesotypic phages result in more complex mesoimmunity phenotypes and virulence profiles. These results highlight that the evolution of immune specificities can be shaped by homotypic and mesotypic interactions and may be more dynamic than previously considered.

Temperate phages encode an immunity system to control lytic gene expression during lysogeny. This gene regulatory circuit consists of multiple interacting genetic elements, and although it is essential for controlling phage growth, it is subject to conflicting evolutionary pressures. During superinfection of a lysogen, the prophage’s circuit interacts with the superinfecting phage’s circuit and prevents lytic growth if the two circuits are closely related. The circuitry is advantageous since it provides the prophage with a defense mechanism, but the circuitry is also disadvantageous since it limits the phage’s host range during superinfection. Evolutionarily related phages have divergent, orthogonal immunity systems that no longer interact and are heteroimmune, but we do not understand how immunity systems evolve new specificities. Here, we use a group of Cluster A mycobacteriophages that exhibit a spectrum of genetic diversity to examine how immunity system evolution impacts superinfection immunity. We show that phages with mesotypic (i.e., genetically related but distinct) immunity systems exhibit asymmetric and incomplete superinfection phenotypes. They form complex immunity networks instead of well-defined immunity groups, and mutations conferring escape (i.e., virulence) from homotypic or mesotypic immunity have various escape specificities. Thus, virulence and the evolution of new immune specificities are shaped by interactions with homotypic and mesotypic immunity systems.

A mycobacteriophage genomics approach to identify novel mycobacteriophage proteins with mycobactericidal properties

Mycobacteriophages that are specific to mycobacteria are sources of various effector proteins that are capable of eliciting bactericidal responses. We describe a genomics approach in combination with bioinformatics to identify mycobacteriophage proteins that are toxic to mycobacteria upon expression. A genomic library comprising phage genome collections was screened for clones capable of killing Mycobacterium smegmatis strain mc²155. We identified four unique clones: clones 45 and 12N (from the mycobacteriophage D29) and clones 66 and 85 (from the mycobacteriophage Che12). The gene products from clones 66 and 45 were identified as Gp49 of the Che12 phage and Gp34 of the D29 phage, respectively. The gene products of the other two clones, 85 and 12N, utilized novel open reading frames (ORFs) coding for synthetic proteins. These four clones (clones 45, 66, 85 and 12N) caused growth defects in M. smegmatis and Mycobacterium bovis upon expression. Clones with Gp49 and Gp34 also induced growth defects in Escherichia coli, indicating that they target conserved host machineries. Their expression induced various morphological changes, indicating that they affected DNA replication and cell division steps. We predicted that Gp34 is a Xis protein that is required in phage DNA excision from the bacterial chromosome. Gp49 is predicted to have an HTH motif with DNA-bending/twisting properties. We suggest that this methodology is useful to identify new phage proteins with the desired properties without laboriously characterizing the individual phages. It is universal and could be applied to other bacteria-phage systems. We speculate that the existence of a virtually unlimited number of phages with unique gene products could offer a cheaper and less hazardous alternative to explore new antimicrobial molecules.

Mycobacteriophage CRB2 defines a new subcluster in mycobacteriophage classification

Mycobacteriophages are viruses -mostly temperates- that infect Mycobacterium smegmatis and sometimes Mycobacterium tuberculosis. Mycobacteriophages are grouped in clusters on the basis of the overall nucleotide sequence homology, being further divided in subclusters as more mycobacteriophage genomes are sequenced and annotated. As part of our on-going screening for novel isolates, we herein report the bioinformatics analysis of CRB2, a mycobacteriophage belonging into the Siphoviridae family that propagates at 30°C. CRB2 has a 72,217 bp genome with a 69.78% GC content that belongs to Cluster B; nucleotide comparison with other B cluster members positions CRB2 as the sole member of a new subcluster, B9, being mycobacteriophage Saguaro (belonging into subcluster B7) its closest relative. Sequencing and annotation of 14 mycobacteriophages isolated by our group has yielded six cluster A members, a singleton, four of the five members of subcluster B6, one of the three reported members of subcluster G4, and CRB2 which defines subcluster B9. Considering the massive mycobacteriophage search performed in USA and the relatively rarity of our phages, we propose that factors other than size of the sampling determine the variability of mycobacteriophage distribution, and thus a world-wide concerted mining would most likely bring extremely rare and yet undiscovered mycobacteriophages.

Mycobacteriophages

Mycobacteriophages are viruses that infect mycobacterial hosts. A large number of mycobacteriophages have been isolated and genomically characterized, providing insights into viral diversity and evolution, as well as fueling development of tools for mycobacterial genetics. Mycobacteriophages have intimate relationships with their hosts and provide insights into the genetics and physiology of the mycobacteria and tools for potential clinical applications such as drug development, diagnosis, vaccines, and potentially therapy.

Characterization and induction of prophages in human gut-associated Bifidobacterium hosts

In the current report, we describe the identification of three genetically distinct groups of prophages integrated into three different chromosomal sites of human gut-associated Bifidobacterium breve and Bifidobacterium longum strains. These bifidobacterial prophages are distantly related to temperate actinobacteriophages of several hosts. Some prophages, integrated within the dnaJ₂ gene, are competent for induction, excision, replication, assembly and lysis, suggesting that they are fully functional and can generate infectious particles, even though permissive hosts have not yet been identified. Interestingly, several of these phages harbor a putative phase variation shufflon (the Rin system) that generates variation of the tail-associated receptor binding protein (RBP). Unlike the analogous coliphage-associated shufflon Min, or simpler Cin and Gin inversion systems, Rin is predicted to use a tyrosine recombinase to promote inversion, the first reported phage-encoded tyrosine-family DNA invertase. The identification of bifidobacterial prophages with RBP diversification systems that are competent for assembly and lysis, yet fail to propagate lytically under laboratory conditions, suggests dynamic evolution of bifidobacteria and their phages in the human gut.

Mycobacteriophage Lysis Enzymes: Targeting the Mycobacterial Cell Envelope

Mycobacteriophages are viruses that specifically infect mycobacteria, which ultimately culminate in host cell death. Dedicated enzymes targeting the complex mycobacterial cell envelope arrangement have been identified in mycobacteriophage genomes, thus being potential candidates as antibacterial agents. These comprise lipolytic enzymes that target the mycolic acid-containing outer membrane and peptidoglycan hydrolases responsive to the atypical mycobacterial peptidoglycan layer. In the recent years, a remarkable progress has been made, particularly on the comprehension of the mechanisms of bacteriophage lysis proteins activity and regulation. Notwithstanding, information about mycobacteriophages lysis strategies is limited and is mainly represented by the studies performed with mycobacteriophage Ms6. Since mycobacteriophages target a specific group of bacteria, which include Mycobacterium tuberculosis responsible for one of the leading causes of death worldwide, exploitation of the use of these lytic enzymes demands a special attention, as they may be an alternative to tackle multidrug resistant tuberculosis. This review focuses on the current knowledge of the function of lysis proteins encoded by mycobacteriophages and their potential applications, which may contribute to increasing the effectiveness of antimycobacterial therapy.

Fluoromycobacteriophages Can Detect Viable Mycobacterium tuberculosis and Determine Phenotypic Rifampicin Resistance in 3-5 Days From Sputum Collection

The World Health Organization (WHO) estimates that 40% of tuberculosis (TB) cases are not diagnosed and treated correctly. Even though there are several diagnostic tests available in the market, rapid, easy, inexpensive detection, and drug susceptibility testing (DST) of Mycobacterium tuberculosis is still of critical importance specially in low and middle-income countries with high incidence of the disease. In this work, we have developed a microscopy-based methodology using the reporter mycobacteriophage mCherry_bombϕ for detection of Mycobacterium spp. and phenotypic determination of rifampicin resistance within just days from sputum sample collection. Fluoromycobacteriophage methodology is compatible with regularly used protocols in clinical laboratories for TB diagnosis and paraformaldehyde fixation after infection reduces biohazard risks with sample analysis by fluorescence microscopy. We have also set up conditions for discrimination between M. tuberculosis complex (MTBC) and non-tuberculous mycobacteria (NTM) strains by addition of p-nitrobenzoic acid (PNB) during the assay. Using clinical isolates of pre-XDR and XDR-TB strains from this study, we tested mCherry_bombΦ for extended DST and we compared the antibiotic resistance profile with those predicted by whole genome sequencing. Our results emphasize the utility of a phenotypic test for M. tuberculosis extended DST. The many attributes of mCherry_bombΦ suggests this could be a useful component of clinical microbiological laboratories for TB diagnosis and since only viable cells are detected this could be a useful tool for monitoring patient response to treatment.

Expression, purification and crystallization of phosphoribosyl transferase from a mycobacteriophage

Tuberculosis (TB) continues to remain a leading cause of death globally. Of particular concern is the emergence and rise in incidence of multidrug-resistant and extremely drug-resistant cases of TB. To counter this threat, it is important to explore alternative therapies, including phage therapy. Phage BTCU-1 specifically infects Mycobacterium spp. and kills the majority of them. Intriguingly, many proteins from the phage do not share high amino-acid sequence identity with proteins from species other than phages. Here, the expression, purification and crystallization of one such protein, a putative phosphoribosyl transferase from phage BTCU-1, is reported. The crystals belonged to space group C2221, with unit-cell parameters a = 59.71, b = 64.42, c = 65.32 Å, α = β = γ = 90°. The crystals diffracted X-rays to 2.2 Å resolution.

Genomic differentiation among wild cyanophages despite widespread horizontal gene transfer

Background

Genetic recombination is a driving force in genome evolution. Among viruses it has a dual role. For genomes with higher fitness, it maintains genome integrity in the face of high mutation rates. Conversely, for genomes with lower fitness, it provides immediate access to sequence space that cannot be reached by mutation alone. Understanding how recombination impacts the cohesion and dissolution of individual whole genomes within viral sequence space is poorly understood across double-stranded DNA bacteriophages (a.k.a phages) due to the challenges of obtaining appropriately scaled genomic datasets.

Results

Here we explore the role of recombination in both maintaining and differentiating whole genomes of 142 wild double-stranded DNA marine cyanophages. Phylogenomic analysis across the 51 core genes revealed ten lineages, six of which were well represented. These phylogenomic lineages represent discrete genotypic populations based on comparisons of intra- and inter- lineage shared gene content, genome-wide average nucleotide identity, as well as detected gaps in the distribution of pairwise differences between genomes. McDonald-Kreitman selection tests identified putative niche-differentiating genes under positive selection that differed across the six well-represented genotypic populations and that may have driven initial divergence. Concurrent with patterns of recombination of discrete populations, recombination analyses of both genic and intergenic regions largely revealed decreased genetic exchange across individual genomes between relative to within populations.

Conclusions

These findings suggest that discrete double-stranded DNA marine cyanophage populations occur in nature and are maintained by patterns of recombination akin to those observed in bacteria, archaea and in sexual eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3286-x) contains supplementary material, which is available to authorized users.

Genetic and life-history traits associated with the distribution of prophages in bacteria

Nearly half of the sequenced bacteria are lysogens and many of their prophages encode adaptive traits. Yet, the variables driving prophage distribution remain undetermined. We identified 2246 prophages in complete bacterial genomes to study the genetic and life-history traits associated with lysogeny. While optimal growth temperatures and average cell volumes were not associated with lysogeny, prophages were more frequent in pathogens and in bacteria with small minimal doubling times. Their frequency also increased with genome size, but only for genomes smaller than 6 Mb. The number of spacers in CRISPR-Cas systems and the frequency of type III systems were anticorrelated with prophage frequency, but lysogens were more likely to encode type I and type II systems. The minimal doubling time was the trait most correlated with lysogeny, followed by genome size and pathogenicity. We propose that bacteria with highly variable growth rates often encounter lower opportunity costs for lysogeny relative to lysis. These results contribute to explain the paucity of temperate phages in certain bacterial clades and of bacterial lysogens in certain environments. They suggest that genetic and life-history traits affect the contributions of temperate phages to bacterial genomes.

The Double-Stranded DNA Virosphere as a Modular Hierarchical Network of Gene Sharing

Virus genomes are prone to extensive gene loss, gain, and exchange and share no universal genes. Therefore, in a broad-scale study of virus evolution, gene and genome network analyses can complement traditional phylogenetics. We performed an exhaustive comparative analysis of the genomes of double-stranded DNA (dsDNA) viruses by using the bipartite network approach and found a robust hierarchical modularity in the dsDNA virosphere. Bipartite networks consist of two classes of nodes, with nodes in one class, in this case genomes, being connected via nodes of the second class, in this case genes. Such a network can be partitioned into modules that combine nodes from both classes. The bipartite network of dsDNA viruses includes 19 modules that form 5 major and 3 minor supermodules. Of these modules, 11 include tailed bacteriophages, reflecting the diversity of this largest group of viruses. The module analysis quantitatively validates and refines previously proposed nontrivial evolutionary relationships. An expansive supermodule combines the large and giant viruses of the putative order "Megavirales" with diverse moderate-sized viruses and related mobile elements. All viruses in this supermodule share a distinct morphogenetic tool kit with a double jelly roll major capsid protein. Herpesviruses and tailed bacteriophages comprise another supermodule, held together by a distinct set of morphogenetic proteins centered on the HK97-like major capsid protein. Together, these two supermodules cover the great majority of currently known dsDNA viruses. We formally identify a set of 14 viral hallmark genes that comprise the hubs of the network and account for most of the intermodule connections.

IMPORTANCE

Viruses and related mobile genetic elements are the dominant biological entities on earth, but their evolution is not sufficiently understood and their classification is not adequately developed. The key reason is the characteristic high rate of virus evolution that involves not only sequence change but also extensive gene loss, gain, and exchange. Therefore, in the study of virus evolution on a large scale, traditional phylogenetic approaches have limited applicability and have to be complemented by gene and genome network analyses. We applied state-of-the art methods of such analysis to reveal robust hierarchical modularity in the genomes of double-stranded DNA viruses. Some of the identified modules combine highly diverse viruses infecting bacteria, archaea, and eukaryotes, in support of previous hypotheses on direct evolutionary relationships between viruses from the three domains of cellular life. We formally identify a set of 14 viral hallmark genes that hold together the genomic network.

Homology modeling, substrate docking, and molecular simulation studies of mycobacteriophage Che12 lysin A

Mycobacteriophages produce lysins that break down the host cell wall at the end of lytic cycle to release their progenies. The ability to lyse mycobacterial cells makes the lysins significant. Mycobacteriophage Che12 is the first reported temperate phage capable of infecting and lysogenising Mycobacterium tuberculosis. Gp11 of Che12 was found to have Chitinase domain that serves as endolysin (lysin A) for Che12. Structure of gp11 was modeled and evaluated using Ramachandran plot in which 98 % of the residues are in the favored and allowed regions. Che12 lysin A was predicted to act on NAG-NAM-NAG molecules in the peptidoglycan of cell wall. The tautomers of NAG-NAM-NAG molecule were generated and docked with lysin A. The stability and binding affinity of lysin A - NAG-NAM-NAG tautomers were studied using molecular dynamics simulations.

Mycobacteriophage putative GTPase-activating protein can potentiate antibiotics

The soaring incidences of infection by antimicrobial resistant (AR) pathogens and shortage of effective antibiotics with new mechanisms of action have renewed interest in phage therapy. This scenario is exemplified by resistant tuberculosis (TB), caused by resistant Mycobacterium tuberculosis. Mycobacteriophage SWU1 A321_gp67 encodes a putative GTPase-activating protein. Mycobacterium smegmatis with gp67 overexpression showed changed colony formation and biofilm morphology and supports the efficacy of streptomycin and capreomycin against Mycobacterium. gp67 down-regulated the transcription of genes involved in cell wall and biofilm development. To our knowledge, this is the first report to show that phage protein in addition to lysin or recombination components can synergize with existing antibiotics. Phage components might represent a promising new clue for better antibiotic potentiators.

Characterization of mycobacteria and mycobacteriophages isolated from compost at the São Paulo Zoo Park Foundation in Brazil and creation of the new mycobacteriophage Cluster U

Background

A large collection of sequenced mycobacteriophages capable of infecting a single host strain of Mycobacterium smegmatis shows considerable genomic diversity with dozens of distinctive types (clusters) and extensive variation within those sharing evident nucleotide sequence similarity. Here we profiled the mycobacterial components of a large composting system at the São Paulo zoo.

Results

We isolated and sequenced eight mycobacteriophages using Mycobacterium smegmatis mc²155 as a host. None of these eight phages infected any of mycobacterial strains isolated from the same materials. The phage isolates span considerable genomic diversity, including two phages (Barriga, Nhonho) related to Subcluster A1 phages, two Cluster B phages (Pops, Subcluster B1; Godines, Subcluster B2), three Subcluster F1 phages (Florinda, Girafales, and Quico), and Madruga, a relative of phage Patience with which it constitutes the new Cluster U. Interestingly, the two Subcluster A1 phages and the three Subcluster F1 phages have genomic relationships indicating relatively recent evolution within a geographically isolated niche in the composting system.

Conclusions

We predict that composting systems such as those used to obtain these mycobacteriophages will be a rich source for the isolation of additional phages that will expand our view of bacteriophage diversity and evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-016-0734-3) contains supplementary material, which is available to authorized users.

Function, expression, specificity, diversity and incompatibility of actinobacteriophage parABS systems

More than 180 individual phages infecting hosts in the phylum Actinobacteria have been sequenced and grouped into Cluster A because of their similar overall nucleotide sequences and genome architectures. These Cluster A phages are either temperate or derivatives of temperate parents, and most have an integration cassette near the centre of the genome containing an integrase gene and attP. However, about 20% of the phages lack an integration cassette, which is replaced by a 1.4 kbp segment with predicted partitioning functions, including plasmid-like parA and parB genes. Phage RedRock forms stable lysogens in Mycobacterium smegmatis in which the prophage replicates at 2.4 copies/chromosome and the partitioning system confers prophage maintenance. The parAB genes are expressed upon RedRock infection of M. smegmatis, but are downregulated once lysogeny is established by binding of RedRock ParB to parS-L, one of two centromere-like sites flanking the parAB genes. The RedRock parS-L and parS-R sites are composed of eight directly repeated copies of an 8 bp motif that is recognized by ParB. The actinobacteriophage parABS cassettes span considerable sequence diversity and specificity, providing a suite of tools for use in mycobacterial genetics.

Vitamin C induced DevR-dependent synchronization of Mycobacterium smegmatis growth and its effect on the proliferation of mycobacteriophage D29

Vitamin C is known to inhibit mycobacterial growth by acting as a hypoxia inducing agent. While investigating how mycobacteriophage growth is influenced by hypoxic conditions induced by vitamin C, using Mycobacterium smegmatis- mycobacteriophage D29 as a model system, it was observed that prior exposure of the host to such conditions resulted in increased burst size of the phage. Vitamin C pre-exposure was also found to induce synchronous growth of the host. A mutant defective in DevR, the response regulator that controls hypoxic responses in mycobacteria, neither supported higher phage bursts nor was it able to undergo synchronized growth following vitamin C pre-exposure, indicating thereby that the two phenomena are interrelated. Further evidence supporting such an interrelationship was obtained from the observation that phage burst sizes varied depending on the stage of synchronous growth that the host cells were in, at the time of infection-higher bursts were observed in the resting/synthetic phases and lower in the dividing ones. The effects were specific in nature as synchronization by an unrelated method, known as 'crowding', did not lead to the same consequence. The results indicate that growth synchronization induced by vitamin C treatment is a DevR-dependent phenomenon which is exploited by mycobacteriophage D29 to grow in larger numbers.

Genetics of Phage Lysis

We have been witnessing an increased interest in bacteriophage studies focused on their use as antibacterial agents to fight pathogenic bacteria. This interest is a consequence of the phages' ability to lyse a bacterial host. Until recently, little was known about the mechanisms used by mycobacteriophages to induce lysis of their complex hosts. However, studies on Ms6-induced lysis have changed this scenario and provided new insights into the mechanisms of bacteriophage-induced lysis. Specific lysis protein genes have been identified in mycobacteriophage genomes, reflecting the particular mycobacterial cell envelope composition. These include enzymes that target mycolic acid-containing lipids and proteins that participate in the secretion of the phage endolysin, functioning as chaperone-like proteins. This chapter focuses on the current knowledge of mycobacteriophage-induced lysis, starting with an overview of phage lysis and basic features of the lysis players.

Dynamics of Mycobacteriophage-Mycobacterial Host Interaction: Evidence for Secondary Mechanisms for Host Lethality

Mycobacteriophages infect mycobacteria, resulting in their death. Therefore, the possibility of using them as therapeutic agents against the deadly mycobacterial disease tuberculosis (TB) is of great interest. To obtain better insight into the dynamics of mycobacterial inactivation by mycobacteriophages, this study was initiated using mycobacteriophage D29 and Mycobacterium smegmatis as the phage-host system. Here, we implemented a goal-oriented iterative cycle of experiments on one hand and mathematical modeling combined with Monte Carlo simulations on the other. This integrative approach lends valuable insight into the detailed kinetics of bacterium-phage interactions. We measured time-dependent changes in host viability during the growth of phage D29 in M. smegmatis at different multiplicities of infection (MOI). The predictions emerging out of theoretical analyses were further examined using biochemical and cell biological assays. In a phage-host interaction system where multiple rounds of infection are allowed to take place, cell counts drop more rapidly than expected if cell lysis is considered the only mechanism for cell death. The phenomenon could be explained by considering a secondary factor for cell death in addition to lysis. Further investigations reveal that phage infection leads to the increased production of superoxide radicals, which appears to be the secondary factor. Therefore, mycobacteriophage D29 can function as an effective antimycobacterial agent, the killing potential of which may be amplified through secondary mechanisms.

Bacteriophages of wastewater foaming-associated filamentous Gordonia reduce host levels in raw activated sludge

Filamentous bacteria are a normal and necessary component of the activated sludge wastewater treatment process, but the overgrowth of filamentous bacteria results in foaming and bulking associated disruptions. Bacteriophages, or phages, were investigated for their potential to reduce the titer of foaming bacteria in a mixed-microbial activated sludge matrix. Foaming-associated filamentous bacteria were isolated from activated sludge of a commercial wastewater treatment plan and identified as Gordonia species by 16S rDNA sequencing. Four representative phages were isolated that target G. malaquae and two un-named Gordonia species isolates. Electron microscopy revealed the phages to be siphophages with long tails. Three of the phages - GordTnk2, Gmala1, and GordDuk1 - had very similar ~76 kb genomes, with >93% DNA identity. These genomes shared limited synteny with Rhodococcus equi phage ReqiDocB7 and Gordonia phage GTE7. In contrast, the genome of phage Gsput1 was smaller (43 kb) and was not similar enough to any known phage to be placed within an established phage type. Application of these four phages at MOIs of 5–15 significantly reduced Gordonia host levels in a wastewater sludge model by approximately 10-fold as compared to non-phage treated reactors. Phage control was observed for nine days after treatment.

Biochemical Characterization of a Mycobacteriophage Derived DnaB Ortholog Reveals New Insight into the Evolutionary Origin of DnaB Helicases

The bacterial replicative helicases known as DnaB are considered to be members of the RecA superfamily. All members of this superfamily, including DnaB, have a conserved C- terminal domain, known as the RecA core. We unearthed a series of mycobacteriophage encoded proteins in which the RecA core domain alone was present. These proteins were phylogenetically related to each other and formed a distinct clade within the RecA superfamily. A mycobacteriophage encoded protein, Wildcat Gp80 that roots deep in the DnaB family, was found to possess a core domain having significant sequence homology (Expect value < 10^-5) with members of this novel cluster. This indicated that Wildcat Gp80, and by extrapolation, other members of the DnaB helicase family, may have evolved from a single domain RecA core polypeptide belonging to this novel group. Biochemical investigations confirmed that Wildcat Gp80 was a helicase. Surprisingly, our investigations also revealed that a thioredoxin tagged truncated version of the protein in which the N-terminal sequences were removed was fully capable of supporting helicase activity, although its ATP dependence properties were different. DnaB helicase activity is thus, primarily a function of the RecA core although additional N-terminal sequences may be necessary for fine tuning its activity and stability. Based on sequence comparison and biochemical studies we propose that DnaB helicases may have evolved from single domain RecA core proteins having helicase activities of their own, through the incorporation of additional N-terminal sequences.

Whole genome comparison of a large collection of mycobacteriophages reveals a continuum of phage genetic diversity

The bacteriophage population is large, dynamic, ancient, and genetically diverse. Limited genomic information shows that phage genomes are mosaic, and the genetic architecture of phage populations remains ill-defined. To understand the population structure of phages infecting a single host strain, we isolated, sequenced, and compared 627 phages of Mycobacterium smegmatis. Their genetic diversity is considerable, and there are 28 distinct genomic types (clusters) with related nucleotide sequences. However, amino acid sequence comparisons show pervasive genomic mosaicism, and quantification of inter-cluster and intra-cluster relatedness reveals a continuum of genetic diversity, albeit with uneven representation of different phages. Furthermore, rarefaction analysis shows that the mycobacteriophage population is not closed, and there is a constant influx of genes from other sources. Phage isolation and analysis was performed by a large consortium of academic institutions, illustrating the substantial benefits of a disseminated, structured program involving large numbers of freshman undergraduates in scientific discovery.

DOI:http://dx.doi.org/10.7554/eLife.06416.001

Viruses are unable to replicate independently. To generate copies of itself, a virus must instead invade a target cell and commandeer that cell's replication machinery. Different viruses are able to invade different types of cell, and a group of viruses known as bacteriophages (or phages for short) replicate within bacteria. The enormous number and diversity of phages in the world means that they play an important role in virtually every ecosystem.

Despite their importance, relatively little is known about how different phage populations are related to each other and how they evolved. Many phages contain their genetic information in the form of strands of DNA. Using genetic sequencing to find out where and how different genes are encoded in the DNA can reveal information about how different viruses are related to each other. These relationships are particularly complicated in phages, as they can exchange genes with other viruses and microbes.

Previous studies comparing the genomes—the complete DNA sequence—of reasonably small numbers of phages that infect the Mycobacterium group of bacteria have found that the phages can be sorted into ‘clusters’ based on similarities in their genes and where these are encoded in their DNA. However, the number of phages investigated so far has been too small to conclude how different clusters are related. Are the clusters separate, or do they form a ‘continuum’ with different genes and DNA sequences shared between different clusters?

Here, Pope, Bowman, Russell et al. compare the individual genomes of 627 bacteriophages that infect the bacterial species Mycobacterium smegmatis. This is by far the largest number of phage genomes analyzed from a single host species. The large number of genomes analyzed allowed a much clearer understanding of the complexity and diversity of these phages to be obtained. The isolation, sequencing and analysis of the hundreds of M. smegmatis bacteriophage genomes was performed by an integrated research and education program, called the Science Education Alliance Phage Hunters Advancing Genomics and Evolutionary Science (SEA-PHAGES) program. This enabled thousands of undergraduate students from different institutions to contribute to the phage discovery and sequencing project, and co-author the report. SEA-PHAGES therefore shows that it is possible to successfully incorporate genuine scientific research into an undergraduate course, and that doing so can benefit both the students and researchers involved.

The results show that while the genomes could be categorized into 28 clusters, the genomes are not completely unrelated. Instead, a spread of diversity is seen, as genes and groups of genes are shared between different clusters. Pope, Bowman, Russell et al. further reveal that the phage population is in a constant state of change, and continuously acquires genes from other microorganisms and viruses.

DOI:http://dx.doi.org/10.7554/eLife.06416.002