|
|
||||||||
1 Division of Medical Microbiology and Genitourinary Medicine, University of Liverpool, Liverpool L69 3GA, UK
2 Regional Adult Cystic Fibrosis Unit, Cardiothoracic Centre, Liverpool L14 3PE, UK
Correspondence
Craig Winstanley
C.Winstanley{at}liv.ac.uk
Received 9 December 2005
Accepted 20 February 2006
Abbreviations: CF, cystic fibrosis; FHA, filamentous haemagglutinin; LES, Liverpool epidemic strain; SSH, suppression subtractive hybridization.
| INTRODUCTION |
|---|
|
|
|---|
A survey of 31 CF centres in England and Wales, in which over 1200 isolates of P. aeruginosa were analysed, has identified the LES as the most common clone, being present in 48 % of CF centres and accounting for 11 % of the isolates (Scott & Pitt, 2004). The strain has also been identified amongst CF isolates in Scotland (Edenborough et al., 2004). The antibiotic susceptibilities of LES isolates vary widely (Panagea et al., 2003; Salunkhe et al., 2005), and the gold-standard typing method, macrorestriction analysis using PFGE, is time-consuming and requires specialist equipment. A PCR-based diagnostic test for the LES (Parsons et al., 2002) has been evaluated in studies with a large number of CF isolates and found to be 100 % in concordance with identification using PFGE (Panagea et al., 2003; Scott & Pitt, 2004). Furthermore, the test has successfully been applied directly to patient sputum samples (Panagea et al., 2003). The PCR assay is used to amplify one of a number of sequences (named PS21) identified previously using suppression subtractive hybridization (SSH), a technique that enables the identification of sequences present in one strain but absent from another (Parsons et al., 2002; Winstanley, 2002). However, more recently, we and others have identified a small number of non-LES isolates that are positive using the diagnostic test (Lewis et al., 2005; F. W. Scott and T. L. Pitt, unpublished data).
P. aeruginosa contains a core genome that possesses a highly conserved backbone comprising the vast majority of the genome and including most of the recognized virulence factors (Ernst et al., 2003; Wolfgang et al., 2003), and an accessory genome (or mobilome) that includes mobile genetic elements and genomic islands (Larbig et al., 2002; Spencer et al., 2003). The accessory genome of the LES is likely to contain not only strain-specific sequences that could be targeted for identification, but also novel genomic islands that may contribute to the extra abilities possessed by this strain. A genome sequence project to identify all such regions has been initiated in collaboration with the Sanger Centre, Hinxton, UK (http://www.sanger.ac.uk/Projects/P_aeruginosa/).
In this study, we report the further characterization of PS21 and two other previously identified subtracted sequences, and report the results of further SSH to identify a much larger number of sequences present in the accessory genome of LES. By using preliminary data from the genome sequence project we further group many of these sequences into islands or gene clusters. Finally, we combine the data to report the development of improved PCR assays to identify the most common epidemic strain of P. aeruginosa in the UK CF community.
| METHODS |
|---|
|
|
|---|
|
Construction and screening of subtraction libraries. Genomic DNA was isolated from P. aeruginosa LES400 and strain PAO1 as described previously (Winstanley & Hart, 2000). SSH was carried out using the CLONTECH PCR-Select Bacterial Genome Subtraction kit (Clontech) as recommended by the supplier, but with a hybridization temperature of 73 °C. In the hybridizations, DNA from the LES was used as the tester and DNA from strain PAO1 was used as the driver. PCR amplicons obtained following SSH were cloned into pGEM-T (Invitrogen). The subtraction library of RsaI fragments thus constructed was screened by sequencing of plasmid DNA extracted from individual clones (Lark Technologies) using M13 forward and reverse vector primers. BLAST searches at the P. aeruginosa PAO1 genome project web site (http://www.pseudomonas.com) were used to determine the presence or absence of sequences in the PAO1 genome. Tester-specific sequences were further analysed using BLASTN and BLASTX searches of the general database using the site http://www.ncbi.nlm.nih.gov. The relative positions of sequences in the LES genome were determined by BLASTN searching the preliminary data available at the site http://www.sanger.ac.uk/cgi-bin/blast/submitblast/p_aeruginosa. Subsequent analysis of downloaded sequence data was carried out using the program ARTEMIS.
PCR amplification screening of strains.
Oligonucleotide primers (Sigma-Genosys) used in PCR assays are listed in Table 2
along with the annealing temperatures used. DNA for PCR amplification was prepared by making a suspension of a few colonies in 200 µl 5 % Chelex-100 (Sigma) solution. After vigorous mixing, the suspension was boiled for 510 min. Following centrifugation, 150 µl of supernatant containing the DNA was removed and stored at 20 °C. Typically, 1 µl of this DNA was used directly in 25 µl volumes containing 1.25 units Taq DNA polymerase (Promega), 1x TaqMaster (Helena Biosciences), 300 nM each primer, 1x Taq buffer, 2.5 mM MgCl2 and 100 µM nucleotides (dATP, dCTP, dGTP, dTTP). Amplifications were carried out in an Eppendorf MasterCycler thermal cycler for 30 cycles consisting of 95 °C (1 min), annealing temperature (1 min) and 72 °C (2 min), with an additional extension time at 72 °C (10 min) following completion of the 30 cycles.
|
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
It has been demonstrated that the insertion of a transposon into the FHA-like genes PA0041, PA2462 and PA4625 results in reduced virulence in a rat model of chronic respiratory infection (Potvin et al., 2003), suggesting a role for each of these genes in pathogenicity. By extending SSH sequence PSA, we found that the LES equivalent of PA2462 contained an alternative C terminus as well as an alternative variable region further upstream. By screening data now available from the LES genome sequence project, we have been able not only to confirm this variation, but also to identify downstream of the PA2462 ORF a region of repeated FHA-like sequences, including an identical truncated C terminus equivalent to the PAO1 C terminus of ORF PA2462. In reporting the genome sequence of N. meningitidis serogroup A strain Z2491, Parkhill et al. (2000) refer to repeat mediated rearrangement of the 3' ends of a gene encoding a surface-exposed FHA-like protein. The sequence downstream of this gene contained complex repeated sequences also found within the gene. Sequencing of another N. meningitidis strain (Z4259, serogroup C) revealed a similar repeat structure, but with additional DNA inserted into the 3' sequence of NMA0688 at one of the repeats. This suggested that the 3' end of the gene can be altered by recombination with alternative 3' ends. Furthermore Klee et al. (2000) have reported that the C terminus of the FHA-like protein identified as similar to the PSA sequence (fhaB of N. meningitidis) is variable between strains. The authors hypothesize that repeated sequences could provide for variation in the C-terminal portion of the FhaB protein. They also report evidence that the fhaB gene is present in two copies in at least some strains of N. meningitidis.
The PA2462 peptide sequence shares significant similarity to Pseudomonas putida KT2440 ORFs PP0168 and PP1449 (Nelson et al., 2002). PP1449 (1508 amino acids) is the better match, and is also identified when the KT2440 genome is interrogated with PA0041 (E value of approximately 1x1046). PP0168 is the largest gene (8682 amino acids) in the KT2440 genome and has the most unusual sequence (Weinel et al., 2002). The protein is rich in asparagine, valine and threonine, but contains no cysteine and is low in arginine. It contains two large threonine-rich repeats: the N-terminal repeat (nine units, each 100 amino acids long) and the C-terminal repeat (29 units, each 219 amino acids long). Although it has much in common with the P. aeruginosa genome, the only virulence-related genes that the genome of the avirulent P. putida KT2440 contains encode putative or confirmed adhesion proteins. Both PP0168 and PP1449 share homology with adhesin genes of enteropathogens and encode adhesins that are essential for seed colonization (Espinosa-Urgel et al., 2000; Nelson et al., 2002).
Thus, FHA-like proteins such as PA2462 may contribute to adhesion and colonization of disparate hosts, whether seed or human, and contain repeat sequences in the C terminus and downstream of the ORF with a probable role in mediating recombination. The biological function of such variations remains unknown, but since there may be a role for such proteins in the colonization of CF patients, further study is certainly merited.
LES subtracted library
Subtraction of the genome of strain PAO1 from that of isolate LES400 produced a library of clones, 106 of which were subjected to plasmid extraction and sequencing. Following the removal of duplicates and non-Tester-specific sequences, the subtraction added 78 sequences to the 14 obtained in a previous study (Parsons et al., 2002). The combined subtracted sequences are summarized in Table 3
. Fifteen of the sequences shared similarity with mobile genetic elements, including bacteriophages D3 and D3112, both of which have been reported previously in P. aeruginosa (Cavenagh & Miller, 1986; Wang et al., 2004). A further six subtracted sequences were associated with transposition, recombination or DNA modification. Amongst the known variable genes matching subtracted sequences were several related to siderophores. We have reported previously that the LES has a type III pyoverdine receptor (Parsons et al., 2002). The SSH identified sequences related to the synthesis of both type III pyoverdine and the corresponding receptor (Smith et al., 2005). In addition, the LES carries genes for serotype O6. Four of the subtracted sequences shared similarity with serotype-related genes. Whilst four subtracted sequences indicated the presence of the known P. aeruginosa genomic islands PAGI-1 (Liang et al., 2001) and PAGI-2 (Larbig et al., 2002), a further five sequences shared similarity with four linked genes from Clostridium thermocellum ATCC 27405, suggesting that the LES may contain a low-G+C-content genomic island related to a gene cluster in this Gram-positive bacterium (Table 3
). The subtraction also identified variations in the pilABC cluster. Such variations between strains of P. aeruginosa have been reported previously (Choi et al., 2002; Kus et al., 2004).
As well as the gene related to mercury resistance associated with PS21, the subtraction identified a putative gene involved in silver resistance. Three subtracted sequences shared similarity with proteins that have roles in the production of polyketides (Table 3
). Two of these sequences matched a regulatory protein (PltR) involved in the synthesis of pyoluteorin, which has antifungal activities (Bender et al., 1999). BLASTX searches for 38 % of the LES subtracted sequences revealed either no significant matches or matches with hypothetical proteins of no known function.
Preliminary searches of the LES genome using SSH sequences
Preliminary searches of the LES genome sequence using BLASTN identified all but two of the SSH sequences. Interestingly, there was no BLASTN match in the LES genome (99.98 % theoretical coverage) for either LESC4 or LESC5, which match adjacent genes from bacteriophage D3112, suggesting that this bacteriophage may be present in isolate LES400 but missing from the genome-sequenced isolate (LESB58). Notably, many of the SSH sequences were found on common contigs containing putative genomic islands. A single contig contained, in order, the subtracted sequences LESC6, LESH12, PS54, LESF10, LESG7, LESE5, PS3, LESF1, LESB11, LESE3, LES11, LESF12, PS1, PS12, LESD8, LESB7 and LES4 (cluster 1). Further analysis indicated that sequence LESC6 is located at one edge of a putative genomic island incorporating all of the other SSH sequences with the exception of LES4, which lies beyond the putative island. The diagnostic probe (PS21) appears to lie in a cargo region of the previously identified island PAGI-2 (Klockgether et al., 2004; Larbig et al., 2002) that also contains LESB9, the putative silver-efflux-related sequence (cluster 2). Twelve of the SSH sequences (in order: LESH1, LESE4, LESB4, LESD11, LESG2, LESF9, LESD7, PS32, LES5, LESA8, LESE6 and LESB2) were clustered on a different contig within a second large putative island containing a number of genes related to C. thermocellum (cluster 3). There were three separate contigs containing multiple SSH sequences, some of which were bacteriophage related. One contained, in order, LES1, LESA9, LESC1, LESB6, LES6 and LESG5 (cluster 4); a second contained, in order, LESA11, LESH3, LESG11, LESE2, LESG3, LESA2, LESC11, LESC10, LESH5, LESA12 and LESF3 (cluster 5); and a third contained, in order, LESC3, LESE12, LESE7, LESF8, LESF6 and LESG9 (cluster 6). There were also several other contigs containing more than one SSH sequence. Not surprisingly, one of these linked the type III pyoverdine-synthesis-related sequences LESH4 and LES8.
Distribution of subtracted sequences
The distribution of 25 LES sequences amongst a panel of strains was assessed by PCR amplification assays using oligonucleotide primers designed from subtracted sequences (Table 4
). The sequences selected included multiple representatives for each of the six major clusters identified by preliminary analysis of the genome sequence, sequence LESD9 (which shared a contig with LESB5, PS15, LESG10 and PS2), sequence LESB3 (which shared a contig with LESH11), and the type III pyoverdine-synthesis-related sequences LESH4 and LES8, and LESA3 (which shared a contig with LESA1).
|
Cluster 2, including PS21 and the silverefflux-pump-related SSH sequence LESB9, was present in all LES isolates. All four PS21 false-positive strains were also PCR positive for LESB9, confirming the linkage of these two sequences. All SSH sequences tested from cluster 3 were PCR positive for all LES isolates but PCR negative for all other genetically characterized isolates (Table 4
). Further PCR assays for sequences LESE4 and LESF9 were used to screen the uncharacterized non-LES CF isolates from Liverpool and from Clatterbridge Hospital. Three isolates (CH31, CH47 and 49278) were PCR positive for LESE4, and two of these isolates, CH47 and 49278, were also PCR positive for LESF9. All other PS21-negative isolates were PCR negative for each of these sequences. PFGE typing of isolates CH31, CH47 and 49278 confirmed that these did not match the LES pulsotype (data not shown).
The PCR-based distribution profiles for the three cluster 4 sequences were similar but not identical, and indicated some instability. Two LES isolates were PCR negative for all three sequences. Two of the three cluster 5 sequences (LESA11 and LESG11) shared identical distributions and were stable amongst LES and Midlands1 isolates. The third sequence tested from cluster 5 (LESH5) was greater than 50 kb away from the other two. Only LES isolates were PCR positive for this sequence. The three cluster 6 SSH sequences shared nearly identical distributions and were detected in all LES isolates. The other sequences screened were present in both LES and Midlands1 isolates. The Manchester epidemic strain isolates were PCR positive for only one of the sequences (LESA3, matching a putative siderophore receptor) tested (Table 4
).
PCR-based identification of the LES
We have reported previously that the LES strain carries approximately 95 % of the ORFs present in strain PAO1 and that estimations based on PFGE indicate that the LES genome is no bigger than that of PAO1 (Salunkhe et al., 2005). In our previous study (Parsons et al., 2002), we identified 14 LES sequences subtracted from strain PAO1, including five (PS1, PS3, PS12, PS21 and PS54) that were LES specific. All of these, except PS21, were absent from one of the LES isolates screened (isolate 109). After considerably expanding the SSH database and conducting preliminary searches of the LESB58 genome sequence, it is clear that sequences PS1, PS3 and PS54, along with 13 other SSH sequences, are part of an LES-specific island that is associated with some instability. Thus, although we have never found this island in strains other than LES, any test for the LES based on this island may be prone to false negatives. In contrast, PS21, part of a merA gene located in an apparent cargo region of PAGI-2, is stable in the LES, but occurs rarely in other P. aeruginosa strains, making tests based on PS21 theoretically susceptible to rare false positives. Sequence LESF9 forms part of a putative genomic island present in all LES isolates but is, like PS21, occasionally also found in non-LES isolates. In the absence of a definitive LES-specific PCR assay that is positive for all LES isolates tested, we suggest a strategy using a combination of PCR assays for the detection of PS21 and LESF9. Since PS21 has proved to be an effective tool, especially in a setting where LES is widespread (Panagea et al., 2003; our unpublished data), we would propose the continuation of this test. For clinicians and patients alike, it is important that all LES isolates give a positive result. In order to achieve this, in the light of the instability of the LES genome, we may have to accept the possibility of rare false positives. However, if a combination of PCR assays for PS21 and LESF9 is used, then at least with our collection of isolates, only LES isolates are PCR positive for both of these sequences. We would suggest that any isolates that are PCR positive for just one of these sequences should be further tested using a different typing method, such as PFGE or RAPD.
The LES mobilome
The SSH data indicated that bacteriophage DNA is likely to make a significant contribution to the accessory genome of the LES. D3 is a lysogenic bacteriophage with a genome of 56.4 kb (Kropinski, 2000) (accession no. AF165214) implicated previously in serotype conversion (Newton et al., 2001) and transduction (Cavenagh & Miller, 1986) in P. aeruginosa. SSH sequences matching widely separated regions of the D3 sequence indicate the possible presence of this bacteriophage genome within the LES. D3112 is a transposable Mu-like bacteriophage with a genome of 37.6 kb (Wang et al., 2004). We identified LES sequences matching neighbouring genes within the D3112 genome, but we have no further indication concerning the extent of the D3112 DNA within the genome of the LES. Interestingly, these two sequences appear to be missing from the LES isolate chosen for genome sequencing, suggesting either that the bacteriophage has been acquired by isolate LES400 or that it has been lost by isolate LESB58. In addition, the bacteriophage-related cluster 4 was absent from some LES isolates. SSH also indicated that the LES contains a putative tail-fibre protein from another bacteriophage reported elsewhere in P. aeruginosa (Nakayama et al., 2000). The presence of several other subtracted sequences with matches sharing lower identities with bacteriophage-related proteins, many clustered on to one of three contigs, suggests that there may be additional whole (or fragments of) bacteriophage genomes within the LES. The details will become clear when the ongoing genome sequence project is further advanced.
| ACKNOWLEDGEMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
C. Winstanley, M. G.I. Langille, J. L. Fothergill, I. Kukavica-Ibrulj, C. Paradis-Bleau, F. Sanschagrin, N. R. Thomson, G. L. Winsor, M. A. Quail, N. Lennard, et al. Newly introduced genomic prophage islands are critical determinants of in vivo competitiveness in the Liverpool Epidemic Strain of Pseudomonas aeruginosa Genome Res., January 1, 2009; 19(1): 12 - 23. [Abstract] [Full Text] [PDF] |
||||
![]() |
K Mohan, J L Fothergill, J Storrar, M J Ledson, C Winstanley, and M J Walshaw Transmission of Pseudomonas aeruginosa epidemic strain from a patient with cystic fibrosis to a pet cat Thorax, September 1, 2008; 63(9): 839 - 840. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Kukavica-Ibrulj, A. Bragonzi, M. Paroni, C. Winstanley, F. Sanschagrin, G. A. O'Toole, and R. C. Levesque In Vivo Growth of Pseudomonas aeruginosa Strains PAO1 and PA14 and the Hypervirulent Strain LESB58 in a Rat Model of Chronic Lung Infection J. Bacteriol., April 15, 2008; 190(8): 2804 - 2813. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Tralau, S. Vuilleumier, C. Thibault, B. J. Campbell, C. A. Hart, and M. A. Kertesz Transcriptomic Analysis of the Sulfate Starvation Response of Pseudomonas aeruginosa J. Bacteriol., October 1, 2007; 189(19): 6743 - 6750. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. H. M. Smart, F. W. Scott, E. A. Wright, M. J. Walshaw, C. A. Hart, T. L. Pitt, and C. Winstanley Development of a diagnostic test for the Midlands 1 cystic fibrosis epidemic strain of Pseudomonas aeruginosa. J. Med. Microbiol., August 1, 2006; 55(Pt 8): 1085 - 1091. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | J MED MICROBIOL | MICROBIOLOGY | J GEN VIROL | ALL SGM JOURNALS |