Adaptive mgl-regulatory mutations and genetic diversity evolving in glucose-limited Escherichia coli populations

The mutational adaptation of E. coli to low glucose concentrations was studied in chemostats over 280 generations of growth. All members of six independent populations acquired increased fitness through the acquisition of mutations at the mgl locus, increasing the binding protein-dependent transport of glucose. These mutations provided a strong fitness advantage (up to 10-fold increase in glucose affinity) and were present in most isolates after 140 generations. mgl constitutivity in some isolates was caused by base substitution, short duplication, small deletion and IS1 insertion in the 1041 bp gene encoding the repressor of the mgl system, mglD (galS). But an unexpectedly large proportion of mutations were located in the short mgl operator sequence (mglO), and the majority of mutations were in mglO after 280 generations of selection. The adaptive mglO substitutions in several independent populations were at exactly the positions conserved in the two 8 bp half-sites of the mgl operator, with the nature of the base changes also completely symmetrical. Either mutations were directed to the operator or the particular operator mutations had a selective advantage under glucose limitation. Indeed, isolates carrying mglO mutations showed greater rates of transport for glucose and galactose at low concentrations than those carrying mglD null mutations. mglO mutations avoid cross-talk by members of the GalR-Lacl repressor family, reducing transporter expression and providing a competitive advantage in a glucose-limited environment. Another interesting aspect of these results was that each adapted population acquired multiple mgl alleles, with several populations containing at least six different mgl-regulatory mutations co-existing after 200 generations. The diversity of mutations in the mglO/mglD region, generally in combination with mutations at other loci regulating glucose uptake (malT, mlc, ptsG), provided evidence for multiple clones in each population. Increased fitness was accompanied by the generation of genetic diversity and not the evolution of a single winner clone, as predicted by the periodic selection model of bacterial populations.     

The influence of cellular physiology on the initiation of mutational pathways in Escherichia coli populations

The factors affecting the direction of evolutionary pathways and the reproducibility of adaptive responses were investigated under closely related but non-identical conditions. Replicate chemostat cultures of Escherichia coli were compared when adapting to partial or severe glucose limitation. Four independent populations used a reproducible sequence of early mutational changes under both conditions, with rpoS mutations always occurring first before mgl. However, there were interesting differences in the timing of mutational sweeps: rpoS mutations appeared in a clock-like fashion under both partial and severe glucose limitation, while mgl sweeps arose under both conditions but at different times. Interestingly, malT and mlc mutations appeared only under severe limitation. Even though the ancestors were genotypically identical, the semi-differentiated properties of bacteria growing with mild or severe glucose limitation sent the populations in characteristic directions. Mutation supply and the fitness contribution of mutations were estimated and demonstrated to be potential influences in the choice of particular adaptation pathways under severe and mild glucose limitation. Predicting all the mutations fixed in adapting populations is beyond our current understanding of evolutionary processes, but the interplay between ancestor physiology and the initiation of adaptation pathways is demonstrated and definable in bacterial populations.     

Experimental analysis of molecular events during mutational periodic selections in bacterial evolution

A fundamental feature of bacterial evolution is a succession of adaptive mutational sweeps when fitter mutants take over a population. To understand the processes involved in mutational successions, Escherichia coli continuous cultures were analyzed for changes at two loci where mutations provide strong transport advantages to fitness under steady-state glucose limitation. Three separate sweeps, observed as classic periodic selection events causing a change in the frequency of neutral mutations (in fhuA causing phage T5 resistance), were identified with changes at particular loci. Two of the sweeps were associated with a reduction in the frequency of neutral mutations and the concurrent appearance of at least 13 alleles at the mgl or mlc loci, respectively. These mgl and mlc polymorphisms were of many mutational types, so were not the result of a mutator or directed mutation event. The third sweep observed was altogether distinct and involved hitchhiking between T5 resistance and advantageous mgl mutations. Moreover, the hitchhiking event coincided with an increase in mutation rates, due to the transient appearance of a strong mutator in the population. The spectrum of mgl mutations among mutator isolates was distinct and due to mutS. The mutator-associated periodic selection also resulted in mgl and fhuA polymorphism in the sweeping population. These examples of periodic selections maintained significant genotypic diversity even in a rapidly evolving culture, with no individual "winner clone" or genotype purging the population.     

OmpF changes and the complexity of Escherichia coli adaptation to prolonged lactose limitation

The disaccharide lactose has no specific diffusion pathway across the outer membrane of Escherichia coli. At least three classes of spontaneous mutation affecting outer membrane permeability arose with each of three independent E. coli populations adapting to prolonged lactose limitation in chemostats. Both structural and regulatory mutations affecting OmpF porin predominated in isolates after 210-280 generations of culture. Six types of ompF mutation were found, including in-frame deletions and substitutions at Arg82 and Asp113, all affecting the channel constriction residues of OmpF. Isolates had increased susceptibility to antibiotics and were affected in the quantity of OmpF, LamB and OmpA proteins. A minimum of three or four mutations was evident in all isolates after 280 generations in a lactose-limited environment, in addition to lac mutations defined in previous studies.     

Enrichment and elimination of mutY mutators in Escherichia coli populations

The kinetics of mutator sweeps was followed in two independent populations of Escherichia coli grown for up to 350 generations in glucose-limited continuous culture. A rapid elevation of mutation rates was observed in both populations within 120-150 generations, as was apparent from major increases in the proportion of the populations with unselected mutations in fhuA. The increase in mutation rates was due to sweeps by mutY mutators. In both cultures, the enrichment of mutators resulted from hitchhiking with identified beneficial mutations increasing fitness under glucose limitation; mutY hitchhiked with mgl mutations in one culture and ptsG in the other. In both cases, mutators were enriched to constitute close to 100% of the population before a periodic selection event reduced the frequency of unselected mutations and mutators in the cultures. The high proportion of mutators persisted for 150 generations in one population but began to be eliminated within 50 generations in the other. The persistence of mutator, as well as experimental data showing that mutY bacteria were as fit as near-isogenic mutY(+) bacteria in competition experiments, suggest that mutator load by deleterious mutations did not explain the rapidly diminishing proportion of mutators in the populations. The nonmutators sweeping out mutators were also unlikely to have arisen by reversion or antimutator mutations; the mutY mutations were major deletions in each case and the bacteria sweeping out mutators contained intact mutY. By following mgl allele frequencies in one population, we discovered that mutators were outcompeted by bacteria that had rare mgl mutations previously as well as additional beneficial mutation(s). The pattern of appearance of mutY, but not its elimination, conforms to current models of mutator sweeps in bacterial populations. A mutator with a narrow mutational spectrum like mutY may be lost if the requirement for beneficial mutations is for changes other than GC --> TA transversions. Alternatively, epistatic interactions between mutator mutation and beneficial mutations need to be postulated to explain mutator elimination.     

Derepression of LamB protein facilitates outer membrane permeation of carbohydrates into Escherichia coli under conditions of nutrient stress.

The level of LamB protein in the outer membrane of Escherichia coli was derepressed in the absence of a known inducer (maltodextrins) under carbohydrate-limiting conditions in chemostats. LamB protein contributed to the ability of the bacteria to remove sugar from glucose-limited chemostats, and well-characterized lamB mutants with reduced stability constants for glucose were less growth competitive under glucose limitation than those with wild-type affinity. In turn, wild-type bacteria were less growth competitive than lamB mutants with enhanced sugar affinity. In contrast to an earlier report, we found that LamB- bacteria were less able to compete in carbohydrate-limited chemostats (with glucose, lactose, arabinose, or glycerol as the carbon and energy sources) when mixed with LamB+ bacteria. The transport Km for [14C]glucose was affected by the presence or affinity of LamB, but only in chemostat-grown bacteria, with their elevated LamB levels. The pattern of expression of LamB and the advantage it confers for growth on low concentrations of carbohydrates are consistent with a wider role in sugar permeation than simply maltosaccharide transport, and hence the well-known maltoporin activity of LamB is but one facet of its role as the general glycoporin of E. coli. A corollary of these findings is that OmpF/OmpC porins, present at high levels in carbon-limited bacteria, do not provide sufficient permeability to sugars or even glycerol to support high growth rates at low concentrations. Hence, the sugar-binding site of LamB protein is an important contributor to the permeability of the outer membrane to carbohydrates in habitats with low extracellular nutrient concentrations.     

Allele-specific malE mutations that restore interactions between maltose-binding protein and the inner-membrane components of the maltose transport system

Active accumulation of maltose and maltodextrins by Escherichia coli depends on an outer-membrane protein. LamB, a periplasmic maltose-binding protein (MalE, MBP) and three inner-membrane proteins, MalF, MalG and MalK. MalF and MalG are integral transmembrane proteins, while MalK is associated with the inner aspect of the cytoplasmic membrane via an interaction with MalG. Previously we have shown that MBP is essential for movement of maltose across the inner membrane. We have taken advantage of malF and malG mutants in which MBP interacts improperly with the membrane proteins. We describe the properties of malE mutations in which a proper interaction between MBP and defective MalF and MalG proteins has been restored. We found that these malE suppressor mutations are able to restore transport activity in an allele-specific manner. That is, a given malE mutation restores transport activity to different extents in different malF and malG mutants. Since both malF and malG mutations could be suppressed by allele-specific malE suppressors, we propose that, in wild-type bacteria, MBP interacts with sites on both MalF and MalG during active transport. The locations of different malE suppressor mutations indicate specific regions on MBP that are important for interacting with MalF and MalG.     

Osmoregulation of the maltose regulon in Escherichia coli.

The maltose regulon consists of four operons that direct the synthesis of proteins required for the transport and metabolism of maltose and maltodextrins. Expression of the mal genes is induced by maltose and maltodextrins and is dependent on a specific positive regulator, the MalT protein, as well as on the cyclic AMP-catabolite gene activator protein complex. In the absence of an exogenous inducer, expression of the mal regulon was greatly reduced when the osmolarity of the growth medium was high; maltose-induced expression was not affected, and malTc-dependent expression was only weakly affected. Mutants lacking MalK, a cytoplasmic membrane protein required for maltose transport, expressed the remaining mal genes at a high level, presumably because an internal inducer of the mal system accumulated; this expression was also strongly repressed at high osmolarity. The repression of mal regulon expression at high osmolarity was not caused by reduced expression of the malT, envZ, or crp gene or by changes in cellular cyclic AMP levels. In strains carrying mutations in genes encoding amylomaltase (malQ), maltodextrin phosphorylase (malP), amylase (malS), or glycogen (glg), malK mutations still led to elevated expression at low osmolarity. The repression at high osmolarity no longer occurred in malQ mutants, however, provided that glycogen was present.     

Effects of chemical contaminants on genetic diversity in natural populations: implications for biomonitoring and ecotoxicology

The conservation of genetic diversity has emerged as one of the central issues in conservation biology. Although researchers in the areas of evolutionary biology, population management, and conservation biology routinely investigate genetic variability in natural populations, only a handful of studies have addressed the effects of chemical contamination on population genetics. Chemical contamination can cause population reduction by the effects of somatic and heritable mutations, as well as non-genetic modes of toxicity. Stochastic processes in small populations, increased mutation load, and the phenomenon of mutational meltdown are compounding factors that cause reduced fitness and accelerate the process of population extirpation. Although the original damage caused by chemical contaminants is at the molecular level, there are emergent effects at the level of populations, such as the loss of genetic diversity, that are not predictable based solely on knowledge of the mechanism of toxicity of the chemical contaminants. Therefore, the study of evolutionary toxicology, which encompasses the population-genetic effects of environmental contaminants, should be an important focus of ecotoxicology. This paper reviews the issues surrounding the genetic effects of pollution, summarizes the technical approaches that can be used to address these issues, and provides examples of studies that have addressed some of them.     

The malX malY operon of Escherichia coli encodes a novel enzyme II of the phosphotransferase system recognizing glucose and maltose and an enzyme abolishing the endogenous induction of the maltose system.

Mutants lacking MalK, a subunit of the binding protein-dependent maltose-maltodextrin transport system, constitutively express the maltose genes. A second site mutation in malI abolishes the constitutive expression. The malI gene (at 36 min on the linkage map) codes for a typical repressor protein that is homologous to the Escherichia coli LacI, GalR, or CytR repressor (J. Reidl, K. R?misch, M. Ehrmann, and W. Boos, J. Bacteriol. 171:4888-4899, 1989). We now report that MalI regulates an adjacent and divergently oriented operon containing malX and malY. MalX encodes a protein with a molecular weight of 56,654, and the deduced amino acid sequence of MalX exhibits 34.9% identity to the enzyme II of the phosphototransferase system for glucose (ptsG) and 32.1% identity to the enzyme II for N-acetylglucosamine (nagE). When constitutively expressed, malX can complement a ptsG ptsM double mutant for growth on glucose. Also, a delta malE malT(Con) strain that is unable to grow on maltose due to its maltose transport defect becomes Mal+ after introduction of malI::Tn10 and the plasmid carrying malX. MalX-mediated transport of glucose and maltose is likely to occur by facilitated diffusion. We conclude that malX encodes a phosphotransferase system enzyme II that can recognize glucose and maltose as substrates even though these sugars may not represent the natural substrates of the system. The second gene in the operon, malY, encodes a protein of 43,500 daltons. Its deduced amino acid sequence exhibits weak homology to aminotransferase sequences. The presence of plasmid-encoded MalX alone was sufficient for complementing growth on glucose in a ptsM ptsG glk mutant, and the plasmid-encoded MalY alone was sufficient to abolish the constitutivity of the mal genes in a malK mutant. The overexpression of malY in a strain that is wild type with respect to the maltose genes strongly interferes with growth on maltose. This is not the case in a malT(Con) strain that expresses the mal genes constitutively. We conclude that malY encodes an enzyme that degrades the inducer of the maltose system or prevents its synthesis.     

Distribution of mannosamine and mannosaminuronic acid among cell walls of Bacillus species

Some Escherichia coli K-12 lamB mutants, those producing reduced amounts of LamB protein (one-tenth the wild type amount), grow normally on dextrins but transport maltose when present at a concentration of 1 microM at about one-tenth the normal rate. lamB Dex- mutants were found as derivatives of these strains. These Dex- mutants are considerably impaired in the transport of maltose at low concentrations (below 10 microM), and they have a structurally altered LamB protein which is impaired in its interaction with phages lambda and K10 but still interacts with a lambda host range mutant lambda hh*. The Dex- mutants are double lamB mutants carrying one mutation, already present in the parental strains, that reduces LamB synthesis and a second that alters LamB structure. The secondary mutations, present in different independent Dex- mutants, are clustered in the same region of the lamB gene. Dex+ revertants were isolated and analyzed: when the altered LamB protein is made in wild-type amount, due to a reversion of the first mutation, the phenotype reverts to Dex+. However, these Dex+ revertants are still very significantly impaired in maltose transport at low concentrations (below 10 microM).     

Epistasis in polygenic traits and the evolution of genetic architecture under stabilizing selection

We consider the effects of epistasis in a polygenic trait in the balance of mutation and stabilizing selection. The main issues are the genetic variation maintained in equilibrium and the evolution of the mutational effect distribution. The model assumes symmetric mutation and a continuum of alleles at all loci. Epistasis is modeled proportional to pairwise products of the single-locus effects. A general analytical formalism is developed. Assuming linkage equilibrium, we derive results for the equilibrium mutation load and the genetic and mutational variance in the house of cards and the Gaussian approximation. The additive genetic variation maintained in mutation-selection balance is reduced by any pattern of the epistatic interactions. The mutational variance, in contrast, is often increased. Large differences in mutational effects among loci emerge, and a negative correlation among (standard mean) locus mutation effects and mutation rates is predicted. Contrary to the common view since Waddington, we find that stabilizing selection in general does not lead to canalization of the trait. We propose that canalization as a target of selection instead occurs at the genic level. Here, primarily genes with a high mutation rate are buffered, often at the cost of decanalization of other genes. An intuitive interpretation of this view is given in the discussion.     

Stabilizing selection,mutational bias,and the evolution of sex*

Stabilizing selection around a fixed phenotypic optimum is expected to disfavor sexual reproduction, since asexually reproducing organisms can maintain a higher fitness at equilibrium, while sex disrupts combinations of compensatory mutations. This conclusion rests on the assumption that mutational effects on phenotypic traits are unbiased, that is, mutation does not tend to push phenotypes in any particular direction. In this article, we consider a model of stabilizing selection acting on an arbitrary number of polygenic traits coded by bialellic loci, and show that mutational bias may greatly reduce the mean fitness of asexual populations compared with sexual ones in regimes where mutations have weak to moderate fitness effects. Indeed, mutation and drift tend to push the population mean phenotype away from the optimum, this effect being enhanced by the low effective population size of asexual populations. In a second part, we present results from individual‐based simulations showing that positive rates of sex are favored when mutational bias is present, while the population evolves toward complete asexuality in the absence of bias. We also present analytical (QLE) approximations for the selective forces acting on sex in terms of the effect of sex on the mean and variance in fitness among offspring.     

Mutation and evolution of microsatellite loci in Neurospora

The patterns of mutation and evolution at 13 microsatellite loci were studied in the filamentous fungal genus Neurospora. First, a detailed investigation was performed on five microsatellite loci by sequencing each microsatellite, together with its nonrepetitive flanking regions, from a set of 147 individuals from eight species of Neurospora. To elucidate the genealogical relationships among microsatellite alleles, repeat number was mapped onto trees constructed from flanking-sequence data. This approach allowed the potentially convergent microsatellite mutations to be placed in the evolutionary context of the less rapidly evolving flanking regions, revealing the complexities of the mutational processes that have generated the allelic diversity conventionally assessed in population genetic studies. In addition to changes in repeat number, frequent substitution mutations within the microsatellites were detected, as were substitutions and insertion/deletions within the flanking regions. By comparing microsatellite and flanking-sequence divergence, clear evidence of interspecific allele length homoplasy and microsatellite mutational saturation was observed, suggesting that these loci are not appropriate for inferring phylogenetic relationships among species. In contrast, little evidence of intraspecific mutational saturation was observed, confirming the utility of these loci for population-level analyses. Frequency distributions of alleles within species were generally consistent with the stepwise mutational model. By comparing variation within species at the microsatellites and the flanking-sequence, estimated microsatellite mutation rates were approximately 2500 times greater than mutation rates of flanking DNA and were consistent with estimates from yeast and fruit flies. A positive relationship between repeat number and variance in repeat number was significant across three genealogical depths, suggesting that longer microsatellite alleles are more mutable than shorter alleles. To test if the observed patterns of microsatellite variation and mutation could be generalized, an additional eight microsatellite loci were characterized and sequenced from a subset of the same Neurospora individuals.     

Coexistence and Within-Host Evolution of Diversified Lineages of Hypermutable Pseudomonas aeruginosa in Long-term Cystic Fibrosis Infections

The advent of high-throughput sequencing techniques has made it possible to follow the genomic evolution of pathogenic bacteria by comparing longitudinally collected bacteria sampled from human hosts. Such studies in the context of chronic airway infections by Pseudomonas aeruginosa in cystic fibrosis (CF) patients have indicated high bacterial population diversity. Such diversity may be driven by hypermutability resulting from DNA mismatch repair system (MRS) deficiency, a common trait evolved by P. aeruginosa strains in CF infections. No studies to date have utilized whole-genome sequencing to investigate within-host population diversity or long-term evolution of mutators in CF airways. We sequenced the genomes of 13 and 14 isolates of P. aeruginosa mutator populations from an Argentinian and a Danish CF patient, respectively. Our collection of isolates spanned 6 and 20 years of patient infection history, respectively. We sequenced 11 isolates from a single sample from each patient to allow in-depth analysis of population diversity. Each patient was infected by clonal populations of bacteria that were dominated by mutators. The in vivo mutation rate of the populations was ∼100 SNPs/year–∼40-fold higher than rates in normo-mutable populations. Comparison of the genomes of 11 isolates from the same sample showed extensive within-patient genomic diversification; the populations were composed of different sub-lineages that had coexisted for many years since the initial colonization of the patient. Analysis of the mutations identified genes that underwent convergent evolution across lineages and sub-lineages, suggesting that the genes were targeted by mutation to optimize pathogenic fitness. Parallel evolution was observed in reduction of overall catabolic capacity of the populations. These findings are useful for understanding the evolution of pathogen populations and identifying new targets for control of chronic infections.     

Optimal bacteriophage mutation rates for phage therapy

The mutability of bacteriophages offers a particular advantage in the treatment of bacterial infections not afforded by other antimicrobial therapies. When phage-resistant bacteria emerge, mutation may generate phage capable of exploiting and thus limiting population expansion among these emergent types. However, while mutation potentially generates beneficial variants, it also contributes to a genetic load of deleterious mutations. Here, we model the influence of varying phage mutation rate on the efficacy of phage therapy. All else being equal, phage types with historical mutation rates of approximately 0.1 deleterious mutations per genome per generation offer a reasonable balance between beneficial mutational diversity and deleterious mutational load. We determine that increasing phage inoculum density can undesirably increase the peak density of a mutant bacterial class by limiting the in situ production of mutant phage variants. For phage populations with minimal genetic load, engineering mutation rate increases beyond the mutation-selection balance optimum may provide even greater protection against emergent bacterial types, but only with very weak selective coefficients for de novo deleterious mutations (below approximately 0.01). Increases to the mutation rate beyond the optimal value at mutation-selection balance may therefore prove generally undesirable.     

Evaluating the within-host fitness effects of mutations fixed during virus adaptation to different ecotypes of a new host

The existence of genetic variation for resistance in host populations is assumed to be essential to the spread of an emerging virus. Models predict that the rate of spread slows down with the increasing frequency and higher diversity of resistance alleles in the host population. We have been using the experimental pathosystem Arabidopsis thaliana—tobacco etch potyvirus (TEV) to explore the interplay between genetic variation in host''s susceptibility and virus diversity. We have recently shown that TEV populations evolving in A. thaliana ecotypes that differ in susceptibility to infection gained within-host fitness, virulence and infectivity in a manner compatible with a gene-for-gene model of host–parasite interactions: hard-to-infect ecotypes were infected by generalist viruses, whereas easy-to-infect ecotypes were infected by every virus. We characterized the genomes of the evolved viruses and found cases of host-driven convergent mutations. To gain further insights in the mechanistic basis of this gene-for-gene model, we have generated all viral mutations individually as well as in specific combinations and tested their within-host fitness effects across ecotypes. Most of these mutations were deleterious or neutral in their local ecotype and only a very reduced number had a host-specific beneficial effect. We conclude that most of the mutations fixed during the evolution experiment were so by drift or by selective sweeps along with the selected driver mutation. In addition, we evaluated the ruggedness of the underlying adaptive fitness landscape and found that mutational effects were mostly multiplicative, with few cases of significant epistasis.