An intron in the nuclear ribosomal DNA of Didymium iridis codes for a group I ribozyme and a novel ribozyme that cooperate in self-splicing

Mobile group I introns at nuclear rDNA position L2066 harbor sense and antisense homing endonuclease genes intervened by spliceosomal introns

Background

Mobile group I introns encode homing endonucleases that confer intron mobility initiated by a double-strand break in the intron-lacking allele at the site of insertion. Nuclear ribosomal DNA of some fungi and protists contain mobile group I introns harboring His-Cys homing endonuclease genes (HEGs). An intriguing question is how protein-coding genes embedded in nuclear ribosomal DNA become expressed. To address this gap of knowledge we analyzed nuclear L2066 group I introns from myxomycetes and ascomycetes.

Results

A total of 34 introns were investigated, including two identified mobile-type introns in myxomycetes with HEGs oriented in sense or antisense directions. Intriguingly, both HEGs are interrupted by spliceosomal introns. The intron in Didymium squamulosum, which harbors an antisense oriented HEG, was investigated in more detail. The group I intron RNA self-splices in vitro, thus generating ligated exons and full-length intron circles. The intron HEG is expressed in vivo in Didymium cells, which involves removal of a 47-nt spliceosomal intron (I-47) and 3′ polyadenylation of the mRNA. The D. squamulosum HEG (lacking the I-47 intron) was over-expressed in E. coli, and the corresponding protein was purified and shown to confer endonuclease activity. The homing endonuclease was shown to cleave an intron-lacking DNA and to produce a pentanucleotide 3′ overhang at the intron insertion site.

Conclusions

The L2066 family of nuclear group I introns all belong to the group IE subclass. The D. squamulosum L2066 intron contains major hallmarks of a true mobile group I intron by encoding a His-Cys homing endonuclease that generates a double-strand break at the DNA insertion site. We propose a potential model to explain how an antisense HEG becomes expressed from a nuclear ribosomal DNA locus.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-022-00280-4.

Targeted editing and evolution of engineered ribosomes in vivo by filtered editing

Genome editing technologies introduce targeted chromosomal modifications in organisms yet are constrained by the inability to selectively modify repetitive genetic elements. Here we describe filtered editing, a genome editing method that embeds group 1 self-splicing introns into repetitive genetic elements to construct unique genetic addresses that can be selectively modified. We introduce intron-containing ribosomes into the E. coli genome and perform targeted modifications of these ribosomes using CRISPR/Cas9 and multiplex automated genome engineering. Self-splicing of introns post-transcription yields scarless RNA molecules, generating a complex library of targeted combinatorial variants. We use filtered editing to co-evolve the 16S rRNA to tune the ribosome's translational efficiency and the 23S rRNA to isolate antibiotic-resistant ribosome variants without interfering with native translation. This work sets the stage to engineer mutant ribosomes that polymerize abiological monomers with diverse chemistries and expands the scope of genome engineering for precise editing and evolution of repetitive DNA sequences.

A Phylogenetic Approach to Structural Variation in Organization of Nuclear Group I Introns and Their Ribozymes

Nuclear group I introns are restricted to the ribosomal DNA locus where they interrupt genes for small subunit and large subunit ribosomal RNAs at conserved sites in some eukaryotic microorganisms. Here, the myxomycete protists are a frequent source of nuclear group I introns due to their unique life strategy and a billion years of separate evolution. The ribosomal DNA of the myxomycete Mucilago crustacea was investigated and found to contain seven group I introns, including a direct repeat-containing intron at insertion site S1389 in the small subunit ribosomal RNA gene. We collected, analyzed, and compared 72 S1389 group IC1 introns representing diverse myxomycete taxa. The consensus secondary structure revealed a conserved ribozyme core, but with surprising sequence variations in the guanosine binding site in segment P7. Some S1389 introns harbored large extension sequences in the peripheral region of segment P9 containing direct repeat arrays. These repeats contained up to 52 copies of a putative internal guide sequence motif. Other S1389 introns harbored homing endonuclease genes in segment P1 encoding His-Cys proteins. Homing endonuclease genes were further interrupted by small spliceosomal introns that have to be removed in order to generate the open reading frames. Phylogenetic analyses of S1389 intron and host gene indicated both vertical and horizontal intron transfer during evolution, and revealed sporadic appearances of direct repeats, homing endonuclease genes, and guanosine binding site variants among the myxomycete taxa.

Use of a Lariat Capping Ribozyme to Study Cap Function In Vivo

A lariat cap is a naturally occurring substitute of a conventional mRNA cap and is found in a particular genomic setting in a few eukaryotic microorganisms. It is installed by the lariat capping ribozyme acting in cis. In principle, any RNA molecule in any organism can be equipped with a lariat cap in vivo when expressed downstream of a lariat capping ribozyme. Lariat capping is thus a versatile tool for studying the importance of the 5' end structure of RNA molecules. In this chapter, we present protocols to validate the presence of the lariat cap and measure the efficiency of in vivo cleavage by the lariat capping ribozyme.

Conformational adaptation of UNCG loops upon crowding

If the A-form helix is the major structural motif found in RNA, the loops that cap them constitute the second most important family of motifs. Among those, two are overrepresented, GNRA and UNCG tetraloops. Recent surveys of RNA structures deposited in the PDB show that GNRA and UNCG tetraloops can adopt tertiary folds that are very different from their canonical conformations, characterized by the presence of a U-turn of a Z-turn, respectively. Crystallographic data from both a lariat-capping (LC) ribozyme and a group II intron ribozyme reveal that a given UUCG tetraloop can adopt a distinct fold depending on its structural environment. Specifically, when the crystal packing applies relaxed constraints on the loop, the canonical Z-turn conformation is observed. In contrast, a highly packed environment induces "squashing" of the tetraloop by distorting its sugar-phosphate backbone in a specific way that expels the first and fourth nucleobases out of the loop, and falls in van der Waals distance of the last base pair of the helix, taking the place of the pair formed between the first and fourth residues in Z-turn loops. The biological relevance of our observations is supported by the presence of similarly deformed loops in the highly packed environment of the ribosome and in a complex between a dsRNA and a RNase III. The finding that Z-turn loops change conformation under higher molecular packing suggests that, in addition to their demonstrated role in stabilizing RNA folding, they may contribute to the three-dimensional structure of RNA by mediating tertiary interactions with distal residues.

Integrated evolution of ribosomal RNAs, introns, and intron nurseries

The initial components of ribosomes first appeared more than 3.8 billion years ago during a time when many types of RNAs were evolving. While modern ribosomes are complex molecular machines consisting of rRNAs and proteins, they were assembled during early evolution by the association and joining of small functional RNA units. Introns may have provided the means to ligate many of these pieces together. All four classes of introns (group I, group II, spliceosomal, and archaeal) are present in many rRNA gene loci over a broad phylogenetic range. A survey of rRNA intron sequences across the three major life domains suggests that some of the classes of introns may have diverged from one another within rRNA gene loci. Analyses of rRNA sequences revealed self-splicing group I and group II introns are present in ancestral regions of the SSU (small subunit) and LSU (large subunit), whereas spliceosomal and archaeal introns appeared in sections of the rRNA that evolved later. Most classes of introns increased in number for approximately 1 billion years. However, their frequencies are low in the most recently evolved regions added to the SSU and LSU rRNAs. Furthermore, many of the introns appear to have been in the same locations for billions of years, suggesting an ancient origin for these sequences. In this Perspectives paper, I reviewed and analyzed rRNA intron sequences, locations, structural characteristics, and splicing mechanisms; and suggest that rRNA gene loci may have served as evolutionary nurseries for intron formation and diversification.

Accumulation of Stable Full-Length Circular Group I Intron RNAs during Heat-Shock

Group I introns in nuclear ribosomal RNA of eukaryotic microorganisms are processed by splicing or circularization. The latter results in formation of full-length circular introns without ligation of the exons and has been proposed to be active in intron mobility. We applied qRT-PCR to estimate the copy number of circular intron RNA from the myxomycete Didymium iridis. In exponentially growing amoebae, the circular introns are nuclear and found in 70 copies per cell. During heat-shock, the circular form is up-regulated to more than 500 copies per cell. The intron harbours two ribozymes that have the potential to linearize the circle. To understand the structural features that maintain circle integrity, we performed chemical and enzymatic probing of the splicing ribozyme combined with molecular modeling to arrive at models of the inactive circular form and its active linear counterpart. We show that the two forms have the same overall structure but differ in key parts, including the catalytic core element P7 and the junctions at which reactions take place. These differences explain the relative stability of the circular species, demonstrate how it is prone to react with a target molecule for circle integration and thus supports the notion that the circular form is a biologically significant molecule possibly with a role in intron mobility.

Convergent evolution of twintron-like configurations: One is never enough

Introns inserted within introns are commonly referred to as twintrons, however the original definition for twintron implied that splicing of the external member of the twintron could only proceed upon splicing of the internal member. This review examines the various types of twintron-like arrangements that have been reported and assigns them to either nested or twintron categories that are subdivided further into subtypes based on differences of their mode of splicing. Twintron-like arrangements evolved independently by fortuitous events among different types of introns but once formed they offer opportunities for the evolution of new regulatory strategies and/or novel genetic elements.

RNA circularization strategies in vivo and in vitro

In the plenitude of naturally occurring RNAs, circular RNAs (circRNAs) and their biological role were underestimated for years. However, circRNAs are ubiquitous in all domains of life, including eukaryotes, archaea, bacteria and viruses, where they can fulfill diverse biological functions. Some of those functions, as for example playing a role in the life cycle of viral and viroid genomes or in the maturation of tRNA genes, have been elucidated; other putative functions still remain elusive. Due to the resistance to exonucleases, circRNAs are promising tools for in vivo application as aptamers, trans-cleaving ribozymes or siRNAs. How are circRNAs generated in vivo and what approaches do exist to produce ring-shaped RNAs in vitro? In this review we illustrate the occurrence and mechanisms of RNA circularization in vivo, survey methods for the generation of circRNA in vitro and provide appropriate protocols.

Molecular characterization of a new member of the lariat capping twin-ribozyme introns

BACKGROUND

Twin-ribozyme introns represent a complex class of mobile group I introns that harbour a lariat capping (LC) ribozyme and a homing endonuclease gene embedded in a conventional self-splicing group I ribozyme (GIR2). Twin-ribozyme introns have so far been confined to nucleolar DNA in Naegleria amoeboflagellates and the myxomycete Didymium iridis.

RESULTS

We characterize structural organization, catalytic properties and molecular evolution of a new twin-ribozyme intron in Allovahlkampfia (Heterolobosea). The intron contains two ribozyme domains with different functions in ribosomal RNA splicing and homing endonuclease mRNA maturation. We found Allovahlkampfia GIR2 to be a typical group IC1 splicing ribozyme responsible for addition of the exogenous guanosine cofactor (exoG), exon ligation and circularization of intron RNA. The Allovahlkampfia LC ribozyme, by contrast, represents an efficient self-cleaving ribozyme that generates a small 2',5' lariat cap at the 5' end of the homing endonuclease mRNA, and thus contributes to intron mobility.

CONCLUSIONS

The discovery of a twin-ribozyme intron in a member of Heterolobosea expands the distribution pattern of LC ribozymes. We identify a putative regulatory RNA element (AP2.1) in the Allovahlkampfia LC ribozyme that involves homing endonuclease mRNA coding sequences as an important structural component.

Speciation of a group I intron into a lariat capping ribozyme

The lariat-capping (LC) ribozyme is a natural ribozyme isolated from eukaryotic microorganisms. Despite apparent structural similarity to group I introns, the LC ribozyme catalyzes cleavage by a 2',5' branching reaction, leaving the 3' product with a 3-nt lariat cap that functionally substitutes for a conventional mRNA cap in the downstream pre-mRNA encoding a homing endonuclease. We describe the crystal structures of the precleavage and postcleavage LC ribozymes, which suggest that structural features inherited from group I ribozymes have undergone speciation due to profound changes in molecular selection pressure, ultimately giving rise to an original branching ribozyme family. The structures elucidate the role of key elements that regulate the activity of the LC ribozyme by conformational switching and suggest a mechanism by which the signal for branching is transmitted to the catalytic core. The structures also show how conserved interactions twist residues, forming the lariat to join chemical groups involved in branching.

Gene regulation by structured mRNA elements

The precise temporal and spatial coordination of gene activity, based on the integration of internal and external signals, is crucial for the accurate functioning of all biological processes. Although the basic principles of gene expression were established some 60 years ago, recent research has revealed a surprising complexity in the control of gene activity. Many of these gene regulatory mechanisms occur at the level of the mRNA, including sophisticated gene control tasks mediated by structured mRNA elements. We now know that mRNA folds can serve as highly specific receptors for various types of molecules, as exemplified by metabolite-binding riboswitches, and interfere with pro- and eukaryotic gene expression at the level of transcription, translation, and RNA processing. Gene regulation by structured mRNA elements comprises versatile strategies including self-cleaving ribozymes, RNA-folding-mediated occlusion or presentation of cis-regulatory sequences, and sequestration of trans-acting factors including other RNAs and proteins.

Nuclear group I introns in self-splicing and beyond

Group I introns are a distinct class of RNA self-splicing introns with an ancient origin. All known group I introns present in eukaryote nuclei interrupt functional ribosomal RNA genes located in ribosomal DNA loci. The discovery of the Tetrahymena intron more than 30 years ago has been essential to our understanding of group I intron catalysis, higher-order RNA structure, and RNA folding, but other intron models have provided information about the biological role. Nuclear group I introns appear widespread among eukaryotic microorganisms, and the plasmodial slime molds (myxomycetes) contain an abundance of self-splicing introns. Here, we summarize the main conclusions from previous work on the Tetrahymena intron on RNA self-splicing catalysis as well as more recent work on myxomycete intron biology. Group I introns in myxomycetes that represent different evolutionary stages, biological roles, and functional settings are discussed.

Homing endonucleases: from basics to therapeutic applications

Homing endonucleases (HE) are double-stranded DNAses that target large recognition sites (12-40 bp). HE-encoding sequences are usually embedded in either introns or inteins. Their recognition sites are extremely rare, with none or only a few of these sites present in a mammalian-sized genome. However, these enzymes, unlike standard restriction endonucleases, tolerate some sequence degeneracy within their recognition sequence. Several members of this enzyme family have been used as templates to engineer tools to cleave DNA sequences that differ from their original wild-type targets. These custom HEs can be used to stimulate double-strand break homologous recombination in cells, to induce the repair of defective genes with very low toxicity levels. The use of tailored HEs opens up new possibilities for gene therapy in patients with monogenic diseases that can be treated ex vivo. This review provides an overview of recent advances in this field.

A group I self-splicing intron in the flagellin gene of the thermophilic bacterium Geobacillus stearothermophilus

A group I intron that can be spliced in vivo and in vitro was identified in the flagellin gene of the thermophilic bacterium Geobacillus stearothermophilus. We also found one or two intervening sequences (IVS) of flagellin genes in five additional bacterial species. Furthermore, we report the presence of these sequences in two sites of a highly conserved region in the flagellin gene.

Natural and unnatural ribozymes: back to the primordial RNA world

We review natural and in vitro selected ribozymes, for which combined studies could provide us with both insight into the functions performed by ancient RNA molecules in a primitive RNA world and a hypothesis about evolutionary steps that led to the contemporary world.

A conformational switch in the DiGIR1 ribozyme involved in release and folding of the downstream I-DirI mRNA

DiGIR1 is a group I-like cleavage ribozyme found as a structural domain within a nuclear twin-ribozyme group I intron. DiGIR1 catalyzes cleavage by branching at an Internal Processing Site (IPS) leading to formation of a lariat cap at the 5'-end of the 3'-cleavage product. The 3'-cleavage product is subsequently processed into an mRNA encoding a homing endonuclease. By analysis of combinations of 5'- and 3'-deletions, we identify a hairpin in the 5'-UTR of the mRNA (HEG P1) that is formed by conformational switching following cleavage. The formation of HEG P1 inhibits the reversal of the branching reaction, thus giving it directionality. Furthermore, the release of the mRNA is a consequence of branching rather than hydrolytic cleavage. A model is put forward that explains the release of the I-DirI mRNA with a lariat cap and a structured 5'-UTR as a direct consequence of the DiGIR1 branching reaction. The role of HEG P1 in GIR1 branching is reminiscent of that of hairpin P-1 in splicing of the Tetrahymena rRNA group I intron and illustrates a general principle in RNA-directed RNA processing.

Recognition of a common rDNA target site in archaea and eukarya by analogous LAGLIDADG and His-Cys box homing endonucleases

The presence of a homing endonuclease gene (HEG) within a microbial intron or intein empowers the entire element with the ability to invade genomic targets. The persistence of a homing endonuclease lineage depends in part on conservation of its DNA target site. One such rDNA sequence has been invaded both in archaea and in eukarya, by LAGLIDADG and His–Cys box homing endonucleases, respectively. The bases encoded by this target include a universally conserved ribosomal structure, termed helix 69 (H69) in the large ribosomal subunit. This region forms the ‘B2a’ intersubunit bridge to the small ribosomal subunit, contacts bound tRNA in the A- and P-sites, and acts as a trigger for ribosome disassembly through its interactions with ribosome recycling factor. We have determined the DNA-bound structure and specificity profile of an archaeal LAGLIDADG homing endonuclease (I-Vdi141I) that recognizes this target site, and compared its specificity with the analogous eukaryal His–Cys box endonuclease I-PpoI. These homodimeric endonuclease scaffolds have arrived at similar specificity profiles across their common biological target and analogous solutions to the problem of accommodating conserved asymmetries within the DNA sequence, but with differences at individual base pairs that are fine-tuned to the sequence conservation of archaeal versus eukaryal ribosomes.

Toward predicting self-splicing and protein-facilitated splicing of group I introns

In the current era of massive discoveries of noncoding RNAs within genomes, being able to infer a function from a nucleotide sequence is of paramount interest. Although studies of individual group I introns have identified self-splicing and nonself-splicing examples, there is no overall understanding of the prevalence of self-splicing or the factors that determine it among the >2300 group I introns sequenced to date. Here, the self-splicing activities of 12 group I introns from various organisms were assayed under six reaction conditions that had been shown previously to promote RNA catalysis for different RNAs. Besides revealing that assessing self-splicing under only one condition can be misleading, this survey emphasizes that in vitro self-splicing efficiency is correlated with the GC content of the intron (>35% GC was generally conductive to self-splicing), and with the ability of the introns to form particular tertiary interactions. Addition of the Neurospora crassa CYT-18 protein activated splicing of two nonself-splicing introns, but inhibited the second step of self-splicing for two others. Together, correlations between sequence, predicted structure and splicing begin to establish rules that should facilitate our ability to predict the self-splicing activity of any group I intron from its sequence.

Splicing and evolution of an unusually small group I intron

Introns are common in the rRNA gene loci of fungal genomes, but biochemical studies to investigate splicing are rare. Here, self-splicing of a very small (67 nucleotide) group I intron is demonstrated. The PaSSU intron (located within the rRNA small subunit gene of Phialophora americana) splices in vitro under group I intron conditions. Most group I ribozymes contain pairing regions P1-P10, with a conserved G.U pair at the 5' splice site, and a G at the 3' intron border. The PaSSU intron contains only P1, P7, and P10. While it contains the G.U pair at the 5' splice, a U is found at the 3' end of the intron instead of a G. Phylogenetic analysis places it within subgroup IC1, whose members are found in the nuclear rRNA genes of fungi. The structural elements are similar to those in the centermost regions of other group I introns. Its size can be explained by a single large deletion that removed P2 through much of P9. Part of the original P9 region has assumed the function of P7. Its small size and genealogy makes it an excellent model to study RNA catalysis and evolution.

Expression of protein-coding genes embedded in ribosomal DNA

Ribosomal DNA (rDNA) is a specialised chromosomal location that is dedicated to high-level transcription of ribosomal RNA genes. Interestingly, rDNAs are frequently interrupted by parasitic elements, some of which carry protein genes. These are non-LTR retrotransposons and group II introns that encode reverse transcriptase-like genes, and group I introns and archaeal introns that encode homing endonuclease genes (HEGs). Although rDNA-embedded protein genes are widespread in nuclei, organelles and bacteria, there is surprisingly little information available on how these genes are expressed. Exceptions include a handful of HEGs from group I introns. Recent studies have revealed unusual and essential roles of group I and group I-like ribozymes in the endogenous expression of HEGs. Here we discuss general aspects of rDNA-embedded protein genes and focus on HEG expression from group I introns in the nucleolus.

Patterns of group I intron presence in nuclear SSU rDNA of the Lichen family Parmeliaceae

Group I introns are commonly reported within nuclear SSU ribosomal DNA of eukaryotic micro-organisms, especially in lichen-forming fungi. We have studied the primary and secondary structure of 70 new nuclear SSU rDNA group I introns of Parmeliaceae (Ascomycota: Lecanorales) and compared them with those available in databases, covering more than 60 species. The analyzed samples of Parmeliaceae fell into two groups, one having an intron at the 1506 site and another lacking this one but having another at the 1516 or 1521 position. Introns at the 1521 position seem to be transposed from 1516 sites. Introns at the 1516 position were similar in structure to ones previously reported at this site and known from other lecanoralean fungi, while those at the 1506 position showed structural differences and no similar introns are known from related fungi. The study of the distribution of group I introns within a large monophyletic ensemble of fungi has revealed an unexpected correlation between intron types and ecological and geographical parameters. The introns at the 1516 position occurred in mainly arctic, boreal, and temperate lichens, while those at position 1506 were present in mainly tropical and subtropical to oceanic mild-temperate taxa. Further, the 1516 introns occurred in genera with few distributed species that could represent older taxa, while the 1506 ones were mainly in species-rich genera that could be of recent speciation, as many species have wide distribution areas. The transition between two different environments has been accompanied by a change in introns gained and lost.

A new RNA branching activity: the GIR1 ribozyme

The formation of lariat intermediates during the first step of splicing of group II introns and spliceosomal introns is a well-studied fundamental reaction in molecular biology. Apart from this prominent example, there are surprisingly few occurrences of branched nucleotides or even 2',5'-phosphodiester bonds in biology. We recently described a new ribozyme, the GIR1 branching ribozyme, which catalyzes the formation of a tiny lariat that caps an mRNA. This new example together with work on artificial branching ribozymes and deoxyribozymes shows that branching is facile and points to the possibility that branching reactions could be more prevalent than previously recognized.

Rapid PCR-based method for detection and differentiation of Didymiaceae and Physaraceae (myxomycetes) in environmental samples

Ecological studies of myxomycetes have been limited by the absence of universal cultivation techniques and the lack of life stage independent identification methods. We designed a novel PCR primer pair for the specific amplification of small subunit ribosomal RNA gene of Didymiaceae and Physaraceae. The primers produced amplicons from 192 fruiting body samples belonging to 10 genera. Twenty-four samples yielded longer fragments and sequence analysis revealed the presence of intron(s). As for the exonic regions, while sequence heterogeneities within a single species/varietas/forma were frequently observed, identical sequences were obtained only from identical species/varietas. The effectiveness of this primer pair in the analysis of morphologically unidentifiable samples was confirmed with the applications to samples of environmental plasmodium/sclerotium and soil. Denaturing gradient gel electrophoresis analysis was also tested with the soil samples. The results presented here demonstrate this PCR-based method can facilitate further ecological studies of Physaraceae and Didymiaceae in the environment.

An mRNA is capped by a 2', 5' lariat catalyzed by a group I-like ribozyme

Twin-ribozyme introns are formed by two ribozymes belonging to the group I family and occur in some ribosomal RNA transcripts. The group I-like ribozyme, GIR1, liberates the 5' end of a homing endonuclease messenger RNA in the slime mold Didymium iridis. We demonstrate that this cleavage occurs by a transesterification reaction with the joining of the first and the third nucleotide of the messenger by a 2',5'-phosphodiester linkage. Thus, a group I-like ribozyme catalyzes an RNA branching reaction similar to the first step of splicing in group II introns and spliceosomal introns. The resulting short lariat, by forming a protective 5' cap, might have been useful in a primitive RNA world.

Capping by branching: a new ribozyme makes tiny lariats

The number of naturally occurring RNA enzymes has just been expanded by the discovery of a new branching ribozyme. But this ribozyme has unexpected relatives: group I introns.

Reverse splicing of a mobile twin-ribozyme group I intron into the natural small subunit rRNA insertion site

A mobile group I intron containing two ribozyme domains and a homing endonuclease gene (twin-ribozyme intron organization) can integrate by reverse splicing into the small subunit rRNA of bacteria and yeast. The integration is sequence-specific and corresponds to the natural insertion site (homing site) of the intron. The reverse splicing is independent of the homing endonuclease gene, but is dependent on the group I splicing ribozyme domain. The observed distribution of group I introns in nature can be explained by horizontal transfer between natural homing sites by reverse splicing and subsequent spread in populations by endonuclease-dependent homing.

Higher-Order Phylogeny of Plasmodial Slime Molds (Myxogastria) Based on Elongation Factor 1-A and Small Subunit rRNA Gene Sequences

The Myxogastria are common soil microorganisms with a life cycle comprised of a plasmodial trophic stage and large fruiting bodies generally visible with the unaided eye. Until now, their classification has been based exclusively on a combination of morphological, ultrastructural, and developmental characters. Our study is the first attempt to examine phylogenetic relationships among these taxa using molecular data. Partial small-subunit ribosomal RNA and/or elongation factor 1-alpha gene sequences were obtained from eleven, mostly field-collected species representing the five orders of Myxogastria. Nineteen sequences were obtained and subjected to phylogenetic analysis together with 10 sequences available from GenBank. Separate and combined analyses of the two data sets support the division of Myxogastria into three distinct groups. The most basal clade consists of the Echinosteliales, an order considered to have affinities with Protostelia. The three species examined possess unpigmented or slightly pigmented spores. The second group consists of Liceales and Trichiales, taxa characterized by the presence of clear, but pigmented, spores. The third group consists of the two remaining orders, Physarales and Stemonitales, both possessing dark spores. This suggests that spore pigmentation is an evolutionarily conservative character in myxogastrians, and that the simple morphology of echinostelids is not a derived feature.

The recent transfer of a homing endonuclease gene

The myxomycete Didymium iridis (isolate Panama 2) contains a mobile group I intron named Dir.S956-1 after position 956 in the nuclear small subunit (SSU) rRNA gene. The intron is efficiently spread through homing by the intron-encoded homing endonuclease I-DirI. Homing endonuclease genes (HEGs) usually spread with their associated introns as a unit, but infrequently also spread independent of introns (or inteins). Clear examples of HEG mobility are however sparse. Here, we provide evidence for the transfer of a HEG into a group I intron named Dir.S956-2 that is inserted into the SSU rDNA of the Costa Rica 8 isolate of D.iridis. Similarities between intron sequences that flank the HEG and rDNA sequences that flank the intron (the homing endonuclease recognition sequence) suggest that the HEG invaded the intron during the recent evolution in a homing-like event. Dir.S956-2 is inserted into the same SSU site as Dir.S956-1. Remarkably, the two group I introns encode distantly related splicing ribozymes with phylogenetically related HEGs inserted on the opposite strands of different peripheral loop regions. The HEGs are both interrupted by small spliceosomal introns that must be removed during RNA maturation.

Site-specific reverse splicing of a HEG-containing group I intron in ribosomal RNA

The wide, but scattered distribution of group I introns in nature is a result of two processes; the vertical inheritance of introns with or without losses, and the occasional transfer of introns across species barriers. Reversal of the group I intron self-splicing reaction, termed reverse splicing, coupled with reverse transcription and genomic integration potentially mediate an RNA-based intron mobility pathway. Compared to the well characterized endonuclease-mediated intron homing, reverse splicing is less specific and represents a likely explanation for many intron transpositions into new genomic sites. However, the frequency and general role of an RNA-based mobility pathway in the spread of natural group I introns is still unclear. We have used the twin-ribozyme intron (Dir.S956-1) from the myxomycete Didymium iridis to test how a mobile group I intron containing a homing endonuclease gene (HEG) selects between potential insertion sites in the small subunit (SSU) rRNA in vitro, in Escherichia coli and in yeast. Surprisingly, the results show a site-specific RNA-based targeting of Dir.S956-1 into its natural (S956) SSU rRNA site. Our results suggest that reverse splicing, in addition to the established endonuclease-mediated homing mechanism, potentially accounts for group I intron spread into the homologous sites of different strains and species.

The molecular evolution and structural organization of self-splicing group I introns at position 516 in nuclear SSU rDNA of myxomycetes

Group I introns are relatively common within nuclear ribosomal DNA of eukaryotic microorganisms, especially in myxomycetes. Introns at position S516 in the small subunit ribosomal RNA gene are particularly common, but have a sporadic occurrence in myxomycetes. Fuligo septica, Badhamia gracilis, and Physarum flavicomum, all members of the family Physaraceae, contain related group IC1 introns at this site. The F. septica intron was studied at the molecular level and found to self-splice as naked RNA and to generate full-length intron RNA circles during incubation. Group I introns at position S516 appear to have a particularly widespread distribution among protists and fungi. Secondary structural analysis of more than 140 S516 group I introns available in the database revealed five different types of organization, including IC1 introns with and without His-Cys homing endonuclease genes, complex twin-ribozyme introns, IE introns, and degenerate group I-like introns. Both intron structural and phylogenetic analyses indicate a multiple origin of the S516 introns during evolution. The myxomycete introns are related to S516 introns in the more distantly related brown algae and Acanthamoeba species. Possible mechanisms of intron transfer both at the RNA- and DNA-levels are discussed in order to explain the observed widespread, but scattered, phylogenetic distribution.

Long-term evolution of the S788 fungal nuclear small subunit rRNA group I introns

More than 1000 group I introns have been identified in fungal rDNA. Little is known, however, of the splicing and secondary structure evolution of these ribozymes. Here, we use a combination of comparative and biochemical methods to address the evolution and splicing of a vertically inherited group I intron found at position 788 in the fungal small subunit (S) rRNA. The ancestral state of the S788 intron contains a highly conserved core and an extended P5 domain typical of IC1 introns. In contrast, the more derived introns have lost most of P5, and have an accelerated divergence rate within the core region with three functionally important substitutions that unambiguously separate them from the ancestral pool. Of 14 S788 group I introns that were tested for splicing, five, all of the ancestral type, were able to self-splice and produced intron RNA circles in vitro. The more derived S788 introns did not self-splice, and potentially rely on fungal-specific factors to facilitate splicing. In summary, we demonstrate one possible fate of vertically inherited group I introns, the loss of secondary structure elements, lessened selective constraints in the intron core, and ultimately, dependence on host-mediated splicing.

Hydrolytic cleavage by a group I intron ribozyme is dependent on RNA structures not important for splicing

DiGIR2 is the group I splicing-ribozyme of the mobile twin-ribozyme intron Dir.S956-1, present in Didymium nuclear ribosomal DNA. DiGIR2 is responsible for intron excision, exon ligation, 3'-splice site hydrolysis, and full-length intron RNA circle formation. We recently reported that DiGIR2 splicing (intron excision and exon ligation) competes with hydrolysis and subsequent full-length intron circularization. Here we present experimental evidence that hydrolysis at the 3'-splice site in DiGIR2 is dependent on structural elements within the P9 subdomain not involved in splicing. Whereas the GCGA tetra-loop in P9b was found to be important in hydrolytic cleavage, probably due to tertiary RNA-RNA interactions, the P9.2 hairpin structure was found to be essential for hydrolysis. The most important positions in P9.2 include three adenosines in the terminal loop (L9.2) and a consensus kink-turn motif in the proximal stem. We suggest that the L9.2 adenosines and the kink-motif represent key regulatory elements in the splicing and hydrolytic reaction pathways.

Twelve Group I introns in the same pre-rRNA transcript of the myxomycete Fuligo septica: RNA processing and evolution

The ribosomal DNA region of the myxomycete Fuligo septica was investigated and found to contain 12 group I introns (four in the small subunit and eight in the large subunit ribosomal RNAs). We have performed molecular and phylogenetic analyses to provide insight into intron structure and function, intron-host biology, and intron origin and evolution. The introns vary in size from 398 to 943 nt, all lacking detectable open reading frames. Secondary structure models revealed considerable structural diversity, but all, except one (subclass IB), represent the common group IC1 intron subclass. In vitro splicing analysis revealed that 10 of the 12 introns were able to self-splice as naked RNA, but all 12 introns were able to splice out from the precursor rRNA in vivo as evaluated by reverse transcription PCR analysis on total F. septica RNA. Furthermore, RNA processing analyses in vitro and in vivo showed that 10 of 12 introns perform hydrolytic cleavage at the 3' splice site, as well as intron circularization. Full-length intron RNA circles were detected in vivo. The order of splicing was analyzed by a reverse transcription PCR approach on cellular RNA, but no strict order of intron excision could be detected. Phylogenetic analysis indicated that most Fuligo introns were distantly related to each other and were independently gained in ribosomal DNA during evolution.

Active self-splicing group I introns in 23S rRNA genes of hyperthermophilic bacteria, derived from introns in eukaryotic organelles

Group I introns are common in the 23 rRNA genes of mitochondria and chloroplasts. Often, they encode "homing endonucleases," which target highly conserved gene sequences and drive interorganellar intron mobility, even across species and genus lines. Most bacterial 23S rRNA genes show these same endonuclease-sensitive target sequences. However, only two bacterial 23S rRNA genes are known to contain group I introns: that of Simkania negevensis [Everett, K. D., Kahane, S., Bush, R. M. & Friedman, M. G. (1999) J. Bacteriol. 181, 4734-4740], where the intron is not spliced and probably limits growth, and that of Coxiella burnetii [Seshadri, R., et al. (2003) Proc. Natl. Acad. Sci. USA 100, 5455-5460], where no direct evidence of splicing exists. Both bacteria are intracellular parasites and might have acquired introns from eukaryotic hosts. Here we provide direct evidence for splicing, and evolutionary evidence for mobility, of group I introns in the 23S rRNA genes of several free-living hyperthermophilic bacteria of the genus Thermotoga. These bacteria do not live closely with eukaryotes, but phylogenetic analyses suggest that their introns were also acquired from eukaryotic (probably algal) organelles. In vivo, their introns must be spliced at temperatures approaching 90 degrees C, making them the most thermostable natural ribozymes so far described. We demonstrate that at least some of these introns can also self-splice in vitro.

The group I-like ribozyme DiGIR1 mediates alternative processing of pre-rRNA transcripts in Didymium iridis

During starvation induced encystment, cells of the myxomycete Didymium iridis accumulate a 7.5-kb RNA that is the result of alternative processing of pre-rRNA. The 5' end corresponds to an internal processing site cleaved by the group I-like ribozyme DiGIR1, located within the twin-ribozyme intron Dir.S956-1. The RNA retains the majority of Dir.S956-1 including the homing endonuclease gene and a small spliceosomal intron, the internal transcribed spacers ITS1 and ITS2, and the large subunit rRNA lacking its two group I introns. The formation of this RNA implies cleavage by DiGIR1 in a new RNA context, and presents a new example of the cost to the host of intron load. This is because the formation of the 7.5-kb RNA is incompatible with the formation of functional ribosomal RNA from the same transcript. In the formation of the 7.5-kb RNA, DiGIR1 catalysed cleavage takes place without prior splicing performed by DiGIR2. This contrasts with the processing order leading to mature rRNA and I-DirI mRNA in growing cells, suggesting an interplay between the two ribozymes of a twin-ribozyme intron.

DiGIR1 and NaGIR1: naturally occurring group I-like ribozymes with unique core organization and evolved biological role

The group I-like ribozyme GIR1 is a unique example of a naturally occurring ribozyme with an evolved biological function. GIR1 generates the 5'-end of a nucleolar encoded messenger RNA involved in intron mobility. GIR1 is found as a cis-cleaving ribozyme within two very different rDNA group I introns (twin-ribozyme introns) in distantly related organisms. The Didymium GIR1 (DiGIR1) and Naegleria GIR1 (NaGIR1) share fundamental features in structural organization and reactivity, and display significant differences when compared to the related group I splicing ribozymes. GIR1 lacks the characteristic P1 segment present in all group I splicing ribozymes, it has a novel core organization, and it catalyses two site-specific hydrolytic cleavages rather than splicing. DiGIR1 and NaGIR1 appear to have originated from eubacterial group I introns in order to fulfil a common biological challenge: the expression of a protein encoding gene in a nucleolar context.

Characterization of the self-splicing products of two complex Naegleria LSU rDNA group I introns containing homing endonuclease genes

The two group I introns Nae.L1926 and Nmo.L2563, found at two different sites in nuclear LSU rRNA genes of Naegleria amoebo-flagellates, have been characterized in vitro. Their structural organization is related to that of the mobile Physarum intron Ppo.L1925 (PpLSU3) with ORFs extending the L1-loop of a typical group IC1 ribozyme. Nae.L1926, Nmo.L2563 and Ppo.L1925 RNAs all self-splice in vitro, generating ligated exons and full-length intron circles as well as internal processed excised intron RNAs. Formation of full-length intron circles is found to be a general feature in RNA processing of ORF-containing nuclear group I introns. Both Naegleria LSU rDNA introns contain a conserved polyadenylation signal at exactly the same position in the 3' end of the ORFs close to the internal processing sites, indicating an RNA polymerase II-like expression pathway of intron proteins in vivo. The intron proteins I-NaeI and I-NmoI encoded by Nae.L1926 and Nmo.L2563, respectively, correspond to His-Cys homing endonucleases of 148 and 175 amino acids. I-NaeI contains an additional sequence motif homologous to the unusual DNA binding motif of three antiparallel beta sheets found in the I-PpoI endonuclease, the product of the Ppo.L1925 intron ORF.

Flanking sequences with an essential role in hydrolysis of a self-cleaving group I-like ribozyme

DiGIR1 is a group I-like ribozyme derived from the mobile twin ribozyme group I intron DiSSU1 in the nuclear ribosomal DNA of the myxomycete Didymium iridis. This ribozyme is responsible for intron RNA processing in vitro and in vivo at two internal sites close to the 5'-end of the intron endo-nuclease open reading frame and is a unique example of a group I ribozyme with an evolved biological function. DiGIR1 is the smallest functional group I ribozyme known from nature and has an unusual core organization including the 6 bp P15 pseudoknot. Here we report results of functional and structural analyses that identify RNA elements critical for hydrolysis outside the DiGIR1 ribozyme core moiety. Results from deletion analysis, disruption/compensation mutagenesis and RNA structure probing analysis all support the existence of two new segments, named P2 and P2.1, involved in the hydrolysis of DiGIR1. Significant decreases in the hydrolysis rate, k (obs), were observed in disruption mutants involving both segments. These effects were restored by compensatory base pairing mutants. The possible role of P2 is to tether the ribozyme core, whereas P2.1 appears to be more directly involved in catalysis.

Functional alpha-fragment of beta-galactosidase can be expressed from the mobile group I intron PpLSU3 embedded in yeast pre-ribosomal RNA derived from the chromosomal rDNA locus

PpLSU3, a mobile group I intron found in the ribo-somal RNA genes of Physarum polycephalum, encodes the I-PpoI homing endonuclease. This enzyme represents one of the rare cases in nature where a protein is expressed from an RNA polymerase I transcript. Our previous results showed that the full length intron, but not a further processed species, is the messenger for I-PpoI, implying a role of the untranslated region (UTR) in gene expression. To study the function of the 3'-UTR in expression of the endonuclease and in splicing of the intron, we replaced the I-PpoI gene in PpLSU3 with the gene for the alpha-fragment of Escherichia coli beta-galactosidase, and then integrated this chimeric intron into all the chromosomal rDNA repeats of yeast. The resulting cells synthesized functional alpha-fragment, as evidenced by a complementation assay analogous to that used in E.coli. The beta-galactosidase activity thus provides an unusual and potentially valuable readout for Pol I transcription from chromosomal rDNA. This is the first example in which a eucaryotic homing endonuclease gene has been successfully replaced by a heterologous gene. Using deletion mutagenesis and a novel randomization approach with the alpha-fragment as a reporter, we found that a small segment of the 3'-UTR dramatically influences both splicing and protein expression.

In vivo expression of the nucleolar group I intron-encoded I-dirI homing endonuclease involves the removal of a spliceosomal intron

The Didymium iridis DiSSU1 intron is located in the nuclear SSU rDNA and has an unusual twin-ribozyme organization. One of the ribozymes (DiGIR2) catalyses intron excision and exon ligation. The other ribozyme (DiGIR1), which along with the endonuclease-encoding I-DirI open reading frame (ORF) is inserted in DiGIR2, carries out hydrolysis at internal processing sites (IPS1 and IPS2) located at its 3' end. Examination of the in vivo expression of DiSSU1 shows that after excision, DiSSU1 is matured further into the I-DirI mRNA by internal DiGIR1-catalysed cleavage upstream of the ORF 5' end, as well as truncation and polyadenylation downstream of the ORF 3' end. A spliceosomal intron, the first to be reported within a group I intron and the rDNA, is removed before the I-DirI mRNA associates with the polysomes. Taken together, our results imply that DiSSU1 uses a unique combination of intron-supplied ribozyme activity and adaptation to the general RNA polymerase II pathway of mRNA expression to allow a protein to be produced from the RNA polymerase I-transcribed rDNA.

I-PpoI, the endonuclease encoded by the group I intron PpLSU3, is expressed from an RNA polymerase I transcript

PpLSU3, a mobile group I intron in the rRNA genes of Physarum polycephalum, also can home into yeast chromosomal ribosomal DNA (rDNA) (D. E. Muscarella and V. M. Vogt, Mol. Cell. Biol. 13:1023-1033, 1993). By integrating PpLSU3 into the rDNA copies of a yeast strain temperature sensitive for RNA polymerase I, we have shown that the I-PpoI homing endonuclease encoded by PpLSU3 is expressed from an RNA polymerase I transcript. We have also developed a method to integrate mutant forms of PpLSU3 as well as the Tetrahymena intron TtLSU1 into rDNA, by expressing I-PpoI in trans. Analysis of I-PpoI expression levels in these mutants, along with subcellular fractionation of intron RNA, strongly suggests that the full-length excised intron RNA, but not RNAs that are further cleaved, serves as or gives rise to the mRNA.

Group I twintrons: genetic elements in myxomycete and schizopyrenid amoeboflagellate ribosomal DNAs

Protists are unicellular eukaryotes which represent a significant fraction of the global biodiversity. The myxomycete Didymium and the schizopyrenid amoeboflagellate Naegleria are distantly related protists. However, we have noted several striking similarities in life cycle, cell morphology, and ribosomal DNA organization between these organisms. Both have multicopy nuclear extrachromosomal ribosomal DNAs. Here the small subunit ribosomal RNA genes are interrupted by an optional group I twintron, a novel category among the group I introns. Group I twintrons are mobile self-splicing introns of 1.3-1.4 kb in size, with a complex organization at the RNA level. A group I twintron consists of two distinct ribozymes (catalytic RNAs) with different functions in RNA processing, and an open reading frame encoding a functional homing endonuclease--all with prospects of application as molecular tools in biotechnology. Updated RNA secondary structure models of group I twintrons, as well as an example of in vitro ribozyme activity, are presented. We suggest that the group I twintrons have been independently established in myxomycetes and schizopyrenid amoeboflagellates by horizontal gene transfer due to a combination of the phagocytotic behavior in natural environments and the extrachromosomal multicopy nature of ribosomal DNA.

Multiple origins of fungal group I introns located in the same position of nuclear SSU rRNA gene

The archiascomycetous fungus Protomyces pachydermus has two group I introns within the nuclear small subunit (nSSU) rRNA gene. One of these introns has an internal open reading frame (ORF) that encodes a predicted protein of 228 amino acid residues. On the other hand, Protomyces macrosporus has two group I introns that insert at the same positions as P. pachydermus, which have no ORF. Each alignment was constructed with Protomyces group I introns located in the same position and other introns retrieved by the BLAST Search. Each phylogenetic tree based on the alignment shows that Protomyces introns are monophyletic but the relationships among fungal introns do not reflect on the fungal phylogeny. Therefore, it is suggested that two different horizontal transfers of group I introns occurred at the early stage of Protomyces species diversification.

Crystallization and preliminary X-ray studies of I-PpoI: a nuclear, intron-encoded homing endonuclease from Physarum polycephalum

The homing endonuclease I-PpoI is encoded by an optional third intron, Pp LSU 3, found in nuclear, extrachromosomal copies of the Physarum polycephalum 26S rRNA gene. This endonuclease promotes the lateral transfer or "homing" of its encoding intron by recognizing and cleaving a partially symmetric, 15 bp homing site in 26S rDNA alleles that lack the Pp LSU 3 intron. The open reading frame encoding I-PpoI has been subcloned, and the endonuclease has been overproduced in E. coli. Purified recombinant I-PpoI has been co-crystallized with a 21 bp homing site DNA duplex. The crystals belong to space group P3(1)21, with unit cell dimensions a = b = 114 A, c = 89 A. The results of initial X-ray diffraction experiments indicate that the asymmetric unit contains an enzyme homodimer and one duplex DNA molecule, and that the unit cell has a specific volume of 3.4 A3/dalton. These experiments also provide strong evidence that I-PpoI contains several bound zinc ions as part of its structure.

New loop-loop tertiary interactions in self-splicing introns of subgroup IC and ID: a complete 3D model of the Tetrahymena thermophila ribozyme

BACKGROUND

Group I introns self-splice via two consecutive trans-esterification reactions in the presence of guanosine cofactor and magnesium ions. Comparative sequence analysis has established that a catalytic core of about 120 nucleotides is conserved in all known group I introns. This core is generally not sufficient for activity, however, and most self-splicing group I introns require non-conserved peripheral elements to stabilize the complete three-dimensional (3D) structure. The physico-chemical properties of group I introns make them excellent systems for unraveling the structural basis of the RNA-RNA interactions responsible for promoting the self-assembly of complex RNAs.

RESULTS

We present phylogenetic and experimental evidence for the existence of three additional tertiary base pairings between hairpin loops within peripheral components of subgroup IC1 and ID introns. Each of these new long range interactions, called P13, P14 and P16, involves a terminal loop located in domain 2. Although domains 2 of IC and ID introns share very strong sequence similarity, their terminal loops interact with domains 5 and 9 (subgroup IC1) and domain 6 (subgroup ID). Based on these tertiary contacts, comparative sequence analysis, and published experimental results such as Fe(II)-EDTA protection patterns, we propose 3D models for two entire group I introns, the subgroup IC1 intron in the large ribosomal precursor RNA of Tetrahymena thermophila and the SdCob.1 subgroup ID intron found in the cytochrome b gene of Saccharomyces douglasii.

CONCLUSIONS

Three-dimensional models of group I introns belonging to four different subgroups are now available. They all emphasize the modular and hierarchical organization of the architecture of group I introns and the widespread use of base-pairings between terminal hairpin loops for stabilizing the folded and active structures of large and complex RNA molecules.