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.

Group II introns: versatile ribozymes and retroelements.

Group II introns are catalytic RNAs (ribozymes) and retroelements found in the genomes of bacteria, archaebacteria, and organelles of some eukaryotes. The prototypical retroelement form consists of a structurally conserved RNA and a multidomain reverse transcriptase protein, which interact with each other to mediate splicing and mobility reactions. A wealth of biochemical, cross-linking, and X-ray crystal structure studies have helped to reveal how the two components cooperate to carry out the splicing and mobility reactions. In addition to the standard retroelement form, group II introns have evolved into derivative forms by either losing specific splicing or mobility characteristics, or becoming functionally specialized. Of particular interest are the eukaryotic derivatives-the spliceosome, spliceosomal introns, and non-LTR retroelements-which together make up approximately half of the human genome. On a practical level, the properties of group II introns have been exploited to develop group II intron-based biotechnological tools. WIREs RNA 2016, 7:341-355. doi: 10.1002/wrna.1339 For further resources related to this article, please visit the WIREs website.

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.

Evolution of group II introns.

Present in the genomes of bacteria and eukaryotic organelles, group II introns are an ancient class of ribozymes and retroelements that are believed to have been the ancestors of nuclear pre-mRNA introns. Despite long-standing speculation, there is limited understanding about the actual pathway by which group II introns evolved into eukaryotic introns. In this review, we focus on the evolution of group II introns themselves. We describe the different forms of group II introns known to exist in nature and then address how these forms may have evolved to give rise to spliceosomal introns and other genetic elements. Finally, we summarize the structural and biochemical parallels between group II introns and the spliceosome, including recent data that strongly support their hypothesized evolutionary relationship.

Novel RNA structural features of an alternatively splicing group II intron from Clostridium tetani.

Group II introns are ribozymes in bacterial and organellar genomes that function as self-splicing introns and as retroelements. Previously, we reported that the group II intron C.te.I1 of Clostridium tetani alternatively splices in vivo to produce five distinct coding mRNAs. Accurate fusion of upstream and downstream reading frames requires a shifted 5' splice site located 8 nt upstream of the usual 5' GUGYG motif. This site is specified by the ribozyme through an altered intron/exon-binding site 1 (IBS1-EBS1) pairing. Here we use mutagenesis and self-splicing assays to investigate in more detail the significance of the structural features of the C.te.I1 ribozyme. The shifted 5' splice site is shown to be affected by structures in addition to IBS1-EBS1, and unlike other group II introns, C.te.I1 appears to require a spacer between IBS1 and the GUGYG motif. In addition, the mechanism of 3' exon recognition is modified from the ancestral IIB mechanism to a IIA-like mechanism that appears to be longer than the typical single base-pair interaction and may extend up to 4 bp. The novel ribozyme properties that have evolved for C.te.I1 illustrate the plasticity of group II introns in adapting new structural and catalytic properties that can be utilized to affect gene expression.

Alternative splicing of a group II intron in a surface layer protein gene in Clostridium tetani.

Group II introns are ribozymes and retroelements found in bacteria, and are thought to have been the ancestors of nuclear pre-mRNA introns. Whereas nuclear introns undergo prolific alternative splicing in some species, group II introns are not known to carry out equivalent reactions. Here we report a group II intron in the human pathogen Clostridium tetani, which undergoes four alternative splicing reactions in vivo. Together with unspliced transcript, five mRNAs are produced, each encoding a distinct surface layer protein isoform. Correct fusion of exon reading frames requires a shifted 5' splice site located 8 nt upstream of the canonical boundary motif. The shifted junction is accomplished by an altered IBS1-EBS1 pairing between the intron and 5' exon. Growth of C. tetani under a variety of conditions did not result in large changes in alternative splicing levels, raising the possibility that alternative splicing is constitutive. This work demonstrates a novel type of gene organization and regulation in bacteria, and provides an additional parallel between group II and nuclear pre-mRNA introns.

A pipeline of programs for collecting and analyzing group II intron retroelement sequences from GenBank.

BACKGROUND: Accurate and complete identification of mobile elements is a challenging task in the current era of sequencing, given their large numbers and frequent truncations. Group II intron retroelements, which consist of a ribozyme and an intron-encoded protein (IEP), are usually identified in bacterial genomes through their IEP; however, the RNA component that defines the intron boundaries is often difficult to identify because of a lack of strong sequence conservation corresponding to the RNA structure. Compounding the problem of boundary definition is the fact that a majority of group II intron copies in bacteria are truncated. RESULTS: Here we present a pipeline of 11 programs that collect and analyze group II intron sequences from GenBank. The pipeline begins with a BLAST search of GenBank using a set of representative group II IEPs as queries. Subsequent steps download the corresponding genomic sequences and flanks, filter out non-group II introns, assign introns to phylogenetic subclasses, filter out incomplete and/or non-functional introns, and assign IEP sequences and RNA boundaries to the full-length introns. In the final step, the redundancy in the data set is reduced by grouping introns into sets of ≥95% identity, with one example sequence chosen to be the representative. CONCLUSIONS: These programs should be useful for comprehensive identification of group II introns in sequence databases as data continue to rapidly accumulate.

Structural variation and uniformity among tetraloop-receptor interactions and other loop-helix interactions in RNA crystal structures.

Tetraloop-receptor interactions are prevalent structural units in RNAs, and include the GAAA/11-nt and GNRA-minor groove interactions. In this study, we have compiled a set of 78 nonredundant loop-helix interactions from X-ray crystal structures, and examined them for the extent of their sequence and structural variation. Of the 78 interactions in the set, only four were classical GAAA/11-nt motifs, while over half (48) were GNRA-minor groove interactions. The GNRA-minor groove interactions were not a homogeneous set, but were divided into five subclasses. The most predominant subclass is characterized by two triple base pair interactions in the minor groove, flanked by two ribose zipper contacts. This geometry may be considered the "standard" GNRA-minor groove interaction, while the other four subclasses are alternative ways to form interfaces between a minor groove and tetraloop. The remaining 26 structures in the set of 78 have loops interacting with mostly idiosyncratic receptors. Among the entire set, a number of sequence-structure correlations can be identified, which may be used as initial hypotheses in predicting three-dimensional structures from primary sequences. Conversely, other sequence patterns are not predictive; for example, GAAA loop sequences and GG/CC receptors bind to each other with three distinct geometries. Finally, we observe an example of structural evolution in group II introns, in which loop-receptor motifs are substituted for each other while maintaining the larger three-dimensional geometry. Overall, the study gives a more complete view of RNA loop-helix interactions that exist in nature.

Database for bacterial group II introns.

The Database for Bacterial Group II Introns (http://webapps2.ucalgary.ca/~groupii/index.html#) provides a catalogue of full-length, non-redundant group II introns present in bacterial DNA sequences in GenBank. The website is divided into three sections. The first section provides general information on group II intron properties, structures and classification. The second and main section lists information for individual introns, including insertion sites, DNA sequences, intron-encoded protein sequences and RNA secondary structure models. The final section provides tools for identification and analysis of intron sequences. These include a step-by-step guide to identify introns in genomic sequences, a local BLAST tool to identify closest intron relatives to a query sequence, and a boundary-finding tool that predicts 5' and 3' intron-exon junctions in an input DNA sequence. Finally, selected intron data can be downloaded in FASTA format. It is hoped that this database will be a useful resource not only to group II intron and RNA researchers, but also to microbiologists who encounter these unexpected introns in genomic sequences.

A reverse transcriptase-related protein mediates phage resistance and polymerizes untemplated DNA in vitro.

Reverse transcriptases (RTs) are RNA-dependent DNA polymerases that usually function in the replication of selfish DNAs such as retrotransposons and retroviruses. Here, we have biochemically characterized a RT-related protein, AbiK, which is required for abortive phage infection in the Gram-positive bacterium Lactococcus lactis. In vitro, AbiK does not exhibit the properties expected for an RT, but polymerizes long DNAs of 'random' sequence, analogous to a terminal transferase. Moreover, the polymerized DNAs appear to be covalently attached to the AbiK protein, presumably because an amino acid serves as a primer. Mutagenesis experiments indicate that the polymerase activity resides in the RT motifs and is essential for phage resistance in vivo. These results establish a novel biochemical property and a non-replicative biological role for a polymerase.

Group II introns: mobile ribozymes that invade DNA.

Group II introns are mobile ribozymes that self-splice from precursor RNAs to yield excised intron lariat RNAs, which then invade new genomic DNA sites by reverse splicing. The introns encode a reverse transcriptase that stabilizes the catalytically active RNA structure for forward and reverse splicing, and afterwards converts the integrated intron RNA back into DNA. The characteristics of group II introns suggest that they or their close relatives were evolutionary ancestors of spliceosomal introns, the spliceosome, and retrotransposons in eukaryotes. Further, their ribozyme-based DNA integration mechanism enabled the development of group II introns into gene targeting vectors ("targetrons"), which have the unique feature of readily programmable DNA target specificity.

A broadscale phylogenetic analysis of group II intron RNAs and intron-encoded reverse transcriptases.

Group II introns are self-splicing RNAs that are frequently assumed to be the ancestors of spliceosomal introns. They are widely distributed in bacteria and are also found in organelles of plants, fungi, and protists. In this study, we present a broadscale phylogenetic analysis of group II introns using sequence data from both the conserved RNA structure and the intron-encoded reverse transcriptase (RT). Two similar phylogenies are estimated for the RT open reading frame (ORF), based on either amino acid or nucleotide sequence, whereas one phylogeny is produced for the RNA. In making these estimates, we confronted nearly all the classic challenges to phylogenetic inference, including positional saturation, base composition heterogeneity, short internodes with low support, and sensitivity to taxon sampling. Although the major lineages are well-defined, robust resolution of topology is not possible between these lineages. The approximately unbiased (AU) and Shimodaira-Hasegawa topology tests indicated that the RT ORF and RNA ribozyme data sets are in significant conflict under a variety of models, revealing the possibility of imperfect coevolution between group II introns and their intron-encoded ORFs. The high level of sequence divergence, large timescale, and limited number of alignable characters in our study are representative of many RTs and group I introns, and our results suggest that phylogenetic analyses of any of these sequences could suffer from the same sources of error and instability identified in this study.

A diversity of uncharacterized reverse transcriptases in bacteria.

Retroelements are usually considered to be eukaryotic elements because of the large number and variety in eukaryotic genomes. By comparison, reverse transcriptases (RTs) are rare in bacteria, with only three characterized classes: retrons, group II introns and diversity-generating retroelements (DGRs). Here, we present the results of a bioinformatic survey that aims to define the landscape of RTs across eubacterial, archaeal and phage genomes. We identify and categorize 1021 RTs, of which the majority are group II introns (73%). Surprisingly, a plethora of novel RTs are found that do not belong to characterized classes. The RTs have 11 domain architectures and are classified into 20 groupings based on sequence similarity, phylogenetic analyses and open reading frame domain structures. Interestingly, group II introns are the only bacterial RTs to exhibit clear evidence for independent mobility, while five other groups have putative functions in defense against phage infection or promotion of phage infection. These examples suggest that additional beneficial functions will be discovered among uncharacterized RTs. The study lays the groundwork for experimental characterization of these highly diverse sequences and has implications for the evolution of retroelements.

Group II introns in eubacteria and archaea: ORF-less introns and new varieties.

Group II introns are a major class of ribozymes found in bacteria, mitochondria, and plastids. Many introns contain reverse transcriptase open reading frames (ORFs) that confer mobility to the introns and allow them to persist as selfish DNAs. Here, we report an updated compilation of group II introns in Eubacteria and Archaea comprising 234 introns. One new phylogenetic class is identified, as well as several specialized lineages. In addition, we undertake a detailed search for ORF-less group II introns in bacterial genomes in order to find undiscovered introns that either entirely lack an ORF or encode a novel ORF. Unlike organellar group II introns, we find only a handful of ORF-less introns in bacteria, suggesting that if a substantial number exist, they must be divergent from known introns. Together, these results highlight the retroelement character of bacterial group II introns, and suggest that their long-term survival is dependent upon retromobility.

A three-dimensional model of a group II intron RNA and its interaction with the intron-encoded reverse transcriptase.

Group II introns are self-splicing ribozymes believed to be the ancestors of spliceosomal introns. Many group II introns encode reverse transcriptases that promote both RNA splicing and intron mobility to new genomic sites. Here we used a circular permutation and crosslinking method to establish 16 intramolecular distance relationships within the mobile Lactococcus lactis Ll.LtrB-DeltaORF intron. Using these new constraints together with 13 established tertiary interactions and eight published crosslinks, we modeled a complete three-dimensional structure of the intron. We also used the circular permutation strategy to map RNA-protein interaction sites through fluorescence quenching and crosslinking assays. Our model provides a comprehensive structural framework for understanding the function of group II ribozymes, their natural structural variations, and the mechanisms by which the intron-encoded protein promotes RNA splicing and intron mobility. The model also suggests an arrangement of active site elements that may be conserved in the spliceosome.

Insertion of group II intron retroelements after intrinsic transcriptional terminators.

Mobile DNAs use many mechanisms to minimize damage to their hosts. Here we show that a subclass of group II introns avoids host damage by inserting directly after transcriptional terminator motifs in bacterial genomes (stem-loops followed by Ts). This property contrasts with the site-specific behavior of most group II introns, which insert into homing site sequences. Reconstituted ribonucleo protein particles of the Bacillus halodurans intron B.h.I1 are shown to reverse-splice into DNA targets in vitro but require the DNA to be single-stranded and fold into a stem-loop analogous to the RNA structure that forms during transcription termination. Recognition of this DNA stem-loop motif accounts for in vivo target specificity. Insertion after terminators is a previously unrecognized strategy for a selfish DNA because it prevents interruption of coding sequences and restricts expression of the mobile DNA after integration.

Self-splicing of a group IIC intron: 5' exon recognition and alternative 5' splicing events implicate the stem-loop motif of a transcriptional terminator.

Bacterial IIC introns are a newly recognized subclass of group II introns whose ribozyme properties have not been characterized in detail. IIC introns are typically located downstream of transcriptional terminator motifs (inverted repeat followed by T's) or other inverted repeats in bacterial genomes. Here we have characterized the self-splicing activity of a IIC intron, B.h.I1, from Bacillus halodurans. B.h.I1 self-splices in vitro through hydrolysis to produce linear intron, but interestingly, additional unexpected products were formed that were highly dependent on ionic conditions. These products were determined to represent alternative splicing events at the 5' junction and cleavages throughout the RNA transcript. The alternative splicing and cleavage events occurred at cryptic splice sites containing stem-loop and IBS1 motifs, suggesting that the 5' exon is recognized by both elements. These results provide the first example of a group II intron that uses 5' splice sites nonadjacent to the ribozyme structure. Furthermore, the data suggest that IIC introns differ from IIA and IIB introns with respect to 5' exon definition, and that the terminator stem-loop substitutes in part for the missing IBS2-EBS2 (intron and exon binding sites 2) interaction.

Origin and evolution of the chloroplast trnK (matK) intron: a model for evolution of group II intron RNA structures.

The trnK intron of plants encodes the matK open reading frame (ORF), which has been used extensively as a phylogenetic marker for classification of plants. Here we examined the evolution of the trnK intron itself as a model for group II intron evolution in plants. Representative trnK intron sequences were compiled from species spanning algae to angiosperms, and four introns were newly sequenced. Phylogenetic analyses showed that the matK ORFs belong to the ML (mitochondrial-like) subclass of group II intron ORFs, indicating that they were derived from a mobile group II intron of the class. RNA structures of the introns were folded and analyzed, which revealed progressive RNA structural deviations and degenerations throughout plant evolution. The data support a model in which plant organellar group II introns were derived from bacterial-like introns that had "standard" RNA structures and were competent for self-splicing and mobility and that subsequently the ribozyme structures degenerated to ultimately become dependent upon host-splicing factors. We propose that the patterns of RNA structure evolution seen for the trnK intron will apply to the other group II introns in plants.

Mobile group II introns.

Mobile group II introns, found in bacterial and organellar genomes, are both catalytic RNAs and retrotransposable elements. They use an extraordinary mobility mechanism in which the excised intron RNA reverse splices directly into a DNA target site and is then reverse transcribed by the intron-encoded protein. After DNA insertion, the introns remove themselves by protein-assisted, autocatalytic RNA splicing, thereby minimizing host damage. Here we discuss the experimental basis for our current understanding of group II intron mobility mechanisms, beginning with genetic observations in yeast mitochondria, and culminating with a detailed understanding of molecular mechanisms shared by organellar and bacterial group II introns. We also discuss recently discovered links between group II intron mobility and DNA replication, new insights into group II intron evolution arising from bacterial genome sequencing, and the evolutionary relationship between group II introns and both eukaryotic spliceosomal introns and non-LTR-retrotransposons. Finally, we describe the development of mobile group II introns into gene-targeting vectors, "targetrons," which have programmable target specificity.

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

Bordetella bacteriophages generate diversity in a gene that specifies host tropism. This microevolutionary adaptation is produced by a genetic element that combines the basic retroelement life cycle of transcription, reverse transcription and integration with site-directed, adenine-specific mutagenesis. Central to this process is a reverse transcriptase-mediated exchange between two repeats; one serving as a donor template (TR) and the other as a recipient of variable sequence information (VR). Here we describe the genetic basis for diversity generation. The directionality of information transfer is determined by a 21-base-pair sequence present at the 3' end of VR. On the basis of patterns of marker transfer in response to variant selective pressures, we propose that a TR reverse transcript is mutagenized, integrated into VR as a single non-coding strand, and then partially converted to the parental VR sequence. This allows the diversity-generating system to minimize variability to the subset of bases under selection. Using the Bordetella phage cassette as a signature, we have identified numerous related elements in diverse bacteria. These elements constitute a new family of retroelements with the potential to confer selective advantages to their host genomes.