Analysis of class I introns in a mitochondrial plasmid associated with senescence of Podospora anserina reveals extraordinary resemblance to the Tetrahymena ribosomal intron

The Cryphonectria parasitica mitochondrial rns gene: plasmid-like elements, introns and homing endonucleases

The mt-rns gene of Cryphonectria parasitica is 9872bp long and includes two group I and two group II introns. An analysis of intronic protein-encoding sequences revealed that LAGLIDADG ORFs, which usually are associated with group I introns, were transferred at least twice into group II introns. A plasmid-like mitochondrial element (plME) that appears in high amounts in previously mutagen-induced mit1 and mit2 hypovirulent mutants of the Ep155 standard virulent strain of C. parasitica was found to be derived from a short region of the mt-rns gene, including the exon 1 and most of the first intron. The plME is a 4.2-kb circular, multimeric DNA and an autonomously-replicating mtDNA fragment. Although sexual transmission experiments indicate that the plME does not directly cause hypovirulence, its emergence is one manifestation of the many complex molecular and genetic events that appear to underlie this phenotype.

The sporadic occurrence of a group I intron-like element in the mtDNA rnl gene of Ophiostoma novo-ulmi subsp. americana

The presence of group I intron-like elements within the U7 region of the mtDNA large ribosomal subunit RNA gene (rnl) was investigated in strains of Ophiostoma novo-ulmi subsp. americana from Canada, Europe and Eurasia, and in selected strains of O. ips, O. minus, O. piceae, O. ulmi, and O. himal-ulmi. This insertion is of interest as it has been linked previously to the generation of plasmid-like mtDNA elements in diseased strains of O. novo-ulmi. Among 197 O. novo-ulmi subsp. americana strains tested, 61 contained a 1.6kb insertion within the rnl-U7 region and DNA sequence analysis suggests the presence of a group I intron (IA1 type) that encodes a potential double motif LAGLIDADG homing endonuclease-like gene (HEG). Phylogenetic analysis of rnl-U7 intron encoded HEG-like elements supports the view that double motif HEGs originated from a duplication event of a single-motif HEG followed by a fusion event that combined the two copies into one open reading frame (ORF). The data also show that rnl-U7 intron encoded ORFs belong to a clade that includes ORFs inserted into different types of group I introns, e.g. IB, ID, IC3, IA1, present within a variety of different mtDNA genes, such as the small ribosomal subunit RNA gene (rns), apo-cytochrome b gene (cob), NADH dehydrogenase subunit 5 (nad5), cytochrome oxidase subunit 1 gene (coxI), and ATPase subunit 9 gene (atp9). We also compared the occurrence of the rnl-U7 intron in our collection of 227 strains with the presence of the rnl-U11 group I intron and concluded that the U7 intron appears to be an optional element and the U11 intron is probably essential among the strains tested.

Optional mitochondrial introns and evidence for a homing-endonuclease gene in the mtDNA rnl gene in Ophiostoma ulmi s. lat

Strains of Ophiostoma ulmi, O. novo-ulmi subsp. americana, O. novo-ulmi subsp. novo-ulmi and O. himal-ulmi were examined for optional introns/insertions within the following mitochondrial genes: small subunit RNA gene (rns), large ribosomal subunit gene (rnl) and the cytochrome oxidase subunit I gene (coxI). Insertions were noted in the rns and coxI genes in strains of O. ulmi, the less aggressive species, but absent in strains of the more aggressive O. novo-ulmi subsp. americana. Strains of all species examined had a group I intron present in the U11 region of the mitochondrial-rnl gene. In all but two strains of O. novo-ulmi subsp. americana, this rnl-U11 intron was about 1.5 kb in length whereas a 2.6 kb version of this element was present in all strains representing O. ulmi, O. novo-ulmi subsp. novo-ulmi, and Ophiostoma himal-ulmi. Irrespective of size, this intron based on RNA folds is a class IA1 group I intron and it encodes a putative ORF for the rps3 ribosomal protein. The size variation of the rnl-U11 intron was examined in detail for two strains of O. novo-ulmi subsp. americana and sequence data suggests the presence of a complex ORF within the 2.6 kb version of this intron; here a homing endonuclease-like gene has been inserted in frame and fused to the carboxyl-terminus of the putative rps3 coding region. The mitochondrial optional introns/insertions in combination with nuclear markers might be useful in distinguishing among the various species and subspecies of the O. ulmi s. lat. complex.

Comparative mitochondrial genomics in zygomycetes: bacteria-like RNase P RNAs, mobile elements and a close source of the group I intron invasion in angiosperms

To generate data for comparative analyses of zygomycete mitochondrial gene expression, we sequenced mtDNAs of three distantly related zygomycetes, Rhizopus oryzae, Mortierella verticillata and Smittium culisetae. They all contain the standard fungal mitochondrial gene set, plus rnpB, the gene encoding the RNA subunit of the mitochondrial RNase P (mtP-RNA) and rps3, encoding ribosomal protein S3 (the latter lacking in R.oryzae). The mtP-RNAs of R.oryzae and of additional zygomycete relatives have the most eubacteria-like RNA structures among fungi. Precise mapping of the 5′ and 3′ termini of the R.oryzae and M.verticillata mtP-RNAs confirms their expression and processing at the exact sites predicted by secondary structure modeling. The 3′ RNA processing of zygomycete mitochondrial mRNAs, SSU-rRNA and mtP-RNA occurs at the C-rich sequence motifs similar to those identified in fission yeast and basidiomycete mtDNAs. The C-rich motifs are included in the mature transcripts, and are likely generated by exonucleolytic trimming of RNA 3′ termini. Zygomycete mtDNAs feature a variety of insertion elements: (i) mtDNAs of R.oryzae and M.verticillata were subject to invasions by double hairpin elements; (ii) genes of all three species contain numerous mobile group I introns, including one that is closest to an intron that invaded angiosperm mtDNAs; and (iii) at least one additional case of a mobile element, characterized by a homing endonuclease insertion between partially duplicated genes [Paquin,B., Laforest,M.J., Forget,L., Roewer,I., Wang,Z., Longcore,J. and Lang,B.F. (1997) Curr. Genet., 31, 380–395]. The combined mtDNA-encoded proteins contain insufficient phylogenetic signal to demonstrate monophyly of zygomycetes.

Intronic GIY-YIG endonuclease gene in the mitochondrial genome of Podospora curvicolla: evidence for mobility

Endonuclease genes encoded in invasive introns are themselves supposed to be mobile elements which, during evolution, have colonized pre-existing introns converting them into invasive elements. This hypothesis is supported by numerous data concerning the LAGLI-DADG subclass of intronic endonucleases. Less is known about the GIY-YIG ORFs which constitute another family of endonucleases. In this paper we describe the presence of one optional GIY-YIG ORF in the second intron of the mitochondrial cytochrome b gene in the fungus Podospora curvicolla. We show that this GIY-YIG ORF is efficiently transferred from an ORF-containing intron to an ORF-less allele. We also show that the products of both the GIY-YIG ORF and the non-canonical LAGLI-DADG-GIY-YIG ORF, which is generated by its integration, have endonuclease activities which recognize and cut the insertion site of the optional sequence. This constitutes the first direct evidence for potential mobility of an intronic GIY-YIG endonuclease. We discuss the role that such a mobile sequence could have played during evolution.

Solution conformation of a five-nucleotide RNA bulge loop from a group I intron

We present the solution conformation, determined by NMR spectroscopy, of a five-nucleotide RNA bulge loop. The bulge interrupts the stem of a 25-nucleotide RNA hairpin, and its sequence and flanking sequences are those of a conserved bulge from a Group I intron. The secondary structure of the bulge loop in the hairpin context is that predicted by the secondary structure prediction algorithm of Zuker. It differs, however, from the secondary structure deduced from sequence covariation of the bulge in the context of the functionally folded Group I introns and observed in the crystal structure of an independently folding domain of the Group I intron from Tetrahymena thermophila. This difference represents an exception to the heierarchical model of RNA folding in which preformed elements of secondary structure interact to form a tertiary structure. The three-dimensional structure of the bulge loop is characterized by discontinuous base stacking. Adjacent adenines stack with each other and with the flanking double helices. However, the position of the central uracil is not well defined by NOE distance constraints and is a point of discontinuity in the base stacking.

Protein footprinting approach to mapping DNA binding sites of two archaeal homing enzymes: evidence for a two-domain protein structure

The archaeal intron-encoded homing enzymes I-PorI and I-DmoI belong to a family of endonucleases that contain two copies of a characteristic LAGLIDADG motif. These endonucleases cleave their intron- or intein- alleles site-specifically, and thereby facilitate homing of the introns or inteins which encode them. The protein structure and the mechanism of DNA recognition of these homing enzymes is largely unknown. Therefore, we examined these properties of I-PorI and I-DmoI by protein footprinting. Both proteins were susceptible to proteolytic cleavage within regions that are equidistant from each of the two LAGLIDADG motifs. When complexed with their DNA substrates, a characteristic subset of the exposed sites, located in regions immediately after and 40-60 amino acids after each of the LAGLIDADG motifs, were protected. Our data suggest that the enzymes are structured into two, tandemly repeated, domains, each containing both the LAGLIDADG motif and two putative DNA binding regions. The latter contains a potentially novel DNA binding motif conserved in archaeal homing enzymes. The results are consistent with a model where the LAGLIDADG endonucleases bind to their non-palindromic substrates as monomeric enzymes, with each of the two domains recognizing one half of the DNA substrate.

Sequence effects on RNA bulge-induced helix bending and a conserved five-nucleotide bulge from the group I introns

Bulge loops introduce bends in RNA double helices. Thus, a role for bulge loops in the tertiary folding of RNA is to orient helical elements. The location, size, and sequence of a five-nucleotide bulge are conserved in many of the self-splicing group I introns. We have used gel electrophoretic analysis of helix bending to test the hypothesis that this bulge loop is conserved to control the angle between the flanking helices. Interruption of an RNA duplex by the five-nucleotide bulge of the group I intron from Tetrahymena thermophila results in an electrophoretically retarded species, indicative of bending by the bulge. However, mutation of conserved bases in the bulge has a small effect on the retardation, suggesting that the average induced bend angle is not strongly dependent on the conserved sequence. Electrophoretic analysis of a mixture of bulged duplexes containing all five-nucleotide bulges reveals that most five-nucleotide bulge sequences induce bends that are similar to the bend induced by the conserved bulge. We have calibrated relative electrophoretic mobilities with bends of known magnitude, and characterized the distribution of bulge sequences among bend angles. Though the entire range of bend angles induced by different five-nucleotide bulges is from approximately 45 degrees to 75 degrees, most ( > 85%) five-nucleotide bulge loops induce bends between 65 degrees and 75 degrees. We have identified several of the anomalous five-nucleotide bulge sequences that induce bends of magnitude smaller than 65 degrees. They are generally, though not universally, pyrimidine-rich.

Primary and secondary structure analyses of the rDNA group-I introns of the Zygnematales (Charophyta)

The Zygnematales (Charophyta) contain a group-I intron (subgroupIC1) within their nuclear-encoded small subunit ribosomal DNA (SSU rDNA) coding region. This intron, which is inserted after position 1506 (relative to the SSU rDNA of Escherichia coli), is proposed to have been vertically inherited since the origin of the Zygnematales approximately 350-400 million years ago. Primary and secondary structure analyses were carried out to model group-I intron evolution in the Zygnematales. Secondary structure analyses support genetic data regarding sequence conservation within regions known to be functionally important for in vitro self-splicing of group-I introns. Comparisons of zygnematalean group-I intron secondary structures also provided some new insights into sequences that may have important roles in in vivo RNA splicing. Sequence analyses showed that sequence divergence rates and the nucleotide compositions of introns and coding regions within any one taxon varied widely, suggesting that the "1506" group-I introns and rDNA coding regions in the Zygnematales evolve independently.

Nucleotide-sequence analysis indicates that a DNA plasmid in a diseased isolate of Ophiostoma novo-ulmi is derived by recombination between two long repeat sequences in the mitochondrial large subunit ribosomal RNA gene

The nucleotide sequence of a mitochondrial plasmid (2234 bp) in a diseased isolate of Ophiostoma novo-ulmi, and sequences of the mitochondrial DNA that overlap and flank the plasmid end-points, have been determined. The plasmid was shown to be derived from the O. novo-ulmi mitochondrial large subunit ribosomal RNA gene and contained most of intron 1, the whole of exon 2, and probably the first part of intron 2. Within intron 1 there is an open reading frame with the potential to encode a 323 amino-acid polypeptide which contained dodecapeptide sequences typical of RNA maturases and DNA endonucleases. The endpoints of the plasmid in the mtDNA were located within two 90-bp direct imperfect repeat sequences, one of which comprised the last 7 bp of exon 1 and the first 83 bp of intron 1 whilst the other comprised the last 7 bp of exon 2 and the first 83 bp of intron 2. It is proposed that the Ld plasmid was generated by intramolecular recombination between these two repeats with the crossover point probably within the last 15 bp.

A group-I intron in the mitochondrial small subunit ribosomal RNA gene of Sclerotinia sclerotiorum

A 1,380-bp intervening sequence within the mitochondrial small subunit ribosomal RNA (mt SSU rRNA) gene of the fungus Sclerotinia sclerotiorum has been sequenced and identified as a group-I intron. This is the first report of an intron in the mt SSU rRNA gene. The intron shows close similarity in secondary structure to the subgroup-IC2 introns from Podospora (ND3i1, ND5i2, and COIi5) and Neurospora (ND5i1). The intron has an open reading frame (ORF) that encodes a putative protein of 420 amino acids which contains two copies of the LAGLI-DADG motif. The ORF belongs to a family of ORFs identified in Podospora (ND3i1, ND4Li1, ND4Li2, ND5i2, and COIi5) and Neurospora (ND5i1). The putative 420-aa polypeptide is also similar to a site-specific endonuclease in the chloroplast large subunit ribosomal RNA (LSU rRNA) gene of the green alga Chlamydomonas eugametos. In each clone of S. sclerotiorum examined, including several clones which were sampled over a 3-year period from geographically separated sites, all isolates either had the intron or lacked the intron within the mt SSU rRNA gene. Screening by means of Southern hybridization and PCR amplification detected the intron in the mt SSU rRNA genes of S. minor, S. trifoliorum and Sclerotium cepivorum, but not in other members of the Sclerotiniaceae, such as Botrytis anamorphs of Botryotinia spp., or in other ascomycetous and basidiomycetous fungi.

Phylogenetic evidence for horizontal transfer of an intervening sequence between species in a spirochete genus

The 23S rRNA genes (rrl genes) of some strains of certain species of the spirochete genus Leptospira carry an intervening sequence (IVS) of 485 to 759 bases flanked by terminal inverted repeat and encoding an open reading frame for a putative protein of over 120 amino acids. The structure and the sporadic distribution of the IVS suggest that it might be a mobile element that can be horizontally transferred within or between species. Phylogenetic hypotheses based on the sequences for six IVS open reading frames from various species were compared with hypotheses constructed by using DNA sequences from the 16S rRNA gene (rrs), which is not closely linked to rrl in this genus. The predicted phylogenies for the IVS and rrs differed in a major respect: one strain that claded with L. weillii in the tree based on the rrs data claded with L. noguchi in the tree based on the IVS data. Neither set of data supported a tree in which this strain was constrained to be in the same clade as was supported by the other set of data. This result indicates a probable horizontal transfer of the IVS from a recent ancestor of L. noguchi to a recent ancestor of one of the L. weillii strains. This observation is the first indication of horizontal transfer of elements encoded on the chromosomes of spirochetes.

Transposition of a group II intron

Among mobile genetic elements, self-splicing introns are of particular interest. They belong to either group I or group II depending on their three-dimensional structure. Homing, the systematic intron invasion of an intronless gene when it encounters its homologous intron-bearing allele, is the only means for intron mobility so far demonstrated. It depends on the activity of the intron-encoded protein and is very specific for the acceptor site. Intron transposition, the transfer of an intron to a novel site, predicted on the basis of phylogenetic studies and in vitro reverse-splicing experiments, has been proposed to be responsible for evolutionary intron spreading. Here we present results from polymerase chain reaction experiments consistent with transposition of a group II intron. This event is proposed to account for the site-specific deletion in the mitochondrial chromosome of the fungus Podospora anserina that is associated with the premature death syndrome and might also be involved in the senescence process affecting this species.

The single group-I intron in the chloroplast rrnL gene of Chlamydomonas humicola encodes a site-specific DNA endonuclease (I-ChuI)

The single group-I intron (ChLSU.1) in the chloroplast (cp) large subunit rRNA-encoding gene (rrnL) of the green alga Chlamydomonas humicola is located at a position at which no introns have previously been characterized in other systems. In the present study, the nucleotide (nt) sequence of this 1118-bp intron was found to contain an internal open reading frame (ORF) that potentially encodes a basic protein of 218 amino acid residues. The putative C. humicola protein features two copies of the LAGLI-DADG motif and is part of the family of intron-encoded proteins comprising the endonucleases (ENases), I-SceI, I-SceIV and I-CsmI. Expression of the ChLSU.1 intron ORF in vitro in the presence of a 260-bp DNA fragment containing the exon 1-2 junction of an intronless version of the C. humicola rrnL resulted in specific cleavage of the DNA fragment very close to the intron insertion site. This novel intron-encoded ENase, designated I-ChuI, was also shown to generate a staggered cut with 4-nt (CTCG) 3'-OH overhangs 2 bp downstream from the intron insertion site.

Group I introns in the liverwort mitochondrial genome: the gene coding for subunit 1 of cytochrome oxidase shares five intron positions with its fungal counterparts

The complete nucleotide sequence of the mitochondrial DNA (mtDNA) from a liverwort, Marchantia polymorpha, contains thirty-two introns. Twenty-five of these introns possess the characteristic secondary structures and consensus sequences of group II introns. The remaining seven are group I introns, six of which happen to interrupt the gene coding for subunit 1 of cytochrome oxidase (cox1). Interestingly, the insertion sites of one group II and four group I introns in the cox1 gene coincide with those of the respective fungal mitochondrial interns. Moreover, comparison of the four group I introns with their fungal counterparts shows that group I introns inserted at identical genomic sites in different organisms are indeed related to one another, in terms of the peptide sequences generated from the complete or fragmental ORFs encoded by these introns. At the same time, the liverwort introns turned out to be more divergent from their fungal cognates than the latter are from one another. We therefore conclude that vertical transmission from a common ancestor organism is the simplest explanation for the presence of cognate introns in liverwort and fungal mitochondrial genomes.

Imported sequences in the mitochondrial yeast genome identified by nucleotide linguistics

In addition to universally appearing mitochondrial (mt) genes, origins of replication and transcription start regions typical of all mt genome variants of the yeast Saccharomyces cerevisiae, the mt genomes of some of the strains contain variable sequences. These sequences are apparently largely dispensable. They are mainly composed of group-I and -II introns and intergenic open reading frames (ORFs). Many of the introns contain ORFs, some of which were shown by genetic and biochemical means to be involved in splicing and transposition of the mt introns. Some of the optional sequences are hypothesized to be mobile genetic elements. Nucleotide (nt) sequences of the mt genome of S. cerevisiae were examined by analyzing occurrences of oligodeoxyribonucleotide (oligo) 'words'. This linguistic technique had been found to be sensitive to both function and origin of the sequence [Pietrokovski et al., J. Biomol. Struct. Dyn. 7 (1990) 1251-1268]. A clear difference is found between the oligo vocabularies of the optional and basic yeast mt sequences. The difference is mainly located in protein coding segments of the optional sequences which contain conserved amino acid motifs, characteristic of intronic and intergenic ORFs. The use of nt linguistics to detect the sequence dissimilarity and its causes in yeast mitochondria provides fast and straightforward results, identifying the intronic and intergenic ORFs as DNA sequences of foreign, non-mt origin.

The mitochondrial tyrosyl-tRNA synthetase of Podospora anserina is a bifunctional enzyme active in protein synthesis and RNA splicing

The Neurospora crassa mitochondrial tyrosyl-tRNA synthetase (mt tyrRS), which is encoded by the nuclear gene cyt-18, functions not only in aminoacylation but also in the splicing of group I introns. Here, we isolated the cognate Podospora anserina mt tyrRS gene, designated yts1, by using the N. crassa cyt-18 gene as a hybridization probe. DNA sequencing of the P. anserina gene revealed an open reading frame (ORF) of 641 amino acids which has significant similarity to other tyrRSs. The yts1 ORF is interrupted by two introns, one near its N terminus at the same position as the single intron in the cyt-18 gene and the other downstream in a region corresponding to the nucleotide-binding fold. The P. anserina yts1+ gene transformed the N. crassa cyt-18-2 mutant at a high frequency and rescued both the splicing and protein synthesis defects. Furthermore, the YTS1 protein synthesized in Escherichia coli was capable of splicing the N. crassa mt large rRNA intron in vitro. Together, these results indicate that YTS1 is a bifunctional protein active in both splicing and protein synthesis. The P. anserina YTS1 and N. crassa CYT-18 proteins share three blocks of amino acids that are not conserved in bacterial or yeast mt tyrRSs which do not function in splicing. One of these blocks corresponds to the idiosyncratic N-terminal domain shown previously to be required for splicing activity of the CYT-18 protein. The other two are located in the putative tRNA-binding domain toward the C terminus of the protein and also appear to be required for splicing. Since the E. coli and yeast mt tyrRSs do not function in splicing, the adaptation of the Neurospora and Podospora spp. mt tyrRSs to function in splicing most likely occurred after the divergence of their common ancestor from yeast.

The mitochondrial tyrosyl-tRNA synthetase of Podospora anserina is a bifunctional enzyme active in protein synthesis and RNA splicing

The Neurospora crassa mitochondrial tyrosyl-tRNA synthetase (mt tyrRS), which is encoded by the nuclear gene cyt-18, functions not only in aminoacylation but also in the splicing of group I introns. Here, we isolated the cognate Podospora anserina mt tyrRS gene, designated yts1, by using the N. crassa cyt-18 gene as a hybridization probe. DNA sequencing of the P. anserina gene revealed an open reading frame (ORF) of 641 amino acids which has significant similarity to other tyrRSs. The yts1 ORF is interrupted by two introns, one near its N terminus at the same position as the single intron in the cyt-18 gene and the other downstream in a region corresponding to the nucleotide-binding fold. The P. anserina yts1+ gene transformed the N. crassa cyt-18-2 mutant at a high frequency and rescued both the splicing and protein synthesis defects. Furthermore, the YTS1 protein synthesized in Escherichia coli was capable of splicing the N. crassa mt large rRNA intron in vitro. Together, these results indicate that YTS1 is a bifunctional protein active in both splicing and protein synthesis. The P. anserina YTS1 and N. crassa CYT-18 proteins share three blocks of amino acids that are not conserved in bacterial or yeast mt tyrRSs which do not function in splicing. One of these blocks corresponds to the idiosyncratic N-terminal domain shown previously to be required for splicing activity of the CYT-18 protein. The other two are located in the putative tRNA-binding domain toward the C terminus of the protein and also appear to be required for splicing. Since the E. coli and yeast mt tyrRSs do not function in splicing, the adaptation of the Neurospora and Podospora spp. mt tyrRSs to function in splicing most likely occurred after the divergence of their common ancestor from yeast.

Recognition of the folding consensus in RNA secondary structures by the topological-filtering method

Functionally homologous RNA sequences can substantially diverge in their primary sequences but it can be reasonably assumed that they are related in their higher-degree structures. The problem to find such structures and simultaneously satisfy as far as possible the free-energy-minimization criterion, is considered here in two aspects. Firstly a quantitative measure of the folding consensus among secondary structures is defined, translating each structure into a linear representation and using the correlation theorem to compare them. Secondly an algorithm for the parallel search for secondary structures according to the free-energy-minimization criterion, but with a filtering action on the basis of the folding consensus measure is presented. The method is tested on groups of RNA sequences different in origin and in functions, for which proposals of homologous secondary structures based on experimental data exist. A comparison of the results with a blank consisting of a search on the basis of the free energy minimization alone is always performed. In these tests the method shows its ability in obtaining, from different sequences, secondary structures characterized by a high-folding consensus measure also when lower free energy but not homologous structures are possible. Two applications are also shown. The first demonstrates the transfer of experimental data available for one sequence, to a functionally related and therefore homologous one. The second application is the possibility of using a topological probe in the search for precise structural motifs.

Nucleotide sequence of the COX1 gene in Kluyveromyces lactis mitochondrial DNA: evidence for recent horizontal transfer of a group II intron

The cytochrome oxidase subunit 1 gene (COX1) in K. lactis K8 mtDNA spans 8,826 bp and contains five exons (termed E1-E5) totalling 1,602 bp that show 88% nucleotide base matching and 91% amino acid homology to the equivalent gene in S. cerevisiae. The four introns (termed K1 cox1.1-1.4) contain open reading frames encoding proteins of 786, 333, 319 and 395 amino acids respectively that potentially encode maturase enzymes. The first intron belongs to group II whereas the remaining three are group I type B. Introns K1 cox1.1, 1.3, and 1.4 are found at identical locations to introns Sc cox1.2, 1.5 a, and 1.5 b respectively from S. cerevisiae. Horizontal transfer of an intron between recent progenitors of K. lactis and S. cerevisiae is suggested by the observation that K1 cox1.1 and Sc cox1.2 show 96% base matching. Sequence comparisons between K1 cox1.3/Sc cox1.5 a and K1 cox1.4/Sc cox1.5 b suggest that these introns are likely to have been present in the ancestral COX1 gene of these yeasts. Intron K1 cox1.2 is not found in S. cerevisiae and appears at an unique location in K. lactis. A feature of the DNA sequences of the group I introns K1 cox1.2, 1.3, and 1.4 is the presence of 11 GC-rich clusters inserted into both coding and noncoding regions. Immediately downstream of the COX1 gene is the ATPase subunit 8 gene (A8) that shows 82.6% base matching to its counterpart in S. cerevisiae mtDNA.

The phage T4 nrdB intron: a deletion mutant of a version found in the wild

Bacteriophage T4 possesses three self-splicing group I introns. Two of the three introns are mobile elements; the third, in the gene encoding a subunit of the phage nucleotide reductase (nrdB), is not mobile. Because intron mobility offers a reasonable explanation for the paradoxical occurrence of large intervening sequences in a space-efficient eubacterial phage, it is puzzling that the nrdB intron is not mobile like its compatriots. We have discovered a larger nrdB intron in a closely related phage, and we infer from comparative sequence data that the T4 intron is a deletion mutant derived from this larger intron. This larger nrdB intron encodes an open reading frame of 269 codons, which we have cloned and overexpressed. The overexpressed protein shows a dsDNA endonuclease activity specific for the intronless nrdB gene, typical of mobile introns. Thus, we believe that all three introns of T4 are or were mobile "infectious introns" and that they have entered into and been maintained in the phage population by virtue of this efficient mobility.

Incipient mitochondrial evolution in yeasts. II. The complete sequence of the gene coding for cytochrome b in Saccharomyces douglasii reveals the presence of both new and conserved introns and discloses major differences in the fixation of mutations in evolution

We have determined the complete sequence of the mitochondrial gene coding for cytochrome b in Saccharomyces douglasii. The gene is 6310 base-pairs long and is interrupted by four introns. The first one (1311 base-pairs) belongs to the group ID of secondary structure, contains a fragment open reading frame with a characteristic GIY ... YIG motif, is absent from Saccharomyces cerevisiae and is inserted in the same site in which introns 1 and 2 are inserted in Neurospora crassa and Podospora anserina, respectively. The next three S. douglasii introns are homologous to the first three introns of S. cerevisiae, are inserted at the same positions and display various degrees of similarity ranging from an almost complete identity (intron 2 and 4) to a moderate one (intron 3). We have compared secondary structures of intron RNAs, and nucleotide and amino acid sequences of cytochrome b exons and intron open reading frames in the two Saccharomyces species. The rules that govern fixation of mutations in exon and intron open reading frames are different: the relative proportion of mutations occurring in synonymous codons is low in some introns and high in exons. The overall frequency of mutations in cytochrome b exons is much smaller than in nuclear genes of yeasts, contrary to what has been found in vertebrates, where mitochondrial mutations are more frequent. The divergence of the cytochrome b gene is modular: various parts of the gene have changed with a different mode and tempo of evolution.

Six group I introns and three internal transcribed spacers in the chloroplast large subunit ribosomal RNA gene of the green alga Chlamydomonas eugametos

The chloroplast large subunit rRNA gene of Chlamydomonas eugametos and its 5' flanking region encoding tRNA(Ile) (GAU) and tRNA(Ala) (UGC) have been sequenced. The DNA sequence data along with the results of a detailed RNA analysis disclosed two unusual features of this green algal large subunit rRNA gene: (1) the presence of six group I introns (CeLSU.1-CeLSU.6) whose insertion positions have not been described previously, and (2) the presence of three short internal transcribed spacers that are post-transcriptionally excised to yield four rRNA species of 280, 52, 810 and 1720 nucleotides, positioned in this order (5' to 3') in the primary transcript. Together, these RNA species can assume a secondary structure that is almost identical to that proposed for the 23 S rRNA of Escherichia coli. All three internal transcribed spacers map to variable regions of primary sequence and/or potential secondary structure, whereas all six introns lie within highly conserved regions. The first three introns are inserted within the sequence encoding the 810 nucleotide rRNA species and map within domain II of the large subunit rRNA structure; the remaining introns, found in the sequence encoding the 1720 nucleotide rRNA species, lie within either domain IV or V, as is the case for all other large subunit rDNA introns that have been documented to date. CeLSU.5 and CeLSU.6 each contain a long open reading frame (ORF) of more than 200 codons. While the CeLSU.6 ORF is not related to any known ORFs, the CeLSU.5 ORF belongs to a family of ORFs that have been identified in Podospora and Neurospora mitochondrial group I introns. The finding that a polymorphic marker showing unidirectional gene conversion during crosses between C. eugametos and Chlamydomonas moewusii is located within the CeLSU.5 ORF makes it likely that this intron is a mobile element and that its ORF encodes a site-specific endonuclease promoting the transfer of the intron DNA sequence.

The mitochondrial genome of fission yeast: inability of all introns to splice autocatalytically, and construction and characterization of an intronless genome

In this paper we report the inability of four group I introns in the gene encoding subunit I of cytochrome c oxidase (cox1) and the group II intron in the apocytochrome b gene (cob) to splice autocatalytically. Furthermore we present the characterization of the first cox1 intron in the mutator strain anar-14 and the construction and characterization of strains with intronless mitochondrial genomes. We provide evidence that removal of introns at the DNA level (termed DNA splicing) is dependent on an active RNA maturase. Finally we demonstrate that the absence of introns does not abolish homologous mitochondrial recombination.

Modelling of the three-dimensional architecture of group I catalytic introns based on comparative sequence analysis

Alignment of the 87 available sequences of group I self-splicing introns reveals numerous instances of covariation between distant sites. Some of these covariations cannot be ascribed to historical coincidences or the known secondary structure of group I introns, and are, therefore, best explained as reflecting tertiary contacts. With the help of stereochemical modelling, we have taken advantage of these novel interactions to derive a three-dimensional model of the conserved core of group I introns. Two noteworthy features of that model are its extreme compactness and the fact that all of the most evolutionarily conserved residues happen to converge around the two helices that constitute the substrate of the core ribozyme and the site that binds the guanosine cofactor necessary for self-splicing. Specific functional implications are discussed, both with regard to the way the substrate helices are recognized by the core and possible rearrangements of the introns during the self-splicing process. Concerning potential long-range interactions, emphasis is put on the possible recognition of two consecutive purines in the minor groove of a helix by a GAAA or related terminal loop.

Mutational analysis of conserved nucleotides in a self-splicing group I intron

We have constructed all single base substitutions in almost all of the highly conserved residues of the Tetrahymena self-splicing intron. Mutation of highly conserved residues almost invariably leads to loss of enzymatic activity. In many cases, activity could be regained by making additional mutations that restored predicted base-pairings; these second site suppressors in general confirm the secondary structure derived from phylogenetic data. At several positions, our suppression data can be most readily explained by assuming non-Watson-Crick base-pairings. In addition to the requirements imposed by the secondary structure, the sequence of the intron is constrained by "negative interactions", the exclusion of particular nucleotide sequences that would form undesirable secondary structures. A comparison of genetic and phylogenetic data suggests sites that may be involved in tertiary structural interactions.

The complete DNA sequence of the mitochondrial genome of Podospora anserina

The complete 94,192 bp sequence of the mitochondrial genome from race s of Podospora anserina is presented (1 kb = 10(3) base pairs). Three regions unique to race A are also presented bringing the size of this genome to 100,314 bp. Race s contains 31 group I introns (33 in race A) and 2 group II introns (3 in race A). Analysis shows that the group I introns can be categorized according to families both with regard to secondary structure and their open reading frames. All identified genes are transcribed from the same strand. Except for the lack of ATPase 9, the Podospora genome contains the same genes as its fungal counterparts, N. crassa and A. nidulans. About 20% of the genome has not yet been identified. DNA sequence studies of several excision-amplification plasmids demonstrate a common feature to be the presence of short repeated sequences at both termini with a prevalence of GGCGCAAGCTC.

Mitochondrial DNA sequence analysis of the cytochrome oxidase subunit II gene from Podospora anserina. A group IA intron with a putative alternative splice site

A 5 kb region of the 95 kb mitochondrial genome of Podospora anserina race s has been mapped and sequenced (1 kb = 10(3) base-pairs). This DNA region is continuous with the sequence for the ND4L and ND5 gene complex in the accompanying paper. We show that this sequence contains the gene for cytochrome oxidase subunit II (COII). This gene is 4 kb in length and is interrupted by a subgroup IB intron (1267 base-pairs (bp) in length) and a subgroup IA intron (1992 bp in length). This group IA intron has a long open reading frame (ORF; 472 amino acid residues) discontinuous with the upstream exon sequence. A putative alternative splice site is present, which brings the ORF into phase with the 5' exon sequence. The 5'- and 3'-flanking regions of the COII gene contain G + C-rich palindromic sequences that resemble similar sequences flanking many Neurospora crassa mitochondrial genes.

DNA sequence analysis of the mitochondrial ND4L-ND5 gene complex from Podospora anserina. Duplication of the ND4L gene within its intron

A 15 kb region of the 100 kb mitochondrial genome of Podospora anserina has been mapped and sequenced (1 kb = 10(3) base-pairs). The genes for ND4L and ND5 are identified as contiguous genes with overlapping termination and initiation codons. In race A (101 kb) the gene for ND4L (4.3 kb) has a gene duplication within an intron including a second subgroup IC intron. Race s (95 kb) lacks this second gene complex. Each intron has the identical 5' exon boundary. Secondary structure analysis showed that the closest relative of the second intron is the first intron itself. The open reading frames of the two introns are also closely related to each other as well as to their counterpart in the ND4L gene of Neurospora crassa. The 9.9 kb ND5 gene starts immediately at the termination codon of ND4L and is split by two group IB introns, one group IC intron and one group II intron. The group II intron is closely related to other group II introns although its open reading frame sequence similarity with retroviral reverse transcriptase appears to be more divergent. The similarities in secondary structure and open reading frames for these six introns are discussed.

DNA sequence analysis of the 24.5 kilobase pair cytochrome oxidase subunit I mitochondrial gene from Podospora anserina: a gene with sixteen introns

The DNA sequence of a 26.7 Kilobase pair (10(3) base pairs = 1 Kb) region of the mitochondrial genomes of races s and A from Podospora anserina was determined. Within this region, the 24.5 Kb cytochrome oxidase subunit I gene was located and its exon sequences determined by computer analysis comparisons with other fungal genes. The Podospora COI gene was interrupted by two group II introns (one in race s) and fourteen group I introns ranging in size from about 2.2 Kb to 404 bp. Earlier studies on secondary structure analysis, as well as comparison of their open reading frames (ORFs), showed that the two group II introns were closely related. The fourteen group I introns were representatives of three subgroupings (IB, C and a new category, subgroup ID). Two of these group I introns were separated by just a single exon codon. The analysis of all these introns is discussed in comparison with other fungal introns as well as with the known Podospora anserina introns.

DNA sequence analysis of the apocytochrome b gene of Podospora anserina: a new family of intronic open reading frame

The 5,969 bp (base pair) DNA sequence of the apocytochrome b mitochondrial (mt) gene of race A Podospora anserina was located in a 8.5 Kbp region. This gene contained a 2,499 bp subgroup IB and a 1,306 bp subgroup ID intron as well as a 990 bp subgroup IB intron which is present in race A but not race s. The large subgroup IB intron and the race A specific IB intron both contained potential alternate splice sites which brought their open reading frames into phase with their upstream exon sequences. All three introns were compared with regard to their secondary structures and open reading frames to the other 30 group I introns in Podospora anserina, as well as to other fungal introns. We detected a new family of intronic ORFs comprising seven P. anserina introns, several N. crassa introns, as well as the T4td bacteriophage intron. Sequence similarities to intron-encoded endonucleases were noteworthy. The DNA sequences reported here and in the accompanying paper complete the analysis of race s and race A mitochondrial DNA.

Group I introns as mobile genetic elements: facts and mechanistic speculations--a review

Group I introns form a structural and functional group of introns with widespread but irregular distribution among very diverse organisms and genetic systems. Evidence is now accumulating that several group I introns are mobile genetic elements with properties similar to those originally described for the omega system of Saccharomyces cerevisiae: mobile group I introns encode sequence-specific double-strand (ds) endoDNases, which recognize and cleave intronless genes to insert a copy of the intron by a ds-break repair mechanism. This mechanism results in: the efficient propagation of group I introns into their cognate sites; their maintenance at the site against spontaneous loss; and, perhaps, their transposition to different sites. The spontaneous loss of group I introns occurs with low frequency by an RNA-mediated mechanism. This mechanism eliminates introns defective for mobility and/or for RNA splicing. Mechanisms of intron acquisition and intron loss must create an equilibrium, which explains the irregular distribution of group I introns in various genetic systems. Furthermore, the observed distribution also predicts that horizontal transfer of intron sequences must occur between unrelated species, using vectors yet to be discovered.

Comparative and functional anatomy of group II catalytic introns--a review

The 70 published sequences of group II introns from fungal and plant mitochondria and plant chloroplasts are analyzed for conservation of primary sequence, secondary structure and three-dimensional base pairings. Emphasis is put on structural elements with known or suspected functional significance with respect to self-splicing: the exon-binding and intron-binding sites, the bulging A residue involved in lariat formation, structural domain V and two isolated base pairs, one of them involving the last intron nucleotide and the other one, the first nt of the 3' exon. Separate sections are devoted to the 29 group II-like introns from Euglena chloroplasts and to the possible relationship of catalytic group II introns to nuclear premessenger introns. Alignments of all available sequences of group II introns are provided in the APPENDIX.

Improved predictions of secondary structures for RNA

The accuracy of computer predictions of RNA secondary structure from sequence data and free energy parameters has been increased to roughly 70%. Performance is judged by comparison with structures known from phylogenetic analysis. The algorithm also generates suboptimal structures. On average, the best structure within 10% of the lowest free energy contains roughly 90% of phylogenetically known helixes. The algorithm does not include tertiary interactions or pseudoknots and employs a crude model for single-stranded regions. The only favorable interactions are base pairing and stacking of terminal unpaired nucleotides at the ends of helixes. The excellent performance is consistent with these interactions being the primary interactions determining RNA secondary structure.

Selection of the 3'-splice site in group I introns

A model for selection of 3'-splice sites in splicing of RNA precursors containing group I introns is presented. The key feature of this model is a newly identified tertiary interaction between the catalytic core of the intron and the 3'-splice site. This tertiary pairing would bring the 3'-splice site into the core of the intron, which is known to contain RNA sequences and structures essential for catalyzing the splicing reactions. The proposed tertiary interaction can coexist with P10, a pairing between 3'-exon sequences and the 'internal guide sequence' near the 5'-end of the intron. The model predicts that three RNA-RNA interactions are important in selection of 3'-splice sites: (i) binding of intron sequences with the core; (ii) pairing of exon sequences with the internal guide sequence; and (iii) binding of the terminal guanosine to an unknown site within the core.

Two group I introns with long internal open reading frames in the chloroplast psbA gene of Chlamydomonas moewusii

We report the nucleotide sequence of the chloroplast psbA gene encoding the 32 kilodalton protein of photosystem II from Chlamydomonas moewusii. Like its land plant homologues, this green algal protein consists of 353 amino acids. The C. moewusii psbA gene is composed of three exons containing 252, 11 and 90 codons and of two group I introns containing 2363 and 1807 nucleotides. Each of the introns features an internal open reading frame (ORF) that potentially encodes a basic protein of more than 300 residues. The primary sequences of the putative intron-encoded proteins are unrelated and none of them shares conserved elements with any of the proteins predicted from the group I intron sequences published so far. The first C. moewusii intron is inserted at the same position as the fourth intron of the psbA gene from Chlamydomonas reinhardtii; the second intron lies at a novel site downstream of this position. On the basis of their RNA secondary structures, the C. moewusii introns 1 and 2 can be assigned to subgroups IA and IB, respectively. However, intron 1 is not typical of subgroup IA introns, its most unusual feature being the location of the ORF in the "loop L5" region. To our knowledge, this is the first time that an ORF is located in this region of the group I intron structure.

Reverse self-splicing of the tetrahymena group I intron: implication for the directionality of splicing and for intron transposition

Using short oligoribonucleotides as ligated exon substrates, we show that splicing of the Tetrahymena rRNA group I intron is fully reversible in vitro. Incubation of ligated exon RNA with linear intron produces a molecule in which the splice site sequences of the precursor are reformed. Reversal of self-splicing is favored by high RNA concentration, high magnesium and temperature, and the absence of guanosine. 5' exon sequences that can pair with the internal guide sequence of the intron are required, whereas 3' exon sequences are not essential. Integration of the intron into ligated exon substrates that have the ability to form stem-loop structures is reduced at least one order of magnitude over short, unstructured substrates. We propose that the formation of these structures helps drive splicing in the forward direction. We also show that the Tetrahymena intron can integrate into a beta-globin transcript. This has implications for transposition of group I introns.

An optional group I intron between the chloroplast small subunit rRNA genes of Chlamydomonas moewusii and C. eugametos

We report the presence of a 402 bp group I intron in the chloroplast small subunit (SSU) rRNA gene of Chlamydomonas moewusii. The intron in inserted within the highly conserved '530 loop', at a site corresponding to positions 531-532 of the E. coli 16rRNA. Residues surrounding the insertion site almost certainly play an important role in ribosomal proofreading function as they proved to be protected by tRNAs in E. coli 16S rRNA (Moazed and Noller 1986; Stern et al. 1986). The C. moewusii intron revealed a secondary structure model which differs substantially from those of the typical subgroup IA and IB introns. This model, however, shows striking similarities with the structures of the C. reinhardtii chloroplast 23S rRNA gene intron (Rochaix et al. 1985), the S. cerevisiae mitochondrial COB3 intron (Holl et al. 1985) and the three introns of phage T4 in the nrdB, td and sunY genes (Shub et al. 1988). The SSU rRNA gene intron is absent from C. eugametos, an alga that is interfertile with C. moewusii. The presence/absence of the intron account for a 390 bp restriction fragment length polymorphism between the two algal SSU rRNA genes, a polymorphic locus that is strictly co-inherited with a tightly linked streptomycin resistance mutation (sr-2) in interspecific hybrids between the two algae.

DNA sequence and secondary structures of the large subunit rRNA coding regions and its two class I introns of mitochondrial DNA from Podospora anserina

DNA sequence analysis has shown that the gene coding for the mitochondrial (mt) large subunit ribosomal RNA (rRNA) from Podospora anserina is interrupted by two class I introns. The coding region for the large subunit rRNA itself is 3715 bp and the two introns are 1544 (r1) and 2404 (r2) bp in length. Secondary structure models for the large subunit rRNA were constructed and compared with the equivalent structure from Escherichia coli 23S rRNA. The two structures were remarkably similar despite an 800-base difference in length. The additional bases in the P. anserina rRNA appear to be mostly in unstructured regions in the 3' part of the RNA. Secondary structure models for the two introns show striking similarities with each other as well as with the intron models from the equivalent introns in Saccharomyces cerevisiae, Neurospora crassa, and Aspergillus nidulans. The long open reading frames in each intron are different from each other, however, and the nucleotide sequence similarity diverges as it proceeds away from the core structure. Each intron is located within regions of the large subunit rRNA gene that are highly conserved in both sequence and structure. Computer analysis showed that the open reading frame for intron r1 contained a common maturase-like polypeptide. The open reading frames of intron r2 appeared to be chimeric, displaying high sequence similarity with the open reading frames in the r1 and ATPase 6 introns of N. crassa.

Molecular genetics of group I introns: RNA structures and protein factors required for splicing--a review

In vivo and in vitro genetic techniques have been widely used to investigate the structure-function relationships and requirements for splicing of group-I introns. Analyses of group-I introns from extremely diverse genetic systems, including fungal mitochondria, protozoan nuclei, and bacteriophages, have yielded results which are complementary and highly consistent. In vivo genetic studies of fungal mitochondrial systems have served to identify cis-acting sequences within mitochondrial introns, and trans-acting protein products of mitochondrial and nuclear genes which are important for splicing, and to show that some mitochondrial introns are mobile genetic elements. In vitro genetic studies of the self-splicing intron within the Tetrahymena thermophila nuclear large ribosomal RNA precursor (Tetrahymena LSU intron) have been used to examine essential and nonessential RNA sequences and structures in RNA-catalyzed splicing. In vivo and in vitro genetic analysis of the intron within the bacteriophage T4 td gene has permitted the detailed examination of mutant phenotypes by analyzing splicing in vivo and self-splicing in vitro. The genetic studies combined with phylogenetic analysis of intron structure based on comparative nucleotide sequence data [Cech 73 (1988) 259-271] and with biochemical data obtained from in vitro splicing experiments have resulted in significant advances in understanding the biology and chemistry of group-I introns.