Examination of Prokaryotic Multipartite Genome Evolution through Experimental Genome Reduction

Many bacteria carry two or more chromosome-like replicons. This occurs in pathogens such as Vibrio cholerea and Brucella abortis as well as in many N₂-fixing plant symbionts including all isolates of the alfalfa root-nodule bacteria Sinorhizobium meliloti. Understanding the evolution and role of this multipartite genome organization will provide significant insight into these important organisms; yet this knowledge remains incomplete, in part, because technical challenges of large-scale genome manipulations have limited experimental analyses. The distinct evolutionary histories and characteristics of the three replicons that constitute the S. meliloti genome (the chromosome (3.65 Mb), pSymA megaplasmid (1.35 Mb), and pSymB chromid (1.68 Mb)) makes this a good model to examine this topic. We transferred essential genes from pSymB into the chromosome, and constructed strains that lack pSymB as well as both pSymA and pSymB. This is the largest reduction (45.4%, 3.04 megabases, 2866 genes) of a prokaryotic genome to date and the first removal of an essential chromid. Strikingly, strains lacking pSymA and pSymB (ΔpSymAB) lost the ability to utilize 55 of 74 carbon sources and various sources of nitrogen, phosphorous and sulfur, yet the ΔpSymAB strain grew well in minimal salts media and in sterile soil. This suggests that the core chromosome is sufficient for growth in a bulk soil environment and that the pSymA and pSymB replicons carry genes with more specialized functions such as growth in the rhizosphere and interaction with the plant. These experimental data support a generalized evolutionary model, in which non-chromosomal replicons primarily carry genes with more specialized functions. These large secondary replicons increase the organism''s niche range, which offsets their metabolic burden on the cell (e.g. pSymA). Subsequent co-evolution with the chromosome then leads to the formation of a chromid through the acquisition of functions core to all niches (e.g. pSymB).     

Dinucleotide compositional analysis of Sinorhizobium meliloti using the genome signature: distinguishing chromosomes and plasmids

The symbiotic N₂-fixing α-proteobacterium Sinorhizobium meliloti has three replicons: a circular chromosome (3.7 Mb) and two smaller replicons, pSymA (1.4 Mb) and pSymB (1.7 Mb). Sequence analysis has revealed that an essential gene is carried on pSymB, which brings into question whether pSymB should be considered a chromosome or a plasmid. Based on the criterion that essential genes define a chromosome, several species have been shown to have multiple chromosomes. Many of these species are part of the α subdivision of the Proteobacteria family. Here, additional justification is presented for designating the pSymB replicon as a chromosome. It is shown that chromosomes within a species share a more similar dinucleotide composition, or genome signature, than plasmids do with the host chromosome(s). Dinucleotide signatures were determined for each of the S. meliloti replicons, and, consistent with the suggestion that pSymB is a chromosome, it is shown that the pSymB signature more closely resembles that of the S. meliloti chromosome, while the pSymA signature is typical of other α-proteobacterial plasmids. Electronic Publication     

Sequence Analysis of the 144-Kilobase Accessory Plasmid pSmeSM11a, Isolated from a Dominant Sinorhizobium meliloti Strain Identified during a Long-Term Field Release Experiment

The genome of Sinorhizobium meliloti type strain Rm1021 consists of three replicons: the chromosome and two megaplasmids, pSymA and pSymB. Additionally, many indigenous S. meliloti strains possess one or more smaller plasmids, which represent the accessory genome of this species. Here we describe the complete nucleotide sequence of an accessory plasmid, designated pSmeSM11a, that was isolated from a dominant indigenous S. meliloti subpopulation in the context of a long-term field release experiment with genetically modified S. meliloti strains. Sequence analysis of plasmid pSmeSM11a revealed that it is 144,170 bp long and has a mean G+C content of 59.5 mol%. Annotation of the sequence resulted in a total of 160 coding sequences. Functional predictions could be made for 43% of the genes, whereas 57% of the genes encode hypothetical or unknown gene products. Two plasmid replication modules, one belonging to the repABC replicon family and the other belonging to the plasmid type A replicator region family, were identified. Plasmid pSmeSM11a contains a mobilization (mob) module composed of the type IV secretion system-related genes traG and traA and a putative mobC gene. A large continuous region that is about 42 kb long is very similar to a corresponding region located on S. meliloti Rm1021 megaplasmid pSymA. Single-base-pair deletions in the homologous regions are responsible for frameshifts that result in nonparalogous coding sequences. Plasmid pSmeSM11a carries additional copies of the nodulation genes nodP and nodQ that are responsible for Nod factor sulfation. Furthermore, a tauD gene encoding a putative taurine dioxygenase was identified on pSmeSM11a. An acdS gene located on pSmeSM11a is the first example of such a gene in S. meliloti. The deduced acdS gene product is able to deaminate 1-aminocyclopropane-1-carboxylate and is proposed to be involved in reducing the phytohormone ethylene, thus influencing nodulation events. The presence of numerous insertion sequences suggests that these elements mediated acquisition of accessory plasmid modules.     

The complete genome sequence of the dominant Sinorhizobium meliloti field isolate SM11 extends the S. meliloti pan-genome

Isolates of the symbiotic nitrogen-fixing species Sinorhizobium meliloti usually contain a chromosome and two large megaplasmids encoding functions that are absolutely required for the specific interaction of the microsymbiont with corresponding host plants leading to an effective symbiosis. The complete genome sequence, including the megaplasmids pSmeSM11c (related to pSymA) and pSmeSM11d (related to pSymB), was established for the dominant, indigenous S. meliloti strain SM11 that had been isolated during a long-term field release experiment with genetically modified S. meliloti strains. The chromosome, the largest replicon of S. meliloti SM11, is 3,908,022 bp in size and codes for 3785 predicted protein coding sequences. The size of megaplasmid pSmeSM11c is 1,633,319 bp and it contains 1760 predicted protein coding sequences whereas megaplasmid pSmeSM11d is 1,632,395 bp in size and comprises 1548 predicted coding sequences. The gene content of the SM11 chromosome is quite similar to that of the reference strain S. meliloti Rm1021. Comparison of pSmeSM11c to pSymA of the reference strain revealed that many gene regions of these replicons are variable, supporting the assessment that pSymA is a major hot-spot for intra-specific differentiation. Plasmids pSymA and pSmeSM11c both encode unique genes. Large gene regions of pSmeSM11c are closely related to corresponding parts of Sinorhizobium medicae WSM419 plasmids. Moreover, pSmeSM11c encodes further novel gene regions, e.g. additional plasmid survival genes (partition, mobilisation and conjugative transfer genes), acdS encoding 1-aminocyclopropane-1-carboxylate deaminase involved in modulation of the phytohormone ethylene level and genes having predicted functions in degradative capabilities, stress response, amino acid metabolism and associated pathways. In contrast to Rm1021 pSymA and pSmeSM11c, megaplasmid pSymB of strain Rm1021 and pSmeSM11d are highly conserved showing extensive synteny with only few rearrangements. Most remarkably, pSmeSM11b contains a new gene cluster predicted to be involved in polysaccharide biosynthesis. Compilation of the S. meliloti SM11 genome sequence contributes to an extension of the S. meliloti pan-genome.     

Spatiotemporal choreography of chromosome and megaplasmids in the Sinorhizobium meliloti cell cycle

A considerable share of bacterial species maintains multipartite genomes. Precise coordination of genome replication and segregation with cell growth and division is vital for proliferation of these bacteria. The α‐proteobacterium Sinorhizobium meliloti possesses a tripartite genome composed of one chromosome and the megaplasmids pSymA and pSymB. Here, we investigated the spatiotemporal pattern of segregation of these S. meliloti replicons at single cell level. Duplication of chromosomal and megaplasmid origins of replication occurred spatially and temporally separated, and only once per cell cycle. Tracking of FROS (fluorescent repressor operator system)‐labelled origins revealed a strict temporal order of segregation events commencing with the chromosome followed by pSymA and then by pSymB. The repA2B2C2 region derived from pSymA was sufficient to confer the spatiotemporal behaviour of this megaplasmid to a small plasmid. Altering activity of the ubiquitous prokaryotic replication initiator DnaA, either positively or negatively, resulted in an increase in replication initiation events or G1 arrest of the chromosome only. This suggests that interference with DnaA activity does not affect replication initiation control of the megaplasmids.     

Presence of 16S rRNA genes in multiple replicons in Azospirillum brasilense

Rhizosphere and endophytic Azospirillum brasilense isolates recovered from sugarcane plants and the reference strains Sp7 and Cd were analyzed for plasmid occurrence. All of the 26 A. brasilense isolates analyzed harbored from five to eight replicons. Several strains contained small plasmids from 45 to 70 kb, but all of the isolates harbored other plasmids ranging from 100 to 290 kb and two megareplicons of approximately 1700 and over 1800 kb. Most of the strains contained a replicon with a size of either 570 or 630 kb, and another large 910- or 980-kb replicon. The 1700-kb megareplicon and some others around 600 kb strongly hybridized to 16S rDNA genes, while the 910- or 980-kb replicons hybridized only slightly. This suggests that the A. brasilense genome is composed of multiple minichromosomes instead of a single circular chromosome. The apparent genome complexity of A. brasilense deserves to eventually be resolved by complete genome sequencing.     

Genome sequence of the saprophyte Leptospira biflexa provides insights into the evolution of Leptospira and the pathogenesis of leptospirosis

Leptospira biflexa is a free-living saprophytic spirochete present in aquatic environments. We determined the genome sequence of L. biflexa, making it the first saprophytic Leptospira to be sequenced. The L. biflexa genome has 3,590 protein-coding genes distributed across three circular replicons: the major 3,604 chromosome, a smaller 278-kb replicon that also carries essential genes, and a third 74-kb replicon. Comparative sequence analysis provides evidence that L. biflexa is an excellent model for the study of Leptospira evolution; we conclude that 2052 genes (61%) represent a progenitor genome that existed before divergence of pathogenic and saprophytic Leptospira species. Comparisons of the L. biflexa genome with two pathogenic Leptospira species reveal several major findings. Nearly one-third of the L. biflexa genes are absent in pathogenic Leptospira. We suggest that once incorporated into the L. biflexa genome, laterally transferred DNA undergoes minimal rearrangement due to physical restrictions imposed by high gene density and limited presence of transposable elements. In contrast, the genomes of pathogenic Leptospira species undergo frequent rearrangements, often involving recombination between insertion sequences. Identification of genes common to the two pathogenic species, L. borgpetersenii and L. interrogans, but absent in L. biflexa, is consistent with a role for these genes in pathogenesis. Differences in environmental sensing capacities of L. biflexa, L. borgpetersenii, and L. interrogans suggest a model which postulates that loss of signal transduction functions in L. borgpetersenii has impaired its survival outside a mammalian host, whereas L. interrogans has retained environmental sensory functions that facilitate disease transmission through water.     

Why Genes Evolve Faster on Secondary Chromosomes in Bacteria

In bacterial genomes composed of more than one chromosome, one replicon is typically larger, harbors more essential genes than the others, and is considered primary. The greater variability of secondary chromosomes among related taxa has led to the theory that they serve as an accessory genome for specific niches or conditions. By this rationale, purifying selection should be weaker on genes on secondary chromosomes because of their reduced necessity or usage. To test this hypothesis we selected bacterial genomes composed of multiple chromosomes from two genera, Burkholderia and Vibrio, and quantified the evolutionary rates (dN and dS) of all orthologs within each genus. Both evolutionary rate parameters were faster among orthologs found on secondary chromosomes than those on the primary chromosome. Further, in every bacterial genome with multiple chromosomes that we studied, genes on secondary chromosomes exhibited significantly weaker codon usage bias than those on primary chromosomes. Faster evolution and reduced codon bias could in turn result from global effects of chromosome position, as genes on secondary chromosomes experience reduced dosage and expression due to their delayed replication, or selection on specific gene attributes. These alternatives were evaluated using orthologs common to genomes with multiple chromosomes and genomes with single chromosomes. Analysis of these ortholog sets suggested that inherently fast-evolving genes tend to be sorted to secondary chromosomes when they arise; however, prolonged evolution on a secondary chromosome further accelerated substitution rates. In summary, secondary chromosomes in bacteria are evolutionary test beds where genes are weakly preserved and evolve more rapidly, likely because they are used less frequently.     

The Glycobiome of the Rumen Bacterium Butyrivibrio proteoclasticus B316T Highlights Adaptation to a Polysaccharide-Rich Environment

Determining the role of rumen microbes and their enzymes in plant polysaccharide breakdown is fundamental to understanding digestion and maximising productivity in ruminant animals. Butyrivibrio proteoclasticus B316^T is a Gram-positive, butyrate-forming rumen bacterium with a key role in plant polysaccharide degradation. The 4.4Mb genome consists of 4 replicons; a chromosome, a chromid and two megaplasmids. The chromid is the smallest reported for all bacteria, and the first identified from the phylum Firmicutes. B316 devotes a large proportion of its genome to the breakdown and reassembly of complex polysaccharides and has a highly developed glycobiome when compared to other sequenced bacteria. The secretion of a range of polysaccharide-degrading enzymes which initiate the breakdown of pectin, starch and xylan, a subtilisin family protease active against plant proteins, and diverse intracellular enzymes to break down oligosaccharides constitute the degradative capability of this organism. A prominent feature of the genome is the presence of multiple gene clusters predicted to be involved in polysaccharide biosynthesis. Metabolic reconstruction reveals the absence of an identifiable gene for enolase, a conserved enzyme of the glycolytic pathway. To our knowledge this is the first report of an organism lacking an enolase. Our analysis of the B316 genome shows how one organism can contribute to the multi-organism complex that rapidly breaks down plant material in the rumen. It can be concluded that B316, and similar organisms with broad polysaccharide-degrading capability, are well suited to being early colonizers and degraders of plant polysaccharides in the rumen environment.     

Combined Analysis of Variation in Core,Accessory and Regulatory Genome Regions Provides a Super-Resolution View into the Evolution of Bacterial Populations

The use of whole-genome phylogenetic analysis has revolutionized our understanding of the evolution and spread of many important bacterial pathogens due to the high resolution view it provides. However, the majority of such analyses do not consider the potential role of accessory genes when inferring evolutionary trajectories. Moreover, the recently discovered importance of the switching of gene regulatory elements suggests that an exhaustive analysis, combining information from core and accessory genes with regulatory elements could provide unparalleled detail of the evolution of a bacterial population. Here we demonstrate this principle by applying it to a worldwide multi-host sample of the important pathogenic E. coli lineage ST131. Our approach reveals the existence of multiple circulating subtypes of the major drug–resistant clade of ST131 and provides the first ever population level evidence of core genome substitutions in gene regulatory regions associated with the acquisition and maintenance of different accessory genome elements.     

Complete genome sequencing of Agrobacterium sp. H13-3, the former Rhizobium lupini H13-3, reveals a tripartite genome consisting of a circular and a linear chromosome and an accessory plasmid but lacking a tumor-inducing Ti-plasmid

Agrobacterium sp. H13-3, formerly known as Rhizobium lupini H13-3, is a soil bacterium that was isolated from the rhizosphere of Lupinus luteus. The isolate has been established as a model system for studying novel features of flagellum structure, motility and chemotaxis within the family Rhizobiaceae. The complete genome sequence of Agrobacterium sp. H13-3 has been established and the genome structure and phylogenetic assignment of the organism was analysed. For de novo sequencing of the Agrobacterium sp. H13-3 genome, a combined strategy comprising 454-pyrosequencing on the Genome Sequencer FLX platform and PCR-based amplicon sequencing for gap closure was applied. The finished genome consists of three replicons and comprises 5,573,770 bases. Based on phylogenetic analyses, the isolate could be assigned to the genus Agrobacterium biovar I and represents a genomic species G1 strain within this biovariety. The highly conserved circular chromosome (2.82 Mb) of Agrobacterium sp. H13-3 mainly encodes housekeeping functions characteristic for an aerobic, heterotrophic bacterium. Agrobacterium sp. H13-3 is a motile bacterium driven by the rotation of several complex flagella. Its behaviour towards external stimuli is regulated by a large chemotaxis regulon and a total of 17 chemoreceptors. Comparable to the genome of Agrobacterium tumefaciens C58, Agrobacterium sp. H13-3 possesses a linear chromosome (2.15 Mb) that is related to its reference replicon and features chromosomal and plasmid-like properties. The accessory plasmid pAspH13-3a (0.6 Mb) is only distantly related to the plasmid pAtC58 of A. tumefaciens C58 and shows a mosaic structure. A tumor-inducing Ti-plasmid is missing in the sequenced strain H13-3 indicating that it is a non-virulent isolate.     

A Personal Reflection on the Replicon Theory: From R1 Plasmid to Replication Timing Regulation in Human Cells

Fifty years after the Replicon Theory was originally presented, detailed mechanistic insight into prokaryotic replicons has been obtained and rapid progress is being made to elucidate the more complex regulatory mechanisms of replicon regulation in eukaryotic cells. Here, I present my personal perspectives on how studies of model replicons have contributed to our understanding of the basic mechanisms of DNA replication as well as the evolution of replication regulation in human cells. I will also discuss how replication regulation contributes to the stable maintenance of the genome and how disruption of replication regulation leads to human diseases.     

Draft genome sequence of Sinorhizobium meliloti RU11/001, a model organism for flagellum structure,motility and chemotaxis

Sinorhizobium meliloti of the order Rhizobiales is a symbiotic nitrogen-fixing bacterium nodulating plants of the genera Medicago, Trigonella and Melilotus and hence is of great agricultural importance. In its free-living state it is motile and capable of modulating its movement patterns in response to chemical attractants. Here, the draft genome consisting of a circular chromosome, the megaplasmids pSymA and pSymB and three accessory plasmids of Sinorhizobium meliloti RU11/001, a model organism for flagellum structure, motility and chemotaxis, is reported.     

Novel DNA Sequences from Natural Strains of the Nitrogen-Fixing Symbiotic Bacterium Sinorhizobium meliloti

Variation in genome size and content is common among bacterial strains. Identifying these naturally occurring differences can accelerate our understanding of bacterial attributes, such as ecological specialization and genome evolution. In this study, we used representational difference analysis to identify potentially novel sequences not present in the sequenced laboratory strain Rm1021 of the nitrogen-fixing bacterium Sinorhizobium meliloti. Using strain Rm1021 as the driver and the type strain of S. meliloti ATCC 9930, which has a genome size ~370 kilobases bigger than that of strain Rm1021, as the tester, we identified several groups of sequences in the ATCC 9930 genome not present in strain Rm1021. Among the 85 novel DNA fragments examined, 55 showed no obvious homologs anywhere in the public databases. Of the remaining 30 sequences, 24 contained homologs to the Rm1021 genome as well as unique segments not found in Rm1021, 3 contained sequences homologous to those published for another S. meliloti strain but absent in Rm1021, 2 contained sequences homologous to other symbiotic nitrogen-fixing bacteria (Rhizobium etli and Bradyrhizobium japonicum), and 1 contained a sequence homologous to a gene in a non-nitrogen-fixing species, Pseudomonas sp. NK87. Using PCR, we assayed the distribution of 12 of the above 85 novel sequences in a collection of 59 natural S. meliloti strains. The distribution varied widely among the 12 novel DNA fragments, from 1.7% to 72.9%. No apparent correlation was found between the distribution of these novel DNA sequences and their genotypes obtained using multilocus enzyme electrophoresis. Our results suggest potentially high rates of gene gain and loss in S. meliloti genomes.     

Comparative genome architecture and dynamics in bacteria

This article reviews current knowledge about the organization of bacterial genomes, of which a number of components (replicons), namely chromosomes, plasmids and prophages, have been well characterized. The historical position of the acceptance of the idea of circularity and unit copy number of replicons in bacterial cells has been readdressed by new methods of genome analysis, particularly pulsed-field gcl electrophoresis, which have facilitated identification of variation in replicon number and distinction whether the replicons are circular or linear DNA structures. Much has also been learnt about the origins of DNA replication in replicons and how they function via the controlling role of specific proteins or RNA. A related aspect is thc problem of how the replication products are stabilized, segregated and partitioned into daughter cells at cell division. Our understanding of replicons has also been improved by application of state-of-the-art computer software methods of comparative DNA and protein sequence analysis. This knowledge has provided insights into the fundamental nature of these processes and their origin and evolution in single-cell and multicellular organisms.