Diversity-generating retroelement homing regenerates target sequences for repeated rounds of codon rewriting and protein diversification

RNA control of reverse transcription in a diversity-generating retroelement

Diversity-generating retroelements (DGRs) create massive protein sequence variation (up to 1030)1 in ecologically diverse microorganisms. A recent survey identified around 31,000 DGRs from more than 1,500 bacterial and archaeal genera, constituting more than 90 environment types2. DGRs are especially enriched in the human gut microbiome2,3 and nano-sized microorganisms that seem to comprise most microbial life and maintain DGRs despite reduced genomes4,5. DGRs are also implicated in the emergence of multicellularity6,7. 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 complementary DNA synthesis to be carried out by a complex of the DGR reverse transcriptase bRT and associated protein Avd. The function of the surrounding RNA was unknown. Here we show through cryogenic electron microscopy that this RNA envelops bRT and lies over the barrel-shaped Avd, forming an intimate ribonucleoprotein. An abundance of essential interactions in the ribonucleoprotein precisely position an RNA homoduplex in the bRT active site for initiation of reverse transcription. Our results explain how the surrounding RNA primes complementary DNA synthesis, promotes processivity, terminates polymerization and strictly limits mutagenesis to specific proteins through mechanisms that are probably conserved in DGRs belonging to distant taxa.

Targeted protein evolution in the gut microbiome by diversity-generating retroelements

Diversity-generating retroelements (DGRs) accelerate evolution by rapidly diversifying variable proteins. The human gastrointestinal microbiota harbors the greatest density of DGRs known in nature, suggesting they play adaptive roles in this environment. We identified >1,100 unique DGRs among human-associated Bacteroides species and discovered a subset that diversify adhesive components of Type V pili and related proteins. We show that Bacteroides DGRs are horizontally transferred across species, that some are highly active while others are tightly controlled, and that they preferentially alter the functional characteristics of ligand-binding residues on adhesive organelles. Specific variable protein sequences are enriched when Bacteroides strains compete with other commensal bacteria in gnotobiotic mice. Analysis of >2,700 DGRs from diverse phyla in mother-infant pairs shows that Bacteroides DGRs are preferentially transferred to vaginally delivered infants where they actively diversify. Our observations provide a foundation for understanding the roles of stochastic, targeted genome plasticity in shaping host-associated microbial communities.

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.

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.

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.

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.

Does over a century of aerobic phage work provide a solid framework for the study of phages in the gut?

Bacterial viruses (bacteriophages, phages) of the gut have increasingly become a focus in microbiome studies, with an understanding that they are likely key players in health and disease. However, characterization of the virome remains largely based on bioinformatic approaches, with the impact of these viromes inferred based on a century of knowledge from aerobic phage work. Studying the phages infecting anaerobes is difficult, as they are often technically demanding to isolate and propagate. In this review, we primarily discuss the phages infecting three well-studied anaerobes in the gut: Bifidobacterium, Clostridia and Bacteroides, with a particular focus on the challenges in isolating and characterizing these phages. We contrast the lessons learned from these to other anaerobic work on phages infecting facultative anaerobes of the gut: Enterococcus and Lactobacillus. Phages from the gut do appear to adhere to the lessons learned from aerobic work, but the additional challenges of working on them has required ingenious new approaches to enable their study. This, in turn, has uncovered remarkable biology likely underpinning phage-host relationships in many stable environments.

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.

DGR mutagenic transposition occurs via hypermutagenic reverse transcription primed by nicked template RNA

Diversity-generating retroelements (DGRs) are molecular evolution machines that facilitate microbial adaptation to environmental changes. Hypervariation occurs via a mutagenic retrotransposition process from a template repeat (TR) to a variable repeat (VR) that results in adenine-to-random nucleotide conversions. Here we show that reverse transcription of the Bordetella phage DGR is primed by an adenine residue in TR RNA and is dependent on the DGR-encoded reverse transcriptase (bRT) and accessory variability determinant (Avd ), but is VR-independent. We also find that the catalytic center of bRT plays an essential role in site-specific cleavage of TR RNA for cDNA priming. Adenine-specific mutagenesis occurs during reverse transcription and does not involve dUTP incorporation, indicating it results from bRT-catalyzed misincorporation of standard deoxyribonucleotides. In vivo assays show that this hybrid RNA-cDNA molecule is required for mutagenic transposition, revealing a unique mechanism of DNA hypervariation for microbial adaptation.

Retroelement-guided protein diversification abounds in vast lineages of Bacteria and Archaea

Major radiations of enigmatic bacteria and archaea with large inventories of uncharacterized proteins are a striking feature of the Tree of Life^1,2,3,4,5. The processes that led to functional diversity in these lineages, which may contribute to a host-dependent lifestyle, are poorly understood. Here we show that diversity-generating retroelements (DGRs), which guide site-specific protein hypervariability^6,7,8, are prominent features of genomically-reduced organisms from the bacterial candidate phyla radiation (CPR) and yet uncultivated phyla belonging to the DPANN archaeal superphylum. From reconstructed genomes we defined monophyletic bacterial and archaeal DGR lineages that expand known DGR range by 120% and reveal a history of horizontal retroelement transfer. Retroelement-guided diversification is further shown to be active in current CPR and DPANN populations, with an assortment of protein targets potentially involved in attachment, defense, and regulation. Based on observations of DGR abundance, function, and evolutionary history, we find that targeted protein diversification is a pronounced trait of CPR and DPANN phyla compared to other bacterial and archaeal phyla. This diversification mechanism may provide CPR and DPANN organisms a versatile tool that could be used for adaptation to a dynamic, host-dependent, existence.

Genomic and Metagenomic Analysis of Diversity-Generating Retroelements Associated with Treponema denticola

Diversity-generating retroelements (DGRs) are genetic cassettes that can produce massive protein sequence variation in prokaryotes. Presumably DGRs confer selective advantages to their hosts (bacteria or viruses) by generating variants of target genes—typically resulting in target proteins with altered ligand-binding specificity—through a specialized error-prone reverse transcription process. The only extensively studied DGR system is from the Bordetella phage BPP-1, although DGRs are predicted to exist in other species. Using bioinformatics analysis, we discovered that the DGR system associated with the Treponema denticola species (a human oral-associated periopathogen) is dynamic (with gains/losses of the system found in the isolates) and diverse (with multiple types found in isolated genomes and the human microbiota). The T. denticola DGR is found in only nine of the 17 sequenced T. denticola strains. Analysis of the DGR-associated template regions and reverse transcriptase gene sequences revealed two types of DGR systems in T. denticola: the ATCC35405-type shared by seven isolates including ATCC35405; and the SP32-type shared by two isolates (SP32 and SP33), suggesting multiple DGR acquisitions. We detected additional variants of the T. denticola DGR systems in the human microbiomes, and found that the SP32-type DGR is more abundant than the ATCC35405-type in the healthy human oral microbiome, although the latter is found in more sequenced isolates. This is the first comprehensive study to characterize the DGRs associated with T. denticola in individual genomes as well as human microbiomes, demonstrating the importance of utilizing both individual genomes and metagenomes for characterizing the elements, and for analyzing their diversity and distribution in human populations.

An Unexplored Diversity of Reverse Transcriptases in Bacteria

Reverse transcriptases (RTs) are usually thought of as eukaryotic enzymes, but they are also present in bacteria and likely originated in bacteria and migrated to eukaryotes. Only three types of bacterial retroelements have been substantially characterized: group II introns, diversity-generating retroelements, and retrons. Recent work, however, has identified a myriad of uncharacterized RTs and RT-related sequences in bacterial genomes, which exhibit great sequence diversity and a range of domain structures. Apart from group II introns, none of these putative RTs show evidence of active retromobility. Instead, available information suggests that they are involved in useful processes, sometimes related to phages or phage resistance. This article reviews our knowledge of both characterized and uncharacterized RTs in bacteria. The range of their sequences and genomic contexts promises the discovery of new biochemical reactions and biological phenomena.

Diversity-generating Retroelements in Phage and Bacterial Genomes

Diversity-generating retroelements (DGRs) are DNA diversification machines found in diverse bacterial and bacteriophage genomes that accelerate the evolution of ligand-receptor interactions. Diversification results from a unidirectional transfer of sequence information from an invariant template repeat (TR) to a variable repeat (VR) located in a protein-encoding gene. Information transfer is coupled to site-specific mutagenesis in a process called mutagenic homing, which occurs through an RNA intermediate and is catalyzed by a unique, DGR-encoded reverse transcriptase that converts adenine residues in the TR into random nucleotides in the VR. In the prototype DGR found in the Bordetella bacteriophage BPP-1, the variable protein Mtd is responsible for phage receptor recognition. VR diversification enables progeny phage to switch tropism, accelerating their adaptation to changes in sequence or availability of host cell-surface molecules for infection. Since their discovery, hundreds of DGRs have been identified, and their functions are just beginning to be understood. VR-encoded residues of many DGR-diversified proteins are displayed in the context of a C-type lectin fold, although other scaffolds, including the immunoglobulin fold, may also be used. DGR homing is postulated to occur through a specialized target DNA-primed reverse transcription mechanism that allows repeated rounds of diversification and selection, and the ability to engineer DGRs to target heterologous genes suggests applications for bioengineering. This chapter provides a comprehensive review of our current understanding of this newly discovered family of beneficial retroelements.

Regulatory RNAs in photosynthetic cyanobacteria

Regulatory RNAs play versatile roles in bacteria in the coordination of gene expression during various physiological processes, especially during stress adaptation. Photosynthetic bacteria use sunlight as their major energy source. Therefore, they are particularly vulnerable to the damaging effects of excess light or UV irradiation. In addition, like all bacteria, photosynthetic bacteria must adapt to limiting nutrient concentrations and abiotic and biotic stress factors. Transcriptome analyses have identified hundreds of potential regulatory small RNAs (sRNAs) in model cyanobacteria such as Synechocystis sp. PCC 6803 or Anabaena sp. PCC 7120, and in environmentally relevant genera such as Trichodesmium, Synechococcus and Prochlorococcus. Some sRNAs have been shown to actually contain μORFs and encode short proteins. Examples include the 40-amino-acid product of the sml0013 gene, which encodes the NdhP subunit of the NDH1 complex. In contrast, the functional characterization of the non-coding sRNA PsrR1 revealed that the 131 nt long sRNA controls photosynthetic functions by targeting multiple mRNAs, providing a paradigm for sRNA functions in photosynthetic bacteria. We suggest that actuatons comprise a new class of genetic elements in which an sRNA gene is inserted upstream of a coding region to modify or enable transcription of that region.

Targeted diversity generation by intraterrestrial archaea and archaeal viruses

In the evolutionary arms race between microbes, their parasites, and their neighbours, the capacity for rapid protein diversification is a potent weapon. Diversity-generating retroelements (DGRs) use mutagenic reverse transcription and retrohoming to generate myriad variants of a target gene. Originally discovered in pathogens, these retroelements have been identified in bacteria and their viruses, but never in archaea. Here we report the discovery of intact DGRs in two distinct intraterrestrial archaeal systems: a novel virus that appears to infect archaea in the marine subsurface, and, separately, two uncultivated nanoarchaea from the terrestrial subsurface. The viral DGR system targets putative tail fibre ligand-binding domains, potentially generating >10¹⁸ protein variants. The two single-cell nanoarchaeal genomes each possess ≥4 distinct DGRs. Against an expected background of low genome-wide mutation rates, these results demonstrate a previously unsuspected potential for rapid, targeted sequence diversification in intraterrestrial archaea and their viruses.

Diversity-generating retroelements (DGRs) are genetic elements that introduce sequence variation within target genes in bacteria and their viruses. Here, Paul et al. report the discovery of DGRs in an archaeal virus and in two archaea from marine and terrestrial subsurface environments, respectively.

The low incidence of diversity-generating retroelements in sequenced genomes

The insertion of a retrotransposable element is usually associated with adverse or, at best, neutral effects on the host. Diversity-generating retroelements (DGRs) are the first elements that seem to offer a direct selective advantage to their phage or prokaryote host by exact replacement of a short, defined region of a host gene with a hypermutated variant. In a previous study, we presented the software DiGReF for identification of DGRs in genome sequences, and compiled the first comprehensive set of diversity-generating retroelements in public databases. We identified 155 elements in more than 6000 prokaryotic and phage genomes, which was a surprisingly low number. In this commentary, we will discuss the low incidence of these elements and speculate about the biological role of bacterial DGRs.

The primary transcriptome of the marine diazotroph Trichodesmium erythraeum IMS101

Blooms of the dinitrogen-fixing marine cyanobacterium Trichodesmium considerably contribute to new nitrogen inputs into tropical oceans. Intriguingly, only 60% of the Trichodesmium erythraeum IMS101 genome sequence codes for protein, compared with ~85% in other sequenced cyanobacterial genomes. The extensive non-coding genome fraction suggests space for an unusually high number of unidentified, potentially regulatory non-protein-coding RNAs (ncRNAs). To identify the transcribed fraction of the genome, here we present a genome-wide map of transcriptional start sites (TSS) at single nucleotide resolution, revealing the activity of 6,080 promoters. We demonstrate that T. erythraeum has the highest number of actively splicing group II introns and the highest percentage of TSS yielding ncRNAs of any bacterium examined to date. We identified a highly transcribed retroelement that serves as template repeat for the targeted mutation of at least 12 different genes by mutagenic homing. Our findings explain the non-coding portion of the T. erythraeum genome by the transcription of an unusually high number of non-coding transcripts in addition to the known high incidence of transposable elements. We conclude that riboregulation and RNA maturation-dependent processes constitute a major part of the Trichodesmium regulatory apparatus.

Bacterial genome instability

Bacterial genomes are remarkably stable from one generation to the next but are plastic on an evolutionary time scale, substantially shaped by horizontal gene transfer, genome rearrangement, and the activities of mobile DNA elements. This implies the existence of a delicate balance between the maintenance of genome stability and the tolerance of genome instability. In this review, we describe the specialized genetic elements and the endogenous processes that contribute to genome instability. We then discuss the consequences of genome instability at the physiological level, where cells have harnessed instability to mediate phase and antigenic variation, and at the evolutionary level, where horizontal gene transfer has played an important role. Indeed, this ability to share DNA sequences has played a major part in the evolution of life on Earth. The evolutionary plasticity of bacterial genomes, coupled with the vast numbers of bacteria on the planet, substantially limits our ability to control disease.

A new topology of the HK97-like fold revealed in Bordetella bacteriophage by cryoEM at 3.5 A resolution

Bacteriophage BPP-1 infects and kills Bordetella species that cause whooping cough. Its diversity-generating retroelement (DGR) provides a naturally occurring phage-display system, but engineering efforts are hampered without atomic structures. Here, we report a cryo electron microscopy structure of the BPP-1 head at 3.5 Å resolution. Our atomic model shows two of the three protein folds representing major viral lineages: jellyroll for its cement protein (CP) and HK97-like ('Johnson') for its major capsid protein (MCP). Strikingly, the fold topology of MCP is permuted non-circularly from the Johnson fold topology previously seen in viral and cellular proteins. We illustrate that the new topology is likely the only feasible alternative of the old topology. β-sheet augmentation and electrostatic interactions contribute to the formation of non-covalent chainmail in BPP-1, unlike covalent inter-protein linkages of the HK97 chainmail. Despite these complex interactions, the termini of both CP and MCP are ideally positioned for DGR-based phage-display engineering. DOI: http://dx.doi.org/10.7554/eLife.01299.001.

Surface display of a massively variable lipoprotein by a Legionella diversity-generating retroelement

Diversity-generating retroelements (DGRs) are a unique family of retroelements that confer selective advantages to their hosts by facilitating localized DNA sequence evolution through a specialized error-prone reverse transcription process. We characterized a DGR in Legionella pneumophila, an opportunistic human pathogen that causes Legionnaires disease. The L. pneumophila DGR is found within a horizontally acquired genomic island, and it can theoretically generate 10(26) unique nucleotide sequences in its target gene, legionella determinent target A (ldtA), creating a repertoire of 10(19) distinct proteins. Expression of the L. pneumophila DGR resulted in transfer of DNA sequence information from a template repeat to a variable repeat (VR) accompanied by adenine-specific mutagenesis of progeny VRs at the 3'end of ldtA. ldtA encodes a twin-arginine translocated lipoprotein that is anchored in the outer leaflet of the outer membrane, with its C-terminal variable region surface exposed. Related DGRs were identified in L. pneumophila clinical isolates that encode unique target proteins with homologous VRs, demonstrating the adaptability of DGR components. This work characterizes a DGR that diversifies a bacterial protein and confirms the hypothesis that DGR-mediated mutagenic homing occurs through a conserved mechanism. Comparative bioinformatics predicts that surface display of massively variable proteins is a defining feature of a subset of bacterial DGRs.

Insights into the physiology and ecology of the brackish-water-adapted Cyanobacterium Nodularia spumigena CCY9414 based on a genome-transcriptome analysis

Nodularia spumigena is a filamentous diazotrophic cyanobacterium that dominates the annual late summer cyanobacterial blooms in the Baltic Sea. But N. spumigena also is common in brackish water bodies worldwide, suggesting special adaptation allowing it to thrive at moderate salinities. A draft genome analysis of N. spumigena sp. CCY9414 yielded a single scaffold of 5,462,271 nucleotides in length on which genes for 5,294 proteins were annotated. A subsequent strand-specific transcriptome analysis identified more than 6,000 putative transcriptional start sites (TSS). Orphan TSSs located in intergenic regions led us to predict 764 non-coding RNAs, among them 70 copies of a possible retrotransposon and several potential RNA regulators, some of which are also present in other N2-fixing cyanobacteria. Approximately 4% of the total coding capacity is devoted to the production of secondary metabolites, among them the potent hepatotoxin nodularin, the linear spumigin and the cyclic nodulapeptin. The transcriptional complexity associated with genes involved in nitrogen fixation and heterocyst differentiation is considerably smaller compared to other Nostocales. In contrast, sophisticated systems exist for the uptake and assimilation of iron and phosphorus compounds, for the synthesis of compatible solutes, and for the formation of gas vesicles, required for the active control of buoyancy. Hence, the annotation and interpretation of this sequence provides a vast array of clues into the genomic underpinnings of the physiology of this cyanobacterium and indicates in particular a competitive edge of N. spumigena in nutrient-limited brackish water ecosystems.

Structure of the essential diversity-generating retroelement protein bAvd and its functionally important interaction with reverse transcriptase

Diversity-generating retroelements (DGRs) are the only known source of massive protein sequence variation in prokaryotes. These elements transfer coding information from a template region (TR) through an RNA intermediate to a protein-encoding variable region. This retrohoming process is accompanied by unique adenine-specific mutagenesis and, in the prototypical BPP-1 DGR, requires a reverse transcriptase (bRT) and an accessory variability determinant (bAvd) protein. To understand the role of bAvd, we determined its 2.69 Å resolution structure, which revealed a highly positively charged pentameric barrel. In accordance with its charge, bAvd bound both DNA and RNA, albeit without a discernable sequence preference. We found that the coding sequence of bAvd functioned as part of TR but identified means to mutate bAvd without affecting TR. This mutational analysis revealed a strict correspondence between retrohoming and interaction of bAvd with bRT, suggesting that the bRT-bAvd complex is important for DGR retrohoming.

Evolution of microbes and viruses: a paradigm shift in evolutionary biology?

When Charles Darwin formulated the central principles of evolutionary biology in the Origin of Species in 1859 and the architects of the Modern Synthesis integrated these principles with population genetics almost a century later, the principal if not the sole objects of evolutionary biology were multicellular eukaryotes, primarily animals and plants. Before the advent of efficient gene sequencing, all attempts to extend evolutionary studies to bacteria have been futile. Sequencing of the rRNA genes in thousands of microbes allowed the construction of the three- domain “ribosomal Tree of Life” that was widely thought to have resolved the evolutionary relationships between the cellular life forms. However, subsequent massive sequencing of numerous, complete microbial genomes revealed novel evolutionary phenomena, the most fundamental of these being: (1) pervasive horizontal gene transfer (HGT), in large part mediated by viruses and plasmids, that shapes the genomes of archaea and bacteria and call for a radical revision (if not abandonment) of the Tree of Life concept, (2) Lamarckian-type inheritance that appears to be critical for antivirus defense and other forms of adaptation in prokaryotes, and (3) evolution of evolvability, i.e., dedicated mechanisms for evolution such as vehicles for HGT and stress-induced mutagenesis systems. In the non-cellular part of the microbial world, phylogenomics and metagenomics of viruses and related selfish genetic elements revealed enormous genetic and molecular diversity and extremely high abundance of viruses that come across as the dominant biological entities on earth. Furthermore, the perennial arms race between viruses and their hosts is one of the defining factors of evolution. Thus, microbial phylogenomics adds new dimensions to the fundamental picture of evolution even as the principle of descent with modification discovered by Darwin and the laws of population genetics remain at the core of evolutionary biology.

Analysis of a comprehensive dataset of diversity generating retroelements generated by the program DiGReF

Background

Diversity Generating Retroelements (DGRs) are genetic cassettes that can introduce tremendous diversity into a short, defined region of the genome. They achieve hypermutation through replacement of the variable region with a strongly mutated cDNA copy generated by the element-encoded reverse transcriptase. In contrast to “selfish” retroelements such as group II introns and retrotransposons, DGRs impart an advantage to their host by increasing its adaptive potential. DGRs were discovered in a bacteriophage, but since then additional examples have been identified in some bacterial genomes.

Results

Here we present the program DiGReF that allowed us to comprehensively screen available databases for DGRs. We identified 155 DGRs which are found in all major classes of bacteria, though exhibiting sporadic distribution across species. Phylogenetic analysis and sequence comparison showed that DGRs move between genomes by associating with various mobile elements such as phages, transposons and plasmids. The DGR cassettes exhibit high flexibility in the arrangement of their components and easily acquire additional paralogous target genes. Surprisingly, the genomic data alone provide new insights into the molecular mechanism of DGRs. Most notably, our data suggest that the template RNA is transcribed separately from the rest of the element.

Conclusions

DiGReF is a valuable tool to detect DGRs in genome data. Its output allows comprehensive analysis of various aspects of DGR biology, thus deepening our understanding of the role DGRs play in prokaryotic genome plasticity, from the global down to the molecular level.

Hypervariable loci in the human gut virome

Genetic variation is critical in microbial immune evasion and drug resistance, but variation has rarely been studied in complex heterogeneous communities such as the human microbiome. To begin to study natural variation, we analyzed DNA viruses present in the lower gastrointestinal tract of 12 human volunteers by determining 48 billion bases of viral DNA sequence. Viral genomes mostly showed low variation, but 51 loci of ∼100 bp showed extremely high variation, so that up to 96% of the viral genomes encoded unique amino acid sequences. Some hotspots of hypervariation were in genes homologous to the bacteriophage BPP-1 viral tail-fiber gene, which is known to be hypermutagenized by a unique reverse-transcriptase (RT)-based mechanism. Unexpectedly, other hypervariable loci in our data were in previously undescribed gene types, including genes encoding predicted Ig-superfamily proteins. Most of the hypervariable loci were linked to genes encoding RTs of a single clade, which we find is the most abundant clade among gut viruses but only a minor component of bacterial RT populations. Hypervariation was targeted to 5'-AAY-3' asparagine codons, which allows maximal chemical diversification of the encoded amino acids while avoiding formation of stop codons. These findings document widespread targeted hypervariation in the human gut virome, identify previously undescribed types of genes targeted for hypervariation, clarify association with RT gene clades, and motivate studies of hypervariation in the full human microbiome.

Self-made phage libraries with heterologous inserts in the Mtd of Bordetella bronchiseptica

Phage display libraries are widely used as tools for identifying, dissecting and optimizing ligands. Development of a simple method to access greater library diversities could expedite and expand the technique. This paper reports progress toward harnessing the naturally occurring diversity generating retroelement used by Bordetella bronchiseptica bacteriophage to alter its tail-fiber protein. Mutagenesis and testing identified four sites amenable to the insertion of <19-residue heterologous peptides within the variable region. Such sites allow auto-generation of peptide libraries surrounded by a scaffold with additional variations. The resultant self-made phage libraries were used successfully for selections targeting anti-FLAG antibody, immobilized metal affinity chromatography microtiter plates and HIV-1 gp41. The reported experiments demonstrate the utility of the major tropism determinant protein of B.bronchiseptica as a natural scaffold for diverse, phage-constructed libraries with heterologous self-made phage libraries.

Target site recognition by a diversity-generating retroelement

Diversity-generating retroelements (DGRs) are in vivo sequence diversification machines that are widely distributed in bacterial, phage, and plasmid genomes. They function to introduce vast amounts of targeted diversity into protein-encoding DNA sequences via mutagenic homing. Adenine residues are converted to random nucleotides in a retrotransposition process from a donor template repeat (TR) to a recipient variable repeat (VR). Using the Bordetella bacteriophage BPP-1 element as a prototype, we have characterized requirements for DGR target site function. Although sequences upstream of VR are dispensable, a 24 bp sequence immediately downstream of VR, which contains short inverted repeats, is required for efficient retrohoming. The inverted repeats form a hairpin or cruciform structure and mutational analysis demonstrated that, while the structure of the stem is important, its sequence can vary. In contrast, the loop has a sequence-dependent function. Structure-specific nuclease digestion confirmed the existence of a DNA hairpin/cruciform, and marker coconversion assays demonstrated that it influences the efficiency, but not the site of cDNA integration. Comparisons with other phage DGRs suggested that similar structures are a conserved feature of target sequences. Using a kanamycin resistance determinant as a reporter, we found that transplantation of the IMH and hairpin/cruciform-forming region was sufficient to target the DGR diversification machinery to a heterologous gene. In addition to furthering our understanding of DGR retrohoming, our results suggest that DGRs may provide unique tools for directed protein evolution via in vivo DNA diversification.

Diversity-generating retroelements function through a unique, reverse transcriptase–mediated “copy and replace” mechanism that enables repeated rounds of protein diversification, selection, and optimization. The ability of DGRs to introduce targeted diversity into protein-coding DNA sequences has the potential to dramatically accelerate the evolution of adaptive traits. The utility of these elements in nature is underscored by their widespread distribution throughout the bacterial domain. Here we define DNA sequences and structures that are necessary and sufficient to direct the diversification machinery to specified target sequences. In addition to providing mechanistic insights into conserved features of DGR activity, our results provide a blueprint for the use of DGRs for a broad range of protein engineering applications.

Conservation of the C-type lectin fold for massive sequence variation in a Treponema diversity-generating retroelement

Anticipatory ligand binding through massive protein sequence variation is rare in biological systems, having been observed only in the vertebrate adaptive immune response and in a phage diversity-generating retroelement (DGR). Earlier work has demonstrated that the prototypical DGR variable protein, major tropism determinant (Mtd), meets the demands of anticipatory ligand binding by novel means through the C-type lectin (CLec) fold. However, because of the low sequence identity among DGR variable proteins, it has remained unclear whether the CLec fold is a general solution for DGRs. We have addressed this problem by determining the structure of a second DGR variable protein, TvpA, from the pathogenic oral spirochete Treponema denticola. Despite its weak sequence identity to Mtd (∼16%), TvpA was found to also have a CLec fold, with predicted variable residues exposed in a ligand-binding site. However, this site in TvpA was markedly more variable than the one in Mtd, reflecting the unprecedented approximate 10(20) potential variability of TvpA. In addition, similarity between TvpA and Mtd with formylglycine-generating enzymes was detected. These results provide strong evidence for the conservation of the formylglycine-generating enzyme-type CLec fold among DGRs as a means of accommodating massive sequence variation.