Department of Biological Sciences, Faculty of Science, National University of Singapore, Science Drive 4, Singapore 117543, Republic of Singapore,1 Department of Immunology, University of Washington, Seattle, Washington 98195-7650,2 Department of Microbiology, University of Barcelona, Barcelona, Spain 08071,3 Department of Microbiology and Immunology, Faculty of Medicine, University of Saskatchewan, Saskatoon, Canada S7N 5E54
Received 29 November 2004/ Accepted 14 February 2005
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The pathogenesis of A. hydrophila is multifactorial. Most studies to date have concentrated on the characterization of a few virulence factors in different animal models by using different strains (10, 34, 42), making it very difficult to evaluate the significance of each gene in the virulence of A. hydrophila. This study aims to identify more putative virulence determinants and characterize them in an integrated manner.
Suppressive subtraction hybridization or genomic subtraction offers a genome-level approach to identifying the genetic differences between virulent and avirulent strains of bacteria (22, 47). In this study, we used genomic subtraction to identify 22 unique DNA fragments encoding 19 putative virulence factors and seven new open reading frames (ORFs), which are frequently present in a group of virulent strains of A. hydrophila. In addition, four genomic islands (GIs), which differ from the rest of the genome in G+C content and which carry mobility-associated genes (integrases or transposes) and putative virulence genes, were found. These are genes expressing the O-antigen and capsule (48, 49), a phage-associated gene cluster, and a type III secretion system (TTSS) gene cluster. Subsequently, these putative virulence-associated genes were located on a preliminary physical map of A. hydrophila PPD134/91 to determine their genetic organization in A. hydrophila, which will provide insight into the molecular basis of pathogenicity of this bacterium.
We further studied the individual or combined roles of these putative virulence genes using a consistent animal model and in the same strains. Knockout mutations were constructed to elucidate the contributions of a number of these virulence genes to A. hydrophila pathogenesis. Of the single-knockout mutations examined, only that in the TTSS (
ascN) significantly affected the 50% lethal doses (LD50s) in a blue gourami infection model. In addition, a triple deletion mutant of S-layer, serine protease and metalloprotease (
ahsA
serA
mepA) increased the LD50 to the same level as the ascN mutant.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
DNA sequencing and sequence analysis.
DNA sequencing was carried out on an Applied Biosystems PRISM 3100 genetic analyzer with an ABI PRISM BigDye Terminator Cycle Sequencing kit (Applied Biosystems). The sequences were edited using the manufacturer's software. Sequence assembly and further editing were carried out with Vector NTI DNA analysis software (InforMax). BLASTN, BLASTP, and BLASTX sequence homology analyses and a protein conserved-domain database analysis (CDD Search) were performed using the BLAST network server of the National Center for Biotechnology Information.
Genome walking and genomic subtraction.
Advantage Polymerase 2 (Clontech) was used for genome walking. Genome Walker libraries were constructed by using five restriction enzymes (DraI, EcoRV, PvuII, ScaI, and StuI). The cycling parameters for genome walking were as follows: 7 cycles of 15 s at 94°C, 4 min at 72°C, 32 cycles of 15 s at 94°C, and 3 min at 67°C. The amplified fragments were cloned into the pGEM-T Easy Vector (Promega), transformed into E. coli JM109-competent cells, and sequenced. Bacterial genome subtraction was performed as described previously by Zhang and coworkers (47).
Construction of defined insertion mutants and deletion mutants.
To obtain defined insertion mutants, oligonucleotides containing flanking XbaI restriction sites (primers sequences available upon request) were used to amplify internal fragments from the targeted genes. The amplified fragment was ligated to pGEM-T Easy vector and transformed into E. coli JM109. The internal fragment was recovered by XbaI restriction digestion and ligated to XbaI-digested and dephosphorylated suicide vector pRE112 (11). The recombinant plasmid was transformed into E. coli MC1061(
pir), selecting for Cmr to isolate the pRE-112-derived plasmids. These Plasmids were transformed into E. coli S17-1(
pir) and transferred by conjugation to A. hydrophila AH-1 (Apr), selecting for Cmr and Apr colonies. The insertion of plasmids into the chromosomes of these mutants was confirmed by PCR using appropriate primers.
Nonpolar deletion mutants were constructed according to the method for suicide plasmid pRE112 (11). Briefly, the targeted gene and at least 300 bp of flanking sequences were amplified with oligonucleotides containing XbaI restriction sites and then cloned into pGEM-T easy vector, followed by inverse PCR using these constructs as templates. The inverse PCR products were purified and self ligated to get the deleted constructs of targeted genes. The deletion constructs were digested with XbaI and ligated to XbaI-digested and dephosphorylated pRE112. This construct was transformed into E. coli S17-1(
pir). The single-crossover mutants were obtained by conjugal transfer into A. hydrophila AH-1. Double-crossover mutants were obtained by selecting against the presence of the sacB gene carried by the vector. Mutants resistant to sucrose were isolated by being plated onto LB-sucrose agar (1% tryptone, 0.5% yeast tract, 1.5% agar, 12% sucrose). The double-crossover mutants were confirmed by PCR.
LD50 studies of fish.
Healthy blue gourami (Trichogaster trichopterus Pallas) were obtained from a commercial fish farm, maintained in well-aerated dechlorinated water at 25 ± 2°C, and acclimatized to laboratory conditions for at least 15 days. Fish were approximately 13 g each and were about 3 months old. Three groups of 10 fish each were injected intramuscularly with 0.1 ml of phosphate-buffered saline-washed bacterial cells adjusted to the required concentrations. LD50 studies were carried out in triplicate for all of the strains. Fish were monitored for mortality for 7 days, and LD50s were calculated by the method of Reed and Muench (31).
Pulse-field gel electrophoresis. Pulse-field gel electrophoresis was carried out as previously described (46). ProMega-Markers Lambda Ladders (Promega) were used as the size markers. Band sizes were calculated from at least three independent separations of restriction fragments.
Statistical analysis.
LD50 data were analyzed using Student's t test. P values of <0.05 were considered significant.
Nucleotide sequence accession numbers.
Thirteen genes of A. hydrophila AH-1 or AH-3 for construction of knockouts were assigned the following accession numbers: aerA (AY442276), bvgA and bvgS (AY841798), f85 (AY442275), flhA1 (AY841793), flhA3 (AY841794), hlyA (AY442273), hup1 (AY442272), hup3 (AY841797), mepA (AY841796), ompAI (AY442271), opdA (AY442274), and serA (AY841795).
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
4) but absent in most of the six avirulent strains (
2) (Table 3). Altogether, 22 unique DNA fragments (F97 and PB45 are in the same DNA fragment) representing about 200 subtracted clones were obtained from two rounds of genomic subtraction and were commonly present in the eight virulent strains (
4) (Table 2) (47). These 22 DNA fragments were assumed to encode potential virulence genes and were subjected to further analysis.
|
|
The other nine genes identified are new putative virulence factors for A. hydrophila and have high similarities to known virulence proteins of other pathogens, including the homologues of a histone-like protein (HU and F52), oligopeptidase A (OpdA and F88/F109/F61/F72), outer membrane protein (OmpAI and F2/F3), VsdC and a type III secretion protein homologue (F11), acinetobactin biosynthesis protein (F20), a two-component BvgA and BvgS system, and arylsulfotransferase.
VsdC is essential for virulence in Salmonella enterica serovar Dublin (18). The VsdC homologue in A. hydrophila contains a VIP2 domain, which belongs to a family of actin-ADP-ribosylating toxins. The presence of a type III secretion protein homologue in A. hydrophila PPD134/91 correlates well with the identification of a TTSS in this strain, as discussed further below. The acinetobactin biosynthesis protein F20 is involved in iron acquisition and may have an important role in the pathogenesis of Acinetobacter baumannii (25). In Bordetella pertussis, the bvgA and bvgS two-component system controls the expression of many virulence genes, including those encoding pili (26) and adenylate cyclase (12). These bvgA and bvgS homologues may thus also be involved in the regulation of pathogenesis in A. hydrophila PPD134/91. On the other hand, arylsulfotransferase has been suggested to be a virulence factor of the Campylobacter-Wolinella family of many oral bacteria and is widely distributed in Campylobacter species (43, 45).
The rest of the eight proteins showed high homology to proteins in other bacteria that have not yet been linked to virulence. They are topoisomerase (PB28/PB35), GGDEF family protein (PB60), a YaiI/YqxD family protein (F89), para-aminobenzoate synthase (F92), a putative exported protein (F93), a transporter transmembrane protein (F99/F106), a transmembrane protein (PB28/PB35), and a reverse transcriptase-like protein (F97). The functions of these proteins in A. hydrophila pathogenesis remain to be clarified.
For the last seven DNA fragments (F32, F34, F58, F85, F87, PA1/PA4/PB62, and PA6/PA98/PB38/PB78), the ORFs encoded had no significant matches with entries in the GenBank database and may represent novel virulence determinants. These presumptive virulence factors will be studied in future experiments.
Identification of a phage-associated genomic island.
The second approach adopted to search for putative virulence genes was to look for pathogenicity island (PAI) markers. PAIs are specific genomic islands with a large unstable chromosomal region that encodes virulence genes, and they are present in most of the pathogenic bacteria (13). Recently, ssrA, a small stable RNA molecule (tmRNA) has been reported to reside in or near the junction point of a mosaic of Salmonella-specific sequences (8) and serves as the insertion site for acquired sequences such as the cryptic phage CP4-57 in E. coli and PAIs in Vibrio cholerae (17) and Dichelobacter nodusus (3).
Primers ssrA-F (5-CAAACGACGAAAACTACGC-3) and ssrA-R (5-GGTACTACATGCTTAGTC) within the conserved ssrA region were designed, and an ssrA gene was identified in A. hydrophila PPD134/91 via a PCR using these primers. Subsequently, a series of genome walking experiments led to the identification of a 23-kb DNA region. Flanked by ssrA at the left end, this region exhibits significantly lower G+C content (49.7%) than the average genome G+C content of A. hydrophila (57 to 63%), which is a common characteristic of GIs (13). It also shows a mosaic distribution of G+C content (Table 4; Fig. 1). Furthermore, two ORFs (ORF1 and ORF2) (Table 4) of this region show high similarities to the insertion element IS1650 and are adjacent to ssrA. The G+C content of thirteen ORFs (Table 4) did not differ significantly from that of A. hydrophila, indicating that they might have been acquired from species with G+C content similar to that of A. hydrophila or that the base composition of these acquired DNA has gradually adapted to the host genome (19).
|
|
The distribution of this putative GI among 14 A. hydrophila strains was surveyed by Southern analysis using probes from four regions of the DNA (Fig. 1; Table 4). The four probes hybridized with most of the eight virulent strains but only hybridized with one out of the six avirulent strains tested (Table 5). The hybridization signals were strong in virulent strains ATCC 7966, PPD134/91, and Xs91-4-1 and the avirulent strain L36 but weak in the other five virulent strains, suggesting some genetic variations among the virulent strains. The results therefore indicate that this putative GI is PPD134/91 specific and/or virulent strain specific. Identification and sequencing of this GI in other virulent strains will help to clarify this issue.
|
Identification of a TTSS gene cluster.
TTSSs have been reported in many gram-negative animal and plant pathogens, including A. hydrophila (6, 9, 46). As reported previously (46), an ascV homologue is present in A. hydrophila PPD134/91. Furthermore, an ascU homologue, near one end of the TTSS in A. hydrophila AH-1, was also identified in A. hydrophila PPD134/91 by PCR (ascU-F, 5-TGGTGATCGCCATCGCCGA-3; ascU-R, 5-GACGGCGCTTGCTCTTGAT-3). The identification of these two TTSS homologues strongly indicated that a TTSS cluster is also present in A. hydrophila PPD134/91. This TTSS cluster was located on the chromosome of A. hydrophila PPD134/91 (Fig. 2). The complete TTSS sequence of A. hydrophila PPD134/91 is not yet available. However, analysis of the TTSS gene cluster of A. hydrophila strain AH-1 revealed an ORF found near ascU showing high homology to the P4-family integrase of a variety of bacteria. Hueck (15) reported that the TTSS may be acquired as intact genetic blocks by horizontal gene transfer during evolution. It is thus reasonable to speculate that the integrase may have been involved in the original mobilization of the TTSS into the chromosome of A. hydrophila. Downstream sequencing of ascU in A. hydrophila PPD134/91 will facilitate the understanding of the evolutionary history of this TTSS. Moreover, complete TTSS gene clusters have been recently identified in the two A. hydrophila strains, AH-3 and SSU (accession no. AY763611) (38), which further suggests that a complete TTSS gene cluster may be also present in PPD134/91.
|
The 22 unique DNA fragments from the two rounds of genomic subtraction as well as the O-antigen and capsule clusters (48, 49), the phage-associated island, and the TTSS gene cluster have been positioned on the physical map by Southern hybridization (Fig. 2). The results showed O-antigen and capsule gene clusters were located at different regions of the chromosome (Fig. 2), while these two clusters were located near each other in E. coli (33). The locations of the TTSS cluster and the F11 fragment, which contains a type III secretion protein homologue, are quite far away from each other. This was not surprising, as the genes encoding a secretion apparatus are usually clustered, while the genes encoding the secreted proteins and their transcriptional regulators are often located in unlinked positions (15).
Ten subtracted clones consisting of seven unique DNA fragments were found to be located in the same PacI-digested fragment C14 (Fig. 2) of strain PPD134/91. However, there is as yet no indication that these virulence genes are located in a PAI. Most of these genes were not clustered in the same PacI-digested fragments in the genomes of seven other virulent strains by Southern analysis (data not shown).
Construction and characterization of mutants.
One strategy to examine whether the putative virulence genes identified in this study are involved in pathogenesis is to construct knockouts and test their virulence in the blue gourami fish model. Previous attempts to introduce plasmids or transposons into A. hydrophila PPD134/91 were not successful. A genetic barrier may exist in strain PPD134/91 and prevent the construction of knockouts in this strain. Therefore, other well-studied pathogenic strains of A. hydrophila such as AH-1 and AH-3 (24, 46) were used as the hosts for the construction of mutants. Fifteen genes in A. hydrophila AH-1 and/or AH-3, which are homologous to their corresponding genes in A. hydrophila PPD134/91, were successfully cloned for the construction of knockouts. In general, insertion mutants were constructed first to examine the effect of the target genes. If the LD50s of insertion mutants increased by at least 0.5 log, deletion mutants were then constructed in the same genes to confirm that the attenuation was not due to polar effects. We also included several deletion mutants such as
ompAI,
ascN, and
bvgA to confirm that there were no differences in the LD50s between insertion and deletion mutations. vsdC, the putative type III secretion protein (F11), and the phage-associated island sequences could not be successfully cloned from strains AH-1 and AH-3, and these genes were thus not characterized in this study. The construction of knockouts in these PPD134/91-specific genes is ongoing and will shed light on their roles in PPD134/91 pathogenesis.
The genes used for construction of mutants can be classified into three groups. Group I includes genes expressing hemolysin (hlyA) and aerolysin (aerA), which have been reported to be associated with the virulence of A. hydrophila (1, 42). Group II includes genes expressing histone-like protein (hup1 for AH-1 and hup3 for AH-3) (36); opdA, expressing ligopeptidase A (27); bvgA, bvgS, and ascN (46); and ompAI, expressing an outer membrane protein (40). These gene products are homologous to known virulence factors of other pathogens. Group III includes genes expressing a novel putative virulence protein (f85) and a recently identified polar flagellar assembly protein homologue, FlhA (flhA1 for AH-1 and flhA3 for AH-3), which appears to be involved in biofilm formation (H. B. Yu and K. Y. Leung, unpublished data).
Our results showed that the LD50s of all the mutants except ascN were comparable to that of the wild-type strains by using blue gourami fish as the infection model. The LD50s of the wild types AH-1 and AH-3 were 105.3 and 106.5, respectively. LD50s of both insertion and deletion mutants of ascN in AH-1 were determined to be 106.3. Wong and coworkers (42) reported that inactivation of hlyA or aerA alone showed no statistically significant attenuation in a suckling mouse model when compared to the wild type. Our virulence assay showed similar results for these genes, and extended this finding to most of the other virulence genes which we analyzed. In all the mutants tested, only ascN exhibited an LD50 1 log higher than that of the wild type-strain AH-1. AscN, a homologue of the YscN of Yersinia species (41), possibly interacts with membrane-bound components of the TTSS apparatus to energize secretion or to provide the energy for the assembly of the secretion apparatus. Disruption or deletion of ascN may thus render the TTSS nonfunctional, leading to the increase in LD50s and consistent with our previous study (46). The fact that most of the virulence genes did not attenuate the mutants suggested that the pathogenesis of A. hydrophila is multifactorial in nature. Disruption of more than one gene or a whole gene cluster as in the case of TTSS appears to be necessary to observe differences in the LD50s. On the other hand, the virulence of A. hydrophila may be also bacterial strain, infection route, and animal model dependent. It has been reported that a single mutation in cytotoxic enterotoxin of A. hydrophila SSU resulted in an
300-fold increase in LD50s by intraperitoneal injection in Swiss-Webster mice (44); a single mutation in elastase of A. hydrophila AG2 resulted in an
100-fold increase in LD50s by intraperitoneal challenge in rainbow trout (5).
In addition, double deletions mutations were made in strain AH-1, such as S-layer and metalloprotease (
ahsA
mepA; LD50 = 105.7), S-layer and serine protease (
ahsA
serA; LD50 = 105.7), and a triple-deletion mutant of S-layer, metalloprotease, and serine protease (
ahsA
serA
mepA; LD50 = 106.5). Different proteases have been shown to be involved in A. hydrophila virulence (5). The double-deletion mutants resulted in about half of a log increase in LD50s, but the attenuation was not statistically significant (P > 0.05). However, fish infected with these double mutants survived 1 to 2 days longer than those infected with the wild type at the same lethal dosage. Statistically significant attenuation was observed for the triple mutant, where the LD50 increased about 1 log. Our results therefore strongly support the idea that virulence factors in A. hydrophila pathogenesis work in a concerted manner and multiple factors are required to produce the observed deleterious effect.
Conclusion.
A variety of putative virulence genes in A. hydrophila have been identified by both genomic subtraction and GI analysis in this study. These include known A. hydrophila virulence genes (encoding hemolysin and aerolysin), as well as other genes showing homologies to known virulence factors, such as bvgA, bvgS, vsdC, and ompAI, which have not yet been examined with A. hydrophila. In addition, the putative virulence gene clusters, such as the presence of a phage-associated island and a TTSS in A. hydrophila PPD134/91, were established. Subsequent positioning of these putative genes and gene clusters on a physical map of strain PPD134/91 has helped us in better understanding the chromosome organization and the molecular mechanisms of pathogenicity in this bacterium.
This is the first report to present a comparative study of different virulence factors in A. hydrophila pathogenesis by constructing knockouts in the same strains and infecting the same infection host. Our virulence assay of these mutants demonstrated that, as is increasingly observed for other pathogens, virulence in A. hydrophila is complex and involves multiple virulence factors, which may work in concert. Construction of more multiple mutants, infection of different animal models and inclusion of other functional assays are under way to elucidate the mechanisms of these putative virulence factors. The putative virulence genes presented in this work will thus form the basis for further investigation of the pathogenesis of A. hydrophila and will be useful for data-mining for the development of effective vaccines, diagnostic and novel therapeutics against animal and human infection caused by motile aeromonads.
| ACKNOWLEDGMENTS |
|---|
We are grateful to Michael Janda from the California Department of Health Services for providing us with some of the A. hydrophila isolates.
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| J. Bacteriol. | Microbiol. Mol. Biol. Rev. | Eukaryot. Cell | All ASM Journals |
|---|