An intron in the thymidylate synthase gene of Bacillus bacteriophage beta 22: evidence for independent evolution of a gene, its group I intron, and the intron open reading frame

An intron endonuclease facilitates interference competition between coinfecting viruses

Introns containing homing endonucleases are widespread in nature and have long been assumed to be selfish elements that provide no benefit to the host organism. These genetic elements are common in viruses, but whether they confer a selective advantage is unclear. In this work, we studied intron-encoded homing endonuclease gp210 in bacteriophage ΦPA3 and found that it contributes to viral competition by interfering with the replication of a coinfecting phage, ΦKZ. We show that gp210 targets a specific sequence in ΦKZ, which prevents the assembly of progeny viruses. This work demonstrates how a homing endonuclease can be deployed in interference competition among viruses and provide a relative fitness advantage. Given the ubiquity of homing endonucleases, this selective advantage likely has widespread evolutionary implications in diverse plasmid and viral competition as well as virus-host interactions.

Genome Study of a Novel Virulent Phage vB_SspS_KASIA and Mu-like Prophages of Shewanella sp. M16 Provides Insights into the Genetic Diversity of the Shewanella Virome

Shewanella is a ubiquitous bacterial genus of aquatic ecosystems, and its bacteriophages are also isolated from aquatic environments (oceans, lakes, ice, and wastewater). In this study, the isolation and characterization of a novel virulent Shewanella phage vB_SspS_KASIA and the identification of three prophages of its host, Shewanella sp. M16, including a mitomycin-inducible Mu-like siphovirus, vB_SspS_MuM16-1, became the starting point for comparative analyses of phages infecting Shewanella spp. and the determination of their position among the known bacterial viruses. A similarity networking analysis revealed the high diversity of Shewanella phages in general, with vB_SspS_KASIA clustering exclusively with Colwellia phage 9A, with which it forms a single viral cluster composed of two separate viral subclusters. Furthermore, vB_SspS_MuM16-1 presented itself as being significantly different from the phages deposited in public databases, expanding the diversity of the known Mu-like phages and giving potential molecular markers for the identification of Mu-like prophages in bacterial genomes. Moreover, the functional analysis performed for vB_SspS_KASIA suggested that, despite the KASIA host, the M16 strain grows better in a rich medium and at 30 °C the phage replication cycle seems to be optimal in restrictive culture conditions mimicking their natural environment, the Zloty Stok gold and arsenic mine.

Long-read viral metagenomics captures abundant and microdiverse viral populations and their niche-defining genomic islands

Marine viruses impact global biogeochemical cycles via their influence on host community structure and function, yet our understanding of viral ecology is constrained by limitations in host culturing and a lack of reference genomes and 'universal' gene markers to facilitate community surveys. Short-read viral metagenomic studies have provided clues to viral function and first estimates of global viral gene abundance and distribution, but their assemblies are confounded by populations with high levels of strain evenness and nucleotide diversity (microdiversity), limiting assembly of some of the most abundant viruses on Earth. Such features also challenge assembly across genomic islands containing niche-defining genes that drive ecological speciation. These populations and features may be successfully captured by single-virus genomics and fosmid-based approaches, at least in abundant taxa, but at considerable cost and technical expertise. Here we established a low-cost, low-input, high throughput alternative sequencing and informatics workflow to improve viral metagenomic assemblies using short-read and long-read technology. The 'VirION' (Viral, long-read metagenomics via MinION sequencing) approach was first validated using mock communities where it was found to be as relatively quantitative as short-read methods and provided significant improvements in recovery of viral genomes. We then then applied VirION to the first metagenome from a natural viral community from the Western English Channel. In comparison to a short-read only approach, VirION: (i) increased number and completeness of assembled viral genomes; (ii) captured abundant, highly microdiverse virus populations, and (iii) captured more and longer genomic islands. Together, these findings suggest that VirION provides a high throughput and cost-effective alternative to fosmid and single-virus genomic approaches to more comprehensively explore viral communities in nature.

Isolation and characterization of a N4-like lytic bacteriophage infecting Vibrio splendidus, a pathogen of fish and bivalves

A novel virulent bacteriophage, vB_VspP_pVa5, infecting a strain of Vibrio splendidus was isolated from a sea-cage aquaculture farm in Greece, and characterized using microbiological methods and genomic analysis. Bacteriophage vB_VspP_pVa5 is a N4-like podovirus with an icosahedral head measuring 85 nm in length and a short non-contractile tail. The phage had a narrow host range infecting only the bacterial host, a latent period of 30 min and a burst size of 24 virions per infected bacterium. Its genome size was 78,145 bp and genomic analysis identified 107 densely-packed genes, 40 of which could be annotated. In addition to the very large virion encapsulated DNA-dependent RNA polymerase which is the signature of the N4-like genus, an interesting feature of the novel phage is the presence of a self-splicing group I intron in the thymidylate synthase gene. A tRNA^Stop interrupted by a ~2.5kb open reading frame–containing area was also identified. The absence of genes related to lysogeny along with the high efficacy observed during in vitro cell lysis trials, indicate that the vB_VspP_pVa5 is a potential candidate component in a bacteriophage cocktail suitable for the biological control of V. splendidus in aquaculture.

Putative type 1 thymidylate synthase and dihydrofolate reductase as signature genes of a novel Bastille-like group of phages in the subfamily Spounavirinae

Background

Spounavirinae viruses have received an increasing interest as tools for the control of harmful bacteria due to their relatively broad host range and strictly virulent phenotype.

Results

In this study, we collected and analyzed the complete genome sequences of 61 published phages, either ICTV-classified or candidate members of the Spounavirinae subfamily of the Myoviridae. A set of comparative analyses identified a distinct, recently proposed Bastille-like phage group within the Spounavirinae. More importantly, type 1 thymidylate synthase (TS1) and dihydrofolate reductase (DHFR) genes were shown to be unique for the members of the proposed Bastille-like phage group, and are suitable as molecular markers. We also show that the members of this group encode beta-lactamase and/or sporulation-related SpoIIIE homologs, possibly questioning their suitability as biocontrol agents.

Conclusions

We confirm the creation of a new genus—the “Bastille-like group”—in Spounavirinae, and propose that the presence of TS1- and DHFR-encoding genes could serve as signatures for the new Bastille-like group. In addition, the presence of metallo-beta-lactamase and/or SpoIIIE homologs in all members of Bastille-like group phages makes questionable their suitability for use in biocontrol.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1757-0) contains supplementary material, which is available to authorized users.

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.

Learning to live together: mutualism between self-splicing introns and their hosts

Group I and II introns can be considered as molecular parasites that interrupt protein-coding and structural RNA genes in all domains of life. They function as self-splicing ribozymes and thereby limit the phenotypic costs associated with disruption of a host gene while they act as mobile DNA elements to promote their spread within and between genomes. Once considered purely selfish DNA elements, they now seem, in the light of recent work on the molecular mechanisms regulating bacterial and phage group I and II intron dynamics, to show evidence of co-evolution with their hosts. These previously underappreciated relationships serve the co-evolving entities particularly well in times of environmental stress.

Better late than early: delayed translation of intron-encoded endonuclease I-TevI is required for efficient splicing of its host group I intron

The td group I intron interrupting the thymidylate synthase (TS) gene of phage T4 is a mobile intron that encodes the homing endonuclease I-TevI. Efficient RNA splicing of the intron is required to restore function of the TS gene, while expression of I-TevI from within the intron is required to initiate intron mobility. Three distinct layers of regulation temporally limit I-TevI expression to late in the T4 infective cycle, yet the biological rationale for stringent regulation has not been tested. Here, we deleted key control elements to deregulate I-TevI expression at early and middle times post T4 infection. Strikingly, we found that deregulation of I-TevI, or of a catalytically inactive variant, generated a thymidine-dependent phenotype that is caused by a reduction in td intron splicing. Prematurely terminating I-TevI translation restores td splicing, full-length TS synthesis, and rescues the thymidine-dependent phenotype. We suggest that stringent translational control of I-TevI evolved to prevent the ribosome from disrupting key structural elements of the td intron that are required for splicing and TS function at early and middle times post T4 infection. Analogous translational regulatory mechanisms in unrelated intron-open reading frame arrangements may also function to limit deleterious consequences on splicing and host gene function.

The Widespread Evolutionary Significance of Viruses

In the last 30 years, the study of virus evolution has undergone a transformation. Originally concerned with disease and its emergence, virus evolution had not been well integrated into the general study of evolution. This chapter reviews the developments that have brought us to this new appreciation for the general significance of virus evolution to all life. We now know that viruses numerically dominate all habitats of life, especially the oceans. Theoretical developments in the 1970s regarding quasispecies, error rates, and error thresholds have yielded many practical insights into virus–host dynamics. The human diseases of HIV-1 and hepatitis C virus cannot be understood without this evolutionary framework. Yet recent developments with poliovirus demonstrate that viral fitness can be the result of a consortia, not one fittest type, a basic Darwinian concept in evolutionary biology. Darwinian principles do apply to viruses, such as with Fisher population genetics, but other features, such as reticulated and quasispecies-based evolution distinguish virus evolution from classical studies. The available phylogenetic tools have greatly aided our analysis of virus evolution, but these methods struggle to characterize the role of virus populations. Missing from many of these considerations has been the major role played by persisting viruses in stable virus evolution and disease emergence. In many cases, extreme stability is seen with persisting RNA viruses. Indeed, examples are known in which it is the persistently infected host that has better survival. We have also recently come to appreciate the vast diversity of phage (DNA viruses) of prokaryotes as a system that evolves by genetic exchanges across vast populations (Chapter 10). This has been proposed to be the “big bang” of biological evolution. In the large DNA viruses of aquatic microbes we see surprisingly large, complex and diverse viruses. With both prokaryotic and eukaryotic DNA viruses, recombination is the main engine of virus evolution, and virus host co-evolution is common, although not uniform. Viral emergence appears to be an unending phenomenon and we can currently witness a selective sweep by retroviruses that infect and become endogenized in koala bears.

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.

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.

Intron-encoded homing endonuclease I-TevI also functions as a transcriptional autorepressor

Customary binding sites of intron-encoded homing endonucleases lie within cognate intronless alleles, at the so-called homing sites. Here, we describe a novel, high-affinity binding site for I-TevI endonuclease, encoded within the group I td intron of phage T4. This site is an operator that overlaps the T4 late promoter, which drives I-TevI expression from within the td intron. I-TevI binds the operator and homing sites with equal affinity, and functions as a transcriptional autorepressor. Distinct sequence and spacing requirements of the catalytic domain result in reduced cleavage activity on operator DNA. Crystallographic studies showed that the overall interactions of the DNA-binding domain with the operator and homing sites are similar, but have some different hydrogen-bonding contacts. We present a model in which the flexibility in protein-DNA interactions allows I-TevI to bind variant intronless alleles to promote intron mobility while facilitating its function in autorepression, and thereby persistence in its host.

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.

The nicking homing endonuclease I-BasI is encoded by a group I intron in the DNA polymerase gene of the Bacillus thuringiensis phage Bastille

Here we describe the discovery of a group I intron in the DNA polymerase gene of Bacillus thuringiensis phage Bastille. Although the intron insertion site is identical to that of the Bacillus subtilis phages SPO1 and SP82 introns, the Bastille intron differs from them substantially in primary and secondary structure. Like the SPO1 and SP82 introns, the Bastille intron encodes a nicking DNA endonuclease of the H-N-H family, I-BasI, with a cleavage site identical to that of the SPO1-encoded enzyme I-HmuI. Unlike I-HmuI, which nicks both intron-minus and intron-plus DNA, I-BasI cleaves only intron-minus alleles, which is a characteristic of typical homing endonucleases. Interestingly, the C-terminal portions of these H-N-H phage endonucleases contain a conserved sequence motif, the intron-encoded endonuclease repeat motif (IENR1) that also has been found in endonucleases of the GIY-YIG family, and which likely comprises a small DNA-binding module with a globular betabetaalphaalphabeta fold, suggestive of module shuffling between different homing endonuclease families.

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.

Group I self-splicing intron in the recA gene of Bacillus anthracis

Self-splicing introns are rarely found in bacteria and bacteriophages. They are classified into group I and II according to their structural features and splicing mechanisms. While the group I introns are occasionally found in protein-coding regions of phage genomes and in several tRNA genes of cyanobacteria and proteobacteria, they had not been found in protein-coding regions of bacterial genomes. Here we report a group I intron in the recA gene of Bacillus anthracis which was initially found by DNA sequencing as an intervening sequence (IVS). By using reverse transcriptase PCR, the IVS was shown to be removable from the recA precursor mRNA for RecA that was being translated in E. coli. The splicing was visualized in vitro with labeled free GTP, indicating that it is a group I intron, which is also implied by its predicted secondary structure. The RecA protein of B. anthracis expressed in E. coli was functional in its ability to complement a recA defect. When recA-negative E. coli cells were irradiated with UV, the Bacillus RecA reduced the UV susceptibility of the recA mutant, regardless of the presence of intron.

Evolution of rbcL group IA introns and intron open reading frames within the colonial Volvocales (Chlorophyceae)

Mobile group I introns sometimes contain an open reading frame (ORF) possibly encoding a site-specific DNA endonuclease. However, previous phylogenetic studies have not clearly deduced the evolutionary roles of the group I intron ORFs. In this paper, we examined the phylogeny of group IA2 introns inserted in the position identical to that of the chloroplast-encoded rbcL coding region (rbcL-462 introns) and their ORFs from 13 strains of five genera (Volvox, Pleodorina, Volvulina, Astrephomene, and Gonium) of the colonial Volvocales (Chlorophyceae) and a related unicellular green alga, Vitreochlamys. The rbcL-462 introns contained an intact or degenerate ORF of various sizes except for the Gonium multicoccum rbcL-462 intron. Partial amino acid sequences of some rbcL-462 intron ORFs exhibited possible homology to the endo/excinuclease amino acid terminal domain. The distribution of the rbcL-462 introns is sporadic in the phylogenetic trees of the colonial Volvocales based on the five chloroplast exon sequences (6021 bp). Phylogenetic analyses of the conserved intron sequences resolved that the G. multicoccum rbcL-462 intron had a phylogenetic position separate from those of other colonial volvocalean rbcL-462 introns, indicating the recent horizontal transmission of the intron in the G. multicoccum lineage. However, the combined data set from conserved intron sequences and ORFs from most of the rbcL-462 introns resolved robust phylogenetic relationships of the introns that were consistent with those of the host organisms. Therefore, most of the extant rbcL-462 introns may have been vertically inherited from the common ancestor of their host organisms, whereas such introns may have been lost in other lineages during evolution of the colonial Volvocales. In addition, apparently higher synonymous substitutions than nonsynonymous substitutions in the rbcL-462 intron ORFs indicated that the ORFs might evolve under functional constraint, which could result in homing of the rbcL-462 intron in cases of spontaneous intron loss. On the other hand, the presence of intact to largely degenerate ORFs of the rbcL-462 introns within the three isolates of Gonium viridistellatum and the rare occurrence of the ORF-lacking rbcL-462 intron suggested that the ORFs might degenerate to result in the spontaneous intron loss during a very short evolutionary time following the loss of the ORF function. Thus, the sporadic distribution of the rbcL-462 introns within the colonial Volvocales can be largely explained by an equilibrium between maintenance of the introns by the intron ORF and spontaneous loss of introns when the introns do not have a functional ORF.

Two self-splicing group I introns in the ribonucleotide reductase large subunit gene of Staphylococcus aureus phage Twort

We have recently described three group I introns inserted into a single gene, orf142, of the staphylococcal bacteriophage Twort and suggested the presence of at least two additional self-splicing introns in this phage genome. Here we report that two previously uncharacterized introns, 429 and 1087 nt in length, interrupt the Twort gene coding for the large subunit of ribonucleotide reductase (nrdE). Reverse transcription-polymerase chain reaction (RT-PCR) of RNA isolated from Staphylococcus aureus after phage infection indicates that the introns are removed from the primary transcript in vivo. Both nrdE introns show sequence similarity to the Twort orf142 introns I2 and I3, suggesting either a common origin of these introns or shuffling of intron structural elements. Intron 2 encodes a DNA endonuclease, I-TwoI, with similarity to homing endonucleases of the HNH family. Like I-HmuI and I-HmuII, intron-encoded HNH endonucleases in Bacillus subtilis phages SPO1 and SP82, I-TwoI nicks only one strand of its DNA recognition sequence. However, whereas I-HmuI and I-HmuII cleave the template strand in exon 2, I-TwoI cleaves the coding strand in exon 1. In each case, the 3' OH created on the cut strand is positioned to prime DNA synthesis towards the intron, suggesting that this reaction contributes to the mechanism of intron homing. Both nrdE introns are inserted in highly conserved regions of the ribonucleotide reductase gene, next to codons for functionally important residues.

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.

Related homing endonucleases I-BmoI and I-TevI use different strategies to cleave homologous recognition sites

A typical homing endonuclease initiates mobility of its group I intron by recognizing DNA both upstream and downstream of the intron insertion site of intronless alleles, preventing the endonuclease from binding and cleaving its own intron-containing allele. Here, we describe a GIY-YIG family homing endonuclease, I-BmoI, that possesses an unusual recognition sequence, encompassing 1 base pair upstream but 38 base pairs downstream of the intron insertion site. I-BmoI binds intron-containing and intronless substrates with equal affinity but can nevertheless discriminate between the two for cleavage. I-BmoI is encoded by a group I intron that interrupts the thymidylate synthase (TS) gene (thyA) of Bacillus mojavensis s87-18. This intron resembles one inserted 21 nucleotides further downstream in a homologous TS gene (td) of Escherichia coli phage T4. I-TevI, the T4 td intron-encoded GIY-YIG endonuclease, is very similar to I-BmoI, but each endonuclease gene is inserted within a different position of its respective intron. Remarkably, I-TevI and I-BmoI bind a homologous stretch of TS-encoding DNA and cleave their intronless substrates in very similar positions. Our results suggest that each endonuclease has independently evolved the ability to distinguish intron-containing from intronless alleles while maintaining the same conserved recognition sequence centered on DNA-encoding active site residues of TS.

IS1414, an Escherichia coli insertion sequence with a heat-stable enterotoxin gene embedded in a transposase-like gene

Enteroaggregative Escherichia coli (EAEC) heat-stable enterotoxin 1 (EAST1) was originally discovered in EAEC but has also been associated with enterotoxigenic E. coli (ETEC). Multiple genomic restriction fragments from each of three ETEC strains of human origin showed homology with an EAST1 gene probe. A single hybridizing fragment was detected on the plasmid of ETEC strain 27D that also encodes heat-stable enterotoxin Ib and colonization factor antigen I. We isolated and characterized this fragment, showing that it (i) carries an allele of astA nearly identical to that originally reported from EAEC 17-2 and (ii) expressed enterotoxic activity. Sequence analysis of the toxin coding region revealed that astA is completely embedded within a 1,209-bp open reading frame (ORF1), whose coding sequence is on the same strand but in the -1 reading frame in reference to the toxin gene. In vitro expression of the predicted M(r)- approximately 46,000 protein product of ORF1 was demonstrated. ORF1 is highly similar to transposase genes of IS285 from Yersinia pestis, IS1356 from Burkholderia cepacia, and ISRm3 from Rhizobium meliloti. It is bounded by 30-bp imperfect inverted repeat sequences and flanked by 8-bp direct repeats. Based on these structural features, pathognomonic of a regular insertion sequence, this element was designated IS1414. Preliminary experiments to show IS1414 translocation were unsuccessful. Overlapping genes of the type suggested by the IS1414 core region have heretofore not been described in bacteria. It seems to offer a most efficient mechanism for intragenomic and horizontal dissemination of EAST1.

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.

Unexpected abundance of self-splicing introns in the genome of bacteriophage Twort: introns in multiple genes, a single gene with three introns, and exon skipping by group I ribozymes

Analysis of RNA that can be labeled with GTP indicates the existence of group I introns in genes of at least three transcriptional classes in the genome of Staphylococcus aureus bacteriophage Twort. A single ORF of 142 amino acids (Orf142) is interrupted by three self-splicing group I introns, providing the first example of a phage gene with multiple intron insertions. Twort Orf142 is encoded in a message that is abundant 15-20 min after infection and is highly similar to a late gene product (Orf8) of the morphologically related Listeria phage A511. The introns in orf142 are spliced in vivo and contain all the conserved features of primary sequence and secondary structure of group I introns in subgroup IA2, which includes the introns in Escherichia coli phage T4 and the Bacillus phages beta22 and SPO1. Introns I2 and I3 in orf142 are highly similar, and their intron insertion sites are closely spaced. The presence of transcripts with a skipped exon between these introns indicates that they may fold into a single active ribozyme resulting in alternative splicing. Alternatively, the cleaved 5' exon preceding I2 may undergo trans splicing to the 3' exon that follows I3. Regardless of the detailed mechanism, these results demonstrate a new means whereby a single gene can give rise to multiple messenger RNAs.

Characterization of rbcL group IA introns from two colonial volvocalean species (Chlorophyceae)

Group I introns were reported for the first time in the large subunit of Rubisco (rbcL) genes, using two colonial green algae, Pleodorina californica and Gonium multicoccum (Volvocales). The rbcL gene of P. californica contained an intron (PIC intron) of 1320 bp harboring an open reading frame (ORF). The G. multicoccum rbcL gene had two ORF-lacking introns of 549 (GM1 intron) and 295 (GM2 intron) base pairs. Based on the conserved nucleotide sequences of the secondary structure, the PIC and GM1 introns were assigned to group IA2 whereas the GM2 intron belonged to group IA1. Southern hybridization analyses of nuclear and chloroplast DNAs indicated that such intron-containing rbcL genes are located in the chloroplast genome. Sequencing RNAs from the two algae revealed that these introns are spliced out during mRNA maturation. In addition, the PIC and GM1 introns were inserted in the same position of the rbcL exons, and phylogenetic analysis of group IA introns indicated a close phylogenetic relationship between the PIC and GM1 introns within the lineage of bacteriophage group IA2 introns. However, P. californica and G. multicoccum occupy distinct clades in the phylogenetic trees of the colonial Volvocales, and the majority of other colonial volvocalean species do not have such introns in the rbcL genes. Therefore, these introns might have been recently inserted in the rbcL genes independently by horizontal transmission by viruses or bacteriophage.

Non-canonical inteins

Previous analyses have shown that inteins (protein splicing elements) employ two structural organizations: the 'canonical' Nintein-Dod-inteinC found in dozens of inteins and a 'non-canonical' Nintein-inteinC described in two inteins, where Nintein at the N-terminus and inteinC at the C-terminus are conserved domains involved in self-splicing and Dod is the Dod DNA endonuclease (DNase). In this study, four non-canonical inteins, each with unique structural features, have been identified using alignment-based Hidden Markov Models. A Nintein-inteinC intein, carrying an unprecedented replacement of the N-terminal catalytic Cys(Ser) by Ala, is described in a putative ATPase encoded by Methanococcus jannaschii . Three replicative proteins of Synechocystis spp. contain inteins with the organizations: (i) Nintein minus X minus inteinC over Dod, where X is an uncharacterized domain and Dod DNase is located in an alternative open reading frame (ORF) being embedded between two novel CG and YK domains; (ii) Nintein-HN-inteinC, where HN stands for phage-like DNase from the EX1H-HX3H family; (iii) Nintein>|<inteinC, where >|< indicates that the intein domains are associated with a disrupted host protein encoded by two spatially separated ORFs. The expression of some of these newly identified inteins may affect the intein hosts. The variety of structural forms of inteins could have evolved through invasion of self-splicing proteases by different mobile DNases or the departure of mobile DNases from canonical inteins.

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

The two putative ribonucleotide reductase subunits of the Bacillus subtilis bacteriophage SPbeta are encoded by the bnrdE and bnrdF genes that are highly similar to corresponding host paralogs, located on the opposite replication arm. In contrast to their bacterial counterparts, bnrdE and bnrdF each are interrupted by a group I intron, efficiently removed in vivo by mRNA processing. The bnrdF intron contains an ORF encoding a polypeptide similar to homing endonucleases responsible for intron mobility, whereas the bnrdE intron has no obvious trace of coding sequence. The downstream bnrdE exon harbors an intervening sequence not excised at the level of the primary transcript, which encodes an in-frame polypeptide displaying all the features of an intein. Presently, this is the only intein identified in bacteriophages. In addition, bnrdE provides an example of a group I intron and an intein coding sequence within the same gene.

Multiplex sequencing of 1.5 Mb of the Mycobacterium leprae genome

The nucleotide sequence of 1.5 Mb of genomic DNA from Mycobacterium leprae was determined using computer-assisted multiplex sequencing technology. This brings the 2.8-Mb M. leprae genome sequence to approximately 66% completion. The sequences, derived from 43 recombinant cosmids, contain 1046 putative protein-coding genes, 44 repetitive regions, 3 tRNAs, and 15 tRNAs. The gene density of one per 1.4 kb is slightly lower than that of Mycoplasma (1.2 kb). Of the protein coding genes, 44% have significant matches to genes with well-defined functions. Comparison of 1157 M. leprae and 1564 Mycobacterium tuberculosis proteins shows a complex mosaic of homologous genomic blocks with up to 22 adjacent proteins in conserved map order. Matches to known enzymatic, antigenic, membrane, cell wall, cell division, multidrug resistance, and virulence proteins suggest therapeutic and vaccine targets. Unusual features of the M. leprae genome include large polyketide synthase (pks) operons, inteins, and highly fragmented pseudogenes.

Rules for RNA recognition of GNRA tetraloops deduced by in vitro selection: comparison with in vivo evolution

Terminal loops with a GNRA consensus sequence are a prominent feature of large self-assembling RNA molecules. In order to investigate tertiary interactions involving GNRA loops, we have devised an in vitro selection system derived from a group I ribozyme. Two selections, destined to isolate RNA sequences that would recognize two of the most widespread loops (GUGA and GAAA), yielded variants of previously identified receptors for those loops, and also some yet unrecognized, high-affinity binders with novel specificities towards members of the GNRA family. By taking advantage of available crystal structures, we have attempted to rationalize these results in terms of RNA-RNA contacts and to expose some of the structural principles that govern GNRA loop-mediated tertiary interactions; the role of loop nucleotide 2 in ensuring specific recognition by receptors is emphasized. More generally, comparison of the products of in vitro and natural selection is shown to provide insights into the mechanisms underlying the in vivo evolution of self-assembling RNA molecules.

Mitochondrial intronic open reading frames in Podospora: mobility and consecutive exonic sequence variations

The mitochondrial genome of 23 wild-type strains belonging to three different species of the filamentous fungus Podospora was examined. Among the 15 optional sequences identified are two intronic reading frames, nad1-i4-orf1 and cox1-i7-orf2. We show that the presence of these sequences was strictly correlated with tightly clustered nucleotide substitutions in the adjacent exon. This correlation applies to the presence or absence of closely related open reading frames (ORFs), found at the same genetic locations, in all the Pyrenomycete genera examined. The recent gain of these optional ORFs in the evolution of the genus Podospora probably account for such sequence differences. In the homoplasmic progeny from heteroplasmons constructed between Podospora strains differing by the presence of these optional ORFs, nad1-i4-orf1 and cox1-i7-orf2 appeared highly invasive. Sequence comparisons in the nad1-i4 intron of various strains of the Pyrenomycete family led us to propose a scenario of its evolution that includes several events of loss and gain of intronic ORFs. These results strongly reinforce the idea that group 1 intronic ORFs are mobile elements and that their transfer, and concomitant modification of the adjacent exon, could participate in the modular evolution of mitochondrial genomes.

A general method for the consecutive integration of single copies of a heterologous gene at multiple locations in the Bacillus subtilis chromosome by replacement recombination

We have devised a two-step procedure by which multiple copies of a heterologous gene can be consecutively integrated into the Bacillus subtilis 168 chromosome without the simultaneous integration of markers (antibiotic resistance). The procedure employs the high level of transformability of B. subtilis 168 strains and makes use of the observation that thymine-auxotrophic mutants of B. subtilis are resistant to the folic acid antagonist trimethoprim (Tmpr), whereas thymine prototrophs are sensitive. First, a thymine-auxotrophic B. subtilis mutant is transformed to prototrophy by integration of a thymidylate synthetase-encoding gene at the desired chromosomal locus. In a second step, the mutant strain is transformed with a DNA fragment carrying the heterologous gene and Tmpr colonies are selected. Approximately 5% of these appear to be thymine auxotrophic and contain a single copy of the heterologous gene at the chromosomal locus previously carrying the thymidylate synthetase-encoding gene. Repetition of the procedure at different locations on the bacterial chromosome allows the isolation of strains carrying multiple copies of the heterologous gene. The method was used to construct B. subtilis strains carrying one, two, and three copies of the Bacillus stearothermophilus branching enzyme gene (glgB) in their genomes.