Evolution of group II introns

Protein-free catalysis of DNA hydrolysis and self-integration by a ribozyme

Group II introns are ancient self-splicing ribozymes and retrotransposons. Though long speculated to have originated before translation, their dependence on intron-encoded proteins for splicing and mobility has cast doubt on this hypothesis. While some group II introns are known to retain part of their catalytic repertoire in the absence of protein cofactors, protein-free complete reverse splicing of a group II intron into a DNA target has never been demonstrated. Here, we demonstrate the complete independence of a group II intron from protein cofactors in all intron-catalyzed reactions. The ribozyme is capable of fully reverse splicing into single-stranded DNA targets in vitro, readily hydrolyzes DNA substrates and is even able to unwind and react with stably duplexed DNA. Our findings make a protein-free origin for group II introns plausible by expanding their known catalytic capabilities beyond what would be needed to survive the transition from RNA to DNA genomes. Furthermore, the intron's capacity to react with both single and double-stranded DNA in conjunction with its expanded sequence recognition may represent a promising starting point for the development of protein-free genomic editing tools.

PCIS1, Encoded by a Pentatricopeptide Protein Co-expressed Gene, Is Required for Splicing of Three Mitochondrial nad Transcripts in Angiosperms

Group II introns are large catalytic RNAs, which reside mainly within genes encoding respiratory complex I (CI) subunits in angiosperms' mitochondria. Genetic and biochemical analyses led to the identification of many nuclear-encoded factors that facilitate the splicing of the degenerated organellar introns in plants. Here, we describe the analysis of the pentatricopeptide repeat (PPR) co-expressed intron splicing-1 (PCIS1) factor, which was identified in silico by its co-expression pattern with many PPR proteins. PCIS1 is well conserved in land plants but has no sequence similarity with any known protein motifs. PCIS1 mutant lines are arrested in embryogenesis and can be maintained by the temporal expression of the gene under the embryo-specific ABI3 promoter. The pABI3::PCIS1 mutant plants display low germination and stunted growth phenotypes. RNA-sequencing and quantitative RT-PCR analyses of wild-type and mutant plants indicated that PCIS1 is a novel splicing cofactor that is pivotal for the maturation of several nad transcripts in Arabidopsis mitochondria. These phenotypes are tightly associated with respiratory CI defects and altered plant growth. Our data further emphasize the key roles of nuclear-encoded cofactors that regulate the maturation and expression of mitochondrial transcripts for the biogenesis of the oxidative phosphorylation system, and hence for plant physiology. The discovery of novel splicing factors other than typical RNA-binding proteins suggests further complexity of splicing mechanisms in plant mitochondria.

Interconnected roles of fungal nuclear- and intron-encoded maturases: at the crossroads of mitochondrial intron splicing

Group I and II introns are large catalytic RNAs (ribozymes) that are frequently encountered in fungal mitochondrial genomes. The discovery of respiratory mutants linked to intron splicing defects demonstrated that for the efficient removal of organellar introns there appears to be a requirement of protein splicing factors. These splicing factors can be intron-encoded proteins with maturase activities that usually promote the splicing of the introns that encode them (cis-acting) and/or nuclear-encoded factors that can promote the splicing of a range of different introns (trans-acting). Compared to plants organellar introns, fungal mitochondrial intron splicing is still poorly explored, especially in terms of the synergy of nuclear factors with intron-encoded maturases that has direct impact on splicing through their association with intron RNA. In addition, nuclear-encoded accessory factors might drive the splicing impetus through translational activation, mitoribosome assembly, and phosphorylation-mediated RNA turnover. This review explores protein-assisted splicing of introns by nuclear and mitochondrial-encoded maturases as a means of mitonuclear interplay that could respond to environmental and developmental factors promoting phenotypic adaptation and potentially speciation. It also highlights key evolutionary events that have led to changes in structure and ATP-dependence to accommodate the dual functionality of nuclear and organellar splicing factors.

Robust genome and cell engineering via in vitro and in situ circularized RNAs

Circularization can improve RNA persistence, yet simple and scalable approaches to achieve this are lacking. Here we report two methods that facilitate the pursuit of circular RNAs (cRNAs): cRNAs developed via in vitro circularization using group II introns, and cRNAs developed via in-cell circularization by the ubiquitously expressed RtcB protein. We also report simple purification protocols that enable high cRNA yields (40-75%) while maintaining low immune responses. These methods and protocols facilitate a broad range of applications in stem cell engineering as well as robust genome and epigenome targeting via zinc finger proteins and CRISPR-Cas9. Notably, cRNAs bearing the encephalomyocarditis internal ribosome entry enabled robust expression and persistence compared with linear capped RNAs in cardiomyocytes and neurons, which highlights the utility of cRNAs in these non-dividing cells. We also describe genome targeting via deimmunized Cas9 delivered as cRNA and a long-range multiplexed protein engineering methodology for the combinatorial screening of deimmunized protein variants that enables compatibility between persistence of expression and immunogenicity in cRNA-delivered proteins. The cRNA toolset will aid research and the development of therapeutics.

Epigenetics and transgenerational inheritance

Epigenetic modifications can alter the function of genes. The epigenetics changes are caused by environmental effects, which lead to chemical modifications of the DNA or the chromatin. The mechanisms involve the influence of small interfering siRNAs on gene silencing. Epigenetic changes normally last only during the life-time of an individual and are erased in embryos and eggs for a naive progeny. The genomes are reprogrammed and the chemical modifications removed to restart the next generation. However, there are mechanisms that allow the genome to escape from such a clearing effect so that modifications can be transmitted to one or more subsequent generations. In the germline of animal cells small RNAs, including piRNAs, have evolved which guarantee a higher degree of fidelity for transmission of genetic information, guarding especially against the detrimental effect caused by transposon activity. piRNA is essential for transposon silencing for survival of a species and protection of subsequent generations. Inactivation of piRNA results in abundant transposon activity and sperm infertility. The effect in humans has been described but is less distinct. Some stress-induced epigenetic changes are transitory in mice and can be reversed by a change of environment or lifestyle.

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.

Root Primordium Defective 1 Encodes an Essential PORR Protein Required for the Splicing of Mitochondria-Encoded Group II Introns and for Respiratory Complex I Biogenesis

Cellular respiration involves complex organellar metabolic activities that are pivotal for plant growth and development. Mitochondria contain their own genetic system (mitogenome, mtDNA), which encodes key elements of the respiratory machinery. Plant mtDNAs are notably larger than their counterparts in Animalia, with complex genome organization and gene expression characteristics. The maturation of the plant mitochondrial transcripts involves extensive RNA editing, trimming and splicing events. These essential processing steps rely on the activities of numerous nuclear-encoded cofactors, which may also play key regulatory roles in mitochondrial biogenesis and function and hence in plant physiology. Proteins that harbor the plant organelle RNA recognition (PORR) domain are represented in a small gene family in plants. Several PORR members, including WTF1, WTF9 and LEFKOTHEA, are known to act in the splicing of organellar group II introns in angiosperms. The AT4G33495 gene locus encodes an essential PORR protein in Arabidopsis, termed ROOT PRIMORDIUM DEFECTIVE 1 (RPD1). A null mutation of At.RPD1 causes arrest in early embryogenesis, while the missense mutant lines, rpd1.1 and rpd1.2, exhibit a strong impairment in root development and retarded growth phenotypes, especially under high-temperature conditions. Here, we further show that RPD1 functions in the splicing of introns that reside in the coding regions of various complex I (CI) subunits (i.e. nad2, nad4, nad5 and nad7), as well as in the maturation of the ribosomal rps3 pre-RNA in Arabidopsis mitochondria. The altered growth and developmental phenotypes and modified respiration activities are tightly correlated with respiratory chain CI defects in rpd1 mutants.

RNA-Small-Molecule Interaction: Challenging the "Undruggable" Tag

RNA targeting, specifically with small molecules, is a relatively new and rapidly emerging avenue with the promise to expand the target space in the drug discovery field. From being "disregarded" as an "undruggable" messenger molecule to FDA approval of an RNA-targeting small-molecule drug Risdiplam, a radical change in perspective toward RNA has been observed in the past decade. RNAs serve important regulatory functions beyond canonical protein synthesis, and their dysregulation has been reported in many diseases. A deeper understanding of RNA biology reveals that RNA molecules can adopt a variety of structures, carrying defined binding pockets that can accommodate small-molecule drugs. Due to its functional diversity and structural complexity, RNA can be perceived as a prospective target for therapeutic intervention. This perspective highlights the proof of concept of RNA-small-molecule interactions, exemplified by targeting of various transcripts with functional modulators. The advent of RNA-oriented knowledge would help expedite drug discovery.

Hydrolytic endonucleolytic ribozyme (HYER) is programmable for sequence-specific DNA cleavage

Ribozymes are catalytic RNAs with diverse functions including self-splicing and polymerization. This work aims to discover natural ribozymes that behave as hydrolytic and sequence-specific DNA endonucleases, which could be repurposed as DNA manipulation tools. Focused on bacterial group II-C introns, we found that many systems without intron-encoded protein propagate multiple copies in their resident genomes. These introns, named HYdrolytic Endonucleolytic Ribozymes (HYERs), cleaved RNA, single-stranded DNA, bubbled double-stranded DNA (dsDNA), and plasmids in vitro. HYER1 generated dsDNA breaks in the mammalian genome. Cryo-electron microscopy analysis revealed a homodimer structure for HYER1, where each monomer contains a Mg2+-dependent hydrolysis pocket and captures DNA complementary to the target recognition site (TRS). Rational designs including TRS extension, recruiting sequence insertion, and heterodimerization yielded engineered HYERs showing improved specificity and flexibility for DNA manipulation.

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.

Exploring the highly reduced spliceosome of Pseudoloma neurophilia

Spliceosomal introns evolved early in eukaryogenesis, originating from self-splicing group II introns that invaded the proto-eukaryotic genome1. Elements of these ribozymes, now called snRNAs (U1, U2, U4, U5, U6), were co-opted to excise these invasive elements. Prior to eukaryotic diversification, the spliceosome is predicted to have accumulated hundreds of proteins2. This early complexification has obscured our understanding of spliceosomal evolution. Reduced systems with few introns and tiny spliceosomes give insights into the plasticity of the splicing reaction and provide an opportunity to study the evolution of the spliceosome3,4. Microsporidia are intracellular parasites possessing extremely reduced genomes that have lost many, and in some instances all, introns5. In the purportedly intron-lacking genome of the microsporidian Pseudoloma neurophilia6, we identified two introns that are spliced at high levels. Furthermore, with only 14 predicted proteins, the P. neurophilia spliceosome could be the smallest known. Intriguingly, the few proteins retained are divergent compared to canonical orthologs. Even the central spliceosomal protein Prp8, which originated from the proteinaceous component of group II introns, is extremely divergent. This is unusual given that Prp8 is highly conserved across eukaryotes, including other microsporidia. All five P. neurophilia snRNAs are present, and all but U2 have diverged extensively, likely resulting from the loss of interacting proteins. Despite this divergence, U1 and U2 are predicted to pair with intron sequences more extensively than previously described. The P. neurophilia spliceosome is retained to splice a mere two introns and, with few proteins and reliance on RNA-RNA interactions, could function in a manner more reminiscent of presumed ancestral splicing.

Structural insights into intron catalysis and dynamics during splicing

The group II intron ribonucleoprotein is an archetypal splicing system with numerous mechanistic parallels to the spliceosome, including excision of lariat introns1,2. Despite the importance of branching in RNA metabolism, structural understanding of this process has remained elusive. Here we present a comprehensive analysis of three single-particle cryogenic electron microscopy structures captured along the splicing pathway. They reveal the network of molecular interactions that specifies the branchpoint adenosine and positions key functional groups to catalyse lariat formation and coordinate exon ligation. The structures also reveal conformational rearrangements of the branch helix and the mechanism of splice site exchange that facilitate the transition from branching to ligation. These findings shed light on the evolution of splicing and highlight the conservation of structural components, catalytic mechanism and dynamical strategies retained through time in premessenger RNA splicing machines.

Transposon-encoded nucleases use guide RNAs to promote their selfish spread

Insertion sequences are compact and pervasive transposable elements found in bacteria, which encode only the genes necessary for their mobilization and maintenance1. IS200- and IS605-family transposons undergo 'peel-and-paste' transposition catalysed by a TnpA transposase2, but they also encode diverse, TnpB- and IscB-family proteins that are evolutionarily related to the CRISPR-associated effectors Cas12 and Cas9, respectively3,4. Recent studies have demonstrated that TnpB and IscB function as RNA-guided DNA endonucleases5,6, but the broader biological role of this activity has remained enigmatic. Here we show that TnpB and IscB are essential to prevent permanent transposon loss as a consequence of the TnpA transposition mechanism. We selected a family of related insertion sequences from Geobacillus stearothermophilus that encode several TnpB and IscB orthologues, and showed that a single TnpA transposase was broadly active for transposon mobilization. The donor joints formed upon religation of transposon-flanking sequences were efficiently targeted for cleavage by RNA-guided TnpB and IscB nucleases, and co-expression of TnpB and TnpA led to substantially greater transposon retention relative to conditions in which TnpA was expressed alone. Notably, TnpA and TnpB also stimulated recombination frequencies, surpassing rates observed with TnpB alone. Collectively, this study reveals that RNA-guided DNA cleavage arose as a primal biochemical activity to bias the selfish inheritance and spread of transposable elements, which was later co-opted during the evolution of CRISPR-Cas adaptive immunity for antiviral defence.

Circular RNA vaccine in disease prevention and treatment

CircRNAs are a class of single-stranded RNAs with covalently linked head-to-tail topology. In the decades since its initial discovery, their biogenesis, regulation, and function have rapidly disclosed, permitting a better understanding and adoption of them as new tools for medical applications. With the development of biotechnology and molecular medicine, artificial circRNAs have been engineered as a novel class of vaccines for disease treatment and prevention. Unlike the linear mRNA vaccine which applications were limited by its instability, inefficiency, and innate immunogenicity, circRNA vaccine which incorporate internal ribosome entry sites (IRESs) and open reading frame (ORF) provides an improved approach to RNA-based vaccination with safety, stability, simplicity of manufacture, and scalability. However, circRNA vaccines are at an early stage, and their optimization, delivery and applications require further development and evaluation. In this review, we comprehensively describe circRNA vaccine, including their history and superiority. We also summarize and discuss the current methodological research for circRNA vaccine preparation, including their design, synthesis, and purification. Finally, we highlight the delivery options of circRNA vaccine and its potential applications in diseases treatment and prevention. Considering their unique high stability, low immunogenicity, protein/peptide-coding capacity and special closed-loop construction, circRNA vaccine, and circRNA-based therapeutic platforms may have superior application prospects in a broad range of diseases.

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.

Transposon-encoded nucleases use guide RNAs to selfishly bias their inheritance

Insertion sequences (IS) are compact and pervasive transposable elements found in bacteria, which encode only the genes necessary for their mobilization and maintenance. IS 200 /IS 605 elements undergo 'peel-and-paste' transposition catalyzed by a TnpA transposase, but intriguingly, they also encode diverse, TnpB- and IscB-family proteins that are evolutionarily related to the CRISPR-associated effectors Cas12 and Cas9, respectively. Recent studies demonstrated that TnpB-family enzymes function as RNA-guided DNA endonucleases, but the broader biological role of this activity has remained enigmatic. Here we show that TnpB/IscB are essential to prevent permanent transposon loss as a consequence of the TnpA transposition mechanism. We selected a family of related IS elements from Geobacillus stearothermophilus that encode diverse TnpB/IscB orthologs, and showed that a single TnpA transposase was active for transposon excision. The donor joints formed upon religation of IS-flanking sequences were efficiently targeted for cleavage by RNA-guided TnpB/IscB nucleases, and co-expression of TnpB together with TnpA led to significantly greater transposon retention, relative to conditions in which TnpA was expressed alone. Remarkably, TnpA and TnpB/IscB recognize the same AT-rich transposon-adjacent motif (TAM) during transposon excision and RNA-guided DNA cleavage, respectively, revealing a striking convergence in the evolution of DNA sequence specificity between collaborating transposase and nuclease proteins. Collectively, our study reveals that RNA-guided DNA cleavage is a primal biochemical activity that arose to bias the selfish inheritance and spread of transposable elements, which was later co-opted during the evolution of CRISPR-Cas adaptive immunity for antiviral defense.

Plant mitochondrial introns as genetic markers - conservation and variation

Plant genomes are comprised of nuclear, plastid and mitochondrial components characterized by different patterns of inheritance and evolution. Genetic markers from the three genomes provide complementary tools for investigations of inheritance, genetic relationships and phenotypic contributions. Plant mitochondrial genomes are challenging for universal marker development because they are highly variable in terms of size, gene order and intergenic sequences and highly conserved with respect to protein-coding sequences. PCR amplification of introns with primers that anneal to conserved, flanking exons is effective for the development of polymorphic nuclear genome markers. The potential for plant mitochondrial intron polymorphisms to distinguish between congeneric species or intraspecific varieties has not been systematically investigated and is possibly constrained by requirements for intron secondary structure and interactions with co-evolved organelle intron splicing factors. To explore the potential for broadly applicable plant mitochondrial intron markers, PCR primer sets based upon conserved sequences flanking 11 introns common to seven angiosperm species were tested across a range of plant orders. PCR-amplified introns were screened for indel polymorphisms among a group of cross-compatible Citrus species and relatives; two Raphanus sativus mitotypes; representatives of the two Phaseolus vulgaris gene pools; and congeneric pairs of Cynodon, Cenchrus, Solanum, and Vaccinium species. All introns were successfully amplified from each plant entry. Length polymorphisms distinguishable by gel electrophoresis were common among genera but infrequent within genera. Sequencing of three introns amplified from 16 entries identified additional short indel polymorphisms and nucleotide substitutions that separated Citrus, Cynodon, Cenchrus and Vaccinium congeners, but failed to distinguish Solanum congeners or representatives of the Phaseolus vulgaris major gene pools. The ability of primer sets to amplify a wider range of plant species' introns and the presence of intron polymorphisms that distinguish congeners was confirmed by in silico analysis. While mitochondrial intron variation is limited in comparison to nuclear introns, these exon-based primer sets provide robust tools for the amplification of mitochondrial introns across a wide range of plant species wherein useful polymorphisms can be identified.

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.

A New EBS2b-IBS2b Base Paring (A-8/T-8) Improved the Gene-Targeting Efficiency of Thermotargetron in Escherichia coli

Thermophilic group II intron is one type of retrotransposon composed of intron RNA and intron-encoded protein (IEP), which can be utilized in gene targeting by harnessing their novel ribozyme-based DNA integration mechanism termed "retrohoming." It is mediated by a ribonucleoprotein (RNP) complex that contains the excised intron lariat RNA and an IEP with reverse transcriptase (RT) activity. The RNP recognizes targeting sites by exon-binding sequences 2 (EBS2)/intron-binding sequences 2 (IBS2), EBS1/IBS1, and EBS3/IBS3 bases pairing. Previously, we developed the TeI3c/4c intron as a thermophilic gene targeting system-Thermotargetron (TMT). However, we found that the targeting efficiency of TMT varies significantly at different targeting sites, which leads to a relatively low success rate. To further improve the success rate and gene-targeting efficiency of TMT, we constructed a Random Gene-targeting Plasmids Pool (RGPP) to analyze the sequence recognition preference of TMT. A new base pairing, located at the -8 site between EBS2/IBS2 and EBS1/IBS1 (named EBS2b-IBS2b), increased the success rate (2.45- to 5.07-fold) and significantly improved gene-targeting efficiency of TMT. A computer algorithm (TMT 1.0), based on the newly discovered sequence recognition roles, was also developed to facilitate the design of TMT gene-targeting primers. The present work could essentially expand the practicalities of TMT in the genome engineering of heat-tolerance mesophilic and thermophilic bacteria. IMPORTANCE The randomized base pairing in the interval of IBS2 and IBS1 of Tel3c/4c intron (-8 and -7 sites) in Thermotargetron (TMT) results in a low success rate and gene-targeting efficiency in bacteria. In the present work, we constructed a randomized gene-targeting plasmids pool (RGPP) to study whether there is a base preference in target sequences. Among all the successful "retrohoming" targets, we found that a new EBS2b-IBS2b base paring (A-8/T-8) significantly increased TMT's gene-targeting efficiency, and the concept is also applicable to other gene targets in redesigned gene-targeting plasmids pool in E. coli. The improved TMT is a promising tool for the genetic engineering of bacteria and could promote metabolic engineering and synthetic biology research in valuable microbes that recalcitrance for genetic manipulation.

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

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

Unprecedented frequency of mitochondrial introns in colonial bilaterians

Animal mitogenomes are typically devoid of introns. Here, we report the largest number of mitochondrial introns ever recorded from bilaterian animals. Mitochondrial introns were identified for the first time from the phylum Bryozoa. They were found in four species from three families (Order Cheilostomatida). A total of eight introns were found in the complete mitogenome of Exechonella vieirai, and five, 17 and 18 introns were found in the partial mitogenomes of Parantropora penelope, Discoporella cookae and Cupuladria biporosa, respectively. Intron-encoded protein domains reverse transcriptase and intron maturase (RVT-IM) were identified in all species. Introns in E. vieirai and P. penelope had conserved Group II intron ribozyme domains V and VI. Conserved domains were lacking from introns in D. cookae and C. biporosa, preventing their further categorization. Putative origins of metazoan introns were explored in a phylogenetic context, using an up-to-date alignment of mitochondrial RVT-IM domains. Results confirmed previous findings of multiple origins of annelid, placozoan and sponge RVT-IM domains and provided evidence for common intron donor sources across metazoan phyla. Our results corroborate growing evidence that some metazoans with regenerative abilities (i.e. placozoans, sponges, annelids and bryozoans) are susceptible to intron integration, most likely via horizontal gene transfer.

Mitochondrion-encoded circular RNAs are widespread and translatable in plants

Nucleus-encoded circular RNAs (ncircRNAs) have been widely detected in eukaryotes, and most circRNA identification algorithms are designed to identify them. However, using these algorithms, few mitochondrion-encoded circRNAs (mcircRNAs) have been identified in plants, and the role of plant mcircRNAs has not yet been addressed. Here, we developed a circRNA identification algorithm, mitochondrion-encoded circRNA identifier, based on common features of plant mitochondrial genomes. We identified 7,524, 9,819, 1,699, 1,821, 1,809, and 5,133 mcircRNAs in maize (Zea mays), Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), tomato (Solanum lycopersicum), cucumber (Cucumis sativus), and grape (Vitis vinifera), respectively. These mcircRNAs were experimentally validated. Plant mcircRNAs had distinct characteristics from ncircRNAs, and they were more likely to be derived from RNA degradation but not intron backsplicing. Alternative circularization was prevalent in plant mitochondria, and most parental genomic regions hosted multiple mcircRNA isoforms, which have homogenous 5' termini but heterogeneous 3' ends. By analysis of mitopolysome and mitoribosome profiling data, 1,463 mcircRNAs bound to ribosomes were detected in maize and Arabidopsis. Further analysis of mass spectrometry-based proteomics data identified 358 mcircRNA-derived polypeptides. Overall, we developed a computational pipeline that efficiently identifies plant mcircRNAs, and we demonstrated mcircRNAs are widespread and translated in plants.

Group II Intron-Encoded Proteins (IEPs/Maturases) as Key Regulators of Nad1 Expression and Complex I Biogenesis in Land Plant Mitochondria

Mitochondria are semi-autonomous organelles that produce much of the energy required for cellular metabolism. As descendants of a bacterial symbiont, most mitochondria harbor their own genetic system (mtDNA/mitogenome), with intrinsic machineries for transcription and protein translation. A notable feature of plant mitochondria involves the presence of introns (mostly group II-type) that reside in many organellar genes. The splicing of the mtRNAs relies on the activities of various protein cofactors, which may also link organellar functions with cellular or environmental signals. The splicing of canonical group II introns is aided by an ancient class of RT-like enzymes (IEPs/maturases, MATs) that are encoded by the introns themselves and act specifically on their host introns. The plant organellar introns are degenerated in structure and are generally also missing their cognate intron-encoded proteins. The factors required for plant mtRNA processing are mostly nuclearly-encoded, with the exception of a few degenerated MATs. These are in particular pivotal for the maturation of NADH-dehydrogenase transcripts. In the following review we provide an update on the non-canonical MAT factors in angiosperm mitochondria and summarize the current knowledge of their essential roles in regulating Nad1 expression and complex I (CI) biogenesis during embryogenesis and early plant life.

Refining Mitochondrial Intron Classification With ERPIN: Identification Based on Conservation of Sequence Plus Secondary Structure Motifs

Mitochondrial genomes—in particular those of fungi—often encode genes with a large number of Group I and Group II introns that are conserved at both the sequence and the RNA structure level. They provide a rich resource for the investigation of intron and gene structure, self- and protein-guided splicing mechanisms, and intron evolution. Yet, the degree of sequence conservation of introns is limited, and the primary sequence differs considerably among the distinct intron sub-groups. It makes intron identification, classification, structural modeling, and the inference of gene models a most challenging and error-prone task—frequently passed on to an “expert” for manual intervention. To reduce the need for manual curation of intron structures and mitochondrial gene models, computational methods using ERPIN sequence profiles were initially developed in 2007. Here we present a refinement of search models and alignments using the now abundant publicly available fungal mtDNA sequences. In addition, we have tested in how far members of the originally proposed sub-groups are clearly distinguished and validated by our computational approach. We confirm clearly distinct mitochondrial Group I sub-groups IA1, IA3, IB3, IC1, IC2, and ID. Yet, IB1, IB2, and IB4 ERPIN models are overlapping substantially in predictions, and are therefore combined and reported as IB. We have further explored the conversion of our ERPIN profiles into covariance models (CM). Current limitations and prospects of the CM approach will be discussed.

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.

Sensitivity of the U.S. Wheat Powdery Mildew Population to Quinone Outside Inhibitor Fungicides and Determination of the Complete Blumeria graminis f. sp. tritici Cytochrome b Gene

Wheat powdery mildew, caused by Blumeria graminis f. sp. tritici, is managed primarily with cultivar resistance and foliar fungicides. Quinone outside inhibitors (QoIs), which target the mitochondrial cytochrome b (cytb) gene, are one of the two main fungicide classes used on wheat. While European populations of B. graminis f. sp. tritici are widely insensitive to QoIs, largely because of the cytb mutation G143A, the QoI sensitivity of the U.S. B. graminis f. sp. tritici population had never been evaluated despite years of QoI use on U.S. wheat. A total of 381 B. graminis f. sp. tritici isolates from 15 central and eastern U.S. states were screened for sensitivity to QoI fungicides pyraclostrobin and picoxystrobin. A modest range of sensitivities was observed, with maximum resistance factors of 11.2 for pyraclostrobin and 5.3 for picoxystrobin. The F129L, G137R, and G143A cytb mutations were not detected in the U.S. B. graminis f. sp. tritici population, nor were mutations identified in the PEWY loop, a key part of the Qo site. Thus, no genetic basis for the observed quantitative variation in QoI sensitivity of U.S. B. graminis f. sp. tritici was identified. Isolate sporulation was weakly negatively associated with reduced QoI sensitivity, suggesting a fitness cost. In the course of the study, the complete B. graminis f. sp. tritici cytb gene sequence was determined for the first time in the isolate 96224 v. 3.16 reference genome. Contrary to previous reports, the gene has an intron that appears to belong to intron group II, which is unusual in fungi. The study was the first QoI sensitivity screening of a large, geographically diverse set of U.S. B. graminis f. sp. tritici isolates, and while the population as a whole remains relatively sensitive, some quantitative loss of efficacy was observed.

Group II intron and repeat-rich red algal mitochondrial genomes demonstrate the dynamic recent history of autocatalytic RNAs

Background

Group II introns are mobile genetic elements that can insert at specific target sequences, however, their origins are often challenging to reconstruct because of rapid sequence decay following invasion and spread into different sites. To advance understanding of group II intron spread, we studied the intron-rich mitochondrial genome (mitogenome) in the unicellular red alga, Porphyridium.

Results

Analysis of mitogenomes in three closely related species in this genus revealed they were 3–6-fold larger in size (56–132 kbp) than in other red algae, that have genomes of size 21–43 kbp. This discrepancy is explained by two factors, group II intron invasion and expansion of repeated sequences in large intergenic regions. Phylogenetic analysis demonstrates that many mitogenome group II intron families are specific to Porphyridium, whereas others are closely related to sequences in fungi and in the red alga-derived plastids of stramenopiles. Network analysis of intron-encoded proteins (IEPs) shows a clear link between plastid and mitochondrial IEPs in distantly related species, with both groups associated with prokaryotic sequences.

Conclusion

Our analysis of group II introns in Porphyridium mitogenomes demonstrates the dynamic nature of group II intron evolution, strongly supports the lateral movement of group II introns among diverse eukaryotes, and reveals their ability to proliferate, once integrated in mitochondrial DNA.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01200-3.

An overview of RNA splicing and functioning of splicing factors in land plant chloroplasts

RNA splicing refers to a process by which introns of a pre-mRNA are excised and the exons at both ends are joined together. Chloroplast introns are inherently self-splicing ribozymes, but over time, they have lost self-splicing ability due to the degeneration of intronic elements. Thus, the splicing of chloroplast introns relies heavily on nuclear-encoded splicing factors, which belong to diverse protein families. Different splicing factors and their shared intron targets are supposed to form ribonucleoprotein particles (RNPs) to facilitate intron splicing. As characterized in a previous review, around 14 chloroplast intron splicing factors were identified until 2010. However, only a few genetic and biochemical evidence has shown that these splicing factors are required for the splicing of one or several introns. The roles of splicing factors are generally believed to facilitate intron folding; however, the precise role of each protein in RNA splicing remains ambiguous. This may be because the precise binding site of most of these splicing factors remains unexplored. In the last decade, several new splicing factors have been identified. Also, several splicing factors were found to bind to specific sequences within introns, which enhanced the understanding of splicing factors. Here, we summarize recent progress on the splicing factors in land plant chloroplasts and discuss their possible roles in chloroplast RNA splicing based on previous studies.

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.

Identification and classification of reverse transcriptases in bacterial genomes and metagenomes

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

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.

Genic Selection Within Prokaryotic Pangenomes

Understanding the evolutionary forces shaping prokaryotic pangenome structure is a major goal of microbial evolution research. Recent work has highlighted that a substantial proportion of accessory genes appear to confer niche-specific adaptations. This work has primarily focused on selection acting at the level of individual cells. Herein, we discuss a lower level of selection that also contributes to pangenome variation: genic selection. This refers to cases where genetic elements, rather than individual cells, are the entities under selection. The clearest examples of this form of selection are selfish mobile genetic elements, which are those that have either a neutral or a deleterious effect on host fitness. We review the major classes of these and other mobile elements and discuss the characteristic features of such elements that could be under genic selection. We also discuss how genetic elements that are beneficial to hosts can also be under genic selection, a scenario that may be more prevalent but not widely appreciated, because disentangling the effects of selection at different levels (i.e., organisms vs. genes) is challenging. Nonetheless, an appreciation for the potential action and implications of genic selection is important to better understand the evolution of prokaryotic pangenomes.

Transposable elements shape the evolution of mammalian development

Transposable elements (TEs) promote genetic innovation but also threaten genome stability. Despite multiple layers of host defence, TEs actively shape mammalian-specific developmental processes, particularly during pre-implantation and extra-embryonic development and at the maternal-fetal interface. Here, we review how TEs influence mammalian genomes both directly by providing the raw material for genetic change and indirectly via co-evolving TE-binding Krüppel-associated box zinc finger proteins (KRAB-ZFPs). Throughout mammalian evolution, individual activities of ancient TEs were co-opted to enable invasive placentation that characterizes live-born mammals. By contrast, the widespread activity of evolutionarily young TEs may reflect an ongoing co-evolution that continues to impact mammalian development.

Organelle Genome Variation in the Red Algal Genus Ahnfeltia (Florideophyceae)

The agarophyte Ahnfeltia (Ahnfeltiales, Rhodophyta) is a globally widespread genus with 11 accepted species names. Two of the most widespread species in this genus, A. plicata and A. fastigiata, may have diverged genetically due to past geographic changes and subsequent geographic isolation. To investigate this genomic and genetic diversity, we generated new plastid (ptDNAs) and mitochondrial genomes (mtDNAs) of these Ahnfeltia species from four different regions (A. plicata - Chile and UK and A. fastigiata - Korea and Oregon). Two architecture variations were found in the Ahnfeltia genomes: in ptDNA of A. fastigiata Oregon, the hypothetical pseudogene region was translocated, likely due to recombination with palindromic repeats or a gene transfer from a red algal plasmid. In mtDNA of A. fastigiata Korea, the composition of the group II intronic ORFs was distinct from others suggesting different scenarios of gain and loss of group II intronic ORFs. These features resulted in genome size differences between the two species. Overall gene contents of organelle genomes of Ahnfeltia were conserved. Phylogenetic analysis using concatenated genes from ptDNAs and mtDNAs supported the monophyly of the Ahnfeltiophycidae. The most probable individual gene trees showed that the Ahnfeltia populations were genetically diversified. These trees, the cox1 haplotype network, and a dN/dS analysis all supported the theory that these Ahnfeltia populations have diversified genetically in accordance with geographic distribution.

Organelle Genomes and Transcriptomes of Nymphaea Reveal the Interplay between Intron Splicing and RNA Editing

Posttranscriptional modifications, including intron splicing and RNA editing, are common processes during regulation of gene expression in plant organelle genomes. However, the intermediate products of intron-splicing, and the interplay between intron-splicing and RNA-editing were not well studied. Most organelle transcriptome analyses were based on the Illumina short reads which were unable to capture the full spectrum of transcript intermediates within an organelle. To fully investigate the intermediates during intron splicing and the underlying relationships with RNA editing, we used PacBio DNA-seq and Iso-seq, together with Illumina short reads genome and transcriptome sequencing data to assemble the chloroplast and mitochondrial genomes of Nymphaea 'Joey Tomocik' and analyze their posttranscriptional features. With the direct evidence from Iso-seq, multiple intermediates partially or fully intron-spliced were observed, and we also found that both cis- and trans-splicing introns were spliced randomly. Moreover, by using rRNA-depleted and non-Oligo(dT)-enrichment strand-specific RNA-seq data and combining direct SNP-calling and transcript-mapping methods, we identified 98 and 865 RNA-editing sites in the plastome and mitogenome of N. 'Joey Tomocik', respectively. The target codon preference, the tendency of increasing protein hydrophobicity, and the bias distribution of editing sites are similar in both organelles, suggesting their common evolutionary origin and shared editing machinery. The distribution of RNA editing sites also implies that the RNA editing sites in the intron and exon regions may splice synchronously, except those exonic sites adjacent to intron which could only be edited after being intron-spliced. Our study provides solid evidence for the multiple intermediates co-existing during intron-splicing and their interplay with RNA editing in organelle genomes of a basal angiosperm.

Organellar Introns in Fungi, Algae, and Plants

Introns are ubiquitous in eukaryotic genomes and have long been considered as ‘junk RNA’ but the huge energy expenditure in their transcription, removal, and degradation indicate that they may have functional significance and can offer evolutionary advantages. In fungi, plants and algae introns make a significant contribution to the size of the organellar genomes. Organellar introns are classified as catalytic self-splicing introns that can be categorized as either Group I or Group II introns. There are some biases, with Group I introns being more frequently encountered in fungal mitochondrial genomes, whereas among plants Group II introns dominate within the mitochondrial and chloroplast genomes. Organellar introns can encode a variety of proteins, such as maturases, homing endonucleases, reverse transcriptases, and, in some cases, ribosomal proteins, along with other novel open reading frames. Although organellar introns are viewed to be ribozymes, they do interact with various intron- or nuclear genome-encoded protein factors that assist in the intron RNA to fold into competent splicing structures, or facilitate the turn-over of intron RNAs to prevent reverse splicing. Organellar introns are also known to be involved in non-canonical splicing, such as backsplicing and trans-splicing which can result in novel splicing products or, in some instances, compensate for the fragmentation of genes by recombination events. In organellar genomes, Group I and II introns may exist in nested intronic arrangements, such as introns within introns, referred to as twintrons, where splicing of the external intron may be dependent on splicing of the internal intron. These nested or complex introns, with two or three-component intron modules, are being explored as platforms for alternative splicing and their possible function as molecular switches for modulating gene expression which could be potentially applied towards heterologous gene expression. This review explores recent findings on organellar Group I and II introns, focusing on splicing and mobility mechanisms aided by associated intron/nuclear encoded proteins and their potential roles in organellar gene expression and cross talk between nuclear and organellar genomes. Potential application for these types of elements in biotechnology are also discussed.

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.

In and out: Evolution of viral sequences in the mitochondrial genomes of legumes (Fabaceae)

Plant specific mitoviruses in the 'genus' Mitovirus (Narnaviridae) and their integrated sequences (non-retroviral endogenous RNA viral elements or NERVEs) have been recently identified in various plant lineages. However, the sparse phylogenetic coverage of complete plant mitochondrial genome (mitogenome) sequences and the non-conserved nature of mitochondrial intergenic regions have hindered comparative studies on mitovirus NERVEs in plants. In this study, 10 new mitogenomes were sequenced from legumes (Fabaceae). Based on comparative genomic analysis of 27 total mitogenomes, we identified mitovirus NERVEs and transposable elements across the family. All legume mitogenomes included NERVEs and total NERVE length varied from ca. 2 kb in the papilionoid Trifolium to 35 kb in the mimosoid Acacia. Most of the NERVE integration sites were in highly variable intergenic regions, however, some were positioned in six cis-spliced mitochondrial introns. In the Acacia mitogenome, there were L1-like transposon sequences including an almost full-length copy with target site duplications (TSDs). The integration sites of NERVEs in four introns showed evidence of L1-like retrotransposition events. Phylogenetic analysis revealed that there were multiple instances of precise deletion of NERVEs between TSDs. This study provides clear evidence that a L1-like retrotransposition mechanism has a long history of contributing to the integration of viral RNA into plant mitogenomes while microhomology-mediated deletion can restore the integration site.

nMAT3 is an essential maturase splicing factor required for holo-complex I biogenesis and embryo development in Arabidopsis thaliana plants

Group-II introns are self-splicing mobile genetic elements consisting of catalytic intron-RNA and its related intron-encoded splicing maturase protein cofactor. Group-II sequences are particularly plentiful within the mitochondria of land plants, where they reside within many critical gene loci. During evolution, the plant organellar introns have degenerated, such as they lack regions that are are required for splicing, and also lost their evolutionary related maturase proteins. Instead, for their splicing the organellar introns in plants rely on different host-acting protein cofactors, which may also provide a means to link cellular signals with respiratory functions. The nuclear genome of Arabidopsis thaliana encodes four maturase-related factors. Previously, we showed that three of the maturases, nMAT1, nMAT2 and nMAT4, function in the excision of different group-II introns in Arabidopsis mitochondria. The function of nMAT3 (encoded by the At5g04050 gene locus) was found to be essential during early embryogenesis. Using a modified embryo-rescue method, we show that nMAT3-knockout plants are strongly affected in the splicing of nad1 introns 1, 3 and 4 in Arabidopsis mitochondria, resulting in complex-I biogenesis defects and altered respiratory activities. Functional complementation of nMAT3 restored the organellar defects and embryo-arrested phenotypes associated with the nmat3 mutant line. Notably, nMAT3 and nMA4 were found to act on the same RNA targets but have no redundant functions in the splicing of nad1 transcripts. The two maturases, nMAT3 and nMAT4 are likely to cooperate together in the maturation of nad1 pre-RNAs. Our results provide important insights into the roles of maturases in mitochondria gene expression and the biogenesis of the respiratory system during early plant life.

Identification of polycistronic transcriptional units and non-canonical introns in green algal chloroplasts based on long-read RNA sequencing data

BACKGROUND

Chloroplasts are important semi-autonomous organelles in plants and algae. Unlike higher plants, the chloroplast genomes of green algal linage have distinct features both in organization and expression. Despite the architecture of chloroplast genome having been extensively studied in higher plants and several model species of algae, little is known about the transcriptional features of green algal chloroplast-encoded genes.

RESULTS

Based on full-length cDNA (Iso-Seq) sequencing, we identified widely co-transcribed polycistronic transcriptional units (PTUs) in the green alga Caulerpa lentillifera. In addition to clusters of genes from the same pathway, we identified a series of PTUs of up to nine genes whose function in the plastid is not understood. The RNA data further allowed us to confirm widespread expression of fragmented genes and conserved open reading frames, which are both important features in green algal chloroplast genomes. In addition, a newly fragmented gene specific to C. lentillifera was discovered, which may represent a recent gene fragmentation event in the chloroplast genome. With the newly annotated exon-intron boundary information, gene structural annotation was greatly improved across the siphonous green algae lineages. Our data also revealed a type of non-canonical Group II introns, with a deviant secondary structure and intronic ORFs lacking known splicing or mobility domains. These widespread introns have conserved positions in their genes and are excised precisely despite lacking clear consensus intron boundaries.

CONCLUSION

Our study fills important knowledge gaps in chloroplast genome organization and transcription in green algae, and provides new insights into expression of polycistronic transcripts, freestanding ORFs and fragmented genes in algal chloroplast genomes. Moreover, we revealed an unusual type of Group II intron with distinct features and conserved positions in Bryopsidales. Our data represents interesting additions to knowledge of chloroplast intron structure and highlights clusters of uncharacterized genes that probably play important roles in plastids.

Screening for potential interaction partners with surface plasmon resonance imaging coupled to MALDI mass spectrometry

We coupled SPR imaging (SPRi) with matrix-assisted laser desorption/ionization mass spectrometry (MALDI MS) to identify new potential RNA binders. Here, we improve this powerful method, especially by optimizing the proteolytic digestion (type of reducing agent, its concentration, and incubation time), to work with complex mixtures, specifically a lysate of the rough mitochondrial fraction from yeast. The advantages of this hyphenated method compared to column-based or separate analyses are (i) rapid and direct visual readout from the SPRi array, (ii) possibility of high-throughput analysis of different interactions in parallel, (iii) high sensitivity, and (iv) no sample loss or contamination due to elution or micro-recovery procedures. The model system used is a catalytically active RNA (group IIB intron from Saccharomyces cerevisiae, Sc.ai5γ) and its cofactor Mss116. The protein supports the RNA folding process and thereby the subsequent excision of the intronic RNA from the coding part. Using the novel approach of coupling SPR with MALDI MS, we report the identification of potential RNA-binding proteins from a crude yeast mitochondrial lysate in a non-targeted approach. Our results show that proteins other than the well-known cofactor Mss116 interact with Sc.ai5γ (Dbp8, Prp8, Mrp13, and Cullin-3), suggesting that the intron folding and splicing are regulated by more than one cofactor in vivo.

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

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

Cloning and Detection of Aptamer-Ribozyme Conjugations

RNA aptamers can be used to target proteins or nucleic acids for therapeutic purposes and are candidates for RNA-mediated gene therapy. Like other small therapeutic RNAs, they can be expressed in cells from DNA templates that include a cellular promoter upstream of the RNA coding sequence. Secondary structures flanking aptamers can be used to enhance the activity or stability of these molecules. Notably, flanking self-cleaving ribozymes to remove extraneous nucleotides included during transcription as well as flanking hairpins to improve RNA stability have been used to increase the effect of therapeutic aptamers. Here we describe the cloning procedure of aptamers containing different flanking secondary structures and methods to compare their expression levels by a northern blot protocol optimized for the detection of small RNA molecules.

Methylation of rRNA as a host defense against rampant group II intron retrotransposition

BACKGROUND

Group II introns are mobile retroelements, capable of invading new sites in DNA. They are self-splicing ribozymes that complex with an intron-encoded protein to form a ribonucleoprotein that targets DNA after splicing. These molecules can invade DNA site-specifically, through a process known as retrohoming, or can invade ectopic sites through retrotransposition. Retrotransposition, in particular, can be strongly influenced by both environmental and cellular factors.

RESULTS

To investigate host factors that influence retrotransposition, we performed random insertional mutagenesis using the ISS1 transposon to generate a library of over 1000 mutants in Lactococcus lactis, the native host of the Ll.LtrB group II intron. By screening this library, we identified 92 mutants with increased retrotransposition frequencies (RTP-ups). We found that mutations in amino acid transport and metabolism tended to have increased retrotransposition frequencies. We further explored a subset of these RTP-up mutants, the most striking of which is a mutant in the ribosomal RNA methyltransferase rlmH, which exhibited a reproducible 20-fold increase in retrotransposition frequency. In vitro and in vivo experiments revealed that ribosomes in the rlmH mutant were defective in the m3Ψ modification and exhibited reduced binding to the intron RNA.

CONCLUSIONS

Taken together, our results reinforce the importance of the native host organism in regulating group II intron retrotransposition. In particular, the evidence from the rlmH mutant suggests a role for ribosome modification in limiting rampant retrotransposition.

The Mitogenomes of Ophiostoma minus and Ophiostoma piliferum and Comparisons With Other Members of the Ophiostomatales

Fungi assigned to the Ophiostomatales are of economic concern as many are blue-stain fungi and some are plant pathogens. The mitogenomes of two blue-stain fungi, Ophiostoma minus and Ophiostoma piliferum, were sequenced and compared with currently available mitogenomes for other members of the Ophiostomatales. Species representing various genera within the Ophiostomatales have been examined for gene content, gene order, phylogenetic relationships, and the distribution of mobile elements. Gene synteny is conserved among the Ophiostomatales but some members were missing the atp9 gene. A genome wide intron landscape has been prepared to demonstrate the distribution of the mobile genetic elements (group I and II introns and homing endonucleases) and to provide insight into the evolutionary dynamics of introns among members of this group of fungi. Examples of complex introns or nested introns composed of two or three intron modules have been observed in some species. The size variation among the mitogenomes (from 23.7 kb to about 150 kb) is mostly due to the presence and absence of introns. Members of the genus Sporothrix sensu stricto appear to have the smallest mitogenomes due to loss of introns. The taxonomy of the Ophiostomatales has recently undergone considerable revisions; however, some lineages remain unresolved. The data showed that genera such as Raffaelea appear to be polyphyletic and the separation of Sporothrix sensu stricto from Ophiostoma is justified.