Evidence for horizontal transfer of the EcoT38I restriction-modification gene to chromosomal DNA by the P2 phage and diversity of defective P2 prophages in Escherichia coli TH38 strains

Going viral: The role of mobile genetic elements in bacterial immunity

Bacteriophages and other mobile genetic elements (MGEs) pose a significant threat to bacteria, subjecting them to constant attacks. In response, bacteria have evolved a sophisticated immune system that employs diverse defensive strategies and mechanisms. Remarkably, a growing body of evidence suggests that most of these defenses are encoded by MGEs themselves. This realization challenges our traditional understanding of bacterial immunity and raises intriguing questions about the evolutionary forces at play. Our review provides a comprehensive overview of the latest findings on the main families of MGEs and the defense systems they encode. We also highlight how a vast diversity of defense systems remains to be discovered and their mechanism of mobility understood. Altogether, the composition and distribution of defense systems in bacterial genomes only makes sense in the light of the ecological and evolutionary interactions of a complex network of MGEs.

The role of rhizosphere phages in soil health

While the One Health framework has emphasized the importance of soil microbiomes for plant and human health, one of the most diverse and abundant groups-bacterial viruses, i.e. phages-has been mostly neglected. This perspective reviews the significance of phages for plant health in rhizosphere and explores their ecological and evolutionary impacts on soil ecosystems. We first summarize our current understanding of the diversity and ecological roles of phages in soil microbiomes in terms of nutrient cycling, top-down density regulation, and pathogen suppression. We then consider how phages drive bacterial evolution in soils by promoting horizontal gene transfer, encoding auxiliary metabolic genes that increase host bacterial fitness, and selecting for phage-resistant mutants with altered ecology due to trade-offs with pathogen competitiveness and virulence. Finally, we consider challenges and avenues for phage research in soil ecosystems and how to elucidate the significance of phages for microbial ecology and evolution and soil ecosystem functioning in the future. We conclude that similar to bacteria, phages likely play important roles in connecting different One Health compartments, affecting microbiome diversity and functions in soils. From the applied perspective, phages could offer novel approaches to modulate and optimize microbial and microbe-plant interactions to enhance soil health.

Mining bacterial NGS data vastly expands the complete genomes of temperate phages

Temperate phages (active prophages induced from bacteria) help control pathogenicity, modulate community structure, and maintain gut homeostasis. Complete phage genome sequences are indispensable for understanding phage biology. Traditional plaque techniques are inapplicable to temperate phages due to their lysogenicity, curbing their identification and characterization. Existing bioinformatics tools for prophage prediction usually fail to detect accurate and complete temperate phage genomes. This study proposes a novel computational temperate phage detection method (TemPhD) mining both the integrated active prophages and their spontaneously induced forms (temperate phages) from next-generation sequencing raw data. Applying the method to the available dataset resulted in 192 326 complete temperate phage genomes with different host species, expanding the existing number of complete temperate phage genomes by more than 100-fold. The wet-lab experiments demonstrated that TemPhD can accurately determine the complete genome sequences of the temperate phages, with exact flanking sites, outperforming other state-of-the-art prophage prediction methods. Our analysis indicates that temperate phages are likely to function in the microbial evolution by (i) cross-infecting different bacterial host species; (ii) transferring antibiotic resistance and virulence genes and (iii) interacting with hosts through restriction-modification and CRISPR/anti-CRISPR systems. This work provides a comprehensively complete temperate phage genome database and relevant information, which can serve as a valuable resource for phage research.

Comparative Genomic Analysis of Rhodococcus equi: An Insight into Genomic Diversity and Genome Evolution

Rhodococcus equi, a member of the Rhodococcus genus, is a gram-positive pathogenic bacterium. Rhodococcus possesses an open pan-genome that constitutes the basis of its high genomic diversity and allows for adaptation to specific niche conditions and the changing host environments. Our analysis further showed that the core genome of R. equi contributes to the pathogenicity and niche adaptation of R. equi. Comparative genomic analysis revealed that the genomes of R. equi shared identical collinearity relationship, and heterogeneity was mainly acquired by means of genomic islands and prophages. Moreover, genomic islands in R. equi were always involved in virulence, resistance, or niche adaptation and possibly working with prophages to cause the majority of genome expansion. These findings provide an insight into the genomic diversity, evolution, and structural variation of R. equi and a valuable resource for functional genomic studies.

Prokaryotic horizontal gene transfer within the human holobiont: ecological-evolutionary inferences, implications and possibilities

The ubiquity of horizontal gene transfer in the living world, especially among prokaryotes, raises interesting and important scientific questions regarding its effects on the human holobiont i.e., the human and its resident bacterial communities considered together as a unit of selection. Specifically, it would be interesting to determine how particular gene transfer events have influenced holobiont phenotypes in particular ecological niches and, conversely, how specific holobiont phenotypes have influenced gene transfer events. In this synthetic review, we list some notable and recent discoveries of horizontal gene transfer among the prokaryotic component of the human microbiota, and analyze their potential impact on the holobiont from an ecological-evolutionary viewpoint. Finally, the human-Helicobacter pylori association is presented as an illustration of these considerations, followed by a delineation of unresolved questions and avenues for future research.

Occurrence, integrity and functionality of AcaML1-like viruses infecting extreme acidophiles of the Acidithiobacillus species complex

General knowledge on the diversity and biology of microbial viruses infecting bacterial hosts from extreme acidic environments lags behind most other econiches. In this study, we analyse the AcaML1 virus occurrence in the taxon, its genetic composition and infective behaviour under standard acidic and SOS-inducing conditions to assess its integrity and functionality. Occurrence analysis in sequenced acidithiobacilli showed that AcaML1-like proviruses are confined to the mesothermophiles Acidithiobacillus caldus and Thermithiobacillus tepidarius. Among A. caldus strains and isolates this provirus had a modest prevalence (30%). Comparative genomic analysis revealed a significant conservation with the T. tepidarius AcaML1-like provirus, excepting the tail genes, and a high conservation of the virus across strains of the A. caldus species. Such conservation extends from the modules architecture to the gene level, suggesting that organization and composition of these viruses are preserved for functional reasons. Accordingly, the AcaML1 proviruses were demonstrated to excise from their host genomes under DNA-damaging conditions triggering the SOS-response and to produce DNA-containing VLPs. Despite this fact, under the conditions evaluated (acidic) the VLPs obtained from A. caldus ATCC 51756 could not produce productive infections of a candidate sensitive strain (#6) nor trigger it lysis.

A Type 3 Prophage of 'Candidatus Liberibacter asiaticus' Carrying a Restriction-Modification System

Prophages, the lysogenic form of bacterial phages, are important genetic entities of 'Candidatus Liberibacter asiaticus' (CLas), a nonculturable α-proteobacterium associated with citrus Huanglongbing. Two CLas prophages have been described, SC1 (NC_019549.1, Type 1) and SC2 (NC_019550.1, Type 2), which involve the lytic cycle and the lysogenic cycle, respectively. To explore the prophage repertoire, 523 CLas DNA samples extracted from leaf petioles of CLas-infected citrus were collected from southern China and surveyed for Type 1 and Type 2 prophages by specific PCR. Eighteen samples were found lacking both prophages. One sample, JXGC, sequenced using Illumina HiSeq, generated >100 million short sequence reads (150 bp per read). Read mapping to known prophage sequences showed a sequence coverage of 46% to SC1 and 50% to SC2. BLAST search using SC1 and SC2 as queries identified three contigs from the JXGC de novo assembly that form a circular P-JXGC-3 (31,449 bp), designated as a new Type 3 prophage. Chromosomal integration of P-JXGC-3 was detected to occur within a helicase gene, resulting in a duplication of this gene. P-JXGC-3 had 36 open reading frames (ORFs), 10 of which were not found in Type 1 or Type 2 prophages, including four genes that encoded a restriction-modification (R-M) system (hsdR, hsdS, hsdM1, and hsdM2). Typed by prophage-specific PCR, the CLas strains in southern China contained all combinations of the three prophage types with the exception of a Type 2-Type 3 combination, suggesting active ongoing prophage-phage interactions. Based on gene annotation, P-JXGC-3 is not capable of reproduction via the lytic cycle. The R-M system was speculated to play a role against Type 1 prophage-phage invasion.

Prophage-mediated defence against viral attack and viral counter-defence

Temperate phages are common, and prophages are abundant residents of sequenced bacterial genomes. Mycobacteriophages are viruses that infect mycobacterial hosts including Mycobacterium tuberculosis and Mycobacterium smegmatis, encompass substantial genetic diversity and are commonly temperate. Characterization of ten Cluster N temperate mycobacteriophages revealed at least five distinct prophage-expressed viral defence systems that interfere with the infection of lytic and temperate phages that are either closely related (homotypic defence) or unrelated (heterotypic defence) to the prophage. Target specificity is unpredictable, ranging from a single target phage to one-third of those tested. The defence systems include a single-subunit restriction system, a heterotypic exclusion system and a predicted (p)ppGpp synthetase, which blocks lytic phage growth, promotes bacterial survival and enables efficient lysogeny. The predicted (p)ppGpp synthetase coded by the Phrann prophage defends against phage Tweety infection, but Tweety codes for a tetrapeptide repeat protein, gp54, which acts as a highly effective counter-defence system. Prophage-mediated viral defence offers an efficient mechanism for bacterial success in host-virus dynamics, and counter-defence promotes phage co-evolution.

Whole genomic DNA sequencing and comparative genomic analysis of Arthrospira platensis: high genome plasticity and genetic diversity

Arthrospira platensis is a multi-cellular and filamentous non-N2-fixing cyanobacterium that is capable of performing oxygenic photosynthesis. In this study, we determined the nearly complete genome sequence of A. platensis YZ. A. platensis YZ genome is a single, circular chromosome of 6.62 Mb in size. Phylogenetic and comparative genomic analyses revealed that A. platensis YZ was more closely related to A. platensis NIES-39 than Arthrospira sp. PCC 8005 and A. platensis C1. Broad gene gains were identified between A. platensis YZ and three other Arthrospira speices, some of which have been previously demonstrated that can be laterally transferred among different species, such as restriction-modification systems-coding genes. Moreover, unprecedented extensive chromosomal rearrangements among different strains were observed. The chromosomal rearrangements, particularly the chromosomal inversions, were analysed and estimated to be closely related to palindromes that involved long inverted repeat sequences and the extensively distributed type IIR restriction enzyme in the Arthrospira genome. In addition, species from genus Arthrospira unanimously contained the highest rate of repetitive sequence compared with the other species of order Oscillatoriales, suggested that sequence duplication significantly contributed to Arthrospira genome phylogeny. These results provided in-depth views into the genomic phylogeny and structural variation of A. platensis, as well as provide a valuable resource for functional genomics studies.

Pectobacterium atrosepticum and Pectobacterium carotovorum Harbor Distinct, Independently Acquired Integrative and Conjugative Elements Encoding Coronafacic Acid that Enhance Virulence on Potato Stems

Integrative and conjugative elements (ICEs) play a central role in the evolution of bacterial virulence, their transmission between bacteria often leading to the acquisition of virulence factors that alter host range or aggressiveness. Much is known about the functions of the virulence determinants that ICEs harbor, but little is understood about the cryptic effects of ICEs on their host cell. In this study, the importance of horizontally acquired island 2 (HAI2), an ICE in the genome of Pectobacterium atrosepticum SCRI1043, was studied using a strain in which the entire ICE had been removed by CRISPR-Cas-mediated genome editing. HAI2 encodes coronafacic acid, a virulence factor that enhances blackleg disease of potato stems caused by P. atrosepticum SCRI1043. As expected, deletion of HAI2 resulted in reduced blackleg symptoms in potato stems. A subsequent screen for HAI2-related ICEs in other strains of the Pectobacterium genus revealed their ubiquitous nature in P. atrosepticum. Yet, HAI2-related ICEs were only detected in the genomes of a few P. carotovorum strains. These strains were notable as blackleg causing strains belonging to two different subspecies of P. carotovorum. Sequence analysis of the ICEs in different strains of both P. atrosepticum and P. carotovorum confirmed that they were diverse and were present in different locations on the genomes of their bacterial host, suggesting that the cfa cluster was probably acquired independently on a number of occasions via chromosomal insertion of related ICEs. Excision assays also demonstrated that the ICEs in both P. atrosepticum and P. carotovorum are mobilized from the host chromosome. Thus, the future spread of these ICEs via lateral gene transfer might contribute to an increase in the prevalence of blackleg-causing strains of P. carotovorum.

The interplay of restriction-modification systems with mobile genetic elements and their prokaryotic hosts

The roles of restriction-modification (R-M) systems in providing immunity against horizontal gene transfer (HGT) and in stabilizing mobile genetic elements (MGEs) have been much debated. However, few studies have precisely addressed the distribution of these systems in light of HGT, its mechanisms and its vectors. We analyzed the distribution of R-M systems in 2261 prokaryote genomes and found their frequency to be strongly dependent on the presence of MGEs, CRISPR-Cas systems, integrons and natural transformation. Yet R-M systems are rare in plasmids, in prophages and nearly absent from other phages. Their abundance depends on genome size for small genomes where it relates with HGT but saturates at two occurrences per genome. Chromosomal R-M systems might evolve under cycles of purifying and relaxed selection, where sequence conservation depends on the biochemical activity and complexity of the system and total gene loss is frequent. Surprisingly, analysis of 43 pan-genomes suggests that solitary R-M genes rarely arise from the degradation of R-M systems. Solitary genes are transferred by large MGEs, whereas complete systems are more frequently transferred autonomously or in small MGEs. Our results suggest means of testing the roles for R-M systems and their associations with MGEs.

Bacteriophage orphan DNA methyltransferases: insights from their bacterial origin, function, and occurrence

Type II DNA methyltransferases (MTases) are enzymes found ubiquitously in the prokaryotic world, where they play important roles in several cellular processes, such as host protection and epigenetic regulation. Three classes of type II MTases have been identified thus far in bacteria which function in transferring a methyl group from S-adenosyl-l-methionine (SAM) to a target nucleotide base, forming N-6-methyladenine (class I), N-4-methylcytosine (class II), or C-5-methylcytosine (class III). Often, these MTases are associated with a cognate restriction endonuclease (REase) to form a restriction-modification (R-M) system protecting bacterial cells from invasion by foreign DNA. When MTases exist alone, which are then termed orphan MTases, they are believed to be mainly involved in regulatory activities in the bacterial cell. Genomes of various lytic and lysogenic phages have been shown to encode multi- and mono-specific orphan MTases that have the ability to confer protection from restriction endonucleases of their bacterial host(s). The ability of a phage to overcome R-M and other phage-targeting resistance systems can be detrimental to particular biotechnological processes such as dairy fermentations. Conversely, as phages may also be beneficial in certain areas such as phage therapy, phages with additional resistance to host defenses may prolong the effectiveness of the therapy. This minireview will focus on bacteriophage-encoded MTases, their prevalence and diversity, as well as their potential origin and function.

Complete genome sequence of temperate bacteriophage AcaML1 from the extreme acidophile Acidithiobacillus caldus ATCC 51756

Development of reproducible genetic tools in the industrially important acidithiobacilli is urgently required. Inducible temperate phages which may be modified in vitro, propagated in suitable hosts, and used to transduce relevant genetic information to other strains and/or species are potentially valuable tools in this field of research. In order to address these current limitations, the genome sequence of an inducible temperate Myoviridae-like bacteriophage from the Acidithiobacillus caldus type strain was annotated and analyzed bioinformatically. Here, we announce the genome sequence of AcaML1 and report major findings from its annotation.

Movement of DNA sequence recognition domains between non-orthologous proteins

Comparisons of proteins show that they evolve through the movement of domains. However, in many cases, the underlying mechanisms remain unclear. Here, we observed the movements of DNA recognition domains between non-orthologous proteins within a prokaryote genome. Restriction-modification (RM) systems, consisting of a sequence-specific DNA methyltransferase and a restriction enzyme, contribute to maintenance/evolution of genomes/epigenomes. RM systems limit horizontal gene transfer but are themselves mobile. We compared Type III RM systems in Helicobacter pylori genomes and found that target recognition domain (TRD) sequences are mobile, moving between different orthologous groups that occupy unique chromosomal locations. Sequence comparisons suggested that a likely underlying mechanism is movement through homologous recombination of similar DNA sequences that encode amino acid sequence motifs that are conserved among Type III DNA methyltransferases. Consistent with this movement, incongruence was observed between the phylogenetic trees of TRD regions and other regions in proteins. Horizontal acquisition of diverse TRD sequences was suggested by detection of homologs in other Helicobacter species and distantly related bacterial species. One of these RM systems in H. pylori was inactivated by insertion of another RM system that likely transferred from an oral bacterium. TRD movement represents a novel route for diversification of DNA-interacting proteins.

Evolutionary genome engineering using a restriction-modification system

Modification of complex microbial cellular processes is often necessary to obtain organisms with particularly favorable characteristics, but such experiments can take many generations to achieve. In the present article, we accelerated the experimental evolution of Escherichia coli populations under selection for improved growth using one of the restriction-modification systems, which have shaped bacterial genomes. This resulted in faster evolutionary changes in both the genome and bacterial growth. Transcriptome/genome analysis at various stages enabled prompt identification of sequential genome rearrangements and dynamic gene-expression changes associated with growth improvement. The changes were related to cell-to-cell communication, the cell death program, as well as mass production and energy consumption. These observed changes imply that improvements in microorganism population growth can be achieved by inactivating the cellular mechanisms regulating fraction of active cells in a population. Some of the mutations were shown to have additive effects on growth. These results open the way for the application of evolutionary genome engineering to generate organisms with desirable properties.

The host-binding domain of the P2 phage tail spike reveals a trimeric iron-binding structure

The adsorption and infection of bacteriophage P2 is mediated by tail fibres and tail spikes. The tail spikes on the tail baseplate are used to irreversibly adsorb to the host cells. Recently, a P2 phage tail-spike protein, gpV, was purified and it was shown that a C-terminal domain, Ser87-Leu211, is sufficient for the binding of gpV to host Escherichia coli membranes [Kageyama et al. (2009), Biochemistry, 48, 10129-10135]. In this paper, the crystal structure of the C-terminal domain of P2 gpV is reported. The structure is a triangular pyramid and looks like a spearhead composed of an intertwined β-sheet, a triple β-helix and a metal-binding region containing iron, calcium and chloride ions.

Domain movement within a gene: a novel evolutionary mechanism for protein diversification

A protein function is carried out by a specific domain localized at a specific position. In the present study, we report that, within a gene, a specific amino acid sequence can move between a certain position and another position. This was discovered when the sequences of restriction-modification systems within the bacterial species Helicobacter pylori were compared. In the specificity subunit of Type I restriction-modification systems, DNA sequence recognition is mediated by target recognition domain 1 (TRD1) and TRD2. To our surprise, several sequences are shared by TRD1 and TRD2 of genes (alleles) at the same locus (chromosomal location); these domains appear to have moved between the two positions. The gene/protein organization can be represented as x-(TRD1)-y-x-(TRD2)-y, where x and y represent repeat sequences. Movement probably occurs by recombination at these flanking DNA repeats. In accordance with this hypothesis, recombination at these repeats also appears to decrease two TRDs into one TRD or increase these two TRDs to three TRDs (TRD1-TRD2-TRD2) and to allow TRD movement between genes even at different loci. Similar movement of domains between TRD1 and TRD2 was observed for the specificity subunit of a Type IIG restriction enzyme. Similar movement of domain between TRD1 and TRD2 was observed for Type I restriction-modification enzyme specificity genes in two more eubacterial species, Streptococcus pyogenes and Mycoplasma agalactiae. Lateral domain movements within a protein, which we have designated DOMO (domain movement), represent novel routes for the diversification of proteins.

Genome comparison and context analysis reveals putative mobile forms of restriction-modification systems and related rearrangements

The mobility of restriction–modification (RM) gene complexes and their association with genome rearrangements is a subject of active investigation. Here we conducted systematic genome comparisons and genome context analysis on fully sequenced prokaryotic genomes to detect RM-linked genome rearrangements. RM genes were frequently found to be linked to mobility-related genes such as integrase and transposase homologs. They were flanked by direct and inverted repeats at a significantly high frequency. Insertion by long target duplication was observed for I, II, III and IV restriction types. We found several RM genes flanked by long inverted repeats, some of which had apparently inserted into a genome with a short target duplication. In some cases, only a portion of an apparently complete RM system was flanked by inverted repeats. We also found a unit composed of RM genes and an integrase homolog that integrated into a tRNA gene. An allelic substitution of a Type III system with a linked Type I and IV system pair, and allelic diversity in the putative target recognition domain of Type IIG systems were observed. This study revealed the possible mobility of all types of RM systems, and the diversity in their mobility-related organization.

The Crp-activated small noncoding regulatory RNA CyaR (RyeE) links nutritional status to group behavior

Small noncoding regulatory RNAs (sRNAs) play a key role in regulating the expression of many genes in Escherichia coli and other bacteria. Many of the sRNAs identified in E. coli bind to mRNAs in an Hfq-dependent manner and stimulate or inhibit translation of the mRNAs. Several sRNAs are regulated by well-studied global regulators. Here, we report characterization of the CyaR (RyeE) sRNA, which was previously identified in a global search for sRNAs in E. coli. We demonstrated that CyaR is positively regulated by the global regulator Crp under conditions in which cyclic AMP levels are high. We showed by using microarray analysis and Northern blotting that several genes are negatively regulated by CyaR, including ompX, encoding a major outer membrane protein; luxS, encoding the autoinducer-2 synthase; nadE, encoding an essential NAD synthetase; and yqaE, encoding a predicted membrane protein with an unknown function. Using translational lacZ fusions to yqaE, ompX, nadE, and luxS, we demonstrated that the negative regulation of these genes by CyaR occurs at the posttranscriptional level and is direct. Different portions of a highly conserved 3' region of CyaR are predicted to pair with sequences near the ribosome binding site of each of these targets; mutations in this sequence affected regulation, and compensatory mutations in the target mRNA restored regulation, confirming that there is direct regulation by the sRNA. These results provide insight into the mechanisms by which Crp negatively regulates genes such as luxS and ompX and provide a link between catabolite repression, quorum sensing, and nitrogen assimilation in E. coli.

Identification of a novel prophage-like gene cluster actively expressed in both virulent and avirulent strains of Leptospira interrogans serovar Lai

DNA microarray analysis was used to compare the differential gene expression profiles between Leptospira interrogans serovar Lai type strain 56601 and its corresponding attenuated strain IPAV. A 22-kb genomic island covering a cluster of 34 genes (i.e., genes LA0186 to LA0219) was actively expressed in both strains but concomitantly upregulated in strain 56601 in contrast to that of IPAV. Reverse transcription-PCR assays proved that the gene cluster comprised five transcripts. Gene annotation of this cluster revealed characteristics of a putative prophage-like remnant with at least 8 of 34 sequences encoding prophage-like proteins, of which the LA0195 protein is probably a putative prophage CI-like regulator. The transcription initiation activities of putative promoter-regulatory sequences of transcripts I, II, and III, all proximal to the LA0195 gene, were further analyzed in the Escherichia coli promoter probe vector pKK232-8 by assaying the reporter chloramphenicol acetyltransferase (CAT) activities. The strong promoter activities of both transcripts I and II indicated by the E. coli CAT assay were well correlated with the in vitro sequence-specific binding of the recombinant LA0195 protein to the corresponding promoter probes detected by the electrophoresis mobility shift assay. On the other hand, the promoter activity of transcript III was very low in E. coli and failed to show active binding to the LA0195 protein in vitro. These results suggested that the LA0195 protein is likely involved in the transcription of transcripts I and II. However, the identical complete DNA sequences of this prophage remnant from these two strains strongly suggests that possible regulatory factors or signal transduction systems residing outside of this region within the genome may be responsible for the differential expression profiling in these two strains.

Analyses of the vrl gene cluster in Desulfococcus multivorans: homologous to the virulence-associated locus of the ovine footrot pathogen Dichelobacter nodosus strain A198

Major parts of the virulence-associated vrl locus known from the gammaproteobacterium Dichelobacter nodosus, the causative agent of ovine footrot, were analyzed in the genome of the sulfate-reducing deltaproteobacterium Desulfococcus multivorans. In the genome of D. multivorans 13 of the 19 vrl genes described for D. nodosus are present and highly conserved with respect to gene sequence and order. The vrl locus and its flanking regions suggest a bacteriophage-mediated transfer into the genome of D. multivorans. Comparative analysis of the deduced Vrl proteins reveals a wide distribution of parts of the virulence-associated vrl locus in distantly related bacteria. Horizontal transfer is suggested as driving mechanism for the circulation of the vrl genes in bacteria. Except for the vrlBMN genes D. multivorans and Desulfovibrio desulfuricans G20 together contain all vrl genes displaying a high degree of similarity. For D. multivorans it could be shown that guanine plus cytosine (GC) content, GC skew, di-, tri- or tetranucleotide distribution did not differ between the vrl locus and its flanking sequences. This could be a hint that the vrl locus originated from a related organism or at least a genome with similar characteristics. The conspicuous high degree of conservation of the analyzed vrl genes may result from a recent transfer event or reflect a function of the vrl genes, which is still unknown and not necessarily disease associated. The latter is supported by the evidence for expression of the vrl genes in D. multivorans, which has not been described as pathogen or to be associated to any disease pattern before.

Identification of a gene encoding a functional reverse transcriptase within a highly variable locus in the P2-like coliphages

The P2-like coliphages are highly similar; the structural genes show at least 96% identity. However, at two loci they have genes believed to be horizontally transferred. We show that the genetic content at the second loci, the TO region, contains six completely different sequences with high AT contents and with different open reading frames. The product of one of them exhibits reverse transcriptase activity and blocks infection of phage T5.

Underrepresentation of short palindromes in Selenomonas ruminantium DNA: evidence for horizontal gene transfer of restriction and modification systems?

Molecular analysis of isolates of the rumen bacterium Selenomonas ruminantium revealed a high variety and frequency of site-specific (restriction) endonucleases. While all known S. ruminantium restriction and modification systems recognize hexanucleotide sequences only, consistently low counts of both 6-bp and 4-bp palindromes were found in DNA sequences of S. ruminantium. Statistical analysis indicated that there is some correlation between the degree of underrepresentation of tetranucleotide words and the number of known restriction endonucleases for a given sequence. Control analysis showed the same correlation in lambda DNA but not in human adenovirus DNA. Based on the data presented, it could be proposed that there is a much higher historical occurrence of restriction and modification systems in S. ruminantium and (or) frequent horizontal gene transfer of restriction and modification gene complexes.

Codon usage comparison of novel genes in clinical isolates of Haemophilus influenzae

A similarity statistic for codon usage was developed and used to compare novel gene sequences found in clinical isolates of Haemophilus influenzae with a reference set of 80 prokaryotic, eukaryotic and viral genomes. These analyses were performed to obtain an indication as to whether individual genes were Haemophilus-like in nature, or if they probably had more recently entered the H.influenzae gene pool via horizontal gene transfer from other species. The average and SD values were calculated for the similarity statistics from a study of the set of all genes in the H.influenzae Rd reference genome that encoded proteins of 100 amino acids or longer. Approximately 80% of Rd genes gave a statistic indicating that they were most like other Rd genes. Genes displaying codon usage statistics >1 SD above this range were either considered part of the highly expressed group of H.influenzae genes, or were considered of foreign origin. An alternative determinant for identifying genes of foreign origin was when the similarity statistics produced a value that was much closer to a non-H.influenzae reference organism than to any of the Haemophilus species contained in the reference set. Approximately 65% of the novel sequences identified in the H.influenzae clinical isolates displayed codon usages most similar to Haemophilus sp. The remaining novel sequences produced similarity statistics closer to one of the other reference genomes thereby suggesting that these sequences may have entered the H.influenzae gene pool more recently via horizontal transfer.

Evolutionary relationships of Fusobacterium nucleatum based on phylogenetic analysis and comparative genomics

Background

The phylogenetic position and evolutionary relationships of Fusobacteria remain uncertain. Especially intriguing is their relatedness to low G+C Gram positive bacteria (Firmicutes) by ribosomal molecular phylogenies, but their possession of a typical gram negative outer membrane. Taking advantage of the recent completion of the Fusobacterium nucleatum genome sequence we have examined the evolutionary relationships of Fusobacterium genes by phylogenetic analysis and comparative genomics tools.

Results

The data indicate that Fusobacterium has a core genome of a very different nature to other bacterial lineages, and branches out at the base of Firmicutes. However, depending on the method used, 35–56% of Fusobacterium genes appear to have a xenologous origin from bacteroidetes, proteobacteria, spirochaetes and the Firmicutes themselves. A high number of hypothetical ORFs with unusual codon usage and short lengths were found and hypothesized to be remnants of transferred genes that were discarded. Some proteins and operons are also hypothesized to be of mixed ancestry. A large portion of the Gram-negative cell wall-related genes seems to have been transferred from proteobacteria.

Conclusions

Many instances of similarity to other inhabitants of the dental plaque that have been sequenced were found. This suggests that the close physical contact found in this environment might facilitate horizontal gene transfer, supporting the idea of niche-specific gene pools. We hypothesize that at a point in time, probably associated to the rise of mammals, a strong selective pressure might have existed for a cell with a Clostridia-like metabolic apparatus but with the adhesive and immune camouflage features of Proteobacteria.