Previous Article | Next Article ![]()
Applied and Environmental Microbiology, February 2008, p. 861-874, Vol. 74, No. 3
0099-2240/08/$08.00+0 doi:10.1128/AEM.02099-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Université de Lyon, Lyon, F-69003, France, Université Lyon 1, Lyon, F-69003, France, CNRS, UMR 5557, Ecologie Microbienne, Villeurbanne F-69622, France, and IFR 41, Villeurbanne F-69622, France,1 Génoscope, Centre National de Séquençage, 2 rue Gaston Crémieux CP 5706, 91057 Evry Cedex, France,2 Centro de Ciencias Genómicas, Universidad Nacional Autónoma de México, AP565-A Cuernavaca, Morelos 62210, México3
Received 13 September 2007/ Accepted 25 November 2007
|
|
|---|
Ab-Cd) genome was determined as a double-stranded DNA circular molecule of 62,337 pb that encodes 95 predicted proteins. Only 14 of the predicted proteins could be assigned functions, some of which were involved in DNA processing, phage morphogenesis, and bacterial lysis. In addition, the
Ab-Cd complete genome was mapped as a prophage on a 570-kb replicon of strain A. brasilense Cd, and a region of 27.3 kb of
Ab-Cd was found to be duplicated on the 130-kb pRhico plasmid previously sequenced from A. brasilense Sp7, the parental strain of A. brasilense Cd. |
|
|---|
1030 phage particles (57). Bacteriophages are ubiquitous in nature, and it has been suggested that they are environmentally important both in controlling the number of bacteria and in facilitating gene transfer. Indeed, bacteriophages represent one of the three major mobile genetic elements that contribute significantly to horizontal gene transfer in bacterial genomes by transduction. It was calculated that 1024 bacteria could be infected per second, which implies that a continuous flow of genetic material occurs between bacteria sharing the same environments (60). Seventy percent of complete bacterial genomes investigated contain prophage sequences of >10 kb (16), and it is believed that well over 50% of bacteria isolated from various environments contain prophages or are lysogenized. The integrated viral DNA, referred to as a prophage, can harbor genes that contribute to increasing the fitness of bacterial cells, such as those encoding phage-encoded virulence factors of some pathogenic bacteria (15). Consequently, the majority of the data concerning prophages come from studies of gram-positive bacteria and gammaproteobacteria of medical and industrial interest (16), but very little information is available concerning phages isolated from commensal or symbiotic bacteria. Plant-beneficial bacteria, including some strains belonging to the genus Azospirillum, are usually referred to as plant growth-promoting rhizobacteria. Members of the genus Azospirillum represent a group of free-living alphaproteobacteria whose promotion of the growth of cereals and grasses leads to an increase of crop yield of up to 30%. The actual benefit deriving from biological nitrogen fixation has been questioned, and plant growth promotion by Azospirillum spp. seems to be due mainly to production of phytohormones (7, 20). The most abundant phytohormone produced is indole-3-acetic acid, allowing an increase in the number of lateral roots and root hairs; this results in a higher absorption of water and minerals from the soil (7). Although many studies of Azospirillum spp. are related to bacterial physiology and ecology, few genomic data are available. It has been shown that the genomes of Azospirillum strains differ in size (from 4.8 to 9.7 Mb) and in architecture, one of the main features being the presence of numerous plasmids (14, 44). Some plasmids were shown to be involved in major genomic rearrangements occurring during phenotypic switching, indicating that genomes of Azospirillum are highly dynamic (63). To date, only two temperate phages were described for this genus. The first one, termed Al-1, was isolated from Brazilian soil and is able to form plaques on A. lipoferum Br17 (24). The second one was characterized after spontaneous or mitomycin C induction in an A. brasilense Sp7 culture (28, 30).
In the frame of investigation of the role of phages in the biology of Azospirillum spp., we first studied their prevalence by examining the ability of 24 Azospirillum strains from four species to release phage particles after mitomycin C treatment. Phage particles were isolated, and morphological studies were performed by transmission electron microscopy (TEM). Phage DNA was then purified, and the molecular relatedness of phages was investigated by Southern hybridization experiments. Finally, the complete genome sequence of the lytic phage
Ab-Cd from A. brasilense Cd strain was obtained and analyzed.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Azospirillum strains used in this study
|
TEM.
The production of phage particles in Azospirillum cultures was induced as described above. Then, bacterial lysate (12 ml) was cleared by two centrifugations at 5,000 x g for 10 min at 20°C and filtration through a 0.45-µm-pore-size filter (type CA; Millipore, Billerica, MA). The phage suspension was then transferred into a polyallomer centrifuge tube (14 by 89 mm; Beckman Coulter, Roissy, France) and centrifuged at 4°C for 5 h at 75,000 x g. The phage pellet was resuspended in 30 µl of suspension buffer without pipetting. A drop of phage suspension was applied to a 300-mesh carbon/formvar-coated grid. The sample was stained with either aqueous 1% phosphotungstic acid or 2% uranyl acetate for 30 s. Excess stain was removed with a filter paper, and the grids were allowed to air dry prior to examination using a Philips CM 120 transmission electron microscope at 100 kV.
PFGE of phage DNA and hybridization experiments.
Phage DNA was extracted from purified particles as described above. A total of 100 ng of undigested DNA was deposited per well of 1% agarose gel (pulsed-field certified agarose; Bio-Rad Laboratories) in 0.5x Tris-borate-EDTA buffer (Euromedex, Mundolsheim, France). Pulsed-field gel electrophoresis (PFGE) was run using a Chef-DRIII system (Bio-Rad Laboratories) with the following parameters: 0.5x Tris-borate-EDTA running buffer at 14°C, 1 s initial switch time, 6 s final switch time, 6 V cm–1, 15 h run time, and 14°C buffer temperature. The gels were stained with ethidium bromide and visualized using Gel-doc 2000 (Bio-Rad Laboratories, Hercules, CA). For Southern blot experiments, phage DNAs were transferred to GeneScreen Plus nylon filters (PerkinElmer, Zaventem, Belgium). Phage DNA randomly labeled with [
-32P]dCTP by use of a Redi-prime DNA labeling system (Amersham plc, Little Chalfont, England) was used as a probe. Southern blotting and hybridization were carried out according to established protocols (53).
DNA sequencing and bioinformatics analyses.
The DNA template used for sequencing the A. brasilense Cd phage was purified from phage particles as described above. The template was subjected to mechanical shearing and cloned as fragments of 10 kb into low-copy-number pCNS home vector (pSU18 modified) (6). Plasmid DNAs were purified and end sequenced using dye-terminator chemistry on ABI3730xl DNA analyzer sequencers. A total of 1,500 sequences were generated from both ends of the clones. The Phred/Phrap/Consed software package (www.phrap.com) was used for sequence assembly and quality assessment (26, 27, 31). Open reading frames (ORFs) in the final genome sequences were predicted using Glimmer 2.0 software (19) and annotated with ARTEMIS 8 software (52). In our analyses, putative ORFs contained AUG (methionine), UUG (leucine), or GUG (valine) as the starting codon. The predicted proteins were compared to those listed in the NCBI protein databases by using the BLASTX algorithm (2). Searches for putative conserved domains of identified proteins were done using the INTERPROSCAN and CD-Search programs (43, 67).
Phylogenetic analysis of functional genes.
The translated amino acid sequences of putative genes encoding terminase (ORF29), a phage major capsid protein (ORF33), lysozyme (ORF53), and integrase (ORF72) were used to construct phylogenetic trees. The BLAST program was used to identify sequences homologous to the sequences of these ORFs in the draft of the A. brasilense Sp245 genome (available at http://genomics.ornl.gov/research/azo). These amino acid sequences were aligned with sequences of viruses found in the GenBank database, using the program CLUSTAL W, and phylogenetic analyses were conducted using MEGA 4.0 software (58). A Blosum 30 matrix was calculated with a gap penalty of 10.0. Distance trees were constructed with the neighbor-joining algorithm, and the bootstrap method was employed with 1,000 replicates to estimate the robustness of the tree topologies.
Plasmid content and prophage localization.
Plasmids were separated by modified Eckhardt agarose gel electrophoresis (33) as applied by Vial and coworkers (63). Plasmid sizes were estimated by comparison with those of A. brasilense Sp7 (14) and Rhizobium etli CFN42 (51). Hybridization experiments were carried out on plasmid profiles as previously described (63). To obtain a probe specific to the A. brasilense Cd phage, several ORFs were first amplified by PCR using template DNA from strain Cd as follows: ORF40 was amplified with primers F4908 (5'-AGAAGATCAACGCCAGCTTC-3') and F4909 (5'-GACTGTGTGCCTGCGTAGAA-3'), ORF50 with primers F5561 (5'-GCAAACCACAATCACCACAG-3') and F5562 (5'-AAGCATCTGGTAACCGTTGG-3'), and ORF53 with primers F4906 (5'-GGGCTGTATCTGACCGCATA-3') and F4907 (5'-CATTGACGGTGGCGTAGAC-3'). The amplification cycle consisted of an initial 5 min at 95°C and 35 cycles of 30 s at 95°C, 30 s at annealing temperature, and 30 s at 72°C followed by a final 7 min extension at 72°C. Then, ORF50 was chosen, randomly prime labeled as described above, and used as a probe to hybridize plasmid profiles. For a probe specific to the pRhico plasmid, the exoC gene was amplified with primers F5559 (5'-CCTTCAGGTCGACGAGATTC-3') and F5560 (5'-GGACTACAACGGCATCAAGA-3').
Nucleotide sequence accession number.
The accession number for the phage DNA sequence obtained in this study is CU468217.
|
|
|---|
![]() View larger version (9K): [in a new window] |
FIG. 1. Effect of mitomycin C treatment on growth of Azospirillum bacteria. (A) A. brasilense strain Cd; (B) A. lipoferum strain B518. Cell growth was followed by measurement of optical densities at 580 nm in a culture treated with 0.1 µg/ml mitomycin C (open circles) and in a control culture without treatment (closed circles). The addition of mitomycin is indicated by arrows.
|
Ab-Cd designates the phage isolated from the strain A. brasilense Cd,
Al-B518 designates the phage isolated from the strain A. lipoferum B518, etc.
![]() View larger version (109K): [in a new window] |
FIG. 2. TEM micrographs of negatively stained bacteriophage particles of Azospirillum bacteria. After 0.1 µg ml–1 mitomycin C was added, phage particles were observed in the lysates of A. brasilense Cd (A), A. brasilense Wb1 (B), and A. lipoferum B518 (C). Bars, 100 nm.
|
Ab-Cd as a probe, hybridization signals were clearly obtained exclusively with DNA of phages
Ab-Cd,
Ab-Sp7, and
Ab-Wb1 (Fig. 3B, lanes 8, 10, and 11), suggesting DNA homology between those three phage genomes. When the
Al-B518 10-kb DNA molecule was used as a probe, hybridization signals were observed only with DNA from smaller phages isolated from A. doebereinerae and A. lipoferum strains (data not shown). In addition, identical restriction patterns were observed for
Ab-Cd and
Ab-Sp7 (Fig. 4), suggesting similar architectures for these two phages. This would indeed be consistent with the genetic relatedness of these two A. brasilense strains, strain Cd being isolated as a scarlet variant of strain Sp7 after inoculation on a plant (25, 59). The restriction pattern of
Ab-Wb1 (Fig. 4, lane 4) was unrelated to those of
Ab-Cd and
Ab-Sp7 despite a strong hybridization signal; this result could be attributed either to the presence of only a few identical genes or to genome recombination. Comparison of restriction profiles of A. lipoferum and A. doebereinerae phage DNAs was not possible, as only smears were obtained whatever the endonuclease used (data not shown).
![]() View larger version (69K): [in a new window] |
FIG. 3. PFGE of DNA molecules isolated from phage particles (A) and hybridization experiments (B). DNA was isolated from the following phage particles: a phage isolated from A. doebereinerae GSF71 and named Ad-GSF71 (lane 1); phages isolated from A. lipoferum 4T ( Al-4T; lane 2), TVV3 ( Al-TVV3; lane 3), B510 ( Al-B510; lane 4), B518 ( Al-B518; lane 5), 4VI ( Al-4VI; lane 6), and 4B ( Al-4B; lane 7); and phages isolated from A. brasilense Wb1 ( Ab-Wb1; lane 8), Sp245 ( Ab-Sp245; lane 9), Cd ( Ab-Cd; lane 10), and Sp7 ( Ab-Sp7; lane 11). Lanes M, 5 kb ladder DNA size standard. The probe used for the hybridization experiment was the DNA molecule extracted from Ab-Cd (B).
|
![]() View larger version (64K): [in a new window] |
FIG. 4. Restriction analysis of phage DNAs. Lane M, DNA digested with HindIII. Phage DNA digested with EcoRI: Ab-Cd (lane 1), Ab-Sp7 (lane 2), Ab-Sp245 (lane 3), and Ab-Wb1 (lane 4).
|
Ab-Cd genome.
Ab-Cd was determined using a shotgun strategy followed by an assembly process with phrap software and annotation with ARTEMIS 8 software (52). The
Ab-Cd genome represents a double-stranded DNA circular molecule of 62,337 pb with a G+C content of 63.93%, which is slightly lower than that described for Azospirillum spp. (69% to 71%) (59). The coding sequence represents 86.35% of the genome, with approximately 1.5 gene per kpb of nucleotide sequence, indicating a close-packed genome. A total of 95 ORFs were predicted to encode proteins larger than 40 amino acids. Searches in databases indicated that 65 (68.42%) ORFs code for hypothetical proteins, 16 (16.8%) ORFs display sequence similarity to proteins with unknown functions, and only 14 (14.73%) ORFs could be assigned functions based on protein homology (Table 2 and Fig. 5). No tRNA genes have been identified. |
View this table: [in a new window] |
TABLE 2. Genome organization of A. brasilense phages Ab-Cd
|
![]() View larger version (22K): [in a new window] |
FIG. 5. Genome map of Ab-Cd and homology with pRhico. The double black circle represents the circular DNA molecule of the Ab-Cd genome. The predicted ORFs are depicted by thin arrows or arrowheads pointing in the direction of transcription and are numbered consecutively (see Table 2). ORFs encoding hypothetical proteins (no homologs in BLAST-P searches) are represented by black thin arrows. The ORFs encoding proteins with assigned functions are represented by colored thin arrows as follows: red arrows, DNA processing; blue, phage morphogenesis; yellow, bacterial lysis; purple, putative phage proteins. Homologies between the Ab-Cd sequence and contigs of pRhico are represented by thick green arrows.
|
Ab-Cd sequence (Fig. 5). Two explanations could account for this finding: (i)
Ab-Cd is integrated as a prophage in the pRhico plasmid (but this would imply that the pRhico sequence has not been entirely described), or (ii)
Ab-Cd is integrated elsewhere in the genome as a prophage and part (18.3%) of
Ab-Cd is duplicated within the pRhico. Subsequent experiments (see below) will allow retaining one of these hypotheses. ORF predictions performed in our study allowed the identification of ORFs similar to those already predicted for contigs 3 and 4 of pRhico, but additional ORFs were found (Table 2). Even if the majority of the putative ORFs identified code for hypothetical proteins or for proteins without assigned functions, some genes encode putative phage proteins that have conserved functional domains and that could be grouped in four classes: proteins involved in DNA processing, proteins involved in phage morphogenesis, proteins related to bacterial lysis, and putative phage proteins. Nine ORFs (ORF10, ORF11, ORF13, ORF19, ORF27, ORF29, ORF57, ORF72, and ORF91) could be ascribed a function related to DNA processing (see Table 2). The deduced protein encoded by ORF10, although displaying 82% identity with conserved bacterial proteins with unassigned functions (COG3750.2), possesses a conserved domain matching the tRNA-binding arm, a domain found in aminoacyl tRNA synthetases. ORF11 could encode a protein showing homology with epoxidase-like protein. A search for conserved domains revealed a putative helix-turn-helix type 3 DNA binding domain (Interpro IPR 001387), identified in proteins belonging to the Lambda repressor-like family; thus, the product encoded by ORF11 could act as a transcriptional repressor like the regulators Cro and cI that control the life cycle of the Lambda bacteriophage. Site-specific DNA methylase domains involved in DNA modification have been identified in the two putative proteins encoded by ORF13 and ORF91. The product encoded by ORF19 exhibits high similarities with deoxycytidylate deaminase-like proteins that catalyze the deamination of dCMP into dUMP. The deduced product of ORF27 is homologous to endonucleases of the HNH family found in several phages (such as LambdaSa2 from Streptococcus agalactiae and bIL170 from Lactococcus spp.). ORF29 is likely to encode the phage terminase large subunit, an enzyme that is part of a large nucleoprotein complex dedicated to packaging. The integrase function allowing the DNA recombination and phage integration process was ascribed to the product encoded by ORF72, showing best identities with alphaproteobacteria like Oceanicola granulosus HTCC2516 (35% identity; RefSeq accession number ZP_01157933) and Magnetospirillum magneticum AMB-1 (37% identity; RefSeq accession number YP_422072), a bacterium taxonomically related to the Azospirillum genera.
Deduced proteins encoded by ORF31, ORF32, ORF33, and ORF50 were revealed to be structural proteins commonly found in phage particles. ORF31 encodes a putative portal protein, the best identity being with a putative portal protein of Sinorhizobium medicae WSM419 (43% identity; GenBank accession number ZP_01412930). The deduced amino acid sequence of ORF32 displays strong similarity to phage prohead proteases that are involved in the processing of the prohead protein. ORF33 may encode the major component of the phage capsid, as it belongs to the family of head proteins of many bacteriophages, such as HK97, phi-105, and P27 (Interpro IPR 006444). The deduced protein of ORF50 shares similarity with phage tail proteins, but no conserved domain could be identified.
As for functions related to bacterial lysis, only one product encoded by ORF53 could be assigned to this process; indeed, this protein revealed similarities to members of the glycoside hydrolase family, i.e., lysozyme (Interpro IPR 002196). Finally, the fourth class of ORFs encoding putative phage proteins includes ORF9, ORF26, ORF28, and ORF30.
Phylogenetic analysis of functional genes.
Phage particle observations have revealed a morphology that is typical of members of the Siphoviridae family of bacteriophages. Information about the genus can be established by analyzing phylogenetic relationships of
Ab-Cd with other functional Caudovirales phages. It has been suggested that phylogenetic analysis should be performed at the level of individual genes in order to resolve viral taxonomical issues (41). Indeed, as bacteriophages are known for horizontal exchange of functional genetic elements drawn from a large shared pool (32, 55), any phage-specific protein could be obtained by a single modular interchange during evolution. Therefore, four proteins of
Ab-Cd, including terminase (ORF29), phage major capsid protein (ORF33), lysozyme (ORF53), and integrase (ORF72), which are commonly present in tailed phages, were chosen to build phylogenetic trees. Homologous protein sequences other than that of the phage major capsid protein have been identified from the draft of the genome of the Sp245 strain and are likely to be components of
Ab-Sp245, as corresponding genes lie in a region containing phage ORFs whose length matches the estimated size of
Ab-Sp245 observed here by PFGE. These sequences were compared with those of functional Caudovirales phages deposited in GenBank in order to construct the phylogenetic trees shown in Fig. 6.
Ab-Cd and
Ab-Sp245 lie on the same branch in the tree based on integrase, lysozyme, and terminase data but are not as closely related as other phages isolated from the same bacterial species, such as Salmonella or Staphylococcus phages, are related to each other. All trees show that
Ab-Cd and
Ab-Sp245 are deeply rooted, suggesting that they are quite different from previously described phages. The low bootstrap values obtained for the deep branches indicate that these phage proteins have diverged to an extent that prevents determining their evolutionary relatedness and therefore prevents classification of Azospirillum brasilense phages in known genera.
![]() View larger version (43K): [in a new window] |
FIG. 6. Unrooted neighbor-joining trees based on the aligned amino acid sequences encoded by four phage genes. Only genome data obtained from virions were used. (A) Terminase; (B) major capsid protein; (C) lysozyme; (D) integrase from Ab-Cd, from Ab-Sp245 (no protein homologous to a major phage protein could be identified in the A. brasilense Sp245 draft genome), and from 20 other phages. Amino acid sequences from A. brasilense phages are highlighted in bold. A Blosum 30 matrix was calculated, with a gap penalty of 10.0. The numbers at the nodes are bootstrap values based on 1,000 resamplings. Bootstrap values less than 50 are not shown. Scale bar, 0.2 amino acid substitutions per residue. M, Myoviridae; S, Siphoviridae; P, Podoviridae; C, unclassified Caudovirales.
|
Ab-Cd prophage.
Ab-Cd and pRhico sequences have revealed that the two entities encompassed regions that were 43.8% (27.3 kb) identical. Given the relatedness between strain A. brasilense Cd and strain A. brasilense Sp7 (strain Cd being isolated as a scarlet variant of strain Sp7 after inoculation on a plant) (25, 59) and the homology between the two phages isolated from those strains (see above), the same genetic organization could be expected in both strains. This implied either that
Ab-Cd was integrated as a prophage in plasmid pRhico or that a region of the prophage was duplicated in the A. brasilense Cd genome. To test this hypothesis, PCR amplifications of ORF40, ORF50, and ORF53 of
Ab-Cd were carried out with DNA isolated from
Ab-Cd,
Ab-Sp7,
Ab-Wb1, and
Ab-Sp245. PCR fragments of the expected size were obtained from DNA isolated from
Ab-Cd and from
Ab-Sp7, sustaining the hypothesis that the
Ab-Cd and
Ab-Sp7 genomes are strongly related (data not shown). No amplification products could be obtained from DNA isolated from
Ab-Wb1 and
Ab-Sp245, suggesting that these ORFs are absent from these phages or have divergent sequences. To determine on which replicons
Ab-Cd and
Ab-Sp7 were integrated, a separation of plasmids was performed (see Materials and Methods) followed by Southern hybridization experiments. First, identical plasmid patterns were observed for strains A. brasilense Cd and Sp7 that coincided with the pattern previously described, with sizes of replicons estimated at 130, 570, 630, 1,700, and >1,800 kb (Fig. 7A) (14). Second, a hybridization experiment performed with exoC, an ORF specific to pRhico (62), as a probe, enabled us to identify pRhico as being the 130-kb plasmid of both the Sp7 and Cd strains (data not shown). Third, ORF50 of
Ab-Cd, when used as a probe, clearly hybridized with the 570-kb plasmid of both Sp7 and Cd strains, but no signal was obtained with the pRhico plasmid (Fig. 7B). These observations indicate that
Ab-Cd and
Ab-Sp7 are not integrated as prophages in pRhico but are present in a larger replicon of 570 kb. These results also imply that a region of 27.3 kb spanning contigs 3 and 4 of pRhico (Fig. 5) is duplicated in the genomes of the two strains and that this duplication is located on the pRhico plasmid.
![]() View larger version (76K): [in a new window] |
FIG. 7. Genome localization of Ab-Cd. Plasmids of strains A. brasilense Cd (lane 1) and Sp7 (lane 2) were separated (A) and hybridized with ORF50 of Ab-Cd (B). The pRhico plasmid is indicated by an asterisk.
|
|
|
|---|
TEM analysis revealed two types of bacteriophages displaying ultrastructural features that are typical of members of the Siphoviridae family. Identical features were previously described for the bacteriophage of A. brasilense Sp7 (30). The two types can be distinguished not only by the sizes of particles but also by the sizes of encapsidated DNA and by the durations of the eclipse period upon induction by mitomycin C. Moreover, each type of bacteriophages seems to be associated with a specific bacterial species, as only "big" phages were evidenced for A. brasilense and "small" phages for A. lipoferum strains. This finding could be related to the usually narrow host range of a given bacteriophage, limited to a bacterial species or strain (18, 30).
Relatedness among Azospirillum phage genomes.
Hybridization experiments suggested DNA homology between the phage genomes of
Ab-Cd,
Ab-Sp7, and
Ab-Wb1. Restriction profiles of
Ab-Cd and
Ab-Sp7 are identical, and, given the relatedness between those two strains, one can speculate that these two phages are identical. The homology between
Ab-Cd and
Ab-Wb1 is less striking, as their restriction profiles display only a few common bands. Genomes of members of the Siphoviridae family are usually organized into functional modules (i.e., multigenic interchangeable elements whose genes are involved in the same process, such as the lysogeny module, tail gene module, etc.), and such modules appear to be conserved among prophages hosted by the same species, as demonstrated by genome analysis of six phages isolated from Streptococcus thermophilus (12, 15). However, the organization of modules can differ from one phage genome to another due to the accumulation of mutations and recombination events, leading to mosaic genomes. Such an organization might be found in the A. brasilense phages, explaining the differences in restriction profiles despite the presence of strong DNA hybridization signals.
Strain A. brasilense Sp245 seems to host a prophage that is not genetically linked to
Ab-Cd, as shown by the absence of hybridization signal. Preliminary analysis of the draft Sp245 genome (available at http://genomics.ornl.gov/research/azo) is in accordance with this finding, as some putative phage ORFs display low or no similarity with the
Ab-Cd sequence obtained in this study. Comparative genomic analysis of 18 Pseudomonas aeruginosa bacteriophages and 27 Staphylococcus aureus bacteriophages showed that some prophages hosted by strains of the same species are highly homologous whereas others display very little homology (38, 39). Therefore, phage diversity within a bacterial species is not restricted to the A. brasilense species. This diversity is even more obvious at the genus level, as prophages hosted by A. lipoferum strains do not share any common sequences with prophages hosted by the A. brasilense species. In contrast, the phage hosted by the A. doebereinerae strain might share common sequences with
Al-B518, and if so, this can probably be attributed to the genetic proximity between those two species (22).
One of the most striking features of the A. lipoferum phages induced in this study is the small size of encapsidated DNA (about 10 kb). Genome sizes of members of the Siphoviridae family usually range from 18 kb to 500 kb (17); indeed, the smallest genome of a Siphoviridae phage described in NCBI is 14,957 pb (phage bIL310 of Lactococcus lactis IL-1403; GenBank accession number AF323671). Attempts to establish the restriction profiles of A. lipoferum phage DNAs by use of various enzymes always resulted in smears; moreover, attempts to clone and sequence those DNAs failed to give data that could be assembled. All these observations seem to imply the presence of prophage-like gene transfer agents rather than of real prophages. In contrast to functional prophages, gene transfer agents typically package bacterial genome fragments and atypically package a portion of their own genome and constitute conspicuous mechanisms of generalized transduction (56). Bioinformatic analysis showed that gene transfer agents are widely distributed within alphaproteobacteria, especially in members of the Rhodobacterales family (40). However, this hypothesis needs to be experimentally validated.
Phage genome sequence.
The complete genome sequence of
Ab-Cd obtained in this study constitutes the first description of an Azospirillum bacteriophage. The process of annotation of the
Ab-Cd genome revealed many ORFs (81 out of 95) encoding either hypothetical proteins or conserved proteins with unknown functions; this result might be attributable to the low number of alphaproteobacteria bacteriophages actually described. Indeed, only three sequences of phages were obtained through complete genome sequencing, i.e., phage
JL001 of a marine alphaproteobacterium, phage SIO1 of Roseobacter SIO67, and phage PBC5 of Sinorhizobium meliloti (42, 50). Of the 14 ORFs of
Ab-Cd to which a function could be attributed, 11 display homology with phage ORFs found in alphaproteobacterium whole genomes but none was homologous to ORFs specific to phages of other bacterial classes. This could be due to high rates of sequence divergence and to the narrow host range of bacteriophages. Interestingly, a function can be attributed to only about 25% of the ORFs even when several sequences of prophages hosted by the same species are available, as illustrated with prophages of S. aureus or Mycobacterium tuberculosis (38, 39, 48). Thus, phages can be considered to be a rich reservoir for new proteins. Some of the ORFs encoding hypothetical proteins could be genes of lysogenic conversion, contributing to increase bacterial fitness; such a case is observed for many bacterial pathogens whose virulence factors are phage encoded (13, 16). However, some of the ORFs identified here encode proteins with functions that are essential in the lytic phage cycle, such as proteins for particle formation (structural proteins, terminase), proteins for bacterial lysis (lysozyme), and proteins for integration/excision (integrase). Using such proteins to resolve the taxonomy of
Ab-Cd proved to be unsuccessful, a result that emphasizes the problems associated with mosaicism of viral genomes when using sequence data for taxonomical purposes (41). The identification of functional modules remained difficult due to the low number of ORFs with an attributed function, but ORF26 to ORF33 could constitute a conserved genetic module encoding structural and maturation proteins.
Phages and dynamics of bacterial genomes and bacterial community.
Ab-Cd is found integrated on a 570-kb plasmid; whether this plasmid is conjugative or not is not known, but this could imply a putative horizontal transfer of phage DNA by conjugation. A 27.3-kb fragment of
Ab-Cd was also recovered within another replicon, the pRhico plasmid, indicating that genomic rearrangements occurred between plasmids and phages, leading to a duplication event. Major genomic rearrangements were previously characterized for several Azospirillum strains undergoing phenotypic switching, demonstrating genome plasticity among members of the Azospirillum genus (63). Bacteriophages could thus account for the dynamics and evolution of Azospirillum genomes and could be major contributors to the genetic individuality of Azospirillum strains, as described previously for other species (15, 37, 47).
Phage particle formation was induced here in vitro by the SOS-inducing agent mitomycin C; it would be of interest to study phage-Azospirillum interactions within the rhizosphere, the natural habitat of Azospirillum spp. Indeed, environmental abiotic factors, such as nutrients and temperature, were shown to induce production of bacteriophages in A. brasilense Sp7 (30). Eukaryote-bacterium interactions may also influence the cycle of phages; for example, a lysogenic avian E. coli pathogen showed an upregulation of phage gene expression when injected into a chicken (21) and phage genes were shown to respond to molecular cues emitted by the mammalian host of the pathogen S. pyogenes (10). Moreover, upregulation of phage-derived genes has been observed during promoter trap studies of Ralstonia solanacearum grown in the presence of tomato plants and of A. brasilense induced by wheat extracts (11, 49). The fact that bacteriophages can be induced in natural conditions has a considerable impact on bacterial densities; such an impact has been demonstrated for viruses in aquatic environments (57, 65, 66) and by artificial release of a Serratia liquefaciens lysogen into soil samples (3). Despite a previous study showing no effect on the declines in the viability of strains A. brasilense Sp7 and Cd in soil inoculated with the Sp7 phage (29), the potential role of bacteriophages in population dynamics of Azospirillum spp. and of rhizospheric bacterial communities cannot be excluded. Also, the presence of prophages might provide immunity to infection by any other bacterial virus. Finally, this work might help optimize the use of bacteriophages as a tool for Azospirillum genetic studies.
This work was supported by a fellowship from the French Ministère de l'Education Nationale de la Recherche et des Nouvelles Technologies awarded to M.B.
Published ahead of print on 7 December 2007. ![]()
|
|
|---|
JL001: insights into its interaction with a sponge-associated alpha-proteobacterium. Appl. Environ. Microbiol. 71:1598-1609.
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»