SegH and Hef: two novel homing endonucleases whose genes replace the mobC and mobE genes in several T4-related phages

Nuclease genes occupy boundaries of genetic exchange between bacteriophages

Homing endonuclease genes (HEGs) are ubiquitous selfish elements that generate targeted double-stranded DNA breaks, facilitating the recombination of the HEG DNA sequence into the break site and contributing to the evolutionary dynamics of HEG-encoding genomes. Bacteriophages (phages) are well-documented to carry HEGs, with the paramount characterization of HEGs being focused on those encoded by coliphage T4. Recently, it has been observed that the highly sampled vibriophage, ICP1, is similarly enriched with HEGs distinct from T4's. Here, we examined the HEGs encoded by ICP1 and diverse phages, proposing HEG-driven mechanisms that contribute to phage evolution. Relative to ICP1 and T4, we found a variable distribution of HEGs across phages, with HEGs frequently encoded proximal to or within essential genes. We identified large regions (> 10kb) of high nucleotide identity flanked by HEGs, deemed HEG islands, which we hypothesize to be mobilized by the activity of flanking HEGs. Finally, we found examples of domain swapping between phage-encoded HEGs and genes encoded by other phages and phage satellites. We anticipate that HEGs have a larger impact on the evolutionary trajectory of phages than previously appreciated and that future work investigating the role of HEGs in phage evolution will continue to highlight these observations.

Nuclease genes occupy boundaries of genetic exchange between bacteriophages

Homing endonuclease genes (HEGs) are ubiquitous selfish elements that generate targeted double-stranded DNA breaks, facilitating the recombination of the HEG DNA sequence into the break site and contributing to the evolutionary dynamics of HEG-encoding genomes. Bacteriophages (phages) are well-documented to carry HEGs, with the paramount characterization of HEGs being focused on those encoded by coliphage T4. Recently, it has been observed that the highly sampled vibriophage, ICP1, is similarly enriched with HEGs distinct from T4’s. Here, we examined the HEGs encoded by ICP1 and diverse phages, proposing HEG-driven mechanisms that contribute to phage evolution. Relative to ICP1 and T4, we found a variable distribution of HEGs across phages, with HEGs frequently encoded proximal to or within essential genes. We identified large regions (> 10kb) of high nucleotide identity flanked by HEGs, deemed HEG islands, which we hypothesize to be mobilized by the activity of flanking HEGs. Finally, we found examples of domain swapping between phage-encoded HEGs and genes encoded by other phages and phage satellites. We anticipate that HEGs have a larger impact on the evolutionary trajectory of phages than previously appreciated and that future work investigating the role of HEGs in phage evolution will continue to highlight these observations.

Comparative genomics of Cp8viruses with special reference to Campylobacter phage vB_CjeM_los1, isolated from a slaughterhouse in Ireland

Campylobacter phage vB_CjeM_Los1 was recently isolated from a slaughterhouse in the Republic of Ireland using the host Campylobacter jejuni subsp. jejuni PT14, and full-genome sequencing and annotation were performed. The genome was found to be 134,073 bp in length and to contain 169 predicted open reading frames. Transmission electron microscopy images of vB_CjeM_Los1 revealed that it belongs to the family Myoviridae, with tail fibres observed in both extended and folded conformations, as seen in T4. The genome size and morphology of vB_CjeM_Los1 suggest that it belongs to the genus Cp8virus, and seven other Campylobacter phages with similar size characteristics have also been fully sequenced. In this work, comparative studies were performed in relation to genomic rearrangements and conservation within each of the eight genomes. None of the eight genomes were found to have undergone internal rearrangements, and their sequences retained more than 98% identity with one another despite the widespread geographical distribution of each phage. Whole-genome phylogenetics were also performed, and clades were shown to be representative of the differing number of tRNAs present in each phage. This may be an indication of lineages within the genus, despite their striking homology.

Evolutionary and biogeographical implications of degraded LAGLIDADG endonuclease functionality and group I intron occurrence in stony corals (Scleractinia) and mushroom corals (Corallimorpharia)

Group I introns and homing endonuclease genes (HEGs) are mobile genetic elements, capable of invading target sequences in intron-less genomes. LAGLIDADG HEGs are the largest family of endonucleases, playing a key role in the mobility of group I introns in a process known as ‘homing’. Group I introns and HEGs are rare in metazoans, and can be mainly found inserted in the COXI gene of some sponges and cnidarians, including stony corals (Scleractinia) and mushroom corals (Corallimorpharia). Vertical and horizontal intron transfer mechanisms have been proposed as explanations for intron occurrence in cnidarians. However, the central role of LAGLIDADG motifs in intron mobility mechanisms remains poorly understood. To resolve questions regarding the evolutionary origin and distribution of group I introns and HEGs in Scleractinia and Corallimorpharia, we examined intron/HEGs sequences within a comprehensive phylogenetic framework. Analyses of LAGLIDADG motif conservation showed a high degree of degradation in complex Scleractinia and Corallimorpharia. Moreover, the two motifs lack the respective acidic residues necessary for metal-ion binding and catalysis, potentially impairing horizontal intron mobility. In contrast, both motifs are highly conserved within robust Scleractinia, indicating a fully functional endonuclease capable of promoting horizontal intron transference. A higher rate of non-synonymous substitutions (K_a) detected in the HEGs of complex Scleractinia and Corallimorpharia suggests degradation of the HEG, whereas lower K_a rates in robust Scleractinia are consistent with a scenario of purifying selection. Molecular-clock analyses and ancestral inference of intron type indicated an earlier intron insertion in complex Scleractinia and Corallimorpharia in comparison to robust Scleractinia. These findings suggest that the lack of horizontal intron transfers in the former two groups is related to an age-dependent degradation of the endonuclease activity. Moreover, they also explain the peculiar geographical patterns of introns in stony and mushroom corals.

Characterization of two polyvalent phages infecting Enterobacteriaceae

Bacteriophages display remarkable genetic diversity and host specificity. In this study, we explore phages infecting bacterial strains of the Enterobacteriaceae family because of their ability to infect related but distinct hosts. We isolated and characterized two novel virulent phages, SH6 and SH7, using a strain of Shigella flexneri as host bacterium. Morphological and genomic analyses revealed that phage SH6 belongs to the T1virus genus of the Siphoviridae family. Conversely, phage SH7 was classified in the T4virus genus of the Myoviridae family. Phage SH6 had a short latent period of 16 min and a burst size of 103 ± 16 PFU/infected cell while the phage SH7 latent period was 23 min with a much lower burst size of 26 ± 5 PFU/infected cell. Moreover, phage SH6 was sensitive to acidic conditions (pH < 5) while phage SH7 was stable from pH 3 to 11 for 1 hour. Of the 35 bacterial strains tested, SH6 infected its S. flexneri host strain and 8 strains of E. coli. Phage SH7 lysed additionally strains of E. coli O157:H7, Salmonella Paratyphi, and Shigella dysenteriae. The broader host ranges of these two phages as well as their microbiological properties suggest that they may be useful for controlling bacterial populations.

Comparative Functional Genomic Analysis of Two Vibrio Phages Reveals Complex Metabolic Interactions with the Host Cell

Sequencing and annotation was performed for two large double stranded DNA bacteriophages, φGrn1 and φSt2 of the Myoviridae family, considered to be of great interest for phage therapy against Vibrios in aquaculture live feeds. In addition, phage–host metabolic interactions and exploitation was studied by transcript profiling of selected viral and host genes. Comparative genomic analysis with other large Vibrio phages was also performed to establish the presence and location of homing endonucleases highlighting distinct features for both phages. Phylogenetic analysis revealed that they belong to the “schizoT4like” clade. Although many reports of newly sequenced viruses have provided a large set of information, basic research related to the shift of the bacterial metabolism during infection remains stagnant. The function of many viral protein products in the process of infection is still unknown. Genome annotation identified the presence of several viral open reading frames (ORFs) participating in metabolism, including a Sir2/cobB (sirtuin) protein and a number of genes involved in auxiliary NAD⁺ and nucleotide biosynthesis, necessary for phage DNA replication. Key genes were subsequently selected for detail study of their expression levels during infection. This work suggests a complex metabolic interaction and exploitation of the host metabolic pathways and biochemical processes, including a possible post-translational protein modification, by the virus during infection.

Genomic analysis of Staphylococcus phage Stau2 isolated from medical specimen

Stau2 is a lytic myophage of Staphylococcus aureus isolated from medical specimen. Exhibiting a broad host range against S. aureus clinical isolates, Stau2 is potentially useful for topical phage therapy or as an additive in food preservation. In this study, Stau2 was firstly revealed to possess a circularly permuted linear genome of 133,798 bp, with low G + C content, containing 146 open reading frames, but encoding no tRNA. The genome is organized into several modules containing genes for packaging, structural proteins, replication/transcription and host-cell-lysis, with the structural proteins and DNA polymerase modules being organized similarly to that in Twort-like phages of Staphylococcus. With the encoded DNA replication genes, Stau2 can possibly use its own system for replication. In addition, analysis in silico found several introns in seven genes, including those involved in DNA metabolism, packaging, and structure, while one of them (helicase gene) is experimentally confirmed to undergo splicing. Furthermore, phylogenetic analysis suggested Stau2 to be most closely related to Staphylococcus phages SA11 and Remus, members of Twort-like phages. The results of sodium dodecyl sulfate polyacrylamide gel electrophoresis showed 14 structural proteins of Stau2 and N-terminal sequencing identified three of them. Importantly, this phage does not encode any proteins which are known or suspected to be involved in toxicity, pathogenicity, or antibiotic resistance. Therefore, further investigations of feasible therapeutic application of Stau2 are needed.

Mapping homing endonuclease cleavage sites using in vitro generated protein

Mapping the precise position of endonucleolytic cleavage sites is a fundamental experimental technique used to describe the function of a homing endonuclease. However, these proteins are often recalcitrant to cloning and over-expression in biological systems because of toxicity induced by spurious DNA cleavage events. In this chapter we outline the steps to successfully express a homing endonuclease in vitro and use this product in nucleotide-resolution cleavage assays.

A bioinformatic analysis of ribonucleotide reductase genes in phage genomes and metagenomes

Background

Ribonucleotide reductase (RNR), the enzyme responsible for the formation of deoxyribonucleotides from ribonucleotides, is found in all domains of life and many viral genomes. RNRs are also amongst the most abundant genes identified in environmental metagenomes. This study focused on understanding the distribution, diversity, and evolution of RNRs in phages (viruses that infect bacteria). Hidden Markov Model profiles were used to analyze the proteins encoded by 685 completely sequenced double-stranded DNA phages and 22 environmental viral metagenomes to identify RNR homologs in cultured phages and uncultured viral communities, respectively.

Results

RNRs were identified in 128 phage genomes, nearly tripling the number of phages known to encode RNRs. Class I RNR was the most common RNR class observed in phages (70%), followed by class II (29%) and class III (28%). Twenty-eight percent of the phages contained genes belonging to multiple RNR classes. RNR class distribution varied according to phage type, isolation environment, and the host’s ability to utilize oxygen. The majority of the phages containing RNRs are Myoviridae (65%), followed by Siphoviridae (30%) and Podoviridae (3%). The phylogeny and genomic organization of phage and host RNRs reveal several distinct evolutionary scenarios involving horizontal gene transfer, co-evolution, and differential selection pressure. Several putative split RNR genes interrupted by self-splicing introns or inteins were identified, providing further evidence for the role of frequent genetic exchange. Finally, viral metagenomic data indicate that RNRs are prevalent and highly dynamic in uncultured viral communities, necessitating future research to determine the environmental conditions under which RNRs provide a selective advantage.

Conclusions

This comprehensive study describes the distribution, diversity, and evolution of RNRs in phage genomes and environmental viral metagenomes. The distinct distributions of specific RNR classes amongst phages, combined with the various evolutionary scenarios predicted from RNR phylogenies suggest multiple inheritance sources and different selective forces for RNRs in phages. This study significantly improves our understanding of phage RNRs, providing insight into the diversity and evolution of this important auxiliary metabolic gene as well as the evolution of phages in response to their bacterial hosts and environments.

T4-related bacteriophage LIMEstone isolates for the control of soft rot on potato caused by 'Dickeya solani'

The bacterium 'Dickeya solani', an aggressive biovar 3 variant of Dickeya dianthicola, causes rotting and blackleg in potato. To control this pathogen using bacteriophage therapy, we isolated and characterized two closely related and specific bacteriophages, vB_DsoM_LIMEstone1 and vB_DsoM_LIMEstone2. The LIMEstone phages have a T4-related genome organization and share DNA similarity with Salmonella phage ViI. Microbiological and molecular characterization of the phages deemed them suitable and promising for use in phage therapy. The phages reduced disease incidence and severity on potato tubers in laboratory assays. In addition, in a field trial of potato tubers, when infected with 'Dickeya solani', the experimental phage treatment resulted in a higher yield. These results form the basis for the development of a bacteriophage-based biocontrol of potato plants and tubers as an alternative for the use of antibiotics.

Activity, specificity and structure of I-Bth0305I: a representative of a new homing endonuclease family

Novel family of putative homing endonuclease genes was recently discovered during analyses of metagenomic and genomic sequence data. One such protein is encoded within a group I intron that resides in the recA gene of the Bacillus thuringiensis 0305ϕ8–36 bacteriophage. Named I-Bth0305I, the endonuclease cleaves a DNA target in the uninterrupted recA gene at a position immediately adjacent to the intron insertion site. The enzyme displays a multidomain, homodimeric architecture and footprints a DNA region of ∼60 bp. Its highest specificity corresponds to a 14-bp pseudopalindromic sequence that is directly centered across the DNA cleavage site. Unlike many homing endonucleases, the specificity profile of the enzyme is evenly distributed across much of its target site, such that few single base pair substitutions cause a significant decrease in cleavage activity. A crystal structure of its C-terminal domain confirms a nuclease fold that is homologous to very short patch repair (Vsr) endonucleases. The domain architecture and DNA recognition profile displayed by I-Bth0305I, which is the prototype of a homing lineage that we term the ‘EDxHD’ family, are distinct from previously characterized homing endonucleases.

Campylobacter jejuni group III phage CP81 contains many T4-like genes without belonging to the T4-type phage group: implications for the evolution of T4 phages

CP81 is a virulent Campylobacter group III phage whose linear genome comprises 132,454 bp. At the nucleotide level, CP81 differs from other phages. However, a number of its structural and replication/recombination proteins revealed a relationship to the group II Campylobacter phages CP220/CPt10 and to T4-type phages. Unlike the T4-related phages, the CP81 genome does not contain conserved replication and virion modules. Instead, the respective genes are scattered throughout the phage genome. Moreover, most genes for metabolic enzymes of CP220/CPt10 are lacking in CP81. On the other hand, the CP81 genome contains nine similar genes for homing endonucleases which may be involved in the attrition of the conserved gene order for the virion core genes of T4-type phages. The phage apparently possesses an unusual modification of C or G bases. Efficient cleavage of its DNA was only achieved with restriction enzymes recognizing pure A/T sites. Uncommonly, phenol extraction leads to a significant loss of CP81 DNA from the aqueous layer, a property not yet described for other phages belonging to the T4 superfamily.

Mobile DNA elements in T4 and related phages

Mobile genetic elements are common inhabitants of virtually every genome where they can exert profound influences on genome structure and function in addition to promoting their own spread within and between genomes. Phage T4 and related phage have long served as a model system for understanding the molecular mechanisms by which a certain class of mobile DNA, homing endonucleases, promote their spread. Homing endonucleases are site-specific DNA endonucleases that initiate mobility by introducing double-strand breaks at defined positions in genomes lacking the endonuclease gene, stimulating repair and recombination pathways that mobilize the endonuclease coding region. In phage T4, homing endonucleases were first discovered as encoded within the self-splicing td, nrdB and nrdD introns of T4. Genomic data has revealed that homing endonucleases are extremely widespread in T-even-like phage, as evidenced by the astounding fact that ~11% of the T4 genome encodes homing endonuclease genes, with most of them located outside of self-splicing introns. Detailed studies of the mobile td intron and its encoded endonuclease, I-TevI, have laid the foundation for genetic, biochemical and structural aspects that regulate the mobility process, and more recently have provided insights into regulation of homing endonuclease function. Here, we summarize the current state of knowledge regarding T4-encoded homing endonucleases, with particular emphasis on the td/I-TevI model system. We also discuss recent progress in the biology of free-standing endonucleases, and present areas of future research for this fascinating class of mobile genetic elements.

Assembly of a fragmented ribonucleotide reductase by protein interaction domains derived from a mobile genetic element

Ribonucleotide reductase (RNR) is a critical enzyme of nucleotide metabolism, synthesizing precursors for DNA replication and repair. In prokaryotic genomes, RNR genes are commonly targeted by mobile genetic elements, including free standing and intron-encoded homing endonucleases and inteins. Here, we describe a unique molecular solution to assemble a functional product from the RNR large subunit gene, nrdA that has been fragmented into two smaller genes by the insertion of mobE, a mobile endonuclease. We show that unique sequences that originated during the mobE insertion and that are present as C- and N-terminal tails on the split NrdA-a and NrdA-b polypeptides, are absolutely essential for enzymatic activity. Our data are consistent with the tails functioning as protein interaction domains to assemble the tetrameric (NrdA-a/NrdA-b)(2) large subunit necessary for a functional RNR holoenzyme. The tails represent a solution distinct from RNA and protein splicing or programmed DNA rearrangements to restore function from a fragmented coding region and may represent a general mechanism to neutralize fragmentation of essential genes by mobile genetic elements.

Phage T4 mobE promotes trans homing of the defunct homing endonuclease I-TevIII

Homing endonucleases are site-specific DNA endonucleases that typically function as mobile genetic elements by introducing a double-strand break (DSB) in genomes that lack the endonuclease, resulting in a unidirectional gene conversion event that mobilizes the homing endonuclease gene and flanking DNA. Here, we characterize phage T4-encoded mobE, a predicted free-standing HNH family homing endonuclease. We show that mobE is promoterless and dependent on upstream transcription for expression, and that an internal intrinsic terminator regulates mobE transcript levels. Crucially, in vivo mapping experiments revealed a MobE-dependent, strand-specific nick in the non-coding strand of the nrdB gene of phage T2. An internal deletion of the predicted HNH catalytic motif of MobE abolishes nicking, and reduces high-frequency inheritance of mobE. Sequence polymorphisms of progeny phage that inherit mobE are consistent with DSB repair pathways. Significantly, we found that mobility of the neighboring I-TevIII, a defunct homing endonuclease encoded within a group I intron interrupting the nrdB gene of phage T4, was dependent on an intact mobE gene. Thus, our data indicate that the stagnant nrdB intron and I-TevIII are mobilized in trans as a consequence of a MobE-dependent gene conversion event, facilitating persistence of genetic elements that have no inherent means of promoting their own mobility.

Genetic insertions and diversification of the PolB-type DNA polymerase (gp43) of T4-related phages

In Escherichia coli phage T4 and many of its phylogenetic relatives, gene 43 consists of a single cistron that encodes a PolB family (PolB-type) DNA polymerase. We describe the divergence of this phage gene and its protein product (gp43) (gene product 43) among 26 phylogenetic relatives of T4 and discuss our observations in the context of diversity among the widely distributed PolB enzymes in nature. In two T4 relatives that grow in Aeromonas salmonicida phages 44RR and 25, gene 43 is fragmented by different combinations of three distinct types of DNA insertion elements: (a) a short intercistronic untranslated sequence (IC-UTS) that splits the polymerase gene into two cistrons, 43A and 43B, corresponding to N-terminal (gp43A) and C-terminal (gp43B) protein products; (b) a freestanding homing endonuclease gene (HEG) inserted between the IC-UTS and the 43B cistron; and (c) a group I intron in the 43B cistron. Phage 25 has all three elements, whereas phage 44RR has only the IC-UTS. We present evidence that (a) the split gene of phage 44RR encodes a split DNA polymerase consisting of a complex between gp43A and gp43B subunits; (b) the putative HEG encodes a double-stranded DNA endonuclease that specifically cleaves intron-free homologues of the intron-bearing 43B site; and (c) the group I intron is a self-splicing RNA. Our results suggest that some freestanding HEGs can mediate the homing of introns that do not encode their own homing enzymes. The results also suggest that different insertion elements can converge on a polB gene and evolve into a single integrated system for lateral transfer of polB genetic material. We discuss the possible pathways for the importation of such insertion elements into the genomes of T4-related phages.

Viruses are essential agents within the roots and stem of the tree of life

In contrast with former definitions of life limited to membrane-bound cellular life forms which feed, grow, metabolise and replicate (i) a role of viruses as genetic symbionts, (ii) along with peripheral phenomena such as cryptobiosis and (iii) the horizontal nature of genetic information acquisition and processing broaden our view of the tree of life. Some researchers insist on the traditional textbook conviction of what is part of the community of life. In a recent review [Moreira, D., Lopez-Garcia, P., 2009. Ten reasons to exclude viruses from the tree of life. Nat. Rev. Microbiol. 7, 306-311.] they assemble four main arguments which should exclude viruses from the tree of life because of their inability to self-sustain and self-replicate, their polyphyly, the cellular origin of their cell-like genes and the volatility of their genomes. In this article we will show that these features are not coherent with current knowledge about viruses but that viral agents play key roles within the roots and stem of the tree of life.

Biochemical and mutagenic analysis of I-CreII reveals distinct but important roles for both the H-N-H and GIY-YIG motifs

Homing endonucleases typically contain one of four conserved catalytic motifs, and other elements that confer tight DNA binding. I-CreII, which catalyzes homing of the Cr.psbA4 intron, is unusual in containing two potential catalytic motifs, H-N-H and GIY-YIG. Previously, we showed that cleavage by I-CreII leaves ends (2-nt 3' overhangs) that are characteristic of GIY-YIG endonucleases, yet it has a relaxed metal requirement like H-N-H enzymes. Here we show that I-CreII can bind DNA without an added metal ion, and that it binds as a monomer, akin to GIY-YIG enzymes. Moreover, cleavage of supercoiled DNA, and estimates of strand-specific cleavage rates, suggest that I-CreII uses a sequential cleavage mechanism. Alanine substitution of a number of residues in the GIY-YIG motif, however, did not block cleavage activity, although DNA binding was substantially reduced in several variants. Substitution of conserved histidines in the H-N-H motif resulted in variants that did not promote DNA cleavage, but retained high-affinity DNA binding-thus identifying it as the catalytic motif. Unlike the non-specific H-N-H colicins, however; substitution of the conserved asparagine substantially reduced DNA binding (though not the ability to promote cleavage). These results indicate that, in I-CreII, two catalytic motifs have evolved to play important roles in specific DNA binding. The data also indicate that only the H-N-H motif has retained catalytic ability.

A likely pathway for formation of mobile group I introns

Mobile group I introns are RNA splicing elements that have been invaded by endonuclease genes. These endonucleases facilitate intron mobility by a unidirectional, duplicative gene-conversion process known as homing [1]. Survival of the invading endonuclease depends upon its ability to promote intron mobility. Therefore, the endonuclease must either quickly change its cleavage specificity to match the site of intron insertion, or it must already be preadapted to cleave this sequence. Here we show that the group I intron in the DNA polymerase gene of T7-like bacteriophage PhiI is mobile, dependent upon its intronic HNH homing endonuclease gene, I-TslI. We also show that gene 5.3 of phage T3, located adjacent to its intronless DNA polymerase gene, is a homologous homing endonuclease gene whose protein product initiates efficient spread of gene 5.3 into empty sites in related phages. Both of these endonucleases cleave intronless DNA polymerase genes at identical positions. This shared feature between an intronic and free-standing endonuclease is unprecedented. Based on this evidence, we propose that introns and their homing endonucleases evolve separately to target the same highly conserved sequences, uniting afterwards to create a composite mobile element.

An RNA hairpin sequesters the ribosome binding site of the homing endonuclease mobE gene

Previous transcript mapping of the bacteriophage Aeh1 nrd operon revealed a predicted RNA hairpin upstream of the homing endonuclease mobE gene. We enzymatically mapped the hairpin, showing that the mobE ribosome binding site is sequestered. Cloning of the hairpin upstream of lacZ resulted in reduced beta-galactosidase activity, consistent with translational regulation.

Amino acid residues in the GIY-YIG endonuclease II of phage T4 affecting sequence recognition and binding as well as catalysis

Phage T4 endonuclease II (EndoII), a GIY-YIG endonuclease lacking a carboxy-terminal DNA-binding domain, was subjected to site-directed mutagenesis to investigate roles of individual amino acids in substrate recognition, binding, and catalysis. The structure of EndoII was modeled on that of UvrC. We found catalytic roles for residues in the putative catalytic surface (G49, R57, E118, and N130) similar to those described for I-TevI and UvrC; in addition, these residues were found to be important for substrate recognition and binding. The conserved glycine (G49) and arginine (R57) were essential for normal sequence recognition. Our results are in agreement with a role for these residues in forming the DNA-binding surface and exposing the substrate scissile bond at the active site. The conserved asparagine (N130) and an adjacent proline (P127) likely contribute to positioning the catalytic domain correctly. Enzymes in the EndoII subfamily of GIY-YIG endonucleases share a strongly conserved middle region (MR, residues 72 to 93, likely helical and possibly substituting for heterologous helices in I-TevI and UvrC) and a less strongly conserved N-terminal region (residues 12 to 24). Most of the conserved residues in these two regions appeared to contribute to binding strength without affecting the mode of substrate binding at the catalytic surface. EndoII K76, part of a conserved NUMOD3 DNA-binding motif of homing endonucleases found to overlap the MR, affected both sequence recognition and catalysis, suggesting a more direct involvement in positioning the substrate. Our data thus suggest roles for the MR and residues conserved in GIY-YIG enzymes in recognizing and binding the substrate.

Convergent mechanisms of genome evolution of large and giant DNA viruses

We have taken advantage of the availability of the genome sequences of a collection of large and giant viruses infecting bacteria (T4 family) and eukaryotes (NCLDV group) to assess some of the evolutionary forces which might have shaped their genomes. Despite having apparently different ancestors, these two groups of viruses are affected by convergent evolutionary forces. Both types of virus probably originated from a simple and ancient viral ancestor with a small subset of 30-35 genes encoding replication and structural proteins. The genome size and diversity of the descendants most likely grew progressively by (i) lineage-specific gene duplications, (ii) lateral gene transfers of cellular genes and (iii) accretion of diverse families of mobile genetic elements. These results argue against the hypothesis that giant viruses derive from a regressive cell.

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.

A functional homing endonuclease in the Bacillus anthracis nrdE group I intron

The essential Bacillus anthracis nrdE gene carries a self-splicing group I intron with a putative homing endonuclease belonging to the GIY-YIG family. Here, we show that the nrdE pre-mRNA is spliced and that the homing endonuclease cleaves an intronless nrdE gene 5 nucleotides (nt) upstream of the intron insertion site, producing 2-nt 3' extensions. We also show that the sequence required for efficient cleavage spans at least 4 bp upstream and 31 bp downstream of the cleaved coding strand. The position of the recognition sequence in relation to the cleavage position is as expected for a GIY-YIG homing endonuclease. Interestingly, nrdE genes from several other Bacillaceae were also susceptible to cleavage, with those of Bacillus cereus, Staphylococcus epidermidis (nrdE1), B. anthracis, and Bacillus thuringiensis serovar konkukian being better substrates than those of Bacillus subtilis, Bacillus lichenformis, and S. epidermidis (nrdE2). On the other hand, nrdE genes from Lactococcus lactis, Escherichia coli, Salmonella enterica serovar Typhimurium, and Corynebacterium ammoniagenes were not cleaved. Intervening sequences (IVSs) residing in protein-coding genes are often found in enzymes involved in DNA metabolism, and the ribonucleotide reductase nrdE gene is a frequent target for self-splicing IVSs. A comparison of nrdE genes from seven gram-positive low-G+C bacteria, two bacteriophages, and Nocardia farcinica showed five different insertion sites for self-splicing IVSs within the coding region of the nrdE gene.

Multiple controls regulate the expression of mobE, an HNH homing endonuclease gene embedded within a ribonucleotide reductase gene of phage Aeh1

Mobile genetic elements have the potential to influence the expression of genes surrounding their insertion site upon invasion of a genome. Here, we examine the transcriptional organization of a ribonucleotide reductase operon (nrd) that has been invaded by an HNH family homing endonuclease, mobE. In Aeromonas hydrophila phage Aeh1, mobE has inserted into the large-subunit gene (nrdA) of aerobic ribonucleotide reductase (RNR), splitting it into two smaller genes, nrdA-a and nrdA-b. This gene organization differs from that in phages T4, T6, RB2, RB3, RB15, and LZ7, where mobE is inserted in the nrdA-nrdB intergenic region. We present evidence that the expression of Aeh1 mobE is regulated by transcriptional, posttranscriptional, and translational controls. An Aeh1-specific late promoter drives expression of mobE, but strikingly the mobE transcript is processed internally at an RNase E-like site. We also identified a putative stem-loop structure upstream of mobE that sequesters the mobE ribosome binding site, presumably acting to down regulate MobE translation. Moreover, our transcriptional analyses indicate that the surrounding nrd genes of phage Aeh1 are expressed by a different strategy than are the corresponding phage T4 genes and that transcriptional readthrough is the only mechanism by which the promoterless Aeh1 nrdB gene is expressed. We suggest that the occurrence of multiple layers of control to limit the expression of mobE to late in the Aeh1 infection cycle is an adaptation of Aeh1 to reduce any effects on expression of the surrounding nrd genes early in phage infection when RNR function is critical.

Insertion of a homing endonuclease creates a genes-in-pieces ribonucleotide reductase that retains function

In bacterial and phage genomes, coding regions are sometimes interrupted by self-splicing introns or inteins, which can encode mobility-promoting homing endonucleases. Homing endonuclease genes are also found free-standing (not intron- or intein-encoded) in phage genomes where they are inserted in intergenic regions. One example is the HNH family endonuclease, mobE, inserted between the large (nrdA) and small (nrdB) subunit genes of aerobic ribonucleotide reductase (RNR) of T-even phages T4, RB2, RB3, RB15, and LZ7. Here, we describe an insertion of mobE into the nrdA gene of Aeromonas hydrophila phage Aeh1. The insertion creates a unique genes-in-pieces arrangement, where nrdA is split into two independent genes, nrdA-a and nrdA-b, each encoding cysteine residues that correspond to the active-site residues of uninterrupted NrdA proteins. Remarkably, the mobE insertion does not inactivate NrdA function, although the insertion is not a self-splicing intron or intein. We copurified the NrdA-a, NrdA-b, and NrdB proteins as complex from Aeh1-infected cells and also showed that a reconstituted complex has RNR activity. Class I RNR activity in phage Aeh1 is thus assembled from separate proteins that interact to form a composite active site, demonstrating that the mobE insertion is phenotypically neutral in that its presence as an intervening sequence does not disrupt the function of the surrounding gene.

Plasticity of the gene functions for DNA replication in the T4-like phages

We have completely sequenced and annotated the genomes of several relatives of the bacteriophage T4, including three coliphages (RB43, RB49 and RB69), three Aeromonas salmonicida phages (44RR2.8t, 25 and 31) and one Aeromonas hydrophila phage (Aeh1). In addition, we have partially sequenced and annotated the T4-like genomes of coliphage RB16 (a close relative of RB43), A. salmonicida phage 65, Acinetobacter johnsonii phage 133 and Vibrio natriegens phage nt-1. Each of these phage genomes exhibited a unique sequence that distinguished it from its relatives, although there were examples of genomes that are very similar to each other. As a group the phages compared here diverge from one another by several criteria, including (a) host range, (b) genome size in the range between approximately 160 kb and approximately 250 kb, (c) content and genetic organization of their T4-like genes for DNA metabolism, (d) mutational drift of the predicted T4-like gene products and their regulatory sites and (e) content of open-reading frames that have no counterparts in T4 or other known organisms (novel ORFs). We have observed a number of DNA rearrangements of the T4 genome type, some exhibiting proximity to putative homing endonuclease genes. Also, we cite and discuss examples of sequence divergence in the predicted sites for protein-protein and protein-nucleic acid interactions of homologues of the T4 DNA replication proteins, with emphasis on the diversity in sequence, molecular form and regulation of the phage-encoded DNA polymerase, gp43. Five of the sequenced phage genomes are predicted to encode split forms of this polymerase. Our studies suggest that the modular construction and plasticity of the T4 genome type and several of its replication proteins may offer resilience to mutation, including DNA rearrangements, and facilitate the adaptation of T4-like phages to different bacterial hosts in nature.