Prediction, identification, and artificial selection of DNA rearrangements in Rhizobium: toward a natural genomic design

Gene amplification mutations originate prior to selective stress in Acinetobacter baylyi

The controversial theory of adaptive amplification states gene amplification mutations are induced by selective environments where they are enriched due to the stress caused by growth restriction on unadapted cells. We tested this theory with three independent assays using an Acinetobacter baylyi model system that exclusively selects for cat gene amplification mutants. Our results demonstrate all cat gene amplification mutant colonies arise through a multistep process. While the late steps occur during selection exposure, these mutants derive from low-level amplification mutant cells that form before growth-inhibiting selection is imposed. During selection, these partial mutants undergo multiple secondary steps generating higher amplification over several days to multiple weeks to eventually form visible high-copy amplification colonies. Based on these findings, amplification in this Acinetobacter system can be explained by a natural selection process that does not require a stress response. These findings have fundamental implications to understanding the role of growth-limiting selective environments on cancer development. We suggest duplication mutations encompassing growth factor genes may serve as new genomic biomarkers to facilitate early cancer detection and treatment, before high-copy amplification is attained.

Unstable Mechanisms of Resistance to Inhibitors of Escherichia coli Lipoprotein Signal Peptidase

Despite increasing evidence suggesting that antibiotic heteroresistance can lead to treatment failure, the significance of this phenomena in the clinic is not well understood, because many clinical antibiotic susceptibility testing approaches lack the resolution needed to reliably classify heteroresistant strains. Here we present G0790, a new globomycin analog and potent inhibitor of the Escherichia coli type II signal peptidase LspA. We demonstrate that in addition to previously known mechanisms of resistance to LspA inhibitors, unstable genomic amplifications containing lspA can lead to modest yet biologically significant increases in LspA protein levels that confer a heteroresistance phenotype.

Clinical development of antibiotics with novel mechanisms of action to kill pathogenic bacteria is challenging, in part, due to the inevitable emergence of resistance. A phenomenon of potential clinical importance that is broadly overlooked in preclinical development is heteroresistance, an often-unstable phenotype in which subpopulations of bacterial cells show decreased antibiotic susceptibility relative to the dominant population. Here, we describe a new globomycin analog, G0790, with potent activity against the Escherichia coli type II signal peptidase LspA and uncover two novel resistance mechanisms to G0790 in the clinical uropathogenic E. coli strain CFT073. Building on the previous finding that complete deletion of Lpp, the major Gram-negative outer membrane lipoprotein, leads to globomycin resistance, we also find that an unexpectedly modest decrease in Lpp levels mediated by insertion-based disruption of regulatory elements is sufficient to confer G0790 resistance and increase sensitivity to serum killing. In addition, we describe a heteroresistance phenotype mediated by genomic amplifications of lspA that result in increased LspA levels sufficient to overcome inhibition by G0790 in culture. These genomic amplifications are highly unstable and are lost after as few as two subcultures in the absence of G0790, which places amplification-containing resistant strains at high risk of being misclassified as susceptible by routine antimicrobial susceptibility testing. In summary, our study uncovers two vastly different mechanisms of resistance to LspA inhibitors in E. coli and emphasizes the importance of considering the potential impact of unstable and heterogenous phenotypes when developing antibiotics for clinical use.

Intracellular Competitions Reveal Determinants of Plasmid Evolutionary Success

Plasmids are autonomously replicating genetic elements that are ubiquitous in all taxa and habitats where they constitute an integral part of microbial genomes. The stable inheritance of plasmids depends on their segregation during cell division and their long-term persistence in a host population is thought to largely depend on their impact on the host fitness. Nonetheless, many plasmids found in nature are lacking a clear trait that is advantageous to their host; the determinants of plasmid evolutionary success in the absence of plasmid benefit to the host remain understudied. Here we show that stable plasmid inheritance is an important determinant of plasmid evolutionary success. Borrowing terminology from evolutionary biology of cellular living forms, we hypothesize that Darwinian fitness is key for the plasmid evolutionary success. Performing intracellular plasmid competitions between non-mobile plasmids enables us to compare the evolutionary success of plasmid genotypes within the host, i.e., the plasmid fitness. Intracellular head-to-head competitions between stable and unstable variants of the same model plasmid revealed that the stable plasmid variant has a higher fitness in comparison to the unstable plasmid. Preemptive plasmid competitions reveal that plasmid fitness may depend on the order of plasmid arrival in the host. Competitions between plasmids characterized by similar stability of inheritance reveal plasmid fitness differences depending on the plasmid-encoded trait. Our results further reveal that competing plasmids can be maintained in coexistence following plasmid fusions that maintain unstable plasmid variants over time. Plasmids are not only useful accessory genetic elements to their host but they are also evolving and replicating entities, similarly to cellular living forms. There is a clear link between plasmid genetics and plasmid evolutionary success – hence plasmids are evolving entities whose fitness is quantifiable.

Prediction and identification of recurrent genomic rearrangements that generate chimeric chromosomes in Saccharomyces cerevisiae

Genomes are dynamic structures. Different mechanisms participate in the generation of genomic rearrangements. One of them is nonallelic homologous recombination (NAHR). This rearrangement is generated by recombination between pairs of repeated sequences with high identity. We analyzed rearrangements mediated by repeated sequences located in different chromosomes. Such rearrangements generate chimeric chromosomes. Potential rearrangements were predicted by localizing interchromosomal identical repeated sequences along the nuclear genome of the Saccharomyces cerevisiae S288C strain. Rearrangements were identified by a PCR-based experimental strategy. PCR primers are located in the unique regions bordering each repeated region of interest. When the PCR is performed using forward primers from one chromosome and reverse primers from another chromosome, the break point of the chimeric chromosome structure is revealed. In all cases analyzed, the corresponding chimeric structures were found. Furthermore, the nucleotide sequence of chimeric structures was obtained, and the origin of the unique regions bordering the repeated sequence was located in the expected chromosomes, using the perfect-match genomic landscape strategy (PMGL). Several chimeric structures were searched in colonies derived from single cells. All of the structures were found in DNA isolated from each of the colonies. Our findings indicate that interchromosomal rearrangements that generate chimeric chromosomes are recurrent and occur, at a relatively high frequency, in cell populations of S. cerevisiae.

Mosaic plasmids are abundant and unevenly distributed across prokaryotic taxa

Mosaic plasmids, plasmids composed of genetic elements from distinct sources, are associated with the spread of antibiotic resistance genes. Transposons are considered the primary mechanism for mosaic plasmid formation, though other mechanisms have been observed in specific instances. The frequency with which mosaic plasmids have been described suggests they may play an important role in plasmid population dynamics. Our survey of the confirmed plasmid sequences available from complete and draft genomes in the RefSeq database shows that 46% of them fit a strict definition of mosaic. Mosaic plasmids are also not evenly distributed over the taxa represented in the database. Plasmids from some genera, including Piscirickettsia and Yersinia, are almost all mosaic, while plasmids from other genera, including Borrelia, are rarely mosaic. While some mosaic plasmids share identical regions with hundreds of others, the median mosaic plasmid only shares with 8 other plasmids. When considering only plasmids from finished genomes (51.6% of the total), mosaic plasmids have significantly higher proportions of transposase and antibiotic resistance genes. Conversely, only 56.6% of mosaic fragments (DNA fragments shared between mosaic plasmids) contain a recognizable transposase gene, and only 1.2% of mosaic fragments are flanked by inverted repeats. Mosaic fragments associated with the IS26 transposase gene are 3.8-fold more abundant than any other sequence shared between mosaic plasmids in the database, though this is at least partly due to overrepresentation of Enterobacteriaceae plasmids. Mosaic plasmids are a complicated trait of some plasmid populations, only partly explained by transposition. Though antibiotic resistance genes led to the identification of many mosaic plasmids, mosaic plasmids are a broad phenomenon encompassing many more traits than just antibiotic resistance. Further research will be required to determine the influence of ecology, host repair mechanisms, conjugation, and plasmid host range on the formation and influence of mosaic plasmids. AUTHOR SUMMARY: Plasmids are extrachromosomal genetic entities that are found in many prokaryotes. They serve as flexible storage for genes, and individual cells can make substantial changes to their characteristics by acquiring, losing, or modifying a plasmid. In some pathogenic bacteria, such as Escherichia coli, antibiotic resistance genes are known to spread primarily on plasmids. By analyzing a database of 8592 plasmid sequences we determined that many of these plasmids have exchanged genes with each other, becoming mosaics of genes from different sources. We next separated these plasmids into groups based on the organism they were isolated from and found that different groups had different fractions of mosaic plasmids. This result was unexpected and suggests that the mechanisms and selective pressures causing mosaic plasmids do not occur evenly over all species. It also suggests that plasmids may provide different levels of potential variation to different species. This work uncovers a previously unrecognized pattern in plasmids across prokaryotes, that could lead to new insights into the evolutionary role that plasmids play.

Genomic Diversity in the Endosymbiotic Bacterium Rhizobium leguminosarum

Rhizobium leguminosarum bv. viciae is a soil α-proteobacterium that establishes a diazotrophic symbiosis with different legumes of the Fabeae tribe. The number of genome sequences from rhizobial strains available in public databases is constantly increasing, although complete, fully annotated genome structures from rhizobial genomes are scarce. In this work, we report and analyse the complete genome of R. leguminosarum bv. viciae UPM791. Whole genome sequencing can provide new insights into the genetic features contributing to symbiotically relevant processes such as bacterial adaptation to the rhizosphere, mechanisms for efficient competition with other bacteria, and the ability to establish a complex signalling dialogue with legumes, to enter the root without triggering plant defenses, and, ultimately, to fix nitrogen within the host. Comparison of the complete genome sequences of two strains of R. leguminosarum bv. viciae, 3841 and UPM791, highlights the existence of different symbiotic plasmids and a common core chromosome. Specific genomic traits, such as plasmid content or a distinctive regulation, define differential physiological capabilities of these endosymbionts. Among them, strain UPM791 presents unique adaptations for recycling the hydrogen generated in the nitrogen fixation process.

The Symbiosome: Legume and Rhizobia Co-evolution toward a Nitrogen-Fixing Organelle?

In legume nodules, symbiosomes containing endosymbiotic rhizobial bacteria act as temporary plant organelles that are responsible for nitrogen fixation, these bacteria develop mutual metabolic dependence with the host legume. In most legumes, the rhizobia infect post-mitotic cells that have lost their ability to divide, although in some nodules cells do maintain their mitotic capacity after infection. Here, we review what is currently known about legume symbiosomes from an evolutionary and developmental perspective, and in the context of the different interactions between diazotroph bacteria and eukaryotes. As a result, it can be concluded that the symbiosome possesses organelle-like characteristics due to its metabolic behavior, the composite origin and differentiation of its membrane, the retargeting of host cell proteins, the control of microsymbiont proliferation and differentiation by the host legume, and the cytoskeletal dynamics and symbiosome segregation during the division of rhizobia-infected cells. Different degrees of symbiosome evolution can be defined, specifically in relation to rhizobial infection and to the different types of nodule. Thus, our current understanding of the symbiosome suggests that it might be considered a nitrogen-fixing link in organelle evolution and that the distinct types of legume symbiosomes could represent different evolutionary stages toward the generation of a nitrogen-fixing organelle.

Ecological dynamics and complex interactions of Agrobacterium megaplasmids

As with many pathogenic bacteria, agrobacterial plant pathogens carry most of their virulence functions on a horizontally transmissible genetic element. The tumor-inducing (Ti) plasmid encodes the majority of virulence functions for the crown gall agent Agrobacterium tumefaciens. This includes the vir genes which drive genetic transformation of host cells and the catabolic genes needed to utilize the opines produced by infected plants. The Ti plasmid also encodes, an opine-dependent quorum sensing system that tightly regulates Ti plasmid copy number and its conjugal transfer to other agrobacteria. Many natural agrobacteria are avirulent, lacking the Ti plasmid. The burden of harboring the Ti plasmid depends on the environmental context. Away from diseased hosts, plasmid costs are low but the benefit of the plasmid is also absent. Consequently, plasmidless genotypes are favored. On infected plants the costs of the Ti plasmid can be very high, but balanced by the opine benefits, locally favoring plasmid bearing cells. Cheating derivatives which do not incur virulence costs but can benefit from opines are favored on infected plants and in most other environments, and these are frequently isolated from nature. Many agrobacteria also harbor an At plasmid which can stably coexist with a Ti plasmid. At plasmid genes are less well characterized but in general facilitate metabolic activities in the rhizosphere and bulk soil, such as the ability to breakdown plant exudates. Examination of A. tumefaciens C58, revealed that harboring its At plasmid is much more costly than harboring it's Ti plasmid, but conversely the At plasmid is extremely difficult to cure. The interactions between these co-resident plasmids are complex, and depend on environmental context. However, the presence of a Ti plasmid appears to mitigate At plasmid costs, consistent with the high frequency with which they are found together.

Genome-wide stochastic adaptive DNA amplification at direct and inverted DNA repeats in the parasite Leishmania

The human parasite Leishmania uses adaptive gene rearrangements and amplification involving repeated sequences on a genome-wide scale as one strategy to adapt to a changing environment.

Gene amplification of specific loci has been described in all kingdoms of life. In the protozoan parasite Leishmania, the product of amplification is usually part of extrachromosomal circular or linear amplicons that are formed at the level of direct or inverted repeated sequences. A bioinformatics screen revealed that repeated sequences are widely distributed in the Leishmania genome and the repeats are chromosome-specific, conserved among species, and generally present in low copy number. Using sensitive PCR assays, we provide evidence that the Leishmania genome is continuously being rearranged at the level of these repeated sequences, which serve as a functional platform for constitutive and stochastic amplification (and deletion) of genomic segments in the population. This process is adaptive as the copy number of advantageous extrachromosomal circular or linear elements increases upon selective pressure and is reversible when selection is removed. We also provide mechanistic insights on the formation of circular and linear amplicons through RAD51 recombinase-dependent and -independent mechanisms, respectively. The whole genome of Leishmania is thus stochastically rearranged at the level of repeated sequences, and the selection of parasite subpopulations with changes in the copy number of specific loci is used as a strategy to respond to a changing environment.

Variations in the copy number of DNA segments account for a substantial amount of genome diversity of most organisms. DNA amplification, a contributor to copy number variation, can occur in response to various stresses or after altered growth conditions, leading to extensive and often reversible genetic variation. DNA amplification in the parasite Leishmania occurs outside the normal chromosomes and arises by DNA rearrangements involving homologous repeated sequences. We show here that such repeated sequences are widespread in the Leishmania genome and that most of the Leishmania genome is subject to stochastic gene rearrangements mediated by these low-copy repeat sequences. Thus, although cells in the population have a common core genome, many individual cells will differ from the rest of the population by carrying one or more distinct extrachromosomal amplicon. Upon selection with either drugs or culture conditions, a subpopulation can emerge where the amplicon copy number per cell increases, and this clone of cells can then expand to dominate the population. We propose that Leishmania uses adaptive gene amplification at a genome-wide scale as one strategy to adapt to a changing environment.

Caffeine junkie: an unprecedented glutathione S-transferase-dependent oxygenase required for caffeine degradation by Pseudomonas putida CBB5

Caffeine and other N-methylated xanthines are natural products found in many foods, beverages, and pharmaceuticals. Therefore, it is not surprising that bacteria have evolved to live on caffeine as a sole carbon and nitrogen source. The caffeine degradation pathway of Pseudomonas putida CBB5 utilizes an unprecedented glutathione-S-transferase-dependent Rieske oxygenase for demethylation of 7-methylxanthine to xanthine, the final step in caffeine N-demethylation. The gene coding this function is unusual, in that the iron-sulfur and non-heme iron domains that compose the normally functional Rieske oxygenase (RO) are encoded by separate proteins. The non-heme iron domain is located in the monooxygenase, ndmC, while the Rieske [2Fe-2S] domain is fused to the RO reductase gene, ndmD. This fusion, however, does not interfere with the interaction of the reductase with N1- and N3-demethylase RO oxygenases, which are involved in the initial reactions of caffeine degradation. We demonstrate that the N7-demethylation reaction absolutely requires a unique, tightly bound protein complex composed of NdmC, NdmD, and NdmE, a novel glutathione-S-transferase (GST). NdmE is proposed to function as a noncatalytic subunit that serves a structural role in the complexation of the oxygenase (NdmC) and Rieske domains (NdmD). Genome analyses found this gene organization of a split RO and GST gene cluster to occur more broadly, implying a larger function for RO-GST protein partners.

The detection of inherent homologous recombination between repeat sequences in H. pylori 26695 by the PCR-based method

Helicobacter pylori infects more than half of the world's population, making it the most widespread infection of bacteria. It has high genetic diversity and has been considered as one of the most variable bacterial species. In the present study, a PCR-based method was used to detect the presence and the relative frequency of homologous recombination between repeat sequences (>500 bp) in H. pylori 26695. All the recombinant structures have been confirmed by sequencing. The inversion generated between inverted repeats showed distinct features from the recombination for duplication or deletion between direct repeats. Meanwhile, we gave the mathematic reasoning of a general formula for the calculation of relative recombination frequency and indicated the conditions for its application. This formula could be extensively applied to detect the frequency of homologous recombination, site-specific recombination, and other types of predictable recombination. Our results should be helpful for better understanding the genome evolution and adaptation of bacteria.

Natural genetic transformation generates a population of merodiploids in Streptococcus pneumoniae

Partial duplication of genetic material is prevalent in eukaryotes and provides potential for evolution of new traits. Prokaryotes, which are generally haploid in nature, can evolve new genes by partial chromosome duplication, known as merodiploidy. Little is known about merodiploid formation during genetic exchange processes, although merodiploids have been serendipitously observed in early studies of bacterial transformation. Natural bacterial transformation involves internalization of exogenous donor DNA and its subsequent integration into the recipient genome by homology. It contributes to the remarkable plasticity of the human pathogen Streptococcus pneumoniae through intra and interspecies genetic exchange. We report that lethal cassette transformation produced merodiploids possessing both intact and cassette-inactivated copies of the essential target gene, bordered by repeats (R) corresponding to incomplete copies of IS861. We show that merodiploidy is transiently stimulated by transformation, and only requires uptake of a ∼3-kb DNA fragment partly repeated in the chromosome. We propose and validate a model for merodiploid formation, providing evidence that tandem-duplication (TD) formation involves unequal crossing-over resulting from alternative pairing and interchromatid integration of R. This unequal crossing-over produces a chromosome dimer, resolution of which generates a chromosome with the TD and an abortive chromosome lacking the duplicated region. We document occurrence of TDs ranging from ∼100 to ∼900 kb in size at various chromosomal locations, including by self-transformation (transformation with recipient chromosomal DNA). We show that self-transformation produces a population containing many different merodiploid cells. Merodiploidy provides opportunities for evolution of new genetic traits via alteration of duplicated genes, unrestricted by functional selective pressure. Transient stimulation of a varied population of merodiploids by transformation, which can be triggered by stresses such as antibiotic treatment in S. pneumoniae, reinforces the plasticity potential of this bacterium and transformable species generally.

Merodiploids are defined as cells possessing a partial duplication of their genetic material, which potentially allows evolution of new genes. Historically, some have been observed in studies of natural genetic transformation. Transformation allows the bacteria to take up foreign DNA and incorporate it into their genome by homology. It is key to the high diversity observed in the human pathogen Streptococcus pneumoniae (the pneumococcus). Here we show that transformation with self DNA generates a population of merodiploids with varied chromosomal duplications, up to almost half a genome in size. We show that formation of merodiploids by transformation only requires uptake of a small donor DNA fragment partially repeated in the chromosome. The donor repeat recombines with an alternative repeat on one arm of a replicating chromosome, whilst the non-repeated part recombines with its complement on the other arm, bridging the two. Subsequent recombination events generate a merodiploid chromosome with the region between the two repeats duplicated. Our results demonstrate that transformation, which is induced by stresses such as antibiotic treatments, transiently increases the ability of a population to form merodiploids. We suggest that creating a variety of merodiploids at a time of stress maximizes the adaptive potential of this pathogen.

Exploration of noncoding sequences in metagenomes

Environment-dependent genomic features have been defined for different metagenomes, whose genes and their associated processes are related to specific environments. Identification of ORFs and their functional categories are the most common methods for association between functional and environmental features. However, this analysis based on finding ORFs misses noncoding sequences and, therefore, some metagenome regulatory or structural information could be discarded. In this work we analyzed 23 whole metagenomes, including coding and noncoding sequences using the following sequence patterns: (G+C) content, Codon Usage (Cd), Trinucleotide Usage (Tn), and functional assignments for ORF prediction. Herein, we present evidence of a high proportion of noncoding sequences discarded in common similarity-based methods in metagenomics, and the kind of relevant information present in those. We found a high density of trinucleotide repeat sequences (TRS) in noncoding sequences, with a regulatory and adaptive function for metagenome communities. We present associations between trinucleotide values and gene function, where metagenome clustering correlate with microorganism adaptations and kinds of metagenomes. We propose here that noncoding sequences have relevant information to describe metagenomes that could be considered in a whole metagenome analysis in order to improve their organization, classification protocols, and their relation with the environment.

Multiple pathways of duplication formation with and without recombination (RecA) in Salmonella enterica

Duplications are often attributed to "unequal recombination" between separated, directly repeated sequence elements (>100 bp), events that leave a recombinant element at the duplication junction. However, in the bacterial chromosome, duplications form at high rates (10(-3)-10(-5)/cell/division) even without recombination (RecA). Here we describe 1800 spontaneous lac duplications trapped nonselectively on the low-copy F'(128) plasmid, where lac is flanked by direct repeats of the transposable element IS3 (1258 bp) and by numerous quasipalindromic REP elements (30 bp). Duplications form at a high rate (10(-4)/cell/division) that is reduced only about 11-fold in the absence of RecA. With and without RecA, most duplications arise by recombination between IS3 elements (97%). Formation of these duplications is stimulated by IS3 transposase (Tnp) and plasmid transfer functions (TraI). Three duplication pathways are proposed. First, plasmid dimers form at a high rate stimulated by RecA and are then modified by deletions between IS3 elements (resolution) that leave a monomeric plasmid with an IS3-flanked lac duplication. Second, without RecA, duplications occur by single-strand annealing of DNA ends generated in different sister chromosomes after transposase nicks DNA near participating IS3 elements. The absence of RecA may stimulate annealing by allowing chromosome breaks to persist. Third, a minority of lac duplications (3%) have short (0-36 bp) junction sequences (SJ), some of which are located within REP elements. These duplication types form without RecA, Tnp, or Tra by a pathway in which the palindromic junctions of a tandem inversion duplication (TID) may stimulate deletions that leave the final duplication.

Rhizobial plasmids — replication, structure and biological role

Soil bacteria, collectively named rhizobia, can establish mutualistic relationships with legume plants. Rhizobia often have multipartite genome architecture with a chromosome and several extrachromosomal replicons making these bacteria a perfect candidate for plasmid biology studies. Rhizobial plasmids are maintained in the cells using a tightly controlled and uniquely organized replication system. Completion of several rhizobial genome-sequencing projects has changed the view that their genomes are simply composed of the chromosome and cryptic plasmids. The genetic content of plasmids and the presence of some important (or even essential) genes contribute to the capability of environmental adaptation and competitiveness with other bacteria. On the other hand, their mosaic structure results in the plasticity of the genome and demonstrates a complex evolutionary history of plasmids. In this review, a genomic perspective was employed for discussion of several aspects regarding rhizobial plasmids comprising structure, replication, genetic content, and biological role. A special emphasis was placed on current post-genomic knowledge concerning plasmids, which has enriched the view of the entire bacterial genome organization by the discovery of plasmids with a potential chromosome-like role.

Rhizobial extrachromosomal replicon variability, stability and expression in natural niches

In bacteria, niche adaptation may be determined by mobile extrachromosomal elements. A remarkable characteristic of Rhizobium and Ensifer (Sinorhizobium) but also of Agrobacterium species is that almost half of the genome is contained in several large extrachromosomal replicons (ERs). They encode a plethora of functions, some of them required for bacterial survival, niche adaptation, plasmid transfer or stability. In spite of this, plasmid loss is common in rhizobia upon subculturing. Rhizobial gene-expression studies in plant rhizospheres with novel results from transcriptomic analysis of Rhizobium phaseoli in maize and Phaseolus vulgaris roots highlight the role of ERs in natural niches and allowed the identification of common extrachromosomal genes expressed in association with plant rootlets and the replicons involved.

Diversity and biogeography of rhizobia isolated from root nodules of Glycine max grown in Hebei Province, China

A total of 215 rhizobial strains were isolated and analyzed with 16S rRNA gene, 16S-23S intergenic spacer, housekeeping genes atpD, recA, and glnII, and symbiotic genes nifH and nodC to understand the genetic diversity of soybean rhizobia in Hebei province, China. All the strains except one were symbiotic bacteria classified into nine genospecies in the genera of Bradyrhizobium and Sinorhizobium. Surveys on the distribution of these rhizobia in different regions showed that Bradyrhizobium japonicum and Bradyrhizobium elkanii strains were found only in neutral to slightly alkaline soils whereas Bradyrhizobium yuanmingense, Bradyrhizobium liaoningense-related strains and strains of five Sinorhizobium genospecies were found in alkaline-saline soils. Correspondence and canonical correspondence analyses on the relationship of rhizobial distribution and their soil characteristics reveal that high soil pH, electrical conductivity, and potassium content favor distribution of the B. yuanmingense and the five Sinorhizobium species but inhibit B. japonicum and B. elkanii. High contents of available phosphorus and organic matters benefit Sinorhizobium fredii and B. liaoningense-related strains and inhibit the others groups mentioned above. The symbiotic gene (nifH and nodC) lineages among B. elkanii, B. japonicum, B. yuanmingense, and Sinorhizobium spp. were observed in the strains, signifying that vertical gene transfer was the main mechanism to maintain these genes in the soybean rhizobia. However, lateral transfer of symbiotic genes commonly in Sinorhizobium spp. and rarely in Bradyrhizobium spp. was also detected. These results showed the genetic diversity, the biogeography, and the soil determinant factors of soybean rhizobia in Hebei province of China.

Mechanism of chromosomal transfer of Enterococcus faecalis pathogenicity island, capsule, antimicrobial resistance, and other traits

The Enterococcus faecalis pathogenicity island (PAI) encodes known virulence traits and >100 additional genes with unknown roles in enterococcal biology. Phage-related integration and excision genes, and direct repeats flanking the island, suggest it moves as an integrative conjugative element (ICE). However, transfer was observed not to require these genes. Transfer only occurred from donors possessing a pheromone responsive-type of conjugative plasmid, and was invariably accompanied by transfer of flanking donor chromosome sequences. Deletion of plasmid transfer functions, including the cis-acting origin of transfer (oriT), abolished movement. In addition to demonstrating PAI movement by a mechanism involving plasmid integration, we observed transfer of a selectable marker placed virtually anywhere on the chromosome. Transfer of this selectable marker was observed to be accompanied by chromosome-chromosome transfer of vancomycin resistance, MLST markers, and capsule genes as well. Plasmid mobilization therefore appears to be a major mechanism for horizontal gene transfer in the evolution of antibiotic resistant E. faecalis strains capable of causing human infection.

Recovery of nonpathogenic mutant bacteria from tumors caused by several Agrobacterium tumefaciens strains: a frequent event?

We have evaluated the interaction that bacterial genotypes and plant hosts have with the loss of pathogenicity in tumors, using seven Agrobacterium tumefaciens strains inoculated on 12 herbaceous and woody hosts. We performed a screening of the agrobacteria present inside the tumors, looking for nonpathogenic strains, and found a high variability of those strains in this niche. To verify the origin of the putative nonpathogenic mutant bacteria, we applied an efficient, reproducible, and specific randomly amplified polymorphic DNA analysis method. In contrast with previous studies, we recovered a very small percentage (0.01%) of nonpathogenic strains that can be considered true mutants. Of 5,419 agrobacterial isolates examined, 662 were nonpathogenic in tomato, although only 7 (from pepper and tomato tumors induced by two A. tumefaciens strains) could be considered to derive from the inoculated strain. Six mutants were affected in the transferred DNA (T-DNA) region; one of them contained IS426 inserted into the iaaM gene, whereas the whole T-DNA region was apparently deleted in three other mutants, and the virulence of the remaining two mutants was fully restored with the T-DNA genes as well. The plasmid profile was altered in six of the mutants, with changes in the size of the Ti plasmid or other plasmids and/or the acquisition of new plasmids. Our results also suggest that the frequent occurrence of nonpathogenic clones in the tumors is probably due to the preferential growth of nonpathogenic agrobacteria, of either endophytic or environmental origin, but different from the bacterial strain inducing the tumor.

Plastic architecture of bacterial genome revealed by comparative genomics of Photorhabdus variants

BACKGROUND

The phenotypic consequences of large genomic architecture modifications within a clonal bacterial population are rarely evaluated because of the difficulties associated with using molecular approaches in a mixed population. Bacterial variants frequently arise among Photorhabdus luminescens, a nematode-symbiotic and insect-pathogenic bacterium. We therefore studied genome plasticity within Photorhabdus variants.

RESULTS

We used a combination of macrorestriction and DNA microarray experiments to perform a comparative genomic study of different P. luminescens TT01 variants. Prolonged culturing of TT01 strain and a genomic variant, collected from the laboratory-maintained symbiotic nematode, generated bacterial lineages composed of primary and secondary phenotypic variants and colonial variants. The primary phenotypic variants exhibit several characteristics that are absent from the secondary forms. We identify substantial plasticity of the genome architecture of some variants, mediated mainly by deletions in the 'flexible' gene pool of the TT01 reference genome and also by genomic amplification. We show that the primary or secondary phenotypic variant status is independent from global genomic architecture and that the bacterial lineages are genomic lineages. We focused on two unusual genomic changes: a deletion at a new recombination hotspot composed of long approximate repeats; and a 275 kilobase single block duplication belonging to a new class of genomic duplications.

CONCLUSION

Our findings demonstrate that major genomic variations occur in Photorhabdus clonal populations. The phenotypic consequences of these genomic changes are cryptic. This study provides insight into the field of bacterial genome architecture and further elucidates the role played by clonal genomic variation in bacterial genome evolution.

Recurrent DNA inversion rearrangements in the human genome

Several lines of evidence suggest that reiterated sequences in the human genome are targets for nonallelic homologous recombination (NAHR), which facilitates genomic rearrangements. We have used a PCR-based approach to identify breakpoint regions of rearranged structures in the human genome. In particular, we have identified intrachromosomal identical repeats that are located in reverse orientation, which may lead to chromosomal inversions. A bioinformatic workflow pathway to select appropriate regions for analysis was developed. Three such regions overlapping with known human genes, located on chromosomes 3, 15, and 19, were analyzed. The relative proportion of wild-type to rearranged structures was determined in DNA samples from blood obtained from different, unrelated individuals. The results obtained indicate that recurrent genomic rearrangements occur at relatively high frequency in somatic cells. Interestingly, the rearrangements studied were significantly more abundant in adults than in newborn individuals, suggesting that such DNA rearrangements might start to appear during embryogenesis or fetal life and continue to accumulate after birth. The relevance of our results in regard to human genomic variation is discussed.

Molecular and population analyses of a recombination event in the catabolic plasmid pJP4

Cupriavidus necator JMP134(pJP4) harbors a catabolic plasmid, pJP4, which confers the ability to grow on chloroaromatic compounds. Repeated growth on 3-chlorobenzoate (3-CB) results in selection of a recombinant strain, which degrades 3-CB better but no longer grows on 2,4-dichlorophenoxyacetate (2,4-D). We have previously proposed that this phenotype is due to a double homologous recombination event between inverted repeats of the multicopies of this plasmid within the cell. One recombinant form of this plasmid (pJP4-F3) explains this phenotype, since it harbors two copies of the chlorocatechol degradation tfd gene clusters, which are essential to grow on 3-CB, but has lost the tfdA gene, encoding the first step in degradation of 2,4-D. The other recombinant plasmid (pJP4-FM) should harbor two copies of the tfdA gene but no copies of the tfd gene clusters. A molecular analysis using a multiplex PCR approach to distinguish the wild-type plasmid pJP4 from its two recombinant forms, was carried out. Expected PCR products confirming this recombination model were found and sequenced. Few recombinant plasmid forms in cultures grown in several carbon sources were detected. Kinetic studies indicated that cells containing the recombinant plasmid pJP4-FM were not selectable by sole carbon source growth pressure, whereas those cells harboring recombinant plasmid pJP4-F3 were selected upon growth on 3-CB. After 12 days of repeated growth on 3-CB, the complete plasmid population in C. necator JMP134 apparently corresponds to this form. However, wild-type plasmid forms could be recovered after growing this culture on 2,4-D, indicating that different plasmid forms can be found in C. necator JMP134 at the population level.

Isolation and characterization of functional insertion sequences of rhizobia

Rhizobia are a group of bacteria that form nodules on the roots of legume host plants. The sequenced genomes of the rhizobia are characterized by the presence of many putative insertion sequences (IS) elements. However, it is unknown whether these IS elements are functional and it is therefore relevant to assess their transposition activity. In this work, several functional insertion sequences belonging to the IS1256, IS3, IS5, IS166, and IS21 families were captured from Rhizobium tropici, Rhizobium sp. NGR234 and Sinorhizobium meliloti, using pGBG1 as a trapping system. In silico analysis shows that homologs of rhizobia mobile elements are present in distantly related genomes, suggesting that Rhizobium IS elements are prone to genetic transfer.

The partitioned Rhizobium etli genome: genetic and metabolic redundancy in seven interacting replicons

We report the complete 6,530,228-bp genome sequence of the symbiotic nitrogen fixing bacterium Rhizobium etli. Six large plasmids comprise one-third of the total genome size. The chromosome encodes most functions necessary for cell growth, whereas few essential genes or complete metabolic pathways are located in plasmids. Chromosomal synteny is disrupted by genes related to insertion sequences, phages, plasmids, and cell-surface components. Plasmids do not show synteny, and their orthologs are mostly shared by accessory replicons of species with multipartite genomes. Some nodulation genes are predicted to be functionally related with chromosomal loci encoding for the external envelope of the bacterium. Several pieces of evidence suggest an exogenous origin for the symbiotic plasmid (p42d) and p42a. Additional putative horizontal gene transfer events might have contributed to expand the adaptive repertoire of R. etli, because they include genes involved in small molecule metabolism, transport, and transcriptional regulation. Twenty-three putative sigma factors, numerous isozymes, and paralogous families attest to the metabolic redundancy and the genomic plasticity necessary to sustain the lifestyle of R. etli in symbiosis and in the soil.

Diversification of DNA sequences in the symbiotic genome of Rhizobium etli

Bacteria of the genus Rhizobium and related genera establish nitrogen-fixing symbioses with the roots of leguminous plants. The genetic elements that participate in the symbiotic process are usually compartmentalized in the genome, either as independent replicons (symbiotic plasmids) or as symbiotic regions or islands in the chromosome. The complete nucleotide sequence of the symbiotic plasmid of Rhizobium etli model strain CFN42, symbiont of the common bean plant, has been reported. To better understand the basis of DNA sequence diversification of this symbiotic compartment, we analyzed the distribution of single-nucleotide polymorphisms in homologous regions from different Rhizobium etli strains. The distribution of polymorphisms is highly asymmetric in each of the different strains, alternating regions containing very few changes with regions harboring an elevated number of substitutions. The regions showing high polymorphism do not correspond with discrete genetic elements and are not the same in the different strains, indicating that they are not hypervariable regions of functional genes. Most interesting, some highly polymorphic regions share exactly the same nucleotide substitutions in more than one strain. Furthermore, in different regions of the symbiotic compartment, different sets of strains share the same substitutions. The data indicate that the majority of nucleotide substitutions are spread in the population by recombination and that the contribution of new mutations to polymorphism is relatively low. We propose that the horizontal transfer of homologous DNA segments among closely related organisms is a major source of genomic diversification.

Genome reduction in the α-Proteobacteria

More than 20 alpha-proteobacterial genomes are currently available. These range in size from 1-9 Mb and represent excellent model systems for evolutionary studies of the organizational features of bacterial genomes. Computational inferences have shown that genome reductions have occurred independently in lineages such as Rickettsia and Bartonella that are associated with intracellular lifestyles. Analyses of these reduced genomes have provided insights into the evolution of vector-borne transmission pathways. Further research into the population biology of bacteria, arthropods and vertebrate hosts will help to refine the biology of host-pathogen interactions and will facilitate the design of vaccines and vector-control programs.

Mechanisms of, and barriers to, horizontal gene transfer between bacteria

Bacteria evolve rapidly not only by mutation and rapid multiplication, but also by transfer of DNA, which can result in strains with beneficial mutations from more than one parent. Transformation involves the release of naked DNA followed by uptake and recombination. Homologous recombination and DNA-repair processes normally limit this to DNA from similar bacteria. However, if a gene moves onto a broad-host-range plasmid it might be able to spread without the need for recombination. There are barriers to both these processes but they reduce, rather than prevent, gene acquisition.

Selection for gene clustering by tandem duplication

In prokaryotic genomes, related genes are frequently clustered in operons and higher-order arrangements that reflect functional context. Organization emerges despite rearrangements that constantly shuffle gene and operon order. Evidence is presented that the tandem duplication of related genes acts as a driving evolutionary force in the origin and maintenance of clusters. Gene amplification can be viewed as a dynamic and reversible regulatory mechanism that facilitates adaptation to variable environments. Clustered genes confer selective benefits via their ability to be coamplified. During evolution, rearrangements that bring together related genes can be selected if they increase the fitness of the organism in which they reside. Similarly, the benefits of gene amplification can prevent the dispersal of existing clusters. Examples of frequent and spontaneous amplification of large genomic fragments are provided. The possibility is raised that tandem gene duplication works in concert with horizontal gene transfer as interrelated evolutionary forces for gene clustering.

Combining mathematical models and statistical methods to understand and predict the dynamics of antibiotic-sensitive mutants in a population of resistant bacteria during experimental evolution

Temporarily discontinuing the use of antibiotics has been proposed as a means to eliminate resistant bacteria by allowing sensitive clones to sweep through the population. In this study, we monitored a tetracycline-sensitive subpopulation that emerged during experimental evolution of E. coli K12 MG1655 carrying the multiresistance plasmid pB10 in the absence of antibiotics. The fraction of tetracycline-sensitive mutants increased slowly over 500 generations from 0.1 to 7%, and loss of resistance could be attributed to a recombination event that caused deletion of the tet operon. To help understand the population dynamics of these mutants, three mathematical models were developed that took into consideration recurrent mutations, increased host fitness (selection), or a combination of both mechanisms (full model). The data were best explained by the full model, which estimated a high mutation frequency (lambda = 3.11 x 10(-5)) and a significant but small selection coefficient (sigma = 0.007). This study emphasized the combined use of experimental data, mathematical models, and statistical methods to better understand and predict the dynamics of evolving bacterial populations, more specifically the possible consequences of discontinuing the use of antibiotics.

Flavonoids induce temporal shifts in gene-expression of nod-box controlled loci in Rhizobium sp. NGR234

Rhizobia, soil bacteria of the Rhizobiales, enter the roots of homologous legumes, where they induce the formation of nitrogen-fixing nodules. Signals emanating from both symbiotic partners control nodule development. Efficient nodulation requires precise, temporal regulation of symbiotic genes. Roots continuously release flavonoids that interact with transcriptional activators of the LysR family. NodD proteins, which are members of this family, act both as sensors of the environment and modulate the expression of genes preceded by conserved promoter sequences called nod-boxes. The symbiotic plasmid of the broad host-range Rhizobium sp. NGR234 caries 19 nod-boxes (NB1 to NB19), all of which were cloned upstream of a lacZ-reporter gene. A flavonoid, daidzein was able to induce 18 of the 19 nod-boxes in a NodD1-dependent manner. Interestingly, induction of four nod-boxes (NB6, NB15, NB16 and NB17) is highly dependent on NodD2 and was delayed in comparison with the others. In turn, NodD2 is involved in the repression of the NB8 nodABCIJnolOnoeI operon. Activation of transcription of nodD2 is also dependent on flavonoids despite the absence of a nod-box like sequence in the upstream promoter region. Mutational analysis showed that syrM 2 (another member of the LysR family), which is controlled by NB19, is also necessary for expression of nodD 2. Thus, NodD1, NodD2 and SyrM2 co-modulate a flavonoid-inducible regulatory cascade that coordinates the expression of symbiotic genes with nodule development.

DNA amplification and rearrangements in archival Salmonella enterica serovar Typhimurium LT2 cultures

Variations in genome size and gene order were observed in archival Salmonella enterica serovar Typhimurium cultures stored for over 40 years. In one strain, microarray analysis revealed a large, stable amplification. PCR analysis of the same strain revealed a genomic duplication that underwent a translocation. Other strains had smaller duplications and deletions. These results demonstrate that storage in stabs over time at room temperature not only allows for further bacterial growth but also may produce an environment that selects for a variety of mutations, including genomic rearrangements.

Natural genomic design in Sinorhizobium meliloti: novel genomic architectures

The complete nucleotide sequence of the genome of Sinorhizobium meliloti, the symbiont of alfalfa, was reported in 2001 by an international consortium of laboratories. The genome comprises a chromosome of 3.65 megabases (Mb) and two megaplasmids, pSymA and pSymB, of 1.35 Mb and 1.68 Mb, respectively. Based on the nucleotide sequence of the whole genome, we designed a pathway of consecutive rearrangements leading to novel genomic architectures. In a first step we obtained derivative strains containing two replicons; in a second step we obtained a strain containing the genetic information in one single replicon of 6.68 MB. From this last architecture we isolated revertants containing two replicons, and from these we could return to the original architecture showing the three replicons. We found that the relative frequency of excision of cointegrated replicons is higher at the site used for the cointegration than at other sites. This might conciliate two apparently opposed facts: the highly dynamic state of genomic architecture in S. meliloti and the common observation that different isolates and derived cellular clones of S. meliloti usually present the architecture of one chromosome and two distinct megaplasmids. Different aspects that must be considered to obtain full advantage of the strategy of natural genomic design are discussed.

Rhizobium etli and Rhizobium gallicum nodulate common bean (Phaseolus vulgaris) in a traditionally managed milpa plot in Mexico: population genetics and biogeographic implications

The stability of the genetic structure of rhizobial populations nodulating Phaseolus vulgaris cultivated in a traditionally managed milpa plot in Mexico was studied over three consecutive years. The set of molecular markers analyzed (including partial rrs, glnII, nifH, and nodB sequences), along with host range experiments, placed the isolates examined in Rhizobium etli bv. phaseoli and Rhizobium gallicum bv. gallicum. Cluster analysis of multilocus enzyme electrophoresis and plasmid profile data separated the two species and identified numerically dominant clones within each of them. Population genetic analyses showed that there was high genetic differentiation between the two species and that there was low intrapopulation differentiation of the species over the 3 years. The results of linkage disequilibrium analyses are consistent with an epidemic genetic structure for both species, with frequent genetic exchange taking place within conspecific populations but not between the R. etli and R. gallicum populations. A subsample of isolates was selected and used for 16S ribosomal DNA PCR-restriction fragment length polymorphism analysis, nifH copy number determination, and host range experiments. Plasmid profiles and nifH hybridization patterns also revealed the occurrence of lateral plasmid transfer among distinct multilocus genotypes within species but not between species. Both species were recovered from nodules of the same plants, indicating that mechanisms other than host, spatial, or temporal isolation may account for the genetic barrier between the species. The biogeographic implications of finding an R. gallicum bv. gallicum population nodulating common bean in America are discussed.

The mosaic structure of the symbiotic plasmid of Rhizobium etli CFN42 and its relation to other symbiotic genome compartments

BACKGROUND

Symbiotic bacteria known as rhizobia interact with the roots of legumes and induce the formation of nitrogen-fixing nodules. In rhizobia, essential genes for symbiosis are compartmentalized either in symbiotic plasmids or in chromosomal symbiotic islands. To understand the structure and evolution of the symbiotic genome compartments (SGCs), it is necessary to analyze their common genetic content and organization as well as to study their differences. To date, five SGCs belonging to distinct species of rhizobia have been entirely sequenced. We report the complete sequence of the symbiotic plasmid of Rhizobium etli CFN42, a microsymbiont of beans, and a comparison with other SGC sequences available.

RESULTS

The symbiotic plasmid is a circular molecule of 371,255 base-pairs containing 359 coding sequences. Nodulation and nitrogen-fixation genes common to other rhizobia are clustered in a region of 125 kilobases. Numerous sequences related to mobile elements are scattered throughout. In some cases the mobile elements flank blocks of functionally related sequences, thereby suggesting a role in transposition. The plasmid contains 12 reiterated DNA families that are likely to participate in genomic rearrangements. Comparisons between this plasmid and complete rhizobial genomes and symbiotic compartments already sequenced show a general lack of synteny and colinearity, with the exception of some transcriptional units. There are only 20 symbiotic genes that are shared by all SGCs.

CONCLUSIONS

Our data support the notion that the symbiotic compartments of rhizobia genomes are mosaic structures that have been frequently tailored by recombination, horizontal transfer and transposition.

Dynamics of genome architecture in Rhizobium sp. strain NGR234

Bacterial genomes are usually partitioned in several replicons, which are dynamic structures prone to mutation and genomic rearrangements, thus contributing to genome evolution. Nevertheless, much remains to be learned about the origins and dynamics of the formation of bacterial alternative genomic states and their possible biological consequences. To address these issues, we have studied the dynamics of the genome architecture in Rhizobium sp. strain NGR234 and analyzed its biological significance. NGR234 genome consists of three replicons: the symbiotic plasmid pNGR234a (536,165 bp), the megaplasmid pNGR234b (>2,000 kb), and the chromosome (>3,700 kb). Here we report that genome analyses of cell siblings showed the occurrence of large-scale DNA rearrangements consisting of cointegrations and excisions between the three replicons. As a result, four new genomic architectures have emerged. Three consisted of the cointegrates between two replicons: chromosome-pNGR234a, chromosome-pNGR234b, and pNGR234a-pNGR234b. The other consisted of a cointegrate of the three replicons (chromosome-pNGR234a-pNGR234b). Cointegration and excision of pNGR234a with either the chromosome or pNGR234b were studied and found to proceed via a Campbell-type mechanism, mediated by insertion sequence elements. We provide evidence showing that changes in the genome architecture did not alter the growth and symbiotic proficiency of Rhizobium derivatives.

Séquences répétées des génomes de Rhizobium sp. NGR234 et Sinorhizobium meliloti : une analyse comparative par séquençage aléatoire

Amongst prokaryotic genomes, those of nitrogen-fixing members of the Rhizobiaceae family are relatively large (6–9 Mb), often include mega-plasmids of 1.5–2 Mb, and contain numerous families of repeated DNA sequences. Although most essential nodulation and nitrogen fixation genes are well characterized, these represent only a small fraction of the DNA content. Little is known about the detailed structure of rhizobial genomes. With the development of sequencing techniques and new bio-informatic tools such studies become possible, however. Using the 2275 shot-gun sequences of ANU265 (a derivative of NGR234 cured of pNGR234a), we have identified numerous families of repeats. Amongst these, the 58-bp-long NGRREP-4 represents the third most abundant DNA sequence after the RIME1 and RIME2 repeats, all of which are also found in Sinorhizobium meliloti. Surprisingly, studies on the distribution of these elements showed that in proportion to its size, the chromosome of NGR234 carries many more RIME modules than pNGR234a or pNGR234b. Together with the presence in NGR234 and S. meliloti 1021 of an insertion sequence (IS) element more conserved than essential nodulation and nitrogen fixation genes, these results give new insights into the origin and evolution of rhizobial genomes.Key words: shot-gun, repeats, BIME.