Comprehensive phylogenetic analysis of bacterial reverse transcriptases

Gut virome-wide association analysis identifies cross-population viral signatures for inflammatory bowel disease

BACKGROUND

The gut virome has been implicated in inflammatory bowel disease (IBD), yet a full understanding of the gut virome in IBD patients, especially across diverse geographic populations, is lacking.

RESULTS

In this study, we conducted a comprehensive gut virome-wide association study in a Chinese cohort of 71 IBD patients (15 with Crohn's disease and 56 with ulcerative colitis) and 77 healthy controls via viral-like particle (VLP) and bulk virome sequencing of their feces. By utilizing an integrated gut virus catalog tailored to the IBD virome, we revealed fundamental alterations in the gut virome in IBD patients. These characterized 139 differentially abundant viral signatures, including elevated phages predicted to infect Escherichia, Klebsiella, Enterococcus_B, Streptococcus, and Veillonella species, as well as IBD-depleted phages targeting Prevotella, Ruminococcus_E, Bifidobacterium, and Blautia species. Remarkably, these viral signatures demonstrated high consistency across diverse populations such as those in Europe and the USA, emphasizing their significance and broad relevance in the disease context. Furthermore, fecal virome transplantation experiments verified that the colonization of these IBD-characterized viruses can modulate experimental colitis in mouse models.

CONCLUSIONS

Building upon these insights into the IBD gut virome, we identified potential biomarkers for prognosis and therapy in IBD patients, laying the foundation for further exploration of viromes in related conditions. Video Abstract.

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.

Harnessing CRISPR-Cas adaptation for RNA recording and beyond

Prokaryotes encode clustered regularly interspaced short palindromic repeat (CRISPR) arrays and CRISPR-associated (Cas) genes as an adaptive immune machinery. CRISPR-Cas systems effectively protect hosts from the invasion of foreign enemies, such as bacteriophages and plasmids. During a process called 'adaptation', non-self-nucleic acid fragments are acquired as spacers between repeats in the host CRISPR array, to establish immunological memory. The highly conserved Cas1-Cas2 complexes function as molecular recorders to integrate spacers in a time course manner, which can subsequently be expressed as crRNAs complexed with Cas effector proteins for the RNAguided interference pathways. In some of the RNA-targeting type III systems, Cas1 proteins are fused with reverse transcriptase (RT), indicating that RT-Cas1-Cas2 complexes can acquire RNA transcripts for spacer acquisition. In this review, we summarize current studies that focus on the molecular structure and function of the RT-fused Cas1-Cas2 integrase, and its potential applications as a directional RNA-recording tool in cells. Furthermore, we highlight outstanding questions for RT-Cas1-Cas2 studies and future directions for RNA-recording CRISPR technologies. [BMB Reports 2024; 57(1): 40-49].

Integrating Phylogenetics With Intron Positions Illuminates the Origin of the Complex Spliceosome

Eukaryotic genes are characterized by the presence of introns that are removed from pre-mRNA by a spliceosome. This ribonucleoprotein complex is comprised of multiple RNA molecules and over a hundred proteins, which makes it one of the most complex molecular machines that originated during the prokaryote-to-eukaryote transition. Previous works have established that these introns and the spliceosomal core originated from self-splicing introns in prokaryotes. Yet, how the spliceosomal core expanded by recruiting many additional proteins remains largely elusive. In this study, we use phylogenetic analyses to infer the evolutionary history of 145 proteins that we could trace back to the spliceosome in the last eukaryotic common ancestor. We found that an overabundance of proteins derived from ribosome-related processes was added to the prokaryote-derived core. Extensive duplications of these proteins substantially increased the complexity of the emerging spliceosome. By comparing the intron positions between spliceosomal paralogs, we infer that most spliceosomal complexity postdates the spread of introns through the proto-eukaryotic genome. The reconstruction of early spliceosomal evolution provides insight into the driving forces behind the emergence of complexes with many proteins during eukaryogenesis.

Novel prokaryotic system employing previously unknown nucleic acids-based receptors

The present study describes a previously unknown universal system that orchestrates the interaction of bacteria with the environment, named the Teazeled receptor system (TR-system). The identical system was recently discovered within eukaryotes. The system includes DNA- and RNA-based molecules named "TezRs", that form receptor's network located outside the membrane, as well as reverse transcriptases and integrases. TR-system takes part in the control of all major aspects of bacterial behavior, such as intra cellular communication, growth, biofilm formation and dispersal, utilization of nutrients including xenobiotics, virulence, chemo- and magnetoreception, response to external factors (e.g., temperature, UV, light and gas content), mutation events, phage-host interaction, and DNA recombination activity. Additionally, it supervises the function of other receptor-mediated signaling pathways. Importantly, the TR-system is responsible for the formation and maintenance of cell memory to preceding cellular events, as well the ability to "forget" preceding events. Transcriptome and biochemical analysis revealed that the loss of different TezRs instigates significant alterations in gene expression and proteins synthesis.

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.

UG/Abi: a highly diverse family of prokaryotic reverse transcriptases associated with defense functions

Reverse transcriptases (RTs) are enzymes capable of synthesizing DNA using RNA as a template. Within the last few years, a burst of research has led to the discovery of novel prokaryotic RTs with diverse antiviral properties, such as DRTs (Defense-associated RTs), which belong to the so-called group of unknown RTs (UG) and are closely related to the Abortive Infection system (Abi) RTs. In this work, we performed a systematic analysis of UG and Abi RTs, increasing the number of UG/Abi members up to 42 highly diverse groups, most of which are predicted to be functionally associated with other gene(s) or domain(s). Based on this information, we classified these systems into three major classes. In addition, we reveal that most of these groups are associated with defense functions and/or mobile genetic elements, and demonstrate the antiphage role of four novel groups. Besides, we highlight the presence of one of these systems in novel families of human gut viruses infecting members of the Bacteroidetes and Firmicutes phyla. This work lays the foundation for a comprehensive and unified understanding of these highly diverse RTs with enormous biotechnological potential.

Structural Analysis of Monomeric RNA-Dependent Polymerases Revisited

In the past few years, our understanding of the RNA virosphere has changed dramatically due to the growth and spurt of metagenomics, exponentially increasing the number of RNA viral sequences, and providing a better understanding of their range of potential hosts. As of today, the only conserved protein among RNA viruses appears to be the monomeric RNA-dependent RNA polymerase. This enzyme belongs to the right-hand DNA-and RNA polymerases, which also includes reverse transcriptases and eukaryotic replicative DNA polymerases. The ubiquity of this protein in RNA viruses makes it a unique evolutionary marker and an appealing broad-spectrum antiviral target. In this work pairwise structural comparisons of viral RdRps and RTs were performed, including tertiary structures that have been obtained in the last few years. The resulting phylogenetic tree shows that the RdRps from (+)ss- and dsRNA viruses might have been recruited several times throughout the evolution of mobile genetic elements. RTs also display multiple evolutionary routes. We have identified a structural core comprising the entire palm, a large moiety of the fingers and the N-terminal helices of the thumb domain, comprising over 300 conserved residues, including two regions that we have named the “knuckles” and the “hypothenar eminence”. The conservation of an helix bundle in the region preceding the polymerase domain confirms that (−)ss and dsRNA Reoviruses’ polymerases share a recent ancestor. Finally, the inclusion of DNA polymerases into our structural analyses suggests that monomeric RNA-dependent polymerases might have diverged from B-family polymerases.

Phages and their satellites encode hotspots of antiviral systems

Bacteria carry diverse genetic systems to defend against viral infection, some of which are found within prophages where they inhibit competing viruses. Phage satellites pose additional pressures on phages by hijacking key viral elements to their own benefit. Here, we show that E. coli P2-like phages and their parasitic P4-like satellites carry hotspots of genetic variation containing reservoirs of anti-phage systems. We validate the activity of diverse systems and describe PARIS, an abortive infection system triggered by a phage-encoded anti-restriction protein. Antiviral hotspots participate in inter-viral competition and shape dynamics between the bacterial host, P2-like phages, and P4-like satellites. Notably, the anti-phage activity of satellites can benefit the helper phage during competition with virulent phages, turning a parasitic relationship into a mutualistic one. Anti-phage hotspots are present across distant species and constitute a substantial source of systems that participate in the competition between mobile genetic elements.

Two co-dominant nitrogen-fixing cyanobacteria demonstrate distinct acclimation and adaptation responses to cope with ocean warming

The globally dominant N₂ -fixing cyanobacteria Trichodesmium and Crocosphaera provide vital nitrogen supplies to subtropical and tropical oceans, but little is known about how they will be affected by long-term ocean warming. We tested their thermal responses using experimental evolution methods during 2 years of selection at optimal (28°C), supra-optimal (32°C) and suboptimal (22°C) temperatures. After several hundred generations under thermal selection, changes in growth parameters, as well as N and C fixation rates, suggested that Trichodesmium did not adapt to the three selection temperature regimes during the 2-year evolution experiment, but could instead rapidly and reversibly acclimate to temperature shifts from 20°C to 34°C. In contrast, over the same timeframe apparent thermal adaptation was observed in Crocosphaera, as evidenced by irreversible phenotypic changes as well as whole-genome sequencing and variant analysis. Especially under stressful warming conditions (34°C), 32°C-selected Crocosphaera cells had an advantage in survival and nitrogen fixation over cell lines selected at 22°C and 28°C. The distinct strategies of phenotypic plasticity versus irreversible adaptation in these two sympatric diazotrophs are both viable ways to maintain fitness despite long-term temperature changes, and so could help to stabilize key ocean nitrogen cycle functions under future warming scenarios.

Distinct Expansion of Group II Introns During Evolution of Prokaryotes and Possible Factors Involved in Its Regulation

Group II introns (G2Is) are ribozymes that have retroelement characteristics in prokaryotes. Although G2Is are suggested to have been an important evolutionary factor in the prokaryote-to-eukaryote transition, comprehensive analyses of these introns among the tens of thousands of prokaryotic genomes currently available are still limited. Here, we developed a bioinformatic pipeline that systematically collects G2Is and applied it to prokaryotic genomes. We found that in bacteria, 25% (447 of 1,790) of the total representative genomes had an average of 5.3 G2Is, and in archaea, 9% (28 of 296) of the total representative genomes had an average of 3.0 G2Is. The greatest number of G2Is per genome was 101 in Arthrospira platensis (phylum Cyanobacteriota). A comprehensive sequence analysis of the intron-encoded protein (IEP) in each G2I sequence was conducted and resulted in the addition of three new IEP classes (U1-U3) to the previous classification. This analysis suggested that about 30% of all IEPs are non-canonical IEPs. The number of G2Is per genome was defined almost at the phylum level, and at least in the following two phyla, Firmicutes, and Cyanobacteriota, the type of IEP was largely associated as a factor in the G2I increase, i.e., there was an explosive increase in G2Is with bacterial C-type IEPs, mainly in the phylum Firmicutes, and in G2Is with CL-type IEPs, mainly in the phylum Cyanobacteriota. We also systematically analyzed the relationship between genomic signatures and the mechanism of these increases in G2Is. This is the first study to systematically characterize G2Is in the prokaryotic phylogenies.

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.

Viruses Defined by the Position of the Virosphere within the Replicator Space

Originally, viruses were defined as miniscule infectious agents that passed through filters that retain even the smallest cells. Subsequently, viruses were considered obligate intracellular parasites whose reproduction depends on their cellular hosts for energy supply and molecular building blocks. However, these features are insufficient to unambiguously define viruses as they are broadly understood today. We outline possible approaches to define viruses and explore the boundaries of the virosphere within the virtual space of replicators and the relationships between viruses and other types of replicators. Regardless of how, exactly, viruses are defined, viruses clearly have evolved on many occasions from nonviral replicators, such as plasmids, by recruiting host proteins to become virion components. Conversely, other types of replicators have repeatedly evolved from viruses. Thus, the virosphere is a dynamic entity with extensive evolutionary traffic across its boundaries. We argue that the virosphere proper, here termed orthovirosphere, consists of a distinct variety of replicators that encode structural proteins encasing the replicators' genomes, thereby providing protection and facilitating transmission among hosts. Numerous and diverse replicators, such as virus-derived but capsidless RNA and DNA elements, or defective viruses occupy the zone surrounding the orthovirosphere in the virtual replicator space. We define this zone as the perivirosphere. Although intense debates on the nature of certain replicators that adorn the internal and external boundaries of the virosphere will likely continue, we present an operational definition of virus that recently has been accepted by the International Committee on Taxonomy of Viruses.

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.

Structural coordination between active sites of a CRISPR reverse transcriptase-integrase complex

CRISPR-Cas systems provide adaptive immunity in bacteria and archaea, beginning with integration of foreign sequences into the host CRISPR genomic locus and followed by transcription and maturation of CRISPR RNAs (crRNAs). In some CRISPR systems, a reverse transcriptase (RT) fusion to the Cas1 integrase and Cas6 maturase creates a single protein that enables concerted sequence integration and crRNA production. To elucidate how the RT-integrase organizes distinct enzymatic activities, we present the cryo-EM structure of a Cas6-RT-Cas1-Cas2 CRISPR integrase complex. The structure reveals a heterohexamer in which the RT directly contacts the integrase and maturase domains, suggesting functional coordination between all three active sites. Together with biochemical experiments, our data support a model of sequential enzymatic activities that enable CRISPR sequence acquisition from RNA and DNA substrates. These findings highlight an expanded capacity of some CRISPR systems to acquire diverse sequences that direct CRISPR-mediated interference.

Bacterial Group II Intron Genomic Neighborhoods Reflect Survival Strategies: Hiding and Hijacking

Group II (gII) introns are mobile retroelements that can spread to new DNA sites through retrotransposition, which can be influenced by a variety of host factors. To determine if these host factors bear any relationship to the genomic location of gII introns, we developed a bioinformatic pipeline wherein we focused on the genomic neighborhoods of bacterial gII introns within their native contexts and sought to determine global relationships between introns and their surrounding genes. We found that, although gII introns inhabit diverse regions, these neighborhoods are often functionally enriched for genes that could promote gII intron retention or proliferation. On one hand, we observe that gII introns are frequently found hiding in mobile elements or after transcription terminators. On the other hand, gII introns are enriched in locations in which they could hijack host functions for their movement, potentially timing expression of the intron with genes that produce favorable conditions for retrotransposition. Thus, we propose that gII intron distributions have been shaped by relationships with their surrounding genomic neighbors.

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.

Diverse enzymatic activities mediate antiviral immunity in prokaryotes

Bacteria and archaea are frequently attacked by viruses and other mobile genetic elements and rely on dedicated antiviral defense systems, such as restriction endonucleases and CRISPR, to survive. The enormous diversity of viruses suggests that more types of defense systems exist than are currently known. By systematic defense gene prediction and heterologous reconstitution, here we discover 29 widespread antiviral gene cassettes, collectively present in 32% of all sequenced bacterial and archaeal genomes, that mediate protection against specific bacteriophages. These systems incorporate enzymatic activities not previously implicated in antiviral defense, including RNA editing and retron satellite DNA synthesis. In addition, we computationally predict a diverse set of other putative defense genes that remain to be characterized. These results highlight an immense array of molecular functions that microbes use against viruses.

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.

CRISPR-Cas Systems and the Paradox of Self-Targeting Spacers

CRISPR-Cas immune systems in bacteria and archaea record prior infections as spacers within each system’s CRISPR arrays. Spacers are normally derived from invasive genetic material and direct the immune system to complementary targets as part of future infections. However, not all spacers appear to be derived from foreign genetic material and instead can originate from the host genome. Their presence poses a paradox, as self-targeting spacers would be expected to induce an autoimmune response and cell death. In this review, we discuss the known frequency of self-targeting spacers in natural CRISPR-Cas systems, how these spacers can be incorporated into CRISPR arrays, and how the host can evade lethal attack. We also discuss how self-targeting spacers can become the basis for alternative functions performed by CRISPR-Cas systems that extend beyond adaptive immunity. Overall, the acquisition of genome-targeting spacers poses a substantial risk but can aid in the host’s evolution and potentially lead to or support new functionalities.

Spacer acquisition from RNA mediated by a natural reverse transcriptase-Cas1 fusion protein associated with a type III-D CRISPR-Cas system in Vibrio vulnificus

The association of reverse transcriptases (RTs) with CRISPR–Cas system has recently attracted interest because the RT activity appears to facilitate the RT-dependent acquisition of spacers from RNA molecules. However, our understanding of this spacer acquisition process remains limited. We characterized the in vivo acquisition of spacers mediated by an RT-Cas1 fusion protein linked to a type III-D system from Vibrio vulnificus strain YJ016, and showed that the adaptation module, consisting of the RT-Cas1 fusion, two different Cas2 proteins (A and B) and one of the two CRISPR arrays, was completely functional in a heterologous host. We found that mutations of the active site of the RT domain significantly decreased the acquisition of new spacers and showed that this RT-Cas1-associated adaptation module was able to incorporate spacers from RNA molecules into the CRISPR array. We demonstrated that the two Cas2 proteins of the adaptation module were required for spacer acquisition. Furthermore, we found that several sequence-specific features were required for the acquisition and integration of spacers derived from any region of the genome, with no bias along the 5′and 3′ends of coding sequences. This study provides new insight into the RT-Cas1 fusion protein-mediated acquisition of spacers from RNA molecules.

Spacer acquisition from RNA mediated by a natural reverse transcriptase-Cas1 fusion protein associated with a type III-D CRISPR-Cas system in Vibrio vulnificus

The association of reverse transcriptases (RTs) with CRISPR-Cas system has recently attracted interest because the RT activity appears to facilitate the RT-dependent acquisition of spacers from RNA molecules. However, our understanding of this spacer acquisition process remains limited. We characterized the in vivo acquisition of spacers mediated by an RT-Cas1 fusion protein linked to a type III-D system from Vibrio vulnificus strain YJ016, and showed that the adaptation module, consisting of the RT-Cas1 fusion, two different Cas2 proteins (A and B) and one of the two CRISPR arrays, was completely functional in a heterologous host. We found that mutations of the active site of the RT domain significantly decreased the acquisition of new spacers and showed that this RT-Cas1-associated adaptation module was able to incorporate spacers from RNA molecules into the CRISPR array. We demonstrated that the two Cas2 proteins of the adaptation module were required for spacer acquisition. Furthermore, we found that several sequence-specific features were required for the acquisition and integration of spacers derived from any region of the genome, with no bias along the 5'and 3'ends of coding sequences. This study provides new insight into the RT-Cas1 fusion protein-mediated acquisition of spacers from RNA molecules.

Large Plasmid Complement Resolved: Complete Genome Sequencing of Lactobacillus plantarum MF1298, a Candidate Probiotic Strain Associated with Unfavorable Effect

Considerable attention has been given to the species Lactobacillus plantarum regarding its probiotic potential. L. plantarum strains have shown health benefits in several studies, and even nonstrain-specific claims are allowed in certain markets. L. plantarum strain MF1298 was considered a candidate probiotic, demonstrating in vitro probiotic properties and the ability to survive passage through the human intestinal tract. However, the strain showed an unfavorable effect on symptoms in subjects with irritable bowel syndrome in a clinical trial. The properties and the genome of this strain are thus of general interest. Obtaining the complete genome of strain MF1298 proved difficult due to its large plasmid complement. Here, we exploit a combination of sequencing approaches to obtain the complete chromosome and plasmid assemblies of MF1298. The Oxford Nanopore Technologies MinION long-read sequencer was particularly useful in resolving the unusually large number of plasmids in the strain, 14 in total. The complete genome sequence of 3,576,440 basepairs contains 3272 protein-encoding genes, of which 315 are located on plasmids. Few unique regions were found in comparison with other L. plantarum genomes. Notably, however, one of the plasmids contains genes related to vitamin B12 (cobalamin) turnover and genes encoding bacterial reverse transcriptases, features not previously reported for L. plantarum. The extensive plasmid information will be important for future studies with this strain.

Multiple origins of reverse transcriptases linked to CRISPR-Cas systems

Prokaryotic genomes harbour a plethora of uncharacterized reverse transcriptases (RTs). RTs phylogenetically related to those encoded by group-II introns have been found associated with type III CRISPR-Cas systems, adjacent or fused at the C-terminus to Cas1. It is thought that these RTs may have a relevant function in the CRISPR immune response mediating spacer acquisition from RNA molecules. The origin and relationships of these RTs and the ways in which the various protein domains evolved remain matters of debate. We carried out a large survey of annotated RTs in databases (198,760 sequences) and constructed a large dataset of unique representative sequences (9,141). The combined phylogenetic reconstruction and identification of the RTs and their various protein domains in the vicinity of CRISPR adaptation and effector modules revealed three different origins for these RTs, consistent with their emergence on multiple occasions: a larger group that have evolved from group-II intron RTs, and two minor lineages that may have arisen more recently from Retron/retron-like sequences and Abi-P2 RTs, the latter associated with type I-C systems. We also identified a particular group of RTs associated with CRISPR-cas loci in clade 12, fused C-terminally to an archaeo-eukaryotic primase (AEP), a protein domain (AE-Prim_S_like) forming a particular family within the AEP proper clade. Together, these data provide new insight into the evolution of CRISPR-Cas/RT systems.

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).

DNA event recorders send past information of cells to the time of observation

While current omics and single cell technologies have enabled measurements of high-resolution molecular snapshots of cells at a large scale, these technologies all require destruction of samples and prevent us from analyzing dynamic changes in molecular profiles, phenotypes, and behaviors of individual cells in a complex system. One possible direction to overcome this issue is the development of a cell-embedded 'event recorder' system, whereby molecular and phenotypic information of a cell(s) can be obtained at the time of observation with their past event information stored in 'heritable polymers' of the same cell. This concept has been demonstrated by many synthetic cellular circuits that monitor and transmit a certain set of environmental and intracellular signals into DNA, and have now been further accelerated by recent CRISPR-related technologies. Notably, the discovery of the RT-Cas1-Cas2 system, which acquires sequences of cellular transcripts into a specific host genomic region, has enabled recording of a broader range of molecular profile histories in the DNA tapes of cells, to understand the dynamics of complex biological processes that cannot be addressed by current technologies.

Origins and evolution of CRISPR-Cas systems

CRISPR-Cas, the bacterial and archaeal adaptive immunity systems, encompass a complex machinery that integrates fragments of foreign nucleic acids, mostly from mobile genetic elements (MGE), into CRISPR arrays embedded in microbial genomes. Transcripts of the inserted segments (spacers) are employed by CRISPR-Cas systems as guide (g)RNAs for recognition and inactivation of the cognate targets. The CRISPR-Cas systems consist of distinct adaptation and effector modules whose evolutionary trajectories appear to be at least partially independent. Comparative genome analysis reveals the origin of the adaptation module from casposons, a distinct type of transposons, which employ a homologue of Cas1 protein, the integrase responsible for the spacer incorporation into CRISPR arrays, as the transposase. The origin of the effector module(s) is far less clear. The CRISPR-Cas systems are partitioned into two classes, class 1 with multisubunit effectors, and class 2 in which the effector consists of a single, large protein. The class 2 effectors originate from nucleases encoded by different MGE, whereas the origin of the class 1 effector complexes remains murky. However, the recent discovery of a signalling pathway built into the type III systems of class 1 might offer a clue, suggesting that type III effector modules could have evolved from a signal transduction system involved in stress-induced programmed cell death. The subsequent evolution of the class 1 effector complexes through serial gene duplication and displacement, primarily of genes for proteins containing RNA recognition motif domains, can be hypothetically reconstructed. In addition to the multiple contributions of MGE to the evolution of CRISPR-Cas, the reverse flow of information is notable, namely, recruitment of minimalist variants of CRISPR-Cas systems by MGE for functions that remain to be elucidated. Here, we attempt a synthesis of the diverse threads that shed light on CRISPR-Cas origins and evolution. This article is part of a discussion meeting issue 'The ecology and evolution of prokaryotic CRISPR-Cas adaptive immune systems'.

Transposons: the agents of antibiotic resistance in bacteria

Transposons are a group of mobile genetic elements that are defined as a DNA sequence. Transposons can jump into different places of the genome; for this reason, they are called jumping genes. However, some transposons are always kept at the insertion site in the genome. Most transposons are inactivated and as a result, cannot move. Transposons are divided into two main groups: retrotransposons (class І) and DNA transposons (class ІІ). Retrotransposons are often found in eukaryotes. DNA transposons can be found in both eukaryotes and prokaryotes. The bacterial transposons belong to the DNA transposons and the Tn family, which are usually the carrier of additional genes for antibiotic resistance. Transposons can transfer from a plasmid to other plasmids or from a DNA chromosome to plasmid and vice versa that cause the transmission of antibiotic resistance genes in bacteria. The treatment of bacterial infectious diseases is difficult because of existing antibiotic resistance that part of this antibiotic resistance is caused by transposons. Bacterial infectious diseases are responsible for the increasing rise in world mortality rate. In this review, transposons and their roles have been studied in bacterial antibiotic resistance, in detail.

Comparison between complete genomes of an isolate of Pseudomonas syringae pv. actinidiae from Japan and a New Zealand isolate of the pandemic lineage

The modern pandemic of the bacterial kiwifruit pathogen Pseudomonas syringae pv actinidiae (Psa) is caused by a particular Psa lineage. To better understand the genetic basis of the virulence of this lineage, we compare the completely assembled genome of a pandemic New Zealand strain with that of the Psa type strain first isolated in Japan in 1983. Aligning the two genomes shows numerous translocations, constrained so as to retain the appropriate orientation of the Architecture Imparting Sequences (AIMs). There are several large horizontally acquired regions, some of which include Type I, Type II or Type III restriction systems. The activity of these systems is reflected in the methylation patterns of the two strains. The pandemic strain carries an Integrative Conjugative Element (ICE) located at a tRNA-Lys site. Two other complex elements are also present at tRNA-Lys sites in the genome. These elements are derived from ICE but have now acquired some alternative secretion function. There are numerous types of mobile element in the two genomes. Analysis of these elements reveals no evidence of recombination between the two Psa lineages.

On the Origin and Evolutionary Relationships of the Reverse Transcriptases Associated With Type III CRISPR-Cas Systems

Reverse transcriptases (RTs) closely related to those encoded by group II introns but lacking the intron RNA structure have been found associated with type III clustered regularly interspaced short palindromic repeats (CRISPR)-Cas systems, a prokaryotic immune system against invading viruses and foreign genetic elements. Two models have been proposed to explain the origin and evolutionary relationships of these RTs: (i) the “single point of origin” model, according to which these RTs originated from a single acquisition event in bacterial, with the various protein domains (RT, RT-Cas1, and Cas6-RT-Cas1 fusions) corresponding to single points in evolution; and (ii) the “various origins” model, according to which, independent acquisition events in different evolutionary episodes led to these fusions. We tested these alternative hypotheses, by analyzing and integrating published datasets of RT sequences associated with CRISPR-Cas systems and inferring phylogenetic trees by maximum likelihood (ML) methods. The RTs studied could be grouped into 13 clades, mostly in bacteria, in which they probably evolved. The various clades appear to form three independent lineages in bacteria and a recent lineage in archaea. Our data show that the Cas6 domain was acquired twice, independently, through RT-Cas1 fusion, in the bacterial lineages. Taken together, there more evidence to support the “various origins” hypothesis.

Using bioinformatic and phylogenetic approaches to classify transposable elements and understand their complex evolutionary histories

In recent years, much attention has been paid to comparative genomic studies of transposable elements (TEs) and the ensuing problems of their identification, classification, and annotation. Different approaches and diverse automated pipelines are being used to catalogue and categorize mobile genetic elements in the ever-increasing number of prokaryotic and eukaryotic genomes, with little or no connectivity between different domains of life. Here, an overview of the current picture of TE classification and evolutionary relationships is presented, updating the diversity of TE types uncovered in sequenced genomes. A tripartite TE classification scheme is proposed to account for their replicative, integrative, and structural components, and the need to expand in vitro and in vivo studies of their structural and biological properties is emphasized. Bioinformatic studies have now become front and center of novel TE discovery, and experimental pursuits of these discoveries hold great promise for both basic and applied science.

Viruses and mobile elements as drivers of evolutionary transitions

The history of life is punctuated by evolutionary transitions which engender emergence of new levels of biological organization that involves selection acting at increasingly complex ensembles of biological entities. Major evolutionary transitions include the origin of prokaryotic and then eukaryotic cells, multicellular organisms and eusocial animals. All or nearly all cellular life forms are hosts to diverse selfish genetic elements with various levels of autonomy including plasmids, transposons and viruses. I present evidence that, at least up to and including the origin of multicellularity, evolutionary transitions are driven by the coevolution of hosts with these genetic parasites along with sharing of ‘public goods’. Selfish elements drive evolutionary transitions at two distinct levels. First, mathematical modelling of evolutionary processes, such as evolution of primitive replicator populations or unicellular organisms, indicates that only increasing organizational complexity, e.g. emergence of multicellular aggregates, can prevent the collapse of the host–parasite system under the pressure of parasites. Second, comparative genomic analysis reveals numerous cases of recruitment of genes with essential functions in cellular life forms, including those that enable evolutionary transitions.

This article is part of the themed issue ‘The major synthetic evolutionary transitions’.

Mobile Genetic Elements and Evolution of CRISPR-Cas Systems: All the Way There and Back

The Clustered Regularly Interspaced Palindromic Repeats (CRISPR)-CRISPR-associated proteins (Cas) systems of bacterial and archaeal adaptive immunity show multifaceted evolutionary relationships with at least five classes of mobile genetic elements (MGE). First, the adaptation module of CRISPR-Cas that is responsible for the formation of the immune memory apparently evolved from a Casposon, a self-synthesizing transposon that employs the Cas1 protein as the integrase and might have brought additional cas genes to the emerging immunity loci. Second, a large subset of type III CRISPR-Cas systems recruited a reverse transcriptase from a Group II intron, providing for spacer acquisition from RNA. Third, effector nucleases of Class 2 CRISPR-Cas systems that are responsible for the recognition and cleavage of the target DNA were derived from transposon-encoded TnpB nucleases, most likely, on several independent occasions. Fourth, accessory nucleases in some variants of types I and III toxin and type VI effectors RNases appear to be ultimately derived from toxin nucleases of microbial toxin-antitoxin modules. Fifth, the opposite direction of evolution is manifested in the recruitment of CRISPR-Cas systems by a distinct family of Tn7-like transposons that probably exploit the capacity of CRISPR-Cas to recognize unique DNA sites to facilitate transposition as well as by bacteriophages that employ them to cope with host defense. Additionally, individual Cas proteins, such as the Cas4 nuclease, were recruited by bacteriophages and transposons. The two-sided evolutionary connection between CRISPR-Cas and MGE fits the "guns for hire" paradigm whereby homologous enzymatic machineries, in particular nucleases, are shuttled between MGE and defense systems and are used alternately as means of offense or defense.

A re-sequencing-based ultra-dense genetic map reveals a gummy stem blight resistance-associated gene in Cucumis melo

The melon (Cucumis melo) genome and genetic maps with hundreds to thousands of single nucleotide polymorphism markers were recently released. However, a high-resolution genetic map was lacking. Gummy stem blight (Gsb) is a destructive disease responsible for considerable economic losses during melon production. We herein describe the development of an ultra-dense genetic map consisting of 12,932 recombination bin markers covering 1,818 cM, with an average distance of 0.17 cM between adjacent tags. A comparison of the genetic maps for melon, watermelon, and cucumber revealed chromosome-level syntenic relationships and recombination events among the three Cucurbitaceae species. Our genetic map was useful for re-anchoring the genome scaffolds of melon. More than 92% assembly was anchored to 12 pseudo-chromosomes and 90% of them were oriented. Furthermore, 1,135 recombination hotspots revealed an unbalanced recombination rate across the melon genome. Genetic analyses of the Gsb-resistant and -susceptible lines indicated the resistance phenotype is mediated by a single dominant gene. We identified Gsb-resistance gene candidates in a 108-kb region on pseudo-chromosome 4. Our findings verify the utility of an ultra-dense genetic map for mapping a gene of interest, and for identifying new disease resistant genes.

Mobile Group II Introns as Ancestral Eukaryotic Elements

The duality of group II introns, capable of carrying out both self-splicing and retromobility reactions, is hypothesized to have played a profound role in the evolution of eukaryotes. These introns likely provided the framework for the emergence of eukaryotic retroelements, spliceosomal introns and other key components of the spliceosome. Group II introns are found in all three domains of life and are therefore considered to be exceptionally successful mobile genetic elements. Initially identified in organellar genomes, group II introns are found in bacteria, chloroplasts, and mitochondria of plants and fungi, but not in nuclear genomes. Although there is no doubt that prokaryotic and organellar group II introns are evolutionary related, there are remarkable differences in survival strategies between them. Furthermore, an evolutionary relationship of group II introns to eukaryotic retroelements, including telomeres, and spliceosomes is unmistakable.

The Reverse Transcriptases Associated with CRISPR-Cas Systems

CRISPR (clustered regularly interspaced short palindromic repeats) and associated proteins (Cas) act as adaptive immune systems in bacteria and archaea. Some CRISPR-Cas systems have been found to be associated with putative reverse transcriptases (RT), and an RT-Cas1 fusion associated with a type III-B system has been shown to acquire RNA spacers in vivo. Nevertheless, the origin and evolutionary relationships of these RTs and associated CRISPR-Cas systems remain largely unknown. We performed a comprehensive phylogenetic analysis of these RTs and associated Cas1 proteins, and classified their CRISPR-Cas modules. These systems were found predominantly in bacteria, and their presence in archaea may be due to a horizontal gene transfer event. These RTs cluster into 12 major clades essentially restricted to particular phyla, suggesting host-dependent functioning. The RTs and associated Cas1 proteins may have largely coevolved. They are, therefore, subject to the same selection pressures, which may have led to coadaptation within particular protein complexes. Furthermore, our results indicate that the association of an RT with a CRISPR-Cas system has occurred on multiple occasions during evolution.

The Reverse Transcriptases Associated with CRISPR-Cas Systems

CRISPR (clustered regularly interspaced short palindromic repeats) and associated proteins (Cas) act as adaptive immune systems in bacteria and archaea. Some CRISPR-Cas systems have been found to be associated with putative reverse transcriptases (RT), and an RT-Cas1 fusion associated with a type III-B system has been shown to acquire RNA spacers in vivo. Nevertheless, the origin and evolutionary relationships of these RTs and associated CRISPR-Cas systems remain largely unknown. We performed a comprehensive phylogenetic analysis of these RTs and associated Cas1 proteins, and classified their CRISPR-Cas modules. These systems were found predominantly in bacteria, and their presence in archaea may be due to a horizontal gene transfer event. These RTs cluster into 12 major clades essentially restricted to particular phyla, suggesting host-dependent functioning. The RTs and associated Cas1 proteins may have largely coevolved. They are, therefore, subject to the same selection pressures, which may have led to coadaptation within particular protein complexes. Furthermore, our results indicate that the association of an RT with a CRISPR-Cas system has occurred on multiple occasions during evolution.

On the Origin of Reverse Transcriptase-Using CRISPR-Cas Systems and Their Hyperdiverse, Enigmatic Spacer Repertoires

Cas1 integrase is the key enzyme of the clustered regularly interspaced short palindromic repeat (CRISPR)-Cas adaptation module that mediates acquisition of spacers derived from foreign DNA by CRISPR arrays. In diverse bacteria, the cas1 gene is fused (or adjacent) to a gene encoding a reverse transcriptase (RT) related to group II intron RTs. An RT-Cas1 fusion protein has been recently shown to enable acquisition of CRISPR spacers from RNA. Phylogenetic analysis of the CRISPR-associated RTs demonstrates monophyly of the RT-Cas1 fusion, and coevolution of the RT and Cas1 domains. Nearly all such RTs are present within type III CRISPR-Cas loci, but their phylogeny does not parallel the CRISPR-Cas type classification, indicating that RT-Cas1 is an autonomous functional module that is disseminated by horizontal gene transfer and can function with diverse type III systems. To compare the sequence pools sampled by RT-Cas1-associated and RT-lacking CRISPR-Cas systems, we obtained samples of a commercially grown cyanobacterium—Arthrospira platensis. Sequencing of the CRISPR arrays uncovered a highly diverse population of spacers. Spacer diversity was particularly striking for the RT-Cas1-containing type III-B system, where no saturation was evident even with millions of sequences analyzed. In contrast, analysis of the RT-lacking type III-D system yielded a highly diverse pool but reached a point where fewer novel spacers were recovered as sequencing depth was increased. Matches could be identified for a small fraction of the non-RT-Cas1-associated spacers, and for only a single RT-Cas1-associated spacer. Thus, the principal source(s) of the spacers, particularly the hypervariable spacer repertoire of the RT-associated arrays, remains unknown.

While the majority of CRISPR-Cas immune systems adapt to foreign genetic elements by capturing segments of invasive DNA, some systems carry reverse transcriptases (RTs) that enable adaptation to RNA molecules. From analysis of available bacterial sequence data, we find evidence that RT-based RNA adaptation machinery has been able to join with CRISPR-Cas immune systems in many, diverse bacterial species. To investigate whether the abilities to adapt to DNA and RNA molecules are utilized for defense against distinct classes of invaders in nature, we sequenced CRISPR arrays from samples of commercial-scale open-air cultures of Arthrospira platensis, a cyanobacterium that contains both RT-lacking and RT-containing CRISPR-Cas systems. We uncovered a diverse pool of naturally occurring immune memories, with the RT-lacking locus acquiring a number of segments matching known viral or bacterial genes, while the RT-containing locus has acquired spacers from a distinct sequence pool for which the source remains enigmatic.

The diversity of mtDNA rns introns among strains of Ophiostoma piliferum, Ophiostoma pluriannulatum and related species

Background

Based on previous studies, it was suspected that the mitochondrial rns gene within the Ophiostomatales is rich in introns. This study focused on a collection of strains representing Ophiostoma piliferum, Ophiostoma pluriannulatum and related species that cause blue-stain; these fungi colonize the sapwood of trees and impart a dark stain. This reduces the value of the lumber. The goal was to examine the mtDNA rns intron landscape for these important blue stain fungi in order to facilitate future annotation of mitochondrial genomes (mtDNA) and to potentially identify mtDNA introns that can encode homing endonucleases which may have applications in biotechnology.

Results

Comparative sequence analysis identified five intron insertion sites among the ophiostomatoid fungi examined. Positions mS379 and mS952 harbor group II introns, the mS379 intron encodes a reverse transcriptase, and the mS952 intron encodes a potential homing endonuclease. Positions mS569, mS1224, and mS1247 have group I introns inserted and these encode intact or eroded homing endonuclease open reading frames (ORF). Phylogenetic analysis of the intron ORFs showed that they can be found in the same insertion site in closely and distantly related species.

Conclusions

Based on the molecular markers examined (rDNA internal transcribed spacers and rns introns), strains representing O. pilifera, O. pluriannulatum and Ophiostoma novae-zelandiae could not be resolved. Phylogenetic studies suggest that introns are gained and lost and that horizontal transfer could explain the presence of related intron in distantly related fungi. With regard to the mS379 group II intron, this study shows that mitochondrial group II introns and their reverse transcriptases may also follow the life cycle previously proposed for group I introns and their homing endonucleases. This consists of intron invasion, decay of intron ORF, loss of intron, and possible reinvasion.

Electronic supplementary material

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

The early events underlying genome evolution in a localized Sinorhizobium meliloti population

BACKGROUND

Population genetic analyses based on genome-wide sequencing data have been carried out for Sinorhizobium medicae and S. meliloti, two closely related bacterial species forming nitrogen-fixing symbioses with plants of the genus Medicago. However, genome coverage was low or the isolates had a broad geographic distribution, making it difficult to interpret the estimated diversity and to unravel the early events underlying population genetic variations and ecological differentiation.

RESULTS

Here, to gain insight into the early genome level variation and diversification within S. meliloti populations, we first used Illumina paired-end reads technology to sequence a new clone of S. meliloti strain GR4, a highly competitive strain for alfalfa nodulation. The Illumina data and the GR4 genome sequence previously obtained with 454 technology were used to generate a high-quality reference genome sequence. We then used Illumina technology to sequence the genomes of 13 S. meliloti isolates representative of the genomic variation within the GR4-type population, obtained from a single field site with a high degree of coverage. The genome sequences obtained were analyzed to determine nucleotide diversity, divergence times, polymorphism and genomic variation. Similar low levels of nucleotide diversity were observed for the chromosome, pSymB and pSymA replicons. The isolates displayed other types of variation, such as indels, recombination events, genomic island excision and the transposition of mobile elements.

CONCLUSIONS

Our results suggest that the GR4-type population has experienced a process of demographic expansion and behaves as a stable genotypic cluster of genome-wide similarity, with most of the genome following a clonal pattern of evolution. Although some of genetic variation detected within the GR4-type population is probably due to genetic drift, others might be important in diversification and environmental adaptation.

Direct CRISPR spacer acquisition from RNA by a natural reverse transcriptase-Cas1 fusion protein

CRISPR systems mediate adaptive immunity in diverse prokaryotes. CRISPR-associated Cas1 and Cas2 proteins have been shown to enable adaptation to new threats in type I and II CRISPR systems by the acquisition of short segments of DNA (spacers) from invasive elements. In several type III CRISPR systems, Cas1 is naturally fused to a reverse transcriptase (RT). In the marine bacterium Marinomonas mediterranea (MMB-1), we showed that a RT-Cas1 fusion protein enables the acquisition of RNA spacers in vivo in a RT-dependent manner. In vitro, the MMB-1 RT-Cas1 and Cas2 proteins catalyze the ligation of RNA segments into the CRISPR array, which is followed by reverse transcription. These observations outline a host-mediated mechanism for reverse information flow from RNA to DNA.

Bacterial Group II Introns: Identification and Mobility Assay

Group II introns are large catalytic RNAs and mobile retroelements that encode a reverse transcriptase. Here, we provide methods for their identification in bacterial genomes and further analysis of their splicing and mobility capacities.

Reverse transcriptase genes are highly abundant and transcriptionally active in marine plankton assemblages

Genes encoding reverse transcriptases (RTs) are found in most eukaryotes, often as a component of retrotransposons, as well as in retroviruses and in prokaryotic retroelements. We investigated the abundance, classification and transcriptional status of RTs based on Tara Oceans marine metagenomes and metatranscriptomes encompassing a wide organism size range. Our analyses revealed that RTs predominate large-size fraction metagenomes (>5 μm), where they reached a maximum of 13.5% of the total gene abundance. Metagenomic RTs were widely distributed across the phylogeny of known RTs, but many belonged to previously uncharacterized clades. Metatranscriptomic RTs showed distinct abundance patterns across samples compared with metagenomic RTs. The relative abundances of viral and bacterial RTs among identified RT sequences were higher in metatranscriptomes than in metagenomes and these sequences were detected in all metatranscriptome size fractions. Overall, these observations suggest an active proliferation of various RT-assisted elements, which could be involved in genome evolution or adaptive processes of plankton assemblage.

Origins and evolution of viruses of eukaryotes: The ultimate modularity

Viruses and other selfish genetic elements are dominant entities in the biosphere, with respect to both physical abundance and genetic diversity. Various selfish elements parasitize on all cellular life forms. The relative abundances of different classes of viruses are dramatically different between prokaryotes and eukaryotes. In prokaryotes, the great majority of viruses possess double-stranded (ds) DNA genomes, with a substantial minority of single-stranded (ss) DNA viruses and only limited presence of RNA viruses. In contrast, in eukaryotes, RNA viruses account for the majority of the virome diversity although ssDNA and dsDNA viruses are common as well. Phylogenomic analysis yields tangible clues for the origins of major classes of eukaryotic viruses and in particular their likely roots in prokaryotes. Specifically, the ancestral genome of positive-strand RNA viruses of eukaryotes might have been assembled de novo from genes derived from prokaryotic retroelements and bacteria although a primordial origin of this class of viruses cannot be ruled out. Different groups of double-stranded RNA viruses derive either from dsRNA bacteriophages or from positive-strand RNA viruses. The eukaryotic ssDNA viruses apparently evolved via a fusion of genes from prokaryotic rolling circle-replicating plasmids and positive-strand RNA viruses. Different families of eukaryotic dsDNA viruses appear to have originated from specific groups of bacteriophages on at least two independent occasions. Polintons, the largest known eukaryotic transposons, predicted to also form virus particles, most likely, were the evolutionary intermediates between bacterial tectiviruses and several groups of eukaryotic dsDNA viruses including the proposed order "Megavirales" that unites diverse families of large and giant viruses. Strikingly, evolution of all classes of eukaryotic viruses appears to have involved fusion between structural and replicative gene modules derived from different sources along with additional acquisitions of diverse genes.

Recent mobility of plastid encoded group II introns and twintrons in five strains of the unicellular red alga Porphyridium

Group II introns are closely linked to eukaryote evolution because nuclear spliceosomal introns and the small RNAs associated with the spliceosome are thought to trace their ancient origins to these mobile elements. Therefore, elucidating how group II introns move, and how they lose mobility can potentially shed light on fundamental aspects of eukaryote biology. To this end, we studied five strains of the unicellular red alga Porphyridium purpureum that surprisingly contain 42 group II introns in their plastid genomes. We focused on a subset of these introns that encode mobility-conferring intron-encoded proteins (IEPs) and found them to be distributed among the strains in a lineage-specific manner. The reverse transcriptase and maturase domains were present in all lineages but the DNA endonuclease domain was deleted in vertically inherited introns, demonstrating a key step in the loss of mobility. P. purpureum plastid intron RNAs had a classic group IIB secondary structure despite variability in the DIII and DVI domains. We report for the first time the presence of twintrons (introns-within-introns, derived from the same mobile element) in Rhodophyta. The P. purpureum IEPs and their mobile introns provide a valuable model for the study of mobile retroelements in eukaryotes and offer promise for biotechnological applications.