When phage, plasmids, and transposons collide: genomic islands, and conjugative- and mobilizable-transposons as a mosaic continuum

Metagenomic analysis of mesopelagic Antarctic plankton reveals a novel deltaproteobacterial group

Phylogenetic screening of 3200 clones from a metagenomic library of Antarctic mesopelagic picoplankton allowed the identification of two bacterial 16S-rDNA-containing clones belonging to the Deltaproteobacteria, DeepAnt-1F12 and DeepAnt-32C6. These clones were very divergent, forming a monophyletic cluster with the environmental sequence GR-WP33-58 that branched at the base of the myxobacteria. Except for the possession of complete rrn operons without associated tRNA genes, DeepAnt-1F12 and DeepAnt-32C6 were very different in gene content and organization. Gene density was much higher in DeepAnt-32C6, whereas nearly one-third of DeepAnt-1F12 corresponded to intergenic regions. Many of the predicted genes encoded by these metagenomic clones were informational (i.e. involved in replication, transcription, translation and related processes). Despite this, a few putative cases of horizontal gene transfer were detected, including a transposase. DeepAnt-1F12 contained one putative gene encoding a long cysteine-rich protein, probably membrane-bound and Ca2+-binding, with only eukaryotic homologues. DeepAnt-32C6 carried some predicted genes involved in metabolic pathways that suggested this organism may be anaerobic and able to ferment and to degrade complex compounds extracellularly.

State of the knowledge of bacterial resistance

Bacteria have adapted a variety of different ways to acquire antibiotic resistance, fostering the rapid development of resistance within a short evolutionary time. The general genetic basis of events leading to and promoting antibiotic resistance formation in bacteria are presented and exemplified by showing the evolution of methicillin, glycopeptide, linezolid, and ketolide resistance in Staphylococcus aureus.

The naphthalene catabolic (nag) genes of Polaromonas naphthalenivorans CJ2: evolutionary implications for two gene clusters and novel regulatory control

Polaromonas naphthalenivorans CJ2, found to be responsible for the degradation of naphthalene in situ at a coal tar waste-contaminated site (C.-O. Jeon et al., Proc. Natl. Acad. Sci. USA 100:13591-13596, 2003), is able to grow on mineral salts agar media with naphthalene as the sole carbon source. Beginning from a 484-bp nagAc-like region, we used a genome walking strategy to sequence genes encoding the entire naphthalene degradation pathway andadditional flanking regions. We found that the naphthalene catabolic genes in P. naphthalenivorans CJ2 were divided into one large and one small gene cluster, separated by an unknown distance. The large gene cluster (nagRAaGHAbAcAdBFCQEDJI'ORF1tnpA) is bounded by a LysR-type regulator (nagR). The small cluster (nagR2ORF2I"KL) is bounded by a MarR-type regulator (nagR2). The catabolic genes of P. naphthalenivorans CJ2 were homologous to many of those of Ralstonia U2, which uses the gentisate pathway to convert naphthalene to central metabolites. However, three open reading frames (nagY, nagM, and nagN), present in Ralstonia U2, were absent. Also, P. naphthalenivorans carries two copies of gentisate dioxygenase (nagI) with 77.4% DNA sequence identity to one another and 82% amino acid identity to their homologue in Ralstonia sp. strain U2. Investigation of the operons using reverse transcription PCR showed that each cluster was controlled independently by its respective promoter. Insertional inactivation and lacZ reporter assays showed that nagR2 is a negative regulator and that expression of the small cluster is not induced by naphthalene, salicylate, or gentisate. Association of two putative Azoarcus-related transposases with the large cluster and one Azoarcus-related putative salicylate 5-hydroxylase gene (ORF2) in the small cluster suggests that mobile genetic elements were likely involved in creating the novel arrangement of catabolic and regulatory genes in P. naphthalenivorans.

Bacteroides fragilis mobilizable transposon Tn5520 requires a 71 base pair origin of transfer sequence and a single mobilization protein for relaxosome formation during conjugation

Tn5520 is the smallest known bacterial mobilizable transposon and was isolated from an antibiotic resistant Bacteroides fragilis clinical isolate. When a conjugation apparatus is provided in trans, Tn5520 is mobilized (transferred) efficiently within, and from, both Bacteroides spp. and Escherichia coli. Only two genes are present on Tn5520; one encodes an integrase, and the other a multifunctional mobilization (Mob) protein BmpH. BmpH is essential for Tn5520 mobility. The focus of this study was to identify the Tn5520 origin of conjugative transfer (oriT) and to study BmpH-oriT binding. We delimited the functional Tn5520 oriT to a 71 bp sequence upstream of the bmpH gene. A plasmid vector harbouring this minimal 71 bp oriT was mobilized at the same frequency as that of intact Tn5520. The minimal oriT contains one 17 bp inverted repeat (IR) sequence. We constructed and tested multiple IR mutants and showed that the IR was essential in its entirety for mobilization. A nick site sequence (5'-GCTAC-3') was also identified within the minimal oriT; this sequence resembled nick sites found in plasmids of Gram positive origin. We further showed that mutation of a highly conserved GC dinucleotide in the nick site sequence completely abolished mobilization. We also purified BmpH and showed that it specifically bound a Tn5520 oriT fragment in electrophoretic mobility shift assays. We also identified non-nick site sequences within the minimal oriT that were essential for mobilization. We hypothesize that transposon-based single Mob protein systems may contribute to efficient gene dissemination from Bacteroides spp., because fewer DNA processing proteins are required for relaxosome formation.

Assessing genetic diversity in plasmids from Escherichia coli and Salmonella enterica using a mixed-plasmid microarray

AIMS

To compare genetic composition of plasmids using microarrays composed of randomly selected fragments of plasmid DNA.

METHODS AND RESULTS

Separate shotgun libraries were constructed from plasmid DNA pooled from Escherichia coli and Salmonella enterica. Cloned fragments were used as probes for microarrays. Plasmid targets were labelled, hybridized overnight, and bound targets were imaged after enzymatic signal amplification. Control hybridizations demonstrated significantly higher signal when probes and targets shared >95% sequence identity. Diagnostic sensitivity and specificity of the assay was 95 and 99%, respectively. Cluster analysis showed close matches for replicate experiments with a high correlation between replicates (r = 0.91) compared with the correlation for nonreplicates (r = 0.09). Analysis of hybridization data from 43 plasmids generated five distinct clusters, two for known serovar-specific plasmids, one for enterohemorrhagic E. coli plasmids, and two for plasmids harboring a recently disseminated antibiotic resistance gene (bla(CMY-2)).

CONCLUSION

Mixed-plasmid microarrays are suitable for comparing genetic content of wild-type plasmids and hybridization results from this study suggest several novel hypotheses about plasmid gene exchange between E. coli and S. enterica.

SIGNIFICANCE AND IMPACT OF STUDY

Mixed-plasmid microarrays permit rapid, low cost analysis and comparison of many plasmids. This ability is critical to understanding the source, fate, and transport of plasmids amongst commensal and pathogenic bacteria.

A hypervariable 130-kilobase genomic region of Magnetospirillum gryphiswaldense comprises a magnetosome island which undergoes frequent rearrangements during stationary growth

Genes involved in magnetite biomineralization are clustered in the genome of the magnetotactic bacterium Magnetospirillum gryphiswaldense. We analyzed a 482-kb genomic fragment, in which we identified an approximately 130-kb region representing a putative genomic "magnetosome island" (MAI). In addition to all known magnetosome genes, the MAI contains genes putatively involved in magnetosome biomineralization and numerous genes with unknown functions, as well as pseudogenes, and it is particularly rich in insertion elements. Substantial sequence polymorphism of clones from different subcultures indicated that this region undergoes frequent rearrangements during serial subcultivation in the laboratory. Spontaneous mutants affected in magnetosome formation arise at a frequency of up to 10(-2) after prolonged storage of cells at 4 degrees C or exposure to oxidative stress. All nonmagnetic mutants exhibited extended and multiple deletions in the MAI and had lost either parts of or the entire mms and mam gene clusters encoding magnetosome proteins. The mutations were polymorphic with respect to the sites and extents of deletions, but all mutations were found to be associated with the loss of various copies of insertion elements, as revealed by Southern hybridization and PCR analysis. Insertions and deletions in the MAI were also found in different magnetosome-producing clones, indicating that parts of this region are not essential for the magnetic phenotype. Our data suggest that the genomic MAI undergoes frequent transposition events, which lead to subsequent deletion by homologous recombination under physiological stress conditions. This can be interpreted in terms of adaptation to physiological stress and might contribute to the genetic plasticity and mobilization of the magnetosome island.

Indications for acquisition of reductive dehalogenase genes through horizontal gene transfer by Dehalococcoides ethenogenes strain 195

The genome of Dehalococcoides ethenogenes strain 195, an anaerobic dehalorespiring bacterium, contains 18 copies of putative reductive dehalogenase genes, including the well-characterized tceA gene, whose gene product functions as the key enzyme in the environmentally important dehalorespiration process. The genome of D. ethenogenes was analyzed using a bioinformatic tool based on the frequency of oligonucleotides. The results in the form of a genomic signature revealed several local disruptions of the host signature along the genome sequence. These fractures represent DNA segments of potentially foreign origin, so-called atypical regions, which may have been acquired by an ancestor through horizontal gene transfer. Most interestingly, 15 of the 18 reductive dehalogenase genes, including the tceA gene, were found to be located in these regions, strongly indicating the foreign nature of the dehalorespiration activity. The GC content and the presence of recombinase genes within some of these regions corroborate this hypothesis. A hierarchical classification of the atypical regions containing the reductive dehalogenase genes indicated that these regions were probably acquired by several gene transfer events.

Elements in the co-evolution of relaxases and their origins of transfer

The central elements in the conjugative mobilization of most plasmids are the relaxase and its cognate origin of transfer (oriT). The relaxase of the plasmid R1162, together with its oriT, belong to a large and widely distributed family of related relaxase/oriT pairs. Several of the properties of these elements are considered for R1162 and for other members of this family with a view to understanding how systems for mobilization might have evolved.

The evolution of IncP catabolic plasmids

The recent adoption of whole plasmid genome sequencing as a routine analytical technique has provided the basis for cataloging the historical events through which plasmids are assembled from the available families of modular plasmid components. Horizontal gene transfer mediated by plasmids plays an important role in the adaptation of bacteria to the presence of specific metabolizable compounds, including man-made chemicals, in the surrounding environment. Bacterial plasmid genome sequence comparisons indicate that plasmids have complex genetic histories resulting from transposition, homologous recombination, and illegitimate recombinational events. Evidence from IncP plasmid genome sequences indicates that cryptic plasmid backbones acquire diverse catabolic pathways through gene capture and horizontal gene transfer.

TheSalmonellagenomic island 1 is an integrative mobilizable element

Salmonella genomic island 1 (SGI1) is a genomic island containing an antibiotic resistance gene cluster identified in several Salmonella enterica serovars. The SGI1 antibiotic resistance gene cluster, which is a complex class 1 integron, confers the common multidrug resistance phenotype of epidemic S. enterica Typhimurium DT104. The SGI1 occurrence in S. enterica serovars Typhimurium, Agona, Paratyphi B, Albany, Meleagridis and Newport indicates the horizontal transfer potential of SGI1. Here, we report that SGI1 could be conjugally transferred from S. enterica donor strains to non-SGI1 S. enterica and Escherichia coli recipient strains where it integrated into the recipient chromosome in a site-specific manner. First, an extrachromosomal circular form of SGI1 was identified by PCR which forms through a specific recombination of the left and right ends of the integrated SGI1. Chromosomal excision of SGI1 was found to require SGI1-encoded integrase which presents similarities to the lambdoid integrase family. Second, the conjugal transfer of SGI1 required the presence of a helper plasmid. The conjugative IncC plasmid R55 could thus mobilize in trans SGI1 which was transferred from the donor to the recipient strains. By this way, the conjugal transfer of SGI1 occurred at a frequency of 10(-5)-10(-6) transconjugants per donor. No transconjugants could be obtained for the SGI1 donor lacking the int integrase gene. Third, chromosomal integration of SGI1 occurred via a site-specific recombination between a 18 bp sequence found in the circular form of SGI1 and a similar 18 bp sequence at the 3' end of thdF gene in the S. enterica and E. coli chromosome. SGI1 appeared to be transmissible only in the presence of additional conjugative functions provided in trans. SGI1 can thus be classified within the group of integrative mobilizable elements (IMEs).

The complete nucleotide sequence of the resistance plasmid R478: defining the backbone components of incompatibility group H conjugative plasmids through comparative genomics

Horizontal transfer of resistance determinants amongst bacteria can be achieved by conjugative plasmid DNA elements. We have determined the complete 274,762 bp sequence of the incompatibility group H (IncH) plasmid R478, originally isolated from the Gram negative opportunistic pathogen Serratia marcescens. This self-transferable extrachromosomal genetic element contains 295 predicted genes, of which 144 are highly similar to coding sequences of IncH plasmids R27 and pHCM1. The regions of similarity among these three IncH plasmids principally encode core plasmid determinants (i.e., replication, partitioning and stability, and conjugative transfer) and we conducted a comparative analysis to define the minimal IncHI plasmid backbone determinants. No resistance determinants are included in the backbone and most of the sequences unique to R478 were contained in a large contiguous region between the two transfer regions. These findings indicate that plasmid evolution occurs through gene acquisition/loss predominantly in regions outside of the core determinants. Furthermore, a modular evolution for R478 was signified by the presence of gene neighbors or operons that were highly related to sequences from a wide range of chromosomal, transposon, and plasmid elements. The conjugative transfer regions are most similar to sequences encoded on SXT, Rts1, pCAR1, R391, and pRS241d. The dual partitioning modules encoded on R478 resemble numerous sequences; including pMT1, pCTX-M3, pCP301, P1, P7, and pB171. R478 also codes for resistance to tetracycline (Tn10), chloramphenicol (cat), kanamycin (aphA), mercury (similar to Tn21), silver (similar to pMG101), copper (similar to pRJ1004), arsenic (similar to pYV), and tellurite (two separate regions similar to IncHI2 ter determinants and IncP kla determinants). Other R478-encoded sequences are related to Tn7, IS26, tus, mucAB, and hok, where the latter is surrounded by insLKJ, and could potentially be involved in post-segregation killing. The similarity to a diverse set of bacterial sequences highlights the ability of horizontally transferable DNA elements to acquire and disseminate genetic traits through the bacterial gene pool.

Comprehensive de novo structure prediction in a systems-biology context for the archaea Halobacterium sp. NRC-1

BACKGROUND

Large fractions of all fully sequenced genomes code for proteins of unknown function. Annotating these proteins of unknown function remains a critical bottleneck for systems biology and is crucial to understanding the biological relevance of genome-wide changes in mRNA and protein expression, protein-protein and protein-DNA interactions. The work reported here demonstrates that de novo structure prediction is now a viable option for providing general function information for many proteins of unknown function.

RESULTS

We have used Rosetta de novo structure prediction to predict three-dimensional structures for 1,185 proteins and protein domains (<150 residues in length) found in Halobacterium NRC-1, a widely studied halophilic archaeon. Predicted structures were searched against the Protein Data Bank to identify fold similarities and extrapolate putative functions. They were analyzed in the context of a predicted association network composed of several sources of functional associations such as: predicted protein interactions, predicted operons, phylogenetic profile similarity and domain fusion. To illustrate this approach, we highlight three cases where our combined procedure has provided novel insights into our understanding of chemotaxis, possible prophage remnants in Halobacterium NRC-1 and archaeal transcriptional regulators.

CONCLUSIONS

Simultaneous analysis of the association network, coordinated mRNA level changes in microarray experiments and genome-wide structure prediction has allowed us to glean significant biological insights into the roles of several Halobacterium NRC-1 proteins of previously unknown function, and significantly reduce the number of proteins encoded in the genome of this haloarchaeon for which no annotation is available.

Evolution of genomic islands by deletion and tandem accretion by site-specific recombination: ICESt1-related elements from Streptococcus thermophilus

The 34 734-bp integrative and potentially conjugative element (putative ICE) ICESt1 has been previously found to be site-specifically integrated in the 3' end of the fda locus of Streptococcus thermophilus CNRZ368. Four types of genomic islands related to ICESt1 are integrated in the same position in seven other strains of S. thermophilus. One of these elements, ICESt3, harbours conjugation and recombination modules closely related to those of ICESt1 and excises by site-specific recombination. Two other types of elements, CIME19258 and CIME302, are flanked by site-specific attachment sites closely related to attL and attR of ICESt1 and ICESt3, whereas Delta CIME308 only possesses a putative attR site; none of these three elements carry complete conjugation and recombination modules. ICESt1 contains a functional internal recombination site, attL', that is almost identical to attL of CIME19258. The recombination between attL' and attR of ICESt1 leads to the excision of the expected circular molecule (putative ICE); a cis-mobilizable element (CIME) flanked by an attL site and an attB' site remains integrated into the 3' end of fda. Furthermore, sequences that could be truncated att sites were found within ICESt1, ICESt3 and CIME302. All together, these data suggest that these genomic islands evolved by deletion and tandem accretion of ICEs and CIMEs resulting from site-specific recombination. A model for this evolution is proposed and its application to other genomic islands is discussed.

Genomic diversity of Burkholderia pseudomallei clinical isolates: subtractive hybridization reveals a Burkholderia mallei-specific prophage in B. pseudomallei 1026b

Burkholderia pseudomallei is the etiologic agent of the disease melioidosis and is a category B biological threat agent. The genomic sequence of B. pseudomallei K96243 was recently determined, but little is known about the overall genetic diversity of this species. Suppression subtractive hybridization was employed to assess the genetic variability between two distinct clinical isolates of B. pseudomallei, 1026b and K96243. Numerous mobile genetic elements, including a temperate bacteriophage designated phi1026b, were identified among the 1026b-specific suppression subtractive hybridization products. Bacteriophage phi1026b was spontaneously produced by 1026b, and it had a restricted host range, infecting only Burkholderia mallei. It possessed a noncontractile tail, an isometric head, and a linear 54,865-bp genome. The mosaic nature of the phi1026b genome was revealed by comparison with bacteriophage phiE125, a B. mallei-specific bacteriophage produced by Burkholderia thailandensis. The phi1026b genes for DNA packaging, tail morphogenesis, host lysis, integration, and DNA replication were nearly identical to the corresponding genes in phiE125. On the other hand, phi1026b genes involved in head morphogenesis were similar to head morphogenesis genes encoded by Pseudomonas putida and Pseudomonas aeruginosa bacteriophages. Consistent with this observation, immunogold electron microscopy demonstrated that polyclonal antiserum against phiE125 reacted with the tail of phi1026b but not with the head. The results presented here suggest that B. pseudomallei strains are genetically heterogeneous and that bacteriophages are major contributors to the genomic diversity of this species. The bacteriophage characterized in this study may be a useful diagnostic tool for differentiating B. pseudomallei and B. mallei, two closely related biological threat agents.

A new Bacteroides conjugative transposon that carries an ermB gene

The erythromycin resistance gene ermB has been found in a variety of gram-positive bacteria. This gene has also been found in Bacteroides species but only in six recently isolated strains; thus, the gene seems to have entered this genus only recently. One of the six Bacteroides ermB-containing isolates, WH207, could transfer ermB to Bacteroides thetaiotaomicron strain BT4001 by conjugation. WH207 was identified as a Bacteroides uniformis strain based on the sequence of its 16S rRNA gene. Results of pulsed-field gel electrophoresis experiments demonstrated that the transferring element was normally integrated into the Bacteroides chromosome. The element was estimated from pulsed-field gel data to be about 100 kb in size. Since the element appeared to be a conjugative transposon (CTn), it was designated CTnBST. CTnBST was able to mobilize coresident plasmids and the circular form of the mobilizable transposon NBU1 to Bacteroides and Escherichia coli recipients. A 13-kb segment that contained ermB was cloned and sequenced. Most of the open reading frames in this region had little similarity at the amino acid sequence level to any proteins in the sequence databases, but a 1,723-bp DNA segment that included a 950-bp segment downstream of ermB had a DNA sequence that was virtually identical to that of a segment of DNA found previously in a Clostridium perfringens strain. This finding, together with the finding that ermB is located on a CTn, supports the hypothesis that CTnBST could have entered Bacteroides from some other genus, possibly from gram-positive bacteria. Moreover, this finding supports the hypothesis that many transmissible antibiotic resistance genes in Bacteroides are carried on CTns.

Horizontal gene transfer and microbial adaptation to xenobiotics: new types of mobile genetic elements and lessons from ecological studies

The characterization of bacteria that degrade organic xenobiotics has revealed that they can adapt to these compounds by expressing 'novel' catabolic pathways. At least some of them appear to have evolved by patchwork assembly of horizontally transmitted genes and subsequent mutations and gene rearrangements. Recent studies have revealed the existence of new types of xenobiotic catabolic mobile genetic elements, such as catabolic genomic islands, which integrate into the chromosome after transfer. The significance of horizontal gene transfer and patchwork assembly for bacterial adaptation to pollutants under real environmental conditions remains uncertain, but recent publications suggest that these processes do occur in a polluted environment.

Structural comparison of the integrative and conjugative elements R391, pMERPH, R997, and SXT

R391 and SXT are members of a group of eleven chromosome-borne conjugative elements found in the gamma-proteobacteria, whose members carry different antibiotic resistance traits. Recent genomic analysis of R391 and SXT revealed a highly conserved 'backbone' encoding integration/excision, conjugative transfer, and regulation functions, augmented by an array of phenotypic traits and transposable elements. In this study, PCR amplification and sequence analysis were employed to investigate the genomic structure of two further MGE of the R391 family, pMERPH (HgR) and R997 (ApR, SmR, SuR). R997 and pMERPH were found to be structurally related to R391 and SXT and share a number of virtually identical regions with them-including putative integration, conjugative transfer, and regulatory determinants-interrupted by variable DNA segments and transposable elements. The presence of a highly conserved backbone in the four elements strongly suggests their origin in a common ancestral element, which itself was a mosaic of sequences related to phages and plasmids. Subsequent genetic recombination and the acquisition of transposable elements resulted in the possession of variable phenotypic traits among the four MGE, and diversification into two distinct lineages, the first one including R391 and pMERPH, the second one containing SXT and R997.

Sequence analysis of the mobile genome island pKLC102 of Pseudomonas aeruginosa C

The Pseudomonas aeruginosa plasmid pKLC102 coexists as a plasmid and a genome island in clone C strains. Whereas the related plasmid pKLK106 reversibly recombines with P. aeruginosa clone K chromosomes at one of the two tRNA(Lys) genes, pKLC102 is incorporated into the tRNA(Lys) gene only close to the pilA locus. Targeting of the other tRNA(Lys) copy in the chromosome is blocked by a 23,395-bp mosaic of truncated PAO open reading frames, transposons, and pKLC102 homologs. Annotation and phylogenetic analysis of the large 103,532-bp pKLC102 sequence revealed that pKLC102 is a hybrid of plasmid and phage origin. The plasmid lineage conferred oriV and genes for replication, partitioning, and conjugation, including a pil cluster encoding type IV thin sex pili and an 8,524-bp chvB glucan synthetase gene that is known to be a major determinant for host tropism and virulence. The phage lineage conferred integrase, att, and a syntenic set of conserved hypothetical genes also observed in the tRNA(Gly)-associated genome islands of P. aeruginosa clone C chromosomes. In subgroup C isolates from patients with cystic fibrosis, pKLC102 was irreversibly fixed into the chromosome by the insertion of the large 23,061-bp class I transposon TNCP23, which is a composite of plasmid, integron, and IS6100 elements. Intramolecular transposition of a copy of IS6100 led to chromosomal inversions and disruption of plasmid synteny. The case of pKLC102 in P. aeruginosa clone C documents the intraclonal evolution of a genome island from a mobile ancestor via a reversibly integrated state to irreversible incorporation and dissipation in the chromosome.

The role of mobile genetic elements in bacterial adaptation to xenobiotic organic compounds

Retrospective studies clearly indicate that mobile genetic elements (MGEs) play a major role in the in situ spread and even de novo construction of catabolic pathways in bacteria, allowing bacterial communities to rapidly adapt to new xenobiotics. The construction of novel pathways seems to occur by an assembly process that involves horizontal gene transfer: different appropriate genes or gene modules that encode different parts of the novel pathway are recruited from phylogenetically related or distant hosts into one single host. Direct evidence for the importance of catabolic MGEs in bacterial adaptation to xenobiotics stems from observed correlations between catabolic gene transfer and accelerated biodegradation in several habitats and from studies that monitor catabolic MGEs in polluted sites.