Introns and intein coding sequence in the ribonucleotide reductase genes of Bacillus subtilis temperate bacteriophage SPbeta

Genomic diversity and comprehensive taxonomical classification of 61 Bacillus subtilis group member infecting bacteriophages, and the identification of ortholog taxonomic signature genes

BACKGROUND

Despite the applications of Bacillus subtilis group species in various sectors, limited information is available regarding their phages. Here, 61 B. subtilis group species-infecting phages (BSPs) were studied for their taxonomic classification considering the genome-size, genomic diversity, and the host, followed by the identification of orthologs taxonomic signature genes.

RESULTS

BSPs have widely ranging genome sizes that can be bunched into groups to demonstrate correlations to family and subfamily classifications. Comparative analysis re-confirmed the existing, BSPs-containing 14 genera and 21 species and displayed inter-genera similarities within existing subfamilies. Importantly, it also revealed the need for the creation of new taxonomic classifications, including 28 species, nine genera, and two subfamilies (New subfamily1 and New subfamily2) to accommodate inter-genera relatedness. Following pangenome analysis, no ortholog shared by all BSPs was identified, while orthologs, namely, the tail fibers/spike proteins and poly-gamma-glutamate hydrolase, that are shared by more than two-thirds of the BSPs were identified. More importantly, major capsid protein (MCP) type I, MCP type II, MCP type III and peptidoglycan binding proteins that are distinctive orthologs for Herelleviridae, Salasmaviridae, New subfamily1, and New subfamily2, respectively, were identified and analyzed which could serve as signatures to distinguish BSP members of the respective taxon.

CONCLUSIONS

In this study, we show the genomic diversity and propose a comprehensive classification of 61 BSPs, including the proposition for the creation of two new subfamilies, followed by the identification of orthologs taxonomic signature genes, potentially contributing to phage taxonomy.

Phage Genome Annotation: Where to Begin and End

With the renewed interest in phage research, coupled with the rising accessibility to affordable sequencing, ever increasing numbers of phage genomes are being sequenced. Therefore, there is an increased need to assemble and annotate phage genomes. There is a plethora of tools and platforms that allow phage genomes to be assembled and annotated. The choice of tools can often be bewildering for those new to phage genome assembly. Here we provide an overview of the assembly and annotation process from obtaining raw reads to genome submission, with worked examples, providing those new to genome assembly and annotation with a guided pathway to genome submission. We focus on the use of open access tools that can be incorporated into workflows to allow easy repetition of steps, highlighting multiple tools that can be used and common pitfalls that may occur.

The life cycle of SPβ and related phages

Phages are viruses of bacteria and are the smallest and most common biological entities in the environment. They can reproduce immediately after infection or integrate as a prophage into their host genome. SPβ is a prophage of the Gram-positive model organism Bacillus subtilis 168, and it has been known for more than 50 years. It is sensitive to dsDNA damage and is induced through exposure to mitomycin C or UV radiation. When induced from the prophage, SPβ requires 90 min to produce and release about 30 virions. Genomes of sequenced related strains range between 128 and 140 kb, and particle-packed dsDNA exhibits terminal redundancy. Formed particles are of the Siphoviridae morphotype. Related isolates are known to infect other B. subtilis clade members. When infecting a new host, SPβ presumably follows a two-step strategy, adsorbing primarily to teichoic acid and secondarily to a yet unknown factor. Once in the host, SPβ-related phages pass through complex lysis-lysogeny decisions and either enter a lytic cycle or integrate as a dormant prophage. As prophages, SPβ-related phages integrate at the host chromosome's replication terminus, and frequently into the spsM or kamA gene. As a prophage, it imparts additional properties to its host via phage-encoded proteins. The most notable of these functional proteins is sublancin 168, which is used as a molecular weapon by the host and ensures prophage maintenance. In this review, we summarise the existing knowledge about the biology of the phage regarding its life cycle and discuss its potential as a research object.

Genomic Characterization of Cyanophage vB_AphaS-CL131 Infecting Filamentous Diazotrophic Cyanobacterium Aphanizomenon flos-aquae Reveals Novel Insights into Virus-Bacterium Interactions

While filamentous cyanobacteria play a crucial role in food web dynamics and biogeochemical cycling of many aquatic ecosystems around the globe, the knowledge regarding the phages infecting them is limited. Here, we describe the complete genome of the virulent cyanophage vB_AphaS-CL131 (here, CL 131), a Siphoviridae phage that infects the filamentous diazotrophic bloom-forming cyanobacterium Aphanizomenon flos-aquae in the brackish Baltic Sea. CL 131 features a 112,793-bp double-stranded DNA (dsDNA) genome encompassing 149 putative open reading frames (ORFs), of which the majority (86%) lack sequence homology to genes with known functions in other bacteriophages or bacteria. Phylogenetic analysis revealed that CL 131 possibly represents a new evolutionary lineage within the group of cyanophages infecting filamentous cyanobacteria, which form a separate cluster from phages infecting unicellular cyanobacteria. CL 131 encodes a putative type V-U2 CRISPR-Cas system with one spacer (out of 10) targeting a DNA primase pseudogene in a cyanobacterium and a putative type II toxin-antitoxin system, consisting of a GNAT family N-acetyltransferase and a protein of unknown function containing the PRK09726 domain (characteristic of HipB antitoxins). Comparison of CL 131 proteins to reads from Baltic Sea and other available fresh- and brackish-water metagenomes and analysis of CRISPR-Cas arrays in publicly available A. flos-aquae genomes demonstrated that phages similar to CL 131 are present and dynamic in the Baltic Sea and share a common history with their hosts dating back at least several decades. In addition, different CRISPR-Cas systems within individual A. flos-aquae genomes targeted several sequences in the CL 131 genome, including genes related to virion structure and morphogenesis. Altogether, these findings revealed new genomic information for exploring viral diversity and provide a model system for investigation of virus-host interactions in filamentous cyanobacteria.IMPORTANCE The genomic characterization of novel cyanophage vB_AphaS-CL131 and the analysis of its genomic features in the context of other viruses, metagenomic data, and host CRISPR-Cas systems contribute toward a better understanding of aquatic viral diversity and distribution in general and of brackish-water cyanophages infecting filamentous diazotrophic cyanobacteria in the Baltic Sea in particular. The results of this study revealed previously undescribed features of cyanophage genomes (e.g., self-excising intein-containing putative dCTP deaminase and putative cyanophage-encoded CRISPR-Cas and toxin-antitoxin systems) and can therefore be used to predict potential interactions between bloom-forming cyanobacteria and their cyanophages.

Metagenome analysis of Russian and Georgian Pyophage cocktails and a placebo-controlled safety trial of single phage versus phage cocktail in healthy Staphylococcus aureus carriers

Bacteriophage therapy is a commonly used treatment for Staphylococcus aureus infections in countries of the former Soviet Union, using both single phages and phage cocktails. The scarce data available on Eastern phage cocktails prompted an investigation into commercially-available Pyophage cocktails from two different manufacturers used to treat skin and wound infections. Comparison of the metagenomic composition of two Pyophage products from Georgia and Russia revealed substantial differences in phage-types targeting Escherichia, Enterococcus, Salmonella, Pseudomonas aeruginosa and Proteus, therefore indicating multiple strategies for composing phage cocktails against these bacterial pathogens. Closely-related Kayvirus-like Myoviruses were, however, a shared component against S. aureus within all products, except for the inclusion of a secondary S. aureus Podovirus in one Microgen cocktail. Metagenomic analysis also revealed the presence of several probable prophage sequences but detected no genetic safety risks in terms of virulence factors or antibiotic resistance genes. The safety of broad-spectrum cocktails was tested by comparing the effects of nasal and oral exposure to Eliava Pyophage, a monospecies counterpart and placebo in healthy human carriers of S. aureus. The lack of adverse effects in any treatment groups supports the clinical safety of S. aureus phages administered as a single phage or as phage cocktail.

Mobile self-splicing introns and inteins as environmental sensors

Self-splicing introns and inteins are often mobile at the level of the genome. Although these RNA and protein elements, respectively, are generally considered to be selfish parasites, group I and group II introns and inteins can be triggered by environmental cues to splice and/or to mobilize. These cues include stressors such as oxidizing agents, reactive oxygen and nitrogen species, starvation, temperature, osmolarity and DNA damage. Their sensitivity to these stimuli leads to a carefully choreographed dance between the mobile element and its host that is in tune with the cellular environment. This responsiveness to a changing milieu provides strong evidence that these diverse, self-splicing mobile elements have adapted to react to prevailing conditions, to the potential advantage of both the element and its host.

Arabidopsis chloroplast mini-ribonuclease III participates in rRNA maturation and intron recycling

RNase III proteins recognize double-stranded RNA structures and catalyze endoribonucleolytic cleavages that often regulate gene expression. Here, we characterize the functions of RNC3 and RNC4, two Arabidopsis thaliana chloroplast Mini-RNase III-like enzymes sharing 75% amino acid sequence identity. Whereas rnc3 and rnc4 null mutants have no visible phenotype, rnc3/rnc4 (rnc3/4) double mutants are slightly smaller and chlorotic compared with the wild type. In Bacillus subtilis, the RNase Mini-III is integral to 23S rRNA maturation. In Arabidopsis, we observed imprecise maturation of 23S rRNA in the rnc3/4 double mutant, suggesting that exoribonucleases generated staggered ends in the absence of specific Mini-III-catalyzed cleavages. A similar phenotype was found at the 3' end of the 16S rRNA, and the primary 4.5S rRNA transcript contained 3' extensions, suggesting that Mini-III catalyzes several processing events of the polycistronic rRNA precursor. The rnc3/4 mutant showed overaccumulation of a noncoding RNA complementary to the 4.5S-5S rRNA intergenic region, and its presence correlated with that of the extended 4.5S rRNA precursor. Finally, we found rnc3/4-specific intron degradation intermediates that are probable substrates for Mini-III and show that B. subtilis Mini-III is also involved in intron regulation. Overall, this study extends our knowledge of the key role of Mini-III in intron and noncoding RNA regulation and provides important insight into plastid rRNA maturation.

Correction: genomic comparison of 93 Bacillus phages reveals 12 clusters, 14 singletons and remarkable diversity

BACKGROUND

The Bacillus genus of Firmicutes bacteria is ubiquitous in nature and includes one of the best characterized model organisms, B. subtilis, as well as medically significant human pathogens, the most notorious being B. anthracis and B. cereus. As the most abundant living entities on the planet, bacteriophages are known to heavily influence the ecology and evolution of their hosts, including providing virulence factors. Thus, the identification and analysis of Bacillus phages is critical to understanding the evolution of Bacillus species, including pathogenic strains.

RESULTS

Whole genome nucleotide and proteome comparison of the 83 extant, fully sequenced Bacillus phages revealed 10 distinct clusters, 24 subclusters and 15 singleton phages. Host analysis of these clusters supports host boundaries at the subcluster level and suggests phages as vectors for genetic transfer within the Bacillus cereus group, with B. anthracis as a distant member. Analysis of the proteins conserved among these phages reveals enormous diversity and the uncharacterized nature of these phages, with a total of 4,442 protein families (phams) of which only 894 (20%) had a predicted function. In addition, 2,583 (58%) of phams were orphams (phams containing a single member). The most populated phams were those encoding proteins involved in DNA metabolism, virion structure and assembly, cell lysis, or host function. These included several genes that may contribute to the pathogenicity of Bacillus strains.

CONCLUSIONS

This analysis provides a basis for understanding and characterizing Bacillus and other related phages as well as their contributions to the evolution and pathogenicity of Bacillus cereus group bacteria. The presence of sparsely populated clusters, the high ratio of singletons to clusters, and the large number of uncharacterized, conserved proteins confirms the need for more Bacillus phage isolation in order to understand the full extent of their diversity as well as their impact on host evolution.

Integration of genomic and proteomic analyses in the classification of the Siphoviridae family

Using a variety of genomic (BLASTN, ClustalW) and proteomic (Phage Proteomic Tree, CoreGenes) tools we have tackled the taxonomic status of members of the largest bacteriophage family, the Siphoviridae. In all over 400 phages were examined and we were able to propose 39 new genera, comprising 216 phage species, and add 62 species to two previously defined genera (Phic3unalikevirus; L5likevirus) grouping, in total, 390 fully sequenced phage isolates. Many of the remainders are orphans which the Bacterial and Archaeal Viruses Subcommittee of the International Committee on Taxonomy of Viruses (ICTV) chooses not to ascribe genus status at the time being.

Cargo capacity of phages and plasmids and other factors influencing horizontal transfers of prokaryote transposable elements

Horizontal transfer of transposable elements (TEs) plays a key role in prokaryote genome evolution. Most TEs do not encode the enzymatic machinery allowing them to transfer between host cells and it is widely assumed in the literature that horizontal transfer of prokaryote TEs is mediated by other mobile genetic elements such as phages and plasmids. In a recent study, we have shown that phages are less tolerant to insertion sequences (IS, the most frequent class of prokaryote TEs) and therefore have a lower cargo capacity than plasmids. Consequently, while our analysis confirmed the crucial role of plasmids as efficient vehicles of IS horizontal transfer, we concluded that phages are unlikely to efficiently shuttle IS elements between prokaryotes. Here, we discuss whether or not the distribution pattern observed for IS elements in phages and plasmids also holds for other TEs, such as transposons and mobile introns. We also further explore various factors that may impact the relative capacity of phages and plasmids to mediate TE horizontal transfer among prokaryotes.

Enigmatic distribution, evolution, and function of inteins

Inteins are mobile genetic elements capable of self-splicing post-translationally. They exist in all three domains of life including in viruses and bacteriophage, where they have a sporadic distribution even among very closely related species. In this review, we address this anomalous distribution from the point of view of the evolution of the host species as well as the intrinsic features of the inteins that contribute to their genetic mobility. We also discuss the incidence of inteins in functionally important sites of their host proteins. Finally, we describe instances of conditional protein splicing. These latter observations lead us to the hypothesis that some inteins have adapted to become sensors that play regulatory roles within their host protein, to the advantage of the organism in which they reside.

Bacillus subtilis class Ib ribonucleotide reductase: high activity and dynamic subunit interactions

The class Ib ribonucleotide reductase (RNR) isolated from Bacillus subtilis was recently purified as a 1:1 ratio of NrdE (α) and NrdF (β) subunits and determined to have a dimanganic-tyrosyl radical (Mn(III)2-Y·) cofactor. The activity of this RNR and the one reconstituted from recombinantly expressed NrdE and reconstituted Mn(III)2-Y· NrdF using dithiothreitol as the reductant, however, was low (160 nmol min(-1) mg(-1)). The apparent tight affinity between the two subunits, distinct from all class Ia RNRs, suggested that B. subtilis RNR might be the protein that yields to the elusive X-ray crystallographic characterization of an "active" RNR complex. We now report our efforts to optimize the activity of B. subtilis RNR by (1) isolation of NrdF with a homogeneous cofactor, and (2) identification and purification of the endogenous reductant(s). Goal one was achieved using anion exchange chromatography to separate apo-/mismetalated-NrdFs from Mn(III)2-Y· NrdF, yielding enzyme containing 4 Mn and 1 Y·/β2. Goal two was achieved by cloning, expressing, and purifying TrxA (thioredoxin), YosR (a glutaredoxin-like thioredoxin), and TrxB (thioredoxin reductase). The success of both goals increased the specific activity to ~1250 nmol min(-1) mg(-1) using a 1:1 mixture of NrdE:Mn(III)2-Y· NrdF and either TrxA or YosR and TrxB. The quaternary structures of NrdE, NrdF, and NrdE:NrdF (1:1) were characterized by size exclusion chromatography and analytical ultracentrifugation. At physiological concentrations (~1 μM), NrdE is a monomer (α) and Mn(III)2-Y· NrdF is a dimer (β2). A 1:1 mixture of NrdE:NrdF, however, is composed of a complex mixture of structures in contrast to expectations.

Group I introns in Staphylococcus bacteriophages

The abundance of group I introns, intragenic RNA sequences capable of self-splicing, in Gram-positive bacteriophage genomes, is illustrated by various new group I introns recently described in Staphylococcus phage genomes. These introns were found to interrupt DNA metabolism genes as well as late genes. These group I introns often code for homing endonucleases, which promote lateral transfer of group I introns, thereby enabling spread through a population. Homing endonucleases encoded by group I introns in Staphylococcus phage genomes were predicted to belong to the GIY-YIG, LAGLIDADG, HNH or EDxHD family of endonucleases. The group I intron distribution in Staphylococcus phage genomes exemplifies the homology between these introns as well as the encoded endonucleases. Despite several suggested functions, the role of group I introns in bacteriophages remains unclear or might be nonexistent. However, transcriptome analysis might provide additional information to elucidate the possible purpose of group I introns in phage genomes.

Romulus and Remus, two phage isolates representing a distinct clade within the Twortlikevirus genus, display suitable properties for phage therapy applications

The renewed interest in controlling Staphylococcus aureus infections using their natural enemies, bacteriophages, has led to the isolation of a limited number of virulent phages so far. These phages are all members of the Twortlikevirus, displaying little variance. We present two novel closely related (95.9% DNA homology) lytic myoviruses, Romulus and Remus, with double-stranded DNA (dsDNA) genomes of 131,333 bp and 134,643 bp, respectively. Despite their relatedness to Staphylococcus phages K, G1, ISP, and Twort and Listeria phages A511 and P100, Romulus and Remus can be proposed as isolates of a new species within the Twortlikevirus genus. A distinguishing feature for these phage genomes is the unique distribution of group I introns compared to that in other staphylococcal myoviruses. In addition, a hedgehog/intein domain was found within their DNA polymerase genes, and an insertion sequence-encoded transposase exhibits splicing behavior and produces a functional portal protein. From a phage therapy application perspective, Romulus and Remus infected approximately 70% of the tested S. aureus isolates and displayed promising lytic activity against these isolates. Furthermore, both phages showed a rapid initial adsorption and demonstrated biofilm-degrading capacity in a proof-of-concept experiment.

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.

Inhibitory effect of prophage SPβ fragments on phage SP10 ribonucleotide reductase function and its multiplication in Bacillus subtilis

Bacteria have evolved various kinds of defense mechanisms against phage infection and multiplication. Analysis of these mechanisms is important for medical and industrial application of phages as well as for their scientific study. Strains of Bacillus subtilis Marburg strain carrying both nonA and nonB mutations are susceptible to the Bacillus phage SP10. The nonB mutation has been shown to have a compromised intrinsic restriction system. The nonA mutation represents the cured state of prophage SPβ whose genome is 135 kb in length and contains 187 ORFs. For this study we investigated the molecular mechanism behind the inhibitory activity of the wild type nonA function against phage SP10 development. The progression of phage-developmental stages was examined in cells harboring wild type nonA, i.e. prophage SPβ. After phage adsorption and DNA injection into host cells, the synthesis of phage specific mRNA proceeded normally. However, phage DNA synthesis was severely inhibited by some effect of wild type nonA. We thus systematically deleted portions of the prophage SPβ region from the B. subtilis genome and the resultant mutant strains were examined as to whether they still retained sufficient wild type nonA functionality to inhibit SP10 phage development. The SPβ region encompassing the bnrdEF gene, which codes for a putative ribonucleotide reductase (RRase), turned out to be responsible for the wild type nonA function. The phage SP10 possesses its own xnrdE gene coding for a putative RRase that complements the temperature-sensitive mutation of the host RRase gene nrdE. This complementation was blocked by an artificially induced transcription from a non-coding strand of the bnrdEF region. It is thus likely that the transcript from the bnrdEF region of SPβ inhibits ribonucleotide reductase function of SP10, resulting in arrest of DNA synthesis during phage SP10 development.

The genome sequence of Bacillus subtilis subsp. spizizenii W23: insights into speciation within the B. subtilis complex and into the history of B. subtilis genetics

The genome sequence of Bacillus subtilis subsp. spizizenii W23 has been determined. The sequence strongly suggests that W23 is a direct descendant of B. subtilis ATCC 6633. W23 shares a 3.6 Mb core genome with the intensively studied model organism B. subtilis subsp. subtilis 168, and gene order within this core has been strongly conserved. Additionally, the W23 genome has 157 accessory (that is, non-core) genome segments that are not found in 168, while the 168 genome has 141 segments not found in W23. The distribution of sequences similar to these accessory segments among other genomes of the B. subtilis species complex shows that those sequences having entered into the phylogeny of the complex more recently tend to be larger and more AT-rich than those having entered earlier. A simple model can account for these observations, in which parasitic or symbiotic DNAs are transferred into the genome and then are reduced in size and modified in base composition during speciation.

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.

Complete genomic sequence of bacteriophage felix o1

Bacteriophage O1 is a Myoviridae A1 group member used historically for identifying Salmonella. Sequencing revealed a single, linear, 86,155-base-pair genome with 39% average G+C content, 131 open reading frames, and 22 tRNAs. Closest protein homologs occur in Erwinia amylovora phage φEa21-4 and Escherichia coli phage wV8. Proteomic analysis indentified structural proteins: Gp23, Gp36 (major tail protein), Gp49, Gp53, Gp54, Gp55, Gp57, Gp58 (major capsid protein), Gp59, Gp63, Gp64, Gp67, Gp68, Gp69, Gp73, Gp74 and Gp77 (tail fiber). Based on phage-host codon differences, 7 tRNAs could affect translation rate during infection. Introns, holin-lysin cassettes, bacterial toxin homologs and host RNA polymerase-modifying genes were absent.

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.

Comparative genomic analysis of 60 Mycobacteriophage genomes: genome clustering, gene acquisition, and gene size

Mycobacteriophages are viruses that infect mycobacterial hosts. Expansion of a collection of sequenced phage genomes to a total of 60-all infecting a common bacterial host-provides further insight into their diversity and evolution. Of the 60 phage genomes, 55 can be grouped into nine clusters according to their nucleotide sequence similarities, 5 of which can be further divided into subclusters; 5 genomes do not cluster with other phages. The sequence diversity between genomes within a cluster varies greatly; for example, the 6 genomes in Cluster D share more than 97.5% average nucleotide similarity with one another. In contrast, similarity between the 2 genomes in Cluster I is barely detectable by diagonal plot analysis. In total, 6858 predicted open-reading frames have been grouped into 1523 phamilies (phams) of related sequences, 46% of which possess only a single member. Only 18.8% of the phams have sequence similarity to non-mycobacteriophage database entries, and fewer than 10% of all phams can be assigned functions based on database searching or synteny. Genome clustering facilitates the identification of genes that are in greatest genetic flux and are more likely to have been exchanged horizontally in relatively recent evolutionary time. Although mycobacteriophage genes exhibit a smaller average size than genes of their host (205 residues compared with 315), phage genes in higher flux average only 100 amino acids, suggesting that the primary units of genetic exchange correspond to single protein domains.

In silico identification of genes in bacteriophage DNA

One of the most satisfying aspects of a genome sequencing project is the identification of the genes contained within it.These are of two types: those which encode tRNAs and those which produce proteins. After a general introduction on the properties of protein-encoding genes and the utility of the Basic Local Alignment Search Tool (BLASTX) to identify genes through homologs, a variety of tools are discussed by their creators. These include for genome annotation: GeneMark, Artemis, and BASys; and, for genome comparisons: Artemis Comparison Tool (ACT), Mauve, CoreGenes, and GeneOrder.

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.

Cyanobacterial ribosomal RNA genes with multiple, endonuclease-encoding group I introns

BACKGROUND

Group I introns are one of the four major classes of introns as defined by their distinct splicing mechanisms. Because they catalyze their own removal from precursor transcripts, group I introns are referred to as autocatalytic introns. Group I introns are common in fungal and protist nuclear ribosomal RNA genes and in organellar genomes. In contrast, they are rare in all other organisms and genomes, including bacteria.

RESULTS

Here we report five group I introns, each containing a LAGLIDADG homing endonuclease gene (HEG), in large subunit (LSU) rRNA genes of cyanobacteria. Three of the introns are located in the LSU gene of Synechococcus sp. C9, and the other two are in the LSU gene of Synechococcus lividus strain C1. Phylogenetic analyses show that these introns and their HEGs are closely related to introns and HEGs located at homologous insertion sites in organellar and bacterial rDNA genes. We also present a compilation of group I introns with homing endonuclease genes in bacteria.

CONCLUSION

We have discovered multiple HEG-containing group I introns in a single bacterial gene. To our knowledge, these are the first cases of multiple group I introns in the same bacterial gene (multiple group I introns have been reported in at least one phage gene and one prophage gene). The HEGs each contain one copy of the LAGLIDADG motif and presumably function as homodimers. Phylogenetic analysis, in conjunction with their patchy taxonomic distribution, suggests that these intron-HEG elements have been transferred horizontally among organelles and bacteria. However, the mode of transfer and the nature of the biological connections among the intron-containing organisms are unknown.

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.

The Conserved Asparagine in the HNH Motif Serves an Important Structural Role in Metal Finger Endonucleases

The HNH motif is a small nucleic acid binding and cleavage module, widespread in metal finger endonucleases in all life kingdoms. Here we studied a non-specific endonuclease, the nuclease domain of ColE7 (N-ColE7), to decipher the role of the conserved asparagine and histidine residues in the HNH motif. We found, using fluorescence resonance energy transfer (FRET) assays, that the DNA hydrolysis activity of H545 N-ColE7 mutants was completely abolished while activities of N560 and H573 mutants varied from 6.9% to 83.2% of the wild-type activity. The crystal structures of three N-ColE7 mutants in complex with the inhibitor Im7, N560A-Im7, N560D-Im7 and H573A-Im7, were determined at a resolution of 1.9 A to 2.2 A. H573 is responsible for metal ion binding in the wild-type protein, as the zinc ion is still partially associated in the structure of H573A, suggesting that H573 plays a supportive role in metal binding. Both N560A and N560D contain a disordered loop in the HNH motif due to the disruption of the hydrogen bond network surrounding the side-chain of residue 560, and as a result, the imidazole ring of the general base residue H545 is tilted slightly and the scissile phosphate is shifted, leading to the large reductions in hydrolysis activities. These results suggest that the highly conserved asparagine in the HNH motif, in general, plays a structural role in constraining the loop in the metal finger structure and keeping the general base histidine and scissile phosphate in the correct position for DNA hydrolysis.

The Bacillus subtilis nrdEF genes, encoding a class Ib ribonucleotide reductase, are essential for aerobic and anaerobic growth

Ribonucleotide reductases (RNRs) are essential for the biosynthesis of the deoxyribonucleoside triphosphates of DNA. Recently, it was proposed that externally supplied deoxyribonucleosides or DNA is required for the growth of Bacillus subtilis under strict anaerobic conditions (M. J. Folmsbee, M. J. McInerney, and D. P. Nagle, Appl. Environ. Microbiol. 70:5252-5257, 2004). Cultivation of B. subtilis on minimal medium in the presence of oxygen indicators in combination with oxygen electrode measurements and viable cell counting demonstrated that growth occurred under strict anaerobic conditions in the absence of externally supplied deoxyribonucleosides. The nrdEF genes encode the only obvious RNR in B. subtilis. A temperature-sensitive nrdE mutant failed to grow under aerobic and anaerobic conditions, indicating that this oxygen-dependent class I RNR has an essential role under both growth conditions. Aerobic growth and anaerobic growth of the nrdE mutant were rescued by addition of deoxynucleotides. The nrd locus consists of an nrdI-nrdE-nrdF-ymaB operon. The 5' end of the corresponding mRNA revealed transcriptional start sites 45 and 48 bp upstream of the translational start of nrdI. Anaerobic transcription of the operon was found to be dependent on the presence of intact genes for the ResDE two-component redox regulatory system. Two potential ResD binding sites were identified approximately 62 bp (site A) and 50 bp (site B) upstream of the transcriptional start sites by a bioinformatic approach. Only mutation of site B eliminated nrd expression. Aerobic transcription was ResDE independent but required additional promoter elements localized between 88 and 275 bp upstream of the transcriptional start.

Subspecies-specific distribution of intervening sequences in the Bacillus subtilis prophage ribonucleotide reductase genes

A collection of 212 gram-positive bacilli isolated from natural habitats was screened for the presence of intervening sequences (introns and intein-coding sequences) in the SPbeta prophage-related ribonucleotide reductase genes bnrdE and bnrdF. Three novel configurations were identified on the basis of the presence of (i) intervening sequences in bnrdE and bnrdF, and (ii) an ORF in the bnrdE-bnrdF spacer. Analysis of the cell wall genetic determinants as well as of the incorporation of radio-labelled glycerol into cell wall allowed newly and previously identified B. subtilis strains with different configurations of bnrdE/bnrdF intervening sequences to be assigned to one of two subspecies. Strains apparently belonging to the subsp. subtilis contain three intervening sequences many of which are associated with the putative homing endonuclease activity. Strains of the subsp. spizizenii contain only one or two ORF-less group I introns. Introns occupying bnrdF are confined to the subspecies subtilis.

Genomic organization and molecular analysis of virulent bacteriophage 2972 infecting an exopolysaccharide-producing Streptococcus thermophilus strain

The Streptococcus thermophilus virulent pac-type phage 2972 was isolated from a yogurt made in France in 1999. It is a representative of several phages that have emerged with the industrial use of the exopolysaccharide-producing S. thermophilus strain RD534. The genome of phage 2972 has 34,704 bp with an overall G+C content of 40.15%, making it the shortest S. thermophilus phage genome analyzed so far. Forty-four open reading frames (ORFs) encoding putative proteins of 40 or more amino acids were identified, and bioinformatic analyses led to the assignment of putative functions to 23 ORFs. Comparative genomic analysis of phage 2972 with the six other sequenced S. thermophilus phage genomes confirmed that the replication module is conserved and that cos- and pac-type phages have distinct structural and packaging genes. Two group I introns were identified in the genome of 2972. They interrupted the genes coding for the putative endolysin and the terminase large subunit. Phage mRNA splicing was demonstrated for both introns, and the secondary structures were predicted. Eight structural proteins were also identified by N-terminal sequencing and/or matrix-assisted laser desorption ionization-time-of-flight mass spectrometry. Detailed analysis of the putative minor tail proteins ORF19 and ORF21 as well as the putative receptor-binding protein ORF20 showed the following interesting features: (i) ORF19 is a hybrid protein, because it displays significant identity with both pac- and cos-type phages; (ii) ORF20 is unique; and (iii) a protein similar to ORF21 of 2972 was also found in the structure of the cos-type phage DT1, indicating that this structural protein is present in both S. thermophilus phage groups. The implications of these findings for phage classification are discussed.

Algal viruses with distinct intraspecies host specificities include identical intein elements

Heterosigma akashiwo virus (HaV) is a large double-stranded DNA virus infecting the single-cell bloom-forming raphidophyte (golden brown alga) H. akashiwo. A molecular phylogenetic sequence analysis of HaV DNA polymerase showed that it forms a sister group with Phycodnaviridae algal viruses. All 10 examined HaV strains, which had distinct intraspecies host specificities, included an intein (protein intron) in their DNA polymerase genes. The 232-amino-acid inteins differed from each other by no more than a single nucleotide change. All inteins were present at the same conserved position, coding for an active-site motif, which also includes inteins in mimivirus (a very large double-stranded DNA virus of amoebae) and in several archaeal DNA polymerase genes. The HaV intein is closely related to the mimivirus intein, and both are apparently monophyletic to the archaeal inteins. These observations suggest the occurrence of horizontal transfers of inteins between viruses of different families and between archaea and viruses and reveal that viruses might be reservoirs and intermediates in horizontal transmissions of inteins. The homing endonuclease domain of the HaV intein alleles is mostly deleted. The mechanism keeping their sequences basically identical in HaV strains specific for different hosts is yet unknown. One possibility is that rapid and local changes in the HaV genome change its host specificity. This is the first report of inteins found in viruses infecting eukaryotic algae.

Insights on evolution of virulence and resistance from the complete genome analysis of an early methicillin-resistant Staphylococcus aureus strain and a biofilm-producing methicillin-resistant Staphylococcus epidermidis strain

Staphylococcus aureus is an opportunistic pathogen and the major causative agent of numerous hospital- and community-acquired infections. Staphylococcus epidermidis has emerged as a causative agent of infections often associated with implanted medical devices. We have sequenced the approximately 2.8-Mb genome of S. aureus COL, an early methicillin-resistant isolate, and the approximately 2.6-Mb genome of S. epidermidis RP62a, a methicillin-resistant biofilm isolate. Comparative analysis of these and other staphylococcal genomes was used to explore the evolution of virulence and resistance between these two species. The S. aureus and S. epidermidis genomes are syntenic throughout their lengths and share a core set of 1,681 open reading frames. Genome islands in nonsyntenic regions are the primary source of variations in pathogenicity and resistance. Gene transfer between staphylococci and low-GC-content gram-positive bacteria appears to have shaped their virulence and resistance profiles. Integrated plasmids in S. epidermidis carry genes encoding resistance to cadmium and species-specific LPXTG surface proteins. A novel genome island encodes multiple phenol-soluble modulins, a potential S. epidermidis virulence factor. S. epidermidis contains the cap operon, encoding the polyglutamate capsule, a major virulence factor in Bacillus anthracis. Additional phenotypic differences are likely the result of single nucleotide polymorphisms, which are most numerous in cell envelope proteins. Overall differences in pathogenicity can be attributed to genome islands in S. aureus which encode enterotoxins, exotoxins, leukocidins, and leukotoxins not found in S. epidermidis.

The 1.2-megabase genome sequence of Mimivirus

We recently reported the discovery and preliminary characterization of Mimivirus, the largest known virus, with a 400-nanometer particle size comparable to mycoplasma. Mimivirus is a double-stranded DNA virus growing in amoebae. We now present its 1,181,404-base pair genome sequence, consisting of 1262 putative open reading frames, 10% of which exhibit a similarity to proteins of known functions. In addition to exceptional genome size, Mimivirus exhibits many features that distinguish it from other nucleocytoplasmic large DNA viruses. The most unexpected is the presence of numerous genes encoding central protein-translation components, including four amino-acyl transfer RNA synthetases, peptide release factor 1, translation elongation factor EF-TU, and translation initiation factor 1. The genome also exhibits six tRNAs. Other notable features include the presence of both type I and type II topoisomerases, components of all DNA repair pathways, many polysaccharide synthesis enzymes, and one intein-containing gene. The size and complexity of the Mimivirus genome challenge the established frontier between viruses and parasitic cellular organisms. This new sequence data might help shed a new light on the origin of DNA viruses and their role in the early evolution of eukaryotes.

Bacterial thymidylate synthase with intein, group II Intron, and distinctive ThyX motifs

The ThyX class of thymidylate synthases was previously characterized by a common ThyX motif, RHRX7S. We report bacterial ThyX sequences having distinctive ThyX motifs, suggesting a more general ThyX motif, R/THRX7-8S. One ThyX sequence has an intein in its ThyX motif that was shown to do protein splicing and a group II intron in its gene, suggesting a hot spot for these self-splicing mobile elements.

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.

DNA Binding and Degradation by the HNH Protein ColE7

The bacterial toxin ColE7 bears an HNH motif which has been identified in hundreds of prokaryotic and eukaryotic endonucleases, involved in DNA homing, restriction, repair, or chromosome degradation. The crystal structure of the nuclease domain of ColE7 in complex with a duplex DNA has been determined at 2.5 A resolution. The HNH motif is bound at the minor groove primarily to DNA phosphate groups at and beyond the 3' side of the scissile phosphate, with little interaction with ribose groups and bases. This result provides a structural basis for sugar- and sequence-independent DNA recognition and the inhibition mechanism by inhibitor Im7, which blocks the substrate binding site but not the active site. Structural comparison shows that two families of endonucleases bind and bend DNA in a similar way to that of the HNH ColE7, indicating that endonucleases containing a "betabetaalpha-metal" fold of active site possess a universal mode for protein-DNA interactions.

Four inteins and three group II introns encoded in a bacterial ribonucleotide reductase gene

A bacterial ribonucleotide reductase gene was found to encode four inteins and three group II introns in the oceanic N2-fixing cyanobacterium Trichodesmium erythraeum. The 13,650-bp ribonucleotide reductase gene is divided into eight extein- or exon-coding sequences that together encode a 768-amino acid mature ribonucleotide reductase protein, with 83% of the gene sequence encoding introns and inteins. The four inteins are encoded on the second half of the gene, and each has conserved sequence motifs for a protein-splicing domain and an endonuclease domain. These four inteins, together with known inteins, define five intein insertion sites in ribonucleotide reductase homologues. Two of the insertion sites are 10 amino acids apart and next to key catalytic residues of the enzyme. Protein-splicing activities of all four inteins were demonstrated in Escherichia coli. The four inteins coexist with three group II introns encoded on the first half of the same gene, which suggests a breakdown of the presumed barrier against intron insertion in this bacterial conserved protein-coding gene.

Bacteriophage T4 genome

Phage T4 has provided countless contributions to the paradigms of genetics and biochemistry. Its complete genome sequence of 168,903 bp encodes about 300 gene products. T4 biology and its genomic sequence provide the best-understood model for modern functional genomics and proteomics. Variations on gene expression, including overlapping genes, internal translation initiation, spliced genes, translational bypassing, and RNA processing, alert us to the caveats of purely computational methods. The T4 transcriptional pattern reflects its dependence on the host RNA polymerase and the use of phage-encoded proteins that sequentially modify RNA polymerase; transcriptional activator proteins, a phage sigma factor, anti-sigma, and sigma decoy proteins also act to specify early, middle, and late promoter recognition. Posttranscriptional controls by T4 provide excellent systems for the study of RNA-dependent processes, particularly at the structural level. The redundancy of DNA replication and recombination systems of T4 reveals how phage and other genomes are stably replicated and repaired in different environments, providing insight into genome evolution and adaptations to new hosts and growth environments. Moreover, genomic sequence analysis has provided new insights into tail fiber variation, lysis, gene duplications, and membrane localization of proteins, while high-resolution structural determination of the "cell-puncturing device," combined with the three-dimensional image reconstruction of the baseplate, has revealed the mechanism of penetration during infection. Despite these advances, nearly 130 potential T4 genes remain uncharacterized. Current phage-sequencing initiatives are now revealing the similarities and differences among members of the T4 family, including those that infect bacteria other than Escherichia coli. T4 functional genomics will aid in the interpretation of these newly sequenced T4-related genomes and in broadening our understanding of the complex evolution and ecology of phages-the most abundant and among the most ancient biological entities on Earth.

Metal ions and phosphate binding in the H-N-H motif: crystal structures of the nuclease domain of ColE7/Im7 in complex with a phosphate ion and different divalent metal ions

H-N-H is a motif found in the nuclease domain of a subfamily of bacteria toxins, including colicin E7, that are capable of cleaving DNA nonspecifically. This H-N-H motif has also been identified in a subfamily of homing endonucleases, which cleave DNA site specifically. To better understand the role of metal ions in the H-N-H motif during DNA hydrolysis, we crystallized the nuclease domain of colicin E7 (nuclease-ColE7) in complex with its inhibitor Im7 in two different crystal forms, and we resolved the structures of EDTA-treated, Zn(2+)-bound and Mn(2+)-bound complexes in the presence of phosphate ions at resolutions of 2.6 A to 2.0 A. This study offers the first determination of the structure of a metal-free and substrate-free enzyme in the H-N-H family. The H-N-H motif contains two antiparallel beta-strands linked to a C-terminal alpha-helix, with a divalent metal ion located in the center. Here we show that the metal-binding sites in the center of the H-N-H motif, for the EDTA-treated and Mg(2+)-soaked complex crystals, were occupied by water molecules, indicating that an alkaline earth metal ion does not reside in the same position as a transition metal ion in the H-N-H motif. However, a Zn(2+) or Mn(2+) ions were observed in the center of the H-N-H motif in cases of Zn(2+) or Mn(2+)-soaked crystals, as confirmed in anomalous difference maps. A phosphate ion was found to bridge between the divalent transition metal ion and His545. Based on these structures and structural comparisons with other nucleases, we suggest a functional role for the divalent transition metal ion in the H-N-H motif in stabilizing the phosphoanion in the transition state during hydrolysis.

The crystal structure of the nuclease domain of colicin E7 suggests a mechanism for binding to double-stranded DNA by the H-N-H endonucleases

The bacterial toxin ColE7 contains an H-N-H endonuclease domain (nuclease ColE7) that digests cellular DNA or RNA non-specifically in target cells, leading to cell death. In the host cell, protein Im7 forms a complex with ColE7 to inhibit its nuclease activity. Here, we present the crystal structure of the unbound nuclease ColE7 at a resolution of 2.1A. Structural comparison between the unbound and bound nuclease ColE7 in complex with Im7, suggests that Im7 is not an allosteric inhibitor that induces backbone conformational changes in nuclease ColE7, but rather one that inhibits by blocking the substrate-binding site. There were two nuclease ColE7 molecules in the P1 unit cell in crystals and they appeared as a dimer related to each other by a non-crystallographic dyad symmetry. Gel-filtration and cross-linking experiments confirmed that nuclease ColE7 indeed formed dimers in solution and that the dimeric conformation was more favored in the presence of double-stranded DNA. Structural comparison of nuclease ColE7 with the His-Cys box homing endonuclease I-PpoI further demonstrated that H-N-H motifs in dimeric nuclease ColE7 were oriented in a manner very similar to that of the betabetaalpha-fold of the active sites found in dimeric I-PpoI. A mechanism for the binding of double-stranded DNA by dimeric H-N-H nuclease ColE7 is suggested.

Sequence analysis of the lactococcal bacteriophage bIL170: insights into structural proteins and HNH endonucleases in dairy phages

The complete 31754 bp genome of bIL170, a virulent bacteriophage of Lactococcus lactis belonging to the 936 group, was analysed. Sixty-four ORFs were predicted and the function of 16 of them was assigned by significant homology to proteins in databases. Three putative homing endonucleases of the HNH family were found in the early region. An HNH endonuclease with zinc-binding motif was identified in the late cluster, potentially being part of the same functional module as terminase. Three putative structural proteins were analysed in detail and show interesting features among dairy phages. Notably, gpl12 (putative fibre) and gpl20 (putative baseplate protein) of bIL170 are related by at least one of their domains to a number of multi-domain proteins encoded by lactococcal or streptococcal phages. A 110- to 150-aa-long hypervariable domain flanked by two conserved motifs of about 20 aa was identified. The analysis presented here supports the participation of some of these proteins in host-range determination and suggests that specific adsorption to the host may involve a complex multi-component system. Divergences in the genome of phages of the 936 group, that may have important biological properties, were noted. Insertions/deletions of units of one or two ORFs were the main source of divergence in the early clusters of the two entirely sequenced phages, bIL170 and sk1. An exchange of fragments probably affected the regions containing the putative origin of replication. It led to the absence in bIL170 of the direct repeats recognized in sk1 and to the presence of different ORFs in the ori region. Shuffling of protein domains affected the endolysin (putative cell-wall binding part), as well as gpl12 and gpl20.

The genome of bacteriophage phiKZ of Pseudomonas aeruginosa

Bacteriophage phiKZ is a giant virus that efficiently infects Pseudomonas aeruginosa strains pathogenic to human and, therefore, it is attractive for phage therapy. We present here the complete phiKZ genome sequence and a preliminary analysis of its genome structure. The 280,334 bp genome is a linear, circularly permutated and terminally redundant, A+T-rich double-stranded DNA molecule. The phiKZ DNA has no detectable sequence homology to other viruses and microorganisms, and it does not contain NotI, PstI, SacI, SmaI, XhoI, and XmaIII endonuclease restriction sites. The genome has 306 open reading frames (ORFs) varying in size from 50 to 2237 amino acid residues. According to the orientation of transcription, ORFs are apparently organized into clusters and most have a clockwise direction. The phiKZ genome also encodes six tRNAs specific for Met (AUG), Asn (AAC), Asp (GAC), Leu (TTA), Thr (ACA), and Pro (CCA). A putative promoter sequence containing a TATATTAC block was identified. Most potential stem-loop transcription terminators contain the tetranucleotide UUCG loops. Some genes may be assigned as phage-encoded RNA polymerase subunits. Only 59 phiKZ gene products exhibit similarity to proteins of known function from a diversity of organisms. Most of these conserved gene products, such as dihydrofolate reductase, ribonucleoside diphosphate reductase, thymidylate synthase, thymidylate kinase, and deoxycytidine triphosphate deaminase are involved in nucleotide metabolism. However, no virus-encoded DNA polymerase, DNA replication-associated proteins, or single-stranded DNA-binding protein were found based on amino acid homology, and they may therefore be strongly divergent from known homologous proteins. Fifteen phiKZ gene products show homology to proteins of pathogenic organisms, including Mycobacterium tuberculosis, Haemophilus influenzae, Listeria sp., Rickettsia prowazakeri, and Vibrio cholerae that must be considered before using this phage as a therapeutic agent. The phiKZ coat contains at least 40 polypeptides, and several proteins are cleaved during virus assembly in a way similar to phage T4. Eleven phiKZ-encoded polypeptides are related to proteins of other bacteriphages that infect a variety of hosts. Among these are four gene products that contain a putative intron-encoded endonuclease harboring the H-N-H motif common to many double-stranded DNA phages. These observations provide evidence that phages infecting diverse hosts have had access to a common genetic pool. However, limited homology on the DNA and protein levels indicates that bacteriophage phiKZ represents an evolutionary distinctive branch of the Myoviridae family.

Analysis of transcription of the Staphylococcus aureus aerobic class Ib and anaerobic class III ribonucleotide reductase genes in response to oxygen

Staphylococcus aureus is a gram-positive facultative aerobe that can grow in the absence of oxygen by fermentation or by using an alternative electron acceptor. To investigate the mechanism by which S. aureus is able to adapt to changes in oxygen concentration, we analyzed the transcriptional regulation of genes that encode the aerobic class Ib and anaerobic class III ribonucleotide reductase (RNR) systems that are responsible for the synthesis of deoxyribonucleotides needed for DNA synthesis. The S. aureus class Ib RNR nrdIEF and class III RNR nrdDG genes and their regulatory regions were cloned and sequenced. Inactivation of the nrdDG genes showed that the class III RNR is essential for anaerobic growth. Inhibition of aerobic growth by hydroxyurea showed that the class Ib RNR is an oxygen-dependent enzyme. Northern blot analysis and primer extension analysis demonstrated that transcription of class III nrdDG genes is regulated by oxygen concentration and was at least 10-fold higher under anaerobic than under aerobic conditions. In contrast, no significant effect of oxygen concentration was found on the transcription of class Ib nrdIEF genes. Disruption or deletion of S. aureus nrdDG genes caused up to a fivefold increase in nrdDG and nrdIEF transcription under anaerobic conditions but not under aerobic conditions. Similarly, hydroxyurea, an inhibitor of the class I RNRs, resulted in increased transcription of class Ib and class III RNR genes under aerobic conditions. These findings establish that transcription of class Ib and class III RNR genes is upregulated under conditions that cause the depletion of deoxyribonucleotide. Promoter analysis of class Ib and class III RNR operons identified several inverted-repeat elements that may account for the transcriptional response of the nrdIEF and nrdDG genes to oxygen.

Ribonucleotide reductase genes of Bacillus prophages: a refuge to introns and intein coding sequences

The ribonucleotide reductase gene tandem bnrdE/bnrdF in SPbeta-related prophages of different Bacillus spp. isolates presents different configurations of intervening sequences, comprising one to three of six non-homologous splicing elements. Insertion sites of group I introns and intein DNA are clustered in three relatively short segments encoding functionally important domains of the ribonucleotide reductase. Comparison of the bnrdE homologs reveals mutual exclusion of a group I intron and an intein coding sequence flanking the codon that specifies a conserved cysteine. In vivo splicing was demonstrated for all introns. However, for two of them a part of the mRNA precursor molecules remains unspliced. Intergenic bnrdE-bnrdF regions are unexpectedly long, comprising between 238 and 541 nt. The longest encodes a putative polypeptide related to HNH homing endonucleases.

Protein-splicing intein: Genetic mobility, origin, and evolution

Intein is the protein equivalent of intron and has been discovered in increasing numbers of organisms and host proteins. A self-splicing intein catalyzes its own removal from the host protein through a posttranslational process of protein splicing. A mobile intein displays a site-specific endonuclease activity that confers genetic mobility to the intein through intein homing. Recent findings of intein structure and the mechanism of protein splicing illuminated how inteins work and yielded clues regarding intein's origin, spread, and evolution. Inteins can evolve into new structures and new functions, such as split inteins that do trans-splicing. The structural basis of intein function needs to be identified for a full understanding of the origin and evolution of this marvelous genetic element.

A chimeric ribozyme in clostridium difficile combines features of group I introns and insertion elements

CdlSt1, a DNA insertion of 1975 bp, was identified within tcdA-C34, the enterotoxin gene of the Clostridium difficile isolate C34. Located in the catalytic domain A1-C34, Cd/St1 combines features of two genetic elements. Within the first 434 nt structures characteristic for group I introns were found; encoding the two transposase-like proteins tlpA and tlpB nucleotides 435-1975 represent the remainder of a IS605-like insertion element. We show that the entire CdlSt1 is accurately spliced from tcdA-C34 primary transcripts and that purified TcdA-C34 toxin is of regular size and catalytic activity. A search for CdlSt1-related sequences demonstrates that the element is widespread in toxinogenic and non-toxinogenic C. difficile strains, indicating the mobility of CdlSt1. In strain C34, we characterize 10 CdlSt1 variants; all are highly homologous to CdlSt1 (> 93% identity), integrated in bacterial open reading frames (ORFs), show the typical composite structure of CdlSt1 and are precisely spliced from their primary transcripts. CdlSt1-like chimeric ribozymes appear to combine the invasiveness of an insertion element with the splicing ability of a group I intron, rendering transposition harmless for the interrupted gene.

Widespread distribution of a group I intron and its three deletion derivatives in the lysin gene of Streptococcus thermophilus bacteriophages

Of 62 Streptococcus thermophilus bacteriophages isolated from various ecological settings, half contain a lysin gene interrupted by a group IA2 intron. Phage mRNA splicing was demonstrated. Five phages possess a variant form of the intron resulting from three distinct deletion events located in the intron-harbored open reading frame (orf 253). The predicted orf 253 gene sequence showed a significantly lower GC content than the surrounding intron and lysin gene sequences, and the predicted protein shared a motif with endonucleases found in phages from both gram-positive and gram-negative bacteria. A comparison of the phage lysin genes revealed a clear division between intron-containing and intron-free alleles, leading to the establishment of a 14-bp consensus sequence associated with intron possession. The conserved intron was not found elsewhere in the phage or S. thermophilus bacterial genomes. Folding of the intron RNA revealed secondary structure elements shared with other phage introns: first, a 38-bp insertion between regions P3 and P4 that can be folded into two stem-loop structures (shared with introns from Bacillus phage SPO1 and relatives); second, a conserved P7.2 region (shared with all phage introns); third, the location of the stop codon from orf 253 in the P8 stem (shared with coliphage T4 and Bacillus phage SPO1 introns); fourth, orf 253, which has sequence similarity with the H-N-H motif of putative endonuclease genes found in introns from Lactococcus, Lactobacillus, and Bacillus phages.

Evolution and horizontal transfer of dUTPase-encoding genes in viruses and their hosts

dUTPase is a ubiquitous and essential enzyme responsible for regulating cellular levels of dUTP. The dut gene exists as single, tandemly duplicated, and tandemly triplicated copies. Crystallized single-copy dUTPases have been shown to assemble as homotrimers. dUTPase is encoded as an auxiliary gene in a number of virus genomes. The origin of viral dut genes has remained unresolved since their initial discovery. A comprehensive analysis of dUTPase amino acid sequence relationships was performed to explore the evolutionary dynamics of dut in viruses and their hosts. Our data set, comprised of 24 host and 51 viral sequences, includes representative sequences from available eukaryotes, archaea, eubacteria cells, and viruses, including herpesviruses. These amino acid sequences were aligned by using a hidden Markov model approach developed to align divergent data. Known secondary structures from single-copy crystals were mapped onto the aligned duplicate and triplicate sequences. We show how duplicated dUTPases might fold into a monomer, and we hypothesize that triplicated dUTPases also assemble as monomers. Phylogenetic analysis revealed at least five viral dUTPase sequence lineages in well-supported monophyletic clusters with eukaryotic, eubacterial, and archaeal hosts. We have identified all five as strong examples of horizontal transfer as well as additional potential transfer of dut genes among eubacteria, between eubacteria and viruses, and between retroviruses. The evidence for horizontal transfers is particularly interesting since eukaryotic dut genes have introns, while DNA virus dut genes do not. This implies that an intermediary retroid agent facilitated the horizontal transfer process between host mRNA and DNA viruses.