Review | Open | Published:
Evaluating genome dynamics: the constraints on rearrangements within bacterial genomes
Genome Biologyvolume 1, Article number: reviews0006.1 (2000)
Inversions and translocations distinguish the genomes of closely related bacterial species, but most of these rearrangements preserve the relationship between the rearranged fragments and the axis of chromosome replication. Within species, such rearrangements are found less frequently, except in the case of clinical isolates of human pathogens, where rearrangements are very frequent
The evolution of biological diversity through genetic change raises questions about how much variation one can expect among closely related genomes. Some answers are emerging from the application of two different technologies to comparative genomics of bacteria. One, complete genome sequencing, is providing detailed blueprints of one or a few examples of each genome of interest. The second, physical mapping by pulsed-field gel electrophoresis (PFGE), is providing 'skeleton' views of large numbers of closely related genomes. Together, these technologies are providing insights into the dynamics of genome plasticity that are both detailed and broad. At this early stage in comparative genomics, the main generalizations that are emerging concerning rearrangements in bacterial genome organization are as follows. First, large chromosome inversions and translocations are common, even between closely related species. Second, chromosome inversions are usually symmetric around the axis of DNA replication. Third, chromosomal rearrangements are less common within species, but a dramatic increase in the frequency of inversions and translocation seems to be associated with the ability of bacteria to infect eukaryotic hosts, possibly reflecting a bacterial response to the challenges posed by the immune system.
The underlying causes of rearrangement
The bacterial RecA protein is required for damage to chromosomes - in particular chromosome breaks - to be repaired, and it acts by using a duplicate copy of the damaged sequence as a template for repair. The template is normally the homologous sequence on a sister chromosome, but when sequences are present in multiple copies within a genome, RecA can promote recombination between paralogs. Such recombination events can result in rearrangements in the order of genes on the chromosome . Thus, recombination between repeated sequences that are in the same orientation as each other (direct repeats) can result in tandem duplication of the region bounded by the repeat sequences (Figure 1a). These duplications are usually unstable, unless maintained by selection or by the accumulation of sufficient mutations to avoid subsequent recombination. Some large duplications have been stabilized during evolution . Recombination between direct repeats can also result in deletion of the intervening sequence, generating a fragment that can potentially insert back into the genome at the site of another copy of the repeat, thus generating a translocation (Figure 1b). Recombination between repeat sequences that are in the opposite orientation to each other (inverse repeats) can result in inversion of the intervening sequence (Figure 1c).
To these phenomena one can add the acquisition of new DNA by horizontal transfer from another genome, which in addition to introducing new genetic information may upset the stability of a genome and trigger other compensating rearrangements. The sequences involved in RecA-mediated rearrangements are usually long repeats, such as rRNA operons, transposons and IS (insertion sequence) elements. Sequences as short as 10-100 nucleotides can also be substrates for homologous recombination, but this is usually limited to sequences in close proximity to one another .
Another outcome of homologous recombination between repeated sequences is gene conversion - homogenization of the sequences within a gene family, to prevent the divergence of repeated sequences - which maintains the sequence similarity required for RecA-mediated recombination over evolutionary time scales . RecA activity is thus a double-edged sword: it is needed to maintain the chromosome integrity required for completing replication, but it also promotes rearrangements within the genome. An interesting exception is that Buchnera lacks a recA gene but compensates by having over 100 copies of its entire genome per cell .
Historical background and new techniques
Microbial genome analysis has its origins in the intensive laboratory analysis of just a few bacterial strains, in particular Escherichia coli K-12 and Salmonella typhimurium LT2 (proper name: Salmonella enterica serovar Typhimurium). These species had a common ancestor approximately 140 million years ago . Their genetic maps are almost identical in organization and their major phenotypic differences can be explained by the horizontal acquisition of DNA segments into one or other of the species. With the exception of a large inversion around the terminus of replication, the two genomes seem to be stable in organization and to be diverging by the accumulation of point mutations. A high frequency of recombination in the terminus region is related to the mechanism of chromosome separation after replication, and different inversions around the terminus are found in other closely related bacteria . The general conclusion drawn from these early comparisons was that bacterial genomes were stable in organization. Within the past few years, however, the application of new technologies to the analysis of the genomes of a wide variety of bacterial species has challenged this view.
The two most important experimental techniques for comparative analysis of genome organization have been whole-genome sequencing and physical mapping of genome organization. Whole-genome sequencing provides complete information on a genome, facilitating many types of analysis including comparative analysis of genome organization. The TIGR Microbial Database  currently lists dozens of completed bacterial genome sequences and over one hundred in progress. In most cases a few genomes from each species are being sequenced. The availability of a genome sequence, while not essential, also facilitates the physical analysis of that genome and of related genomes. Physical mapping, most often by PFGE in conjunction with restriction digestion and probing for specific sequences, is used to generate skeleton structures of genome organization and is suited to screening large numbers of strains. The most informative comparisons, in terms of evaluating genome dynamics, are those made between close relatives - either sister species or isolates within a species.
Genome rearrangements at three levels of comparison
Comparisons between genome arrangements in related bacteria have been made at several levels: interspecific, intraspecific (serovars, or immunologically detectable variants, and biovars, or biochemically detectable variants within a single species), and within presumed clonal populations. The phylogenetic relationships between many of the bacterial species referred to in the following discussion, based on an analysis of their 16S rDNA, are shown in Figure 2.
Complete genome sequence data for three pairs of related species - Vibrio cholera - Escherichia coli; Streptococcus pneumoniae - Streptococcus pyrogenes; and Mycobacterium tuberculosis - Mycobacterium leprae - has been used to compare the positions of conserved sequences within each genome. The comparisons reveal in each case a distinct X-shaped pattern in scatterplots, suggestive of large chromosomal inversions that reverse the genomic sequence symmetrically around the axis of replication . Similarly, Chlamydia trachomatis and Chlamydia pneumoniae differ by multiple large inversions, apparently oriented around the axis of replication. In addition, the region around the terminus of replication in these two species is subject to a high rate of reorganization . Table 1 summarizes all the interspecific rearrangements referred to in this section.
As discussed earlier, another important type of genome rearrangement is translocation. Comparative genomics using whole-genome sequences shows that the genomes of the close relatives Mycoplasma genitalium and Mycoplasma pneumoniae  can be subdivided into six segments, which are ordered differently in the two species. Within each segment the order of genes is conserved, and the increased size of the M. pneumoniae chromosome is due mainly to gene duplications. In both species there is strong uniformity of direction of transcription, and this direction is not changed by the translocations .
Other interspecific genome comparisons reveal that inversions, translocations and deletions typically distinguish closely related species. Thus, Neisseria meningitidis and Neisseria gonorrhoeae differ by multiple translocations and/or inversions of blocks of genetic markers within a 500 kilobase region . Mycobacterium tuberculosis differs from the attenuated vaccine strain Mycobacterium bovis BCG Pasteur in carrying two tandem duplications in its chromosome, of 29 kb and 36 kb respectively . In addition, clinical isolates of Mycobacterium tuberculosis differ from each other in having deletions of up to several kilobases, probably linked to homologous recombination between multiple copies of IS6110 in the genome . Finally, the chromosomal locations of 30 putative orthologs between Bacillus subtilis and Bacillus cereus are arranged in an apparently random manner , similar to what is seen when comparing the genomes of very distantly related organisms. Also, within B. subtilis several variants created by X-ray mutagenesis have large inversions (1,700-1,900 kb) around either the axis of replication, a 100 kb translocation, and smaller duplications and deletions . In conclusion, in every case where interspecific comparisons have been made, clear evidence of large chromosomal rearrangements has been found.
Comparisons within species, including comparisons between different serovars and biovars of the same species, reveal that rearrangements are less common within than between species (summarized in Table 2). Using PFGE after restriction digestion targeted to cut the genome within the conserved and repetitive rRNA operons, 'rrn genomic skeletons' were established for isolates of many serovars of Salmonella enterica . The order of fragments, which is ABCDEFG in S. typhimurium and E. coli K-12, is conserved in most Salmonella serovars, most of which are host-generalists. In S. typhi and S. paratyphi C (which have human hosts), however, and in S. pullorum and S. gallinarum (which have fowl hosts), these fragments are rearranged. Thus, of 127 natural isolates of S. typhi examined, 21 different genome orders were found, all postulated to be due to inversions and translocations with end-points in rrn operons . A feature of these rearrangements is that the distance from the origin of chromosome replication is well conserved, as is the direction of transcription relative to the direction of chromosome replication. A similar PFGE analysis has been made of Shigella, the human pathogenic form of E. coli . This showed that of the four traditional Shigella subgroups (often referred to as species), S. boydii and S. sonnei had chromosomal arrangements identical to E. coli K-12, while S. dysenteriae and S. flexneri had different large rearrangements . Interestingly, the Shiga toxin genes on the S. dysenteriae chromosome are bracketed by IS600 sequences, and increased toxin production is caused by tandem amplification via recombination between the IS600 elements .
The genome sequence of the strain Lactococcus lactis IL1403 has been determined  and physical genome maps have been created for several Lactococcus lactis strains (subspecies lactis and subspecies cremoris strains) and strains of the related Streptococcus thermophilus [20,21]. Within each group, strains were similar with the exception of the L. Lactis subspecies cremoris, where different strains were polymorphic, in part due to an inversion of half the chromosome. This inversion is due to a homologous recombination event between two defective copies of IS905 and does not alter the symmetry of the replication origin and terminus, oriC and terC . Comparison of the physical and genetic maps of strains representing two serovars of Leptospira interrogans suggests that at least two inversions in the large replicon distinguish their genomes .
Brucella is a Gram-negative bacterium pathogenic for animals. The genus is divided into six species and numerous biovars. Physical maps of the genomes of reference strains in each species show a high conservation of restriction sites and the presence of two chromosomes. The exception is a large inversion in the small chromosome of B. abortus. But physical mapping of the genomes of the four biovars within one of these species, Brucella suis, reveals differences in both chromosome number and size. These differences can be explained by rearrangements due to homologous recombination between the three rrn loci in the genome . It is proposed that the ancestor of Brucella had a single chromosome and that recombination, probably between rrn genes, led to the creation of two chromosomes.
Multiple chromosomes are also found in other bacteria, mostly within the proteobacteriaceae, including Rhodobacter sphaeroides, Leptospira interrogans, Rhizobium spp., Burkholderia cepacia, Agrobacterium spp., and Ochrobactrum anthropi . There is no evidence that the presence of multiple chromosomes in these genomes is related to a common phylogeny, since it is not always shared by all members of a genus (for example, the Rhodobacter genus) or even by all strains of the same species (for example, Brucella suis). In Streptomyces, the most common rearrangements found are sequence and length variations in the terminal inverted repeats (TIRs) at the ends of the linear replicons. This variation is due to homologous recombination between repetitive sequences and results in amplifications, deletions, a high frequency of spontaneous mutations, and the transfer of sequences between different chromosome arms . The exchange of telomeric regions has also been described for linear replicons in the unrelated bacterium Borrelia burgdorferi . Instability may not be a particular feature of linearity, however, because when Streptomyces chromosomes are circularized they remain unstable in these regions . The lack of housekeeping genes in a large region at each end of the Streptomyces linear chromosomes probably permits the detection of deletions at a high frequency because they do not affect cell viability under laboratory conditions.
Clinical and clonal populations
The closest relatives that have been subjected to genome organization analysis are presumed clonal derivatives associated with clinical infections (summarized in Table 3). Bordetella pertussis strains from a whooping cough outbreak in Canada were subjected to restriction-enzyme genome mapping. Among 70 isolates, presumed to be descended from the same starting clone, 14 different types were found (distinguished by restriction fragment length polymorphism, RFLP). Representatives of these 14 types were further analyzed and shown to have 11 different genome orders, due in each case to large chromosomal inversions . A similar analysis among different laboratory strains also revealed frequent large inversions . B. pertussis carries about 100 copies of the 1 kilobase insertion sequence IS481, providing many targets for homologous recombination. The positions of the origin and terminus of replication are unknown in Bordetella, but because almost all of the inversions are around the same axis, it seems likely that they are in fact symmetric about the origin-terminus axis, as has been observed in other species. The genomes of clinical isolates of Pseudomonas aeruginosa from cystic fibrosis patients were analyzed by PFGE, revealing that 50% of them have large chromosomal inversions . Most of these inversions are approximately symmetric about the replication axis. It is not known if these rearrangements confer any advantage on the strains in colonizing the lung habitat.
PFGE analysis of the genomes of 30 Neisseria meningitidis epidemic strains belonging to the ET-5 complex, isolated from various parts of the world over a period of 20 years, revealed 10 different types, including some with genome order rearrangements . A striking feature of N. meningiditis revealed by complete genome sequencing [32,33] is the presence of hundreds of repetitive elements that could contribute to genome rearrangements important for evasion of the host immune system. Finally, within a defined lineage of N. gonorrhoeae strains with pilin variations, an inversion of more than one third of the chromosome was found . The end points of the inversion are within a multicopy gene family involved in pilin production.
Constraints on the frequency of rearrangement
Rearrangement could be constrained by the number and size of repetitive sequences in a genome. Most bacterial genomes contain multiple copies of some highly expressed genes (such as rrn genes) or have multiple copies of insertion sequences, however. While the relative positions of these sequences could influence which rearrangements occur most frequently, recombination between short repeats and the mobility of IS elements should increase the variety of rearrangements. A second potential constraint is the rate of recombination, but experimental data from S. typhimurium show that rates of recombination between long repeat sequences are at least as high as nucleotide substitution rates . For example, the rate of inversion between the tuf genes, approximately 10-8 per cell per generation, is equal to the rate of nucleotide substitution within the same genes [3,35]. If one can generalize from this, then recombination can rearrange genome organization as fast as genomes diverge by nucleotide substitution. A third possible constraint is that the fitness of the rearranged genomes is in general reduced. Indeed, inversions that reverse the orientation of sequences on either side of the replication terminus of S. typhimurium and E. coli usually occur very infrequently, or make the bacteria very unfit . Inversions that do not alter the replication axis (that is, inversions that do not change the distance of genes from the origin of replication, or their orientation relative to the direction of replication) may be the least disruptive in terms of fitness, but this has never been rigorously tested. In conclusion, fitness costs may be an important constraint on the fixation of genome rearrangements in bacteria, but there are very few relevant measurements.
Comparing these theoretical constraints on genome rearrangements with the data from natural bacterial isolates, several patterns emerge. One is that inversions and translocations are very common between even closely related species (compare Figure 2 and Table 1). This presumably reflects the frequency of the rearrangement events, on a time scale of tens to hundreds of millions of years. In general, however, the variety of large rearrangements found is quite limited. Almost all rearrangements conserve the axis of chromosome replication, suggesting that this is important for fitness on the evolutionary time scale. Below the species level, within and between serovars/biovars, there are interesting differences in genome stability (Table 2). Some genomes are stable while others, like that of S. typhi, have many different arrangements. Even more striking is the data from clinical isolates that are probably very closely related, separated from one another on a time scale of months to a few years (Table 3). In such cases, a very high frequency of rearranged genomes is found. One thing that S. typhi and the clinical isolates of Neisseria, Pseudomonas and Bordetella have in common is that they have all encountered and survived the human immune system. It is tempting to speculate that such encounters select for variants of the infecting strain. Thus, if an invading pathogen population is targeted by the immune system, bacteria within that population with genome rearrangements may be sufficiently different in phenotype to escape and establish an infection. Variation generated by genome rearrangements has several advantages for a bacterial population invading a complex environment: it occurs at a high frequency, it is reversible, and it can simultaneously alter the expression pattern of many genes.
In conclusion, the main constraint on genome rearrangements on an evolutionary time scale may be bacterial fitness, in particular associated with the global regulation of gene expression patterns and the orderly and efficient replication of the genome. In particular complex environments, however, such as those encountered on invading an eukaryotic host, bacterial fitness may be positively associated, at least on a short time scale, with the generation of genome rearrangements.
Hughes D: Impact of homologous recombination on genome organization and stability. In Organization of the Prokaryotic Genome. Edited by Charlebois RL. Washington: ASM Press;. 1999, 109-128.
Wolfe KH, Shields DC: Molecular evidence for an ancient duplication of the entire yeast genome. Nature. 1997, 387: 708-713. 10.1038/42711.
Abdulkarim F, Hughes D: Homologous recombination between the tuf genes of Salmonella typhimurium. J Mol Biol. 1996, 260: 506-522. 10.1006/jmbi.1996.0418.
Shigenobu S, Watanabe H, Hattori M, Sakaki Y, Ishikawa H: Genome sequence of the endocellular bacterial symbiont of aphids Buchnera sp. APS. Nature. 2000, 407: 81-86. 10.1016/S0168-9002(97)01373-9.
Ochman H, Wilson AC: Evolution in bacteria: evidence for a universal substitution rate in cellular genomes. J Mol Evol. 1987, 26: 74-86.
TIGR Microbial Database. [http://www.tigr.org/tdb/mdb/mdbcomplete.html]
Eisen JA, Heidelberg JF, White O, Salzberg SL: Evidence for symmetric chromosomal inversions around the replication origin in bacteria. Genome Biology. 2000, 1: research0011.1-0011.9. 10.1186/gb-2000-1-6-research0011.
Read TD, Brunham RC, Shen C, Gill SR, Heidelberg JF, White O, Hickey EK, Peterson J, Utterback T, Berry K, et al: Genome sequences of Chlamydia trachomatis MoPn and Chlamydia pneumoniae AR39. Nucleic Acids Res. 2000, 28: 1397-1406. 10.1093/nar/28.6.1397.
Himmelreich R, Plagens H, Hilbert H, Reiner B, Herman R: Comparative analysis of the genomes of the bacteria Mycoplasma pneumoniae and Mycoplasma genitalium. Nucleic Acids Res. 1997, 25: 701-712. 10.1093/nar/25.4.701.
Dempsey JAF, Wallace AB, Cannon JG: The physical map of the chromosome of a serogroup A strain of Neisseria meningitidis shows complex rearrangements relative to the chromosomes of the two mapped strains of the closely related species N. gonorrhoeae. J Bacteriol. 1995, 177: 6390-6400.
Brosch R, Gordon SV, Buchrieser C, Pym AS, Garnier T, Cole ST: Comparative genomics uncovers large tandem chromosomal duplications in Mycobacterium bovis BCG Pasteur. Yeast. 2000, 17: 111-123. 10.1002/1097-0061(20000630)17:2<111::AID-YEA17>3.0.CO;2-G.
Fang Z, Doig C, Kenna DT, Smittipat N, Palittapongarnim P, Watt B, Forbes KJ: IS6110-mediated deletions of wild-type chromosomes of Mycobacterium tuberculosis. J Bacteriol. 1999, 181: 1014-1020.
Økstad OA, Hegna I, Lindback T, Rishovd AL, Kolstø A-B: Genome organization is not conserved between Bacillus cereus and Bacillus subtilis. Microbiology. 1999, 145: 621-631.
Itaya M: Physical map of the Bacillis subtilis 166 genome: evidence for the inversion of an approximately 1900 kb continuous DNA segment, the translocation of an approximately 100 kb segment and the duplication of a 5 kb segment. Microbiology. 1997, 143: 3723-3732.
Sanderson KE, Liu S-L: Chromosomal rearrangements in enteric bacteria. Electrophoresis. 1998, 19: 569-572.
Pupo GM, Lan R, Reeves PR: Multiple independent origins of Shigella clones of Escherichia coli and convergent evolution of many of their characteristics. Proc Natl Acad Sci USA. 2000, 97: 10567-10572. 10.1073/pnas.180094797.
Shu S, Setianingrum E, Zhao L, Li Z-Y, Xu H-X, Kawamura Y, Ezaki T: I-CeuI fragment analysis of the Shigella species: evidence for large-scale chromosome rearrangement in S. dysenteriae and S. flexneri. FEMS Microbiol Lett. 2000, 182: 93-98. 10.1016/S0378-1097(99)00575-3.
McDonough MA, Butterton JR: Spontaneous tandem amplification and deletion of the Shiga toxin operon in Shigella dysenteriae 1. Mol Microbiol. 1999, 34: 1058-1069. 10.1046/j.1365-2958.1999.01669.x.
Bolotin A, Mauger S, Malarme K, Ehrlich SD, Sorokin A: Low-redundancy sequencing of the entire Lactococcus lactis IL1403 genome. Antonie Van Leeuwenhoek. 1999, 76: 27-76. 10.1023/A:1002048720611.
Davidson BE, Kordias N, Dobos M, Hillier AJ: Genomic organization of lactic acid bacteria. Antonie Van Leeuwenhoek. 1996, 70: 161-183.
Le Bourgeois P, Daveran-Mingot M-L, Ritzenthaler P: Genome plasticity among related Lactococcus strains: identification of genetic events associated with macrorestriction polymorphisms. J Bacteriol. 2000, 182: 2481-2491. 10.1128/JB.182.9.2481-2491.2000.
Daveran-Mingot M-L, Campo N, Ritzenthaler P, Le Bourgeois P: A natural large chromosomal inversion in Lactococcus lactis is mediated by homologous recombination between two insertion sequences. J Bacteriol. 1998, 180: 4834-4842.
Zuerner RL, Herrmann JL, Saint Girons I: Comparison of genetic maps for two Leptospira interrogans serovars provides evidence for two chromosomes and interspecies heterogeneity. J Bacteriol. 1993, 175: 5445-5451.
Jumas-Bilak E, Michaux-Charachon S, Bourg G, O'Callaghan D, Ramuz M: Differences in chromosome number in the genus Brucella. Mol Microbiol. 1998, 27: 99-106. 10.1046/j.1365-2958.1998.00661.x.
Fischer G, Wenner T, Decaris B, Leblond P: Chromosomal arm replacement generates a high level of intraspecific polymorphism in the terminal inverted repeats of the linear chromosomal DNA of Streptomyces ambofaciens. Proc Natl Acad Sci USA. 1998, 95: 14296-14301. 10.1073/pnas.95.24.14296.
Casjens S, Palmer N, van Vugt R, Huang WM, Stevensen B, Rosa P, Lathigra R, Sutton G, Peterson J, Dodson RJ, et al: A bacterial genome in flux: the twelve linear and nine circular extra-chromosomal DNAs in an infectious isolate of the Lyme disease spirochete Borrelia burgdorferi. Mol Microbiol. 2000, 35: 490-516. 10.1046/j.1365-2958.2000.01698.x.
Fischer G, Decaris B, Leblond P: Occurrence of deletions, associated with genetic instability in Streptomyces ambofaciens, is independent of the linearity of the chromosomal DNA. J Bacteriol. 1997, 179: 4553-4558.
Stibitz S, Yang M-S: Genomic plasticity in natural populations of Bordetella pertussis. J Bacteriol. 1999, 181: 5512-5515.
Stibitz S, Yang M-S: Genomic fluidity of Bordetella pertussis assessed by a new method for chromosomal mapping. J Bacteriol. 1997, 179: 5820-5826.
Römling U, Schmidt KD, Tümmler B: Large genome rearrangements discovered by the detailed analysis of 21 Pseudomonas aeruginosa clone C isolates found in environment and disease habitats. J Mol Biol. 1997, 271: 386-404. 10.1006/jmbi.1997.1186.
Frøholm LO, Kolstø A-B, Berner J-M, Caugant DA: Genomic rearrangements in Neisseria meningitidis strains of the ET-5 complex. Curr Microbiol. 2000, 40: 372-379. 10.1007/s002840010073.
Parkhill J, Achtman M, James KD, Bentley SD, Churcher C, Klee SR, Morelli G, Basham D, Brown D, Chillingworth T, et al: Complete DNA sequence of a serogroup A strain of Neisseria meningitidis Z2491. Nature. 2000, 404: 502-506. 10.1038/35006655.
Tettelin H, Saunders NJ, Heidelberg J, Jeffries AC, Nelson KE, Eisen JA, Ketchum KA, Hood DW, Peden JF, Dodson RJ, et al: Complete genome sequence of Neisseria meningitidis serogroup B strain MC58. Science. 2000, 287: 1809-1815. 10.1126/science.287.5459.1809.
Gibbs CP, Meyer TF: Genome plasticity in Neisseria gonorrhoeae. FEMS Microbiol Lett. 1996, 145: 173-179. 10.1016/S0378-1097(96)00404-1.
Hughes D: Co-evolution of the tuf genes links gene conversion with the generation of chromosomal inversions. J Mol Biol. 2000, 297: 355-364. 10.1006/jmbi.2000.3587.
My research is supported by the Swedish Natural Sciences Research Council. I thank Ivica Tamas for his assistance in building the phylogenetic tree.