Previous Article | Next Article ![]()
Journal of Bacteriology, August 2006, p. 5722-5730, Vol. 188, No. 16
0021-9193/06/$08.00+0 doi:10.1128/JB.01950-05
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Department of Chemistry and Biomolecular Sciences, Macquarie University, Sydney, NSW 2109, Australia,1 Genome Atlantic, Dalhousie University, 5859 University Avenue, Halifax, Nova Scotia B3H 4H7, Canada,2 Department of Biological Sciences, Macquarie University, Sydney, NSW 2109, Australia,3 Department of Microbiology, Institute of Biology, Sølvgade 83H, DK 1307, Copenhagen K, Denmark4
Received 20 December 2005/ Accepted 8 June 2006
|
|
|---|
|
|
|---|
Class 1 integrons are generally embedded in mobile elements, including numerous plasmids and transposons. Consequently, unlike most other integron classes, they are amenable to lateral gene transfer (LGT) (22, 28) and, when present on a chromosome, they are presumed to have transposed there (27, 47). In contrast, "chromosomal" integrons are capable of recruiting a diverse assortment of novel mobile genes but are, nonetheless, confined to defined phylogenetic lineages (8, 17, 31, 46).
Antibiotic resistance as a problem arose soon after the onset of the clinical use of antibiotics. After that time, resistance genes were commonly observed to be present on mobile elements, such as plasmids and transposons; before that time, they appear to have been relatively rare on these elements (18). The antibiotic era also coincided with a time of rapid appearance of class 1 integrons carrying antibiotic resistance genes, with these integrons appearing in a number of independent locations in different transposons and plasmids (28, 48). The simplest explanation for this rapid appearance in diverse clinical isolates and on many different types of mobile elements is that the class 1 integrons were moving by transposition. If class 1 integrons were originally derived from a chromosomal integron, then association with transposition genes that facilitate their mobility presumably occurred in the preantibiotic era or very early after the onset of the antibiotic era.
Class 1 integrons are characterized by several features (Fig. 1). These include the presence of a 5' conserved segment (CS). This region includes the DNA integrase gene, intI1, and the attI1 site, the latter being the region at which gene cassettes are inserted (33). The outer boundary of this segment is defined by a 25-bp sequence (IRi), which is present as an inverted repeat with respect to another sequence (IRt) located downstream of genes inserted at attI1. IRi and IRt and sequence immediately adjacent within the transposon include identifiable transposase binding sites necessary for transposition (20). These inverted repeats therefore facilitate the movement of class 1 integrons by transposition and define the limits of the mobile unit that includes the class 1 integron in a structural sense. Transposition of a class 1 integron requires four transposition genes of the type found in Tn402 and its relatives (23). These transposition genes, when present, are located downstream of the genes inserted at attI1. Tn402 is an example of both an active transposon (19) and a class 1 integron (39). Transposition is relatively target specific, targeting plasmid and transposon resolution (res) sites (19, 22, 23, 36, 47).
![]() View larger version (17K): [in a new window] |
FIG. 1. General structure of extant class 1 integrons. The element at the top indicates the structure of a class 1 integron which is transposable as a consequence of possessing a complete Tn402-like transposition system. This is presumed to comprise the structure of the ancestral mobile element from which most class 1 integrons, as represented in the bottom diagram, are derived. The filled black rectangle indicates the attI1 site. qacE is a quaternary ammonium compound resistance determinant. Other features are described in the text.
|
Whether the linkage of site-specific recombination functions with Tn402 transposition functions was an ancient event or one that occurred in historical times, it is nonetheless the case that it has become a powerful vehicle for the spread of resistance genes. The site-specific recombination system has allowed individual class 1 integrons to capture a diverse array of resistance determinants. The Tn402-like transposition genes have, in turn, facilitated the incorporation of this element into a myriad of other mobile elements. These other elements can comprise plasmids into which the mobile class 1 integron has transposed directly or other transposons into which the class 1 integron/transposon has inserted (40). The latter situation, an "integron within a transposon within a transposon," is perhaps best exemplified by the capture of a class 1 integron with Tn402 transposition functions by the Tn21 family of transposons (28). As noted above, all class 1 integrons found to date possess at least some evidence of Tn402-like transposon features and/or sequences that make up the 3'-CS (34, 35, 39). Other types of transposition genes, such as those in Tn21, are not universally present. Consequently, the acquisition of such additional transposition functions had to occur through events postdating the linkage of the Tn402-like genes to the site-specific recombination system.
Here, we have screened forest soil and lake sediment environments removed from clinical settings for bacteria containing class 1 integrons. Sample culturing was done in a way that does not bias toward the recovery of antibiotic-resistant bacteria. All environments contained isolates that were positive for the intI1 gene. Class 1 integrons were characterized for four phylogenetically distinct bacteria. While they were identical or nearly identical to class 1 integrons from clinical settings, all lacked the Tn402-like transposition genes that are presumed to have contributed to the dissemination of this class into pathogenic and commensal bacteria in historical times. In addition, none included antibiotic resistance cassettes. The broad distribution of class 1 integrons in this environment, and the evidence for dispersal to numerous locations, demonstrates that extensive transfer of these integrons can occur across very disparate environments. This dispersal may have occurred via a process independent of that which brought about the dispersal of the class 1 integrons described up until now.
|
|
|---|
Strain recovery and characterization. Soil or shallow freshwater sediment was serially diluted in 100 mM sodium phosphate buffer (pH 7.0) and plated to plate count agar (PCA) medium (5 g tryptone, 2.5 g yeast extract, 1 g dextrose and 12 g agar per liter). Plates were incubated at 25°C for 5 days, after which 200 individual colonies from each environment were picked into 100 µl of PCA broth in sterile microtiter trays. Trays were incubated at 25°C for 48 h.
Crude DNA was prepared from all bacterial isolates (at least 180 for each environment) by harvesting 10 µl from each well into 30 µl sterile water, heating to 99°C for 10 min, and centrifuging to pellet cell debris. The 16S rRNA gene was also amplified for samples that tested positive for intI1 (24). The nitrogen fixation gene nifH was amplified from the DNA of isolate MUL2G9 using primers 19B and 407B (11). Amplification mixes consisted of 5 µl of template DNA, 200 nM dNTP, 50 pmol of each primer, 2 mM MgCl2, and 1 U of Red Hot DNA polymerase in the reaction buffer supplied with the enzyme. After initial thermal denaturation, PCR was performed for 30 cycles of 94°C for 1 min, 60°C for 45 s, and 72°C for 80 s. All putative gene products were sequenced using an ABI Prism 377 (PE Biosystems) and BigDye v3.1 chemistry. Table 1 displays a list of the PCR primers used to characterize and sequence class 1 integrons in this study.
|
View this table: [in a new window] |
TABLE 1. Primers for the amplification and sequencing of class 1 integrons
|
40 kb was cut out of the gel, purified from the agarose, ligated to the fosmid vector, packaged in phage capsids, and used to infect Escherichia coli as described in the CopyControl kit manual. Four hundred eighty colonies from each of the resulting libraries were picked from agar plates and used to inoculate 96-well blocks containing 0.5 ml of LB broth with 12.5 µg/ml of chloramphenicol (to select for the presence of fosmids) in each well. These cultures were grown overnight, and glycerol stocks were made by mixing 70 µl of culture from each well with 30 µl of 50% glycerol in a 96-well plate. Fifteen microliters was also sampled from each well and pooled by row in Eppendorf tubes (180 µl per tube). Tubes were centrifuged to remove the LB broth, and cells were resuspended in 25 µl of pH 7.0 Tris-EDTA buffer. One microliter of resuspended cells was added as template to PCRs using primers HS463A and HS464 to screen for the presence of intI1. All wells from a row of the block that was positive were then screened individually. Clones testing positive in the second round of screening were restreaked on LB agar plates containing 12.5 µg/ml of chloramphenicol and were then used to inoculate a liquid culture for extraction of pure fosmid DNA (according to the CopyControl kit manual). One to six positive clones were obtained for three of the four 480 clone libraries, and no positives were obtained for the Acidovorax sp. MUL2A2 library. Sequence assembly. The intI-containing fosmid clone from the Azoarcus communis MUL2G9 strain was subcloned using a TOPO Shotgun subcloning kit (Invitrogen). Subclones were sequenced to obtain eight times coverage of the complete insert. Fragments were assembled using phred, phrap, and consed (http://www.phrap.org/phredphrapconsed.html) (6, 7, 10), which were also used to assess the quality of the sequencing. This yielded a 33-kb contig that included the intI gene with 28 kb of flanking DNA on one side and another 4 kb of sequence on the other side. A set of PCR primers targeting the end of the contig was designed and used to rescreen the library. Both ends of library clones that were positive with this set of primers were sequenced to determine overlap with the original 33-kb contig. A fosmid clone with minimal overlap (2 to 5 kb) was selected to extend the sequence of the contig by walking, which provided an additional 3 kb of sequence.
The integrase-containing fragments found for the Acidovorax sp. MUL2G8 and Burkholderiales bacterium MUL2G11 fosmid libraries were not completely sequenced. Instead, their integrons were sequenced by walking away from the integrase gene in both directions. Sequence (
5 kb) was obtained for each strain, including the complete integron as well as some flanking DNA.
Sequence alignments and phylogenetic analysis.
16S rRNA gene sequences for the MUL2A2, MUL2G8, MUL2G9, and MUL2G11 isolates were used as BLASTN queries to retrieve closely related sequences from GenBank and identify their taxonomic affiliations. Representative 16S rRNA genes from the taxonomic groups to which the isolates belonged were retrieved from GenBank and aligned to the sequences generated in this study using CLUSTALW (50). The alignment was subsequently edited manually to remove ambiguous characters. The 16S rRNA gene trees were constructed with PAUP* 4.04b, applying the heuristic search option and using the tree bisection-reconnection branch-swapping algorithm. Maximum likelihood was used as the tree reconstruction method, with the nucleotide substitution model (general time reversible), gamma rate parameter
, proportion of invariable sites, and nucleotide frequencies determined using MODELTEST (37). The confidence of each node was determined by building a consensus tree of 100 maximum likelihood trees from bootstrap pseudoreplicates of the original data set.
Amino acid alignments of ORFs encoded in the intI1-containing genomic DNA fragments from the MUL2G8, MUL2G9, and MUL2G11 isolates were constructed using CLUSTALW and edited manually to remove ambiguous characters. Maximum likelihood phylogenetic analyses were applied to these datasets using PROML with the JTT amino acid substitution matrix, a rate heterogeneity model with gamma-distributed rates over four categories with the
parameter being estimated using TREE-PUZZLE, global rearrangements, and randomized input order of sequences (10 jumbles). Bootstrap support values represent a consensus (obtained using CONSENSE) of 100 Fitch-Margoliash distance trees (obtained using PUZZLEBOOT and FITCH) from pseudo-replicates (obtained using SEQBOOT) of the original alignment. The settings of PUZZLEBOOT were the same as those used for PROML, except that no global rearrangements and randomized input order of sequences are available with this program. PROML, CONSENSE, FITCH, and SEQBOOT are from the PHYLIP package version 3.6a, available at http://evolution.genetics.washington.edu/phylip.html. TREE-PUZZLE and PUZZLEBOOT can be obtained from the programs' website, http://www.tree-puzzle.de.
Nucleotide sequence accession numbers. The DNA sequences determined in this study have been submitted to the GenBank nucleotide sequence database and have been assigned accession numbers DQ372709 to DQ372715.
|
|
|---|
![]() View larger version (15K): [in a new window] |
FIG. 2. Best maximum likelihood trees for the 16S rRNA genes of the class 1 integron bearing isolates from Lake Yerbury sediments and their relatives. (A) Phylogeny of representative Burkholderiales and isolates MUL2A2, MUL2G8, and MUL2G11. (B) Phylogeny of representative Rhodocyclales and isolate MUL2G9.
|
Characterization of a class 1 integron in Azoarcus communis strain MUL2G9. The isolate most extensively sequenced was MUL2G9, with sequence including a region of approximately 28 kb beyond the outer boundary of the 5'-CS and approximately 5 kb beyond attI1. The integron was missing the first 98 bases of the 5'-CS (Fig. 3). This deletion end point is within one base pair of the insertion point of the insertion sequence ISPa7 present in some class 1 integrons isolated from Pseudomonas aeruginosa (1, 42, 49). In these cases, the insertion element is an insertion into a conventional class 1 integron structure, as the first 99 bases are still present. In contrast, in MUL2G9, the first 98 bases of the 5'-CS are absent. These first 98 bases (measured from the outer boundary of IRi) (Fig. 3) that are lost include all of IRi and the two strong TnpA binding domains necessary for transposition (20). The region immediately adjacent to the 5'-CS in MUL2G9 includes a putative DNA integrase/recombinase (Fig. 4) with another 18 identified ORFs following. Of these 18 ORFs, seven predicted products showed very high matches to proteins located in the sequenced genome of the Azoarcus sp. EbN1, the closest relative of A. communis for which substantial sequence information is available (38) (Fig. 2B). All seven corresponding genes are located on the chromosome of EbN1, and at least three are housekeeping genes, including a sigma-54 transcriptional activator, an acyl coenzyme A synthetase, and an acetone carboxylase beta subunit (see database entry DQ372711). We conclude that the class 1 integron identified here is in the chromosome of A. communis MUL2G9.
![]() View larger version (17K): [in a new window] |
FIG. 3. Left boundary of 5'-CS for MUL2G8, MUL2G9, and MUL2G11. The 25-base sequence comprising IRi in MUL2G11 and previously described class 1 integrons is shown in boldface type and identifies the outer boundary of the 5'-CS, as previously defined. Dots indicate positions of sequence identity with MUL2G11. The vertical arrow indicates the point of insertion of ISPa7 in some class 1 integrons. The underlined three bases are the stop codon (complementary strand) for intI1.
|
![]() View larger version (11K): [in a new window] |
FIG. 4. Structure of class 1 integrons and surrounding regions in MUL2G8, MUL2G9, and MUL2G11. Horizontal dashed lines indicate nonsequenced regions and are not to scale. Regions separated by horizontal dashed lines on the far left and right of diagram are derived by end sequencing using vector-specific primers. Vertical dashed lines delineate regions of DNA or amino acid identity at the levels indicated. Functions assigned to coding regions are putative, based on database matches as indicated in the text. 59-be, 59-base element; tran., trans-acting; acyl-CoA, acyl coenzyme A.
|
Characterization of other class 1 integrons from lake sediment isolates. Sequence derived from a MUL2G8 fosmid clone also confirmed a class 1 integron (Fig. 4). Further, the integron showed structural similarities to the element in MUL2G9. First, it was missing the first 98 bases of the 5'-CS, with the breakpoint at the same location as seen for MUL2G9 (Fig. 3). While the DNA sequence beyond the breakpoint was unrelated to the corresponding sequence in MUL2G9, it was nonetheless the case that the predicted product of the gene immediately adjacent to the 5'-CS was a DNA integrase/recombinase homologous to the one found at the same position in MUL2G9. The amino acid identity between the two proteins was 72% over a region spanning the 165 C-terminal amino acids of these predicted proteins (Fig. 4). Thus, although the two insertion events are independent, a common underlying mechanism of insertion is implied. Two ORFs were found immediately beyond attI1. The first of these is an ORF of unknown function that was followed by a putative transposase. The ORFs were not followed by an identifiable 59-base element, implying that neither gene was part of a gene cassette. The putative transposase was not closely related to the transposase adjacent to the cassette array in the MUL2G9 integron.
In addition to lacking Tn402 transposition functions, the MUL2G8 class 1 integron also lacked IRi and IRt. Overall, the sequence surrounding the MUL2G8 integron suggests that it is located in its host's chromosome. Sequence comparisons of the 5'-CS present revealed that it was identical to the equivalent sequence in a number of class 1 integrons that include In2 of Tn21. In contrast, within the region common to the MUL2G8 and MUL2G9 integrons, 13 differences could be identified. In MUL2G9, these differences are represented by 2 substitutions in attI1 and 11 in intI1 (6 nonsynonymous and 5 synonymous) compared to MUL2G8 and In2. The 5'-CS sequence seen in MUL2G9 is not seen in extant class 1 integrons and is relatively divergent in the context of this region. Given the difference between MUL2G8 and MUL2G9, these two independent capture events are likely to be derived from different parent integrons if the events occurred in historical time.
Sequence analysis of a MUL2G11 fosmid clone revealed a class 1 integron that includes a complete 5'-CS (Fig. 3). That is, the first 98 bases of this region were present. The MUL2G11 5'-CS was also identical to the 5'-CS in In2. The sequence immediately beyond IRi in MUL2G11included a DNA invertase/recombinase beyond which was DNA sequence that was closely related (>90%) to sequence found for several IncP1-ß plasmids (Fig. 4). It was therefore concluded that this integron is located on a plasmid. The sequence beyond the attI1 site was identical to the equivalent region in the MUL2G8 integron. This sequence identity extended at least through the transposase (tnpA) gene, at which point the MUL2G8 clone ends (Fig. 4). Thus, while the locations of the class 1 integrons present in MUL2G8 and MUL2G11 are likely to be different, unlike the MUL2G9 example, they would appear to share a common immediate ancestor. Additional information on the MUL2G11 clone was derived by end sequencing from vector sequence. About 600 bases were sequenced from the right end (Fig. 4) of this clone. Interestingly, this region was about 89% identical to a region of Tn402 that includes the end of tniA and the beginning of tniB (Fig. 5).
![]() View larger version (9K): [in a new window] |
FIG. 5. Best maximum likelihood tree for part of the tniAB genes from the Acidovorax sp. MUL2G8 isolate and various transposons. The numbers in parentheses are accession numbers for the tniAB gene sequences.
|
|
View this table: [in a new window] |
TABLE 2. Determination of the structure of class 1 integrons of bacterial isolates from Lake Yerbury sediments by PCR amplification with primers targeting specific regions
|
|
|
|---|
The features in common for the two types of integrons are an intI gene, an associated attI site, and a Pc, although the last feature has not been experimentally demonstrated in every case. These features collectively allow mobile gene cassettes to be captured by integrons and any cassette-associated genes to be expressed. Beyond these common features, however, the two integron types differ in several aspects. The nature of their respective cassette arrays is their most notable distinguishing feature. Integrons associated with other mobile elements almost exclusively possess gene cassettes that carry antibiotic resistance genes (40); collectively, resistance determinants to most of the clinically important antibiotics are found in class 1 integrons. A typical multicassette array in a class 1 integron contains less than six or so cassettes (40). Cassette arrays of chromosomal integrons, recognized as being generally confined to phylogenetic lines, have different characteristics. The first of these is that the genes within cassettes are very diverse (8, 17, 46, 51). For most, a function cannot be ascribed. Where a potentially attributable function is found, it is almost always unrelated to antibiotic resistance. A second feature of chromosomal cassette arrays is that they can be quite large, in excess of 200 cassettes in the case of some Vibrio species (29).
The distinction between integron types by the use of the descriptors "chromosomal" and "mobile" can be useful, but it does not alter the fact that both types are functionally equivalent in that they capture mobile gene cassettes by site-specific recombination. It has also been shown that cassette arrays containing antibiotic resistance genes can be found in chromosomal integrons and that such arrays within integrons can be a source of resistance genes for class 1 integrons (45).
Chromosomal integrons are very diverse and are broadly distributed phylogenetically, so it is clear that these elements have been in bacteria for a long period of evolutionary time (17, 31, 32, 46). One hypothesis to explain the origin of the mobilizable class 1 integron is that it is derived from a chromosomal integron that became mobile by the acquisition of Tn402-like transposition functions. This, together with the selection for antibiotic resistance as a consequence of human intervention, allowed the dispersal of class 1 integrons and thereby greatly contributed to the clinical problem that we see today.
Is the above hypothesis correct? Given the abundance and distribution of chromosomal integrons, the mobilization of such an integron by the addition of transposition genes is a possibility. However, we observed with at least two cases (MUL2G8 and MUL2G9) that class 1 integrons can be found devoid of Tn402-like transposition functions and need not be associated with antibiotic resistance (see Table 3) but can nonetheless be broadly dispersed among bacteria more likely to be closely associated with plants than animals. The same is nearly true for the remaining two (MUL2G11 and MUL2A2), with the exception that IRi and two TnpA binding sites are present. Our data can be interpreted in one of two ways. The first of these is that the class 1 integrons isolated here were derived from Tn402-like class 1 integrons that are beginning to radiate out of pathogenic and commensal bacteria. Implicit in this hypothesis is that this mobilization leads to the loss of most, or all, of the Tn402-like transposition functions and to the loss of antibiotic resistance genes. Further, the potential radiation of class 1 integrons out of pathogenic bacteria has probably occurred in historical times, that is, in the last 60 years or so, corresponding to the advent of the antibiotic era, since it was only at that time that they became prevalent. A second explanation for the observations made here is that class 1 integrons were already being mobilized prior to the acquisition of the Tn402-like transposition system and prior to the antibiotic era. Since the class 1 integrons found in the sampled lake sediment ecosystem have identical, or nearly identical, counterparts in multidrug resistant strains, it is not possible to say which way the lateral transfer events occurred. Intuitively, the second explanation is more likely. We have identified unusual class 1 integrons with respect to their lack of association with Tn402 transposition functions in at least four distinct bacteria. Also, the associated 5'-CS sequences are not all the same, suggesting either a common ancestor that predates the historical use on antibiotics or recent integration events that involve distinct immediate ancestors. It seems unlikely that multiple events involving simultaneous loss of the many features summarized in Table 3 could occur in historical times. Other anecdotal evidence also hints at the possibility that class 1 integrons may have been widespread in bacteria from sediment and soil environments from a time that predates the widespread use of antibiotics. In a metagenomic survey for integrons from the tailings of an abandoned gold mine in Colorado, 12 examples of class 1 integrons were identified, as determined by the sequencing of a region internal to the intI1 gene (32). It is not known how many of these samples were independent. In addition, we do not have sequence context information or know what the host organism(s) was. Although the possibility of recent anthropogenic disturbance cannot be ruled out, it is nonetheless an intriguing observation, since the mine itself has been abandoned for many decades (D. R. Nemergut, personal communication).
|
View this table: [in a new window] |
TABLE 3. Structural differences between previously identified class 1 integrons and those isolated in this study
|
Putting aside the question of whether class 1 integrons have moved into pathogenic and commensal bacteria or become dispersed from them, the widespread presence of this integron class in soil and sediment bacteria is remarkable. Although some degree of antibiotic contamination of the sampled ecosystem studied here is possible, it is nonetheless relatively remote from commensal populations, and it is certainly the case that the class 1-containing bacteria recovered are examples of neither commensal bacteria nor human pathogens. Nitrogen-fixing Azoarcus sp. strains, in particular, of which MUL2G9 is likely an example (it harbors a nifH gene), are best adapted to the highly specialized environment of the root interior of flood-tolerant grasses with large aerenchymatic air spaces in mature roots (41). The site from which MUL2G9 was recovered is consistent with this environment type in that it was sediment from a lake that had extensive reed growth. The environment type in which Azoarcus is found includes the rhizosphere and is an environment with a diverse microbial ecosystem and where the potential for LGT may be expected to be high. Indeed, the rhizosphere has been identified as an environment of prodigious rates of LGT (3); the sequenced Azoarcus strain EbN1, for example, shows considerable evidence of comprising an unusually plastic genome (38).
The fact that a single gene capture system is commonly present in commensals, pathogens, the rhizosphere, and aquatic environments further underscores the potential ease with which antibiotic resistance genes can transit from the common organisms that manufacture antibiotics to organisms detrimental to humans. In metagenomic studies, antibiotic resistance genes can be recovered from the general environment (43). The presence of class 1 integrons in disparate environmental niches and hosts effectively increases the pool of genes, including those encoding antibiotic resistance, that are accessible to these elements. Finally, it is noteworthy that one of the grasses most amenable to colonization by nitrogen-fixing Azoarcus species is rice (30). The presence of class 1 integrons in bacteria intimately associated with a crop production plant consumed worldwide arguably creates a major conduit for the frequent influx of new mobile genes into bacteria in close association with humans and other animals at a scale not previously appreciated.
We thank W. Ford Doolittle for helpful discussions.
|
|
|---|
plasmid pTB11 from a waste-water treatment plant reveals a highly conserved backbone, a Tn402-like integron and other transposable elements. Plasmid 53:218-238.[CrossRef][Medline]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»