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

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.

Birds Kept in the German Zoo "Tierpark Berlin" Are a Common Source for Polyvalent Yersinia pseudotuberculosis Phages

Yersinia pseudotuberculosis is an important animal pathogen, particularly for birds, rodents, and monkeys, which is also able to infect humans. Indeed, an increasing number of reports have been published on zoo animals that were killed by this species. One option to treat diseased animals is the application of strictly lytic (virulent) phages. However, thus far relatively few phages infecting Y. pseudotuberculosis have been isolated and characterized. To determine the prevalence of Y. pseudotuberculosis phages in zoo animals, fecal samples of birds and some primates, maras, and peccaries kept in the Tierpark Berlin were analyzed. Seventeen out of 74 samples taken in 2013 and 2017 contained virulent phages. The isolated phages were analyzed in detail and could be allocated to three groups. The first group is composed of 10 T4-like phages (PYps2T taxon group: Myoviridae; Tevenvirinae; Tequatrovirus), the second group (PYps23T taxon group: Chaseviridae; Carltongylesvirus; Escherichia virus ST32) consists of five phages encoding a podovirus-like RNA polymerase that is related to an uncommon genus of myoviruses (e.g., Escherichia coli phage phiEcoM-GJ1), while the third group is comprised of two podoviruses (PYps50T taxon group: Autographiviridae; Studiervirinae; Berlinvirus) which are closely related to T7. The host range of the isolated phages differed significantly. Between 5.5 and 86.7% of 128 Y. pseudotuberculosis strains belonging to 20 serotypes were lysed by each phage. All phages were additionally able to lyse Y. enterocolitica B4/O:3 strains, when incubated at 37°C. Some phages also infected Y. pestis strains and even strains belonging to other genera of Enterobacteriaceae. A cocktail containing two of these phages would be able to lyse almost 93% of the tested Y. pseudotuberculosis strains. The study indicates that Y. pseudotuberculosis phages exhibiting a broad-host range can be isolated quite easily from zoo animals, particularly birds.

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.

Phage T4 endonuclease SegD that is similar to group I intron endonucleases does not initiate homing of its own gene

Homing endonucleases are a group of site-specific endonucleases that initiate homing, a nonreciprocal transfer of its own gene into a new allele lacking this gene. This work describes a novel phage T4 endonuclease, SegD, which is homologous to the GIY-YIG family of homing endonucleases. Like other T4 homing endonucleases SegD recognizes an extended, 16bp long, site, cleaves it asymmetrically to form 3'-protruding ends and digests both unmodified DNA and modified T-even phage DNA with similar efficiencies. Surprisingly, we revealed that SegD cleavage site was identical in the genomes of segD^- and segD⁺ phages. We found that segD gene was expressed during the T4 developmental cycle. Nevertheless, endonuclease SegD was not able to initiate homing of its own gene as well as genetic recombination between phages in its site inserted into the rII locus.

MegaTevs: single-chain dual nucleases for efficient gene disruption

Targeting gene disruptions in complex genomes relies on imprecise repair by the non-homologous end-joining DNA pathway, creating mutagenic insertions or deletions (indels) at the break point. DNA end-processing enzymes are often co-expressed with genome-editing nucleases to enhance the frequency of indels, as the compatible cohesive ends generated by the nucleases can be precisely repaired, leading to a cycle of cleavage and non-mutagenic repair. Here, we present an alternative strategy to bias repair toward gene disruption by fusing two different nuclease active sites from I-TevI (a GIY-YIG enzyme) and I-OnuI E2 (an engineered meganuclease) into a single polypeptide chain. In vitro, the MegaTev enzyme generates two double-strand breaks to excise an intervening 30-bp fragment. In HEK 293 cells, we observe a high frequency of gene disruption without co-expression of DNA end-processing enzymes. Deep sequencing of disrupted target sites revealed minimal processing, consistent with the MegaTev sequestering the double-strand breaks from the DNA repair machinery. Off-target profiling revealed no detectable cleavage at sites where the I-TevI CNNNG cleavage motif is not appropriately spaced from the I-OnuI binding site. The MegaTev enzyme represents a small, programmable nuclease platform for extremely specific genome-engineering applications.

The monomeric GIY-YIG homing endonuclease I-BmoI uses a molecular anchor and a flexible tether to sequentially nick DNA

The GIY-YIG nuclease domain is found within protein scaffolds that participate in diverse cellular pathways and contains a single active site that hydrolyzes DNA by a one-metal ion mechanism. GIY-YIG homing endonucleases (GIY-HEs) are two-domain proteins with N-terminal GIY-YIG nuclease domains connected to C-terminal DNA-binding and they are thought to function as monomers. Using I-BmoI as a model GIY-HE, we test mechanisms by which the single active site is used to generate a double-strand break. We show that I-BmoI is partially disordered in the absence of substrate, and that the GIY-YIG domain alone has weak affinity for DNA. Significantly, we show that I-BmoI functions as a monomer at all steps of the reaction pathway and does not transiently dimerize or use sequential transesterification reactions to cleave substrate. Our results are consistent with the I-BmoI DNA-binding domain acting as a molecular anchor to tether the GIY-YIG domain to substrate, permitting rotation of the GIY-YIG domain to sequentially nick each DNA strand. These data highlight the mechanistic differences between monomeric GIY-HEs and dimeric or tetrameric GIY-YIG restriction enzymes, and they have implications for the use of the GIY-YIG domain in genome-editing applications.

Frequent gain and loss of introns in fungal cytochrome b genes

In this study, all available cytochrome b (Cyt b) genes from the GOBASE database were compiled and the evolutionary dynamics of the Cyt b gene introns was assessed. Cyt b gene introns were frequently present in the fungal kingdom and some lower plants, but generally absent or rare in Chromista, Protozoa, and Animalia. Fungal Cyt b introns were found at 35 positions in Cyt b genes and the number of introns varied at individual positions from a single representative to 32 different introns at position 131, showing a wide and patchy distribution. Many homologous introns were present at the same position in distantly related species but absent in closely related species, suggesting that introns of the Cyt b genes were frequently lost. On the other hand, highly similar intron sequences were observed in some distantly related species rather than in closely related species, suggesting that these introns were gained independently, likely through lateral transfers. The intron loss-and-gain events could be mediated by transpositions that might have occurred between nuclear and mitochondria. Southern hybridization analysis confirmed that some introns contained repetitive sequences and might be transposable elements. An intron gain in Botryotinia fuckeliana prevented the development of QoI fungicide resistance, suggesting that intron loss-and-gain events were not necessarily beneficial to their host organisms.

Separate origins of group I introns in two mitochondrial genes of the katablepharid Leucocryptos marina

Mitochondria are descendants of the endosymbiotic α-proteobacterium most likely engulfed by the ancestral eukaryotic cells, and the proto-mitochondrial genome should have been severely streamlined in terms of both genome size and gene repertoire. In addition, mitochondrial (mt) sequence data indicated that frequent intron gain/loss events contributed to shaping the modern mt genome organizations, resulting in the homologous introns being shared between two distantly related mt genomes. Unfortunately, the bulk of mt sequence data currently available are of phylogenetically restricted lineages, i.e., metazoans, fungi, and land plants, and are insufficient to elucidate the entire picture of intron evolution in mt genomes. In this work, we sequenced a 12 kbp-fragment of the mt genome of the katablepharid Leucocryptos marina. Among nine protein-coding genes included in the mt genome fragment, the genes encoding cytochrome b and cytochrome c oxidase subunit I (cob and cox1) were interrupted by group I introns. We further identified that the cob and cox1 introns host open reading frames for homing endonucleases (HEs) belonging to distantly related superfamilies. Phylogenetic analyses recovered an affinity between the HE in the Leucocryptos cob intron and two green algal HEs, and that between the HE in the Leucocryptos cox1 intron and a fungal HE, suggesting that the Leucocryptos cob and cox1 introns possess distinct evolutionary origins. Although the current intron (and intronic HE) data are insufficient to infer how the homologous introns were distributed to distantly related mt genomes, the results presented here successfully expanded the evolutionary dynamism of group I introns in mt genomes.

The genome and proteome of a Campylobacter coli bacteriophage vB_CcoM-IBB_35 reveal unusual features

Background

Campylobacter is the leading cause of foodborne diseases worldwide. Bacteriophages (phages) are naturally occurring predators of bacteria, ubiquitous in the environment, with high host specificity and thus considered an appealing option to control bacterial pathogens. Nevertheless for an effective use of phages as antimicrobial agents, it is important to understand phage biology which renders crucial the analysis of phage genomes and proteomes. The lack of sequence data from Campylobacter phages adds further importance to these studies.

Methods

vB_CcoM-IBB_35 is a broad lytic spectrum Myoviridae Campylobacter phage with high potential for therapeutic use. The genome of this phage was obtained by pyrosequencing and the sequence data was further analyzed. The proteomic analysis was performed by SDS-PAGE and Mass spectrometry.

Results and conclusions

The DNA sequence data of vB_CcoM-IBB_35 consists of five contigs for a total of 172,065 bp with an average GC content of 27%. Attempts to close the gaps between contigs were unsuccessful since the DNA preparations appear to contain substances that inhibited Taq and ϕ29 polymerases. From the 210 identified ORFs, around 60% represent proteins that were not functionally assigned. Homology exists with members of the Teequatrovirinae namely for T4 proteins involved in morphogenesis, nucleotide metabolism, transcription, DNA replication and recombination. Tandem mass spectrometric analysis revealed 38 structural proteins as part of the mature phage particle.

Conclusions

Genes encoding proteins involved in the carbohydrate metabolism along with several incidences of gene duplications, split genes with inteins and introns have been rarely found in other phage genomes yet are found in this phage. We identified the genes encoding for tail fibres and for the lytic cassette, this later, expressing enzymes for bacterial capsular polysaccharides (CPS) degradation, which has not been reported before for Campylobacter phages.

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.

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

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

Genomes of the T4-related bacteriophages as windows on microbial genome evolution

The T4-related bacteriophages are a group of bacterial viruses that share morphological similarities and genetic homologies with the well-studied Escherichia coli phage T4, but that diverge from T4 and each other by a number of genetically determined characteristics including the bacterial hosts they infect, the sizes of their linear double-stranded (ds) DNA genomes and the predicted compositions of their proteomes. The genomes of about 40 of these phages have been sequenced and annotated over the last several years and are compared here in the context of the factors that have determined their diversity and the diversity of other microbial genomes in evolution. The genomes of the T4 relatives analyzed so far range in size between ~160,000 and ~250,000 base pairs (bp) and are mosaics of one another, consisting of clusters of homology between them that are interspersed with segments that vary considerably in genetic composition between the different phage lineages. Based on the known biological and biochemical properties of phage T4 and the proteins encoded by the T4 genome, the T4 relatives reviewed here are predicted to share a genetic core, or "Core Genome" that determines the structural design of their dsDNA chromosomes, their distinctive morphology and the process of their assembly into infectious agents (phage morphogenesis). The Core Genome appears to be the most ancient genetic component of this phage group and constitutes a mere 12-15% of the total protein encoding potential of the typical T4-related phage genome. The high degree of genetic heterogeneity that exists outside of this shared core suggests that horizontal DNA transfer involving many genetic sources has played a major role in diversification of the T4-related phages and their spread to a wide spectrum of bacterial species domains in evolution. We discuss some of the factors and pathways that might have shaped the evolution of these phages and point out several parallels between their diversity and the diversity generally observed within all groups of interrelated dsDNA microbial genomes in nature.

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.

Conservation of intron and intein insertion sites: implications for life histories of parasitic genetic elements

BACKGROUND

Inteins and introns are genetic elements that are removed from proteins and RNA after translation or transcription, respectively. Previous studies have suggested that these genetic elements are found in conserved parts of the host protein. To our knowledge this type of analysis has not been done for group II introns residing within a gene. Here we provide quantitative statistical support from an analyses of proteins that host inteins, group I introns, group II introns and spliceosomal introns across all three domains of life.

RESULTS

To determine whether or not inteins, group I, group II, and spliceosomal introns are found preferentially in conserved regions of their respective host protein, conservation profiles were generated and intein and intron positions were mapped to the profiles. Fisher's combined probability test was used to determine the significance of the distribution of insertion sites across the conservation profile for each protein. For a subset of studied proteins, the conservation profile and insertion positions were mapped to protein structures to determine if the insertion sites correlate to regions of functional activity. All inteins and most group I introns were found to be preferentially located within conserved regions; in contrast, a bacterial intein-like protein, group II and spliceosomal introns did not show a preference for conserved sites.

CONCLUSIONS

These findings demonstrate that inteins and group I introns are found preferentially in conserved regions of their respective host proteins. Homing endonucleases are often located within inteins and group I introns and these may facilitate mobility to conserved regions. Insertion at these conserved positions decreases the chance of elimination, and slows deletion of the elements, since removal of the elements has to be precise as not to disrupt the function of the protein. Furthermore, functional constrains on the targeted site make it more difficult for hosts to evolve immunity to the homing endonuclease. Therefore, these elements will better survive and propagate as molecular parasites in conserved sites. In contrast, spliceosomal introns and group II introns do not show significant preference for conserved sites and appear to have adopted a different strategy to evade loss.

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.

Survey of group I and group II introns in 29 sequenced genomes of the Bacillus cereus group: insights into their spread and evolution

Group I and group II introns are different catalytic self-splicing and mobile RNA elements that contribute to genome dynamics. In this study, we have analyzed their distribution and evolution in 29 sequenced genomes from the Bacillus cereus group of bacteria. Introns were of different structural classes and evolutionary origins, and a large number of nearly identical elements are shared between multiple strains of different sources, suggesting recent lateral transfers and/or that introns are under a strong selection pressure. Altogether, 73 group I introns were identified, inserted in essential genes from the chromosome or newly described prophages, including the first elements found within phages in bacterial plasmids. Notably, bacteriophages are an important source for spreading group I introns between strains. Furthermore, 77 group II introns were found within a diverse set of chromosomal and plasmidic genes. Unusual findings include elements located within conserved DNA metabolism and repair genes and one intron inserted within a novel retroelement. Group II introns are mainly disseminated via plasmids and can subsequently invade the host genome, in particular by coupling mobility with host cell replication. This study reveals a very high diversity and variability of mobile introns in B. cereus group strains.

Phage T4 SegB protein is a homing endonuclease required for the preferred inheritance of T4 tRNA gene region occurring in co-infection with a related phage

Homing endonucleases initiate nonreciprocal transfer of DNA segments containing their own genes and the flanking sequences by cleaving the recipient DNA. Bacteriophage T4 segB gene, which is located in a cluster of tRNA genes, encodes a protein of unknown function, homologous to homing endonucleases of the GIY-YIG family. We demonstrate that SegB protein is a site-specific endonuclease, which produces mostly 3′ 2-nt protruding ends at its DNA cleavage site. Analysis of SegB cleavage sites suggests that SegB recognizes a 27-bp sequence. It contains 11-bp conserved sequence, which corresponds to a conserved motif of tRNA TψC stem-loop, whereas the remainder of the recognition site is rather degenerate. T4-related phages T2L, RB1 and RB3 contain tRNA gene regions that are homologous to that of phage T4 but lack segB gene and several tRNA genes. In co-infections of phages T4 and T2L, segB gene is inherited with nearly 100% of efficiency. The preferred inheritance depends absolutely on the segB gene integrity and is accompanied by the loss of the T2L tRNA gene region markers. We suggest that SegB is a homing endonuclease that functions to ensure spreading of its own gene and the surrounding tRNA genes among T4-related phages.

Comparison of genomes of three Xanthomonas oryzae bacteriophages

BACKGROUND

Xp10 and OP1 are phages of Xanthomonas oryzae pv. oryzae (Xoo), the causative agent of bacterial leaf blight in rice plants, which were isolated in 1967 in Taiwan and in 1954 in Japan, respectively. We recently isolated the Xoo phage Xop411.

RESULTS

The linear Xop411 genome (44,520 bp, 58 ORFs) sequenced here is 147 bp longer than that of Xp10 (60 ORFs) and 735 bp longer than that of OP1 (59 ORFs). The G+C contents of OP1 (51%) and Xop411 and Xp10 (52% each) are less than that of the host (65%). The 9-bp 3'-overhangs (5'-GGACAGTCT-3') in Xop411 and Xp10 are absent from OP1. More of the deduced Xop411 proteins share higher degrees of identity with Xp10 than with OP1 proteins, while the right end of the genomes of Xp10 and OP1, containing all predicted promoters, share stronger homology. Xop411, Xp10, and OP1 contain 8, 7, and 6 freestanding HNH endonuclease genes, respectively. These genes can be classified into five groups depending on their possession of the HNH domain (HNN or HNH type) and/or AP2 domain in intact or truncated forms. While the HNN-AP2 type endonuclease genes dispersed in the genome, the HNH type endonuclease genes, each with a unique copy, were located within the same genome context. Mass spectrometry and N-terminal sequencing showed nine Xop411 coat proteins, among which three were identified, six were assigned as coat proteins (4) and conserved phage proteins (2) in Xp10. The major coat protein, in which only the N-terminal methionine is removed, appears to exist in oligomeric forms containing 2 to 6 subunits. The three phages exhibit different patterns of domain duplication in the N-terminus of the tail fiber, which are involved in determination of the host range. Many short repeated sequences are present in and around the duplicated domains.

CONCLUSION

Geographical separation may have confined lateral gene transfer among the Xoo phages. The HNN-AP2 type endonucleases were more likely to transfer their genes randomly in the genome and may degenerate after successful transmission. Some repeated sequences may be involved in duplication/loss of the domains in the tail fiber genes.

The molecular evolution and structural organization of group I introns at position 1389 in nuclear small subunit rDNA of myxomycetes

The number of nuclear group I introns from myxomycetes is rapidly increasing in GenBank as more rDNA sequences from these organisms are being sequenced. They represent an interesting and complex group of intervening sequences because several introns are mobile (or inferred to be mobile) and many contain large and unusual insertions in peripheral loops. Here we describe related group I introns at position 1389 in the small subunit rDNA of representatives from the myxomycete family Didymiaceae. Phylogenetic analyses support a common origin and mainly vertical inheritance of the intron. All S1389 introns from the Didymiaceae belong to the IC1 subclass of nuclear group I introns. The central catalytic core region of about 100 nt appears divergent in sequence composition even though the introns reside in closely related species. Furthermore, unlike the majority of group I introns from myxomycetes the S1389 introns do not self-splice as naked RNA in vitro under standard conditions, consistent with a dependence on host factors for folding or activity. Finally, the myxomycete S1389 introns are exclusively found within the family Didymiaceae, which suggests that this group I intron was acquired after the split between the families Didymiaceae and Physaraceae.

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.

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.

Evolution of Pleopsidium (Lichenized Ascomycota) S943 Group I Introns and the Phylogeography of an Intron-Encoded Putative Homing Endonuclease

The sporadic distribution of nuclear group I introns among different fungal lineages can be explained by vertical inheritance of the introns followed by successive losses, or horizontal transfers from one lineage to another through intron homing or reverse splicing. Homing is mediated by an intron-encoded homing endonuclease (HE) and recent studies suggest that the introns and their associated HE gene (HEG) follow a recurrent cyclical model of invasion, degeneration, loss, and reinvasion. The purpose of this study was to compare this model to the evolution of HEGs found in the group I intron at position S943 of the nuclear ribosomal DNA of the lichen-forming fungus Pleopsidium. Forty-eight S943 introns were found in the 64 Pleopsidium samples from a worldwide screen, 22 of which contained a full-length HEG that encodes a putative 256-amino acid HE, and 2 contained HE pseudogenes. The HEGs are divided into two closely related types (as are the introns that encode them) that differ by 22.6% in their nucleotide sequences. The evolution of the Pleopsidium intron-HEG element shows strong evidence for a cyclical model of evolution. The intron was likely acquired twice in the genus and then transmitted via two or three interspecific horizontal transfers. Close geographical proximity plays an important role in intron-HEG horizontal transfer because most of these mobile elements were found in Europe. Once acquired in a lineage, the intron-HEG element was also vertically transmitted, and occasionally degenerated or was lost.

Self-splicing of the bacteriophage T4 group I introns requires efficient translation of the pre-mRNA in vivo and correlates with the growth state of the infected bacterium

Bacteriophage T4 contains three self-splicing group I introns in genes in de novo deoxyribonucleotide biosynthesis (in td, coding for thymidylate synthase and in nrdB and nrdD, coding for ribonucleotide reductase). Their presence in these genes has fueled speculations that the introns are retained within the phage genome due to a possible regulatory role in the control of de novo deoxyribonucleotide synthesis. To study whether sequences in the upstream exon interfere with proper intron folding and splicing, we inhibited translation in T4-infected bacteria as well as in bacteria containing recombinant plasmids carrying the nrdB intron. Splicing was strongly reduced for all three T4 introns after the addition of chloramphenicol during phage infection, suggesting that the need for translating ribosomes is a general trait for unperturbed splicing. The splicing of the cloned nrdB intron was markedly reduced in the presence of chloramphenicol or when translation was hindered by stop codons inserted in the upstream exon. Several exon regions capable of forming putative interactions with nrdB intron sequences were identified, and the removal or mutation of these exon regions restored splicing efficiency in the absence of translation. Interestingly, splicing of the cloned nrdB intron was also reduced as cells entered stationary phase and splicing of all three introns was reduced upon the T4 infection of stationary-phase bacteria. Our results imply that conditions likely to be frequently encountered by natural phage populations may limit the self-splicing efficiency of group I introns. This is the first time that environmental effects on bacterial growth have been linked to the regulation of splicing of phage introns.

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.

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

T4 contains two groups of genes with similarity to homing endonucleases, the seg-genes (similarity to endonucleases encoded by group I introns) containing GIY-YIG motifs and the mob-genes (similarity to mobile endonucleases) containing H-N-H motifs. The four seg-genes characterized to date encode homing endonucleases with cleavage sites close to their respective gene loci while none of the mob-genes have been shown to cleave DNA. Of 18 phages screened, only T4 was found to have mobC while mobE genes were found in five additional phages. Interestingly, three phages encoded a seg-like gene (hereby called segH) with a GIY-YIG motif in place of mobC. An additional phage has an unrelated gene called hef (homing endonuclease-like function) in place of the mobE gene. The gene products of both novel genes displayed homing endonuclease activity with cleavage site specificity close to their respective genes. In contrast to intron encoded homing endonucleases, both SegH and Hef can cleave their own DNA as well as DNA from phages without the genes. Both segH and mobE (and most likely hef) can home between phages in mixed infections. We discuss why it might be a selective advantage for phage freestanding homing endonucleases to cleave both HEG-containing and HEG-less genomes.

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.

A self-splicing group I intron in DNA polymerase genes of T7-like bacteriophages

Group I introns are inserted into genes of a wide variety of bacteriophages of gram-positive bacteria. However, among the phages of enteric and other gram-negative proteobacteria, introns have been encountered only in phage T4 and several of its close relatives. Here we report the insertion of a self-splicing group I intron in the coding sequence of the DNA polymerase genes of PhiI and W31, phages that are closely related to T7. The introns belong to subgroup IA2 and both contain an open reading frame, inserted into structural element P6a, encoding a protein belonging to the HNH family of homing endonucleases. The introns splice efficiently in vivo and self-splice in vitro under mild conditions of ionic strength and temperature. We conclude that there is no barrier for maintenance of group I introns in phages of proteobacteria.

Homing endonuclease genes: the rise and fall and rise again of a selfish element

Homing endonuclease genes (HEGs) are selfish genetic elements that spread by first cleaving chromosomes that do not contain them and then getting copied across to the broken chromosome as a byproduct of the repair process. The success of this strategy will depend on the opportunities for homing--in other words, the frequency with which HEG(+) and HEG(-) chromosomes come into contact--which varies widely among host taxa. HEGs are also unusual in that the selection pressure for endonuclease function disappears if they become fixed in a population, which makes them susceptible to degeneration and imposes a need for regular horizontal transmission between species. HEGs will be selected to reduce the harm done to the host organism, and this is expected to influence the evolution of their sequence specificity and maturase functions. HEGs may also be domesticated by their hosts, and are currently being put to human uses.