The Evolution of Insertion Sequences within Enteric Bacteria

To identify mechanisms that influence the evolution of bacterial transposons, DNA sequence variation was evaluated among homologs of insertion sequences IS1, IS3 and IS30 from natural strains of Escherichia coli and related enteric bacteria. The nucleotide sequences within each class of IS were highly conserved among E. coli strains, over 99.7% similar to a consensus sequence. When compared to the range of nucleotide divergence among chromosomal genes, these data indicate high turnover and rapid movement of the transposons among clonal lineages of E. coli. In addition, length polymorphism among IS appears to be far less frequent than in eukaryotic transposons, indicating that nonfunctional elements comprise a smaller fraction of bacterial transposon populations than found in eukaryotes. IS present in other species of enteric bacteria are substantially divergent from E. coli elements, indicating that IS are mobilized among bacterial species at a reduced rate. However, homologs of IS1 and IS3 from diverse species provide evidence that recombination events and horizontal transfer of IS among species have both played major roles in the evolution of these elements. IS3 elements from E. coli and Shigella show multiple, nested, intragenic recombinations with a distantly related transposon, and IS1 homologs from diverse taxa reveal a mosaic structure indicative of multiple recombination and horizontal transfer events.     

Escherichia coli as a model system to study DNA repair genes of eukaryotic organisms

The bacteria Escherichia coli has been widely employed in studies of eukaryotic DNA repair genes. Several eukaryotic genes have been cloned by functional complementation of mutant lineages of E. coli. We examined the similarities and differences among bacterial and eukaryotic DNA repair systems. Based on these data, we examined tools used for gene cloning and functional studies of DNA repair in eukaryotes, using this bacterial system as a model.     

Acquisition of the rfb-gnd cluster in evolution of Escherichia coli O55 and O157

The rfb region specifies the structure of lipopolysaccharide side chains that comprise the diverse gram-negative bacterial somatic (O) antigens. The rfb locus is adjacent to gnd, which is a polymorphic gene encoding 6-phosphogluconate dehydrogenase. To determine if rfb and gnd cotransfer, we sequenced gnd in five O55 and 13 O157 strains of Escherichia coli. E. coli O157:H7 has a gnd allele (allele A) that is only 82% identical to the gnd allele (allele D) of closely related E. coli O55:H7. In contrast, gnd alleles of E. coli O55 in distant lineages are >99.9% identical to gnd allele D. Though gnd alleles B and C in E. coli O157 that are distantly related to E. coli O157:H7 are more similar to allele A than to allele D, there are nucleotide differences at 4 to 6% of their sites. Alleles B and C can be found in E. coli O157 in different lineages, but we have found allele A only in E. coli O157 belonging to the DEC5 lineage. DNA 3' to the O55 gnd allele in diverse E. coli lineages has sequences homologous to tnpA of the Salmonella enterica serovar Typhimurium IS200 element, E. coli Rhs elements (including an H-rpt gene), and portions of the O111 and O157 rfb regions. We conclude that rfb and gnd cotransferred into E. coli O55 and O157 in widely separated lineages and that recombination was responsible for recent antigenic shifts in the emergence of pathogenic E. coli O55 and O157.     

Genome rearrangement distances and gene order phylogeny in gamma-Proteobacteria

Genome rearrangements have been studied in 30 gamma-proteobacterial complete genomes by comparing the order of a reduced set of genes on the chromosome. This set included those genes fulfilling several characteristics, the main ones being that an ortholog was present in every genome and that none of them had been acquired by horizontal gene transfer. Genome rearrangement distances were estimated based on either the number of breakpoints or the minimal number of inversions separating two genomes. Breakpoint and inversion distances were highly correlated, indicating that inversions were the main type of rearrangement event in gamma-Proteobacteria. In general, the progressive increase in sequence-based distances between genome pairs was associated with the increase in their rearrangement-based distances but with several groups of distances not following this pattern. Compared with free-living enteric bacteria, the lineages of Pasteurellaceae were evolving, on average, to relatively higher rates of between 2.02 and 1.64, while the endosymbiotic bacterial lineages of Buchnera aphidicola and Wigglesworthia glossinidia were evolving at moderately higher rates of 1.38 and 1.35, respectively. Because we know that the rearrangement rate in the Bu. aphidicola lineage was close to zero during the last 100-150 Myr of evolution, we deduced that a much higher rate took place in the first period of lineage evolution after the divergence of the Escherichia coli lineage. On the other hand, the lineage of the endosymbiont Blochmannia floridanus did present an almost identical rate to free-living enteric bacteria, indicating that the increase in the genome rearrangement rate is not a general change associated with bacterial endosymbiosis. Phylogenetic reconstruction based on rearrangement distances showed a different topology from the one inferred by sequence information. This topology broke the proposed monophyly of the three endosymbiotic lineages and placed Bl. floridanus as a closer relative to E. coli than Yersinia pestis. These results indicate that the phylogeny of these insect endosymbionts is still an open question that will require the development of specific phylogenetic methods to confirm whether the sisterhood of the three endosymbiotic lineages is real or a consequence of a long-branch attraction phenomenon.     

Heterogeneity of genome sizes among natural isolates of Escherichia coli.

Comparisons of the genetic maps of Escherichia coli K-12 and Salmonella typhimurium LT2 suggest that the size and organization of bacterial chromosomes are highly conserved. Employing pulsed-field gel electrophoresis, we have estimated the extent of variation in genome size among 14 natural isolates of E. coli. The BlnI and NotI restriction fragment patterns were highly variable among isolates, and genome sizes ranged from 4,660 to 5,300 kb, which is several hundred kilobases larger than the variation detected between enteric species. Genome size differences increase with the evolutionary genetic distance between lineages of E. coli, and there are differences in genome size among the major subgroups of E. coli. In general, the genomes of natural isolates are larger than those of laboratory strains, largely because of the fact that laboratory strains were derived from the subgroup of E. coli with the smallest genomes.     

Rapid identification of Escherichia coli safety and laboratory strain lineages based on Multiplex-PCR

Escherichia coli K-12, B, C and W strains are the most frequently used bacterial safety and laboratory strains. Lineage-specific DNA fragments were detected by microplate subtractive hybridization and utilized to create a fast differentiation method using a single PCR reaction to differentiate clearly the four lineages and separate them from pathogenic variants. The method has been evaluated on a comprehensive selection of widely used laboratory strains and a variety of pathogenic E. coli representatives. In addition, in silico analysis on all available E. coli genomes and the genomes of the close relatives Shigella and Salmonella confirmed the reliability of the proposed method. A fast identification and differentiation of E. coli safety strains by Multiplex-PCR is a useful tool for researchers and companies to check and monitor their reference stocks.     

Mhc-DQB repertoire variation in hominoid and Old World primate species.

Comparison of 87 distinct Mhc-DQB sequences, obtained from 13 primate species, demonstrates that five out of eight trans-species Mhc-DQB allele lineages are at least 30 million years old and predate divergence of hominoid and Old World primate species. One lineage may be much older because its members are not only traced back in higher primates, but also are present in a New World primate species. Comparing Mhc-DQB repertoire variation in distinct species, allows one to pinpoint when certain polymorphisms were lost or gained in primate evolution. Heterogeneity observed among members of trans-species Mhc-DQB allele lineages can be explained in major part by point mutations, whereas intraexonic crossing-over is a potent mechanism in generating new allele lineages. The stability of Mhc-DQB polymorphisms is influenced by selective forces because distinct allele lineages appear to have accumulated nucleotide substitutions and amino acid replacements at different rates.     

Chromosomal Homology and Molecular Organization of Muller''s Elements D and E in the Drosophila Repleta Species Group

Thirty-three DNA clones containing protein-coding genes have been used for in situ hybridization to the polytene chromosomes of two Drosophila repleta group species, D. repleta and D. buzzatii. Twenty-six clones gave positive results allowing the precise localization of 26 genes and the tentative identification of another nine. The results were fully consistent with the currently accepted chromosomal homologies and in no case was evidence for reciprocal translocations or pericentric inversions found. Most of the genes mapped to chromosomes 2 and 4 that are homologous, respectively, to chromosome arms 3R and 3L of D. melanogaster (Muller's elements E and D). The comparison of the molecular organization of these two elements between D. melanogaster and D. repleta (two species that belong to different subgenera and diverged some 62 million years ago) showed an extensive reorganization via paracentric inversions. Using a maximum likelihood procedure, we estimated that 130 paracentric inversions have become fixed in element E after the divergence of the two lineages. Therefore, the evolution rate for element E is approximately one inversion per million years. This value is comparable to previous estimates of the rate of evolution of chromosome X and yields an estimate of 4.5 inversions per million years for the whole Drosophila genome.     

Site-directed mutagenesis of the active site loop of the rhodanese-like domain of the human molybdopterin synthase sulfurase MOCS3. Major differences in substrate specificity between eukaryotic and bacterial homologs

Sequence alignments of human molybdopterin synthase sulfurase, MOCS3, showed that the N-terminal domain is homologous to Escherichia coli MoeB, whereas the C-terminal domain is homologous to rhodanese-like proteins. Previous studies showed that the activity of the separately purified rhodanese-like domain of MOCS3 displayed 1000-fold lower activity in comparison to bovine rhodanese with thiosulfate as sulfur source. When the six amino acid active site loop of MOCS3 rhodanese-like domain was exchanged with the loop found in bovine rhodanese, thiosulfate:cyanide sulfurtransferase activity was increased 165-fold. Site-directed mutagenesis of each individual residue of the active site loop of the MOCS3 rhodanese-like domain showed that the charge of the last amino acid determines thiosulfate sulfurtransferase activity. Replacing Asp417 by threonine resulted in 90-fold increased activity, whereas replacing it by arginine increased the activity 470-fold. Using a fully defined in vitro system containing precursor Z, MOCS2A, E. coli MoaE, E. coli MoeB, Mg-ATP, MOCS3 rhodanese-like domain, and thiosulfate, it was shown that sulfur transfer to MOCS2A was also affected by the alterations, but not as drastically. Our studies revealed that in humans and most eukaryotes thiosulfate is not the physiologic sulfur donor for MOCS3, whereas in bacterial homologs, which have an arginine at the last position of the active site loop, thiosulfate can be used as a sulfur source for molybdenum cofactor biosynthesis. The phylogenetic analysis of MoeB homologs showed that eukaryotic homologs are of bacterial origin. Furthermore, it could be shown that an MoeB homolog named MoeZ, where the dual CXXC zinc-binding motif of the MoeB domain is not present, arose independently several times during evolution.     

Solution structure of the bacterial chemotaxis adaptor protein CheW from Escherichia coli

The bacterial chemotaxis adaptor protein CheW physically links the chemoreceptors (MCPs) and the histidine kinase CheA. Extensive investigations using bacterium Escherichia coli have established the central role of CheW in the MCP-modulated activation of CheA. Here we report the solution structure of CheW from E. coli determined by NMR spectroscopy. The results show that E. coli CheW shares an overall fold with previously reported structure of CheW from Thermotoga maritima, whereas local conformational deviations are observed. In particular, the C-terminal alpha-helix is considerably longer in E. coli CheW and appears to shrink the active binding pocket with CheA. Our study provides the structural basis for further investigations in E. coli chemotaxis.     

Life history influences rates of climatic niche evolution in flowering plants

Across angiosperms, variable rates of molecular substitution are linked with life-history attributes associated with woody and herbaceous growth forms. As the number of generations per unit time is correlated with molecular substitution rates, it is expected that rates of phenotypic evolution would also be influenced by differences in generation times. Here, we make the first broad-scale comparison of growth-form-dependent rates of niche evolution. We examined the climatic niches of species on large time-calibrated phylogenies of five angiosperm clades and found that woody lineages have accumulated fewer changes per million years in climatic niche space than related herbaceous lineages. Also, climate space explored by woody lineages is consistently smaller than sister lineages composed mainly of herbaceous taxa. This pattern is probably linked to differences in the rate of climatic niche evolution. These results have implications for niche conservatism; in particular, the role of niche conservatism in the distribution of plant biodiversity. The consistent differences in the rate of climatic niche evolution also emphasize the need to incorporate models of phenotypic evolution that allow for rate heterogeneity when examining large datasets.     

Temporal dynamics of bacterial aging and rejuvenation

Single-celled organisms dividing by binary fission were thought not to age [1-4]. A 2005 study by Stewart et al. [5] reversed the dogma by demonstrating that Escherichia coli were susceptible to aging. A follow-up study by Wang et al. [6] countered those results by demonstrating that E. coli cells trapped in microfluidic devices are able to sustain robust growth without aging. The present study reanalyzed these conflicting data by applying a population genetic model for aging in bacteria [7]. Our reanalysis showed that in E. coli, as predicted by the model, (1) aging and rejuvenation occurred simultaneously in a population; (2) lineages receiving sequentially the maternal old pole converged to a stable attractor state; (3) lineages receiving sequentially the maternal new pole converged to an equivalent but separate attractor state; (4) cells at the old pole attractor had a longer doubling time than ones at the new pole attractor; and (5) the robust growth state identified by Wang et al. corresponds to our predicted attractor for lineages harboring the maternal old pole. Thus, the previous data, rather than opposing each other, together provide strong evidence for bacterial aging.     

COMPARATIVE ANALYSIS OF MORPHOLOGICAL DIVERSITY: DOES DISPARITY ACCUMULATE AT THE SAME RATE IN TWO LINEAGES OF CENTRARCHID FISHES?

Evolutionary lineages differ with regard to the variety of forms they exhibit. We investigated whether comparisons of morphological diversity can be used to identify differences in ecological diversity in two sister clades of centrarchid fishes. Species in the Lepomis clade (sunfishes) feed on a wider range of prey items than species in the Micropterus clade (black basses). We quantified disparity in morphology of the feeding apparatus as within-clade variance on principal components and found that Lepomis exhibits 4.4 and 7.4 times more variance than Micropterus on the first two principal components. However, lineages are expected to diversify morphologically and ecologically given enough time, and this pattern could have arisen due to differences in the amount of time each clade has had to accumulate variance. Despite being sister groups, the age of the most recent common ancestor of Lepomis is approximately 14.6 million years ago and its lineages have a total length of 86.4 million years while the age of the most recent common ancestor of Micropterus is only about 8.4 million years ago, and it has a total branch length of 42.9 million years. We used the Brownian motion model of character evolution to test the hypothesis that time of independent evolution of each clade's lineages accounts for differences in morphological disparity and determined that the rates of evolution of the first two principal components are 4.4 and 7.7 times greater in Lepomis. Thus, time and phylogeny do not account for the differences in morphological disparity observed in Lepomis and Micropterus, and other diversity-promoting mechanisms should be investigated.     

Evidence for genetic drift in endosymbionts (Buchnera): analyses of protein-coding genes.

Buchnera, the bacterial endosymbionts of aphids, undergo severe population bottlenecks during maternal transmission through their hosts. Previous studies suggest an increased effect of drift within these strictly asexual, small populations, resulting in an increased fixation of slightly deleterious mutations. This study further explores sequence evolution in Buchnera using three approaches. First, patterns of codon usage were compared across several homologous Escherichia coli and Buchnera loci, in order to test the prediction that selection for the use of optimal codons is less effective in small populations. A chi 2-based measure of codon bias was developed to adjust for the overall A + T richness of silent positions in the endosymbionts. In contrast to E. coli homologues, adaptive codon bias across Buchnera loci is markedly low, and patterns of codon usage lack a strong relationship with gene expression level. These data suggest that codon usage in Buchnera has been shaped largely by mutational pressure and drift rather than by selection for translational efficiency. One exception to the overall lack of bias is groEL, which is known to be constitutively overexpressed in Buchnera and other endosymbionts. Second, relative-rate tests show elevated rates of sequence evolution of numerous protein-coding loci across Buchnera, compared to E. coli. Finally, consistently higher ratios of nonsynonymous to synonymous substitutions in Buchnera loci relative to the enteric bacteria strongly suggest the accumulation of nonsynonymous substitutions in endosymbiont lineages. Combined, these results suggest a decreased effectiveness of purifying selection in purging endosymbiont populations of slightly deleterious mutations, particularly those affecting codon usage and amino acid identity.     

Low-temperature recovery strategies for the isolation of bacteria from ancient permafrost sediments

Permafrost represents a unique ecosystem that has allowed the prolonged survival of certain bacterial lineages at subzero temperatures. To better understand the permafrost microbial community, it is important to identify isolation protocols that optimize the recovery of genetically diverse bacterial lineages. We have investigated the impact of different low-temperature isolation protocols on recovery of aerobic bacteria from northeast Siberian permafrost of variable geologic origin and frozen for 5000 to 3 million years. Low-nutrient media enhanced the quantitative recovery of bacteria, whereas the isolation of diverse morphotypes was maximized on rich media. Cold enrichments done directly in natural, undisturbed permafrost led not only to recovery of increased numbers of bacteria but also to isolation of genotypes not recovered by means of liquid low-temperature enrichments. On the other hand, direct plating and growth at 4°C also led to recovery of diverse genotypes, some of which were not recovered following enrichment. Strains recovered from different permafrost samples were predominantly oligotrophic and non-spore-forming but were otherwise variable from each other in terms of a number of bacteriological characteristics. Our data suggest that a combination of isolation protocols from different permafrost samples should be used to establish a culture-based survey of the different bacterial lineages in permafrost. Received: December 2, 1999 / Accepted: January 18, 2000     

Evolutionary rate variation in Old World monkeys

We analysed over 8 million base pairs of bacterial artificial chromosome-based sequence alignments of four Old World monkeys and the human genome. Our findings are as follows. (i) Genomic divergences among several Old World monkeys mirror those between well-studied hominoids. (ii) The X-chromosome evolves slower than autosomes, in accord with ‘male-driven evolution’. However, the degree of male mutation bias is lower in Old World monkeys than in hominoids. (iii) Evolutionary rates vary significantly between lineages. The baboon branch shows a particularly slow molecular evolution. Thus, lineage-specific evolutionary rate variation is a common theme of primate genome evolution. (iv) In contrast to the overall pattern, mutations originating from DNA methylation exhibit little variation between lineages. Our study illustrates the potential of primates as a model system to investigate genome evolution, in particular to elucidate molecular mechanisms of substitution rate variation.     

6S RNA is a widespread regulator of eubacterial RNA polymerase that resembles an open promoter

6S RNA is an abundant noncoding RNA in Escherichia coli that binds to sigma70 RNA polymerase holoenzyme to globally regulate gene expression in response to the shift from exponential growth to stationary phase. We have computationally identified >100 new 6S RNA homologs in diverse eubacterial lineages. Two abundant Bacillus subtilis RNAs of unknown function (BsrA and BsrB) and cyanobacterial 6Sa RNAs are now recognized as 6S homologs. Structural probing of E. coli 6S RNA and a B. subtilis homolog supports a common secondary structure derived from comparative sequence analysis. The conserved features of 6S RNA suggest that it binds RNA polymerase by mimicking the structure of DNA template in an open promoter complex. Interestingly, the two B. subtilis 6S RNAs are discoordinately expressed during growth, and many proteobacterial 6S RNAs could be cotranscribed with downstream homologs of the E. coli ygfA gene encoding a putative methenyltetrahydrofolate synthetase. The prevalence and robust expression of 6S RNAs emphasize their critical role in bacterial adaptation.     

Quantitative approach in the study of adhesion of lactic acid bacteria to intestinal cells and their competition with enterobacteria

To describe the phenomena of bacterial adhesion to intestinal cells and the competition for adhesion between bacteria, mathematical equations based on a simple dissociation process involving a finite number of bacterial receptors on intestinal cell surface were developed. The equations allow the estimation of the maximum number of Lactobacillus sp. and Escherichia coli cells that can adhere to Caco-2 cells and intestinal mucus; they also characterize the affinity of the bacteria to Caco-2 cells and intestinal and fecal mucus and the theoretical adhesion ratio of two bacteria present in a mixed suspension. The competition for adhesion between Lactobacillus rhamnosus GG and E. coli TG1 appeared to follow the proposed kinetics, whereas the competition between Lactobacillus casei Shirota and E. coli TG1 may involve multiple adhesion sites or a soluble factor in the culture medium of the former. The displacement of the adhered Lactobacillus by E. coli TG1 seemed to be a rapid process, whereas the displacement of E. coli TG1 by the Lactobacillus took more than an hour.     

Evidence that a major determinant for the identity of a transfer RNA is conserved in evolution

We observed recently that a single G3.U70 base pair in the amino acid acceptor stem of an Escherichia coli alanine tRNA is a major determinant for its identity. Inspection of tRNA sequences shows that G3.U70 is unique to alanine in E. coli and is present in eucaryotic cytoplasmic alanine tRNAs. We show here that single nucleotide changes of G3.U70 to A3.U70 or to G3.C70 eliminate in vitro aminoacylation of an insect and of a human alanine tRNA by the respective homologous synthetase. Compared to the influence of G3.U70, other sequence variations in tRNAAla have a relatively small effect on aminoacylation by the insect and human enzymes. In addition, while these eucaryotic tRNAs have nucleotide differences from E. coli alanine tRNA, they are heterologously charged only with alanine when expressed in E. coli. The results indicate a functional role for G3.U70 that is conserved in evolution. They also suggest that the sequence differences between E. coli and the eucaryotic alanine tRNAs at sites other than the conserved G3.U70 do not create major determinants for recognition by any other bacterial enzyme.     

The Cockroach Origin of the Termite Gut Microbiota: Patterns in Bacterial Community Structure Reflect Major Evolutionary Events

Termites digest wood and other lignocellulosic substrates with the help of their intestinal microbiota. While the functions of the symbionts in the digestive process are slowly emerging, the origin of the bacteria colonizing the hindgut bioreactor is entirely unknown. Recently, our group discovered numerous representatives of bacterial lineages specific to termite guts in a closely related omnivorous cockroach, but it remains unclear whether they derive from the microbiota of a common ancestor or were independently selected by the gut environment. Here, we studied the bacterial gut microbiota in 34 species of termites and cockroaches using pyrotag analysis of the 16S rRNA genes. Although the community structures differed greatly between the major host groups, with dramatic changes in the relative abundances of particular bacterial taxa, we found that the majority of sequence reads belonged to bacterial lineages that were shared among most host species. When mapped onto the host tree, the changes in community structure coincided with major events in termite evolution, such as acquisition and loss of cellulolytic protists and the ensuing dietary diversification. UniFrac analysis of the core microbiota of termites and cockroaches and construction of phylogenetic tree of individual genus level lineages revealed a general host signal, whereas the branching order often did not match the detailed phylogeny of the host. It remains unclear whether the lineages in question have been associated with the ancestral cockroach since the early Cretaceous (cospeciation) or are diet-specific lineages that were independently acquired from the environment (host selection).