Physical mapping of bacterial genomes

Reference-free Structural Variant Detection in Microbiomes via Long-read Coassembly Graphs

Bacterial genome dynamics are vital for understanding the mechanisms underlying microbial adaptation, growth, and their broader impact on host phenotype. Structural variants (SVs), genomic alterations of 10 base pairs or more, play a pivotal role in driving evolutionary processes and maintaining genomic heterogeneity within bacterial populations. While SV detection in isolate genomes is relatively straightforward, metagenomes present broader challenges due to absence of clear reference genomes and presence of mixed strains. In response, our proposed method rhea, forgoes reference genomes and metagenome-assembled genomes (MAGs) by encompassing a single metagenome coassembly graph constructed from all samples in a series. The log fold change in graph coverage between subsequent samples is then calculated to call SVs that are thriving or declining throughout the series. We show rhea to outperform existing methods for SV and horizontal gene transfer (HGT) detection in two simulated mock metagenomes, which is particularly noticeable as the simulated reads diverge from reference genomes and an increase in strain diversity is incorporated. We additionally demonstrate use cases for rhea on series metagenomic data of environmental and fermented food microbiomes to detect specific sequence alterations between subsequent time and temperature samples, suggesting host advantage. Our innovative approach leverages raw read patterns rather than references or MAGs to include all sequencing reads in analysis, and thus provide versatility in studying SVs across diverse and poorly characterized microbial communities for more comprehensive insights into microbial genome dynamics.

Bacterial genome sequencing

For over 30 yr, the Sanger method has been the standard for DNA sequencing. Instruments have been developed and improved over time to increase throughput, but they always relied on the same technology. Today, we are facing a revolution in DNA sequencing with many drastically different platforms that have become or will soon become available on the market. We review a number of sequencing technologies and provide examples of applications. We also discuss the impact genomics and new DNA sequencing approaches have had on various fields of biological research.

Molecular characterisation of the metallo-beta-lactamase genes in imipenem-resistant Gram-negative bacteria from a university hospital in southern Taiwan

In this study, 260 non-replicate imipenem-resistant Gram-negative bacteria isolated between January 2002 and December 2006 were subjected to a screening test for detection of metallo-beta-lactamase (MBL) using the Etest containing imipenem and ethylene diamine tetra-acetic acid (EDTA). Among the 260 strains, 123 (47.3%) appeared to produce MBL. Of these 123 strains, 113 (91.9%) were found by polymerase chain reaction (PCR) to carry MBL genes of types blaVIM-2, blaVIM-3, blaVIM-11 (blaVIM-11a), blaIMP-8 and novel blaIMP-24. One strain of Serratia marcescens harboured two MBL genes (blaVIM-11 and blaIMP-8) simultaneously. Of the 123 strains, 116 strains (94.3%) carrying the intI1 gene and 21 strains carrying integron-associated blaVIM-3, blaVIM-11 and blaIMP-8 genes were identified among Acinetobacter baumannii, Pseudomonas aeruginosa, Acinetobacter haemolyticus and S. marcescens. Using pulsed-field gel electrophoresis (PFGE) and Southern hybridisation with the blaVIM gene probe for I-CeuI-digested genomic DNA, P. aeruginosa 9527 strain harboured two class 1 integron-associated MBL genes in the chromosome, including blaVIM-3-orf2-aacA4 and novel bla(VIM-3)-orf2-aacA4-aadB-aacA4. This is the first description of the blaVIM-11 gene spreading among P. aeruginosa and A. baumannii strains in southern Taiwan. This finding suggests that clinical spread of this blaVIM-11 gene is a matter of great concern for carbapenem resistance in southern Taiwan.

Diversity of genome structure in Salmonella enterica serovar Typhi populations

The genomes of most strains of Salmonella and Escherichia coli are highly conserved. In contrast, all 136 wild-type strains of Salmonella enterica serovar Typhi analyzed by partial digestion with I-CeuI (an endonuclease which cuts within the rrn operons) and pulsed-field gel electrophoresis and by PCR have rearrangements due to homologous recombination between the rrn operons leading to inversions and translocations. Recombination between rrn operons in culture is known to be equally frequent in S. enterica serovar Typhi and S. enterica serovar Typhimurium; thus, the recombinants in S. enterica serovar Typhi, but not those in S. enterica serovar Typhimurium, are able to survive in nature. However, even in S. enterica serovar Typhi the need for genome balance and the need for gene dosage impose limits on rearrangements. Of 100 strains of genome types 1 to 6, 72 were only 25.5 kb off genome balance (the relative lengths of the replichores during bidirectional replication from oriC to the termination of replication [Ter]), while 28 strains were less balanced (41 kb off balance), indicating that the survival of the best-balanced strains was greater. In addition, the need for appropriate gene dosage apparently selected against rearrangements which moved genes from their accustomed distance from oriC. Although rearrangements involving the seven rrn operons are very common in S. enterica serovar Typhi, other duplicated regions, such as the 25 IS200 elements, are very rarely involved in rearrangements. Large deletions and insertions in the genome are uncommon, except for deletions of Salmonella pathogenicity island 7 (usually 134 kb) from fragment I-CeuI-G and 40-kb insertions, possibly a prophage, in fragment I-CeuI-E. The phage types were determined, and the origins of the phage types appeared to be independent of the origins of the genome types.

Physical map of the linear chromosome of Streptomyces hygroscopicus 10-22 deduced by analysis of overlapping large chromosomal deletions

The chromosomal DNA of Streptomyces hygroscopicus 10-22, a derivative of strain 5102-6, was digested with several restriction endonucleases and analyzed by pulsed-field gel electrophoresis (PFGE). Digestions with AseI gave 11 fragments with a total length of ca. 7.36 Mb. The AseI sites were mapped by analysis of overlapping chromosomal deletions in different mutants and confirmed by Southern hybridizations using partially digested genome fragments and linking cosmids as probes. PFGE analysis of DNA with and without proteinase K treatment, together with the hybridization results, suggested a linear organization with terminal proteins and large terminal inverted repeats. Some deletion mutants had circular chromosomes.

A physical map of the syringomycin and syringopeptin gene clusters localized to an approximately 145-kb DNA region of Pseudomonas syringae pv. syringae strain B301D

Genetic and phenotypic mapping of an approximately 145-kb DraI fragment of Pseudomonas syringae pv. syringae strain B301D determined that the syringomycin (syr) and syringopeptin (syp) gene clusters are localized to this fragment. The syr and syp gene clusters encompass approximately 55 kb and approximately 80 kb, respectively. Both phytotoxins are synthesized by a thiotemplate mechanism of biosynthesis, requiring large multienzymatic proteins called peptide synthetases. Genes encoding peptide synthetases were identified within the syr and syp gene clusters, accounting for 90% of the DraI fragment. In addition, genes encoding regulatory and secretion proteins were localized to the DraI fragment. In particular, the salA gene, encoding a regulatory element responsible for syringomycin production and lesion formation in P. syringae pv. syringae strain B728a, was localized to the syr gene cluster. A putative ATP-binding cassette (ABC) transporter homolog was determined to be physically located in the syp gene cluster, but phenotypically affects production of both phytotoxins. Preliminary size estimates of the syr and syp gene clusters indicate that they represent two of the largest nonribosomal peptide synthetase gene clusters. Together, the syr and syp gene clusters encompass approximately 135 kb of DNA and may represent a genomic island in P. syringae pv. syringae that contributes to virulence in plant hosts.

General method of rapid Smith/Birnstiel mapping adds for gap closure in shotgun microbial genome sequencing projects: application to Pseudomonas putida KT2440

A physical mapping strategy has been developed to verify and accelerate the assembly and gap closure phase of a microbial genome shotgun-sequencing project. The protocol was worked out during the ongoing Pseudomonas putida KT2440 genome project. A macro-restriction map was constructed by linking probe hybridisation of SwaI- or I-CeuI-restricted chromosomes to serve as a backbone for the quick quality control of sequence and contig assemblies. The library of PCR-generated SwaI linking probes was derived from the sequence assembly after 3- and 6-fold genome coverage. In order to support gap closure in regions with ambiguous assemblies such as the repetitive sequence of the seven ribosomal operons, high-resolution Smith/Birnstiel maps were generated by Southern hybridisation of pulsed-field gel electrophoresis-separated rare-cutter complete/frequent-cutter partial digestions with rare-cutter fragment end probes. Overall 1.5 Mb of the 6.1 Mb P.putida KT2440 genome has been subjected to high-resolution physical mapping in order to align assemblies generated from shotgun sequencing.

A physical map of the megaplasmid pHG1, one of three genomic replicons in Ralstonia eutropha H16

We have used pulsed field gel electrophoresis and megabase DNA techniques to investigate the basic genomic organization of Ralstonia eutropha H16, and to construct a physical map of its indigenous megaplasmid pHG1. This Gram-negative, soil-dwelling bacterium is a facultative chemolithoautotroph and a denitrifier. In the absence of organic substrates it can grow on H2 as its sole energy source and CO2 as its sole source of carbon. Under anaerobic conditions it can utilize nitrate as a terminal electron acceptor, whereby dinitrogen is released. Essential genetic determinants of the enzyme systems responsible for these metabolic processes are linked to the 0.44-Mb conjugative megaplasmid pHG1. Aside from pHG1, the genome of R. eutropha H16 is comprised of two circular chromosomes measuring 4.1 and 2.9 Mb, adding up to a total genome size of 7.1 Mb. An estimated five copies of rDNA are distributed on the two chromosomes. A macrorestriction map of pHG1 was derived for the endonucleases DraI and XbaI. Hybridization studies showed that genes for anaerobic metabolism are located on all three genomic replicons.

Genome plasticity among related ++Lactococcus strains: identification of genetic events associated with macrorestriction polymorphisms

The genomic diversity of nine strains of the Lactococcus lactis subsp. cremoris (NCDO712, NCDO505, NCDO2031, NCDO763, MMS36, C2, LM0230, LM2301, and MG1363) was studied by macrorestriction enzyme analysis using pulsed-field gel electrophoresis. These strains were considered adequate for the investigation of genomic plasticity because they have been described as belonging to the same genetic lineage. Comparison of ApaI and SmaI genome fingerprints of each strain revealed the presence of several macrorestriction fragment length polymorphisms (RFLPs), despite a high degree of similarity of the generated restriction patterns. The physical map of the MG1363 chromosome was used to establish a genome map of the other strains and allocate the RFLPs to five regions. Southern hybridization analysis correlated the polymorphic regions with genetic events such as chromosomal inversion, integration of prophage DNA, and location of the transposon-like structures carrying conjugative factor or oligopeptide transport system.

A physical map of the chromosome of Xanthomonas campestris pv. phaseoli var. fuscans BXPF65

A physical map of the Xanthomonas campestris pv. phaseoli var. fuscans BXPF65 chromosome was constructed by PFGE and Southern hybridization. The chromosome was 3938 kb and was composed of 39 XbaI restriction fragments. The size is close to that estimated for X. campestris pv. campestris based on DNA renaturation kinetics. Macrorestriction fragments containing genes for a 16S rRNA, pectate lyase, two-component regulatory system and hrp cluster were located on the physical map. Each of these probes hybridized to single macrorestriction fragments except the 16S rRNA probe, which hybridized to five fragments located in different areas on the chromosome map.

The origin of eukaryotes: the difference between prokaryotic and eukaryotic cells

Eukaryotes have long been thought to have arisen by evolving a nucleus, endomembrane, and cytoskeleton. In contrast, it was recently proposed that the first complex cells, which were actually proto-eukaryotes, arose simultaneously with the acquisition of mitochondria. This so-called symbiotic association hypothesis states that eukaryotes emerged when some ancient anaerobic archaebacteria (hosts) engulfed respiring alpha-proteobacteria (symbionts), which evolved into the first energy-producing organelles. Therefore, the intracellular compartmentalization of the energy-converting metabolism that was bound originally to the plasma membrane appears to be the key innovation towards eukaryotic genome and cellular organization. The novel energy metabolism made it possible for the nucleotide synthetic apparatus of cells to be no longer limited by subsaturation with substrates and catalytic components. As a consequence, a considerable increase has occurred in the size and complexity of eukaryotic genomes, providing the genetic basis for most of the further evolutionary changes in cellular complexity. On the other hand, the active uptake of exogenous DNA, which is general in bacteria, was no longer essential in the genome organization of eukaryotes. The mitochondrion-driven scenario for the first eukaryotes explains the chimera-like composition of eukaryotic genomes as well as the metabolic and cellular organization of eukaryotes.

Physical and genetic map of the genome of Vibrio parahaemolyticus: presence of two chromosomes in Vibrio species

We constructed a physical map of the genomic DNA (5.1 Mb) for Vibrio parahaemolyticus strain AQ4673 by combining 17 adjacent NotI fragments. This map shows two circular replicons of 3.2 and 1.9 Mb. Pulsed-field gel electrophoresis (PFGE) of undigested genomic DNA revealed two bands of corresponding sizes. Analysis both by NotI digestion and by Southern blot of the two isolated bands confirmed the existence of two replicons. The presence of genes for 16S rRNA on both the replicons indicates that the replicons are chromosomes rather than megaplasmids. The two bands were also seen after PFGE of undigested genomic DNA of V. parahaemolyticus strains other than AQ4673, and of strains belonging to other Vibrio species, such as V. vulnificus, V. fluvialis and various serovars and biovars of V. cholerae. It is noteworthy that V. cholerae O1 strain 569B, a classical biovar, was also shown to have two replicons of 2.9 and 1.2 Mb, which does not agree with a physical map proposed in a previous study. Our results suggest that a two-replicon structure is common throughout Vibrio species.

The diverse and dynamic structure of bacterial genomes

Bacterial genome sizes, which range from 500 to 10,000 kbp, are within the current scope of operation of large-scale nucleotide sequence determination facilities. To date, 8 complete bacterial genomes have been sequenced, and at least 40 more will be completed in the near future. Such projects give wonderfully detailed information concerning the structure of the organism's genes and the overall organization of the sequenced genomes. It will be very important to put this incredible wealth of detail into a larger biological picture: How does this information apply to the genomes of related genera, related species, or even other individuals from the same species? Recent advances in pulsed-field gel electrophoretic technology have facilitated the construction of complete and accurate physical maps of bacterial chromosomes, and the many maps constructed in the past decade have revealed unexpected and substantial differences in genome size and organization even among closely related bacteria. This review focuses on this recently appreciated plasticity in structure of bacterial genomes, and diversity in genome size, replicon geometry, and chromosome number are discussed at inter- and intraspecies levels.

Combined physical and genetic map of the Pseudomonas putida KT2440 chromosome

A combined physical and genetic map of the Pseudomonas putida KT2440 genome was constructed from data obtained by pulsed-field gel electrophoresis techniques (PFGE) and Southern hybridization. Circular genome size was estimated at 6.0 Mb by adding the sizes of 19 SwaI, 9 PmeI, 6 PacI, and 6 I-CeuI fragments. A complete physical map was achieved by combining the results of (i) analysis of PFGE of the DNA fragments resulting from digestion of the whole genome with PmeI, SwaI, I-CeuI, and PacI as well as double digestion with combinations of these enzymes and (ii) Southern hybridization analysis of the whole wild-type genome digested with different enzymes and hybridized against a series of probes obtained as cloned genes from different pseudomonads of rRNA group I and Escherichia coli, as P. putida DNA obtained by PCR amplification based on sequences deposited at the GenBank database, and by labeling of macrorestriction fragments of the P. putida genome eluted from agarose gels. As an alternative, 10 random mini-Tn5-Km mutants of P. putida KT2440 were used as a source of DNA, and the band carrying the mini-Tn5 in each mutant was identified after PFGE of a series of complete chromosomal digestions and hybridization with the kanamycin resistance gene of the mini-Tn5 as a probe. We established a circular genome map with an average resolution of 160 kb. Among the 63 genes located on the genetic map were key markers such as oriC, 6 rrn loci (rnnA to -F), recA, ftsZ, rpoS, rpoD, rpoN, and gyrB; auxotrophic markers; and catabolic genes for the metabolism of aromatic compounds. The genetic map of P. putida KT2440 was compared to those of Pseudomonas aeruginosa PAO1 and Pseudomonas fluorescens SBW25. The chromosomal backbone revealed some similarity in gene clustering among the three pseudomonads but differences in physical organization, probably as a result of intraspecific rearrangements.

Construction of a physical and preliminary genetic map of Aeromonas hydrophila JMP636

A physical and preliminary genetic map of the Aeromonas hydrophila JMP636 chromosome has been constructed. The topology of the genome was predicted to be circular as chromosomal DNA did not migrate from the origin during PFGE unless linearized by S1 nuclease. Cleavage of the chromosome with PacI and PmeI produced 23 and 14 fragments, respectively, and enabled calculation of the genome size at 4.5 Mb. Digestion of the chromosome with I-CeuI produced 10 fragments, indicating that 10 rrl (23S) genes were likely to be present. Hybridizations between DNA fragments generated with PacI, PmeI and I-CeuI were used to initially determine the relationship between these segments. To accurately map genes previously characterized from JMP636, the suicide vector pJP5603 was modified to introduce restriction sites for PacI and PmeI, producing pJP9540. Following cloning of genes into this vector and recombinational insertion into the JMP636 chromosome, PacI and PmeI cleavage determined the location of genes within macrorestriction fragments with the additional bands produced forming hybridization probes. From the data generated, it was possible to form a physical map comprising all the fragments produced by PacI and PmeI, and assign the contig of I-CeuI fragments on this map. The preliminary genetic map defines the location of six loci for degradative enzymes previously characterized from JMP636, while the locations of the 10 sets of ribosomal genes were assigned with less accuracy from hybridization data.

Physical and genetic map of the Pasteurella multocida A:1 chromosome

A physical and genetic map of the Pasteurella multocida A:1 genome was generated by using the restriction enzymes ApaI, CeuI, and NotI. The positions of 23 restriction sites and 32 genes, including 5 rrn operons, were localized on the 2.35-Mbp single circular chromosome. This report presents the first genetic and physical map for this genus.

Physical and genetic chromosomal maps of Streptococcus agalactiae, serotypes II and III; rRNA operon organization

A detailed analysis of two Streptococcus agalactiae (group B streptococcus, GBS) strains was performed by pulsed field gel electrophoresis (PFGE). Digestion of the chromosomal DNA with SmaI and SgrAI endonucleases, followed by separation and analysis of fragments by PFGE was carried out. Physical chromosomal maps of serotype II/(alpha + beta) and III/alpha strains of S. agalactiae were constructed. The GBS genome size was estimated to be 2200 kb. Sixteen GBS genes were used as probes and were located on the restriction maps of both strains by DNA-DNA hybridization. Six copies of ribosomal operons were found in the genome of the analyzed strains. Significant differences in the restriction patterns of chromosomal DNA and DNA-DNA hybridization between the two analyzed strains were detected so that DNA restriction patterns may be used to trace outbreaks of disease. The overall GBS chromosomal organization as determined is fairly conserved.

A natural large chromosomal inversion in Lactococcus lactis is mediated by homologous recombination between two insertion sequences

Comparative analysis of chromosomal macrorestriction polymorphism of the two closely related Lactococcus lactis subsp. cremoris strains MG1363 and NCDO763 revealed the presence of a large inversion covering half of the genome. To determine what kind of genetic element could be implicated in this rearrangement, the two inversion junctions of MG1363 and NCDO763 chromosomes were cloned and characterized. Nucleotide sequence analysis showed the presence of one copy of the lactococcal IS905 element in each junction. Each copy of this element contained the same nucleotide mutation that inactivates the putative transposase. Comparison of the sequences surrounding the insertion sequence demonstrated that the large inversion arose from a single-step homologous recombination event between the two defective copies of the IS905 element. The large inversion presumably conferred no selective disadvantage on strain NCDO763 because this rearrangement did not alter the oriC-terC symmetry of the chromosome and the local genetic environment.

Bacterial artificial chromosome (BAC) library as a tool for physical mapping of the archaeon Methanosarcina thermophila TM-1

We have used a variety of methods to characterize the genome of the archaeon Methanosarcina thermophila TM-1. Pulsed-field gel analysis indicates a genome size of 2.8 Mb. We have constructed a bacterial artificial chromosome (BAC) library of M. thermophila and have used it to generate physical maps for this organism. The library is made up of 384 clones with an average insert size of 58 kb representing 8.0 genome equivalents. The utility of the library for low-resolution physical mapping was shown by identifying NotI linking clones and using these to order the NotI macrorestriction fragments of M. thermophila into a 2.8 Mb map. Hybridization of nine single copy genes and a 16S rRNA sequence to these macrorestriction fragments forms the basis for the first genetic map in this organism. High-resolution physical maps, consisting of overlapping clones, have been created using HindIII fingerprints of BAC clones. In this way, we identified a minimal path of five clones that span a 270 kb NotI fragment. The ease of manipulating BAC clones makes the BAC system an excellent choice for the construction of low-resolution and high-resolution physical and genetic maps of archaeal genomes. It also provides a substrate for future genome-sequencing efforts.

Mycobacterial genome structure

Genome maps have been constructed for the mycobacterial pathogens Mycobacterium leprae and Mycobacterium tuberculosis, as well as for the attenuated vaccine strain Mycobacterium bovis BCG Pasteur. While the chromosomes of M. tuberculosis and M. bovis BCG Pasteur show extensive conservation at the gross level, comparison with M. leprae revealed a high degree of diversification, with a mosaic-like pattern apparent. The ordered libraries of M. tuberculosis and M. leprae produced during the course of these studies played a central role in the genome sequencing projects of these two bacilli, showing the utility of this approach for systematic sequencing of bacterial genomes.

Genome encyclopedias and their use for comparative analysis of Rhodobacter capsulatus strains

This paper consists of two components: the use of gene encyclopedias in genomic studies and Rhodobacter capsulatus genome project. A survey of vectors used for encyclopedia construction includes a brief discussion of their relative advantages and limitations. Projects employing various methods of encyclopedia assembly including the comparison of restriction patterns, restriction maps, linking by hybridization, oligonucleotide fingerprinting, sequence tagged site (STS) fingerprinting and encyclopedias derived from genetic maps are listed and briefly described. The R. capsulatus SB 1003 genome project started with the construction of its cosmid encyclopedia, which comprises 192 cosmids covering the chromosome and the 134 kbp plasmid in strain SB 1003, with the exact map coordinates of each cosmid. In a pilot sequencing study, several cosmids were individually subcloned using the vector M13mp18 and merged into one 189 kbp contig. About 160 open reading frames (ORFs) identified by the CodonUse program were subjected to similarity searches. The biological functions of eighty ORFs could be assigned reliably using the WIT (what is there) genome investigation environment. Eighty percent of these recognizable ORFs were organized in functional clusters, which simplified assignment decisions and increased the strength of the predictions. A set of 26 genes for cobalamin biosynthesis, genes for polyhydroxyalkanoic acid metabolism, DNA replication and recombination, and DNA gyrase were among those identified. Recently, another 1.2 Mbp genome fragment of the Rhodobacter genome was sequenced using a slightly modified approach. These results together with some genome investigation tools, have been placed at our web site (http://capsulapedia.uchicago.edu). The sequence of R. capsulatus is expected to be completed by summer 1998. A project to construct a systematic set of deletion strains of R. capsulatus in order to assign functions to unknown ORFs has been started. Preliminary data demonstrate the extreme convenience of the unique gene transfer agent (GTA) system to perform such work.

Cloning of prokaryotic genomes in yeast artificial chromosomes: application to the population genetics of Pseudomonas aeruginosa

Yeast artificial chromosomes (YACs) can accommodate large inserts and hence should be attractive tools for intra- and interspecies comparisons of bacterial genomes. YAC libraries were constructed from size-selected partial digests of human and Pseudomonas aeruginosa PAO DNA and SpeI-restricted PAO DNA. Whereas YACs from human DNA had an average size of 350 kilobase pairs (kbp), a P. aeruginosa sequence larger than 120 kbp was absent or truncated in the eukaryotic host. Coligation occurred for YACs smaller than 40 kbp, but stable YACs with 40-120 kbp large inserts of P. aeruginosa DNA were obtained in high yield. SpeI-restricted chromosomes from 97 P. aeruginosa strains representing 47 genotypes were hybridized with stable YACs from three equidistant regions of the PAO genome. The low complexity of hybridizing bands demonstrated that the analyzed 100 kbp sequence contigs were stably maintained in most P. aeruginosa isolates from both disease and environmental habitats.

Chromosomal rearrangements in enteric bacteria

Early genetic studies showed conservation of gene order in the enteric bacteria. Two recent methods using pulsed field gel electrophoresis (PFGE) to determine the physical map of the genome are: (i) partial digestion with the endonuclease I-CeuI, which digests the DNA of bacteria in the rrn operon for rRNA (ribosomal RNA), thus establishing the "rrn genomic skeleton" (the size in kbp of the intervals between rRNA operons); (ii) analysis of XbaI and B1nI sites within Tn10 insertions in the chromosome. The order of I-CeuI fragments, which is ABCDEFG in S. typhimurium LT2 and E. coli K-12, was found to be conserved in most Salmonella species, most of which grow in many hosts (host-generalists). However, in S. typhi, S. paratyphi C, S. gallinarum, and S. pullorum, species which are host-specialized, these fragments are rearranged, due to homologous recombination between the rrn operons, resulting in translocations and inversions. Inversions and translocations not involving the rrn operons are seldom detected except for inversions over the TER (termination of replication) region. Additive genetic changes (due to lateral transfer resulting in insertion of nonhomologous DNA) have resulted in "loops" containing blocks of DNA which provide new genes to specific strains, thus driving rapid evolution of new traits.

Leptospira genomics

The bacterial species Leptospira interrogans (sensu stricto) has a complex genome containing two circular chromosomal replicons. Comparative analysis of the larger chromosome reveals a fluid genetic organization with many large rearrangements differentiating two closely related strains. In the present study new genes were identified by partial sequence analysis of randomly cloned fragments of L. interrogans DNA. These genes were localized in regions of the genome by nucleic acid hybridization with DNA fragments separated by pulsed-field gel electrophoresis. The resulting genetic maps provide improved resolution for each strain and provide evidence for additional chromosomal rearrangements. Insertion elements may be involved in recombination events, as several are near regions of the chromosome that have undergone rearrangement.

Differential genome analysis of bacteria by genomic subtractive hybridization and pulsed field gel electrophoresis

A comprehensive analysis of the differences between the genomes of two closely related bacterial strains should give insight into the molecular basis of their individual phenotypic and genotypic characteristics. Here we present an integrative approach including two different strategies for the thorough investigation of genomic divergence. We have combined two techniques including genomic subtractive hybridization and comparative genome mapping by pulsed field gel electrophoresis (PFGE) techniques. The subtractive method for which a protocol is given herein results in the production of a library of specific DNA sequence tags present only in one strain, while the construction of macrorestriction maps of the bacterial chromosomes yields data about the overall genome organization and the arrangement and distance of gene loci. Comparison of the physical and genetic maps and determination of the map positions of the strain-specific DNA sequences reveals gross chromosomal modifications, insertions or deletions of additional genetic material, and transpositional events. The further investigation of the strain-specific regions yields information about the nature and origin of the acquired DNA and their influence on the evolution of the individual bacterial genome. The two methods were applied to differential genome analysis of clonal divergence in Pseudomonas aeruginosa choosing two clone C isolates from diverse habitats.

Chromosome sizes and phylogenetic relationships between serotypes of Actinobacillus pleuropneumoniae

The genome size of Actinobacillus pleuropneumoniae was determined by pulsed field gel electrophoresis of AscI and ApaI digested chromosomal DNA. The genome size of the type strain 4074T (serotype 1) was determined to be 2404 +/- 40 kb. The chromosome sizes for the reference strains of the other serotypes range between 2.3 and 2.4 Mb. The restriction pattern profiles of AscI, ApaI and NheI digested chromosomes showed a high degree of polymorphism among the different serotype reference strains and allowed their discrimination. The analysis of the macrorestriction pattern polymorphism revealed phylogenetic relationships between the different serotype reference strains which reflect to some extent groups of serotypes known to cross-react serologically. In addition, different pulsed fields gel electrophoresis patterns also revealed heterogeneity in the chromosomal structure among different field strains of serotypes 1, 5a, and 5b, while strains of serotype 9 originating from most distant geographical places showed homogeneous ApaI patterns in pulsed field gel electrophoresis.

Ordered cloned DNA map of the genome of Vibrio cholerae 569B and localization of genetic markers

By using a low-resolution macrorestriction map as the foundation (R. Majumder et al., J. Bacteriol. 176:1105-1112, 1996), an ordered cloned DNA map of the 3.2-Mb chromosome of the hypertoxinogenic strain 569B of Vibrio cholerae has been constructed. A cosmid library the size of about 4,000 clones containing more than 120 Mb of V. cholerae genomic DNA (40-genome equivalent) was generated. By combining landmark analysis and chromosome walking, the cosmid clones were assembled into 13 contigs covering about 90% of the V. cholerae genome. A total of 92 cosmid clones were assigned to the genome and to regions defined by NotI, SfiI, and CeuI macrorestriction maps. Twenty-seven cloned genes, 9 rrn operons, and 10 copies of a repetitive DNA sequence (IS1004) have been positioned on the ordered cloned DNA map.

Physical and genetic map of the genome of Staphylococcus carnosus TM300

A genome map of Staphylococcus carnosus TM300, an important micro-organism in the food industry and long used as a starter culture, was constructed by pulsed-field gel electrophoresis of DNA fragments obtained after digestion with NotI, SfiI and ApaI. The size of the chromosome was estimated to be 2590 kb. The fragments were assembled into a physical map using a combination of complementary methods including multiple and partial digests of genomic DNA, hybridization with homologous gene probes, and cross-Southern hybridization. Fifteen genes or gene clusters were positioned on the physical map by Southern hybridization analysis. The map provides a basis for further analysis of the S. carnosus chromosome.

Macrorestriction analysis of Desulfurella acetivorans and Desulfurella multipotens

The genomes of the phylogenetically and physiologically unique bacteria Desulfurella acetivorans DSM 5264T and D. multipotens DSM 8415T were characterized and compared by pulsed field gel electrophoresis (PFGE). Macrorestriction patterns made of large PFGE separated DNA fragments were generated by digesting the genomic DNAs of both strains with the rare cutting restriction endonucleases ApaI, AscI, EagI, RsrII, SacII, SalI as well as with the intron encoded endonuclease I-CeuI. The sum of calculated fragment sizes from digests of the first six enzymes yielded estimates for the chromosome sizes of D. acetivorans with a mean of 1939.0 +/- 26.0 kb and for D. multipotens with a mean of 1864.0 +/- 23.0 kb. Within the patterns obtained from EagI and RsrII cleavages the apparent differences could be attributed to DNA insertion or deletion and to point mutation. The single, circular chromosomes of the two strains contain two copies of 23S rRNA genes each. Different extrachromosomal elements were detected in both strains.

Differences in chromosome number and genome rearrangements in the genus Brucella

We have studied the genomic structure and constructed the SpeI, PacI and I-CeuI restriction maps of the four biovars of the pathogenic bacterium Brucella suis. B. suis biovar 1 has two chromosomes of 2.1 Mb and 1.15 Mb, similar to those of the other Brucella species: B. melitensis, B. abortus, B. ovis and B. neotomae. Two chromosomes were also observed in the genome of B. suis biovars 2 and 4, but with sizes of 1.85 Mb and 1.35 Mb, whereas only one chromosome with a size of 3.1 Mb was found in B. suis biovar 3. We show that the differences in chromosome size and number can be explained by rearrangements at chromosomal regions containing the three rrn genes. The location and orientation of these genes confirmed that these rearrangements are due to homologous recombination at the rrn loci. This observation allows us to propose a scheme for the evolution of the genus Brucella in which the two chromosome-containing strains can emerge from an hypothetical ancestor with a single chromosome, which is probably similar to that of B. suis biovar 3. As the genus Brucella is certainly monospecific, this is the first time that differences in chromosome number have been observed in strains of the same bacterial species.

Low-resolution sequencing of Rhodobacter sphaeroides 2.4.1T: chromosome II is a true chromosome

The photosynthetic bacterium Rhodobacter sphaeroides 2.4.1T has two chromosomes, CI (approximately 3.0 Mb) and CII (approximately 0.9 Mb). In this study a low-redundancy sequencing strategy was adopted to analyse 23 out of 47 cosmids from an ordered CII library. The sum of the lengths of these 23 cosmid inserts was approximately 495 kb, which comprised approximately 417 kb of unique DNA. A total of 1145 sequencing runs was carried out, with each run generating 559 +/- 268 bases of sequence to give approximately 640 kb of total sequence. After editing, approximately 2.8% bases per run were estimated to be ambiguous. After the removal of vector and Escherichia coli sequences, the remaining approximately 565 kb of R. sphaeroides sequences were assembled, generating approximately 291 kb of unique sequences. BLASTX analysis of these unique sequences suggested that approximately 131 kb (45% of the unique sequence) had matches to either known genes, or database ORFs of hypothetical or unknown function (dORFs). A total of 144 strong matches to the database was found; 101 of these matches represented genes encoding a wide variety of functions, e.g. amino acid biosynthesis, photosynthesis, nutrient transport, and various regulatory functions. Two rRNA operons (rrnB and rrnC) and five tRNAs were also identified. The remaining 160 kb of DNA sequence which did not yield database matches was then analysed using CODONPREFERENCE from the GCG package. This analysis suggested that 122 kb (42% of the total unique DNA sequence) could encode putative ORFs (pORFs), with the remaining 38 kb (13%) possibly representing non-coding intergenic DNA. From the data so far obtained, CII does not appear to be specialized for encoding any particular metabolic function, physiological state or growth condition. These data suggest that CII contains genes which are functionally as diverse as those found on any other bacterial chromosome and also contains sequences (pORFs), which may prove to be unique to this organism.

Construction of the temperature-sensitive vectors pLUCH80 and pLUCH88 for delivery of Tn917::NotI/SmaI and use of these vectors to derive a circular map of Listeria monocytogenes Scott A, a serotype 4b isolate

A physical map of Listeria monocytogenes Scott A was generated by the pulsed-field technique of contour-clamped-homogeneous-electric-field (CHEF) electrophoresis. The circular genome of this serotype 4b strain contains 12 AscI fragments (38 to 790 kb), 5 NotI fragments (55 to 1,400 kb), 3 SrfI fragments (110, 1,110, and 2,000 kb), and 2 SfiI fragments (1,320 and 1,920 kb). Summation of individually sized fragments derived by digestion of Scott A genomic DNA with each of these four enzymes provided an average estimated genome length of 3,210 +/- 60 kb. Efforts to assemble the macrorestriction map benefited greatly from the construction and use of pLUCH80 and pLUCH88, temperature-sensitive vectors for delivering transposon Tn917::NotI/SmaI to the chromosome of Scott A. As another component of this study, the positions of four known virulence genes (inlA, mpl, hly, and prf) and three L. monocytogenes-specific sequences (lisM44, lisM51, and lisM52) were localized on the physical map of Scott A by hybridization. Probes prepared from lisM44, lisM51, and the four virulence genes hybridized within a cluster on a 150-kb fragment of the Scott A genome that overlaps part of the NotI-B and AscI-D fragments. The lisM52 probe hybridized with the AscI-F2 (120-kb) fragment of Scott A, which is separated from the NotI-B-AscI-D region by about 300 kb. These results established the first physical and genetic map of a serotype 4b strain of L. monocytogenes and provided further insight on this important food-borne pathogen at the genome level.

Large chromosomal inversions occur in Pseudomonas aeruginosa clone C strains isolated from cystic fibrosis patients

An analysis of the chromosomal structure of 21 Pseudomonas aeruginosa clone C genotypes of various origins was carried out. The circular chromosomes showed various insertions and deletions which did not alter the gene order in 6 environmental and one ear isolate in comparison with strain PAO and strain C. However, the chromosome structure differed by at least one large inversion in 50% of the genotypes (7/14) isolated from cystic fibrosis patients. Recombination endpoints were scattered throughout the chromosome. All but one inversion included the terminus of replication which was displaced by the recombination events. The rearrangements shifted the positions of the rrn operons with respect to the origin of replication, but the transcription of the rrn operons remained in the same direction as replication. The exclusive detection of inversions in isolates from cystic fibrosis lungs, which represent an atypical habitat for P. aeruginosa, supports the theory that peculiar features of this new ecological niche may select, cause or tolerate genomic changes.

Directed introduction of DNA cleavage sites to produce a high-resolution genetic and physical map of the Acinetobacter sp. strain ADP1 (BD413UE) chromosome

The natural transformability of the soil bacterium Acinetobacter sp. ADP1 (BD413UE), formerly classified as A. calcoaceticus, has facilitated previous physiological and biochemical investigations. In the present studies, the natural transformation system was exploited to generated a physical and genetic map of this strain's 3780 +/- 191 kbp circular chromosome. Previously isolated Acinetobacter genes were modified in vitro to incorporate a recognition sequence for the restriction endonuclease NotI. Following transformation of the wild-type strain by the modified DNA, homologous recombination placed each engineered NotI cleavage site at the chromosomal location of the corresponding gene. This allowed precise gene localization and orientation of more than 40 genes relative to a physical map which was constructed with transverse alternating field electrophoresis (TAFE) and Southern hybridization methods. The positions of NotI, AscI and I-CeuI recognition sites were determined, and the latter enzyme identified the presence of seven ribosomal RNA operons. Multiple chromosomal copies of insertion sequence IS1236 were indicated by hybridization. Several of these copies were concentrated in one region of the chromosome in which a spontaneous deletion of approximately 100 kbp occurred. Moreover, contrary to previous reports, ColE1-based plasmids appeared to replicate autonomously in Acinetobacter sp. ADP1.

Physical and genetic map of the chromosome of Dichelobacter nodosus strain A198

A physical map of the chromosome of Dichelobacter nodosus strain A198 was constructed using the restriction endonucleases EagI and StuI. Mapping data indicated the presence of a single, circular chromosome of 1.54 Mb. The three rRNA operons and the virulence related locus (vrl) were precisely positioned at the junctions of EagI and StuI fragments, and their transcriptional orientations were also determined. Other D. nodosus genes were assigned to specific EagI and StuI fragments. Analysis of the resultant map revealed that the putative virulence genes were not clustered on the chromosome which suggests that the D. nodosus virulence determinants have been acquired gradually and that virulence in D. nodosus is an evolving trait.

Comparative typing of Pseudomonas aeruginosa by random amplification of polymorphic DNA or pulsed-field gel electrophoresis of DNA macrorestriction fragments

Eighty-seven strains of Pseudomonas aeruginosa were typed by random amplification of polymorphic DNA (RAPD) and pulsed-field gel electrophoresis (PFGE) of macrorestriction fragments. Stains were clustered on the basis of interpretative criteria as presented previously for the PFGE analysis. Clusters of strains were also defined on the basis of epidemiological data and subsequently reanalyzed by RAPD. It was found that in an RAPD assay employing the enterobacterial repetitive intergenic consensus sequence ERIC2 as a primer, single band differences can be ignored; in this case, clonally related strains could be grouped as effectively and reliably as with PFGE. These data could be corroborated by the use of other primer species. However, some primers either showed reduced resolution or, in contrast, identified DNA polymorphisms beyond epidemiologically and PFGE-defined limits. Apparently, different primers define different windows of genetic variation. It is suggested that criteria for interpretation of the ERIC2 PCR fingerprints can be simple and straightforward: when single band differences are ignored, RAPD-determined grouping of P. aeruginosa is congruent with that obtained by PFGE. Consequently, this implies that RAPD can be used with trust as a first screen in epidemiological characterization of P. aeruginosa. The ability to measure the rate of molecular evolution of the P. aeruginosa genome clearly depends on the choice of restriction enzyme or primer when RAPD or PFGE, respectively, is applied for the detection of DNA polymorphisms.

Highly plastic chromosomal organization in Salmonella typhi

Gene order in the chromosomes of Escherichia coli K-12 and Salmonella typhimurium LT2, and in many other species of Salmonella, is strongly conserved, even though the genera diverged about 160 million years ago. However, partial digestion of chromosomal DNA of Salmonella typhi, the causal organism of typhoid fever, with the endonuclease I-CeuI followed by separation of the DNA fragments by pulsed-field gel electrophoresis showed that the chromosomes of independent wild-type isolates of S. typhi are rearranged due to homologous recombination between the seven rrn genes that code for ribosomal RNA. The order of genes within the I-CeuI fragments is largely conserved, but the order of the fragments on the chromosome is rearranged. Twenty-one different orders of the I-CeuI fragments were detected among the 127 wild-type strains we examined. Duplications and deletions were not found, but transpositions and inversions were common. Transpositions of I-CeuI fragments into sites that do not change their distance from the origin of replication (oriC) are frequently detected among the wild-type strains, but transpositions that move the fragments much further from oriC were rare. This supports the gene dosage hypothesis that genes at different distances from oriC have different gene dosages and, hence, different gene expression, and that during evolution genes become adapted to their specific location; thus, cells with changes in gene location due to transpositions may be less fit. Therefore, gene dosage may be one of the forces that conserves gene order, although its effects seem less strong in S. typhi than in other enteric bacteria. However, both the gene dosage and the genomic balance hypotheses, the latter of which states that the origin (oriC) and terminus (TER) of replication must be separated by 180 degrees C, need further investigation.

Physical mapping of the Borrelia miyamotoi HT31 chromosome in comparison with that of Borrelia turicatae, an etiological agent of tick-borne relapsing fever

We report the construction of physical maps of chromosomes for Borrelia miyamotoi HT31 (a new species of Borrelia) and Borrelia turicatae (relapsing fever agent) by pulsed-field gel electrophoresis of DNA fragments generated by digestion of chromosomal DNA with rare-cutting restriction endonucleases and reciprocal hybridization. The size of the B. miyamotoi HT31 chromosome was calculated to be approximately 925 kilobase pairs, and the chromosome for B. turicatae was estimated to be 951 kilobase pairs. The chromosomes of B. miyamotoi HT31 and B. turicatae consisted of single linear molecules. The locations of several genes have been assigned to the chromosome maps by Southern hybridization by using specific gene probes. Comparison of the genetic maps of the two species of Borrelia provided evidence that the gene order on the chromosomes is quite similar to that of Borrelia burgdorferi sensu lato strains and is highly conserved in the genus Borrelia.

Physical and genetic chromosomal map of an M type 1 strain of Streptococcus pyogenes

A physical map of the chromosome of an M type 1 strain of Streptococcus pyogenes was constructed following digestion with three different restriction enzymes, SmaI, SfiI, and SgrAI, and separation and analysis of fragments by pulsed-field gel electrophoresis. The genome size of this strain was estimated to be 1,920 kb. By employing Southern hybridization and PCR analysis, 36 genes were located on the map.

Comparative genomic analysis of the Haloferax volcanii DS2 and Halobacterium salinarium GRB contig maps reveals extensive rearrangement

Anonymous probes from the genome of Halobacterium salinarium GRB and 12 gene probes were hybridized to the cosmid clones representing the chromosome and plasmids of Halobacterium salinarium GRB and Haloferax volcanii DS2. The order of and pairwise distances between 35 loci uniquely cross-hybridizing to both chromosomes were analyzed in a search for conservation. No conservation between the genomes could be detected at the 15-kbp resolution used in this study. We found distinct sets of low-copy-number repeated sequences in the chromosome and plasmids of Halobacterium salinarium GRB, indicating some degree of partitioning between these replicons. We propose alternative courses for the evolution of the haloarchaeal genome: (i) that the majority of genomic differences that exist between genera came about at the inception of this group or (ii) that the differences have accumulated over the lifetime of the lineage. The strengths and limitations of investigating these models through comparative genomic studies are discussed.

An integrated map of the genome of the tubercle bacillus, Mycobacterium tuberculosis H37Rv, and comparison with Mycobacterium leprae

An integrated map of the genome of the tubercle bacillus, Mycobacterium tuberculosis, was constructed by using a twin-pronged approach. Pulsed-field gel electrophoretic analysis enabled cleavage sites for Asn I and Dra I to be positioned on the 4.4-Mb circular chromosome, while, in parallel, clones from two cosmid libraries were ordered into contigs by means of fingerprinting and hybridization mapping. The resultant contig map was readily correlated with the physical map of the genome via the landmarked restriction sites. Over 165 genes and markers were localized on the integrated map, thus enabling comparisons with the leprosy bacillus, Mycobacterium leprae, to be undertaken. Mycobacterial genomes appear to have evolved as mosaic structures since extended segments with conserved gene order and organization are interspersed with different flanking regions. Repetitive sequences and insertion elements are highly abundant in M. tuberculosis, but the distribution of IS6110 is apparently nonrandom.

Physical map of the genome of Planctomyces limnophilus, a representative of the phylogenetically distinct planctomycete lineage

A physical map of the chromosome of Planctomyces limnophilus DSM 3776T was constructed by pulsed-field gel electrophoresis techniques. A total of 32 cleavage sites for the rare-cutting restriction endonucleases PacI, PmeI, and SwaI were located on the chromosome, which was shown to be circular and approximately 5.2 Mbp in size. An extrachromosomal element was detected but was found not to be cleaved by any of the enzymes used in the analysis of the chromosome. The order of the fragments on the chromosome was determined by hybridization of excised, labelled restriction fragments to Southern blots of pulsed-field gel electrophoresis-separated restriction digests. Seven genetic markers, rrs, rrl, atpD, tuf, gyrB, rpoD, and dnaK, on the chromosome were located by hybridization. Probes for all genetic markers were obtained by PCR. For five of these markers, probes were constructed by PCR with degenerate primers targeting conserved sequences. The arrangement of the genetic markers was compared with that found in other bacteria.

Alignment of a 1.2-Mb chromosomal region from three strains of Rhodobacter capsulatus reveals a significantly mosaic structure

High-resolution physical maps of the genomes of three Rhodobacter capsulatus strains, derived from ordered cosmid libraries, were aligned. The 1.2-Mb segment of the SB1003 genome studied here is adjacent to a 1-Mb region analyzed previously [Fonstein, M., Nikolskaya, T. & Haselkorn, H. (1995) J. Bacteriol. 177, 2368-2372]. Probes derived from the ordered cosmid set of R. capsulatus SB1003 were used to link cosmids from the St. Louis and 2.3.1 strain libraries. Cosmids selected this way did not merge into a single contig but formed several unlinked groups. EcoRV restriction maps of the ordered cosmids were then constructed using lambda terminase and fused to derive fragments of the chromosomal map. In order to link these fragments, their ends were transcribed to produce secondary probes for hybridization to gridded cosmid libraries of the same strains. This linking reduced the number of subcontigs to three for the St. Louis strain and one for the 2.3.1 strain. Hybridization of the same probes back to the ordered cosmid set of SB1003 positioned the subcontigs on the high-resolution physical map of SB1003. The final alignment of the restriction maps shows numerous large and small translocations in this 1.2-Mb chromosomal region of the three Rhodobacter strains. In addition, the chromosomes of the three strains, whose fine-structure maps can now be compared over 2.2 Mb, are seen to contain regions of 15-80 kb in which restriction sites are highly polymorphic, interspersed among regions in which the positions of restriction sites are highly conserved.

The hyperthermophilic bacterium Thermotoga neapolitana possesses two isozymes of the 3-phosphoglycerate kinase/triosephosphate isomerase fusion protein

The phosphoglycerate kinase (pgk), triosephosphate isomerase (tpi), and enolase (eno) genes from Thermotoga neapolitana have been cloned and expressed in Escherichia coli. In high copy number, the pgk gene complemented an E. coli pgk- strain. In T. neapolitana, the pgk and tpi genes appear to be fused and eno is near those genes. Like T. maritima, T. neapolitaná produces phosphoglycerate kinase as both an individual enzyme and a fusion protein with triosephosphate isomerase, and triosephosphate isomerase activity is not found without associated phosphoglycerate kinase activity. Unlike T. maritima, which forms only a 70-kDa fusion protein, T. neapolitana expresses both 73-kDa and 81-kDa isozymes of this fusion protein. These isozymes are present in both T. neapolitana cells and in E. coli cells expressing T. neapolitana genes.