Modularity and diversity of target selectors in Tn7 transposons.

To spread, transposons must integrate into target sites without disruption of essential genes while avoiding host defense systems. Tn7-like transposons employ multiple mechanisms for target-site selection, including protein-guided targeting and, in CRISPR-associated transposons (CASTs), RNA-guided targeting. Combining phylogenomic and structural analyses, we conducted a broad survey of target selectors, revealing diverse mechanisms used by Tn7 to recognize target sites, including previously uncharacterized target-selector proteins found in newly discovered transposable elements (TEs). We experimentally characterized a CAST I-D system and a Tn6022-like transposon that uses TnsF, which contains an inactivated tyrosine recombinase domain, to target the comM gene. Additionally, we identified a non-Tn7 transposon, Tsy, encoding a homolog of TnsF with an active tyrosine recombinase domain, which we show also inserts into comM. Our findings show that Tn7 transposons employ modular architecture and co-opt target selectors from various sources to optimize target selection and drive transposon spread.

Alternating lysis and lysogeny is a winning strategy in bacteriophages due to Parrondo's paradox.

SignificanceBacteriophages, the most widespread reproducing biological entity on Earth, employ two strategies of virus-host interaction: lysis of the host cell and lysogeny whereby the virus genome integrates into the host genome and propagates vertically with it. We present a population model that reveals an effect known as Parrondo's paradox in game theory: Alternating between lysis and lysogeny is a winning strategy for a bacteriophage, even when each strategy individually is at a disadvantage compared with a competing bacteriophage. Thus, evolution of bacteriophages appears to optimize the ratio between the lysis and lysogeny propensities rather than the phage burst size in any individual phase. This phenomenon is likely to be relevant for understanding evolution of other host-parasites systems.

Phage lysis-lysogeny switches and programmed cell death: Danse macabre.

Exploration of immune systems in prokaryotes, such as restriction-modification or CRISPR-Cas, shows that both innate and adaptive systems possess programmed cell death (PCD) potential. The key outstanding question is how the immune systems sense and "predict" infection outcomes to "decide" whether to fight the pathogen or induce PCD. There is a striking parallel between this life-or-death decision faced by the cell and the decision by temperate viruses to protect or kill their hosts, epitomized by the lysis-lysogeny switch of bacteriophage Lambda. Immune systems and temperate phages sense the same molecular inputs, primarily, DNA damage, that determine whether the cell lives or dies. Because temperate (pro)phages are themselves components of prokaryotic genomes, their shared "interests" with the hosts result in coregulation of the lysis-lysogeny switch and immune systems that jointly provide the cell with the decision machinery to probe and predict infection outcomes, answering the life-or-death question.

Seeker: alignment-free identification of bacteriophage genomes by deep learning.

Recent advances in metagenomic sequencing have enabled discovery of diverse, distinct microbes and viruses. Bacteriophages, the most abundant biological entity on Earth, evolve rapidly, and therefore, detection of unknown bacteriophages in sequence datasets is a challenge. Most of the existing detection methods rely on sequence similarity to known bacteriophage sequences, impeding the identification and characterization of distinct, highly divergent bacteriophage families. Here we present Seeker, a deep-learning tool for alignment-free identification of phage sequences. Seeker allows rapid detection of phages in sequence datasets and differentiation of phage sequences from bacterial ones, even when those phages exhibit little sequence similarity to established phage families. We comprehensively validate Seeker's ability to identify previously unidentified phages, and employ this method to detect unknown phages, some of which are highly divergent from the known phage families. We provide a web portal (seeker.pythonanywhere.com) and a user-friendly Python package (github.com/gussow/seeker) allowing researchers to easily apply Seeker in metagenomic studies, for the detection of diverse unknown bacteriophages.

Genomic and ecological attributes of marine bacteriophages encoding bacterial virulence genes.

BACKGROUND: Bacteriophages encode genes that modify bacterial functions during infection. The acquisition of phage-encoded virulence genes is a major mechanism for the rise of bacterial pathogens. In coral reefs, high bacterial density and lysogeny has been proposed to exacerbate reef decline through the transfer of phage-encoded virulence genes. However, the functions and distribution of these genes in phage virions on the reef remain unknown. RESULTS: Here, over 28,000 assembled viral genomes from the free viral community in Atlantic and Pacific Ocean coral reefs were queried against a curated database of virulence genes. The diversity of virulence genes encoded in the viral genomes was tested for relationships with host taxonomy and bacterial density in the environment. These analyses showed that bacterial density predicted the profile of virulence genes encoded by phages. The Shannon diversity of virulence-encoding phages was negatively related with bacterial density, leading to dominance of fewer genes at high bacterial abundances. A statistical learning analysis showed that reefs with high microbial density were enriched in viruses encoding genes enabling bacterial recognition and invasion of metazoan epithelium. Over 60% of phages could not have their hosts identified due to limitations of host prediction tools; for those which hosts were identified, host taxonomy was not an indicator of the presence of virulence genes. CONCLUSIONS: This study described bacterial virulence factors encoded in the genomes of bacteriophages at the community level. The results showed that the increase in microbial densities that occurs during coral reef degradation is associated with a change in the genomic repertoire of bacteriophages, specifically in the diversity and distribution of bacterial virulence genes. This suggests that phages are implicated in the rise of pathogens in disturbed marine ecosystems.

Development and Use of Personalized Bacteriophage-Based Therapeutic Cocktails To Treat a Patient with a Disseminated Resistant Acinetobacter baumannii Infection.

Widespread antibiotic use in clinical medicine and the livestock industry has contributed to the global spread of multidrug-resistant (MDR) bacterial pathogens, including Acinetobacter baumannii We report on a method used to produce a personalized bacteriophage-based therapeutic treatment for a 68-year-old diabetic patient with necrotizing pancreatitis complicated by an MDR A. baumannii infection. Despite multiple antibiotic courses and efforts at percutaneous drainage of a pancreatic pseudocyst, the patient deteriorated over a 4-month period. In the absence of effective antibiotics, two laboratories identified nine different bacteriophages with lytic activity for an A. baumannii isolate from the patient. Administration of these bacteriophages intravenously and percutaneously into the abscess cavities was associated with reversal of the patient's downward clinical trajectory, clearance of the A. baumannii infection, and a return to health. The outcome of this case suggests that the methods described here for the production of bacteriophage therapeutics could be applied to similar cases and that more concerted efforts to investigate the use of therapeutic bacteriophages for MDR bacterial infections are warranted.

Metabolomics of reef benthic interactions reveals a bioactive lipid involved in coral defence.

Holobionts are assemblages of microbial symbionts and their macrobial host. As extant representatives of some of the oldest macro-organisms, corals and algae are important for understanding how holobionts develop and interact with one another. Using untargeted metabolomics, we show that non-self interactions altered the coral metabolome more than self-interactions (i.e. different or same genus, respectively). Platelet activating factor (PAF) and Lyso-PAF, central inflammatory modulators in mammals, were major lipid components of the coral holobionts. When corals were damaged during competitive interactions with algae, PAF increased along with expression of the gene encoding Lyso-PAF acetyltransferase; the protein responsible for converting Lyso-PAF to PAF. This shows that self and non-self recognition among some of the oldest extant holobionts involve bioactive lipids identical to those in highly derived taxa like humans. This further strengthens the hypothesis that major players of the immune response evolved during the pre-Cambrian.

Recruitment of Mobile Genetic Elements for Diverse Cellular Functions in Prokaryotes.

Prokaryotic genomes are replete with mobile genetic elements (MGE) that span a continuum of replication autonomy. On numerous occasions during microbial evolution, diverse MGE lose their autonomy altogether but, rather than being quickly purged from the host genome, assume a new function that benefits the host, rendering the immobilized MGE subject to purifying selection, and resulting in its vertical inheritance. This mini-review highlights the diversity of the repurposed (exapted) MGE as well as the plethora of cellular functions that they perform. The principal contribution of the exaptation of MGE and their components is to the prokaryotic functional systems involved in biological conflicts, and in particular, defense against viruses and other MGE. This evolutionary entanglement between MGE and defense systems appears to stem both from mechanistic similarities and from similar evolutionary predicaments whereby both MGEs and defense systems tend to incur fitness costs to the hosts and thereby evolve mechanisms for survival including horizontal mobility, causing host addiction, and exaptation for functions beneficial to the host. The examples discussed demonstrate that the identity of an MGE, overall mobility and relationship with the host cell (mutualistic, symbiotic, commensal, or parasitic) are all factors that affect exaptation.

Fishing for phages in metagenomes: what do we catch, what do we miss?

Metagenomics and metatranscriptomics have become the principal approaches for discovery of novel bacteriophages and preliminary characterization of their ecology and biology. Metagenomic sequencing dramatically expanded the known diversity of tailed and non-tailed phages with double-stranded DNA genomes and those with single-stranded DNA genomes, whereas metatranscriptomics led to the discovery of thousands of new single-stranded RNA phages. Apart from expanding phage diversity, metagenomics studies discover major novel groups of phages with unique features of genome organization, expression strategy and virus-host interaction, such as the putative order 'crAssvirales', which includes the most abundant human-associated viruses. The continued success of metagenomics hinges on the combination of the most powerful computational methods for phage genome assembly and analysis including harnessing CRISPR spacers for the discovery of novel phages and host assignment. Together, these approaches could make a comprehensive characterization of the earth phageome a realistic goal.

Cargo Genes of Tn7-Like Transposons Comprise an Enormous Diversity of Defense Systems, Mobile Genetic Elements, and Antibiotic Resistance Genes.

Transposition is a major mechanism of horizontal gene mobility in prokaryotes. However, exploration of the genes mobilized by transposons (cargo) is hampered by the difficulty in delineating integrated transposons from their surrounding genetic context. Here, we present a computational approach that allowed us to identify the boundaries of 6,549 Tn7-like transposons. We found that 96% of these transposons carry at least one cargo gene. Delineation of distinct communities in a gene-sharing network demonstrates how transposons function as a conduit of genes between phylogenetically distant hosts. Comparative analysis of the cargo genes reveals significant enrichment of mobile genetic elements (MGEs) nested within Tn7-like transposons, such as insertion sequences and toxin-antitoxin modules, and of genes involved in recombination, anti-MGE defense, and antibiotic resistance. More unexpectedly, cargo also includes genes encoding central carbon metabolism enzymes. Twenty-two Tn7-like transposons carry both an anti-MGE defense system and antibiotic resistance genes, illustrating how bacteria can overcome these combined pressures upon acquisition of a single transposon. This work substantially expands the distribution of Tn7-like transposons, defines their evolutionary relationships, and provides a large-scale functional classification of prokaryotic genes mobilized by transposition. IMPORTANCE Transposons are major vehicles of horizontal gene transfer that, in addition to genes directly involved in transposition, carry cargo genes. However, characterization of these genes is hampered by the difficulty of identification of transposon boundaries. We developed a computational approach for detecting transposon ends and applied it to perform a comprehensive census of the cargo genes of Tn7-like transposons, a large class of bacterial mobile genetic elements (MGE), many of which employ a unique, CRISPR-mediated mechanism of site-specific transposition. The cargo genes encompass a striking diversity of MGE, defense, and antibiotic resistance systems. Unexpectedly, we also identified cargo genes encoding metabolic enzymes. Thus, Tn7-like transposons mobilize a vast repertoire of genes that can have multiple effects on the host bacteria.

Analysis of metagenome-assembled viral genomes from the human gut reveals diverse putative CrAss-like phages with unique genomic features.

CrAssphage is the most abundant human-associated virus and the founding member of a large group of bacteriophages, discovered in animal-associated and environmental metagenomes, that infect bacteria of the phylum Bacteroidetes. We analyze 4907 Circular Metagenome Assembled Genomes (cMAGs) of putative viruses from human gut microbiomes and identify nearly 600 genomes of crAss-like phages that account for nearly 87% of the DNA reads mapped to these cMAGs. Phylogenetic analysis of conserved genes demonstrates the monophyly of crAss-like phages, a putative virus order, and of 5 branches, potential families within that order, two of which have not been identified previously. The phage genomes in one of these families are almost twofold larger than the crAssphage genome (145-192 kilobases), with high density of self-splicing introns and inteins. Many crAss-like phages encode suppressor tRNAs that enable read-through of UGA or UAG stop-codons, mostly, in late phage genes. A distinct feature of the crAss-like phages is the recurrent switch of the phage DNA polymerase type between A and B families. Thus, comparative genomic analysis of the expanded assemblage of crAss-like phages reveals aspects of genome architecture and expression as well as phage biology that were not apparent from the previous work on phage genomics.

Thousands of previously unknown phages discovered in whole-community human gut metagenomes.

BACKGROUND: Double-stranded DNA bacteriophages (dsDNA phages) play pivotal roles in structuring human gut microbiomes; yet, the gut virome is far from being fully characterized, and additional groups of phages, including highly abundant ones, continue to be discovered by metagenome mining. A multilevel framework for taxonomic classification of viruses was recently adopted, facilitating the classification of phages into evolutionary informative taxonomic units based on hallmark genes. Together with advanced approaches for sequence assembly and powerful methods of sequence analysis, this revised framework offers the opportunity to discover and classify unknown phage taxa in the human gut. RESULTS: A search of human gut metagenomes for circular contigs encoding phage hallmark genes resulted in the identification of 3738 apparently complete phage genomes that represent 451 putative genera. Several of these phage genera are only distantly related to previously identified phages and are likely to found new families. Two of the candidate families, "Flandersviridae" and "Quimbyviridae", include some of the most common and abundant members of the human gut virome that infect Bacteroides, Parabacteroides, and Prevotella. The third proposed family, "Gratiaviridae," consists of less abundant phages that are distantly related to the families Autographiviridae, Drexlerviridae, and Chaseviridae. Analysis of CRISPR spacers indicates that phages of all three putative families infect bacteria of the phylum Bacteroidetes. Comparative genomic analysis of the three candidate phage families revealed features without precedent in phage genomes. Some "Quimbyviridae" phages possess Diversity-Generating Retroelements (DGRs) that generate hypervariable target genes nested within defense-related genes, whereas the previously known targets of phage-encoded DGRs are structural genes. Several "Flandersviridae" phages encode enzymes of the isoprenoid pathway, a lipid biosynthesis pathway that so far has not been known to be manipulated by phages. The "Gratiaviridae" phages encode a HipA-family protein kinase and glycosyltransferase, suggesting these phages modify the host cell wall, preventing superinfection by other phages. Hundreds of phages in these three and other families are shown to encode catalases and iron-sequestering enzymes that can be predicted to enhance cellular tolerance to reactive oxygen species. CONCLUSIONS: Analysis of phage genomes identified in whole-community human gut metagenomes resulted in the delineation of at least three new candidate families of Caudovirales and revealed diverse putative mechanisms underlying phage-host interactions in the human gut. Addition of these phylogenetically classified, diverse, and distinct phages to public databases will facilitate taxonomic decomposition and functional characterization of human gut viromes. Video abstract.

The Missing Tailed Phages: Prediction of Small Capsid Candidates.

Tailed phages are the most abundant and diverse group of viruses on the planet. Yet, the smallest tailed phages display relatively complex capsids and large genomes compared to other viruses. The lack of tailed phages forming the common icosahedral capsid architectures T = 1 and T = 3 is puzzling. Here, we extracted geometrical features from high-resolution tailed phage capsid reconstructions and built a statistical model based on physical principles to predict the capsid diameter and genome length of the missing small-tailed phage capsids. We applied the model to 3348 isolated tailed phage genomes and 1496 gut metagenome-assembled tailed phage genomes. Four isolated tailed phages were predicted to form T = 3 icosahedral capsids, and twenty-one metagenome-assembled tailed phages were predicted to form T < 3 capsids. The smallest capsid predicted was a T = 4/3 ≈ 1.33 architecture. No tailed phages were predicted to form the smallest icosahedral architecture, T = 1. We discuss the feasibility of the missing T = 1 tailed phage capsids and the implications of isolating and characterizing small-tailed phages for viral evolution and phage therapy.

Gp4 is a nuclease required for morphogenesis of T4-like bacteriophages.

An essential step in the morphogenesis of tailed bacteriophages is the joining of heads and tails to form infectious virions. Our understanding of the maturation of complete virus particles remains incomplete. Through an unknown mechanism, phage T4 gene product 4 (gp4) plays an essential role in the head-tail joining step of T4-like phages. Alignment of T4 gp4 homologs identified a type II restriction endonuclease motif. Purified gp4 from both T4 and a marine T4-like bacteriophage, YC, have non-specific nuclease activity in vitro. Mutation of a single conserved amino acid residue in the endonuclease fold of T4 and YC gp4 abrogates nuclease activity. When expressed in trans, the wild type T4 gp4, but neither the mutated T4 protein nor the YC homolog, rescues a T4 gene 4 amber mutant phage. Thus the nuclease activity appears essential for morphogenesis, potentially by cleaving packaged DNA to enable the joining of heads to tails.

A Distinct Contractile Injection System Gene Cluster Found in a Majority of Healthy Adult Human Microbiomes.

Many commensal bacteria antagonize each other or their host by producing syringe-like secretion systems called contractile injection systems (CIS). Members of the Bacteroidales family have been shown to produce only one type of CIS-a contact-dependent type 6 secretion system that mediates bacterium-bacterium interactions. Here, we show that a second distinct cluster of genes from Bacteroidales bacteria from the human microbiome may encode yet-uncharacterized injection systems that we term Bacteroidales injection systems (BIS). We found that BIS genes are present in the gut microbiomes of 99% of individuals from the United States and Europe and that BIS genes are more prevalent in the gut microbiomes of healthy individuals than in those individuals suffering from inflammatory bowel disease. Gene clusters similar to that of the BIS mediate interactions between bacteria and diverse eukaryotes, like amoeba, insects, and tubeworms. Our findings highlight the ubiquity of the BIS gene cluster in the human gut and emphasize the relevance of the gut microbiome to the human host. These results warrant investigations into the structure and function of the BIS and how they might mediate interactions between Bacteroidales bacteria and the human host or microbiome.IMPORTANCE To engage with host cells, diverse pathogenic bacteria produce syringe-like structures called contractile injection systems (CIS). CIS are evolutionarily related to the contractile tails of bacteriophages and are specialized to puncture membranes, often delivering effectors to target cells. Although CIS are key for pathogens to cause disease, paradoxically, similar injection systems have been identified within healthy human microbiome bacteria. Here, we show that gene clusters encoding a predicted CIS, which we term Bacteroidales injection systems (BIS), are present in the microbiomes of nearly all adult humans tested from Western countries. BIS genes are enriched within human gut microbiomes and are expressed both in vitro and in vivo Further, a greater abundance of BIS genes is present within healthy gut microbiomes than in those humans with with inflammatory bowel disease (IBD). Our discovery provides a potentially distinct means by which our microbiome interacts with the human host or its microbiome.

Cystic Fibrosis Rapid Response: Translating Multi-omics Data into Clinically Relevant Information.

Pulmonary exacerbations are the leading cause of death in cystic fibrosis (CF) patients. To track microbial dynamics during acute exacerbations, a CF rapid response (CFRR) strategy was developed. The CFRR relies on viromics, metagenomics, metatranscriptomics, and metabolomics data to rapidly monitor active members of the viral and microbial community during acute CF exacerbations. To highlight CFRR, a case study of a CF patient is presented, in which an abrupt decline in lung function characterized a fatal exacerbation. The microbial community in the patient's lungs was closely monitored through the multi-omics strategy, which led to the identification of pathogenic shigatoxigenic Escherichia coli (STEC) expressing Shiga toxin. This case study illustrates the potential for the CFRR to deconstruct complicated disease dynamics and provide clinicians with alternative treatments to improve the outcomes of pulmonary exacerbations and expand the life spans of individuals with CF.IMPORTANCE Proper management of polymicrobial infections in patients with cystic fibrosis (CF) has extended their life span. Information about the composition and dynamics of each patient's microbial community aids in the selection of appropriate treatment of pulmonary exacerbations. We propose the cystic fibrosis rapid response (CFRR) as a fast approach to determine viral and microbial community composition and activity during CF pulmonary exacerbations. The CFRR potential is illustrated with a case study in which a cystic fibrosis fatal exacerbation was characterized by the presence of shigatoxigenic Escherichia coli The incorporation of the CFRR within the CF clinic could increase the life span and quality of life of CF patients.

A diversity-generating retroelement encoded by a globally ubiquitous Bacteroides phage.

BACKGROUND: Diversity-generating retroelements (DGRs) are genetic cassettes that selectively mutate target genes to produce hypervariable proteins. First characterized in Bordetella bacteriophage BPP-1, the DGR creates a hypervariable phage tail fiber that enables host tropism switching. Subsequent surveys for DGRs conclude that the majority identified to date are bacterial or archaeal in origin. This work examines bacteriophage and bacterial genomes for novel phage-encoded DGRs. RESULTS: This survey discovered 92 DGRs that were only found in phages exhibiting a temperate lifestyle. The majority of phage-encoded DGRs were identified as prophages in bacterial hosts from the phyla Bacteroidetes, Proteobacteria, and Firmicutes. Sequence reads from these previously unidentified prophages were present in viral metagenomes (viromes), indicating these prophages can produce functional viruses. Five phages possessed hypervariable proteins with structural similarity to the tail fiber of BPP-1, whereas the functions of the remaining DGR target proteins were unknown. A novel temperate phage that harbors a DGR cassette targeting a protein of unknown function was induced from Bacteroides dorei. This phage, here named Bacteroides dorei Hankyphage, lysogenizes 13 different Bacteroides species and was present in 34% and 21% of whole-community metagenomes and human-associated viromes, respectively. CONCLUSIONS: Here, the number of known DGR-containing phages is increased from four to 92. All of these phages exhibit a temperate lifestyle, including a cosmopolitan human-associated phage. Targeted hypervariation by temperate phages may be a ubiquitous mechanism underlying phage-bacteria interaction in the human microbiome.