Distribution, sequence homology, and homing of group I introns among T-even-like bacteriophages: evidence for recent transfer of old introns

Self-splicing group I introns are being found in an increasing number of bacteriophages. Most introns contain an open reading frame coding for a homing endo-nuclease that confers mobility to both the intron and the homing endonuclease gene (HEG). The frequent occurrence of intron/HEG has raised questions whether group I introns are spread via horizontal transfer between phage populations. We have determined complete sequences for the known group I introns among T-even-like bacteriophages together with sequences of the intron-containing genes td, nrdB, and nrdD from phages with and without introns. A previously uncharacterized phage isolate, U5, is shown to contain all three introns, the only phage besides T4 found with a "full set" of these introns. Sequence analysis of td and nrdB genes from intron-containing and intronless phages provides evidence that recent horizontal transmission of introns has occurred among the phages. The fact that several of the HEGs have suffered deletions rendering them non-functional implies that the homing endonucleases are of no selective advantage to the phage and are rapidly degenerating and probably dependent upon frequent horizontal transmissions for maintenance within the phage populations. Several of the introns can home to closely related intronless phages during mixed infections. However, the efficiency of homing varies and is dependent on homology in regions flanking the intron insertion site. The occurrence of optional genes flanking the respective intron-containing gene can strongly affect the efficiency of homing. These findings give further insight into the mechanisms of propagation and evolution of group I introns among the T-even-like bacteriophages.     

A latent intron-encoded maturase is also an endonuclease needed for intron mobility

Some yeast mitochondrial introns encode proteins that promote either splicing (maturases) or intron propagation via gene conversion (the fit1 endonuclease). We surveyed introns in the coxl gene for their ability to engage in gene conversion and found that the group I intron, al4 alpha, was efficiently transmitted to genes lacking it. An endonucleolytic cleavage is detectable in recipient DNA molecules near the site of intron insertion in vivo and in vitro. Conversion is dependent on an intact al4 alpha open reading frame. This intron product is a latent maturase, but these data show that it is also a potent endonuclease involved in recombination. Dual function proteins that cleave DNA and facilitate RNA splicing may have played a pivotal role in the propagation and tolerance of introns.     

The inconsistent distribution of introns in the T-even phages indicates recent genetic exchanges.

Group I self-splicing introns are present in the td, nrdB and sunY genes of bacteriophage T4. We previously reported that whereas the td intron is present in T2, T4 and T6, the nrdB intron is present in T4 only. These studies, which argue in favor of introns as mobile genetic elements, have been extended by defining the distribution of all three T4 introns in a more comprehensive collection of T2, T4 and T6 isolates. The three major findings are as follows: First, all three introns are inconsistently distributed throughout the T-even phage family. Second, different T2 isolates have different intron complements, with T2H and T2L having no detectable introns. Third, the intron open reading frames are inherited or lost as a unit with their respective flanking intron core elements. Furthermore, exon sequences flanking sites where introns are inserted in the T4 td, sunY and nrdB genes were determined for all the different T-even isolates studied. Six of eighteen residues surrounding the junction sequences are identical. In contrast, a comprehensive comparison of exon sequences in intron plus and intron minus variants of the sunY gene indicate that sequence changes are concentrated around the site of intron occurrence. This apparent paradox may be resolved by hypothesizing that the recombination events responsible for intron acquisition or loss require a consensus sequence, while these same events result in sequence heterogeneity around the site.     

Intron mobility in phage T4 is dependent upon a distinctive class of endonucleases and independent of DNA sequences encoding the intron core: mechanistic and evolutionary implications

Although mobility of the phylogenetically widespread group I introns appears to be mechanistically similar, the phage T4 intron-encoded endonucleases that promote mobility of the td and sunY introns are different from their eukaryotic counterparts. Most notably, they cleave at a distance from the intron insertion sites. The td enzyme was shown to cleave 23-26 nt 5' and the sunY endonuclease 13-15 nt 3' to the intron insertion site to generate 3-nt or 2-nt 3'-OH extensions, respectively. The absolute coconversion of exon markers between the distant cleavage and insertion sites is consistent with the double-strand-break repair model for intron mobility. As a further critical test of the model we have demonstrated that the mobility event is independent of DNA sequences that encode the catalytic intron core structure. Thus, in derivatives in which the lacZ or kanR coding sequences replace the intron, these marker genes are efficiently inserted into intron-minus alleles when the cognate endonuclease is provided in trans. The process is therefore endonuclease-dependent, rather than dependent on the intron per se. These findings, which imply that the endonucleases rather than the introns themselves were the primordial mobile elements, are incorporated into a model for the evolution of mobile introns.     

Purification and substrate specificity of a T4 phage intron-encoded endonuclease.

The T4 phage td intron-encoded endonuclease (I-Tev I) cleaves the intron-deleted td gene (td delta I) 23 nucleotides upstream of the intron insertion site on the noncoding strand and 25 nucleotides upstream of this site on the coding strand, to generate a 2-base hydroxyl overhang in the 3' end of each DNA strand. I-Tev I-157, a truncated form in which slightly more than one third (88 residues) of the endonuclease is deleted, was purified to homogeneity and shown to possess endonuclease activity similar to that of I-TEV I, the full-length enzyme (245 residues). The minimal length of the td delta I gene that was cleaved by I-Tev I and I-Tev I-157 has been determined to be exactly 39 basepairs, from -27 (upstream in exon1) to +12 (downstream in exon2) relative to the intron insertion site. Similar to the full-length endonuclease, I-Tev I-157 cuts the intronless thymidylate synthase genes from such diverse organisms as Escherichia coli, Lactobacillus casei and the human. The position and nature of the in vitro endonucleolytic cut in these genes are homologous to those in td delta I. Point mutational analysis of the td delta I substrate based on the deduced consensus nucleotide sequence has revealed a very low degree of specificity on either side of the cleavage site, for both the full-length and truncated I-TEV I.     

A site-specific endonuclease and co-conversion of flanking exons associated with the mobile td intron of phage T4

The product of the td intron open reading frame (ORF) of phage T4 is required for high-frequency transfer of the intervening sequence from intron-plus (In+) to intron-minus (In-) alleles. In vivo studies have demonstrated that the td ORF product targets cleavage of td In- DNA, and that cleavage is correlated with intron inheritance [Quirk et al., Cell 56 (1989) 455-465]. In the present study we show by in vitro synthesis of the td intron ORF product, that the protein possesses endonuclease activity and efficiently cleaves double-stranded DNA at or near the site of intron integration. In addition, we demonstrate that intron insertion is accompanied by co-conversion of the flanking exon sequences. Co-conversion of markers within 50 nt surrounding the site of intron insertion occurred at a high frequency (80-100%), and decreased at greater distance from the intervening sequence. Co-conversion may provide a mechanism for maintaining exon-intron RNA contacts required for accurate splicing of the relocated intron. Cleavage of target DNA by an intron endonuclease and co-conversion of flanking exon sequences are both features associated with mobile introns of eukaryotes, indicating a common mechanism for intron transfer in the eukaryotic and prokaryotic kingdoms.     

Multiple self-splicing introns in bacteriophage T4: evidence from autocatalytic GTP labeling of RNA in vitro

RNA from T4-infected cells yielded multiple end-labeled species when incubated with alpha-32P-GTP under self-splicing conditions. One of these corresponds to the previously identified intron from the td gene of T4, while others appear to represent additional group I introns in T4. Two loci distinct from the td gene were found to hybridize to a mixed alpha-32P-GTP-labeled T4 RNA probe. These mapped in or near the unlinked genes nrdB and nrdC. A fragment from the nrdB region that contains the intron has been cloned and shown to generate characteristic group I splice products with RNA synthesized in vivo and in vitro. Multiple introns, and the prospect that these occur within several genes in the same metabolic pathway, suggest a possible regulatory role for splicing in T4.     

The evolution of homing endonuclease genes and group I introns in nuclear rDNA

Group I introns are autonomous genetic elements that can catalyze their own excision from pre-RNA. Understanding how group I introns move in nuclear ribosomal (r)DNA remains an important question in evolutionary biology. Two models are invoked to explain group I intron movement. The first is termed homing and results from the action of an intron-encoded homing endonuclease that recognizes and cleaves an intronless allele at or near the intron insertion site. Alternatively, introns can be inserted into RNA through reverse splicing. Here, we present the sequences of two large group I introns from fungal nuclear rDNA, which both encode putative full-length homing endonuclease genes (HEGs). Five remnant HEGs in different fungal species are also reported. This brings the total number of known nuclear HEGs from 15 to 22. We determined the phylogeny of all known nuclear HEGs and their associated introns. We found evidence for intron-independent HEG invasion into both homologous and heterologous introns in often distantly related lineages, as well as the "switching" of HEGs between different intron peripheral loops and between sense and antisense strands of intron DNA. These results suggest that nuclear HEGs are frequently mobilized. HEG invasion appears, however, to be limited to existing introns in the same or neighboring sites. To study the intron-HEG relationship in more detail, the S943 group I intron in fungal small-subunit rDNA was used as a model system. The S943 HEG is shown to be widely distributed as functional, inactivated, or remnant ORFs in S943 introns.     

The Agaricus bisporus cox1 gene: the longest mitochondrial gene and the largest reservoir of mitochondrial group i introns

In eukaryotes, introns are located in nuclear and organelle genes from several kingdoms. Large introns (up to 5 kbp) are frequent in mitochondrial genomes of plant and fungi but scarce in Metazoa, even if these organisms are grouped with fungi among the Opisthokonts. Mitochondrial introns are classified in two groups (I and II) according to their RNA secondary structure involved in the intron self-splicing mechanism. Most of these mitochondrial group I introns carry a "Homing Endonuclease Gene" (heg) encoding a DNA endonuclease acting in transfer and site-specific integration ("homing") and allowing intron spreading and gain after lateral transfer even between species from different kingdoms. Opposed to this gain mechanism, is another which implies that introns, which would have been abundant in the ancestral genes, would mainly evolve by loss. The importance of both mechanisms (loss and gain) is matter of debate. Here we report the sequence of the cox1 gene of the button mushroom Agaricus bisporus, the most widely cultivated mushroom in the world. This gene is both the longest mitochondrial gene (29,902 nt) and the largest group I intron reservoir reported to date with 18 group I and 1 group II. An exhaustive analysis of the group I introns available in cox1 genes shows that they are mobile genetic elements whose numerous events of loss and gain by lateral transfer combine to explain their wide and patchy distribution extending over several kingdoms. An overview of intron distribution, together with the high frequency of eroded heg, suggests that they are evolving towards loss. In this landscape of eroded and lost intron sequences, the A. bisporus cox1 gene exhibits a peculiar dynamics of intron keeping and catching, leading to the largest collection of mitochondrial group I introns reported to date in a Eukaryote.     

Characterization of I-Ppo, an intron-encoded endonuclease that mediates homing of a group I intron in the ribosomal DNA of Physarum polycephalum.

A novel and only recently recognized class of enzymes is composed of the site-specific endonucleases encoded by some group I introns. We have characterized several aspects of I-Ppo, the endonuclease that mediates the mobility of intron 3 in the ribosomal DNA of Physarum polycephalum. This intron is unique among mobile group I introns in that it is located in nuclear DNA. We found that I-Ppo is encoded by an open reading frame in the 5' half of intron 3, upstream of the sequences required for self-splicing of group I introns. Either of two AUG initiation codons could start this reading frame, one near the beginning of the intron and the other in the upstream exon, leading to predicted polypeptides of 138 and 160 amino acid residues. The longer polypeptide was the major form translated in vitro in a reticulocyte extract. From nuclease assays of proteins synthesized in vitro with partially deleted DNAs, we conclude that both polypeptides possess endonuclease activity. We also have expressed I-Ppo in Escherichia coli, using a bacteriophage T7 RNA polymerase expression system. The longer polypeptide also was the predominant form made in this system. It showed enzymatic activity in bacteria in vivo, as demonstrated by the cleavage of a plasmid carrying the target site. Like several other intron-encoded endonucleases, I-Ppo makes a four-base staggered cut in its ribosomal DNA target sequence, very near the site where intron 3 becomes integrated in crosses of intron 3-containing and intron 3-lacking Physarum strains.     

Homing endonuclease I-TevIII: dimerization as a means to a double-strand break

Homing endonucleases are unusual enzymes, capable of recognizing lengthy DNA sequences and cleaving site-specifically within genomes. Many homing endonucleases are encoded within group I introns, and such enzymes promote the mobility reactions of these introns. Phage T4 has three group I introns, within the td, nrdB and nrdD genes. The td and nrdD introns are mobile, whereas the nrdB intron is not. Phage RB3 is a close relative of T4 and has a lengthier nrdB intron. Here, we describe I-TevIII, the H–N–H endonuclease encoded by the RB3 nrdB intron. In contrast to previous reports, we demonstrate that this intron is mobile, and that this mobility is dependent on I-TevIII, which generates 2-nt 3′ extensions. The enzyme has a distinct catalytic domain, which contains the H–N–H motif, and DNA-binding domain, which contains two zinc fingers required for interaction with the DNA substrate. Most importantly, I-TevIII, unlike the H–N–H endonucleases described so far, makes a double-strand break on the DNA homing site by acting as a dimer. Through deletion analysis, the dimerization interface was mapped to the DNA-binding domain. The unusual propensity of I-TevIII to dimerize to achieve cleavage of both DNA strands underscores the versatility of the H–N–H enzyme family.     

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.     

The I-CeuI endonuclease recognizes a sequence of 19 base pairs and preferentially cleaves the coding strand of the Chlamydomonas moewusii chloroplast large subunit rRNA gene.

The I-CeuI endonuclease is a member of the growing family of homing endonucleases that catalyse mobility of group I introns by making a double-strand break at the homing site of these introns in cognate intronless alleles during genetic crosses. In a previous study, we have shown that a short DNA fragment of 26 bp, encompassing the homing site of the fifth intron in the Chlamydomonas eugametos chloroplast large subunit rRNA gene (Ce LSU.5), was sufficient for I-CeuI recognition and cleavage. Here, we report the recognition sequence of the I-CeuI endonuclease, as determined by random mutagenesis of nucleotide positions adjacent to the I-CeuI cleavage site. Single-base substitutions that completely abolish endonuclease activity delimit a 15-bp sequence whereas those that reduce the cleavage rate define a 19-bp sequence that extends from position -7 to position +12 with respect to the Ce LSU.5 intron insertion site. As the other homing endonucleases that have been studied so far, the I-CeuI endonuclease recognizes a non-symmetric degenerate sequence. The top strand of the recognition sequence is preferred for I-CeuI cleavage and the bottom strand most likely determines the rate of double-strand breaks.     

Fungal origin by horizontal transfer of a plant mitochondrial group I intron in the chimeric coxI gene of Peperomia

We present phylogenetic evidence that a group I intron in an angiosperm mitochondrial gene arose recently by horizontal transfer from a fungal donor species. A 1,716-bp fragment of the mitochondrial coxI gene from the angiosperm Peperomia polybotrya was amplified via the polymerase chain reaction and sequenced. Comparison to other coxI genes revealed a 966-bp group I intron, which, based on homology with the related yeast coxI intron aI4, potentially encodes a 279-amino-acid site-specific DNA endonuclease. This intron, which is believed to function as a ribozyme during its own splicing, is not present in any of 19 coxI genes examined from other diverse vascular plant species. Phylogenetic analysis of intron origin was carried out using three different tree-generating algorithms, and on a variety of nucleotide and amino acid data sets from the intron and its flanking exon sequences. These analyses show that the Peperomia coxI gene intron and exon sequences are of fundamentally different evolutionary origin. The Peperomia intron is more closely related to several fungal mitochondrial introns, two of which are located at identical positions in coxI, than to identically located coxI introns from the land plant Marchantia and the green alga Prototheca. Conversely, the exon sequence of this gene is, as expected, most closely related to other angiosperm coxI genes. These results, together with evidence suggestive of co-conversion of exonic markers immediately flanking the intron insertion site, lead us to conclude that the Peperomia coxI intron probably arose by horizontal transfer from a fungal donor, using the double-strand-break repair pathway. The donor species may have been one of the symbiotic mycorrhizal fungi that live in close obligate association with most plants. Correspondence to: J.C. Vaughn     

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.     

Unconventional GIY-YIG homing endonuclease encoded in group I introns in closely related strains of the Bacillus cereus group

Several group I introns have been previously found in strains of the Bacillus cereus group at three different insertion sites in the nrdE gene of the essential nrdIEF operon coding for ribonucleotide reductase. Here, we identify an uncharacterized group IA intron in the nrdF gene in 12 strains of the B. cereus group and show that the pre-mRNA is efficiently spliced. The Bacillus thuringiensis ssp. pakistani nrdF intron encodes a homing endonuclease, denoted I-BthII, with an unconventional GIY-(X)₈-YIG motif that cleaves an intronless nrdF gene 7 nt upstream of the intron insertion site, producing 2-nt 3′ extensions. We also found four additional occurrences of two of the previously reported group I introns in the nrdE gene of 25 sequenced B. thuringiensis and one B. cereus strains, and one non-annotated group I intron at a fourth nrdE insertion site in the B. thuringiensis ssp. Al Hakam sequenced genome. Two strains contain introns in both the nrdE and the nrdF genes. Phylogenetic studies of the nrdIEF operon from 39 strains of the B. cereus group suggest several events of horizontal gene transfer for two of the introns found in this operon.     

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.     

Naegleria nucleolar introns contain two group I ribozymes with different functions in RNA splicing and processing.

We have characterized the structural organization and catalytic properties of the large nucleolar group I introns (NaSSU1) of the different Naegleria species N. jamiesoni, N. andersoni, N. italica, and N. gruberi. NaSSU1 consists of three distinct RNA domains: an open reading frame encoding a homing-type endonuclease, and a small group I ribozyme (NaGIR1) inserted into the P6 loop of a second group I ribozyme (NaGIR2). The two ribozymes have different functions in RNA splicing and processing. NaGIR1 is an unusual self-cleaving group I ribozyme responsible for intron processing at two internal sites (IPS1 and IPS2), both close to the 5' end of the open reading frame. This processing is hypothesized to lead to formation of a messenger RNA for the endonuclease. Structurally, NaGIR2 is a typical group IC1 ribozyme, catalyzing intron excision and exon ligation reactions. NaGIR2 is responsible for circularization of the excised intron, a reaction that generates full-length RNA circles of wild-type intron. Although it is only distantly related in primary sequence, NaSSU1 RNA has a predicted organization and function very similar to that of the mobile group I intron DiSSU1 of Didymium, the only other group I intron known to encode two ribozymes. We propose that these twin-ribozyme introns define a distinct category of group I introns with a conserved structural organization and function.