Increased production of L-serine in Escherichia coli through Adaptive Laboratory Evolution

L-serine is a promising building block biochemical with a high theoretical production yield from glucose. Toxicity of L-serine is however prohibitive for high-titer production in E. coli. Here, E. coli lacking L-serine degradation pathways was evolved for improved tolerance by gradually increasing L-serine concentration from 3 to 100 g/L using adaptive laboratory evolution (ALE). Genome sequencing of isolated clones revealed multiplication of genetic regions, as well as mutations in thrA, thereby showing a potential mechanism of serine inhibition. Additional mutations were evaluated by MAGE combined with amplicon sequencing, revealing role of rho, lrp, pykF, eno, and rpoB on tolerance and fitness in minimal medium. Production using the tolerant strains resulted in 37 g/L of L-serine with a 24% mass yield. The resulting titer is similar to the highest production reported for any organism thereby highlighting the potential of ALE for industrial biotechnology.     

Genetic Architecture of Intrinsic Antibiotic Susceptibility

Background

Antibiotic exposure rapidly selects for more resistant bacterial strains, and both a drug''s chemical structure and a bacterium''s cellular network affect the types of mutations acquired.

Methodology/Principal Findings

To better characterize the genetic determinants of antibiotic susceptibility, we exposed a transposon-mutagenized library of Escherichia coli to each of 17 antibiotics that encompass a wide range of drug classes and mechanisms of action. Propagating the library for multiple generations with drug concentrations that moderately inhibited the growth of the isogenic parental strain caused the abundance of strains with even minor fitness advantages or disadvantages to change measurably and reproducibly. Using a microarray-based genetic footprinting strategy, we then determined the quantitative contribution of each gene to E. coli''s intrinsic antibiotic susceptibility. We found both loci whose removal increased general antibiotic tolerance as well as pathways whose down-regulation increased tolerance to specific drugs and drug classes. The beneficial mutations identified span multiple pathways, and we identified pairs of mutations that individually provide only minor decreases in antibiotic susceptibility but that combine to provide higher tolerance.

Conclusions/Significance

Our results illustrate that a wide-range of mutations can modulate the activity of many cellular resistance processes and demonstrate that E. coli has a large mutational target size for increasing antibiotic tolerance. Furthermore, the work suggests that clinical levels of antibiotic resistance might develop through the sequential accumulation of chromosomal mutations of small individual effect.     

Protein Homeostasis Imposes a Barrier on Functional Integration of Horizontally Transferred Genes in Bacteria

Horizontal gene transfer (HGT) plays a central role in bacterial evolution, yet the molecular and cellular constraints on functional integration of the foreign genes are poorly understood. Here we performed inter-species replacement of the chromosomal folA gene, encoding an essential metabolic enzyme dihydrofolate reductase (DHFR), with orthologs from 35 other mesophilic bacteria. The orthologous inter-species replacements caused a marked drop (in the range 10–90%) in bacterial growth rate despite the fact that most orthologous DHFRs are as stable as E.coli DHFR at 37°C and are more catalytically active than E. coli DHFR. Although phylogenetic distance between E. coli and orthologous DHFRs as well as their individual molecular properties correlate poorly with growth rates, the product of the intracellular DHFR abundance and catalytic activity (k _cat/K_M), correlates strongly with growth rates, indicating that the drop in DHFR abundance constitutes the major fitness barrier to HGT. Serial propagation of the orthologous strains for ~600 generations dramatically improved growth rates by largely alleviating the fitness barriers. Whole genome sequencing and global proteome quantification revealed that the evolved strains with the largest fitness improvements have accumulated mutations that inactivated the ATP-dependent Lon protease, causing an increase in the intracellular DHFR abundance. In one case DHFR abundance increased further due to mutations accumulated in folA promoter, but only after the lon inactivating mutations were fixed in the population. Thus, by apparently distinguishing between self and non-self proteins, protein homeostasis imposes an immediate and global barrier to the functional integration of foreign genes by decreasing the intracellular abundance of their products. Once this barrier is alleviated, more fine-tuned evolution occurs to adjust the function/expression of the transferred proteins to the constraints imposed by the intracellular environment of the host organism.     

Enhancing protein perdeuteration by experimental evolution of Escherichia coli K‐12 for rapid growth in deuterium‐based media

Deuterium is a natural low abundance stable hydrogen isotope that in high concentrations negatively affects growth of cells. Here, we have studied growth of Escherichia coli MG1655, a wild‐type laboratory strain of E. coli K‐12, in deuterated glycerol minimal medium. The growth rate and final biomass in deuterated medium is substantially reduced compared to cells grown in ordinary medium. By using a multi‐generation adaptive laboratory evolution‐based approach, we have isolated strains that show increased fitness in deuterium‐based growth media. Whole‐genome sequencing identified the genomic changes in the obtained strains and show that there are multiple routes to genetic adaptation to growth in deuterium‐based media. By screening a collection of single‐gene knockouts of nonessential genes, no specific gene was found to be essential for growth in deuterated minimal medium. Deuteration of proteins is of importance for NMR spectroscopy, neutron protein crystallography, neutron reflectometry, and small angle neutron scattering. The laboratory evolved strains, with substantially improved growth rate, were adapted for recombinant protein production by T7 RNA polymerase overexpression systems and shown to be suitable for efficient production of perdeuterated soluble and membrane proteins for structural biology applications.     

Chromosomal position shift of a regulatory gene alters the bacterial phenotype

Recent studies strongly suggest that in bacterial cells the order of genes along the chromosomal origin-to-terminus axis is determinative for regulation of the growth phase-dependent gene expression. The prediction from this observation is that positional displacement of pleiotropic genes will affect the genetic regulation and hence, the cellular phenotype. To test this prediction we inserted the origin-proximal dusB-fis operon encoding the global regulator FIS in the vicinity of replication terminus on both arms of the Escherichia coli chromosome. We found that the lower fis gene dosage in the strains with terminus-proximal dusB-fis operons was compensated by increased fis expression such that the intracellular concentration of FIS was homeostatically adjusted. Nevertheless, despite unchanged FIS levels the positional displacement of dusB-fis impaired the competitive growth fitness of cells and altered the state of the overarching network regulating DNA topology, as well as the cellular response to environmental stress, hazardous substances and antibiotics. Our finding that the chromosomal repositioning of a regulatory gene can determine the cellular phenotype unveils an important yet unexplored facet of the genetic control mechanisms and paves the way for novel approaches to manipulate bacterial physiology.     

Evolution of a regulated operon in the laboratory

The evolution of new metabolic functions is being studied in the laboratory using the EBG system of E. coli as a model system. It is demonstrated that the evolution of lactose utilization by lacZ deletion strains requires a series of structural and regulatory gene mutations. Two structural gene mutations act to increase the activity of ebg enzyme toward lactose, and to permit ebg enzyme to convert lactose into allolactose, an inducer of the lac operon. A regulatory mutation increases the sensitivity of the ebg repressor to lactose, and permits sufficient ebg enzyme activity for growth. The resulting fully evolved ebg operon regulates its own expression, and also regulates the synthesis of the lactose permease.     

Variation in the Fitness Effects of Mutations with Population Density and Size in Escherichia coli

The fitness effects of mutations are context specific and depend on both external (e.g., environment) and internal (e.g., cellular stress, genetic background) factors. The influence of population size and density on fitness effects are unknown, despite the central role population size plays in the supply and fixation of mutations. We addressed this issue by comparing the fitness of 92 Keio strains (Escherichia coli K12 single gene knockouts) at comparatively high (1.2×10⁷ CFUs/mL) and low (2.5×10² CFUs/mL) densities, which also differed in population size (high: 1.2×10⁸; low: 1.25×10³). Twenty-eight gene deletions (30%) exhibited a fitness difference, ranging from 5 to 174% (median: 35%), between the high and low densities. Our analyses suggest this variation among gene deletions in fitness responses reflected in part both gene orientation and function, of the gene properties we examined (genomic position, length, orientation, and function). Although we could not determine the relative effects of population density and size, our results suggest fitness effects of mutations vary with these two factors, and this variation is gene-specific. Besides being a mechanism for density-dependent selection (r-K selection), the dependence of fitness effects on population density and size has implications for any population that varies in size over time, including populations undergoing evolutionary rescue, species invasions into novel habitats, and cancer progression and metastasis. Further, combined with recent advances in understanding the roles of other context-specific factors in the fitness effects of mutations, our results will help address theoretical and applied biological questions more realistically.     

Proteomic analysis of a multi-resistant clinical Escherichia coli isolate of unknown genomic background

Horizontal transfer of gene clusters occurs in Escherichia coli (E. coli), which could lead to evolution of new pathovars and improve survival fitness. However, this genetic event results in genomic plasticity which is a hindrance for proteomic characterization of strains with unknown genetic backgrounds. To characterize such isolate with many specific genetic variations we used the recently in-house designed MSMSpdbb software which merges protein databases from several sources of E. coli including type strains and other commensal and pathogenic isolates. We selected a multidrug resistant clinical isolate in order to check the capacity of our approach to identify selected protein markers. From the 1596 identified proteins, we found important virulence factors such as IutA, OmpA, TraT and selected enzymes conferring antibiotic resistance, such as CTX-M-15 (Extended-Spectrum Beta Lactamase - ESBL) and AAC(6′)-Ib-cr (to aminoglycoside + fluoroquinolone). In addition, we compared the protein identifications with E. coli gene annotation and found that 27% of the proteins identified in the present study corresponded to the pan-genome of E. coli species and are only present in a subset of strains. This demonstrates the ability of our approach to characterize the proteome of bacterial strains with complex genomic plasticity even without its genomic information.     

Oxidative phosphorylation in Escherichia coli K-12: The genetic and biochemical characterisation of a strain carrying a mutation in the uncB gene

1. A mutant strain of Escherichia coli unable to grow with succinate as sole carbon source was isolated. This mutant was found to carry a mutation in a gene (designated uncB) mapping at about minute 73.5 on the E. coli chromosome and close to the uncA gene which is probably the structural gene for (Mg²⁺,Ca²⁺)-stimulated ATPase.2. The uncB401 allele was transduced into two other strains of E. coli and the transductants compared with the parent strains.3. Strains carrying the uncB401 allele have low aerobic growth yields when grown on limiting concentrations of glucose, but unlike mutations in the uncA gene, mutations in the uncB gene do not impair anaerobic growth on a glucose-mineral salts medium.4. Oxidase activities in membranes from the normal strains and strains carrying the uncB401 allele were similar.5. Measurement of P/O ratios indicated that a mutation in the uncB gene causes uncoupling of phosphorylation associated with electron transport with d-lactate as substrate.6. (Mg²⁺,Ca²⁺)-stimulated ATPase activities in the normal strains and in strains carrying the uncB401 allele are similar.7. Estimation of the energy-linked and non-energy-linked transhydrogenase activities in membrane preparations from both the normal and mutant strains indicated that the protein affected by a mutation in the uncB gene is essential for the functioning of the ATP-dependent energy-linked transhydrogenase.8. It is concluded that two proteins, specified by the uncA and uncB genes, are essential for phosphorylation coupled to d-lactate oxidation and also for the energy-linked transhydrogenase activity using ATP as the energy source.     

Plasmidic qnrA3 enhances Escherichia coli fitness in absence of antibiotic exposure

The widespread presence of plasmid-mediated quinolone resistance determinants, particularly qnr genes, has become a current issue. By protecting DNA-gyrase from quinolones, Qnr proteins confer a low level quinolone resistance that is not sufficient to explain their emergence. Since Qnr proteins were hypothesized to act as DNA-binding protein regulators, qnr genes could have emerged by providing a selective advantage other than antibiotic resistance. We investigated host fitness of Escherichia coli isogenic strains after acquisition of the qnrA3 gene, inserted either alone onto a small plasmid (pBR322), or harbored on a large conjugative native plasmid, pHe96(qnrA3) found in a clinical isolate. The isogenic strains were derived from the susceptible E. coli CFT073, a virulent B2 group strain known to infect bladder and kidneys in a mouse model of pyelonephritis. In vitro experiments included growth analysis by automatic spectrophotometry and flow cytometry, and competitions with CFU enumeration. In vivo experiments included infection with each strain and pairwise competitions in absence of antimicrobial exposure. As controls for our experiments we used mutations known to reduce fitness (rpsL K42N mutation) or to enhance fitness (tetA deletion in pBR322). E. coli CFT073 transformed with pBRAM(PBR322-qnrA3) had significantly higher maximal OD than E. coli CFT073 transformed with pBR322 or pBR322ΔtetA, and in vivo competitions were more often won by the qnrA3 carrying strain (24 victories vs. 9 loss among 42 competitions, p = 0.001). In contrast, when pHe96(qnrA3) was introduced by conjugation in E. coli CFT073, it exerted a fitness cost shown by an impaired growth observed in vitro and in vivo and a majority of lost competitions (33/35, p<0.0001). In conclusion, qnrA3 acquisition enhanced bacterial fitness, which may explain qnr emergence and suggests a regulation role of qnr. However, fitness was reduced when qnrA3 was inserted onto multidrug-resistant plasmids and this can slow down its dissemination without antibiotic exposure.     

Evolution of E. coli on [U-13C]Glucose Reveals a Negligible Isotopic Influence on Metabolism and Physiology

¹³C-Metabolic flux analysis (¹³C-MFA) traditionally assumes that kinetic isotope effects from isotopically labeled compounds do not appreciably alter cellular growth or metabolism, despite indications that some biochemical reactions can be non-negligibly impacted. Here, populations of Escherichia coli were adaptively evolved for ~1000 generations on uniformly labeled ¹³C-glucose, a commonly used isotope for ¹³C-MFA. Phenotypic characterization of these evolved strains revealed ~40% increases in growth rate, with no significant difference in fitness when grown on either labeled (¹³C) or unlabeled (¹²C) glucose. The evolved strains displayed decreased biomass yields, increased glucose and oxygen uptake, and increased acetate production, mimicking what is observed after adaptive evolution on unlabeled glucose. Furthermore, full genome re-sequencing revealed that the key genetic changes underlying these phenotypic alterations were essentially the same as those acquired during adaptive evolution on unlabeled glucose. Additionally, glucose competition experiments demonstrated that the wild-type exhibits no isotopic preference for unlabeled glucose, and the evolved strains have no preference for labeled glucose. Overall, the results of this study indicate that there are no significant differences between ¹²C and ¹³C-glucose as a carbon source for E. coli growth.     

Investigation of plasmid-mediated resistance in E. coli isolated from healthy and diarrheic sheep and goats

Escherichia coli is zoonotic bacteria and the emergence of antimicrobial-resistant strains becomes a critical issue in both human and animal health globally. This study was therefore aimed to investigate the plasmid-mediated resistance in E. coli strains isolated from healthy and diarrheic sheep and goats. A total of 234 fecal samples were obtained from 157 sheep (99 healthy and 58 diarrheic) and 77 goats (32 healthy and 45 diarrheic) for the isolation and identification of E. coli. Plasmid DNA was extracted using the alkaline lysis method. Phenotypic antibiotic susceptibility profiles were determined against the three classes of antimicrobials, which resistance is mediated by plasmids (Cephalosporins, Fluoroquinolone, and Aminoglycosides) using the disc-diffusion method. The frequency of plasmid-mediated resistance genes was investigated by PCR. A total of 159 E. coli strains harbored plasmids. The isolates antibiogram showed different patterns of resistance in both healthy and diarrheic animals. A total of (82; 51.5%) E. coli strains were multidrug-resistant. rmtB gene was detected in all Aminoglycoside-resistant E. coli, and the ESBL-producing E. coli possessed different CTX-M genes. Similarly, fluoroquinolone-resistant E. coli possessed different qnr genes. On the analysis of the gyrB gene sequence of fluoroquinolone-resistant E. coli, multiple point mutations were revealed. In conclusion, a high prevalence of E. coli with high resistance patterns to antimicrobials was revealed in the current study, in addition to a wide distribution of their resistance determinants. These findings highlight the importance of sheep and goats as reservoirs for the dissemination of MDR E. coli and resistance gene horizontal transfer.     

Rapid Mapping of Conditional and Auxotrophic Mutations in Escherichia coli K-12

The approximate genetic map locations of auxotrophic and conditional lethal mutations of Escherichia coli can be rapidly determined with replica plating techniques. A set of patches of 15 streptomycin-sensitive (Str^S) Hfr strains with points of origin distributed around the map is replica plated onto a recombinant-selective plate with a lawn of Str^R cells which carry an unmapped mutation. The map interval defined by the Hfr points of origin which are closest to the mutant locus is seen by the presence or absence of heavy patches of recombinants produced by transfer of early wild-type genes from the Hfrs. An alternative method is to replicate patches of different mutant strains (100 per plate) onto Hfr lawns; in this case more than 1,000 different mutants can be mapped in a single experiment in a few days. In this way, many types of mutations with similar phenotypes can be grouped as to approximate location on the genetic map. For ordering mutations within groups, the same replica plating methods can be used to cross F-prime derivatives of mutants with other mutants of the same group. Relative merits of these and other mapping methods of E. coli are discussed.     

Genome-Wide Gene Expression Analysis of Bordetella pertussis Isolates Associated with a Resurgence in Pertussis: Elucidation of Factors Involved in the Increased Fitness of Epidemic Strains

Bordetella pertussis (B. pertussis) is the causative agent of whooping cough, which is a highly contagious disease in the human respiratory tract. Despite vaccination since the 1950s, pertussis remains the most prevalent vaccine-preventable disease in developed countries. A recent resurgence pertussis is associated with the expansion of B. pertussis strains with a novel allele for the pertussis toxin (ptx) promoter ptxP3 in place of resident ptxP1 strains. The recent expansion of ptxP3 strains suggests that these strains carry mutations that have increased their fitness. Compared to the ptxP1 strains, ptxP3 strains produce more Ptx, which results in increased virulence and immune suppression. In this study, we investigated the contribution of gene expression changes of various genes on the increased fitness of the ptxP3 strains. Using genome-wide gene expression profiling, we show that several virulence genes had higher expression levels in the ptxP3 strains compared to the ptxP1 strains. We provide the first evidence that wildtype ptxP3 strains are better colonizers in an intranasal mouse infection model. This study shows that the ptxP3 mutation and the genetic background of ptxP3 strains affect fitness by contributing to the ability to colonize in a mouse infection model. These results show that the genetic background of ptxP3 strains with a higher expression of virulence genes contribute to increased fitness.     

The basis of antagonistic pleiotropy in hfq mutations that have opposite effects on fitness at slow and fast growth rates

Mutations beneficial in one environment may cause costs in different environments, resulting in antagonistic pleiotropy. Here, we describe a novel form of antagonistic pleiotropy that operates even within the same environment, where benefits and deleterious effects exhibit themselves at different growth rates. The fitness of hfq mutations in Escherichia coli affecting the RNA chaperone involved in small-RNA regulation is remarkably sensitive to growth rate. E. coli populations evolving in chemostats under nutrient limitation acquired beneficial mutations in hfq during slow growth (0.1 h⁻¹) but not in populations growing sixfold faster. Four identified hfq alleles from parallel populations were beneficial at 0.1 h⁻¹ and deleterious at 0.6 h⁻¹. The hfq mutations were beneficial, deleterious or neutral at an intermediate growth rate (0.5 h⁻¹) and one changed from beneficial to deleterious within a 36 min difference in doubling time. The benefit of hfq mutations was due to the greater transport of limiting nutrient, which diminished at higher growth rates. The deleterious effects of hfq mutations at 0.6 h⁻¹ were less clear, with decreased viability a contributing factor. The results demonstrate distinct pleiotropy characteristics in the alleles of the same gene, probably because the altered residues in Hfq affected the regulation of expression of different genes in distinct ways. In addition, these results point to a source of variation in experimental measurement of the selective advantage of a mutation; estimates of fitness need to consider variation in growth rate impacting on the magnitude of the benefit of mutations and on their fitness distributions.     

Genetic Characterization of Escherichia coli Populations from Host Sources of Fecal Pollution by Using DNA Fingerprinting

Escherichia coli isolates were obtained from common host sources of fecal pollution and characterized by using repetitive extragenic palindromic (REP) PCR fingerprinting. The genetic relationship of strains within each host group was assessed as was the relationship of strains among different host groups. Multiple isolates from a single host animal (gull, human, or dog) were found to be identical; however, in some of the animals, additional strains occurred at a lower frequency. REP PCR fingerprint patterns of isolates from sewage (n = 180), gulls (n = 133), and dairy cattle (n = 121) were diverse; within a host group, pairwise comparison similarity indices ranged from 98% to as low as 15%. A composite dendrogram of E. coli fingerprint patterns did not cluster the isolates into distinct host groups but rather produced numerous subclusters (approximately >80% similarity scores calculated with the cosine coefficient) that were nearly exclusive for a host group. Approximately 65% of the isolates analyzed were arranged into host-specific groups. Comparable results were obtained by using enterobacterial repetitive intergenic consensus PCR and pulsed-field gel electrophoresis (PFGE), where PFGE gave a higher differentiation of closely related strains than both PCR techniques. These results demonstrate that environmental studies with genetic comparisons to detect sources of E. coli contamination will require extensive isolation of strains to encompass E. coli strain diversity found in host sources of contamination. These findings will assist in the development of approaches to determine sources of fecal pollution, an effort important for protecting water resources and public health.     

The Molecular and Genetic Basis of Repeatable Coevolution between Escherichia coli and Bacteriophage T3 in a Laboratory Microcosm

The objective of this study was to determine the genomic changes that underlie coevolution between Escherichia coli B and bacteriophage T3 when grown together in a laboratory microcosm. We also sought to evaluate the repeatability of their evolution by studying replicate coevolution experiments inoculated with the same ancestral strains. We performed the coevolution experiments by growing Escherichia coli B and the lytic bacteriophage T3 in seven parallel continuous culture devices (chemostats) for 30 days. In each of the chemostats, we observed three rounds of coevolution. First, bacteria evolved resistance to infection by the ancestral phage. Then, a new phage type evolved that was capable of infecting the resistant bacteria as well as the sensitive bacterial ancestor. Finally, we observed second-order resistant bacteria evolve that were resistant to infection by both phage types. To identify the genetic changes underlying coevolution, we isolated first- and second-order resistant bacteria as well as a host-range mutant phage from each chemostat and sequenced their genomes. We found that first-order resistant bacteria consistently evolved resistance to phage via mutations in the gene, waaG, which codes for a glucosyltransferase required for assembly of the bacterial lipopolysaccharide (LPS). Phage also showed repeatable evolution, with each chemostat producing host-range mutant phage with mutations in the phage tail fiber gene T3p48 which binds to the bacterial LPS during adsorption. Two second-order resistant bacteria evolved via mutations in different genes involved in the phage interaction. Although a wide range of mutations occurred in the bacterial waaG gene, mutations in the phage tail fiber were restricted to a single codon, and several phage showed convergent evolution at the nucleotide level. These results are consistent with previous studies in other systems that have documented repeatable evolution in bacteria at the level of pathways or genes and repeatable evolution in viruses at the nucleotide level. Our data are also consistent with the expectation that adaptation via loss-of-function mutations is less constrained than adaptation via gain-of-function mutations.