Gene duplication in bacteria: alteration of gene dosage by sister-chromosome exchanges

Three concurrent mechanisms generate gene copy number variation and transient antibiotic heteroresistance

Heteroresistance is a medically relevant phenotype where small antibiotic-resistant subpopulations coexist within predominantly susceptible bacterial populations. Heteroresistance reduces treatment efficacy across diverse bacterial species and antibiotic classes, yet its genetic and physiological mechanisms remain poorly understood. Here, we investigated a multi-resistant Klebsiella pneumoniae isolate and identified three primary drivers of gene dosage-dependent heteroresistance for several antibiotic classes: tandem amplification, increased plasmid copy number, and transposition of resistance genes onto cryptic plasmids. All three mechanisms imposed fitness costs and were genetically unstable, leading to fast reversion to susceptibility in the absence of antibiotics. We used a mouse gut colonization model to show that heteroresistance due to elevated resistance-gene dosage can result in antibiotic treatment failures. Importantly, we observed that the three mechanisms are prevalent among Escherichia coli bloodstream isolates. Our findings underscore the necessity for treatment strategies that address the complex interplay between plasmids, resistance cassettes, and transposons in bacterial populations.

Large-scale duplication events underpin population-level flexibility in tRNA gene copy number in Pseudomonas fluorescens SBW25

The complement of tRNA genes within a genome is typically considered to be a (relatively) stable characteristic of an organism. Here, we demonstrate that bacterial tRNA gene set composition can be more flexible than previously appreciated, particularly regarding tRNA gene copy number. We report the high-rate occurrence of spontaneous, large-scale, tandem duplication events in laboratory populations of the bacterium Pseudomonas fluorescens SBW25. The identified duplications are up to ∼1 Mb in size (∼15% of the wildtype genome) and are predicted to change the copy number of up to 917 genes, including several tRNA genes. The observed duplications are inherently unstable: they occur, and are subsequently lost, at extremely high rates. We propose that this unusually plastic type of mutation provides a mechanism by which tRNA gene set diversity can be rapidly generated, while simultaneously preserving the underlying tRNA gene set in the absence of continued selection. That is, if a tRNA set variant provides no fitness advantage, then high-rate segregation of the duplication ensures the maintenance of the original tRNA gene set. However, if a tRNA gene set variant is beneficial, the underlying duplication fragment(s) may persist for longer and provide raw material for further, more stable, evolutionary change.

Evolution of new functions de novo and from preexisting genes

How the enormous structural and functional diversity of new genes and proteins was generated (estimated to be 10(10)-10(12) different proteins in all organisms on earth [Choi I-G, Kim S-H. 2006. Evolution of protein structural classes and protein sequence families. Proc Natl Acad Sci 103: 14056-14061] is a central biological question that has a long and rich history. Extensive work during the last 80 years have shown that new genes that play important roles in lineage-specific phenotypes and adaptation can originate through a multitude of different mechanisms, including duplication, lateral gene transfer, gene fusion/fission, and de novo origination. In this review, we focus on two main processes as generators of new functions: evolution of new genes by duplication and divergence of pre-existing genes and de novo gene origination in which a whole protein-coding gene evolves from a noncoding sequence.

GenoBase: comprehensive resource database of Escherichia coli K-12

Comprehensive experimental resources, such as ORFeome clone libraries and deletion mutant collections, are fundamental tools for elucidation of gene function. Data sets by omics analysis using these resources provide key information for functional analysis, modeling and simulation both in individual and systematic approaches. With the long-term goal of complete understanding of a cell, we have over the past decade created a variety of clone and mutant sets for functional genomics studies of Escherichia coli K-12. We have made these experimental resources freely available to the academic community worldwide. Accordingly, these resources have now been used in numerous investigations of a multitude of cell processes. Quality control is extremely important for evaluating results generated by these resources. Because the annotation has been changed since 2005, which we originally used for the construction, we have updated these genomic resources accordingly. Here, we describe GenoBase (http://ecoli.naist.jp/GB/), which contains key information about comprehensive experimental resources of E. coli K-12, their quality control and several omics data sets generated using these resources.

High fitness costs and instability of gene duplications reduce rates of evolution of new genes by duplication-divergence mechanisms

An important mechanism for generation of new genes is by duplication-divergence of existing genes. Duplication-divergence includes several different submodels, such as subfunctionalization where after accumulation of neutral mutations the original function is distributed between two partially functional and complementary genes, and neofunctionalization where a new function evolves in one of the duplicated copies while the old function is maintained in another copy. The likelihood of these mechanisms depends on the longevity of the duplicated state, which in turn depends on the fitness cost and genetic stability of the duplications. Here, we determined the fitness cost and stability of defined gene duplications/amplifications on a low copy number plasmid. Our experimental results show that the costs of carrying extra gene copies are substantial and that each additional kilo base pairs of DNA reduces fitness by approximately 0.15%. Furthermore, gene amplifications are highly unstable and rapidly segregate to lower copy numbers in absence of selection. Mathematical modeling shows that the fitness costs and instability strongly reduces the likelihood of both sub- and neofunctionalization, but that these effects can be offset by positive selection for novel beneficial functions.

Genomic plasticity enables phenotypic variation of Pseudomonas syringae pv. tomato DC3000

Whole genome sequencing revealed the presence of a genomic anomaly in the region of 4.7 to 4.9 Mb of the Pseudomonas syringae pv. tomato (Pst) DC3000 genome. The average read depth coverage of Pst DC3000 whole genome sequencing results suggested that a 165 kb segment of the chromosome had doubled in copy number. Further analysis confirmed the 165 kb duplication and that the two copies were arranged as a direct tandem repeat. Examination of the corresponding locus in Pst NCPPB1106, the parent strain of Pst DC3000, suggested that the 165 kb duplication most likely formed after the two strains diverged via transposition of an ISPsy5 insertion sequence (IS) followed by unequal crossing over between ISPsy5 elements at each end of the duplicated region. Deletion of one copy of the 165 kb region demonstrated that the duplication facilitated enhanced growth in some culture conditions, but did not affect pathogenic growth in host tomato plants. These types of chromosomal structures are predicted to be unstable and we have observed resolution of the 165 kb duplication to single copy and its subsequent re-duplication. These data demonstrate the role of IS elements in recombination events that facilitate genomic reorganization in P. syringae.

Dynamics of gene duplication in the genomes of chlorophyll d-producing cyanobacteria: implications for the ecological niche

Gene duplication may be an important mechanism for the evolution of new functions and for the adaptive modulation of gene expression via dosage effects. Here, we analyzed the fate of gene duplicates for two strains of a novel group of cyanobacteria (genus Acaryochloris) that produces the far-red light absorbing chlorophyll d as its main photosynthetic pigment. The genomes of both strains contain an unusually high number of gene duplicates for bacteria. As has been observed for eukaryotic genomes, we find that the demography of gene duplicates can be well modeled by a birth-death process. Most duplicated Acaryochloris genes are of comparatively recent origin, are strain-specific, and tend to be located on different genetic elements. Analyses of selection on duplicates of different divergence classes suggest that a minority of paralogs exhibit near neutral evolutionary dynamics immediately following duplication but that most duplicate pairs (including those which have been retained for long periods) are under strong purifying selection against amino acid change. The likelihood of duplicate retention varied among gene functional classes, and the pronounced differences between strains in the pool of retained recent duplicates likely reflects differences in the nutrient status and other characteristics of their respective environments. We conclude that most duplicates are quickly purged from Acaryochloris genomes and that those which are retained likely make important contributions to organism ecology by conferring fitness benefits via gene dosage effects. The mechanism of enhanced duplication may involve homologous recombination between genetic elements mediated by paralogous copies of recA.

DNA topology of highly transcribed operons in Salmonella enterica serovar Typhimurium

Bacteria differ from eukaryotes by having the enzyme DNA gyrase, which catalyses the ATP-dependent negative supercoiling of DNA. Negative supercoils are essential for condensing chromosomes into an interwound (plectonemic) and branched structure known as the nucleoid. Topo-1 removes excess supercoiling in an ATP-independent reaction and works with gyrase to establish a topological equilibrium where supercoils move within 10 kb domains bounded by stochastic barriers along the sequence. However, transcription changes the stochastic pattern by generating supercoil diffusion barriers near the sites of gene expression. Using supercoil-dependent Tn3 and γδ resolution assays, we studied DNA topology upstream, downstream and across highly transcribed operons. Whenever two Res sites flanked efficiently transcribed genes, resolution was inhibited and the loss in recombination efficiency was proportional to transcription level. Ribosomal RNA operons have the highest transcription rates, and resolution assays at the rrnG and rrnH operons showed inhibitory levels 40-100 times those measured in low-transcription zones. Yet, immediately upstream and downstream of RNA polymerase (RNAP) initiation and termination sites, supercoiling characteristics were similar to poorly transcribed zones. We present a model that explains why RNAP blocks plectonemic supercoil movement in the transcribed track and suggests how gyrase and TopA control upstream and downstream transcription-driven supercoiling.

Fitness effects of replichore imbalance in Salmonella enterica

A fitness cost due to imbalanced replichores has been proposed to provoke chromosome rearrangements in Salmonella enterica serovars. To determine the impact of replichore imbalance on fitness, the relative fitness of isogenic Salmonella strains containing transposon-held duplications of various sizes and at various chromosomal locations was determined. Although duplication of certain genes influenced fitness, a replichore imbalance of up to 16° did not affect fitness.

Duplication frequency in a population of Salmonella enterica rapidly approaches steady state with or without recombination

Tandem duplications are among the most common mutation events. The high loss rate of duplication suggested that the frequency of duplications in a bacterial population (1/1000) might reflect a steady state dictated by relative rates of formation (k(F)) and loss (k(L)). This possibility was tested for three genetic loci. Without homologous recombination (RecA), duplication loss rate dropped essentially to zero, but formation rate decreased only slightly and a steady state was still reached rapidly. Under all conditions, steady state was reached faster than predicted by formation and loss rates alone. A major factor in determining steady state proved to be the fitness cost, which can exceed 40% for some genomic regions. Depending on the region tested, duplications reached 40-98% of the steady-state frequency within 30 generations-approximately the growth required for a single cell to produce a saturated overnight culture or form a large colony on solid medium (10(9) cells). Long-term bacterial populations are stably polymorphic for duplications of every region of their genome. These polymorphisms contribute to rapid genetic adaptation by providing frequent preexisting mutations that are beneficial whenever imposed selection favors increases in some gene activity. While the reported results were obtained with the bacterium Salmonella enterica, the genetic implications seem likely to be of broader biological relevance.

Adaptive expression of foreign genes in the clonal variants of bacteria: from proteomics to clinical application

Clonal variants of bacteria are able to colonize environmental niches and patients. The factors, that determine the interplay between the colonization of diverse habitats and adaptation, are acquired through horizontal gene transfer. Elucidation of mechanisms, which lead to the prevalence of dominant bacterial clones in patients and the environment, requires the knowledge of complex phenotypes. It was found in the genomes of most bacteria, that upon a conserved chromosomal backbone there were regions of plasticity achieved by insertions, deletions and rearrangements of genomic islands and islets as well as large chromosomal inversions. However, it had been shown that environmental and clinical isolates are indistinguishable in certain pathogenic and biodegradative properties. For example, clonal variants of Pseudomonas aeruginosa exhibit convergent phenotypes despite the presence of numerous DNA insertions in the genome. Apart from this feature, expression of a few genes from the acquired genetic material is important for niche-based adaptation of this organism. Protein expression patterns at the cellular and sub-cellular levels showed common virulence factors and novel drug targets among clonal variants of bacteria. This review will give a short overview on proteomics of different clonal variants of bacteria with respect to clinical applications.

Diverse phenotypes resulting from polyphosphate kinase gene (ppk1) inactivation in different strains of Helicobacter pylori

Connections among biochemical pathways should help buffer organisms against environmental stress and affect the pace and trajectory of genome evolution. To explore these ideas, we studied consequences of inactivating the gene for polyphosphate kinase 1 (ppk1) in strains of Helicobacter pylori, a genetically diverse gastric pathogen. The PPK1 enzyme catalyzes synthesis of inorganic polyphosphate (poly P), a reservoir of high-energy phosphate bonds with multiple roles. Prior analyses in less-fastidious microbes had implicated poly P in stress resistance, motility, and virulence. In our studies, ppk1 inactivation caused the expected near-complete absence of poly P (>250-fold decrease) but had phenotypic effects that differed markedly among unrelated strains: (i) poor initial growth on standard brain heart infusion agar (five of six strains tested); (ii) weakened colonization of mice (4 of 5 strains); (iii) reduced growth on Ham's F-12 agar, a nutritionally limiting medium (8 of 11 strains); (iv) heightened susceptibility to metronidazole (6 of 17 strains); and (v) decreased motility in soft agar (1 of 13 strains). Complementation tests confirmed that the lack of growth of one Deltappk1 strain on F-12 agar and the inability to colonize mice of another were each due to ppk1 inactivation. Thus, the importance of ppk1 to H. pylori differed among strains and the phenotypes monitored. We suggest that quantitative interactions, as seen here, are common among genes that affect metabolic pathways and that H. pylori's high genetic diversity makes it well suited for studies of such interactions, their underlying mechanisms, and their evolutionary consequences.

Diversity of genome structure in Salmonella enterica serovar Typhi populations

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

Concerted evolution in the repeats of an immunomodulating cell surface protein, SOWgp, of the human pathogenic fungi Coccidioides immitis and C. posadasii

Genome dynamics that allow pathogens to escape host immune responses are fundamental to our understanding of host-pathogen interactions. Here we present the first population-based study of the process of concerted evolution in the repetitive domain of a protein-coding gene. This gene, SOWgp, encodes the immunodominant protein in the parasitic phase of the human pathogenic fungi Coccidioides immitis and C. posadasii. We sequenced the entire gene from strains representing the geographic ranges of the two Coccidioides species. By using phylogenetic and genetic distance analyses we discovered that the repetitive part of SOWgp evolves by concerted evolution, predominantly by the mechanism of unequal crossing over. We implemented a mathematical model originally developed for multigene families to estimate the rate of homogenization and recombination of the repetitive array, and the results indicate that the pattern of concerted evolution is a result of homogenization of repeat units proceeding at a rate close to the nucleotide point mutation rate. The release of the SOWgp molecules by the pathogen during proliferation may mislead the host: we speculate that the pathogen benefits from concerted evolution of repeated domains in SOWgp by an enhanced ability to misdirect the host's immune system.

Genomic rearrangements at rrn operons in Salmonella

Most Salmonella serovars are general pathogens that infect a variety of hosts. These "generalist" serovars cause disease in many animals from reptiles to mammals. In contrast, a few serovars cause disease only in a specific host. Host-specific serovars can cause a systemic, often fatal disease in one species yet remain avirulent in other species. Host-specific Salmonella frequently have large genomic rearrangements due to recombination at the ribosomal RNA (rrn) operons while the generalists consistently have a conserved chromosomal arrangement. To determine whether this is the result of an intrinsic difference in recombination frequency or a consequence of lifestyle difference between generalist and host-specific Salmonella, we determined the frequency of rearrangements in vitro. Using lacZ genes as portable regions of homology for inversion analysis, we found that both generalist and host-specific serovars of Salmonella have similar tolerances to chromosomal rearrangements in vitro. Using PCR and genetic selection, we found that generalist and host-specific serovars also undergo rearrangements at rrn operons at similar frequencies in vitro. These observations indicate that the observed difference in genomic stability between generalist and host-specific serovars is a consequence of their distinct lifestyles, not intrinsic differences in recombination frequencies.

Atypical 16S rRNA gene copies in Ochrobactrum intermedium strains reveal a large genomic rearrangement by recombination between rrn copies

Ochrobactrum intermedium is an opportunistic human pathogen belonging to the alpha 2 subgroup of proteobacteria. The 16S rDNA sequences of nine O. intermedium isolates from a collection of clinical and environmental isolates exhibited a 46-bp insertion at position 187, which was present in only one sequence among the 82 complete or partial 16S rDNA sequences of Ochrobactrum spp. available in data banks. Reverse transcription-PCR experiments showed that the 46-bp insertion remained in the 16S rRNA. The inserted sequence folded into a stem-loop structure, which took place in and prolonged helix H184 of the 16S rRNA molecule. Helix H184 has been described as conserved in length among eubacteria, suggesting the idiosyncratic character of the 46-bp insertion. Pulsed-field gel electrophoresis experiments showed that seven of the clinical isolates carrying the 46-bp insertion belonged to the same clone. Insertion and rrn copy numbers were determined by hybridization and I-CeuI digestion. In the set of clonal isolates, the loss of two insertion copies revealed the deletion of a large genomic fragment of 150 kb, which included one rrn copy; deletion occurred during the in vivo evolution of the clone. Determination of the rrn skeleton suggested that the large genomic rearrangement occurred during events involving homologous recombination between rrn copies. The loss of insertion copies suggested a phenomenon of concerted evolution among heterogeneous rrn copies.

Rapid approach to determine rrn arrangement in Salmonella serovars

A PCR method was developed by which to rapidly and accurately determine the rrn arrangement of Salmonella enterica serovars. Primers were designed to the genomic regions flanking each of the seven rrn operons. PCR analysis using combinations of these primers will distinguish each of the possible arrangements of the rrn skeleton.

Genetic mapping by duplication segregation in Salmonella enterica

MudP and MudQ elements were used to induce duplications in Salmonella enterica by formation of a triple crossover between two transduced fragments and the host chromosome. The large size (36 kb) of MudP and MudQ is a favorable trait for duplication formation, probably because homology length is a limiting factor for the central crossover. Additional requirements are a multiplicity of infection of 2 or higher in the infecting phage suspensions (which reflects the need of two transduced fragments) and an exponentially growing recipient (which reflects the need of a chromosome replication fork). We describe a set of 11 strains of S. enterica, each carrying a chromosomal duplication with known endpoints. The collection covers all the Salmonella chromosome except the terminus. For mapping, a dominant marker (e.g., a transposon insertion in or near the locus to be mapped) is transduced into the 11-strain set. Several transductants from each cross are grown nonselectively, and haploid segregants are scored for the presence of the marker. If all the segregants contain the transduced marker, it maps outside the duplication interval. If the marker is found only in a fraction of the segregants, it maps within the duplicated region.

Complete genome sequence of Neisseria meningitidis serogroup B strain MC58

The 2,272,351-base pair genome of Neisseria meningitidis strain MC58 (serogroup B), a causative agent of meningitis and septicemia, contains 2158 predicted coding regions, 1158 (53.7%) of which were assigned a biological role. Three major islands of horizontal DNA transfer were identified; two of these contain genes encoding proteins involved in pathogenicity, and the third island contains coding sequences only for hypothetical proteins. Insights into the commensal and virulence behavior of N. meningitidis can be gleaned from the genome, in which sequences for structural proteins of the pilus are clustered and several coding regions unique to serogroup B capsular polysaccharide synthesis can be identified. Finally, N. meningitidis contains more genes that undergo phase variation than any pathogen studied to date, a mechanism that controls their expression and contributes to the evasion of the host immune system.

Resurgent Vibrio cholerae O139: rearrangement of cholera toxin genetic elements and amplification of rrn operon

The unprecedented genesis of a novel non-O1 Vibrio cholerae strain belonging to serogroup O139, which caused an epidemic in late 1992 in the Indian subcontinent, and its subsequent displacement by El Tor O1 vibrios after 18 months initiated a renewed investigation of the aspects of the organism that are related to pathogenesis. The reappearance of V. cholerae O139 with altered antibiotic sensitivity compared to O139 Bengal (O139B) in late 1996 has complicated the epidemiological scenario of V. cholerae and has necessitated an examination of possible rearrangements in the genome underlying such rapid changes in the phenotypic traits. With a view to investigating whether the phenotypic changes that have occurred are associated with alteration in the genome, the genome of the resurgent V. cholerae O139 (O139R) strains were examined. Pulsed-field gel electrophoresis analysis of NotI- and SfiI-digested genomic DNA of O139R isolates showed restriction fragment length polymorphism including in the cholera toxin (CTX) genetic element locus and with O139B isolates. Analyses of the organization of the CTX genetic elements in O139R strains showed that in contrast to two copies of the elements connected by two direct-repeat sequences (RS) in most of the genomes of O139B isolates, the genomes of all O139R strains examined, except strain AS192, have three such elements connected by a single RS. While the RS present in the upstream of the CTX genetic elements in the genome of O139R is of O139B origin, the RS connecting the cores of the elements has several new restriction sites and has lost the BglII site which is supposed to be conserved in all O1 strains and O139B. The endonuclease I-CeuI, which has sites only in the rrn operons in the genomes of all organisms examined so far, has 10 sites in the genomes of O139R strains, compared to 9 in the genomes of O139B strains. The recent isolates of V. cholerae O139 have thus gained one rrn operon. This variation in the number of rrn operons within a serogroup has not been reported for any other organism. The results presented in this report suggest that like the pathogenic El Tor O1 strains, the genomes of O139 strains are undergoing rapid alterations.

Homologous recombination between rrn operons rearranges the chromosome in host-specialized species of Salmonella

Partial digestion with I-CeuI, which digests bacterial DNA at the gene coding for the large subunit rRNA, established the rrn genomic skeleton (the distance in kb between rRNA operons) in 56 strains of Salmonella, from Salmonella Reference B (SARB) set. All had seven I-CeuI sites, indicating seven rrn operons. The order of I-CeuI fragments was ABCDEFG in S. typhimurium LT2 and in 31 other species, mostly host-generalists; in S. typhi, S. paratyphi C, S. gallinarum, and S. pullorum (host-specialized species), these fragments are rearranged, due to homologous recombination between the rrn operons. Rearrangements, such as inversions and translocations not involving the rrn operons, are rare. I-CeuI fragments of some species are larger than the norm, suggesting the insertion of unique blocks of DNA by lateral transfer from other species.

Chromosomal rearrangements in enteric bacteria

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

Conserved structure of IS200 elements in Salmonella

Sequence analysis of three IS200 elements (two from Salmonella typhimurium, one from Salmonella abortusovis) reveals a highly conserved structure, with a length of 707-708 bp and absence of terminal repeats. IS200 contains an open-reading-frame (ORF) which potentially encodes a peptide of 151 amino acids, with a putative ribosome-binding-site properly placed upstream of the ORF. A potential RNA stem-loop structure that might occlude the ribosome-binding-site of the ORF is also found. Another conserved trait is a potential RNA hairpin which resembles a Rho-independent transcription terminator, located near one end of IS200. The junctions between IS200 and host DNA sequences are A+T-rich. Upon insertion, IS200 duplicates 1-2 bp of host DNA sequences. The observation that IS200 elements characterized as 'hops' are roughly identical to those residing in the Salmonella genome suggests that IS200 transposition is unlikely to generate inactive copies. If such is the case and many or all IS200 elements are active, the extremely low frequency of IS200 transposition may reflect the normal behavior of the element.

Highly plastic chromosomal organization in Salmonella typhi

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

The chromosome of Salmonella paratyphi A is inverted by recombination between rrnH and rrnG

Salmonella paratyphi A, a human-adapted bacterial pathogen, causes paratyphoid enteric fever. We established the genome map of strain ATCC 9150 by the use of four endonucleases, XbaI, I-CeuI, AvrII (= BlnI), and SpeI, which generated 27, 7, 19, and 38 fragments, respectively; the sum of the fragments in each case indicates a genome size of ca. 4,600 kb. With phage P22, we transduced Tn10 insertions in known genes from Salmonella typhimurium LT2 to S. paratyphi A ATCC 9150 and located these insertions on the S. paratyphi A chromosome through the XbaI and AvrII sites in Tn10 and through the increased size of the SpeI fragment bearing a Tn10. Compared with the maps of other Salmonella species, the S. paratyphi A genomic map showed two major differences: (i) an insertion of about 100 kb of DNA between rrnH/G and proB and (ii) an inversion of half the genome between rrnH and rrnG, postulated to be due to homologous recombination between the rrn genes. We propose that during the evolution of S. paratyphi A, the first rearrangement event was the 100-kb insertion, which disrupted the chromosomal balance between oriC and the termination of replication, forcing the rrnH/G inversion to restore the balance. The insertion and the inversion are both present in all 10 independent wild-type S. paratyphi A strains tested.

Suppression of the pleiotropic effects of HisH and HisF overproduction identifies four novel loci on the Salmonella typhimurium chromosome: osmH, sfiW, sfiX, and sfiY

Insertion mutations that suppress some or all the pleiotropic effects of HisH and HisF overproduction were obtained by using transposons Tn10dTet and Tn10dCam. All suppressor mutations proved to be recessive, indicating that their effects were caused by loss of function; thus, the suppressors identify genes that are necessary to trigger the pleiotropic response when HisH and HisF are overproduced. Genetic mapping of the suppressor mutations identifies four novel loci on the Salmonella typhimurium genetic map. Mutations in osmH (min 49) behave as general suppressors that abolish all manifestations of the pleiotropic response. Mutations in sfiY (min 83) suppress cell division inhibition and thermosensitivity but not osmosensitivity. Mutations that suppress only cell division inhibition define another locus, sfiX (min 44). A fourth novel locus, sfiW (min 19), is also involved in cell division inhibition. The phenotype of sfiW mutations is in turn pleiotropic: they suppress cell division inhibition, make S. typhimurium unable to grow in minimal media, and cause slow growth and abnormal colony and cell shape. The inability of sfiW mutants to grow in minimal medium cannot be relieved by any known nutritional requirement or by the use of carbon sources other than glucose. The hierarchy of suppressor phenotypes and the existence of epistatic effects among suppressor mutations suggest a pathway-like model for the Hisc pleiotropic response.

Rearrangements in the genome of the bacterium Salmonella typhi

We have determined the genomic map of the bacterium Salmonella typhi Ty2, the causal organism of typhoid fever, by using pulsed-field gel electrophoresis. Digestion of the Ty2 genome with endonucleases Xba I, Bln I, and Ceu I yielded 33, 26, and 7 fragments, respectively, that were placed in order on a circular chromosome of 4780 kb. Transposon Tn10 was inserted in specific genes of Salmonella typhimurium and transduced into S. typhi, and thus, the positions of 37 S. typhi genes were located through the Xba I and Bln I sites of the Tn10. Gene order on chromosomes of Escherichia coli K-12 and S. typhimurium LT2 is remarkably conserved; however, the gene order in S. typhi Ty2 is different, suggesting it has undergone major genomic rearrangements during its evolution. These rearrangements include inversions and transpositions in the 7 DNA fragments between the seven rrn operons for rRNA (postulated to be due to homologous recombination in these rrn genes), another inversion that covers the replication terminus region (resembling inversions found in other enteric bacteria), and at least three insertions, one as large as 118 kb. Partial digestion of genomic DNA with the intron-encoded endonuclease I-Ceu I, which cuts only in rrn genes, shows chromosomal rearrangements, apparently due to homologous recombination in the rrn genes, that were detected in all wild-type strains of S. typhi tested. These rearrangements may have been selected to compensate for the insertions that otherwise would have altered the locations of genes with respect to the origin and terminus of replication. These observations are relevant to our view of the evolution of the bacterial genome and may be significant in the virulence of S. typhi.

Mutations in rpoA affect expression of anaerobically regulated genes in Salmonella typhimurium

oxrB8, a mutation that diminishes the anaerobic induction of pepT and other anaerobically regulated, oxrA (fnr)-dependent Salmonella typhimurium genes, is an allele of rpoA, the gene for the alpha subunit of RNA polymerase. Four additional rpoA mutations that affect anaerobic pepT expression have been isolated after localized mutagenesis of the rpoA region. All but one of these rpoA mutations appear to have relatively specific effects on genes that require the OxrA (FNR) protein, a positive transcriptional regulator of a family of anaerobically expressed genes. All of these mutations lead to amino acid substitutions in the C-terminal region of the alpha subunit. These results taken with a number of previous observations suggest a role for the alpha subunit in the interaction between RNA polymerase and positive transcriptional regulatory proteins. They also suggest that the C-terminal region of alpha is important for these interactions.

beta-Alanine auxotrophy associated with dfp, a locus affecting DNA synthesis in Escherichia coli

Strains containing the conditional-lethal dfp-707 mutation, which have a defect in DNA synthesis at 42 degrees C, were found to require either pantothenate or its precursor, beta-alanine, for growth at 30 degrees C. The auxotrophy and conditional lethality were corevertible. Through localized mutagenesis of the dfp-pyrE region of Escherichia coli, another mutation, dfp-1, was obtained. It conferred the auxotrophy but not the conditional lethality of dfp-707. Complementation analysis, performed with a set of plasmid-borne deletion and insertion mutations, revealed a correspondence between the complementation of each mutant phenotype and the production of the dfp gene product, previously identified as a 45-kilodalton flavoprotein. The dfp mutants had a normal level of aspartate-1-decarboxylase, which is the only enzyme known to produce beta-alanine in E. coli and which is specified by the distant panD gene. A prototrophic pseudorevertant of a dfp-1 strain was found to have retained the dfp mutation, to be genetically unstable, and to have an elevated level of aspartate-1-decarboxylase, suggesting that it had acquired a duplication of panD. It is not known what steps in pantothenate or DNA metabolism are affected by the mutant dfp product or how its flavin moiety may be involved.

N-terminal methionine-specific peptidase in Salmonella typhimurium

Crude extracts of a multiply peptidase-deficient strain of Salmonella typhimurium contain an aminopeptidase that specifically removes N-terminal methionine from peptides. This activity shows pronounced specificity for the peptide's second amino acid. Methionine is removed from peptides with alanine, threonine, or glycine in this position but not when the second amino acid is leucine or methionine. The activity is stimulated by Co2+ and is inhibited by EDTA. Mutations that lead to overproduction (up to 30-fold) of the activity have been obtained by selecting for growth on Met-Gly-Gly as a methionine source. These mutations map at approximately 3 map units, phage P22 cotransducible with leu. The overproducer mutations are dominant to wild type, and duplication of the wild-type allele of the locus leads to a gene dosage effect on peptidase levels. This suggests that the locus of the overproducer mutations may be the structural gene for the peptidase. NaDodSO4/PAGE shows an increased level of a single protein (34 kDa) in the overproducer mutant. This protein is highly enriched in a purified preparation of the peptidase. The specificity of this enzyme suggests that it is involved in the cleavage of methionine from newly synthesized peptide chains. This activity can specifically remove methionine from the N terminus of a completed protein. Treatment of purified, unprocessed (N-terminal methionine) interleukin 1 beta with the purified peptidase results in removal of N-terminal methionine with no additional alterations. N-terminal processing of at least this protein can occur after translation is complete. We propose to call this enzyme peptidase M (methionine-specific aminopeptidase).

Mapping and spacer identification of rRNA operons of Salmonella typhimurium

The rRNA operons of Salmonella typhimurium have been characterized with respect to their map position, orientation, and type of tRNA spacer. One of the seven rrn operons was found to be linked to pheA and another was found to be linked to aroE. This information, together with published information about the other five rrn operons, shows that S. typhimurium and Escherichia coli are essentially identical in terms of the number, the map position, and the orientation of all seven operons. S. typhimurium and E. coli were also similar in that four of the rrn spacer regions code for tRNAGlu2 and three code for tRNAAla1B. However, the two species differed in that rrnD coded for tRNAGlu2 and rrnB coded for tRNAAla1B in S. typhimurium. This is the opposite of the arrangement in E. coli. We have tabulated the coordinates of the BamHI and PstI sites flanking six of the S. typhimurium rrn genes and present revisions for the coordinates of some of the E. coli sites.

Aspartate-specific peptidases in Salmonella typhimurium: mutants deficient in peptidase E

The only dipeptide found to serve as a leucine source for a Salmonella strain lacking peptidases N, A, B, D, P, and Q was alpha-L-aspartyl-L-leucine. A peptidase (peptidase E) that specifically hydrolyzes Asp-X peptides was identified and partially purified from cell extracts. The enzyme (molecular weight, 35,000) is inactive toward dipeptides with N-terminal asparagine or glutamic acid. Mutants (pepE) lacking this enzyme were isolated by screening extracts for loss of the activity. Genetic mapping placed the pepE locus at 91.5 map units and established the gene order metA pepE zja-861::Tn5 malB. Duplications of the pepE locus showed a gene dosage effect on levels of peptidase E, suggesting that pepE is the structural gene for this enzyme. Mutations in pepE resulted in the loss of the ability to grow on Asp-Pro as a proline source but did not affect utilization of other dipeptides with N-terminal aspartic acid. Loss of peptidase E did not cause a detectable impairment in protein degradation. Two other peptidases present in cell extracts of mutants lacking peptidases N, A, B, D, P, Q, and E also hydrolyze many Asp-X dipeptides.

recA-dependent recombination between rRNA operons generates type II F' plasmids

The formation of type II F' ilv cya metE plasmids from the Salmonella typhimurium Hfr strain SA722 occurs by general recombination between repeated rrn.

In situ transposon replacement and isolation of a spontaneous tandem genetic duplication

Using a specialized transducing P1 phage carrying an insertion of Tn5-132, an insertion of Tn5-wt in the chromosome of Myxococcus xanthus, which codes for resistance to kanamycin, can be replaced with one of Tn5-132, which codes for resistance to tetracycline. That Tn5-132 in the daughter is inserted at the same location in the chromosome as Tn5-wt was in the parent was shown by a variety of physical and genetic tests. Southern blot hybridizations of restriction digests of daughter and parent DNAs probed for sequences homologous to Tn5 show that the physical location is the same. When KmR was transduced from the parent to the TcR daughter by the generalized transducing myxophage Mx4 or Mx8, all the transductants were TcS. Likewise, when the daughter was used as donor, TcR transductants of its KmR parent were KmS. Flanking markers that were linked to KmR in the parent were linked to TcR in the daughter. Spontaneous tandem genetic duplications of portions of bacterial chromosomes can be trapped by transducing a selectable marker from a donor to a recipient that has a different selectable marker at the same genetic location and selecting transductants with both markers. Using Tc-replacement, this technique can be applied to any region of the chromosome. We used it to isolate a spontaneous tandem duplication of part of the M. xanthus chromosome. The duplication was characterized by Southern blot hybridizations probed for Tn5-homologous DNA. It was also shown to be unstable by quantitation of loss of drug resistance. Transduction of the novel joint led to reconstruction of the duplication in the recipient strain. All these tests gave results consistent with the proposed structure. The methods described here are applicable to any bacterium into which transposons can be introduced, and for which some means of genetic exchange is available.

Salmonella typhimurium mutations affecting utilization of L-leucine beta-naphthylamide

L-leucine-beta-naphthylamide (LNA) will support growth of a leucine auxotroph of Salmonella typhimurium. Utilization of this compound depends on the presence in the cells of active peptidase N. Selection for improved growth on a suboptimal concentration of LNA yields mutants some of which contain elevated levels of peptidase N. The properties of these strains indicate that they carry tandem genetic duplication of the pepN locus: they show rec-dependent genetic instability; they contain an approximately doubled level of the pepN gene product; neighboring chromosomal loci are also duplicated; and, the mutants occur with a greatly diminished frequency in rec- strains. When selection for improved growth on LNA is applied to a rec- strain, the mutants obtained do not contain duplications. These strains appear to contain lesions in the pepN gene that lead to the production of a peptidase N with altered substrate specificity.

Chromosomal locations of the genes for rRNA in Escherichia coli K-12

Chromosomal locations of the seven rRNA operons in Escherichia coli K-12 were studied by digesting DNA from various merodiploid strains with SalI restriction enzyme followed by Southern gel analysis with 32P-labeled 23S rRNA as a probe. The seven unique SalI DNA fragments revealed in the autoradiograms were first correlated to the seven rRNA operons previously isolated as hybrid plasmids or transducing phages. The chromosomal locations of six (rrnA, B, C, D, E, and G) of the seven isolated operons were confirmed by increased gene dosage demonstrated in autoradiograms after Southern gel analysis of DNA from relevant merodiploid strains. The gene dosage analysis showed that the location of the remaining operon (now called rrnH) is between metD and proA. No evidence was obtained for the presence of rrnF, which was previously reported to map between aroB and malA. The chromosomal location of rrnH was confirmed by P1 transduction in the following way: a DNA fragment adjacent to rrnH was cloned into pBR322; the resulting hybrid plasmid was integrated at the homologous region of the chromosome of a polA mutant; and the ampicillin resistance marker originally carried by pBR322 was then used for mapping of the nearby rrnH by P1 transduction. A close linkage of rrnH to metD (about 60% cotransduction) was observed, and the data were consistent with the order metD-rrnH-proA. Thus, mapping of all seven rRNA operons has been completed. The present study has also determined the orientation of rrnG and rrnH and demonstrated that the direction of transcription of all the rRNA operons is identical to that of DNA replication.

Transposition of a chromosomal segment bounded by redundant rRNA genes into other rRNA genes in Escherichia coli

We have constructed several mutants of Escherichia coli which have the chromosomal segment between the directly repeated rrnB and rrnE genes deleted from the normal position and transposed into another one of the seven redundant rRNA genes. We have examples where the transposition has been into rrnC, rrnD, rrnG, and rrnH. Included in the evidence for each of these transpositions was the finding that each transposition specifically affected a different one of the seven BamHI-PstI restriction nuclease fragments known to correspond to the seven rrn genes. The transposition mutants were generally healthy, but sensitive mixed-growth experiments revealed that most of them grew somewhat more slowly than the parental control in rich medium. The maximal detrimental effect was a 4 to 5% reduction in growth rate when the transposition of the rrnB-rrnE segment was into rrnG. We have found that a rrnF gene, reported by others to be linked to malA, does not exist in our standard strain, a derivative of Cavalli Hfr. Instead of rrnF, we identified a new rrn gene, rrnH, which mapped near min 5.

Inversions between ribosomal RNA genes of Escherichia coli

It might be anticipated that the presence of redundant but oppositely oriented sequences in a chromosome could allow inversion of the intervening material through homologous recombination. For example, the ribosomal RNA gene rrnD of Escherichia coli has the opposite orientation fro rrnB and rrnE and is separated from these genes by roughly 20% of the chromosome. Starting with a derivative of Cavalli Hfr, we have constructed mutants that have an inversion of the segment between rrnD and either rrnB or rrnE. These mutants are generally quite viable but do exhibit a slight reduction in growth rate relative to the parental strain. A major line of laboratory E. coli, W3110 and its derivatives, also has an inversion between rrnD and rrnE, probably created directly by a recombinational event between these highly homologous genes.

Spontaneous tandem genetic duplications in Salmonella typhimurium arise by unequal recombination between rRNA (rrn) cistrons

A method is described to detect and measure the frequency of spontaneous tandem genetic duplications located throughout the Salmonella genome. The method is based on the ability of duplication-containing strains to inherit two selectable alleles of a single gene during generalized transductional crosses. One allele of the gene carries an insertion of the translocatable tetracycline-resistance element Tn10; the other allele is a wild-type copy of that gene. Using this technique, we have measured the frequency of tandem duplications at 38 chromosomal sites and the amount of material included in 199 independent duplications. These results suggest that, in one region of the chromosome, tandem duplications are particularly frequent events. Such duplications have end points within rRNA (rrn) cistrons and probably arise by unequal cross-over between these dispersed repeated sequences. Spontaneously duplications of this type are harbored by as much as 3% of the bacterial population. Preliminary evidence suggests that such duplications may play a significant regulatory role under conditions of rapid growth. Our analysis has suggested the position on the genome of an additional rRNA cistron.

Regulation of Tn5 transposition in Salmonella typhimurium

The drug-resistance element Tn5 transposes with high frequency immediately after entry into a cell. Establishment of Tn5 within a cell results in a decrease in this transposition frequency. This phenomenon resembles "zygotic induction" of repressible operons and prophages. Evidence is presented that Tn5 transposition is under negative control by a factor encoded within the element itself. Established Tn5 elements (that contain point mutations inactivating the resistance gene) are able to inhibit transposition of an incoming Tn5 element by a factor of 12- to 70-fold. Several deletion derivatives of Tn5 lack the ability to inhibit transposition.

UGA suppressor that maps within a cluster of ribosomal protein genes

A suppressor of UGA mutations (supU) maps near or within a cluster of ribosomal protein genes at 72 min on the Salmonella typhimurium genetic map. The suppressor is relatively inefficient, and its activity is abolished by rpsL (formerly strA) mutations. The suppressor is dominant to a wild-type supU allele. The map position of this suppressor suggests that it may owe its activity to an alteration of ribosome structure.

Deletion of a ribosomal ribonucleic acid operon in Escherichia coli

One (rrnE) of the seven operons which codes for ribosomal ribonucleic acid in Escherichia coli was deleted. No significant change in phenotype was observed even under maximum laboratory growth conditions.

Involvement of ribosomal ribonucleic acid operons in Salmonella typhimurium chromosomal rearrangements

As part of our efforts to understand factors influencing chromosomal organization and rearrangements, we studied a family of Salmonella typhimurium tandum duplication mutants. We found that the duplications were originally generated by unequal recombination between pairs of similarly oriented ribosomal ribonucleic acid operons (rrn). This demonstration involved the physical isolation of the duplicated material as circular deoxyribonucleic acid excised from the duplication. The four rrn operons involved embraced the ilv pur D segment of the chromosome and occurred at positions closely analogous to those previously observed for Escherichia coli. The interval between rrnC and rrnA of S. typhimurium was similar in size to that of E. coli (43 versus 39 kilobases), as was the interval between rrnB and rrnE (94 versus 91 kilobases). The rrnA-to-rrnB interval of S. typhimurium, however, was 155 kilobases, substantially greater than the 126 kilobases observed for E. coli.

Rapid duplication and loss of genes coding for the alpha chains of hemoglobin

Rapid cycles of gene duplication and loss appear to have been going on in the region coding for the alpha chain of adult hemoglobin. This is inferred from restriction endonuclease analysis of the alpha gene region in five species of apes, whose common ancestor lived about 10 million years ago. Because all five species resemble humans in having duplicate alpha genes, the duplicate state of this region is probably at least as old as the common ancestor of all these species. However, the alpha polypeptides within these species are about 10 times more alike than is expected for 10 million years of divergent evolution. Thus, the alpha polypeptides within each species have been evolving in concert. Changes in gene number have also taken place in the apes. Whereas the predominant number of alpha genes per chromosome is two for most species, it is three for chimpanzees. Concerted evolution appears also to have occurred, but far more slowly, in the region coding for the adult beta-like chains of hemoglobin. Consideration of the structural differences between the two regions leads to the hypothesis that the lengths of the noncoding regions are important determinants of the rates at which genes are gained and lost by intergenic recombination.