Previous Article | Next Article ![]()
Journal of Bacteriology, December 2008, p. 8053-8064, Vol. 190, No. 24
0021-9193/08/$08.00+0 doi:10.1128/JB.00834-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
Neil Molyneaux,1
Kristine M. Hujer,2
Heather Lavender,1
Jennifer J. Jamison,3
Ian J. MacDonald,3
Kristienna M. Martin,3
Thomas Russo,3,4,5,7
Anthony A. Campagnari,3,7
Andrea M. Hujer,2
Robert A. Bonomo,2,6 and
Steven R. Gill3,8
Department of Genetics,1 Departments of Pharmacology and Molecular Biology and Microbiology, Case Western Reserve University School of Medicine, Cleveland, Ohio,6 Louis Stokes Cleveland Department of Veterans Affairs Medical Center, Cleveland, Ohio,2 New York State Center of Excellence in Bioinformatics and Life Sciences,3 Veterans Administration Western New York Healthcare System,4 Department of Medicine,5 Department of Microbiology & Immunology,7 Department of Oral Biology, The State University of New York at Buffalo, Buffalo, New York8
Received 16 June 2008/ Accepted 6 October 2008
|
|
|---|
|
|
|---|
A significant nosocomial outbreak of MDR A. baumannii in the United States occurred in the period 2003-2005 in a military treatment facility caring for injured service personnel and civilians (Walter Reed Army Medical Center [WRAMC]) (19, 40). Molecular typing of isolates from this outbreak revealed eight major clone types, and about 60% of the isolates were related to three pan-European types (5, 52), suggesting multiple independent origins (40). Examination of specific resistance determinants in the WRAMC isolates demonstrated considerable variability in the composition of resistance genes within each clone type and similar patterns across certain divergent clone types. Thus, genetic relatedness was a poor predictor of the MDR phenotype. This led to the hypothesis that there exist multiple independent genetic mechanisms leading to MDR in A. baumannii. A large cluster of antibiotic resistance genes and mobile genetic elements is present as an 86-kb "resistance island" (RI) in the A. baumannii AYE genome (15). The RI is not present in the genomes of the drug-susceptible A. baumannii isolates ATCC 17978 (43) and SDF (51), suggesting that it is an important contributor to the MDR phenotype.
In an effort to understand the significance of the RI and the integration of genetic factors that permit A. baumannii to successfully infect patients in multiple clinical settings, we undertook a whole-genome sequence analysis of six clinical isolates in comparison to a nonhuman A. baumannii isolate (SDF) and an environmental Acinetobacter strain, Acinetobacter baylyi ADP1. We compared three MDR isolates to three drug-susceptible isolates in order to explore common features of their evolution and the diversity of mechanisms by which resistance is demonstrated. Our goals were (i) to determine the repertoire and phylogeny of resistance determinants among recent clinical isolates, (ii) to see how these "successful" clinical isolates differ from pathogenic and nonpathogenic strains/species already examined, (iii) to assess the function and purpose of the RI, and (iv) to evaluate the importance and impact of insertion sequences (IS) in the evolution of the resistance phenotype. Our results show that the RI is variable in composition and is only one contributor to the MDR phenotype. Based on closely matched MDR and drug-susceptible isolates, it appears that the MDR phenotype is largely due to the acquisition of a series of genes encoding drug-inactivating enzymes that are located on mobile genetic elements. Moreover, we identified a series of genes that may define the ability of A. baumannii to live in association with the human host and describe extensive diversity in genes that contribute to the lipopolysaccharide barrier.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Selected genetic and phenotypic properties of A. baumannii genomes
|
Genome sequencing and annotation. AB307-0294 and AB900 were sequenced by pyrosequencing, and AB0057 was sequenced by a combination of pyrosequencing and traditional Sanger shotgun sequencing. Finished genome sequences with no gaps were completed for AB307-0294 and AB0057. Genome annotation was performed at The Institute for Genomic Research (TIGR) and later at the J. C. Venter Institute, Rockville, MD, and the Institute for Genome Sciences, Baltimore, MD, using the TIGR Annotation Engine (www.jcvi.org/ems/research/projects/annotation_service/overview/) followed by manual curation. Details of the sequencing, assembly, gap filling, and annotation methods are given in the supplemental material.
PCR. PCR was performed using the primers listed in Table S1 in the supplemental material and as described by Hujer et al. (19).
Nucleotide sequence accession numbers. The genome sequences of AB0057, AB307-0294, and AB900 have been deposited in the DDBJ/EMBL/GenBank databases under the accession numbers CP001182-3, CP001172, and ABXK00000000, respectively.
|
|
|---|
Despite acquisition from geographically diverse regions, AYE, AB0057, and AB307-0294 are remarkably similar at the DNA level, averaging over 99.9% nucleotide sequence identity in orthologous regions (Fig. 1). The other isolates are each
98% identical to one another. Interestingly, the >50-year-old isolate (ATCC 17978) and the louse isolate (SDF) are not more divergent than any of the other isolates at the sequence identity level (Fig. 1A). The relationship of isolates to one another based on shared gene content differs from the relationships inferred based on analysis of conserved genes, reflecting a considerable contribution of lateral gene transfer in the evolution of A. baumannii (Fig. 1B). SDF has a much smaller genome than the other A. baumannii isolates, presumably as a result of extensive IS-mediated deletion events (51). The extent of divergence in gene content in SDF and its questionable role as a human pathogen have prompted us to focus most of our analysis on comparisons among the six human clinical isolates.
![]() View larger version (9K): [in a new window] |
FIG. 1. Phylograms of A. baumannii isolates. The relationships of A. baumannii isolates were inferred based on the similarity of shared sequences that have most likely been retained in each genome from a common ancestor (A) and based on the extent of unique gene content in each isolate (B). (A) A maximum likelihood tree was constructed using dnaml from the PHYLIP package, based on the concatenated DNA sequences of 1,942 ORFs that were reciprocal best matches in each genome. Bar, 0.01 substitution per nucleotide. (B) A tree based on the shared content of accessory genes was constructed as described previously (44). Bar, 10% difference in gene content.
|
![]() View larger version (18K): [in a new window] |
FIG. 2. GO categories of core and A. baumannii-specific genes. GO categories were inferred for AB0057 genes by using the JCVI Annotation Engine. The percentage of genes in each GO category was determined for core A. baumannii genes (those present in all six clinical isolates), Acinetobacter core genes (A. baumannii core genes that are also present in ADP1), and pan-A. baumannii accessory genes (core genes absent from ADP1), using WEGO (58). Percentages indicate the fractions of gene groups annotated with each GO category. Stars indicate GO categories that are significantly underrepresented (open) or overrepresented (filled) among pan-A. baumannii accessory genes compared to core genes (P < 0.05).
|
Adaptation to the human host. A. baylyi ADP1 is a naturally transformable soil isolate. We hypothesize that genes found in the A. baumannii isolates but not in ADP1 may contribute to features of growth in association with a human host. We found 475 genes that are shared by all six A. baumannii clinical isolates and that are missing from the ADP1 genome (pan-A. baumannii accessory genes [see Table S2 in the supplemental material]). Hence, we propose that these genes may play an important role in human adaptation. Reinforcing the significance of these genes for growth in association with humans, 279 (59%) have been lost in the SDF strain, which was isolated from a human body louse.
Gene annotations for each genome were classified using gene ontology (GO) terms to group related genes by molecular function, biological role, and cellular compartment (2). The pan-A. baumannii accessory genes are enriched in the categories of transport activity and transcription factor activity compared with core genes shared with ADP1 (Fig. 2). In contrast, GO categories for metabolism and cellular processes were overrepresented among core genes compared with pan-A. baumannii accessory genes. Remarkably, 59 of the pan-A. baumannii accessory genes are predicted to be transcription factors based on sequence similarity. This suggests a more extensive regulatory program in A. baumannii than in ADP1. Transporter genes comprise 22% of the pan-A. baumannii accessory genes, compared to 10% of all genes. Transport systems are identifiable for urea, taurine, oligopeptides, and histidine. The presence of additional transport systems for amino acids and other nitrogen-containing compounds presumably reflects their greater availability in humans than in soil. Several amino acid utilization enzymes were also found, including arginine N-succinyltransferase and histidine decarboxylase.
Two hundred sixty-one of the 475 pan-A. baumannii accessory genes (55%) are present in clusters of at least four adjacent genes. As noted above, not all clusters should be presumed to be present due to lateral transfer into the A. baumannii ancestor. Although by definition none have orthologs in ADP1, 40% have orthologs in Pseudomonas aeruginosa PAO1 and may have been lost in the ADP1 lineage.
Two clusters are worthy of further note. One is a cluster of genes responsible for the biosynthesis of a homoserine lactone that is involved in quorum sensing (AB57_0151 to AB57_0161). The AbaI protein (AB57_0151) has been shown to be a homoserine lactone synthase which is involved in quorum sensing and is required for biofilm formation (33). The second cluster (AB57_2565 to AB57_2574) is a chaperone/usher pilus assembly gene cluster that has also been shown to be involved in adherence to mammalian cells and in biofilm formation (50).
Variation among A. baumannii isolates. As shown for Haemophilus influenzae, P. aeruginosa, Streptococcus pneumoniae, and other bacterial pathogens, there is considerable variation in the DNA content of each isolate (18, 29, 41, 42). The accessory genome of A. baumannii comprises 2,649 annotated genes that are present in one to five of the six isolates. This number will surely grow as additional isolates are characterized. The number of unique genes ranges from 80 (AB307-0294) to 528 (ATCC 17978).
Genomic regions of unique content (islands) have been noted in each of the reported A. baumannii genomes, but the availability of seven complete genomes enables a much more rigorous evaluation of the contents and distribution of these islands. A variety of terms have been applied to clustered groups of genes in A. baumannii, including (putative) alien islands, pathogenicity islands (PAIs), and RIs (11, 15, 23, 39, 43). Based on the gene content, the latter term seems appropriate for a subset of islands with a clear concentration of genes associated with resistance to several classes of antibiotics. However, the function of most islands is not proven, so we prefer the simple term genomic islands.
In AB0057, 60% of lineage-specific genes are within clusters of at least 10 genes located adjacent to one another. Some genomic islands are likely the result of lateral gene transfer, such as prophage insertions, which in general are the largest genomic islands in each isolate. The RI described below is comprised largely of genes that were likely acquired by lateral transfer. It is also clear that not all genomic islands were derived by lateral gene transfer, however. Some likely represent regions of an ancestral Pseudomonadales genome that have been lost in ADP1, for example, a leucine metabolism cluster at AB57_1597 to AB57_1604 that has orthologs in both P. aeruginosa and another closely related organism, Psychrobacter cryohalolentis (AB0057 numbers are given, but this cluster is shared by all of the A. baumannii isolates except for SDF).
Each genome contains a considerable number of unique genes, as 46% of accessory genes are isolate specific. A disproportionate number of these are hypothetical and are likely cryptic prophage regions. Of the 1,228 total genes that are annotated for only one of the six clinical isolates, 412 fall within prophage regions and 487 others are annotated as encoding hypothetical proteins (73% in the two categories combined). The phage regions are not orthologous, suggesting recent insertion and/or rapid loss and a large pool of potential bacteriophage genomes. One prophage was found in both AB0057 and ACICU, albeit at different chromosome locations in the two genomes (AB57_1230 to AB57_1300 and ACICU_02150 to ACICU_02228). The significance of these two independent insertion events in the context of the natural history of these isolates is unclear.
The excess of hypothetical proteins in the accessory genome could be due to several causes. First, incorrectly annotated (wrong frame, strand, or region) open reading frames (ORFs) that do not code for bona fide proteins will not match ORFs from the other isolates. Second, incomplete or inaccurate annotation can result in the unintended omission of identification of a set of orthologous proteins. Finally, frameshift errors in the DNA sequence can result in unpaired orthologs. The first two groups of hypothetical proteins are difficult to quantify. The number of frameshifted annotations was examined by comparison of the predicted protein set from each isolate to the genomes of the other isolates, using TBLASTN. Instances where the alignment was apportioned across two different frames were counted. The number of potentially frameshifted annotations varied from 21 in AYE to 162 in ATCC 17978 (see Table S3 in the supplemental material).
A. baumannii is thought of as a cousin of another human pathogen, P. aeruginosa, and 65% of the core genes have orthologs in P. aeruginosa PAO1. The accessory genome of P. aeruginosa has been shown to be inserted into limited locations around the genome, defined as regions of genome plasticity (29). In contrast, there does not seem to be a similar preference for strain-specific insertion locations in A. baumannii, with the exception of the location of the RI.
RIs of variable composition. An interesting feature of many drug-resistant genomes is the presence of a genomic island or RI containing a collection of genes encoding proteins related to antibiotic inactivation and efflux. RIs can be plasmid carried (20, 22) or located on the bacterial chromosome (15), and antibiotic resistance genes are usually interspersed with mobile genetic elements (11).
Comparison of the sequence and structure of the genomic island in each A. baumannii isolate revealed a series of insertion/deletion events at this locus (Fig. 3A). The 86-kb island in AYE (AbaR1) is the largest, with 90 annotated genes, including genes related to the inactivation of β-lactams, aminoglycosides, chloramphenicol, rifampin, and tetracycline. These resistance determinants are intermixed with a complex set of partial and complete integrons and transposons (15). The 49-kb AB0057 RI (AbaR3) is largely a subset of AbaR1. It contains eight genes associated with antibiotic resistance. Unique sequences in AbaR3 include a blaTEM gene that is associated with a Tn3 transposon and a small cluster of genes, including two that encode a DNA topoisomerase and a single-strand binding protein that are most similar to proteins from a broad-host-range plasmid (49). The ACICU RI contains the distal portion of the AYE/AB0057 RI, into which an ISAba1-flanked blaOXA-20 gene has been inserted (AbaR2). In ATCC 17978, only the initial and terminal segments are present, and these do not contain resistance determinants. An unrelated 37.7-kb insertion is present at the location of the RI in AB900. The AB900 island appears to have been transferred laterally based on the source of the best protein database matches. It contains a copper resistance operon and a high concentration of hypothetical proteins. ADP1 and AB307-0294 preserve the ancestral gene order across the location of the RI insertion. Given the close relationship between AYE and AB307-0294, it is likely that the larger AYE RI is more like an original insertion, with AB0057 and ACICU derived by series of subsequent deletion, insertion, and rearrangement events. Another interpretation, however, is that the RI represents a series of independent insertion events; this view is supported by the fact that ACICU and AYE/AB0057 are representatives of divergent European clone types II and I, respectively.
![]() View larger version (19K): [in a new window] |
FIG. 3. RI composition. (A) Modules that comprise components of the RI in sequenced A. baumannii genomes are designated by the letters A to J. A bold bar indicates that the region is present in an isolate. (B) Representative genes present in each module are listed. (C) The presence of various RI components in additional A. baumannii isolates was tested by PCR. *, PCR results were inferred based on a BLAST search of each genome with primer sequences; a, these primer pairs flank the island location, so successful PCR amplification (+) indicates that the island, an insertion that is too long for PCR amplification, is absent. AB0057 was reported previously as AB057 (19). Abbreviations: AMP, ampicillin; CIP, ciprofloxacin; CAZ, ceftazidime; FEP, cefepime; AMK, amikacin; TOB, tobramycin; SAM, ampicillin-sulbactam; MEM, meropenem; IPM, imipenem. Hujer, data obtained from reference 19.
|
A distinguishing feature of AbaR1 and AbaR3 in the MDR isolates AYE and AB0057 is the presence of a pair of ISPpu12 IS elements that flank the resistance gene portion of the island (Fig. 3A, fragments D to J). In the second AB0057 island (AbaR4), the uspA gene, encoding a universal stress protein, is present intact, while this gene is interrupted at the corresponding position in the main RI (AbaR3) by the ISPpu12-flanked insertion. The average G+C content of ORFs between the ISPpu12 elements is 57%, compared to 39% for segments A to C and K of the RI and the rest of the genome, suggesting a recent acquisition of the ISPpu12-flanked cassette by lateral transfer. Many of the predicted proteins in segments E to I are nearly identical to segments of a broadly distributed integron (45). In addition, as noted above, AbaR3 in AB0057 contains genes for a plasmid-derived DNA topoisomerase that may contribute to RI mobility. Only one ISPpu12 element is present in ACICU, perhaps reflecting deletion of one element along with a significant fraction of the 5' end of the island, including the 5'-flanking sequence.
How much of the MDR phenotype can be explained by the RI? The RI is a dramatic feature of the MDR isolates AB0057 and AYE. A significant concentration of genes involved in antibiotic resistance is also present in an orthologous location in ACICU. It is tempting to speculate that the MDR phenotype is correlated with the presence of some variant of the RI. To test this hypothesis, we attempted to amplify segments of the RI from several drug-susceptible and MDR isolates (Fig. 3B). The basic structure of the RI, as found in AB0057 and AYE, appears to be restricted to closely related clones, as defined by PCR-electrospray ionization mass spectrometry (PCR/ESI-MS) typing (12). In other PCR/ESI-MS types, however, MDR isolates typically exhibit the presence of some of the resistance genes but not other markers of the RI per se, such as the ISPpu12 elements that flank the resistance gene region. It is possible that these genes either are not associated with an RI or are associated with islands with an alternate structure and/or location. One example is the presence of aadB in several isolates despite the absence of blaVEB-1, which is adjacent to aadB in the AYE RI, indicating that this segment of the RI is not widespread (Fig. 3C, RIN primers). The RI is completely missing in AB307-0294 and a subset of other isolates that are largely susceptible (Fig. 3C, RIH primers). An insertion is present at this genomic location in a large majority of isolates. The genes present may be quite different, though, since completely unrelated sequences are present at this genome location in AB900 and SDF.
In the three MDR isolates, AB0057, AYE, and ACICU, a significant fraction of genes responsible for antibiotic resistance are located in mobile genetic elements, including the RIs in each isolate and a plasmid in ACICU. In AB0057, a second gene encoding ADC (blaADC-40) is present, flanked by ISAba1 elements. The blaADC-40 gene was likely obtained by lateral transfer because it contains a 15-bp insertion relative to most other blaADC genes as well as four amino acid substitutions and several silent changes compared with the chromosomal copy.
The corollary issue, i.e., how many resistance genes are not associated with mobile genetic elements, is equally compelling because it may provide a sense of the innate resources in the genome that are available to counter the antibiotic threat. Resistance to fluoroquinolones is mediated by mutations in gyrA and parC, which are not present on mobile elements. Resistance to aminoglycosides and β-lactams may be mediated by a combination of drug-inactivating enzymes, efflux pumps, and changes in the bacterial outer membrane that reduce drug uptake. Three multidrug efflux pumps, encoded by adeABC, adeIJK, and abeM, have been shown to contribute to drug resistance in A. baumannii (9, 27, 32, 46).
The very close relationship of the drug-susceptible isolate AB307-0294 with the MDR isolates AB0057 and AYE enables another perspective on the acquisition of resistance determinants. Virtually all of the AB307-0294 genome aligns with the AB0057 and AYE genomes, and the average percent identity between orthologs is >99.7%. The genes that are present in AYE and AB0057 but not in AB307-0294 are almost exclusively related to the RIs or to prophage or IS element insertions (see Table S4 in the supplemental material). Only
3% of the AB307-0294 genome is not represented in AB0057 or AYE; this is virtually all in two prophage clusters. The significance of the genetic similarity among these three isolates rests with the observation that AB307-0294 retains all of the transport, metabolism, and regulatory capacities of the MDR isolates and yet is drug susceptible.
The adeABC and adeIJK genes encoding drug efflux pumps are present and the encoded proteins are 100% identical in sequence in all three isolates. The two-component regulatory proteins AdeR and AdeS are also present in AB307-0294, with only a single conservative amino acid substitution in each protein differing from the AB0057 sequence. In fact, the entire adeRSABC and adeIJK gene regions together have only three nucleotide differences in 12,436 bp between AB307-0294 and AB0057 or AYE, suggesting that both their regulation and function are similar among the three isolates. It has been suggested that the AdeIJK proteins contribute to intrinsic but not acquired antibiotic resistance in A. baumannii (9), and these proteins are also identical in AB307-0294, AB0057, and AYE. Regulators of AdeIJK expression have not been described, and it is likely that additional proteins that control the adeABC and adeIJK operons remain to be found.
There are 15 genomic islands in ACICU annotated as being involved in drug resistance (20). For most of these islands, the designation is based on the presence of one or more genes annotated as transporters of the DMT superfamily (21), and their substrates have not been verified experimentally. Most of these genes are also present in the drug-susceptible isolates AB307-0294, AB900, and ATCC 17978 and thus do not seem to be indicative of antibiotic susceptibility status.
In addition to the RIs and prophage regions, there are two additional clusters of genes in the MDR isolates AB0057 and AYE that are absent from AB307-0294: these are a copper resistance cluster and a cluster that includes a glycine/betaine transport operon. Betaines serve as organic osmolytes for protection against osmotic stress, drought, high salinity, or high temperature (57). At least some portion of the copper resistance cluster is present in ATCC 17978 and AB900, and the cluster is absent from ACICU, so these genes are unlikely to contribute directly to antibiotic resistance. The role, if any, of the cluster including the glycine/betaine transport genes in resistance and/or pathogenicity remains to be determined.
IS elements represent another potential source of variability between AB307-0294 and the MDR isolates. AB307-0294 does not contain ISAba1 elements, while AB0057 carries 9 copies and AYE carries 21 copies. Four of the nine ISAba1 elements in AB0057 flank β-lactamase genes that appear to have been transferred laterally into this isolate (blaADC-40 and blaOXA-23). These ISAba1 elements are presumably involved in overexpression of the bla genes from outward-facing promoters in the IS element (17). Two AB0057 ISAba1 insertions are polymorphic and are thus unlikely to be related to drug resistance. One ISAba1 element appears to inactivate the uup gene, in which mutations have been reported to increase the precise excision of transposons (31). The final two insertions are in or adjacent to hypothetical protein genes. None of the ISAba1 insertion locations are shared between AB0057 and AYE, suggesting that any impact on expression of adjacent genes at the insertion sites is not a primary feature of drug resistance.
Other features of antibiotic resistance. The RI carries genes that encode resistance to chloramphenicol, an antibiotic that is not commonly used in contemporary acute care settings in developed countries. This suggests that the RI is not a recent acquisition. Further work will be required to assess the genetic stability of the RI.
All tested isolates of A. baumannii are resistant to aztreonam, which specifically targets penicillin-binding protein 3 (PBP-3) (8). The basis for this resistance has yet to be determined but could be related to the considerable divergence of PBP-3 in Acinetobacter. The AB0057 PBP-3 protein is only 40% identical to the orthologous proteins from the closest non-Acinetobacter relative, Psychrobacter cryohalolentis K5, and to P. aeruginosa PAO1.
All isolates of A. baumannii displayed high-level resistance to nitrofurantoin. Nitrofurantoin is a prodrug that is reduced by the bacterial enzyme nitrofuran reductase to form highly reactive intermediates that in turn produce DNA damage. Although present in the close Acinetobacter relative Psychrobacter sp. strain PRwf-1 and many other gammaproteobacteria, this enzyme is absent from ADP1 and all A. baumannii isolates examined to date.
Plasmids were identified during genome sequencing of AB0057 (1), ATCC17978 (2), AYE (4), SDF (3), and ACICU (2), but only in ACICU is there a direct contribution to antibiotic resistance. Two copies of blaOXA-58 are present on a plasmid in ACICU and are flanked by insertion sequences. No plasmids were found after accounting for all genome sequence reads from AB900 and AB307-0294.
Evaluation of esv genes. Smith et al. identified 16 genes that, when mutated, resulted in attenuation of the virulence of A. baumannii ATCC 17978 in two invertebrate models of infection (43). They also identified 28 potential PAIs based on analysis of base composition and gene content. Surprisingly, only three of the attenuated virulence mutations mapped to PAIs, and the two PAIs containing mutations are prophage regions. Eleven of the 16 so-called esv genes (for ethanol-stimulated virulence) have orthologs in the nonpathogenic organism Acinetobacter baylyi ADP1, and 4 of the remaining 5 are absent from all five recent clinical isolates of A. baumannii discussed here and from six additional tested isolates (data not shown). This suggests that they are not ubiquitous contributors to clinical infections. The one remaining esv gene (esvB; A1S_1232) encodes an HxlR family helix-turn-helix transcription factor of unknown specificity. esvB is also absent from the SDF strain. Further work to determine the role of this gene in pathogenicity of A. baumannii will be required, preferably with a mammalian model system.
Diversity at the O-antigen biosynthetic cluster. The outer surfaces of gram-negative bacteria are comprised of a lipopolysaccharide layer that serves as a barrier between the outer membrane and the environment. Considerable intraspecies variability has been observed in the composition of the O-antigen side chains of the lipopolysaccharide, and this contributes to antigenic variability among isolates. For P. aeruginosa, at least 11 highly divergent gene clusters encoding proteins involved in O-antigen biosynthesis at a conserved genome location have been reported (37). A similar arrangement is found in A. baumannii. Among the seven isolates analyzed here, there are six highly divergent sets of genes located at the O-antigen cluster. In each genome, the cluster is located at an orthologous position (between ACIAD0069 and ACIAD1016 in ADP1) and contains 17 to 30 genes. The sequences of AB307-0294 and AYE are identical across this region. Testing of the O-antigen gene content of other isolates by PCR showed that a subset of PCR/ESI-MS type 15 isolates share the AB0057 gene set, while other isolates match the AB307-0294/AYE gene set (Table 2). Among the closely related isolates defined by PCR/ESI-MS type 15, the amplification pattern is quite consistent across the gene cluster, and all type 15 isolates match either the AB0057 or AYE/AB307-0294 pattern (Table 2). With more diverse isolates, it is less common to see amplification with both of the tested primer pairs for each O-antigen cluster type, implying that the entire cluster is not always present intact (Table 3). Nonetheless, 18 of 20 isolates in the diverse set of isolates tested by PCR were able to amplify at least one of the O-antigen cluster genes (Table 4). Thus, it remains a possibility that a relatively small set of O-antigen gene sets is present among A. baumannii isolates. Previous characterization of the O-antigen diversity among isolates from the A. baumannii-A. calcoaceticus complex by use of monoclonal antibodies showed considerable variation in reactivity, indicating multiple distinct lipopolysaccharide modifications, even among closely related isolates (34).
|
View this table: [in a new window] |
TABLE 2. Distribution of O-antigen genes with first primer set
|
|
View this table: [in a new window] |
TABLE 3. Distribution of O-antigen genes with second primer set
|
|
View this table: [in a new window] |
TABLE 4. O-antigen genes tested by PCR
|
Two two-component systems that are likely to be involved in controlling resistance to heavy metals were found, including cusSR (30) (AB57_0660-1) and a less-well-characterized pair encoded by AB57_2550-1. The cusSR genes are absent from AB307-0294 and ACICU, along with the likely target genes involved in copper resistance. The cusSR genes are also missing from ADP1, while the likely target genes involved in copper resistance are present. In ATCC 17978, the regulatory and likely target genes are present adjacent to one another, but the gene cluster is located in a nonsyntenic position relative to flanking genes in AB0057. The predicted CusS and CusR proteins have their best (non-A. baumannii) matches to proteins from several betaproteobacteria, suggesting that they are derived from a lateral transfer event into an A. baumannii ancestor following the split from A. baylyi ADP1.
At least one component of 11 of the 19 two-component systems that are present in the clinical isolates is absent from SDF, suggesting a reduced capacity to respond to environmental changes in this organism.
|
|
|---|
Among 74 isolates from WRAMC, resistance to the aminoglycosides amikacin and tobramycin was always associated with the presence of at least one aminoglycoside-modifying enzyme (AME) (19). In addition, a recent study of A. baumannii isolates in New York City showed that resistance to cefepime and tigecycline could be ascribed to increased expression of AdeABC in some isolates, while resistance to other cephalosporins and aminoglycosides was mediated by antimicrobial drug-inactivating enzymes and fluoroquinolone resistance was always associated with gyrase gene mutations (4). It is possible that upregulation of efflux pumps may be found to be sufficient to confer resistance to certain classes of antibiotics in other isolates. For example, one WRAMC isolate is ciprofloxacin resistant, despite carrying gyrA and parC alleles associated with drug susceptibility (54, 55), suggesting an alternative mechanism of resistance. Taken together, the extensive similarities between AB307-0294 and the MDR isolates imply that while drug efflux may be a contributor to clinically relevant antibiotic resistance, it is not the primary mechanism of MDR.
The RI is a dominant genomic feature in AYE and AB0057, exhibiting a large concentration of laterally transferred genes associated with antimicrobial resistance intermixed with a diverse set of transposon- and integron-related sequences that have presumably facilitated the transfer of this region of the genome into A. baumannii. The importance of the RI in contributing to the MDR phenotype is considerable in these isolates, with genes encoding resistance to cephalosporins, carbapenems, and aminoglycosides. In each of the MDR isolates, however, additional loci contribute to drug resistance. IS-flanked β-lactamase genes are found outside the RIs in AB0057 and ACICU, and cephalosporin resistance in AYE is likely due to an ISAba1 insertion adjacent to the blaADC-9 gene. Furthermore, the RI found in AYE and AB0057, and to a lesser extent, ACICU, does not appear to be present in all MDR isolates based on PCR amplification of several genes located throughout the RI. Thus, it is likely that alternative genetic arrangements of resistance genes will be found as additional isolates are sequenced.
Quinolone resistance is primarily mediated by mutations in the housekeeping genes gyrA and parC, which are not RI associated. The bla genes (e.g., encoding TEM, OXA, VEB, and CTX-M) are frequently found in locations outside the RI, even in isolates that contain the RI. In ACICU, two copies of blaOXA-58 are present on a plasmid, and in AB0057, blaOXA-23 is located on the second island copy, ISAbaR4. In the MDR isolates analyzed to date, AMEs are associated with the RI. It remains to be determined whether AMEs are associated with RIs of novel structure in isolates that do not contain an RI analogous to those described in AYE, AB0057, and ACICU. In summary, evidence to support the notion that the RI is the sole, or even primary, genetic basis for clinically relevant MDR is still needed.
It is curious that the sequence of the >50-year-old isolate ATCC 17978 is not more divergent than those of other isolates, suggesting that although A. baumannii isolates differ by
2% at the sequence level, there may be considerable stability in each clade. AYE and AB0057 were isolated more than a decade after AB307-0294 yet have diverged <0.1% at the nucleotide level compared with the older isolate. SDF also differs from the other A. baumannii isolates by
2% at the DNA level, although it has also undergone a considerable reduction in gene number, presumably driven by the action of a large number of IS elements (51). It was argued that SDF was acquired by the louse directly from the bloodstream of an infected human in the course of a blood meal, since the louse gut is sterile (24). An alternative explanation is that SDF represents an A. baumannii lineage that is adapted for growth in the louse rather than in humans. Additional epidemiological surveys will be necessary to determine whether SDF-like A. baumannii isolates can be found in humans.
Antimicrobial resistance in A. baumannii arises by a combination of genetic changes and selective pressure. The following three classes of genetic alteration contribute to the acquisition of a resistance phenotype: lateral gene transfer, gene amplification (by duplication and/or increased expression), and mutation of genes or their promoters (resulting in loss of function, enhancement of function, or gain of novel function). Each of these mechanisms is represented in A. baumannii. Further studies that address the balance among these strategies, the genetic and genomic architectures that support them, the role of mobile genetic elements, and the pace of genetic change will be facilitated by the sequences presented here and those of additional isolates that represent further examples of MDR and susceptible isolates of clinical significance.
We thank the J. C. Venter Institute for making available the JCVI Annotation Service, which provided us with automatic annotation data and the manual annotation tool Manatee. We thank Nicholas Ornston for helpful discussions regarding the metabolic capacity of Acinetobacter species. We thank David Craft of WRAMC for providing A. baumannii isolates.
Published ahead of print on 17 October 2008. ![]()
Supplemental material for this article may be found at http://jb.asm.org/. ![]()
Present address: Scripps Healthcare, San Diego, CA. ![]()
|
|
|---|
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»