Short-term sequence evolution and vertical inheritance of the Naegleria twin-ribozyme group I intron
https://doi.org/10.1186/1471-2148-6-39
© Wikmark et al; licensee BioMed Central Ltd. 2006
Received: 28 December 2005
Accepted: 02 May 2006
Published: 02 May 2006
Abstract
Background
Ribosomal DNA of several species of the free-living Naegleria amoeba harbors an optional group I intron within the nuclear small subunit ribosomal RNA gene. The intron (Nae.S516) has a complex organization of two ribozyme domains (NaGIR1 and NaGIR2) and a homing endonuclease gene (NaHEG). NaGIR2 is responsible for intron excision, exon ligation, and full-length intron RNA circularization, reactions typical for nuclear group I intron ribozymes. NaGIR1, however, is essential for NaHEG expression by generating the 5' end of the homing endonuclease messenger RNA. Interestingly, this unusual class of ribozyme adds a lariat-cap at the mRNA.
Results
To elucidate the evolutionary history of the Nae.S516 twin-ribozyme introns we have analyzed 13 natural variants present in distinct Naegleria isolates. Structural variabilities were noted within both the ribozyme domains and provide strong comparative support to the intron secondary structure. One of the introns, present in N. martinezi NG872, contains hallmarks of a degenerated NaHEG. Phylogenetic analyses performed on separate data sets representing NaGIR1, NaGIR2, NaHEG, and ITS1-5.8S-ITS2 ribosomal DNA are consistent with an overall vertical inheritance pattern of the intron within the Naegleria genus.
Conclusion
The Nae.S516 twin-ribozyme intron was gained early in the Naegleria evolution with subsequent vertical inheritance. The intron was lost in the majority of isolates (70%), leaving a widespread but scattered distribution pattern. Why the apparent asexual Naegleria amoebae harbors active intron homing endonucleases, dependent on sexual reproduction for its function, remains a puzzle.
Keywords
Background
Naegleria is a common genus of soil and freshwater free-living amoeba of the vahlkampfiid family [1]. Naegleria apparently lack a sexual reproduction cycle since meiosis never has been observed or proven experimentally. Subsequently, a number of genetically-defined variants have been isolated in nature and proposed as distinct species [1, 2]. A typical Naegleria amoeba cell contains a distinct nucleus with a large and predominant nucleolus containing as much as 3000–5000 copies of an approximately 14-kb ribosomal DNA (rDNA) plasmid [3, 4]. Each rDNA molecule carries a single transcription unit for the ribosomal RNA (rRNA) genes. Some Naegleria isolates have been reported to contain group I intron insertions at conserved sequence sites, both within the small subunit (SSU) and large subunit (LSU) rRNA genes [5]. Introns have been noted at position 516 in SSU rDNA (i.e. a position that is homologous to corresponding position in the E. coli rRNA gene) and at positions 1921, 1926, 1949, and 2563 in LSU rDNA [5–11].
Group I introns are autocatalytic genetic elements carrying a ribozyme domain responsible for the intron self-splicing reaction, and occasionally a homing endonuclease gene (HEG) encoding an endonuclease protein directly involved in intron mobility at the DNA level [12, 13]. A group I splicing ribozyme possess a well-defined three-dimensional structure organized into three functional domains (catalytic domain, folding domain, and substrate domain) by approximately ten paired RNA segments named P1–P10 [14, 15]. The most common and best characterized of the Naegleria rDNA introns is Nae.S516. Group I introns at position 516 in SSU rDNA are relatively common among eukaryotic microorganisms with more than 250 cases reported so far [16, 17], and with both lateral and vertical inheritance patterns compared to that of host rDNA. A widespread distribution and structural diversity among the S516 group I introns have been noted, including several complex introns carrying HEGs [16, 18].
A typical Naegleria S516 intron has a twin-ribozyme organization and represents the most complex class of all group I introns known [19]. Nae.S516 consists of a small group I-like mRNA capping ribozyme (NaGIR1) and a HEG domain, both inserted into the P6b segment of a regular group IC1 splicing ribozyme (NaGIR2). Expression and functional aspects of the Naegleria S516 intron have been reported. The splicing ribozyme (NaGIR2) is responsible for the autocatalytic activity that generates intron excision and exon ligation, as well as full-length intron RNA circle formations [7]. The ability to form full-length intron RNA circles is a general property of nuclear group I introns and could be important in RNA mobility at the RNA level, or even as an intermediate in the expression of the homing endonuclease [11, 20, 21]. The Naegleria HEG (NaHEG) encodes a 245 amino acid protein that belongs to the His-Cys box family of homing endonucleases [22, 23]. The intron endonuclease recognizes and binds to a 19-bp DNA sequence flanking the S516 rDNA site and cleaves the DNA generating a five-nucleotide 3' staggered end [24, 25]. In general, group I intron endonucleases promote intron homing at the DNA-level by generating a double-stranded break in the intron-less target DNA, followed by invasion of the donor intron-containing allele and DNA repair using the intron-containing allele as template [21]. However, sexual mating is the biological framework for nuclear group I intron homing [26, 27] and it is unclear why the apparently asexual Naegleria contains and maintains homing introns.
The expression of the NaHEG is dependent on a functional NaGIR1 ribozyme, which defines the 5' end of the homing endonuclease mRNA by internal processing and modification of the excised Nae.S516 intron [7, 28, 29]. Primer extension analyses of both cellular RNA from Naegleria and in vitro transcribed intron RNA [7, 28] are consistent with the formation of a tiny lariat cap structure between nucleotide 1 and 3 of the messenger, as recently reported in the related DiGIR1 ribozyme [30]. Thus, the NaGIR1 capping ribozyme represents a novel class of ribozymes possessing a new catalytic function, which is reflected in its unique RNA architecture [29, 31].
The complex and unique structural organization of the Naegleria twin-ribozyme intron makes it interesting to investigate the evolutionary origin of the different intron domains, as well as the inheritance pattern within the Naegleria genus. Here we report several new intron variants and have performed sequence and phylogenetic analyses providing new insight into fundamental questions such as intron structure, intron-host biology, and the origin and evolution of intron HEG and ribozyme domains. The Naegleria twin-ribozyme intron serves as an attractive model system in the characterization of evolutionary processes behind a recently gained, but vertically inherited, selfish genetic element.
Results and discussion
Widespread but sporadic distribution of Nae.S516 introns within the Naegleria genus
Distribution of the Nae.S516 introns among Naegleria isolates. NJ phylogenetic tree of the ITS-rDNA based on 287 nucleotide positions and 70 Naegleria isolates. NJ bootstrap values above 50% are shown at the branches. The six clusters of Naegleria isolates (Cluster 1–6) are indicated according to [1,2]. Naegleria isolates known to contain (21 isolates) or lack (49 isolates) the Nae.S516 intron are shown in red or black, respectively. All introns are approximately 1,3 kp in size (corresponding to a twin-ribozyme organization), except N. byersi NG597 which is 375 bp (only NaGIR2). The N. fowleri branch in Cluster 1 represents 21 distinct isolates.
Structural features and sequence variability of intron domains
Consensus structure diagram of the Nae.S516 twin-ribozyme intron. The structure diagram is based on the 13 completely sequenced Nae.S516 introns from distinct natural isolates (see Table 1) and is folded according to previously reported models [7,19] with some modifications. Invariant nucleotide residues are presented as uppercase letters. Red filled circles represents a variable position in one or more intron sequences, and regions with size and structural variations are boxed. Nae.S516 contains the two distinct ribozymes NaGIR1 and NaGIR2, and the homing endonuclease gene NaHEG. IPS, internal processing site; BP, branch point nucleotide (U).
NaGIR1 is a group I-like ribozyme with an evolved biological role in intron NaHEG expression [29, 30], likely to generate a lariat cap-structure at the 5' end of the Naegleria homing endonuclease messenger [30]. The three-dimensional architecture of NaGIR1 is related to that of bacterial tRNA group I intron ribozymes [18, 29, 33], but with a unique catalytic core organization that contains the novel pseudoknot segment P15 [7, 29, 31, 34]. As seen from Figure 2, most of the core nucleotides are highly conserved among the various Naegleria capping ribozymes. Variable regions are almost exclusively located in the terminal loops of P6 and P8, the internal loop junction J5/4, and sequences flanking NaGIR1 and NaHEG.
Alignment of Naegleria homing endonuclease sequences. Identical residues compared to the N. jamiesoni T56E endonuclease sequence are shown by dots and deletions by dashes. Functional important residues involved in zinc binding and catalysis are indicated. Divergent regions within the N. martinezi NG872 sequence due to reading-frame shifts (*) are boxed. For structural comparisons to the I-Ppo I homing endonuclease, see [23].
The splicing ribozyme of Nae.S516 intron is vertically inherited in Naegleria
Phylogeny of NaGIR2 and ITS-rDNA. The NaGIR2 and ITS-rDNA NJ trees were generated from 356 and 415 nt datasets. Bootstrap values from 2000 replicates and Bayesian posterior probability values, all over 50%, are shown at branches. The values are from, in order, NJ, MP, ML, and BAY analyses. See Materials and Methods for more detailed analytical parameters. ITS-rDNA and NaGIR2 phylogenies are congruent, except for N. italica AB-T-F3 (marked in red). See Table 1 for information about sequence accession numbers.
The capping ribozyme NaGIR1 and its downstream NaHEG are evolutionary linked
Phylogeny of NaGIR1 and NaHEG. The NaGIR1 and NaHEG NJ trees were generated from 230 and 747 nt datasets. Bootstrap values from 2000 replicates and Bayesian posterior probability values, all over 50%, are shown at branches. The values are from NJ, MP, ML, and BAY analyses, respectively. See Materials and Methods for more detailed analytical parameters. NaGIR1 and NaHEG phylogenies are congruent, except for N. carteri NG055 (marked in red). See Table 1 for information about sequence accession numbers.
A Naegleria S516 group I intron with only NaHEG or only NaGIR1 insertions has never been observed, suggesting a strong linkage between the domains. Both structural and functional data give further support to a close linkage between the NaGIR1 and NaHEG domains. Jabri and Cech [34] showed that the RNA structure essential for NaGIR1 catalysis includes nucleotide residues within the NaHEG coding region. Functional experiments in yeast conclude that expression of NaHEG, and subsequent endonuclease activity in yeast extracts, is completely dependent on a functional NaGIR1 ribozyme [28]. Thus, the NaGIR1 and NaHEG domains have to be considered as one functional unit within the Nae.S516 intron.
Gain of a L5b GNRA tetraloop in NaGIR2 during Naegleria evolution
One of the best-studied tetraloop receptor interaction is the L5b-P6a tertiary structure in the Tetrahymena group I intron ribozyme [40]. Here, the GAAA loop in L5b specifically binds to the 11-nt receptor motif CCUAAG-UAUGG within the helical stem of P6a by docking into the shallow groove. The L5b-P6a interaction in Tetrahymena is essential for an efficient folding of the P4–P6 domain, and subsequently the folding and activity of the splicing ribozyme.
Structural variations in L5b and P6a of NaGIR2. A) NJ phylogenetic tree of the 13 complete Nae.S516 introns based on 1370 positions. NJ and MP bootstrap values are indicated at the branches. The presence of L5b tetra- (red), penta- (blue), or hexaloops (green) in L5b in NaGIR2 are indicated. All tetra-loops belong to the GNRA-loop family. B) Structure diagrams of P6a (grey boxes) and L9b representing the 13 twin-ribozyme Nae.S516 introns.
The primary function of GNRA tetraloops is to participate in long-range RNA-RNA interactions by specific binding to a receptor motif. A variety of receptor motifs, ranging from 4 to 12 nt, have been recognized experimentally [41–43]. Figure 6B presents secondary structure diagrams of the various NaGIR2 P6a regions and their corresponding L5b loops. A prominent difference in the P6a structure is noted among introns possessing L5b GNRA tetraloops compared to those with penta- or hexaloops. Introns with GNRA loops contain a less tightly base-paired P6a stem with several proposed exposed residues (see Figure 6B) compared to the P5b penta- or hexaloop containing introns (compare N. clarki RU30 and N. carteri NG055). We speculate that the exposed residues in P6a could be involved in RNA-RNA interactions as GNRA receptors. However, these sequences do not fit any known consensus motifs, suggesting that new motifs are yet to be experimentally identified.
Conclusion
Evolutionary aspects on the structural organization of Naegleria twin-ribozyme group I introns have been reported previously [19]. Whereas the NaGIR2 splicing ribozyme appears related to other eukaryote rDNA group IC1 intron [16], the NaGIR1 capping ribozyme has recently evolved from a bacterial tRNA group I introns [29]. Here we present phylogenetic evidence of a vertical inheritance pattern of the Nae.S516 intron in Naegleria that includes each of the domains NaGIR1, NaGIR2, and NaHEG, and corroborates a previous study based on 5 intron sequences [6]. Based on the reported distribution pattern and phylogeny, we propose the following vertical inheritance scenario for the Nae.S516 evolution. 1) A pre-organized twin-ribozyme group I intron was gained in the rDNA early in evolution of the Naegleria genus, but after the Cluster 6 branching (see Figure 1). 2) Once established, the Naegleria intron co-evolved along with its host rDNA by maintaining intron activities including intron splicing, endonuclease mRNA capping, and homing endonuclease DNA cleavage. 3) The intron was subsequently lost (see Figure 1) by sporadic deletions in most isolates (70 %). 4) Degradation of the NaHEG is initiated due to loss of biological function, and subsequent selection pressure (e.g. N. martinezi NG872), resulting in complete deletion of the NaHEG as well as its regulatory NaGIR1 capping ribozyme (e.g. N. byersi NG597) [8]. 5) The remaining introns have to improve and adjust their functions by continuous sequence evolution in order to be maintained within rDNA. Here, a recent gain of a GNRA tetraloop receptor in the P4–P6 domain would facilitate folding of the splicing ribozyme (Figure 6).
What is the biological role of a functional NaHEG in Naegleria S516 introns? There are only two reported examples addressing the biological role of nuclear group I intron HEGs in experimental settings. In sexual matings between intron-containing and intron-lacking strains of either the myxomycetes Physarum polycephalum or Didymium iridis, group I introns were shown to be mobile due to the double-strand-break-repair pathway induced by intron-encoded homing endonucleases [26, 27]. In both cases the homing endonucleases were found to cleave the group I intron lacing alleles in a highly sequence specific manner, resulting in unidirectional transfers of introns into the intron-lacking strains. This process is dependent on a sexual reproduction of the host organism, which is apparently absent in Naegleria. However, the Naegleria intron endonucleases possess hallmarks linked to a function in intron homing. First, sequence comparisons show that the Naegleria enzymes belong to the same His-Cys homing endonuclease family as the known homing endonucleases I-Ppo I and I-Dir I encoded by the mobile Physarum and Didymium group I introns [10, 11, 25]. Second, the Naegleria endonucleases cleave only intron lacking alleles flanking the intron insertion site at the SSU rRNA gene [24, 25]. Finally, artificial expression of the Naegleria endonuclease and its intron in yeast generate intron homing intermediates consistent with a homing endonuclease function [28]. Interestingly, Naegleria may occasionally perform sexual reproduction in nature since Pernin and co-workers [44, 45] reported evidence for genetic exchange in N. lovaniensis, including chromosomal recombination. Both haploid and diploid strains of the N. gruberi NEG isolate have been described based on both amoeba DNA content and UV-sensitivity [46, 47]. Perhaps the observed recombination-like feature of NaGIR1 in N. carteri NG055 (see Figure 5) is a result of rare sexual mating. This possibility remains to be experimentally explored.
Methods
Naegleria strains, DNA amplification, plasmid cloning, and DNA sequencing
The following Naegleria isolates were DNA sequenced at ITS-rDNA and Nae.S516 intron regions in this study: N. clarki RU30 (ITS-rDNA and Nae.S516); N. clarki RU42 (ITS-rDNA and Nae.S516); N. pringsheimi 1D (ITS-rDNA and Nae.S516); N. philippinensis RJTM (ITS-rDNA and Nae.S516); N. carteri NG055 (ITS-rDNA and Nae.S516); Naegleria sp. NG647 (ITS-rDNA and Nae.S516); Naegleria sp. NG358 (ITS-rDNA and Nae.S516); Naegleria sp. NG393 (ITS-rDNA and Nae.S516); Naegleria sp. NG498 (ITS-rDNA and Nae.S516); Naegleria sp. NG169 (ITS-rDNA);Naegleria sp. NG336 (ITS-rDNA); Naegleria sp. NG491 (ITS-rDNA); Naegleria sp. NG492 (ITS-rDNA). The strains without designation are under revision and will be proposed proper species names based on phylogenetic analyses of ITS1-5.8S-ITS2 sequences (De Jonckheere et al. in preparation). A complete list of all 70 Naegleria isolates included in this study is presented in Table 1. DNA samples of Naegleria strains were prepared as described previously [2], dissolved in water, and applied as template in 50 μl standard PCR reactions. Amplified product of interest were plasmid cloned into the pGEM®-T Easy Vector System I (Promega). Individual clones where DNA purified and sequenced with the ABI PRISM BigDyeTerminator Cycle Sequencing Ready Reaction Kit (Perkin-Elmer), running on an ABI Prism 377 system (Perkin-Elmer), or using the sequencing service from MWG Biotech [48]. Two or more individual clones were sequenced from all introns and ITS-rDNA regions analysed. The following oligoprimers were used in Nae.S516 PCR amplifications and DNA sequencing analyses: OP 25 (5'-CTC GAA TTC GCT CTT GGA GCT GGA ATT A-3'), OP 26 (5'-ACG AAG CTT ATT TCT AAG CCT-3'), OP 28 (5'-CAG AGG AGT TTC TTA CCT ATC-3'), OP 131 (5'-AAA CGA ATT CTA TTG ATT AGT AGT-3'), OP 946 (5'-GAA TTG AAA AAG CTT GAT-3'), OP 1200 (5'-AAA CAA ATG CTA TTG ATC A-3'), OP 1201 (5'-GAA CGT CTA GAG ACT ACA CGG-3'), OP 1042 (5'-CGA TTT TCC ATG ATT TGG G-3'), OP 1043 (5'-ATA CCT CAA CAG AGG TCC-3'), OP 1044 (5'-GGA CGT CTA GAG ACT ACA CGG-3'), OP 1045 (5'-TGA TGC ACG TAC GAA TCG GAG C-3'), OP 276 (5'-GGT AAA CAA ATC CCT GTT-3'), OP 823 (5'-TAA CCA TTT TGT ATG GGA-3'). Heteroplasmic rDNA alleles (intron-containing/intron lacing) were not observed. The following oligoprimers were used in ITS-rDNA PCR amplifications and DNA sequencing analyses: OP 918 (5'-AAC CTG CGT AGG GAT CAT TT-3'), OP 919 (5'-TTT CCT CCC CTT ATT AAT AT-3').
Sequence alignment and phylogenetic analysis
Multiple alignment of sequences were performed by using Megalign (Version 5.06) included in the Lasergene package from DNASTAR, Inc [49], Bioedit (Version 7.0.4.1; [50]), manuel refinements. Phylogenetic analyses and non-synonymous to synonymous substitution rates [51] were conducted using MEGA version 2.1 [52], PAUP* (Version 4.0 Beta) [53], and MrBayes (version 3.1) [54, 55]. Trees were built with the methods of neighbor-joining (NJ) using different distance matrixes, maximum parsimony (MP) with the branch and bound search method, as well as Bayesian analyses (BAY) and maximum likelihood (ML) using different evolutionary models. The reliabilities of the tree topologies were evaluated by bootstrapping (NJ, MP, and ML), and posterior probability (BAY).
ITS and intron sequence analyses
Two different data sets of the internal transcribed rDNA spacer region (ITS-rDNA: ITS1-5.8S-ITS2) were used. In analysis with all the 70 Naegleria isolates, only 287 nucleotide positions could unambiguously be aligned due to high sequence variation in ITS2. However, when the analysis was restricted to the 14 intron-containing Naegleria isolates we extended the ITS-rDNA region to 415 nucleotide positions. Based on the multiple sequence alignment a NJ tree was constructed with the Kimura-2 evolutionary model of substitution (K2), with pair wise deletion of gaps and bootstrapped with 2000 replications with a cut-off value of 50%. Similarly, intron trees are constructed with NJ-K2 parameters. The robustness of the tree topologies were tested by the NJ-K2 parameter (first value), MP branch and bound search criteria (second value), and ML with the HKY+G model of substitution selected by Modeltest 3.7 [56], all from 1000 replicates. The last values where constructed by running 1000000 generations of Metropolis-coupled Markov chain Monte Carlo, and trees were sampled every 100 generations (average standard deviation of split frequencies below 0.01). A consensus tree was generated from the 75% last trees to find posterior probabilities (Burn-in value = 2500).
Declarations
Acknowledgements
We thank Peik Haugen and Dag H. Coucheron for discussions. This work was supported by grants to SDJ from the Norwegian Research Council, The Norwegian Cancer Society, and The Aakre Foundation for Cancer Research.
Authors’ Affiliations
References
- De Jonckheere JF: A century of research on the amoebaflagellate genus Naegleria. Acta Protozool. 2002, 41: 309-342.Google Scholar
- De Jonckheere JF: Molecular definition and the ubiquity of species in the genus Naegleria. Protist. 2004, 155: 89-103. 10.1078/1434461000167.View ArticlePubMedGoogle Scholar
- Clark CG, Cross GA: rRNA genes of Naegleria gruberi are carried exclusively on a 14-kilobase-pair plasmid. Mol Cell Biol. 1987, 7: 3027-3031.PubMed CentralView ArticlePubMedGoogle Scholar
- Clark CG, Cross GA: Small-subunit ribosomal RNA sequence from Naegleria gruberi supports the polyphyletic origin of amoebas. Mol Biol Evol. 1988, 5: 512-518.PubMedGoogle Scholar
- Johansen S, Coucheron DH, Haugen P: Group I introns in nuclear ribosomal DNA of Naegleria. Ninth International meeting on the biology and pathogenicity of free-living amoebae proceedings. Edited by: Billot-Bonef S, Cabanes PA, Marciano-Cabral F, Pernin P, Pringuez E. 2001, Paris: John Libby Eurotext, 254-260.Google Scholar
- De Jonckheere JF: Evidence for the ancestral origin of group I introns in the SSUrDNA of Naegleria spp. J Eukaryot Microbiol. 1994, 41: 457-463.View ArticlePubMedGoogle Scholar
- Einvik C, Decatur WA, Embley TM, Vogt VM, Johansen S: Naegleria nucleolar introns contain two group I ribozymes with different functions in RNA splicing and processing. RNA. 1997, 3: 710-720.PubMed CentralPubMedGoogle Scholar
- De Jonckheere JF, Brown S: Loss of the ORF in the SSUrDNA group I intron of one Naegleria lineage. Nucleic Acids Res. 1994, 22: 3925-3927.PubMed CentralView ArticlePubMedGoogle Scholar
- De Jonckheere JF, Brown S: Three different group I introns in the nuclear large subunit ribosomal DNA of the amoeboflagellate Naegleria. Nucleic Acids Res. 1998, 26: 456-461. 10.1093/nar/26.2.456.PubMed CentralView ArticlePubMedGoogle Scholar
- De Jonckheere JF, Brown S: A novel ORF-containing group I intron with His-Cys box in the LSU rDNA of Naegleria. Acta Protozool. 2001, 40: 27-31.Google Scholar
- Haugen P, De Jonckheere DF, Johansen S: Characterization of the self-splicing products of two complex Naegleria LSU rDNA group I introns containing homing endonuclease genes. Eur J Biochem. 2002, 269: 641-1649. 10.1046/j.1432-1327.2002.02802.x.View ArticleGoogle Scholar
- Cech TR: Self-splicing of group I introns. Annu Rev Biochem. 1990, 59: 543-568. 10.1146/annurev.bi.59.070190.002551.View ArticlePubMedGoogle Scholar
- Haugen P, Simon DM, Bhattacharya D: The natural history of group I introns. Trends Genet. 2005, 21: 111-119. 10.1016/j.tig.2004.12.007.View ArticlePubMedGoogle Scholar
- Lehnert V, Jaeger L, Michel F, Westhof E: New loop-loop tertiary interactions in self-splicing introns of subgroup IC and ID: a complete 3D model of the Tetrahymena thermophila ribozyme. Chem Biol. 1996, 3: 993-1009. 10.1016/S1074-5521(96)90166-0.View ArticlePubMedGoogle Scholar
- Guo F, Gooding AR, Cech TR: Structure of the Tetrahymena ribozyme: base triple sandwich and metal ion at the active site. Mol Cell. 2004, 16: 351-362.PubMedGoogle Scholar
- Haugen P, Coucheron DH, Rønning SB, Haugli K, Johansen S: The molecular evolution and structural organization of self-splicing group I introns at position 516 in nuclear SSU rDNA of myxomycetes. J Eukaryot Microbiol. 2003, 50: 283-292. 10.1111/j.1550-7408.2003.tb00135.x.View ArticlePubMedGoogle Scholar
- Comparative RNA Web site. [http://www.rna.icmb.utexas.edu]
- Bhattacharya D, Cannone JJ, Gutell RR: Group I intron lateral transfer between red and brown algal ribosomal RNA. Curr Genet. 2001, 40: 82-90. 10.1007/s002940100227.View ArticleGoogle Scholar
- Einvik C, Elde M, Johansen S: Group I twintrons: genetic elements in myxomycete and schizopyrenid amoeboflagellate ribosomal DNAs. J Biotechnol. 1998, 64: 63-74. 10.1016/S0168-1656(98)00104-7.View ArticlePubMedGoogle Scholar
- Nielsen H, Fiskaa T, Birgisdottir AB, Haugen P, Einvik C, Johansen S: The ability to form full-length intron RNA circles is a general property of nuclear group I introns. RNA. 2003, 9: 1464-1475. 10.1261/rna.5290903.PubMed CentralView ArticlePubMedGoogle Scholar
- Birgisdottir ÅB, Johansen S: Site-specific reverse splicing of a HEG-containing group I intron in ribosomal DNA. Nucleic Acids Res. 2005, 33: 2042-2051. 10.1093/nar/gki341.PubMed CentralView ArticlePubMedGoogle Scholar
- Johansen S, Embley TM, Willassen NP: A family of nuclear homing endonucleases. Nucleic Acids Res. 1993, 21: 4405-PubMed CentralView ArticlePubMedGoogle Scholar
- Chevalier BS, Stoddard BL: Homing endonucleases: structural and functional insight into the catalysts of intron/intein mobility. Nucleic Acids Res. 2001, 29: 3757-3774. 10.1093/nar/29.18.3757.PubMed CentralView ArticlePubMedGoogle Scholar
- Elde M, Haugen P, Willassen NP, Johansen S: I- Nja I, a nuclear intron-encoded homing endonuclease from Naegleria, generates a pentanucleotide 3' cleavage-overhang within a 19 base-pair partially symmetric DNA recognition site. Eur J Biochem. 1999, 259: 281-288. 10.1046/j.1432-1327.1999.00035.x.View ArticlePubMedGoogle Scholar
- Elde M, Willassen NP, Johansen S: Functional characterization of isoschizomeric His-Cys box homing endonucleases from Naegleria. Eur J Biochem. 2000, 267: 7257-7266. 10.1046/j.1432-1327.2000.01862.x.View ArticlePubMedGoogle Scholar
- Muscarella DE, Vogt VM: A mobile group I intron in the nuclear rDNA of Physarum polycephalum. Cell. 1989, 56: 443-454. 10.1016/0092-8674(89)90247-X.View ArticlePubMedGoogle Scholar
- Johansen S, Elde M, Vader A, Haugen P, Haugli K, Haugli F: In vivo mobility of a group I twintron in nuclear ribosomal DNA of the myxomycete Didymium iridis. Mol Microbiol. 1997, 24: 737-745. 10.1046/j.1365-2958.1997.3921743.x.View ArticlePubMedGoogle Scholar
- Decatur WA, Johansen S, Vogt VM: Expression of the Naegleria intron endonuclease is dependent on a functional group I self-cleaving ribozyme. RNA. 2000, 6: 616-627. 10.1017/S1355838200992203.PubMed CentralView ArticlePubMedGoogle Scholar
- Johansen S, Einvik C, Nielsen H: DiGIR1 and NaGIR1. Naturally occurring group I-like ribozymes with unique core organization and evolved biological role. Biochimie. 2002, 84: 905-912. 10.1016/S0300-9084(02)01443-8.View ArticlePubMedGoogle Scholar
- Nielsen H, Westhof E, Johansen S: An mRNA is capped by a 2', 5' lariat catalyzed by a group I-like ribozyme. Science. 2005, 309: 1584-1586. 10.1126/science.1113645.View ArticlePubMedGoogle Scholar
- Einvik C, Nielsen H, Westhof E, Michel F, Johansen S: Group I-like ribozymes with a novel core organization perform obligate sequential hydrolytic cleavages at two processing sites. RNA. 1998, 4: 530-541. 10.1017/S1355838298971758.PubMed CentralView ArticlePubMedGoogle Scholar
- Johansen S, Haugen P: A new nomenclature of group I introns in ribosomal DNA. RNA. 2001, 7: 935-936. 10.1017/S1355838201010500.PubMed CentralView ArticlePubMedGoogle Scholar
- Adams PL, Stahley MR, Kosek AB, Wang J, Strobel SA: Crystal structure of a self-splicing group I intron with both exons. Nature. 2004, 430: 45-50. 10.1038/nature02642.View ArticlePubMedGoogle Scholar
- Jabri E, Cech TR: In vitro selection of the Naegleria GIR1 ribozyme identifies three base changes that dramatically improve activity. RNA. 1998, 4: 1481-1492. 10.1017/S1355838298981237.PubMed CentralView ArticlePubMedGoogle Scholar
- Weiss MA, Narayana N: RNA recognition by arginine-rich peptide motifs. Biopolymers. 1998, 48: 167-180. 10.1002/(SICI)1097-0282(1998)48:2<167::AID-BIP6>3.0.CO;2-8.View ArticlePubMedGoogle Scholar
- Chen Y, Varani G: Protein families and RNA recognition. FEBS J. 2005, 272: 2088-2097. 10.1111/j.1742-4658.2005.04650.x.View ArticlePubMedGoogle Scholar
- Haugen P, Huss VAR, Nielsen H, Johansen S: Complex group I introns in nuclear SSU rDNA of red and green algae: evidence of homing-endonuclease pseudogenes in the Bangiophyceae. Curr Genet. 1999, 36: 345-353. 10.1007/s002940050509.View ArticlePubMedGoogle Scholar
- Haugen P, Reeb V, Lutzoni F, Bhattacharya D: The evolution of homing endonuclease genes and group I introns in nuclear rDNA. Mol Biol Evol. 2004, 21: 129-140. 10.1093/molbev/msh005.View ArticlePubMedGoogle Scholar
- Haugen P, Wikmark OG, Vader A, Coucheron DH, Sjøttem E, Johansen SD: The recent transfer of a homing endonuclease gene. Nucleic Acids Res. 2005, 33: 2734-2741. 10.1093/nar/gki564.PubMed CentralView ArticlePubMedGoogle Scholar
- Cate JH, Gooding AR, Podell E, Zhou K, Golden BL, Kundrot CE, Cech TR, Doudna JA: Crystal structure of a group I ribozyme domain: principles of RNA packing. Science. 1996, 273: 1678-1685.View ArticlePubMedGoogle Scholar
- Jaeger L, Michel F, Westhof E: Involvement of a GNRA tetraloop in long-range RNA tertiary interactions. J Mol Biol. 1994, 236: 1271-1276. 10.1016/0022-2836(94)90055-8.View ArticlePubMedGoogle Scholar
- Costa M, Michel F: Rules for RNA recognition of GNRA tetraloops deduced by in vitro selection: comparison with in vivo evolution. EMBO J. 1997, 16: 3289-3302. 10.1093/emboj/16.11.3289.PubMed CentralView ArticlePubMedGoogle Scholar
- Ikawa Y, Naito D, Aono N, Shiraishi H, Inoue T: A conserved motif in group IC3 introns is a new class of GNRA receptor. Nucleic Acids Res. 1999, 27: 1859-1865. 10.1093/nar/27.8.1859.PubMed CentralView ArticlePubMedGoogle Scholar
- Cariou ML, Pernin P: First evidence for diploidy and genetic recombination in free-living amoebae of the genus Naegleria on the basis of electrophoretic variation. Genetics. 1987, 115: 265-270.PubMed CentralPubMedGoogle Scholar
- Pernin P, Ataya A, Cariou ML: Genetic structure of natural populations of the free-living amoeba, Naegleria lovaniensis. Evidence for sexual reproduction. Heredity. 1992, 68: 173-181.View ArticleGoogle Scholar
- Fulton C: Amebo-flagellates as research partners: the laboratory biology of Naegleria and Tetramitus. Methods in Cell Physiology. Edited by: Prescott DM. 1970, New York: Academic Press, 4: 341-473.Google Scholar
- Fulton C: Naegleria : a research partner for cell and development biology. J Euk Microbiol. 1993, 40: 520-532.View ArticleGoogle Scholar
- MWG Home. [http://www.mwg-biotech.com]
- DNASTAR – Software for molecular biology. [http://www.dnastar.com]
- Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999, 41: 95-98.Google Scholar
- Nei M, Kumar S: Molecular Evolution and Phylogenetics. 2000, Oxform University Press, New YorkGoogle Scholar
- Kumar S, Tamura K, Jakobsen IB, Nei M: MEGA2: Molecular Evolutionary Genetics Analysis software. 2001, Arizona State University, Tempe, Arizona, USAGoogle Scholar
- Swofford DL: PAUP*: Phylogenetic analysis using parsimony (* and other methods). Version 4.0b8. 2002, Sinauer, Sunderland, MAGoogle Scholar
- Huelsenbeck JP, Ronquist F: MrBayes: Bayesian inference of phylogeny. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.View ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP: MRBAYES 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMedGoogle Scholar
- Posada D, Crandall KA: Modeltest: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.View ArticlePubMedGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.