Previous Article | Next Article ![]()
Journal of Clinical Microbiology, January 2005, p. 66-73, Vol. 43, No. 1
0095-1137/05/$08.00+0 doi:10.1128/JCM.43.1.66-73.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Department of Molecular Microbiology and Biotechnology, the George S. Wise Faculty of Life Sciences, Tel-Aviv University, Tel-Aviv, Israel
Received 4 August 2004/ Returned for modification 14 September 2004/ Accepted 21 September 2004
|
|
|---|
|
|
|---|
We have been studying the virulence factors of ExPEC strains involved in human disease as well as in avian colisepticemia, such as serogroups O2 and O78, which constitute two of the major serotypes involved in avian sepsis. The diseases caused by the two serotypes in poultry are symptomatically similar: they start as a respiratory tract infection, which leads to a systemic infection involving the internal organs and sepsis (7, 18, 22). However, several differences between the two serogroups have already been identified, and different virulence factors were found in strains belonging to each serogroup (17). A thorough genetic comparison of these strains could reveal the common mechanisms involved in sepsis in general and in avian sepsis in particular. Human strains of serotype O2 are also associated with human UTIs, while strains of serotype O78 have been known to cause human NBM. Therefore, a comparison of the two serogroups may also reveal genetic factors relevant to the study of human disease.
The molecular procedure of suppression subtractive hybridization (SSH) allows the comparison of two genomes and the detection of specific sequences that are present exclusively in one of them. This technique has been shown to be extremely efficient in detecting genomic sequences and genomic islands that are strain specific (2, 26). In research on pathogenic bacteria, SSH was used to identify virulence factors by comparing pathogenic and nonpathogenic bacterial strains of the same species (15, 29).
In this work, the genome of nonpathogenic laboratory strain K-12 was subtracted from the genome of avian pathogenic strain O2-1772 and, in parallel, from the genome of avian strain O78-9. The results from these two experiments were combined and analyzed. A total of 158 unique sequences that are not present in nonpathogenic strain K-12 were identified. Interestingly, only a few pathogen-specific sequences were shared by the two strains, even though they cause clinically similar diseases.
|
|
|---|
DNA techniques. Isolation of genomic DNA was performed as described previously (23). DNA fragments obtained from the subtractive hybridization were cloned into pGEM-T plasmids of the Promega pGEM-T and pGEM-T EASY Vector Systems kit, according to the instructions of the manufacturer.
Sequence analysis. Sequence analysis was carried out by using National Center for Biotechnology Information algorithms BLASTN, BLASTX, and BLASTP (3).
SSH. The SSH technique was carried out as described previously (8). We used the CLONTECH PCR-Select Bacterial Genome Subtraction kit, as recommended by the manufacturer. Genomic DNAs isolated from strains O2-1772 and O78-9 were used as tester DNAs for two separate hybridization experiments with the driver DNA of strain K-12. Two rounds of PCR were performed with the hybridization products to amplify the specific fragments unique to each of the pathogenic strains.
Confirmation of SSH libraries. Genes found only in one of the SSH libraries were further examined. PCR was performed with the sequences which were found to be virulence related or sequences for which we could not find a known function or homology. PCRs were carried out under nonstringent conditions.
|
|
|---|
Sequence analysis of fragments identified by SSH. A search for sequences homologous to the specific sequences was performed with the BLASTN and BLASTX algorithms (3) in order to further analyze them. The results are summarized in Table 1. The sequences were divided into several categories, as shown in Fig. 1. The first category included nine sequences of putative or verified virulence factors, such as components of iron uptake systems, genes of the type III secretion system (TTSS), and the specific pili of each strain. The second category included sequences classified as potential virulence-associated sequences, according to their predicted conserved domains. These were sequence O2-38, which encodes a predicted ATP-binding protein involved in virulence, and sequence O78-95, which encodes a putative adhesin. The third category included 21 sequences with known or putative functions whose role in pathogenicity is not clear. The fourth category included 34 sequences with no known function that could be divided into two types: fragments with similarity to sequences of unknown function (e.g., conserved hypothetical proteins) and fragments for which no similarity was found.
|
View this table: [in a new window] |
TABLE 1. Specific SSH fragments and their homologies determined with the BLASTN and BLASTX algorithmsa
|
![]() View larger version (20K): [in a new window] |
FIG. 1. Breakdown of the specific sequences of each library into categories according to their homologies.
|
Known PAIs from other E. coli strains. Analysis of the two SSH libraries revealed that many of the specific sequences were homologous to sequences from pathogenicity islands (PAIs) of other pathogenic bacteria, which were presumably acquired by horizontal transfer (13). Interestingly, several sequences of the SSH libraries of the septicemic strains showed similarity to the sequences of E. coli strains involved in UTIs (9) and enterohemorrhagic strain E. coli O157:H7 (20). The data also indicate PAIs which are homologous to those of other bacteria, such as Shigella flexneri (the SRL, SHI-3, and SHE PAIs).
Verification of SSH libraries. Sequences that were present in the subtracted library of only one of the strains were further tested by PCR to verify their absence from the genome of the other strain. Forty-five virulence-related sequences and sequences with unknown function were examined. In general, the PCR results were in agreement with the SSH results, showing that 34 sequences were absent from one of the strains, whereas only 9 sequences not observed in the SSH libraries could be detected by PCR. Although we assume that the SSH library does not contain all the strain-specific sequences, our library appears to represent many large gene clusters which encode complex functions, such as PAIs, as well as single genes.
Comparison of the two SSH libraries. The two SSH libraries verified by PCR have several common characteristics, as can be observed from sequence analysis (Table 1). Most of the sequences showed homology to pathogenic strains of E. coli, as well as to other pathogenic bacteria, most frequently to the phylogenetically closely related species S. flexneri and different serovars of Salmonella. Both libraries include numerous mobility-associated sequences, such as transposases and integrases, along with phage-related sequences, plasmid sequences, and insertion sequence-associated sequences. This implies that both libraries are highly specific for areas of the genome that were acquired horizontally, some of which may be related to virulence.
Components of two major virulence mechanisms were found in both libraries. The first is iron uptake systems: both pathogenic strains were found to possess different iron acquisition systems that were absent from the nonpathogenic strain. Strain O78-9 was found to possess the yersiniabactin system, coded by the Yersinia high-pathogenicity island. Strain O2-1772 has components of the sitABCD system, known from S. flexneri and Salmonella enterica serovar Typhimurium. The second major virulence mechanism is adhesins: both libraries contained components of virulence-related pili. Two minor subunits of the AC/I pilus were found in strain O78-9, and a chaperone involved in the assembly of the P pilus was found in strain O2-1772.
The two SSH libraries were also compared to each other by use of the BLAST algorithm. First, each sequence of the O78-9 library was run by using the BLASTN algorithm. Usually, the sequence showed similarity to a part of an open reading frame (ORF). We selected the complete ORF of the closest homolog and performed a search against the entire library of O2-1772 with the BLAST algorithm. Thus, even if two sequences from the two libraries were homologous to different parts of the same gene, the similarity between them could be determined. Some of the sequences did not show homology within an ORF, and these were run without manipulations against the whole library of O2-1772. We found two adjacent genes of TTSS: eivC of O78-9 and eivA of O2-1772. Both genes had very strong homology to the TTSS of E. coli O157:H7. Sequences O78-151 and O2-37 showed only limited similarity in one domain by use of the BLASTX algorithm and were clearly very distant.
Distribution of specific sequences in other septicemic O2 and O78 strains. The comparison of two strains, one of serotype O2 and one of serotype O78, indicated a high degree of variability in virulence-related sequences between the strains. In order to determine if these sequences were conserved within the same serotype, we performed PCRs to look for the presence of the specific sequences in additional strains of the same serotype. Altogether, six additional strains of each serotype were examined. All the strains were invasive, and about half were of human origin. The results (Fig. 2 and Table 2) indicate that among the 45 sequences examined, only 1 sequence (the putative autotransporter, sequence O78-95) was serogroup specific and appeared in all of the O78 strains but none of the O2 strains. Most other sequences were patchily distributed and were present in only a small subset of strains; therefore, they may have been acquired by fairly recent horizontal gene transfer. Genes such as O2-207, O2-345, and O78-96 probably represent an ancient transfer event, as they were found in all the pathogenic strains tested but not in strain K-12.
![]() View larger version (39K): [in a new window] |
FIG. 2. Results of PCRs of specific sequences, performed with 14 septicemic strains and nonpathogenic strain K-12. Lanes 1 to 16, 1-kb DNA marker, O2-1772, O2-YN, O2-MAN, O2-SA, O2-U33, O2-B18, O2-BEN, O78-9, O78-285, O78-286, O78-287, O78-6, O78-7, O78-63-1, and K-12, respectively. (a) Multiplex PCR with primers specific for the sequences of O2-210, O2-334, O78-18, and O78-55; (b) PCR with primers specific for the putative autotransporter (O78-95 sequence).
|
|
View this table: [in a new window] |
TABLE 2. Distribution of specific sequences in 14 septicemic E. coli strains and nonpathogenic strain K-12a
|
|
|
|---|
We used SSH in order to identify specific sequences of two different serotypes of septicemic E. coli, one of serotype O2 and one of serotype O78. The SSH libraries obtained consisted of 158 sequences, including putative new virulence factors and sequences with homology to PAIs not known to be present in septicemic strains. Although these libraries are substantial, we assume that they do not contain all the unique genes absent from strain K-12. For example, we know of the existence of the yersiniabactin iron uptake system in strain O2-1772 (11), but it was not detected in our SSH library. However, the results of the subsequent PCR experiments indicated that the SSH libraries were reliable: most of the genes absent from the SSH library were also not detected by PCR.
Avian strains of serogroup O2 involved in cellulitis were previously analyzed by SSH (24). However, this study contained very few genes of significance for pathogenesis. Thus, only 62 fragments were obtained in the study of the cellulitis strains, half of which were phage sequences and only 1 of which was similar to PAI genes. Our experiments detected the presence of many more unique sequences (n = 102), many of which were similar to PAIs of various virulent strains. The differences in the results probably reflect the fact that the cellulitis strains are less virulent than the septicemic O2 strain tested in our study. Thus, the more virulent septicemic strains appear to have acquired a larger array of virulence factors.
A comparison of the SSH libraries of strain O2-1772 and strain O78-9 revealed that the diversity between the two strains is very high. These results are consistent with the results from a multilocus sequence type (MLST) analysis showing that the evolutionary distance between these strains is indeed large, as reflected by the sequence divergence of six housekeeping genes with adequate polymorphisms. According to the MLST results, the various O2 and O78 strains (22 in total) are distributed on different branches of the phylogenetic tree (1). This study indicated that clonal division in E. coli O78 and O2 strains is host independent, and strains from the same host show high degrees of variability. It therefore follows that the ability of strains to conquer new niches and infect a new host species depends more on horizontally acquired DNA than on the vertically inherited genotype.
Many of the sequences examined in this study are of genes that have an unknown function or that have no known homologs. Twenty-seven such sequences were present only in the O2-1772 or O78-9 sequence. PCRs performed with additional strains in order to determine if these sequences were conserved within strains of the same serotype indicated that only one sequence was serogroup specific and appeared in all the O78 strains but none of the O2 strains. This sequence (O78-95) may therefore serve as a diagnostic marker for septicemic O78 strains. This is clearly one of the spin-offs of this study that is being investigated further.
The prevalence of sequences associated with genomic plasticity in both libraries supports the assumption that the pathogenic strains evolved by processes involving genome remodeling and the horizontal acquisition of genomic regions from other pathogenic bacteria. In addition to the genes associated with genomic plasticity, many phage sequences were found in both libraries. However, these sequences differ between the two strains, with the exception of lambdoid phages, implying that the two bacterial strains have different susceptibilities to different families of phages. The differences in phage susceptibilities probably reflect differences in restriction modification systems and surface structures, such as pili or capsules, and consequently affect the genetic information that each strain can acquire from these phages.
The common virulence factors, such as the TTSS and the iron uptake systems, could be important to the infection process, as they exist in both strains. The fact that not many virulence factors are common to both strains is unexpected, since phenotypically both strains cause the same disease. Our findings imply that the two strains are using different factors with similar roles in the various stages of the infection process. This possibility is illustrated in Table 3, which compares the virulence factors that are thought to participate in the different stages of infection with strains O2-1772 and O78-9. The two strains have adherence pili, but these are of different types: O2 codes for the P pilus (this work), while O78 codes for the AC/I pilus (6) and for a nonfimbrial putative adhesin (this work). Another virulence factor that varies between the strains is the capsule. While O2 possesses a K1 polysaccharide capsule, O78 lacks such a capsule but is enveloped in a loose layer of extracellular material which is chemically similar to the O antigen (unpublished results).
|
View this table: [in a new window] |
TABLE 3. Virulence factors of strains O2-1772 and O78-9 in the different stages of infection
|
The availability of sequenced genomes and of technologies such as SSH makes it possible to examine the virulence factors through comparative genomics. Previous whole-genome comparisons of E. coli strains have shown great differences in the gene contents of strains. However, none of the comparative genomics of E. coli strains that have been sequenced compared strains that cause the same disease and that target the same host tissues. The only exception is in the two genome sequences of O157:H7. In these strains, the virulence factors and the genes coding for them show a very high degree of similarity. These results are in contrast to our findings for the septicemic O2 and O78 serogroups, which showed a very high level of genome plasticity. An additional and important conclusion from these data is the indication that the large pool of variable virulence genes is accessible to septicemic bacteria, independent of the host. It therefore remains to be determined whether host specificity exists among the septicemic ExPEC strains and which factors control the zoonotic risk.
Comparison of strains causing similar types of diseases is essential for the identification of shared traits which are important for pathogenicity and for study of their conservation, transfer, and evolution.
This work was supported by the Manja and Morris Leigh Chair for Biophysics and Biotechnology, the Israel Center for Emerging Diseases, and European Community project COLIRISK. D.M. and U.G. thank the Constantiner Institute for Molecular Genetics.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»