Tropism switching in Bordetella bacteriophage defines a family of diversity-generating retroelements

The role of noncoding RNAs in bacterial immunity

The evolutionary arms race between bacteria and phages has driven the development of diverse anti-phage defense mechanisms. Recent studies have identified noncoding RNAs (ncRNAs) as key players in bacteria-phage conflicts, including CRISPR-Cas, toxin-antitoxin (TA), and reverse transcriptase (RT)-based defenses; however, our understanding of their roles in immunity is still emerging. In this review, we explore the multifaceted roles of ncRNAs in bacterial immunity, offering insights into their contributions to defense and anti-defense mechanisms, their influence on immune regulatory networks, and potential biotechnological applications. Finally, we highlight key outstanding questions in the field to spark future research directions.

Structural Requirements for Reverse Transcription by a Diversity-generating Retroelement

Diversity-generating retroelements (DGRs) create massive protein sequence variation in ecologically diverse microbes. Variation occurs during reverse transcription of a protein-encoding RNA template coupled to misincorporation at adenosines. In the prototypical Bordetella bacteriophage DGR, the template must be surrounded by upstream and downstream RNA segments for cDNA synthesis by the reverse transcriptase bRT and associated protein Avd. The function of the surrounding RNA was unknown. Cryo-EM revealed that this RNA enveloped bRT and lay over barrel-shaped Avd, forming an intimate ribonucleoprotein (RNP). An abundance of essential interactions between RNA structural elements and bRT-Avd precisely positioned an RNA homoduplex for initiation of cDNA synthesis by cis -priming. Our results explain how the surrounding RNA primes cDNA synthesis, promotes processivity, terminates polymerization, and strictly limits mutagenesis to select proteins through mechanisms that are likely conserved in DGRs from distant taxa.

Mechanisms used for cDNA synthesis and site-specific integration of RNA into DNA genomes by a reverse transcriptase-Cas1 fusion protein

Reverse transcriptase-Cas1 (RT-Cas1) fusion proteins found in some CRISPR systems enable spacer acquisition from both RNA and DNA, but the mechanism of RNA spacer acquisition has remained unclear. Here, we found that Marinomonas mediterranea RT-Cas1/Cas2 adds short 3'-DNA (dN) tails to RNA protospacers, enabling their direct integration into CRISPR arrays as 3'-dN-RNAs or 3'-dN-RNA/cDNA duplexes at rates comparable to similarly configured DNAs. Reverse transcription of RNA protospacers is initiated at 3' proximal sites by multiple mechanisms, including recently described de novo initiation, protein priming with any dNTP, and use of short exogenous or synthesized DNA oligomer primers, enabling synthesis of near full-length cDNAs of diverse RNAs without fixed sequence requirements. The integration of 3'-dN-RNAs or single-stranded DNAs (ssDNAs) is favored over duplexes at higher protospacer concentrations, potentially relevant to spacer acquisition from abundant pathogen RNAs or ssDNA fragments generated by phage defense nucleases. Our findings reveal mechanisms for site-specifically integrating RNA into DNA genomes with potential biotechnological applications.

Targeted hypermutation of putative antigen sensors in multicellular bacteria

Diversity-generating retroelements (DGRs) are used by bacteria, archaea, and viruses as a targeted mutagenesis tool. Through error-prone reverse transcription, DGRs introduce random mutations at specific genomic loci, enabling rapid evolution of these targeted genes. However, the function and benefits of DGR-diversified proteins in cellular hosts remain elusive. We find that 82% of DGRs from one of the major monophyletic lineages of DGR reverse transcriptases are encoded by multicellular bacteria, which often have two or more DGR loci in their genomes. Using the multicellular purple sulfur bacterium Thiohalocapsa sp. PB-PSB1 as an example, we characterized nine distinct DGR loci capable of generating 10282 different combinations of target proteins. With environmental metagenomes from individual Thiohalocapsa aggregates, we show that most of PB-PSB1's DGR target genes are diversified across its biogeographic range, with spatial heterogeneity in the diversity of each locus. In Thiohalocapsa PB-PSB1 and other bacteria hosting this lineage of cellular DGRs, the diversified target genes are associated with NACHT-domain anti-phage defenses and putative ternary conflict systems previously shown to be enriched in multicellular bacteria. We propose that these DGR-diversified targets act as antigen sensors that confer a form of adaptive immunity to their multicellular consortia, though this remains to be experimentally tested. These findings could have implications for understanding the evolution of multicellularity, as the NACHT-domain anti-phage systems and ternary systems share both domain homology and conceptual similarities with the innate immune and programmed cell death pathways of plants and metazoans.

Ultraconserved bacteriophage genome sequence identified in 1300-year-old human palaeofaeces

Bacteriophages are widely recognised as rapidly evolving biological entities. However, knowledge about ancient bacteriophages is limited. Here, we analyse DNA sequence datasets previously generated from ancient palaeofaeces and human gut-content samples, and identify an ancient phage genome nearly identical to present-day Mushuvirus mushu, a virus that infects gut commensal bacteria. The DNA damage patterns of the genome are consistent with its ancient origin and, despite 1300 years of evolution, the ancient Mushuvirus genome shares 97.7% nucleotide identity with its modern counterpart, indicating a long-term relationship between the prophage and its host. In addition, we reconstruct and authenticate 297 other phage genomes from the last 5300 years, including those belonging to unknown families. Our findings demonstrate the feasibility of reconstructing ancient phage genome sequences, thus expanding the known virosphere and offering insights into phage-bacteria interactions spanning several millennia.

Infant gut DNA bacteriophage strain persistence during the first 3 years of life

Bacteriophages are key components of gut microbiomes, yet the phage colonization process in the infant gut remains uncertain. Here, we establish a large phage sequence database and use strain-resolved analyses to investigate DNA phage succession in infants throughout the first 3 years of life. Analysis of 819 fecal metagenomes collected from 28 full-term and 24 preterm infants and their mothers revealed that early-life phageome richness increases over time and reaches adult-like complexity by age 3. Approximately 9% of early phage colonizers, which are mostly maternally transmitted and infect Bacteroides, persist for 3 years and are more prevalent in full-term than in preterm infants. Although rare, phages with stop codon reassignment are more likely to persist than non-recoded phages and generally display an increase in in-frame reassigned stop codons over 3 years. Overall, maternal seeding, stop codon reassignment, host CRISPR-Cas locus prevalence, and diverse phage populations contribute to stable viral colonization.

Tentaclins-A Novel Family of Phage Receptor-Binding Proteins That Can Be Hypermutated by DGR Systems

Diversity-generating retroelements (DGRs) are prokaryotic systems providing rapid modification and adaptation of target proteins. In phages, the main targets of DGRs are receptor-binding proteins that are usually parts of tail structures and the variability of such host-recognizing structures enables phage adaptation to changes on the bacterial host surface. Sometimes, more than one target gene containing a hypermutated variable repeat (VR) can be found in phage DGRs. The role of mutagenesis of two functionally different genes is unclear. In this study, several phage genomes that contain DGRs with two target genes were found in the gut virome of healthy volunteers. Bioinformatics analysis of these genes indicated that they encode proteins with different topology; however, both proteins contain the C-type lectin (C-lec) domain with a hypermutated beta-hairpin on its surface. One of the target proteins belongs to a new family of proteins with a specific topology: N-terminal C-lec domain followed by one or more immunoglobulin domains. Proteins from the new family were named tentaclins after TENTACLe + proteIN. The genes encoding such proteins were found in the genomes of prophages and phages from the gut metagenomes. We hypothesized that tentaclins are involved in binding either to bacterial receptors or intestinal/immune cells.

Mechanisms used for cDNA synthesis and site-specific integration of RNA into DNA genomes by a reverse transcriptase-Cas1 fusion protein

Reverse transcriptase-Cas1 (RT-Cas1) fusion proteins found in some CRISPR systems enable spacer acquisition from both RNA and DNA, but the mechanism of RNA spacer acquisition has remained unclear. Here, we found Marinomonas mediterranea RT-Cas1/Cas2 adds short 3'-DNA (dN) tails to RNA protospacers enabling their direct integration into CRISPR arrays as 3'-dN-RNA/cDNA duplexes or 3'-dN-RNAs at rates comparable to similarly configured DNAs. Reverse transcription of RNA protospacers occurs by multiple mechanisms, including recently described de novo initiation, protein priming with any dNTP, and use of short exogenous or synthesized DNA oligomer primers, enabling synthesis of cDNAs from diverse RNAs without fixed sequence requirements. The integration of 3'-dN-RNAs or single-stranded (ss) DNAs is favored over duplexes at higher protospacer concentrations, potentially relevant to spacer acquisition from abundant pathogen RNAs or ssDNA fragments generated by phage-defense nucleases. Our findings reveal novel mechanisms for site-specifically integrating RNA into DNA genomes with potential biotechnological applications.

Diversity-Generating Retroelements in Prokaryotic Immunity

Adaptive immunity systems found in different organisms fall into two major types. Prokaryotes possess CRISPR-Cas systems that recognize former invaders using memorized (captured) pieces of their DNA as pathogen signatures. Mammals possess a vast repertoire of antibodies and T-cell receptor variants generated in advance. In this second type of adaptive immunity, a pathogen presentation to the immune system specifically activates the cells that express matching antibodies or receptors. These cells proliferate to fight the infection and form the immune memory. The principle of preemptive production of diverse defense proteins for future use can hypothetically take place in microbes too. We propose a hypothesis that prokaryotes employ diversity-generating retroelements to prepare defense proteins against yet-unknown invaders. In this study, we test this hypothesis with the methods of bioinformatics and identify several candidate defense systems based on diversity-generating retroelements.

Community structure and function of microbiomes in polluted stretches of river Yamuna in New Delhi, India, using shotgun metagenomics

The large population residing in the northern region of India surrounding Delhi mostly depends on water of River Yamuna, a tributary of mighty Ganga for agriculture, drinking and various religious activities. However, continuous anthropogenic activities mostly due to pollution mediated by rapid urbanization and industrialization have profoundly affected river microflora and their function thus its health. In this study, potential of whole-genome metagenomics was exploited to unravel the novel consortia of microbiome and their functional potential in the polluted sediments of the river at Delhi. Analysis of high-quality metagenome data from Illumina NextSeq500 revealed substantial differences in composition of microbiota at different sites dominated by Proteobacteria, Bacteroidetes, Firmicutes, Actinobacteria and Chloroflexi phyla. The presence of highly dominant anaerobic bacteria like Dechloromonas aromatica (benzene reducing and denitrifying), Rhodopseudomonas palustris (organic matter reducing), Syntrophus aciditrophicus (fatty acid reducing) and Syntrophobacter fumaroxidans (sulphate reducing) in the polluted river Yamuna signifies the impact of unchecked pollution in declining health of the river ecosystem. A decline in abundance of phages was also noticed along the downstream river Yamuna. Mining of mycobiome reads uncovered plethora of fungal communities (i.e. Nakaseomyces, Aspergillus, Schizosaccharomyces and Lodderomyces) in the polluted stretches due to the availability of higher organic carbon and total nitrogen (%) could be decoded as promising bioindicators of river trophic status. Pathway analysis through KEGG revealed higher abundance of genes involved in energy metabolism (nitrogen and sulphur), methane metabolism, degradation of xenobiotics (Nitrotoluene, Benzoate and Atrazine), two-component system (atoB, cusA and silA) and membrane transport (ABC transporters). Catalase-peroxidase and 4-hydroxybenzoate 3-monooxygenase were the most enriched pollution degrading enzymes in the polluted study sites of river Yamuna. Overall, our results provide crucial insights into microbial dynamics and their function in response to high pollution and could be insightful to the ongoing remediation strategies to clean river Yamuna.

Accelerated Evolution by Diversity-Generating Retroelements

Diversity-generating retroelements (DGRs) create vast amounts of targeted, functional diversity by facilitating the rapid evolution of ligand-binding protein domains. Thousands of DGRs have been identified in bacteria, archaea, and their respective viruses. They are broadly distributed throughout the microbial world, with enrichment observed in certain taxa and environments. The diversification machinery works through a novel mechanism termed mutagenic retrohoming, whereby nucleotide sequence information is copied from an invariant DNA template repeat (TR) into an RNA intermediate, selectively mutagenized at TR adenines during cDNA synthesis by a DGR-encoded reverse transcriptase, and transferred to a variable repeat (VR) region within a variable-protein gene (54). This unidirectional flow of information leaves TR-DNA sequences unmodified, allowing for repeated rounds of mutagenic retrohoming to optimize variable-protein function. DGR target genes are often modular and can encode one or more of a wide variety of discrete functional domains appended to a diversifiable ligand-binding motif. Bacterial variable proteins often localize to cell surfaces, although a subset appear to be cytoplasmic, while phage-encoded DGRs commonly diversify tail fiber-associated receptor-binding proteins. Here, we provide a comprehensive review of the mechanism and consequences of accelerated protein evolution by these unique and beneficial genetic elements. Expected final online publication date for the Annual Review of Microbiology, Volume 76 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

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.

Eco-evolutionary significance of domesticated retroelements in microbial genomes

Since the first discovery of reverse transcriptase in bacteria, and later in archaea, bacterial and archaeal retroelements have been defined by their common enzyme that coordinates diverse functions. Yet, evolutionary refinement has produced distinct retroelements across the tree of microbial life that are perhaps best described in terms of their programmed RNA-a compact sequence that preserves core information for a sophisticated mechanism. From this perspective, reverse transcriptase has been selected as the modular tool for carrying out nature's instructions in various RNA templates. Beneficial retroelements-those that can provide a fitness advantage to their host-evolved to their extant forms in a wide array of microorganisms and their viruses, spanning nearly all habitats. Within each specialized retroelement class, several universal features seem to be shared across diverse taxa, while specific functional and mechanistic insights are based on only a few model retroelement systems from clinical isolates. Currently, little is known about the diversity of cellular functions and ecological significance of retroelements across different biomes. With increasing availability of isolate, metagenome-assembled, and single-amplified genomes, the taxonomic and functional breadth of prokaryotic retroelements is coming into clearer view. This review explores the recently characterized classes of beneficial, yet accessory retroelements of bacteria and archaea. We describe how these specialized mechanisms exploit a form of fixed mobility, whereby the retroelements do not appear to proliferate selfishly throughout the genome. Moreover, we discuss computational approaches for systematic identification of retroelements from vast sequence repositories and highlight recent discoveries in terms of their apparent distribution and ecological significance in nature. Lastly, we present a new perspective on the eco-evolutionary significance of these genetic elements in marine bacteria and demonstrate approaches that enable the characterization of their environmental diversity through metagenomics.

Identification and classification of reverse transcriptases in bacterial genomes and metagenomes

Reverse transcriptases (RTs) are found in different systems including group II introns, Diversity Generating Retroelements (DGRs), retrons, CRISPR-Cas systems, and Abortive Infection (Abi) systems in prokaryotes. Different classes of RTs can play different roles, such as template switching and mobility in group II introns, spacer acquisition in CRISPR-Cas systems, mutagenic retrohoming in DGRs, programmed cell suicide in Abi systems, and recently discovered phage defense in retrons. While some classes of RTs have been studied extensively, others remain to be characterized. There is a lack of computational tools for identifying and characterizing various classes of RTs. In this study, we built a tool (called myRT) for identification and classification of prokaryotic RTs. In addition, our tool provides information about the genomic neighborhood of each RT, providing potential functional clues. We applied our tool to predict RTs in all complete and draft bacterial genomes, and created a collection that can be used for exploration of putative RTs and their associated protein domains. Application of myRT to metagenomes showed that gut metagenomes encode proportionally more RTs related to DGRs, outnumbering retron-related RTs, as compared to the collection of reference genomes. MyRT is both available as a standalone software (https://github.com/mgtools/myRT) and also through a website (https://omics.informatics.indiana.edu/myRT/).

Prokaryotic reverse transcriptases: from retroelements to specialized defense systems

Reverse transcriptases (RTs) catalyze the polymerization of DNA from an RNA template. These enzymes were first discovered in RNA tumor viruses in 1970, but it was not until 1989 that they were found in prokaryotes as a key component of retrons. Apart from RTs encoded by the 'selfish' mobile retroelements known as group II introns, prokaryotic RTs are extraordinarily diverse, but their function has remained elusive. However, recent studies have revealed that different lineages of prokaryotic RTs, including retrons, those associated with CRISPR-Cas systems, Abi-like RTs and other yet uncharacterized RTs, are key components of different lines of defense against phages and other mobile genetic elements. Prokaryotic RTs participate in various antiviral strategies, including abortive infection (Abi), in which the infected cell is induced to commit suicide to protect the host population, adaptive immunity, in which a memory of previous infection is used to build an efficient defense, and other as yet unidentified mechanisms. These prokaryotic enzymes are attracting considerable attention, both for use in cutting-edge technologies, such as genome editing, and as an emerging research topic. In this review, we discuss what is known about prokaryotic RTs, and the exciting evidence for their domestication from retroelements to create specialized defense systems.

Interactions between bacterial and phage communities in natural environments

We commonly acknowledge that bacterial viruses (phages) shape the composition and evolution of bacterial communities in nature and therefore have important roles in ecosystem functioning. This view stems from studies in the 1990s to the first decade of the twenty-first century that revealed high viral abundance, high viral diversity and virus-induced microbial death in aquatic ecosystems as well as an association between collapses in bacterial density and peaks in phage abundance. The recent surge in metagenomic analyses has provided deeper insight into the abundance, genomic diversity and spatio-temporal dynamics of phages in a wide variety of ecosystems, ranging from deep oceans to soil and the mammalian digestive tract. However, the causes and consequences of variations in phage community compositions remain poorly understood. In this Review, we explore current knowledge of the composition and evolution of phage communities, as well as their roles in controlling the population and evolutionary dynamics of bacterial communities. We discuss the need for greater ecological realism in laboratory studies to capture the complexity of microbial communities that thrive in natural environments.

The Human Gut Phageome: Origins and Roles in the Human Gut Microbiome

The investigation of the microbial populations of the human body, known as the microbiome, has led to a revolutionary field of science, and understanding of its impacts on human development and health. The majority of microbiome research to date has focussed on bacteria and other kingdoms of life, such as fungi. Trailing behind these is the interrogation of the gut viruses, specifically the phageome. Bacteriophages, viruses that infect bacterial hosts, are known to dictate the dynamics and diversity of bacterial populations in a number of ecosystems. However, the phageome of the human gut, while of apparent importance, remains an area of many unknowns. In this paper we discuss the role of bacteriophages within the human gut microbiome. We examine the methods used to study bacteriophage populations, how this evolved over time and what we now understand about the phageome. We review the phageome development in infancy, and factors that may influence phage populations in adult life. The role and action of the phageome is then discussed at both a biological-level, and in the broader context of human health and disease.

Bacterial death and TRADD-N domains help define novel apoptosis and immunity mechanisms shared by prokaryotes and metazoans

Several homologous domains are shared by eukaryotic immunity and programmed cell-death systems and poorly understood bacterial proteins. Recent studies show these to be components of a network of highly regulated systems connecting apoptotic processes to counter-invader immunity, in prokaryotes with a multicellular habit. However, the provenance of key adaptor domains, namely those of the Death-like and TRADD-N superfamilies, a quintessential feature of metazoan apoptotic systems, remained murky. Here, we use sensitive sequence analysis and comparative genomics methods to identify unambiguous bacterial homologs of the Death-like and TRADD-N superfamilies. We show the former to have arisen as part of a radiation of effector-associated α-helical adaptor domains that likely mediate homotypic interactions bringing together diverse effector and signaling domains in predicted bacterial apoptosis- and counter-invader systems. Similarly, we show that the TRADD-N domain defines a key, widespread signaling bridge that links effector deployment to invader-sensing in multicellular bacterial and metazoan counter-invader systems. TRADD-N domains are expanded in aggregating marine invertebrates and point to distinctive diversifying immune strategies probably directed both at RNA and retroviruses and cellular pathogens that might infect such communities. These TRADD-N and Death-like domains helped identify several new bacterial and metazoan counter-invader systems featuring underappreciated, common functional principles: the use of intracellular invader-sensing lectin-like (NPCBM and FGS), transcription elongation GreA/B-C, glycosyltransferase-4 family, inactive NTPase (serving as nucleic acid receptors), and invader-sensing GTPase switch domains. Finally, these findings point to the possibility of multicellular bacteria-stem metazoan symbiosis in the emergence of the immune/apoptotic systems of the latter.

Ecology and molecular targets of hypermutation in the global microbiome

Changes in the sequence of an organism's genome, i.e., mutations, are the raw material of evolution. The frequency and location of mutations can be constrained by specific molecular mechanisms, such as diversity-generating retroelements (DGRs). DGRs have been characterized from cultivated bacteria and bacteriophages, and perform error-prone reverse transcription leading to mutations being introduced in specific target genes. DGR loci were also identified in several metagenomes, but the ecological roles and evolutionary drivers of these DGRs remain poorly understood. Here, we analyze a dataset of >30,000 DGRs from public metagenomes, establish six major lineages of DGRs including three primarily encoded by phages and seemingly used to diversify host attachment proteins, and demonstrate that DGRs are broadly active and responsible for >10% of all amino acid changes in some organisms. Overall, these results highlight the constraints under which DGRs evolve, and elucidate several distinct roles these elements play in natural communities.

Novel Siphoviridae Bacteriophages Infecting Bacteroides uniformis Contain Diversity Generating Retroelement

Intestinal phages are abundant and important components of gut microbiota, yet the isolated and characterized representatives that infect abundant gut bacteria are sparse. Here we describe the isolation of human intestinal phages infecting Bacteroidesuniformis. Bacteroides is one of the most common bacterial groups in the global human gut microbiota; however, to date not many Bacteroides specific phages are known. Phages isolated in this study belong to a novel viral genus, Bacuni, within the Siphoviridae family. Their genomes encode diversity-generating retroelements (DGR), which were shown in other bacteriophages to promote phage adaptation to rapidly changing environmental conditions and to broaden their host range. Three isolated phages showed 99.83% genome identity but one of them infected a distinct B. uniformis strain. The tropism of Bacuni phages appeared to be dependent on the interplay of DGR mediated sequence variations of gene encoding putative phage fimbrial tip proteins and mutations in host genes coding for outer-membrane proteins. We found prophages with up to 85% amino acid similarity over two-thirds of the Bacuni phage genome in the B. acidifaciens and Prevotella sp. genomes. Despite the abundance of Bacteroides within the human microbiome, we found Bacuni phages only in a limited subset of published gut metagenomes.

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.

Balance between promiscuity and specificity in phage λ host range

As hosts acquire resistance to viruses, viruses must overcome that resistance to re-establish infectivity, or go extinct. Despite the significant hurdles associated with adapting to a resistant host, viruses are evolutionarily successful and maintain stable coevolutionary relationships with their hosts. To investigate the factors underlying how pathogens adapt to their hosts, we performed a deep mutational scan of the region of the λ tail fiber tip protein that mediates contact with the receptor on λ's host, Escherichia coli. Phages harboring amino acid substitutions were subjected to selection for infectivity on wild type E. coli, revealing a highly restrictive fitness landscape, in which most substitutions completely abrogate function. A subset of positions that are tolerant of mutation in this assay, but diverse over evolutionary time, are associated with host range expansion. Imposing selection for phage infectivity on three λ-resistant hosts, each harboring a different missense mutation in the λ receptor, reveals hundreds of adaptive variants in λ. We distinguish λ variants that confer promiscuity, a general ability to overcome host resistance, from those that drive host-specific infectivity. Both processes may be important in driving adaptation to a novel host.

Determinants of adenine-mutagenesis in diversity-generating retroelements

Diversity-generating retroelements (DGRs) vary protein sequences to the greatest extent known in the natural world. These elements are encoded by constituents of the human microbiome and the microbial ‘dark matter’. Variation occurs through adenine-mutagenesis, in which genetic information in RNA is reverse transcribed faithfully to cDNA for all template bases but adenine. We investigated the determinants of adenine-mutagenesis in the prototypical Bordetella bacteriophage DGR through an in vitro system composed of the reverse transcriptase bRT, Avd protein, and a specific RNA. We found that the catalytic efficiency for correct incorporation during reverse transcription by the bRT-Avd complex was strikingly low for all template bases, with the lowest occurring for adenine. Misincorporation across a template adenine was only somewhat lower in efficiency than correct incorporation. We found that the C6, but not the N1 or C2, purine substituent was a key determinant of adenine-mutagenesis. bRT-Avd was insensitive to the C6 amine of adenine but recognized the C6 carbonyl of guanine. We also identified two bRT amino acids predicted to nonspecifically contact incoming dNTPs, R74 and I181, as promoters of adenine-mutagenesis. Our results suggest that the overall low catalytic efficiency of bRT-Avd is intimately tied to its ability to carry out adenine-mutagenesis.

Systematic prediction of genes functionally associated with bacterial retrons and classification of the encoded tripartite systems

Bacterial retrons consist of a reverse transcriptase (RT) and a contiguous non-coding RNA (ncRNA) gene. One third of annotated retrons carry additional open reading frames (ORFs), the contribution and significance of which in retron biology remains to be determined. In this study we developed a computational pipeline for the systematic prediction of genes specifically associated with retron RTs based on a previously reported large dataset representative of the diversity of prokaryotic RTs. We found that retrons generally comprise a tripartite system composed of the ncRNA, the RT and an additional protein or RT-fused domain with diverse enzymatic functions. These retron systems are highly modular, and their components have coevolved to different extents. Based on the additional module, we classified retrons into 13 types, some of which include additional variants. Our findings provide a basis for future studies on the biological function of retrons and for expanding their biotechnological applications.

The State of the Art in Biodefense Related Bacterial Pathogen Detection Using Bacteriophages: How It Started and How It's Going

Accurate pathogen detection and diagnosis is paramount in clinical success of treating patients. There are two general paradigms in pathogen detection: molecular and immuno-based, and phage-based detection is a third emerging paradigm due to its sensitivity and selectivity. Molecular detection methods look for genetic material specific for a given pathogen in a sample usually by polymerase chain reaction (PCR). Immuno-methods look at the pathogen components (antigens) by antibodies raised against that pathogen specific antigens. There are different variations and products based on these two paradigms with advantages and disadvantages. The third paradigm at least for bacterial pathogen detection entails bacteriophages specific for a given bacterium. Sensitivity and specificity are the two key parameters in any pathogen detection system. By their very nature, bacteriophages afford the best sensitivity for bacterial detection. Bacteria and bacteriophages form the predator-prey pair in the evolutionary arms race and has coevolved over time to acquire the exquisite specificity of the pair, in some instances at the strain level. This specificity has been exploited for diagnostic purposes of various pathogens of concern in clinical and other settings. Many recent reviews focus on phage-based detection and sensor technologies. In this review, we focus on a very special group of pathogens that are of concern in biodefense because of their potential misuse in bioterrorism and their extremely virulent nature and as such fall under the Centers for Disease and Prevention (CDC) Category A pathogen list. We describe the currently available phage methods that are based on the usual modalities of detection from culture, to molecular and immuno- and fluorescent methods. We further highlight the gaps and the needs for more modern technologies and sensors drawing from technologies existing for detection and surveillance of other pathogens of clinical relevance.

Bacterial Retrons Function In Anti-Phage Defense

Retrons are bacterial genetic elements comprised of a reverse transcriptase (RT) and a non-coding RNA (ncRNA). The RT uses the ncRNA as template, generating a chimeric RNA/DNA molecule in which the RNA and DNA components are covalently linked. Although retrons were discovered three decades ago, their function remained unknown. We report that retrons function as anti-phage defense systems. The defensive unit is composed of three components: the RT, the ncRNA, and an effector protein. We examined multiple retron systems and show that they confer defense against a broad range of phages via abortive infection. Focusing on retron Ec48, we show evidence that it "guards" RecBCD, a complex with central anti-phage functions in bacteria. Inhibition of RecBCD by phage proteins activates the retron, leading to abortive infection and cell death. Thus, the Ec48 retron forms a second line of defense that is triggered if the first lines of defense have collapsed.

Role of diversity-generating retroelements for regulatory pathway tuning in cyanobacteria

BACKGROUND

Cyanobacteria maintain extensive repertoires of regulatory genes that are vital for adaptation to environmental stress. Some cyanobacterial genomes have been noted to encode diversity-generating retroelements (DGRs), which promote protein hypervariation through localized retrohoming and codon rewriting in target genes. Past research has shown DGRs to mainly diversify proteins involved in cell-cell attachment or viral-host attachment within viral, bacterial, and archaeal lineages. However, these elements may be critical in driving variation for proteins involved in other core cellular processes.

RESULTS

Members of 31 cyanobacterial genera encode at least one DGR, and together, their retroelements form a monophyletic clade of closely-related reverse transcriptases. This class of retroelements diversifies target proteins with unique domain architectures: modular ligand-binding domains often paired with a second domain that is linked to signal response or regulation. Comparative analysis indicates recent intragenomic duplication of DGR targets as paralogs, but also apparent intergenomic exchange of DGR components. The prevalence of DGRs and the paralogs of their targets is disproportionately high among colonial and filamentous strains of cyanobacteria.

CONCLUSION

We find that colonial and filamentous cyanobacteria have recruited DGRs to optimize a ligand-binding module for apparent function in signal response or regulation. These represent a unique class of hypervariable proteins, which might offer cyanobacteria a form of plasticity to adapt to environmental stress. This analysis supports the hypothesis that DGR-driven mutation modulates signaling and regulatory networks in cyanobacteria, suggestive of a new framework for the utility of localized genetic hypervariation.

Widespread transfer of mobile antibiotic resistance genes within individual gut microbiomes revealed through bacterial Hi-C

The gut microbiome harbors a ‘silent reservoir’ of antibiotic resistance (AR) genes that is thought to contribute to the emergence of multidrug-resistant pathogens through horizontal gene transfer (HGT). To counteract the spread of AR, it is paramount to know which organisms harbor mobile AR genes and which organisms engage in HGT. Despite methods that characterize the overall abundance of AR genes in the gut, technological limitations of short-read sequencing have precluded linking bacterial taxa to specific mobile genetic elements (MGEs) encoding AR genes. Here, we apply Hi-C, a high-throughput, culture-independent method, to surveil the bacterial carriage of MGEs. We compare two healthy individuals with seven neutropenic patients undergoing hematopoietic stem cell transplantation, who receive multiple courses of antibiotics, and are acutely vulnerable to the threat of multidrug-resistant infections. We find distinct networks of HGT across individuals, though AR and mobile genes are associated with more diverse taxa within the neutropenic patients than the healthy subjects. Our data further suggest that HGT occurs frequently over a several-week period in both cohorts. Whereas most efforts to understand the spread of AR genes have focused on pathogenic species, our findings shed light on the role of the human gut microbiome in this process.

Linking antibiotic resistance (AR) in the gut microbiome with their bacterial hosts remains challenging. Here, the authors apply bacterial Hi-C to map mobile genetic elements in metagenomes, and illustrate that genes are present in more diverse taxa in neutropenic patients than healthy subjects.

First crAss-Like Phage Genome Encoding the Diversity-Generating Retroelement (DGR)

A new crAss-like genome encoding diversity-generating retroelement (DGR) was found in the fecal virome of a healthy volunteer. The genome of the phage referred to as the crAssphage LMMB, belonged to the candidate genus I of the AlphacrAssvirinae subfamily. The DGR-cassette of the crAssphage LMMB contained all the essential elements: the gene encoding reverse transcriptase (RT), the target gene (TG) encoding the tail-collar fiber protein, and variable and template repeats (VR and TR) with IMH (initiation of mutagenic homing) and IMH* sequences at the 3′-end of the VR and TR, respectively. Architecture of the DGR-cassette was TG-VR(IMH)-TR(IMH*)-RT and an accessory variable determinant (avd) was absent from the cassette. Analysis of 91 genomes and genome fragments from genus I of the AlphacrAssvirinae showed that 15 (16%) of the genomes had DGRs with the same architecture as the crAssphage LMMB, while 66 of the genomes contained incomplete DGR-cassettes or some elements of the DGR.

Retrons and their applications in genome engineering

Precision genome editing technologies have transformed modern biology. These technologies have arisen from the redirection of natural biological machinery, such as bacteriophage lambda proteins for recombineering and CRISPR nucleases for eliciting site-specific double-strand breaks. Less well-known is a widely distributed class of bacterial retroelements, retrons, that employ specialized reverse transcriptases to produce noncoding intracellular DNAs. Retrons' natural function and mechanism of genetic transmission have remained enigmatic. However, recent studies have harnessed their ability to produce DNA in situ for genome editing and evolution. This review describes retron biology and function in both natural and synthetic contexts. We also highlight areas that require further study to advance retron-based precision genome editing platforms.

Rumen Virus Populations: Technological Advances Enhancing Current Understanding

The rumen contains a multi-kingdom, commensal microbiome, including protozoa, bacteria, archaea, fungi and viruses, which enables ruminant herbivores to ferment and utilize plant feedstuffs that would be otherwise indigestible. Within the rumen, virus populations are diverse and highly abundant, often out-numbering the microbial populations that they both predate on and co-exist with. To date the research effort devoted to understanding rumen-associated viral populations has been considerably less than that given to the other microbial populations, yet their contribution to maintaining microbial population balance, intra-ruminal microbial lysis, fiber breakdown, nutrient cycling and genetic transfer may be highly significant. This review follows the technological advances which have contributed to our current understanding of rumen viruses and drawing on knowledge from other environmental and animal-associated microbiomes, describes the known and potential roles and impacts viruses have on rumen function and speculates on the future directions of rumen viral research.

Highly regulated, diversifying NTP-dependent biological conflict systems with implications for the emergence of multicellularity

Social cellular aggregation or multicellular organization pose increased risk of transmission of infections through the system upon infection of a single cell. The generality of the evolutionary responses to this outside of Metazoa remains unclear. We report the discovery of several thematically unified, remarkable biological conflict systems preponderantly present in multicellular prokaryotes. These combine thresholding mechanisms utilizing NTPase chaperones (the MoxR-vWA couple), GTPases and proteolytic cascades with hypervariable effectors, which vary either by using a reverse transcriptase-dependent diversity-generating system or through a system of acquisition of diverse protein modules, typically in inactive form, from various cellular subsystems. Conciliant lines of evidence indicate their deployment against invasive entities, like viruses, to limit their spread in multicellular/social contexts via physical containment, dominant-negative interactions or apoptosis. These findings argue for both a similar operational 'grammar' and shared protein domains in the sensing and limiting of infections during the multiple emergences of multicellularity.

Comparison of CRISPR and Marker-Based Methods for the Engineering of Phage T7

With the recent rise in interest in using lytic bacteriophages as therapeutic agents, there is an urgent requirement to understand their fundamental biology to enable the engineering of their genomes. Current methods of phage engineering rely on homologous recombination, followed by a system of selection to identify recombinant phages. For bacteriophage T7, the host genes cmk or trxA have been used as a selection mechanism along with both type I and II CRISPR systems to select against wild-type phage and enrich for the desired mutant. Here, we systematically compare all three systems; we show that the use of marker-based selection is the most efficient method and we use this to generate multiple T7 tail fibre mutants. Furthermore, we found the type II CRISPR-Cas system is easier to use and generally more efficient than a type I system in the engineering of phage T7. These results provide a foundation for the future, more efficient engineering of bacteriophage T7.

The interaction of phages and bacteria: the co-evolutionary arms race

Since the dawn of life, bacteria and phages are locked in a constant battle and both are perpetually changing their tactics to overcome each other. Bacteria use various strategies to overcome the invading phages, including adsorption inhibition, restriction-modification (R/E) systems, CRISPR-Cas (clustered regularly interspaced short palindromic repeats-CRISPR-associated proteins) systems, abortive infection (Abi), etc. To counteract, phages employ intelligent tactics for the nullification of bacterial defense systems, such as accessing host receptors, evading R/E systems, and anti-CRISPR proteins. Intense knowledge about the details of these defense pathways is the basis for their broad utilities in various fields of research from microbiology to biotechnology. Hence, in this review, we discuss some strategies used by bacteria to inhibit phage infections as well as phage tactics to circumvent bacterial defense systems. In addition, the application of these strategies will be described as a lesson learned from bacteria and phage combats. The ecological factors that affect the evolution of bacterial immune systems is the other issue represented in this review.

The "fighting wisdom and bravery" of tailed phage and host in the process of adsorption

In the process of bacteriophage and bacteria struggle, adsorption is the key factor to determine who is the winner. In this paper, the molecular mechanism of tailed bacteriophage recognition and adsorption to host and the strategy of "fighting wisdom and courage" between them are reviewed.

Template-assisted synthesis of adenine-mutagenized cDNA by a retroelement protein complex

Diversity-generating retroelements (DGRs) create unparalleled levels of protein sequence variation through mutagenic retrohoming. Sequence information is transferred from an invariant template region (TR), through an RNA intermediate, to a protein-coding variable region. Selective infidelity at adenines during transfer is a hallmark of DGRs from disparate bacteria, archaea, and microbial viruses. We recapitulated selective infidelity in vitro for the prototypical Bordetella bacteriophage DGR. A complex of the DGR reverse transcriptase bRT and pentameric accessory variability determinant (Avd) protein along with DGR RNA were necessary and sufficient for synthesis of template-primed, covalently linked RNA–cDNA molecules, as observed in vivo. We identified RNA–cDNA molecules to be branched and most plausibly linked through 2′-5′ phosphodiester bonds. Adenine-mutagenesis was intrinsic to the bRT-Avd complex, which displayed unprecedented promiscuity while reverse transcribing adenines of either DGR or non-DGR RNA templates. In contrast, bRT-Avd processivity was strictly dependent on the template, occurring only for the DGR RNA. This restriction was mainly due to a noncoding segment downstream of TR, which specifically bound Avd and created a privileged site for processive polymerization. Restriction to DGR RNA may protect the host genome from damage. These results define the early steps in a novel pathway for massive sequence diversification.

Discovery and characterization of the evolution, variation and functions of diversity-generating retroelements using thousands of genomes and metagenomes

Background

Diversity-generating retroelements (DGRs) are a unique family of retroelements that generate sequence diversity of DNA to benefit their hosts by introducing variations and accelerating the evolution of target proteins. They exist widely in bacteria, archaea, phage and plasmid. However, our understanding about DGRs in natural environments was still very limited.

Results

We developed an efficient computational algorithm to identify DGRs, and applied it to characterize DGRs in more than 80,000 sequenced bacterial genomes as well as more than 4,000 human metagenome datasets. In total, we identified 948 non-redundant DGRs, which expanded the number of known DGRs in bacterial genomes and human microbiomes by about 55%, and provided a much more comprehensive reference for the study of DGRs. Phylogenetic analysis was done for identified DGRs. The putative target genes of DGRs were searched, and the functions of these target genes were investigated with a comprehensive alignment against the nr database.

Conclusions

DGR system is a powerful and universal mechanism to generate diversity. DGR evolution is closely associated with the living environment and their cassette structures. Furthermore, it may impact a wide range of functional processes in addition to receptor-binding. These results significantly improved our understanding about DGRs.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5951-3) contains supplementary material, which is available to authorized users.

Diversity-generating retroelements: natural variation, classification and evolution inferred from a large-scale genomic survey

Diversity-generating retroelements (DGRs) are novel genetic elements that use reverse transcription to generate vast numbers of sequence variants in specific target genes. Here, we present a detailed comparative bioinformatic analysis that depicts the landscape of DGR sequences in nature as represented by data in GenBank. Over 350 unique DGRs are identified, which together form a curated reference set of putatively functional DGRs. We classify target genes, variable repeats and DGR cassette architectures, and identify two new accessory genes. The great variability of target genes implies roles of DGRs in many undiscovered biological processes. There is much evidence for horizontal transfers of DGRs, and we identify lineages of DGRs that appear to have specialized properties. Because GenBank contains data from only 10% of described species, the compilation may not be wholly representative of DGRs present in nature. Indeed, many DGR subtypes are present only once in the set and DGRs of the candidate phylum radiation bacteria, and Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, Nanohaloarchaea archaea, are exceptionally diverse in sequence, with little information available about functions of their target genes. Nonetheless, this study provides a detailed framework for classifying and studying DGRs as they are uncovered and studied in the future.

MyDGR: a server for identification and characterization of diversity-generating retroelements

MyDGR is a web server providing integrated prediction and visualization of Diversity-Generating Retroelements (DGR) systems in query nucleotide sequences. It is built upon an enhanced version of DGRscan, a tool we previously developed for identification of DGR systems. DGR systems are remarkable genetic elements that use error-prone reverse transcriptases to generate vast sequence variants in specific target genes, which have been shown to benefit their hosts (bacteria, archaea or phages). As the first web server for annotation of DGR systems, myDGR is freely available on the web at http://omics.informatics.indiana.edu/myDGR with all major browsers supported. MyDGR accepts query nucleotide sequences in FASTA format, and outputs all the important features of a predicted DGR system, including a reverse transcriptase, a template repeat and one (or more) variable repeats and their alignment featuring A-to-N (N can be C, T or G) substitutions, and VR-containing target gene(s). In addition to providing the results as text files for download, myDGR generates a visual summary of the results for users to explore the predicted DGR systems. Users can also directly access pre-calculated, putative DGR systems identified in currently available reference bacterial genomes and a few other collections of sequences (including human microbiomes).

Beyond cells - The virome in the human holobiont

Viromics, or viral metagenomics, is a relatively new and burgeoning field of research that studies the complete collection of viruses forming part of the microbiota in any given niche. It has strong foundations rooted in over a century of discoveries in the field of virology and recent advances in molecular biology and sequencing technologies. Historically, most studies have deconstructed the concept of viruses into a simplified perception of viral agents as mere pathogens, which demerits the scope of large-scale viromic analyses. Viruses are, in fact, much more than regular parasites. They are by far the most dynamic and abundant entity and the greatest killers on the planet, as well as the most effective geo-transforming genetic engineers and resource recyclers, acting on all life strata in any habitat. Yet, most of this uncanny viral world remains vastly unexplored to date, greatly hindered by the bewildering complexity inherent to such studies and the methodological and conceptual limitations. Viromic studies are just starting to address some of these issues but they still lag behind microbial metagenomics. In recent years, however, higher-throughput analysis and resequencing have rekindled interest in a field that is just starting to show its true potential. In this review, we take a look at the scientific and technological developments that led to the advent of viral and bacterial metagenomics with a particular, but not exclusive, focus on human viromics from an ecological perspective. We also address some of the most relevant challenges that current viral studies face and ponder on the future directions of the field.

Synthetic evolution

The combination of modern biotechnologies such as DNA synthesis, λ red recombineering, CRISPR-based editing and next-generation high-throughput sequencing increasingly enables precise manipulation of genes and genomes. Beyond rational design, these technologies also enable the targeted, and potentially continuous, introduction of multiple mutations. While this might seem to be merely a return to natural selection, the ability to target evolution greatly reduces fitness burdens and focuses mutation and selection on those genes and traits that best contribute to a desired phenotype, ultimately throwing evolution into fast forward.

Development of CRISPR-Cas systems for genome editing and beyond

The development of clustered regularly interspaced short-palindromic repeat (CRISPR)-Cas systems for genome editing has transformed the way life science research is conducted and holds enormous potential for the treatment of disease as well as for many aspects of biotechnology. Here, I provide a personal perspective on the development of CRISPR-Cas9 for genome editing within the broader context of the field and discuss our work to discover novel Cas effectors and develop them into additional molecular tools. The initial demonstration of Cas9-mediated genome editing launched the development of many other technologies, enabled new lines of biological inquiry, and motivated a deeper examination of natural CRISPR-Cas systems, including the discovery of new types of CRISPR-Cas systems. These new discoveries in turn spurred further technological developments. I review these exciting discoveries and technologies as well as provide an overview of the broad array of applications of these technologies in basic research and in the improvement of human health. It is clear that we are only just beginning to unravel the potential within microbial diversity, and it is quite likely that we will continue to discover other exciting phenomena, some of which it may be possible to repurpose as molecular technologies. The transformation of mysterious natural phenomena to powerful tools, however, takes a collective effort to discover, characterize, and engineer them, and it has been a privilege to join the numerous researchers who have contributed to this transformation of CRISPR-Cas systems.

Interaction of phages, bacteria, and the human immune system: Evolutionary changes in phage therapy

Phages and bacteria are known to undergo dynamic and co-evolutionary arms race interactions in order to survive. Recent advances from in vitro and in vivo studies have improved our understanding of the complex interactions between phages, bacteria, and the human immune system. This insight is essential for the development of phage therapy to battle the growing problems of antibiotic resistance. It is also pivotal to prevent the development of phage-resistance during the implementation of phage therapy in the clinic. In this review, we discuss recent progress of the interactions between phages, bacteria, and the human immune system and its clinical application for phage therapy. Proper phage therapy design will ideally produce large burst sizes, short latent periods, broad host ranges, and a low tendency to select resistance.

CRISPR Spacers Indicate Preferential Matching of Specific Virioplankton Genes

Viral infection exerts selection pressure on marine microbes, as virus-induced cell lysis causes 20 to 50% of cell mortality, resulting in fluxes of biomass into oceanic dissolved organic matter. Archaeal and bacterial populations can defend against viral infection using the clustered regularly interspaced short palindromic repeat (CRISPR)-associated (Cas) system, which relies on specific matching between a spacer sequence and a viral gene. If a CRISPR spacer match to any gene within a viral genome is equally effective in preventing lysis, no viral genes should be preferentially matched by CRISPR spacers. However, if there are differences in effectiveness, certain viral genes may demonstrate a greater frequency of CRISPR spacer matches. Indeed, homology search analyses of bacterioplankton CRISPR spacer sequences against virioplankton sequences revealed preferential matching of replication proteins, nucleic acid binding proteins, and viral structural proteins. Positive selection pressure for effective viral defense is one parsimonious explanation for these observations. CRISPR spacers from virioplankton metagenomes preferentially matched methyltransferase and phage integrase genes within virioplankton sequences. These virioplankton CRISPR spacers may assist infected host cells in defending against competing phage. Analyses also revealed that half of the spacer-matched viral genes were unknown, some genes matched several spacers, and some spacers matched multiple genes, a many-to-many relationship. Thus, CRISPR spacer matching may be an evolutionary algorithm, agnostically identifying those genes under stringent selection pressure for sustaining viral infection and lysis. Investigating this subset of viral genes could reveal those genetic mechanisms essential to virus-host interactions and provide new technologies for optimizing CRISPR defense in beneficial microbes.IMPORTANCE The CRISPR-Cas system is one means by which bacterial and archaeal populations defend against viral infection which causes 20 to 50% of cell mortality in the ocean. We tested the hypothesis that certain viral genes are preferentially targeted for the initial attack of the CRISPR-Cas system on a viral genome. Using CASC, a pipeline for CRISPR spacer discovery, and metagenome data from oceanic microbes and viruses, we found a clear subset of viral genes with high match frequencies to CRISPR spacers. Moreover, we observed a many-to-many relationship of spacers and viral genes. These high-match viral genes were involved in nucleotide metabolism, DNA methylation, and viral structure. It is possible that CRISPR spacer matching is an evolutionary algorithm pointing to those viral genes most important to sustaining infection and lysis. Studying these genes may advance the understanding of virus-host interactions in nature and provide new technologies for leveraging CRISPR-Cas systems in beneficial microbes.

Contemporary Phage Biology: From Classic Models to New Insights

Bacteriophages, discovered about a century ago, have been pivotal as models for understanding the fundamental principles of molecular biology. While interest in phage biology declined after the phage "golden era," key recent developments, including advances in phage genomics, microscopy, and the discovery of the CRISPR-Cas anti-phage defense system, have sparked a renaissance in phage research in the past decade. This review highlights recently discovered unexpected complexities in phage biology, describes a new arsenal of phage genes that help them overcome bacterial defenses, and discusses advances toward documentation of the phage biodiversity on a global scale.

Crystal structure of a Thermus aquaticus diversity-generating retroelement variable protein

Diversity-generating retroelements (DGRs) are widely distributed in bacteria, archaea, and microbial viruses, and bring about unparalleled levels of sequence variation in target proteins. While DGR variable proteins share low sequence identity, the structures of several such proteins have revealed the C-type lectin (CLec)-fold as a conserved scaffold for accommodating massive sequence variation. This conservation has led to the suggestion that the CLec-fold may be useful in molecular surface display applications. Thermostability is an attractive feature in such applications, and thus we studied the variable protein of a DGR encoded by a prophage of the thermophile Thermus aquaticus. We report here the 2.8 Å resolution crystal structure of the variable protein from the T. aquaticus DGR, called TaqVP, and confirm that it has a CLec-fold. Remarkably, its variable region is nearly identical in structure to those of several other CLec-fold DGR variable proteins despite low sequence identity among these. TaqVP was found to be thermostable, which appears to be a property shared by several CLec-fold DGR variable proteins. These results provide impetus for the pursuit of the DGR variable protein CLec-fold in molecular display applications.

Genomic diversity, lifestyles and evolutionary origins of DPANN archaea

Archaea-a primary domain of life besides Bacteria-have for a long time been regarded as peculiar organisms that play marginal roles in biogeochemical cycles. However, this picture changed with the discovery of a large diversity of archaea in non-extreme environments enabled by the use of cultivation-independent methods. These approaches have allowed the reconstruction of genomes of uncultivated microorganisms and revealed that archaea are diverse and broadly distributed in the biosphere and seemingly include a large diversity of putative symbiotic organisms, most of which belong to the tentative archaeal superphylum referred to as DPANN. This archaeal group encompasses at least 10 different lineages and includes organisms with extremely small cell and genome sizes and limited metabolic capabilities. Therefore, many members of DPANN may be obligately dependent on symbiotic interactions with other organisms and may even include novel parasites. In this contribution, we review the current knowledge of the gene repertoires and lifestyles of members of this group and discuss their placement in the tree of life, which is the basis for our understanding of the deep microbial roots and the role of symbiosis in the evolution of life on Earth.

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.