Mobile Bacterial Group II Introns at the Crux of Eukaryotic Evolution

WHO elements - A new category of selfish genetic elements at the borderline between homing elements and transposable elements

Homing genetic elements are a form of selfish DNA that inserts into a specific target site in the genome and spreads through the population by a process of biased inheritance. Two well-known types of homing element, called inteins and homing introns, were discovered decades ago. In this review we describe WHO elements, a newly discovered type of homing element that constitutes a distinct third category but is rare, having been found only in a few yeast species so far. WHO elements are inferred to spread using the same molecular homing mechanism as inteins and introns: they encode a site-specific endonuclease that cleaves the genome at the target site, making a DNA break that is subsequently repaired by copying the element. For most WHO elements, the target site is in the glycolytic gene FBA1. WHO elements differ from inteins and homing introns in two fundamental ways: they do not interrupt their host gene (FBA1), and they occur in clusters. The clusters were formed by successive integrations of different WHO elements into the FBA1 locus, the result of an 'arms race' between the endonuclease and its target site. We also describe one family of WHO elements (WHO10) that is no longer specifically associated with the FBA1 locus and instead appears to have become transposable, inserting at random genomic sites in Torulaspora globosa with up to 26 copies per strain. The WHO family of elements is therefore at the borderline between homing genetic elements and transposable elements.

Intron lariat spliceosomes convert lariats to true circles: implications for intron transposition

Rare, full-length circular intron RNAs distinct from lariats have been reported in several species, but their biogenesis is not understood. We envisioned and tested a hypothesis for their formation using Saccharomyces cerevisiae, documenting full-length and novel processed circular RNAs from multiple introns. Evidence implicates a previously undescribed catalytic activity of the intron lariat spliceosome (ILS) in which the 3'-OH of the lariat tail (with optional trimming and adenylation by the nuclear 3' processing machinery) attacks the branch, joining the intron 3' end to the 5' splice site in a 3'-5' linked circle. Human U2 and U12 spliceosomes produce analogous full-length and processed circles. Postsplicing catalytic activity of the spliceosome may promote intron transposition during eukaryotic genome evolution.

Mitochondrial RNA Helicases: Key Players in the Regulation of Plant Organellar RNA Splicing and Gene Expression

Mitochondrial genomes of land plants are large and exhibit a complex mode of gene organization and expression, particularly at the post-transcriptional level. The primary organellar transcripts in plants undergo extensive maturation steps, including endo- and/or exo-nucleolytic cleavage, RNA-base modifications (mostly C-to-U deaminations) and both 'cis'- and 'trans'-splicing events. These essential processing steps rely on the activities of a large set of nuclear-encoded factors. RNA helicases serve as key players in RNA metabolism, participating in the regulation of transcription, mRNA processing and translation. They unwind RNA secondary structures and facilitate the formation of ribonucleoprotein complexes crucial for various stages of gene expression. Furthermore, RNA helicases are involved in RNA metabolism by modulating pre-mRNA maturation, transport and degradation processes. These enzymes are, therefore, pivotal in RNA quality-control mechanisms, ensuring the fidelity and efficiency of RNA processing and turnover in plant mitochondria. This review summarizes the significant roles played by helicases in regulating the highly dynamic processes of mitochondrial transcription, RNA processing and translation in plants. We further discuss recent advancements in understanding how dysregulation of mitochondrial RNA helicases affects the splicing of organellar genes, leading to respiratory dysfunctions, and consequently, altered growth, development and physiology of land plants.

Untacking small RNA profiling and RNA fragment footprinting: Approaches and challenges in library construction

Small RNAs (sRNAs) with sizes ranging from 15 to 50 nucleotides (nt) are critical regulators of gene expression control. Prior studies have shown that sRNAs are involved in a broad range of biological processes, such as organ development, tumorigenesis, and epigenomic regulation; however, emerging evidence unveils a hidden layer of diversity and complexity of endogenously encoded sRNAs profile in eukaryotic organisms, including novel types of sRNAs and the previously unknown post-transcriptional RNA modifications. This underscores the importance for accurate, unbiased detection of sRNAs in various cellular contexts. A multitude of high-throughput methods based on next-generation sequencing (NGS) are developed to decipher the sRNA expression and their modifications. Nonetheless, distinct from mRNA sequencing, the data from sRNA sequencing suffer frequent inconsistencies and high variations emanating from the adapter contaminations and RNA modifications, which overall skew the sRNA libraries. Here, we summarize the sRNA-sequencing approaches, and discuss the considerations and challenges for the strategies and methods of sRNA library construction. The pros and cons of sRNA sequencing have significant implications for implementing RNA fragment footprinting approaches, including CLIP-seq and Ribo-seq. We envision that this review can inspire novel improvements in small RNA sequencing and RNA fragment footprinting in future. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA Processing > Processing of Small RNAs Regulatory RNAs/RNAi/Riboswitches > Biogenesis of Effector Small RNAs.

Autonomous transposons tune their sequences to ensure somatic suppression

Transposable elements (TEs) are a major constituent of human genes, occupying approximately half of the intronic space. During pre-messenger RNA synthesis, intronic TEs are transcribed along with their host genes but rarely contribute to the final mRNA product because they are spliced out together with the intron and rapidly degraded. Paradoxically, TEs are an abundant source of RNA-processing signals through which they can create new introns1, and also functional2 or non-functional chimeric transcripts3. The rarity of these events implies the existence of a resilient splicing code that is able to suppress TE exonization without compromising host pre-mRNA processing. Here we show that SAFB proteins protect genome integrity by preventing retrotransposition of L1 elements while maintaining splicing integrity, via prevention of the exonization of previously integrated TEs. This unique dual role is possible because of L1's conserved adenosine-rich coding sequences that are bound by SAFB proteins. The suppressive activity of SAFB extends to tissue-specific, giant protein-coding cassette exons, nested genes and Tigger DNA transposons. Moreover, SAFB also suppresses LTR/ERV elements in species in which they are still active, such as mice and flies. A significant subset of splicing events suppressed by SAFB in somatic cells are activated in the testis, coinciding with low SAFB expression in postmeiotic spermatids. Reminiscent of the division of labour between innate and adaptive immune systems that fight external pathogens, our results uncover SAFB proteins as an RNA-based, pattern-guided, non-adaptive defence system against TEs in the soma, complementing the RNA-based, adaptive Piwi-interacting RNA pathway of the germline.

A Novel Thermostable and Processive Reverse Transcriptase from a Group II Intron of Anoxybacillus flavithermus

Reverse transcriptases (RTs) are a family of enzymes that synthesize DNA using an RNA template and are involved in retrovirus propagation and telomere lengthening. In vitro, RTs are widely applied in various methods, including RNA-seq, RT-PCR, and RT-LAMP. Thermostable RTs from bacterial group II introns are promising tools for biotechnology due to their higher thermostability, fidelity, and processivity compared to commonly used M-MuLV RT and its mutants. However, the diversity of group II intron-encoded RTs is still understudied. In this work, we biochemically characterized a novel RT from a thermophilic bacterium, Anoxybacillus flavithermus, which was isolated from a hot spring in New Zealand and has an optimal growth temperature of around 60 °C. The cloned RT, named Afl RT, retained approximately 40% of the specific activity after a 45 min incubation at 50 °C. The optimal pH was 8.5, the optimal temperature was between 45 and 50 °C, and Mn2+ ions were found to be an optimal cofactor. The processivity analysis with MS2 phage gRNA (3569 b) demonstrated that Afl RT elongated fully up to 36% of the template molecules. In reverse transcription and RT-qLAMP, the enzyme allowed up to 10 copies of MS2 phage genomic RNA to be detected per reaction. Thus, Afl RT holds great potential for a variety of practical applications that require the use of thermostable and processive RTs.

To Be Mobile or Not: The Variety of Reverse Transcriptases and Their Recruitment by Host Genomes

Reverse transcriptases (RT), or RNA-dependent DNA polymerases, are unorthodox enzymes that originally added a new angle to the conventional view of the unidirectional flow of genetic information in the cell from DNA to RNA to protein. First discovered in vertebrate retroviruses, RTs were since re-discovered in most eukaryotes, bacteria, and archaea, spanning essentially all domains of life. For retroviruses, RTs provide the ability to copy the RNA genome into DNA for subsequent incorporation into the host genome, which is essential for their replication and survival. In cellular organisms, most RT sequences originate from retrotransposons, the type of self-replicating genetic elements that rely on reverse transcription to copy and paste their sequences into new genomic locations. Some retroelements, however, can undergo domestication, eventually becoming a valuable addition to the overall repertoire of cellular enzymes. They can be beneficial yet accessory, like the diversity-generating elements, or even essential, like the telomerase reverse transcriptases. Nowadays, ever-increasing numbers of domesticated RT-carrying genetic elements are being discovered. It may be argued that domesticated RTs and reverse transcription in general is more widespread in cellular organisms than previously thought, and that many important cellular functions, such as chromosome end maintenance, may evolve from an originally selfish process of converting RNA into DNA.

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

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

Structural RNA components supervise the sequential DNA cleavage in R2 retrotransposon

Retroelements are the widespread jumping elements considered as major drivers for genome evolution, which can also be repurposed as gene-editing tools. Here, we determine the cryo-EM structures of eukaryotic R2 retrotransposon with ribosomal DNA target and regulatory RNAs. Combined with biochemical and sequencing analysis, we reveal two essential DNA regions, Drr and Dcr, required for recognition and cleavage. The association of 3' regulatory RNA with R2 protein accelerates the first-strand cleavage, blocks the second-strand cleavage, and initiates the reverse transcription starting from the 3'-tail. Removing 3' regulatory RNA by reverse transcription allows the association of 5' regulatory RNA and initiates the second-strand cleavage. Taken together, our work explains the DNA recognition and RNA supervised sequential retrotransposition mechanisms by R2 machinery, providing insights into the retrotransposon and application reprogramming.

MSP1 encodes an essential RNA-binding pentatricopeptide repeat factor required for nad1 maturation and complex I biogenesis in Arabidopsis mitochondria

Mitochondrial biogenesis relies on nuclearly encoded factors, which regulate the expression of the organellar-encoded genes. Pentatricopeptide repeat (PPR) proteins constitute a major gene family in angiosperms that are pivotal in many aspects of mitochondrial (mt)RNA metabolism (e.g. trimming, splicing, or stability). Here, we report the analysis of MITOCHONDRIA STABILITY/PROCESSING PPR FACTOR1 (MSP1, At4g20090), a canonical PPR protein that is necessary for mitochondrial functions and embryo development. Loss-of-function allele of MSP1 leads to seed abortion. Here, we employed an embryo-rescue method for the molecular characterization of msp1 mutants. Our analyses reveal that msp1 embryogenesis fails to proceed beyond the heart/torpedo stage as a consequence of a nad1 pre-RNA processing defect, resulting in the loss of respiratory complex I activity. Functional complementation confirmed that msp1 phenotypes result from a disruption of the MSP1 gene. In Arabidopsis, the maturation of nad1 involves the processing of three RNA fragments, nad1.1, nad1.2, and nad1.3. Based on biochemical analyses and mtRNA profiles of wild-type and msp1 plants, we concluded that MSP1 facilitates the generation of the 3' terminus of nad1.1 transcript, a prerequisite for nad1 exons a-b splicing. Our data substantiate the importance of mtRNA metabolism for the biogenesis of the respiratory system during early plant life.

Categorizing 161 plant (streptophyte) mitochondrial group II introns into 29 families of related paralogues finds only limited links between intron mobility and intron-borne maturases

Group II introns are common in the two endosymbiotic organelle genomes of the plant lineage. Chloroplasts harbor 22 positionally conserved group II introns whereas their occurrence in land plant (embryophyte) mitogenomes is highly variable and specific for the seven major clades: liverworts, mosses, hornworts, lycophytes, ferns, gymnosperms and flowering plants. Each plant group features "signature selections" of ca. 20-30 paralogues from a superset of altogether 105 group II introns meantime identified in embryophyte mtDNAs, suggesting massive intron gains and losses along the backbone of plant phylogeny. We report on systematically categorizing plant mitochondrial group II introns into "families", comprising evidently related paralogues at different insertion sites, which may even be more similar than their respective orthologues in phylogenetically distant taxa. Including streptophyte (charophyte) algae extends our sampling to 161 and we sort 104 streptophyte mitochondrial group II introns into 25 core families of related paralogues evidently arising from retrotransposition events. Adding to discoveries of only recently created intron paralogues, hypermobile introns and twintrons, our survey led to further discoveries including previously overlooked "fossil" introns in spacer regions or e.g., in the rps8 pseudogene of lycophytes. Initially excluding intron-borne maturase sequences for family categorization, we added an independent analysis of maturase phylogenies and find a surprising incongruence between intron mobility and the presence of intron-borne maturases. Intriguingly, however, we find that several examples of nuclear splicing factors meantime characterized simultaneously facilitate splicing of independent paralogues now placed into the same intron families. Altogether this suggests that plant group II intron mobility, in contrast to their bacterial counterparts, is not intimately linked to intron-encoded maturases.

Reverse Transcriptases: From Discovery and Applications to Xenobiology

Reverse transcriptases are DNA polymerases that can use RNA as a template for DNA synthesis. They thus catalyze the reverse of transcription. Although discovered in 1970, reverse transcriptases are still of great interest and are constantly being further developed for numerous modern research approaches. They are frequently used in biotechnological and molecular diagnostic applications. In this review, we describe the discovery of these fascinating enzymes and summarize research results and applications ranging from molecular cloning, direct virus detection, and modern sequencing methods to xenobiology.

Reverse transcriptase meets DNA, again: Possible roles for transposable elements in host DNA repair

Retrotransposons are selfish genetic elements that encode an enzyme, reverse transcriptase (RT), which converts the element-encoded RNA into DNA prior to or during genomic integration. New studies provide compelling evidence that a bacterial group II intron-like RT has adapted enzymatic activities associated with RTs to function in host DNA repair.

Group II intron-like reverse transcriptases function in double-strand break repair

Bacteria encode reverse transcriptases (RTs) of unknown function that are closely related to group II intron-encoded RTs. We found that a Pseudomonas aeruginosa group II intron-like RT (G2L4 RT) with YIDD instead of YADD at its active site functions in DNA repair in its native host and when expressed in Escherichia coli. G2L4 RT has biochemical activities strikingly similar to those of human DNA repair polymerase θ and uses them for translesion DNA synthesis and double-strand break repair (DSBR) via microhomology-mediated end-joining (MMEJ). We also found that a group II intron RT can function similarly in DNA repair, with reciprocal active-site substitutions showing isoleucine favors MMEJ and alanine favors primer extension in both enzymes. These DNA repair functions utilize conserved structural features of non-LTR-retroelement RTs, including human LINE-1 and other eukaryotic non-LTR-retrotransposon RTs, suggesting such enzymes may have inherent ability to function in DSBR in a wide range of organisms.

Dynamic plastid and mitochondrial genomes in Chaetopeltidales (Chlorophyceae) and characterization of a new chlorophyte taxon

PREMISE

Chaetopeltidales is a poorly characterized order in the Chlorophyceae, with only two plastid and no mitochondrial genomes published. Here we describe a new taxon in Chaetopeltidales, Gormaniella terricola gen. et sp. nov. and characterize both of its organellar genomes.

METHODS

Gormaniella terricola was inadvertently isolated from a surface-sterilized hornwort thallus. Light microscopy was used to characterize its vegetative morphology. Organellar genomes were assembled, annotated, and analyzed using a variety of software packages.

RESULTS

The mitochondrial genome (66,927 bp) represents the first complete mitochondrial genome published for Chaetopeltidales. The chloroplast genome, measuring 428,981 bp, is one of the largest plastid genomes published to date and shares this large size and an incredible number of short, dispersed repeats with the other sequenced chloroplast genomes in Chaetopeltidales. Despite these shared features, the chloroplast genomes of Chaetopeltidales appear to be highly rearranged when compared to one another, with numerous inversions, translocations, and duplications, suggesting a particularly dynamic chloroplast genome. Both the chloroplast and mitochondrial genomes of G. terricola contain a number of mobile group I and group II introns, which appear to have invaded separately. Three of the introns within the mitochondrial genome encode homing endonucleases that are phylogenetically nested within those found in fungi, rather than algae, suggesting a possible case of horizontal gene transfer.

CONCLUSIONS

These results help to shed light on a poorly understood group of algae and their unusual organellar genomes, raising additional questions about the unique patterns of genome evolution within Chaetopeltidales.

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.

The spread of the first introns in proto-eukaryotic paralogs

Spliceosomal introns are a unique feature of eukaryotic genes. Previous studies have established that many introns were present in the protein-coding genes of the last eukaryotic common ancestor (LECA). Intron positions shared between genes that duplicated before LECA could in principle provide insight into the emergence of the first introns. In this study we use ancestral intron position reconstructions in two large sets of duplicated families to systematically identify these ancient paralogous intron positions. We found that 20-35% of introns inferred to have been present in LECA were shared between paralogs. These shared introns, which likely preceded ancient duplications, were wide spread across different functions, with the notable exception of nuclear transport. Since we observed a clear signal of pervasive intron loss prior to LECA, it is likely that substantially more introns were shared at the time of duplication than we can detect in LECA. The large extent of shared introns indicates an early origin of introns during eukaryogenesis and suggests an early origin of a nuclear structure, before most of the other complex eukaryotic features were established.

Does the Pachytene Checkpoint, a Feature of Meiosis, Filter Out Mistakes in Double-Strand DNA Break Repair and as a side-Effect Strongly Promote Adaptive Speciation?

This essay aims to explain two biological puzzles: why eukaryotic transcription units are composed of short segments of coding DNA interspersed with long stretches of non-coding (intron) DNA, and the near ubiquity of sexual reproduction. As is well known, alternative splicing of its coding sequences enables one transcription unit to produce multiple variants of each encoded protein. Additionally, padding transcription units with non-coding DNA (often many thousands of base pairs long) provides a readily evolvable way to set how soon in a cell cycle the various mRNAs will begin being expressed and the total amount of mRNA that each transcription unit can make during a cell cycle. This regulation complements control via the transcriptional promoter and facilitates the creation of complex eukaryotic cell types, tissues, and organisms. However, it also makes eukaryotes exceedingly vulnerable to double-strand DNA breaks, which end-joining break repair pathways can repair incorrectly. Transcription units cover such a large fraction of the genome that any mis-repair producing a reorganized chromosome has a high probability of destroying a gene. During meiosis, the synaptonemal complex aligns homologous chromosome pairs and the pachytene checkpoint detects, selectively arrests, and in many organisms actively destroys gamete-producing cells with chromosomes that cannot adequately synapse; this creates a filter favoring transmission to the next generation of chromosomes that retain the parental organization, while selectively culling those with interrupted transcription units. This same meiotic checkpoint, reacting to accidental chromosomal reorganizations inflicted by error-prone break repair, can, as a side effect, provide a mechanism for the formation of new species in sympatry. It has been a long-standing puzzle how something as seemingly maladaptive as hybrid sterility between such new species can arise. I suggest that this paradox is resolved by understanding the adaptive importance of the pachytene checkpoint, as outlined above.

Recruitment of Mobile Genetic Elements for Diverse Cellular Functions in Prokaryotes

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

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.

Identification of Group II Intron RmInt1 Binding Sites in a Bacterial Genome

RmInt1 is a group II intron encoding a reverse transcriptase protein (IEP) lacking the C-terminal endonuclease domain. RmInt1 is an efficient mobile retroelement that predominantly reverse splices into the transient single-stranded DNA at the template for lagging strand DNA synthesis during host replication, a process facilitated by the interaction of the RmInt1 IEP with DnaN at the replication fork. It has been suggested that group II intron ribonucleoprotein particles bind DNA nonspecifically, and then scan for their correct target site. In this study, we investigated RmInt1 binding sites throughout the Sinorhizobium meliloti genome, by chromatin-immunoprecipitation coupled with next-generation sequencing. We found that RmInt1 binding sites cluster around the bidirectional replication origin of each of the three replicons comprising the S. meliloti genome. Our results provide new evidence linking group II intron mobility to host DNA replication.

Natural circularly permuted group II introns in bacteria produce RNA circles

Group II self-splicing introns are large structured RNAs that remove themselves from transcripts while simultaneously sealing the resulting gaps. Some representatives can subsequently reverse splice into DNA, accounting for their pervasive distribution in bacteria. The catalytically active tertiary structure of each group II intron is assembled from six domains that are arranged in a conserved order. Here, we report structural isomers of group II introns, called CP group II ribozymes, wherein the characteristic order of domains has been altered. Domains five and six, which normally reside at the 3' end of group II introns, instead occupy the 5' end to form circularly permuted variants. These unusual group II intron derivatives are catalytically active and generate large linear branched and small circular RNAs, reaction products that are markedly different from those generated by canonical group II introns. The biological role of CP group II ribozymes is currently unknown.

Targeted insertion of large genetic payloads using cas directed LINE-1 reverse transcriptase

A difficult genome editing goal is the site-specific insertion of large genetic constructs. Here we describe the GENEWRITE system, where site-specific targetable activity of Cas endonucleases is coupled with the reverse transcriptase activity of the ORF2p protein of the human retrotransposon LINE-1. This is accomplished by providing two RNAs: a guide RNA targeting Cas endonuclease activity and an appropriately designed payload RNA encoding the desired insertion. Using E. coli as a simple platform for development and deployment, we show that with proper payload design and co-expression of helper proteins, GENEWRITE can enable insertion of large genetic payloads to precise locations, although with off-target effects, using the described approach. Based upon these results, we describe a potential strategy for implementation of GENEWRITE in more complex systems.

TGIRT-seq Protocol for the Comprehensive Profiling of Coding and Non-coding RNA Biotypes in Cellular, Extracellular Vesicle, and Plasma RNAs

High-throughput RNA sequencing (RNA-seq) has extraordinarily advanced our understanding of gene expression and disease etiology, and is a powerful tool for the identification of biomarkers in a wide range of organisms. However, most RNA-seq methods rely on retroviral reverse transcriptases (RTs), enzymes that have inherently low fidelity and processivity, to convert RNAs into cDNAs for sequencing. Here, we describe an RNA-seq protocol using Thermostable Group II Intron Reverse Transcriptases (TGIRTs), which have high fidelity, processivity, and strand-displacement activity, as well as a proficient template-switching activity that enables efficient and seamless RNA-seq adapter addition. By combining these activities, TGIRT-seq enables the simultaneous profiling of all RNA biotypes from small amounts of starting material, with superior RNA-seq metrics, and unprecedented ability to sequence structured RNAs. The TGIRT-seq protocol for Illumina sequencing consists of three steps: (i) addition of a 3' RNA-seq adapter, coupled to the initiation of cDNA synthesis at the 3' end of a target RNA, via template switching from a synthetic adapter RNA/DNA starter duplex; (ii) addition of a 5' RNA-seq adapter, by using thermostable 5' App DNA/RNA ligase to ligate an adapter oligonucleotide to the 3' end of the completed cDNA; (iii) minimal PCR amplification, to add capture sites and indices for Illumina sequencing. TGIRT-seq for the Illumina sequencing platform has been used for comprehensive profiling of coding and non-coding RNAs in ribodepleted, chemically fragmented cellular RNAs, and for the analysis of intact (non-chemically fragmented) cellular, extracellular vesicle (EV), and plasma RNAs, where it yields continuous full-length end-to-end sequences of structured small non-coding RNAs (sncRNAs), including tRNAs, snoRNAs, snRNAs, pre-miRNAs, and full-length excised linear intron (FLEXI) RNAs. Graphic abstract: Figure 1.Overview of the TGIRT-seq protocol for Illumina sequencing.Major steps are: (1) Template switching from a synthetic R2 RNA/R2R DNA starter duplex with a 1-nt 3' DNA overhang (a mixture of A, C, G, and T residues, denoted N) that base pairs to the 3' nucleotide of a target RNA, and upon initiating reverse transcription by adding dNTPs, seamlessly links an R2R adapter to the 5' end of the resulting cDNA; (2) Ligation of an R1R adapter to the 3' end of the completed cDNA; and (3) Minimal PCR amplification with primers that add Illumina capture sites (P5 and P7) and barcode sequences (indices 5 and 7). The index 7 barcode is required, while the index 5 barcode is optional, to provide unique dual indices (UDIs).

A new RNA-DNA interaction required for integration of group II intron retrotransposons into DNA targets

Mobile group II introns are site-specific retrotransposable elements abundant in bacterial and organellar genomes. They are composed of a large and highly structured ribozyme and an intron-encoded reverse transcriptase that binds tightly to its intron to yield a ribonucleoprotein (RNP) particle. During the first stage of the mobility pathway, the intron RNA catalyses its own insertion directly into the DNA target site. Recognition of the proper target rests primarily on multiple base-pairing interactions between the intron RNA and the target DNA, while the protein makes contacts with only a few target positions by yet-unidentified mechanisms. Using a combination of comparative sequence analyses and in vivo mobility assays we demonstrate the existence of a new base-pairing interaction named EBS2a–IBS2a between the intron RNA and its DNA target site. This pairing adopts a Watson–Crick geometry and is essential for intron mobility, most probably by driving unwinding of the DNA duplex. Importantly, formation of EBS2a–IBS2a also requires the reverse transcriptase enzyme which stabilizes the pairing in a non-sequence-specific manner. In addition to bringing to light a new structural device that allows subgroup IIB1 and IIB2 introns to invade their targets with high efficiency and specificity our work has important implications for the biotechnological applications of group II introns in bacterial gene targeting.

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.

Low-bias ncRNA libraries using ordered two-template relay: Serial template jumping by a modified retroelement reverse transcriptase

Retrotransposons are noninfectious, mobile genetic elements that proliferate in host genomes via an RNA intermediate that is copied into DNA by a reverse transcriptase (RT) enzyme. RTs are important for biotechnological applications involving information capture from RNA since RNA is first converted into complementary DNA for detection or sequencing. Here, we biochemically characterized RTs from two retroelements and uncovered several activities that allowed us to design a streamlined, efficient workflow for determining the inventory of RNA sequences in processed RNA pools. The unique properties of nonretroviral RT activities obviate many technical issues associated with current methods of RNA sequence analysis, with wide applications in research, biotechnology, and diagnostics.

Selfish, non-long terminal repeat (non-LTR) retroelements and mobile group II introns encode reverse transcriptases (RTs) that can initiate DNA synthesis without substantial base pairing of primer and template. Biochemical characterization of these enzymes has been limited by recombinant expression challenges, hampering understanding of their properties and the possible exploitation of their properties for research and biotechnology. We investigated the activities of representative RTs using a modified non-LTR RT from Bombyx mori and a group II intron RT from Eubacterium rectale. Only the non-LTR RT supported robust and serial template jumping, producing one complementary DNA (cDNA) from several templates each copied end to end. We also discovered an unexpected terminal deoxynucleotidyl transferase activity of the RTs that adds nucleotide(s) of choice to 3′ ends of single- and/or double-stranded RNA or DNA. Combining these two types of activity with additional insights about nontemplated nucleotide additions to duplexed cDNA product, we developed a streamlined protocol for fusion of next-generation sequencing adaptors to both cDNA ends in a single RT reaction. When benchmarked using a reference pool of microRNAs (miRNAs), library production by Ordered Two-Template Relay (OTTR) using recombinant non-LTR retroelement RT outperformed all commercially available kits and rivaled the low bias of technically demanding home-brew protocols. We applied OTTR to inventory RNAs purified from extracellular vesicles, identifying miRNAs as well as myriad other noncoding RNAs (ncRNAs) and ncRNA fragments. Our results establish the utility of OTTR for automation-friendly, low-bias, end-to-end RNA sequence inventories of complex ncRNA samples.

Molecular characterization of both transesterification reactions of the group II intron circularization pathway

Group II introns can self-splice from RNA transcripts through branching, hydrolysis and circularization, being released as lariats, linear introns and circles, respectively. In contrast to branching, the circularization pathway is mostly based on assumptions and has been largely overlooked. Here, we address the molecular details of both transesterification reactions of the group II intron circularization pathway in vivo. We show that free E1 is recruited by the intron through base pairing interactions and that released intron circles can generate free E1 by the spliced exon reopening reaction. The first transesterification reaction was found to be induced inaccurately by the 3′OH of the terminal residue of free E1 at the 3′ splice site, producing circularization intermediates with heterogeneous 3′ ends. Nevertheless, specific terminal 3′OH, selected by a molecular ruler, was shown to precisely attack the 5′ splice site and release intron circles with 3′–5′ rather than 2′–5′ bonds at their circularization junction. Our work supports a circularization model where the recruitment of free E1 and/or displacement of cis-E1 induce a conformational change of the intron active site from the pre-5′ to the pre-3′ splice site processing conformation, suggesting how circularization might initiate at the 3′ instead of the 5′ splice site.

Structural basis for template switching by a group II intron-encoded non-LTR-retroelement reverse transcriptase

Reverse transcriptases (RTs) can switch template strands during complementary DNA synthesis, enabling them to join discontinuous nucleic acid sequences. Template switching (TS) plays crucial roles in retroviral replication and recombination, is used for adapter addition in RNA-Seq, and may contribute to retroelement fitness by increasing evolutionary diversity and enabling continuous complementary DNA synthesis on damaged templates. Here, we determined an X-ray crystal structure of a TS complex of a group II intron RT bound simultaneously to an acceptor RNA and donor RNA template–DNA primer heteroduplex with a 1-nt 3′-DNA overhang. The structure showed that the 3′ end of the acceptor RNA binds in a pocket formed by an N-terminal extension present in non–long terminal repeat–retroelement RTs and the RT fingertips loop, with the 3′ nucleotide of the acceptor base paired to the 1-nt 3′-DNA overhang and its penultimate nucleotide base paired to the incoming dNTP at the RT active site. Analysis of structure-guided mutations identified amino acids that contribute to acceptor RNA binding and a phenylalanine residue near the RT active site that mediates nontemplated nucleotide addition. Mutation of the latter residue decreased multiple sequential template switches in RNA-Seq. Our results provide new insights into the mechanisms of TS and nontemplated nucleotide addition by RTs, suggest how these reactions could be improved for RNA-Seq, and reveal common structural features for TS by non–long terminal repeat–retroelement RTs and viral RNA–dependent RNA polymerases.

U5 snRNA Interactions With Exons Ensure Splicing Precision

Imperfect conservation of human pre-mRNA splice sites is necessary to produce alternative isoforms. This flexibility is combined with the precision of the message reading frame. Apart from intron-termini GU_AG and the branchpoint A, the most conserved are the exon-end guanine and +5G of the intron start. Association between these guanines cannot be explained solely by base-pairing with U1 snRNA in the early spliceosome complex. U6 succeeds U1 and pairs +5G in the pre-catalytic spliceosome, while U5 binds the exon end. Current U5 snRNA reconstructions by CryoEM cannot explain the conservation of the exon-end G. Conversely, human mutation analyses show that guanines of both exon termini can suppress splicing mutations. Our U5 hypothesis explains the mechanism of splicing precision and the role of these conserved guanines in the pre-catalytic spliceosome. We propose: (1) optimal binding register for human exons and U5-the exon junction positioned at U5Loop1 C₃₉|C₃₈; (2) common mechanism for base-pairing of human U5 snRNA with diverse exons and bacterial Ll.LtrB intron with new loci in retrotransposition-guided by base pair geometry; and (3) U5 plays a significant role in specific exon recognition in the pre-catalytic spliceosome. Statistical analyses showed increased U5 Watson-Crick pairs with the 5'exon in the absence of +5G at the intron start. In 5'exon positions -3 and -5, this effect is specific to U5 snRNA rather than U1 snRNA of the early spliceosome. Increased U5 Watson-Crick pairs with 3'exon position +1 coincide with substitutions of the conserved -3C at the intron 3'end. Based on mutation and X-ray evidence, we propose that -3C pairs with U2 G₃₁ juxtaposing the branchpoint and the 3'intron end. The intron-termini pair, formed in the pre-catalytic spliceosome to be ready for transition after branching, and the early involvement of the 3'intron end ensure that the 3'exon contacts U5 in the pre-catalytic complex. We suggest that splicing precision is safeguarded cooperatively by U5, U6, and U2 snRNAs that stabilize the pre-catalytic complex by Watson-Crick base pairing. In addition, our new U5 model explains the splicing effect of exon-start +1G mutations: U5 Watson-Crick pairs with exon +2C/+3G strongly promote exon inclusion. We discuss potential applications for snRNA therapeutics and gene repair by reverse splicing.

Intracellular DNA transfer events restricted to the genus Convallaria within the Asparagaceae family: Possible mechanisms and potential as genetic markers for biographical studies

Intracellular gene transfer among plant genomes is a common phenomenon. Due to their high conservation and high plastid membrane integrity, chloroplast (cp) genomes incorporate foreign genetic material very rarely. Convallaria is a small monocotyledonous genus consisting of C. keiskei, C. majalis and C. montana. Here, we characterized, analyzed and identified 3.3 and 3.7 kb of mitochondrial DNA sequences in the plastome (MCP) of C. majalis and C. montana, respectively. We identified 6 bp and 23 bp direct repeats and mitochondrial pseudogenes, with rps3, rps19 and rpl10 identified in the MCP region. Additionally, we developed novel plastid molecular genetic markers to differentiate Convallaria spp. based on 21 populations. BEAST and biogeographical analyses suggested that Convallaria separated into Eurasian and North American lineages during the middle Pliocene and originated in East Asia. Vicariance in the genus was followed by dispersal into Europe and southeastern North America. These analyses indicate that the MCP event was restricted to the genus Convallaria of Asparagaceae, in contrast to similar events that occurred in its common ancestors with other families of land plants. However, further mitochondrial and population studies are necessary to understand the integration of the MCP region and gene flow in the genus Convallaria.

Mitochondrial DNA duplication, recombination, and introgression during interspecific hybridization

mtDNA recombination events in yeasts are known, but altered mitochondrial genomes were not completed. Therefore, we analyzed recombined mtDNAs in six Saccharomyces cerevisiae × Saccharomyces paradoxus hybrids in detail. Assembled molecules contain mostly segments with variable length introgressed to other mtDNA. All recombination sites are in the vicinity of the mobile elements, introns in cox1, cob genes and free standing ORF1, ORF4. The transplaced regions involve co-converted proximal exon regions. Thus, these selfish elements are beneficial to the host if the mother molecule is challenged with another molecule for transmission to the progeny. They trigger mtDNA recombination ensuring the transfer of adjacent regions, into the progeny of recombinant molecules. The recombination of the large segments may result in mitotically stable duplication of several genes.

Revealing the high variability on nonconserved core and mobile elements of Austropuccinia psidii and other rust mitochondrial genomes

Mitochondrial genomes are highly conserved in many fungal groups, and they can help characterize the phylogenetic relationships and evolutionary biology of plant pathogenic fungi. Rust fungi are among the most devastating diseases for economically important crops around the world. Here, we report the complete sequence and annotation of the mitochondrial genome of Austropuccinia psidii (syn. Puccinia psidii), the causal agent of myrtle rust. We performed a phylogenomic analysis including the complete mitochondrial sequences from other rust fungi. The genome composed of 93.299 bp has 73 predicted genes, 33 of which encoded nonconserved proteins (ncORFs), representing almost 45% of all predicted genes. A. psidii mtDNA is one of the largest rust mtDNA sequenced to date, most likely due to the abundance of ncORFs. Among them, 33% were within intronic regions of diverse intron groups. Mobile genetic elements invading intron sequences may have played significant roles in size but not shaping of the rust mitochondrial genome structure. The mtDNAs from rust fungi are highly syntenic. Phylogenetic inferences with 14 concatenated mitochondrial proteins encoded by the core genes placed A. psidii according to phylogenetic analysis based on 18S rDNA. Interestingly, cox1, the gene with the greatest number of introns, provided phylogenies not congruent with the core set. For the first time, we identified the proteins encoded by three A. psidii ncORFs using proteomics analyses. Also, the orf208 encoded a transmembrane protein repressed during in vitro morphogenesis. To the best of our knowledge, we presented the first report of a complete mtDNA sequence of a member of the family Sphaerophragmiacea.

Group II Introns Generate Functional Chimeric Relaxase Enzymes with Modified Specificities through Exon Shuffling at Both the RNA and DNA Level

Group II introns are large self-splicing RNA enzymes with a broad but somewhat irregular phylogenetic distribution. These ancient retromobile elements are the proposed ancestors of approximately half the human genome, including the abundant spliceosomal introns and non-long terminal repeat retrotransposons. In contrast to their eukaryotic derivatives, bacterial group II introns have largely been considered as harmful selfish mobile retroelements that parasitize the genome of their host. As a challenge to this view, we recently uncovered a new intergenic trans-splicing pathway that generates an assortment of mRNA chimeras. The ability of group II introns to combine disparate mRNA fragments was proposed to increase the genetic diversity of the bacterial host by shuffling coding sequences. Here, we show that the Ll.LtrB and Ef.PcfG group II introns from Lactococcus lactis and Enterococcus faecalis respectively can both use the intergenic trans-splicing pathway to catalyze the formation of chimeric relaxase mRNAs and functional proteins. We demonstrated that some of these compound relaxase enzymes yield gain-of-function phenotypes, being significantly more efficient than their precursor wild-type enzymes at supporting bacterial conjugation. We also found that relaxase enzymes with shuffled functional domains are produced in biologically relevant settings under natural expression levels. Finally, we uncovered examples of lactococcal chimeric relaxase genes with junctions exactly at the intron insertion site. Overall, our work demonstrates that the genetic diversity generated by group II introns, at the RNA level by intergenic trans-splicing and at the DNA level by recombination, can yield new functional enzymes with shuffled exons, which can lead to gain-of-function phenotypes.

Exon and protein positioning in a pre-catalytic group II intron RNP primed for splicing

Group II introns are the putative progenitors of nuclear spliceosomal introns and use the same two-step splicing pathway. In the cell, the intron RNA forms a ribonucleoprotein (RNP) complex with the intron-encoded protein (IEP), which is essential for splicing. Although structures of spliced group II intron RNAs and RNP complexes have been characterized, structural insights into the splicing process remain enigmatic due to lack of pre-catalytic structural models. Here, we report two cryo-EM structures of endogenously produced group II intron RNPs trapped in their pre-catalytic state. Comparison of the catalytically activated precursor RNP to its previously reported spliced counterpart allowed identification of key structural rearrangements accompanying splicing, including a remodeled active site and engagement of the exons. Importantly, altered RNA-protein interactions were observed upon splicing among the RNP complexes. Furthermore, analysis of the catalytically inert precursor RNP demonstrated the structural impact of the formation of the active site on RNP architecture. Taken together, our results not only fill a gap in understanding the structural basis of IEP-assisted group II intron splicing, but also provide parallels to evolutionarily related spliceosomal splicing.

Group II intron as cold sensor for self-preservation and bacterial conjugation

Group II introns are self-splicing ribozymes and mobile genetic elements. Splicing is required for both expression of the interrupted host gene and intron retromobility. For the pRS01 plasmid-encoded Lactococcus lactis group II intron, Ll.LtrB, splicing enables expression of the intron's host relaxase protein. Relaxase, in turn, initiates horizontal transfer of the conjugative pRS01 plasmid and stimulates retrotransposition of the intron. Little is known about how splicing of bacterial group II introns is influenced by environmental conditions. Here, we show that low temperatures can inhibit Ll.LtrB intron splicing. Whereas autocatalysis is abolished in the cold, splicing is partially restored by the intron-encoded protein (IEP). Structure profiling reveals cold-induced disruptions of key tertiary interactions, suggesting that a kinetic trap prevents the intron RNA from assuming its native state. Interestingly, while reduced levels of transcription and splicing lead to a paucity of excised intron in the cold, levels of relaxase mRNA are maintained, partially due to diminished intron-mediated mRNA targeting, allowing intron spread by conjugal transfer. Taken together, this study demonstrates not only the intrinsic cold sensitivity of group II intron splicing and the role of the IEP for cold-stress adaptation, but also maintenance of horizontal plasmid and intron transfer under cold-shock.

The circle to lariat ratio of the Ll.LtrB group II intron from Lactococcus lactis is greatly influenced by a variety of biological determinants in vivo

Bacterial group II introns mostly behave as versatile retromobile genetic elements going through distinct cycles of gain and loss. These large RNA molecules are also ribozymes splicing autocatalytically from their interrupted pre-mRNA transcripts by two different concurrent pathways, branching and circularization. These two splicing pathways were shown to release in bacterial cells significant amounts of branched intron lariats and perfect end-to-end intron circles respectively. On one hand, released intron lariats can invade new sites in RNA and/or DNA by reverse branching while released intron circles are dead end spliced products since they cannot reverse splice through circularization. The presence of two parallel and competing group II intron splicing pathways in bacteria led us to investigate the conditions that influence the overall circle to lariat ratio in vivo. Here we unveil that removing a prominent processing site within the Ll.LtrB group II intron, raising growth temperature of Lactococcus lactis host cells and increasing the expression level of the intron-interrupted gene all increased the relative amount of released intron circles compared to lariats. Strengthening and weakening the base pairing interaction between the intron and its upstream exon respectively increased and decreased the overall levels of released intron circles in comparison to lariats. Host environment was also found to impact the circle to lariat ratio of the Ll.LtrB and Ll.RlxA group II introns from L. lactis and the Ef.PcfG intron from Enterococcus faecalis. Overall, our data show that multiple factors significantly influence the balance between released intron circles and lariats in bacterial cells.

Transitions between the steps of forward and reverse splicing of group IIC introns

Group II introns are mobile genetic elements that perform both self-splicing and intron mobility reactions. These ribozymes are comprised of a catalytic RNA core that binds to an intron-encoded protein (IEP) to form a ribonucleoprotein (RNP) complex. Splicing proceeds through two competing reactions: hydrolysis or branching. Group IIC intron ribozymes have a minimal RNA architecture, and splice almost exclusively through hydrolysis in ribozyme reactions. Addition of the IEP allows the splicing reaction to form branched lariat RNPs capable of intron mobility. Here we examine ribozyme splicing, IEP-dependent splicing, and mobility reactions of a group IIC intron from the thermophilic bacterium Thermoanerobacter italicus (Ta.it.I1). We show that Ta.it.I1 is highly active for ribozyme activity, forming linear hydrolytic intron products. Addition of purified IEP switches activity to the canonical lariat forming splicing reaction. We demonstrate that the Ta.it.I1 group IIC intron coordinates the progression of the forward splicing reaction through a π-π' interaction between intron domains II and VI. We further show that branched splicing is supported in the absence of the IEP when the π-π' interaction is mutated. We also investigated the regulation of the two steps of reverse splicing during intron mobility into DNA substrates. Using a fluorescent mobility assay that simultaneously visualizes all steps of intron integration into DNA, we show that completion of reverse splicing is tightly coupled to cDNA synthesis regardless of mutation of the π-π' interaction.

Physicochemical Foundations of Life that Direct Evolution: Chance and Natural Selection are not Evolutionary Driving Forces

The current framework of evolutionary theory postulates that evolution relies on random mutations generating a diversity of phenotypes on which natural selection acts. This framework was established using a top-down approach as it originated from Darwinism, which is based on observations made of complex multicellular organisms and, then, modified to fit a DNA-centric view. In this article, it is argued that based on a bottom-up approach starting from the physicochemical properties of nucleic and amino acid polymers, we should reject the facts that (i) natural selection plays a dominant role in evolution and (ii) the probability of mutations is independent of the generated phenotype. It is shown that the adaptation of a phenotype to an environment does not correspond to organism fitness, but rather corresponds to maintaining the genome stability and integrity. In a stable environment, the phenotype maintains the stability of its originating genome and both (genome and phenotype) are reproduced identically. In an unstable environment (i.e., corresponding to variations in physicochemical parameters above a physiological range), the phenotype no longer maintains the stability of its originating genome, but instead influences its variations. Indeed, environment- and cellular-dependent physicochemical parameters define the probability of mutations in terms of frequency, nature, and location in a genome. Evolution is non-deterministic because it relies on probabilistic physicochemical rules, and evolution is driven by a bidirectional interplay between genome and phenotype in which the phenotype ensures the stability of its originating genome in a cellular and environmental physicochemical parameter-depending manner.

A group II intron-encoded protein interacts with the cellular replicative machinery through the β-sliding clamp

Group II introns are self-splicing mobile genetic retroelements. The spliced intron RNA and the intron-encoded protein (IEP) form ribonucleoprotein particles (RNPs) that recognize and invade specific DNA target sites. The IEP is a reverse transcriptase/maturase that may bear a C-terminal endonuclease domain enabling the RNP to cleave the target DNA strand to prime reverse transcription. However, some mobile introns, such as RmInt1, lack the En domain but nevertheless retrohome efficiently to transient single-stranded DNA target sites at a DNA replication fork. Their mobility is associated with host DNA replication, and they use the nascent lagging strand as a primer for reverse transcription. We searched for proteins that interact with RmInt1 RNPs and direct these RNPs to the DNA replication fork. Co-immunoprecipitation assays suggested that DnaN (the β-sliding clamp), a component of DNA polymerase III, interacts with the protein component of the RmInt1 RNP. Pulldown assays, far-western blots and biolayer interferometry supported this interaction. Peptide binding assays also identified a putative DnaN-interacting motif in the RmInt1 IEP structurally conserved in group II intron IEPs. Our results suggest that intron RNP interacts with the β-sliding clamp of the DNA replication machinery, favouring reverse splicing into the transient ssDNA at DNA replication forks.

Template-switching mechanism of a group II intron-encoded reverse transcriptase and its implications for biological function and RNA-Seq

The reverse transcriptases (RTs) encoded by mobile group II introns and other non-LTR retroelements differ from retroviral RTs in being able to template-switch efficiently from the 5' end of one template to the 3' end of another with little or no complementarity between the donor and acceptor templates. Here, to establish a complete kinetic framework for the reaction and to identify conditions that more efficiently capture acceptor RNAs or DNAs, we used a thermostable group II intron RT (TGIRT; GsI-IIC RT) that can template switch directly from synthetic RNA template/DNA primer duplexes having either a blunt end or a 3'-DNA overhang end. We found that the rate and amplitude of template switching are optimal from starter duplexes with a single nucleotide 3'-DNA overhang complementary to the 3' nucleotide of the acceptor RNA, suggesting a role for nontemplated nucleotide addition of a complementary nucleotide to the 3' end of cDNAs synthesized from natural templates. Longer 3'-DNA overhangs progressively decreased the template-switching rate, even when complementary to the 3' end of the acceptor template. The reliance on only a single bp with the 3' nucleotide of the acceptor together with discrimination against mismatches and the high processivity of group II intron RTs enable synthesis of full-length DNA copies of nucleic acids beginning directly at their 3' end. We discuss the possible biological functions of the template-switching activity of group II intron- and other non-LTR retroelement-encoded RTs, as well as the optimization of this activity for adapter addition in RNA- and DNA-Seq protocols.

Improved TGIRT-seq methods for comprehensive transcriptome profiling with decreased adapter dimer formation and bias correction

Thermostable group II intron reverse transcriptases (TGIRTs) with high fidelity and processivity have been used for a variety of RNA sequencing (RNA-seq) applications, including comprehensive profiling of whole-cell, exosomal, and human plasma RNAs; quantitative tRNA-seq based on the ability of TGIRT enzymes to give full-length reads of tRNAs and other structured small ncRNAs; high-throughput mapping of post-transcriptional modifications; and RNA structure mapping. Here, we improved TGIRT-seq methods for comprehensive transcriptome profiling by rationally designing RNA-seq adapters that minimize adapter dimer formation. Additionally, we developed biochemical and computational methods for remediating 5′- and 3′-end biases, the latter based on a random forest regression model that provides insight into the contribution of different factors to these biases. These improvements, some of which may be applicable to other RNA-seq methods, increase the efficiency of TGIRT-seq library construction and improve coverage of very small RNAs, such as miRNAs. Our findings provide insight into the biochemical basis of 5′- and 3′-end biases in RNA-seq and suggest general approaches for remediating biases and decreasing adapter dimer formation.

DNA cleavage and reverse splicing of ribonucleoprotein particles reconstituted in vitro with linear RmInt1 RNA

The RmInt1 group II intron is an efficient self-splicing mobile retroelement that catalyzes its own excision as lariat, linear and circular molecules. In vivo, the RmInt1 lariat and the reverse transcriptase (IEP) it encodes form a ribonucleoprotein particle (RNP) that recognizes the DNA target for site-specific full intron insertion via a two-step reverse splicing reaction. RNPs containing linear group II intron RNA are generally thought to be unable to complete the reverse splicing reaction. Here, we show that reconstituted in vitro RNPs containing linear RmInt1 ΔORF RNA can mediate the cleavage of single-stranded DNA substrates in a very precise manner with the attachment of the intron RNA to the 3´exon as the first step of a reverse splicing reaction. Notably, we also observe molecules in which the 5´exon is linked to the RmInt1 RNA, suggesting the completion of the reverse splicing reaction, albeit rather low and inefficiently. That process depends on DNA target recognition and can be successful completed by RmInt1 RNPs with linear RNA displaying 5´ modifications.

Group II Intron RNPs and Reverse Transcriptases: From Retroelements to Research Tools

Group II introns, self-splicing retrotransposons, serve as both targets of investigation into their structure, splicing, and retromobility and a source of tools for genome editing and RNA analysis. Here, we describe the first cryo-electron microscopy (cryo-EM) structure determination, at 3.8-4.5 Å, of a group II intron ribozyme complexed with its encoded protein, containing a reverse transcriptase (RT), required for RNA splicing and retromobility. We also describe a method called RIG-seq using a retrotransposon indicator gene for high-throughput integration profiling of group II introns and other retrotransposons. Targetrons, RNA-guided gene targeting agents widely used for bacterial genome engineering, are described next. Finally, we detail thermostable group II intron RTs, which synthesize cDNAs with high accuracy and processivity, for use in various RNA-seq applications and relate their properties to a 3.0-Å crystal structure of the protein poised for reverse transcription. Biological insights from these group II intron revelations are discussed.

Origins and Evolution of the Global RNA Virome

Viruses with RNA genomes dominate the eukaryotic virome, reaching enormous diversity in animals and plants. The recent advances of metaviromics prompted us to perform a detailed phylogenomic reconstruction of the evolution of the dramatically expanded global RNA virome. The only universal gene among RNA viruses is the gene encoding the RNA-dependent RNA polymerase (RdRp). We developed an iterative computational procedure that alternates the RdRp phylogenetic tree construction with refinement of the underlying multiple-sequence alignments. The resulting tree encompasses 4,617 RNA virus RdRps and consists of 5 major branches; 2 of the branches include positive-sense RNA viruses, 1 is a mix of positive-sense (+) RNA and double-stranded RNA (dsRNA) viruses, and 2 consist of dsRNA and negative-sense (-) RNA viruses, respectively. This tree topology implies that dsRNA viruses evolved from +RNA viruses on at least two independent occasions, whereas -RNA viruses evolved from dsRNA viruses. Reconstruction of RNA virus evolution using the RdRp tree as the scaffold suggests that the last common ancestors of the major branches of +RNA viruses encoded only the RdRp and a single jelly-roll capsid protein. Subsequent evolution involved independent capture of additional genes, in particular, those encoding distinct RNA helicases, enabling replication of larger RNA genomes and facilitating virus genome expression and virus-host interactions. Phylogenomic analysis reveals extensive gene module exchange among diverse viruses and horizontal virus transfer between distantly related hosts. Although the network of evolutionary relationships within the RNA virome is bound to further expand, the present results call for a thorough reevaluation of the RNA virus taxonomy.IMPORTANCE The majority of the diverse viruses infecting eukaryotes have RNA genomes, including numerous human, animal, and plant pathogens. Recent advances of metagenomics have led to the discovery of many new groups of RNA viruses in a wide range of hosts. These findings enable a far more complete reconstruction of the evolution of RNA viruses than was attainable previously. This reconstruction reveals the relationships between different Baltimore classes of viruses and indicates extensive transfer of viruses between distantly related hosts, such as plants and animals. These results call for a major revision of the existing taxonomy of RNA viruses.

Bacterial group II introns generate genetic diversity by circularization and trans-splicing from a population of intron-invaded mRNAs

Group II introns are ancient retroelements that significantly shaped the origin and evolution of contemporary eukaryotic genomes. These self-splicing ribozymes share a common ancestor with the telomerase enzyme, the spliceosome machinery as well as the highly abundant spliceosomal introns and non-LTR retroelements. More than half of the human genome thus consists of various elements that evolved from ancient group II introns, which altogether significantly contribute to key functions and genetic diversity in eukaryotes. Similarly, group II intron-related elements in bacteria such as abortive phage infection (Abi) retroelements, diversity generating retroelements (DGRs) and some CRISPR-Cas systems have evolved to confer important functions to their hosts. In sharp contrast, since bacterial group II introns are scarce, irregularly distributed and frequently spread by lateral transfer, they have mainly been considered as selfish retromobile elements with no beneficial function to their host. Here we unveil a new group II intron function that generates genetic diversity at the RNA level in bacterial cells. We demonstrate that Ll.LtrB, the model group II intron from Lactococcus lactis, recognizes specific sequence motifs within cellular mRNAs by base pairing, and invades them by reverse splicing. Subsequent splicing of ectopically inserted Ll.LtrB, through circularization, induces a novel trans-splicing pathway that generates exon 1-mRNA and mRNA-mRNA intergenic chimeras. Our data also show that recognition of upstream alternative circularization sites on intron-interrupted mRNAs release Ll.LtrB circles harboring mRNA fragments of various lengths at their splice junction. Intergenic trans-splicing and alternative circularization both produce novel group II intron splicing products with potential new functions. Overall, this work describes new splicing pathways in bacteria that generate, similarly to the spliceosome in eukaryotes, genetic diversity at the RNA level while providing additional functional and evolutionary links between group II introns, spliceosomal introns and the spliceosome.

Testing the retroelement invasion hypothesis for the emergence of the ancestral eukaryotic cell

Phylogenetic evidence suggests that the invasion and proliferation of retroelements, selfish mobile genetic elements that copy and paste themselves within a host genome, was one of the early evolutionary events in the emergence of eukaryotes. Here we test the effects of this event by determining the pressures retroelements exert on simple genomes. We transferred two retroelements, human LINE-1 and the bacterial group II intron Ll.LtrB, into bacteria, and find that both are functional and detrimental to growth. We find, surprisingly, that retroelement lethality and proliferation are enhanced by the ability to perform eukaryotic-like nonhomologous end-joining (NHEJ) DNA repair. We show that the only stable evolutionary consequence in simple cells is maintenance of retroelements in low numbers, suggesting how retrotransposition rates and costs in early eukaryotes could have been constrained to allow proliferation. Our results suggest that the interplay between NHEJ and retroelements may have played a fundamental and previously unappreciated role in facilitating the proliferation of retroelements, elements of which became the ancestors of the spliceosome components in eukaryotes.