IS186 insertion at a hot spot in the lon promoter as a basis for lon protease deficiency of Escherichia coli B: identification of a consensus target sequence for IS186 transposition

Evolution-assisted engineering of E. coli enables growth on formic acid at ambient CO2 via the Serine Threonine Cycle

Atmospheric CO2 poses a major threat to life on Earth by causing global warming and climate change. On the other hand, it can be considered as a resource that is scalable enough to establish a circular carbon economy. Accordingly, technologies to capture and convert CO2 into reduced one-carbon (C1) compounds (e.g. formic acid) are developing and improving fast. Driven by the idea of creating sustainable bioproduction platforms, natural and synthetic C1-utilization pathways are engineered into industrially relevant microbes. The realization of synthetic C1-assimilation cycles in living organisms is a promising but challenging endeavour. Here, we engineer the Serine Threonine Cycle, a synthetic C1-assimilation cycle in Escherichia coli to achieve growth on formic acid. Our stepwise engineering approach in tailored selection strains combined with adaptive laboratory evolution experiments enabled formatotrophic growth of the organism. Whole genome sequencing and reverse engineering allowed us to determine the key mutations linked to pathway activity. The Serine Threonine Cycle strains created in this work use formic acid as sole carbon and energy source and can grow at ambient CO2 cultivation conditions. This work sets an example for the engineering of complex C1-assimilation cycles in heterotrophic microbes.

Multilevel Systematic Optimization To Achieve Efficient Integrated Expression of Escherichia coli

Genomic integration of heterologous genes is the preferred approach in industrial fermentation-related strains due to the drawbacks associated with plasmid-mediated microbial fermentation, including additional growth burden, genetic instability, and antibiotic contamination. Synthetic biology and genome editing advancements have made gene integration convenient. Integrated expression is extensively used in the field of biomanufacturing and is anticipated to become the prevailing method for expressing recombinant proteins. Therefore, it is pivotal to strengthen the expression of exogenous genes at the genome level. Here, we systematically optimized the integrated expression system of Escherichia coli from 3 aspects. First, the integration site slmA with the highest expression activity was screened out of 18 sites in the ORI region of the E. coli BL21 (DE3) genome. Second, we characterized 16 endogenous promoters in E. coli and combined them with the T7 promoter. A constitutive promoter, Plpp-T7, exhibited significantly higher expression strength than the T7 promoter, achieving a 3.3-fold increase in expression levels. Finally, to further enhance the T7 expression system, we proceeded with overexpression of T7 RNA polymerase at the chassis cell level. The resulting constitutive efficient integrated expression system (CEIES_Ecoli) showed a 2-fold increase in GFP expression compared to the pET3b recombinant plasmid. Therefore, CEIES_Ecoli was applied to the integrated expression of nitrilase and hyaluronidase, achieving stable and efficient enzyme expression, with enzyme activities of 22.87 and 12,195 U·mL-1, respectively, comparable to plasmid levels. Overall, CEIES_Ecoli provides a stable and efficient method of gene expression without the need for antibiotics or inducers, making it a robust tool for synthetic biology, enzyme engineering, and related applications.

Suppression of the Escherichia coli rnpA49 conditionally lethal phenotype by different compensatory mutations

RNase P is an essential enzyme found across all domains of life that is responsible for the 5'-end maturation of precursor tRNAs. For decades, numerous studies have sought to elucidate the mechanisms and biochemistry governing RNase P function. However, much remains unknown about the regulation of RNase P expression, the turnover and degradation of the enzyme, and the mechanisms underlying the phenotypes and complementation of specific RNase P mutations, especially in the model bacterium, Escherichia coli In E. coli, the temperature-sensitive (ts) rnpA49 mutation in the protein subunit of RNase P has arguably been one of the most well-studied mutations for examining the enzyme's activity in vivo. Here, we report for the first time naturally occurring temperature-resistant suppressor mutations of E. coli strains carrying the rnpA49 allele. We find that rnpA49 strains can partially compensate the ts defect via gene amplifications of either RNase P subunit (rnpA49 or rnpB) or by the acquisition of loss-of-function mutations in Lon protease or RNase R. Our results agree with previous plasmid overexpression and gene deletion complementation studies, and importantly suggest the involvement of Lon protease in the degradation and/or regulatory pathway(s) of the mutant protein subunit of RNase P. This work offers novel insights into the behavior and complementation of the rnpA49 allele in vivo and provides direction for follow-up studies regarding RNase P regulation and turnover in E. coli.

Two classes of DNA gyrase inhibitors elicit distinct evolutionary trajectories toward resistance in gram-negative pathogens

Comprehensive knowledge of mechanisms driving the acquisition of antimicrobial resistance is essential for the development of new drugs with minimized resistibility. To gain this knowledge, we combine experimental evolution in a continuous culturing device, the morbidostat, with whole genome sequencing of evolving cultures followed by characterization of drug-resistant isolates. Here, this approach was used to assess evolutionary dynamics of resistance acquisition against DNA gyrase/topoisomerase TriBE inhibitor GP6 in Escherichia coli and Acinetobacter baumannii. The evolution of GP6 resistance in both species was driven by a combination of two classes of mutational events: (i) amino acid substitutions near the ATP-binding site of the GyrB subunit of the DNA gyrase target; and (ii) various mutations and genomic rearrangements leading to upregulation of efflux pumps, species-specific (AcrAB/TolC in E. coli and AdeIJK in A. baumannii) and shared by both species (MdtK). A comparison with the experimental evolution of resistance to ciprofloxacin (CIP), previously performed using the same workflow and strains, revealed fundamental differences between these two distinct classes of compounds. Most notable were non-overlapping spectra of target mutations and distinct evolutionary trajectories that, in the case of GP6, were dominated by upregulation of efflux machinery prior to (or even in lieu) of target modification. Most of the efflux-driven GP6-resistant isolates of both species displayed a robust cross-resistance to CIP, while CIP-resistant clones showed no appreciable increase in GP6-resistance.

Two Classes of DNA Gyrase Inhibitors Elicit Distinct Evolutionary Trajectories Toward Resistance in Gram-Negative Pathogens

Comprehensive knowledge of mechanisms driving the acquisition of antimicrobial resistance is essential for the development of new drugs with minimized resistibility. To gain this knowledge, we combine experimental evolution in a continuous culturing device, the morbidostat, with whole genome sequencing of evolving cultures followed by characterization of drug-resistant isolates. Here, this approach was used to assess evolutionary dynamics of resistance acquisition against DNA gyrase/topoisomerase TriBE inhibitor GP6 in Escherichia coli and Acinetobacter baumannii. The evolution of GP6 resistance in both species was driven by a combination of two classes of mutational events: (i) amino acid substitutions near the ATP-binding site of the GyrB subunit of the DNA gyrase target; and (ii) various mutations and genomic rearrangements leading to upregulation of efflux pumps, species-specific (AcrAB/TolC in E. coli and AdeIJK in A. baumannii) and shared by both species (MdtK). A comparison with the experimental evolution of resistance to ciprofloxacin (CIP), previously performed using the same workflow and strains, revealed fundamental differences between these two distinct classes of compounds. Most notable were non-overlapping spectra of target mutations and distinct evolutionary trajectories that, in the case of GP6, were dominated by upregulation of efflux machinery prior to (or even in lieu) of target modification. Most of efflux-driven GP6-resistant isolates of both species displayed a robust cross-resistance to CIP, while CIP-resistant clones showed no appreciable increase in GP6-resistance.

Effect of chromosomal integration on catalytic performance of a multi-component P450 system in Escherichia coli

Cytochromes P450 are useful biocatalysts in synthetic chemistry and important bio-bricks in synthetic biology. Almost all bacterial P450s require separate redox partners for their activity, which are often expressed in recombinant Escherichia coli using multiple plasmids. However, the application of CRISPR/Cas recombineering facilitated chromosomal integration of heterologous genes which enables more stable and tunable expression of multi-component P450 systems for whole-cell biotransformations. Herein, we compared three E. coli strains W3110, JM109, and BL21(DE3) harboring three heterologous genes encoding a P450 and two redox partners either on plasmids or after chromosomal integration in two genomic loci. Both loci proved to be reliable and comparable for the model regio- and stereoselective two-step oxidation of (S)-ketamine. Furthermore, the CRISPR/Cas-assisted integration of the T7 RNA polymerase gene enabled an easy extension of T7 expression strains. Higher titers of soluble active P450 were achieved in E. coli harboring a single chromosomal copy of the P450 gene compared to E. coli carrying a medium copy pET plasmid. In addition, improved expression of both redox partners after chromosomal integration resulted in up to 80% higher (S)-ketamine conversion and more than fourfold increase in total turnover numbers.

The diversity in antimicrobial resistance of MDR Enterobacteriaceae among Chinese broiler and laying farms and two mcr-1 positive plasmids revealed their resistance-transmission risk

This research aimed to investigate the microbial composition and diversity of antimicrobial resistance genes (ARGs) found in Chinese broiler and layer family poultry farms. We focused on the differences in resistance phenotypes and genotypes of multidrug-resistant Enterobacteriaceae (MDRE) isolated from the two farming environments and the existence and transmissibility of colistin resistance gene mcr-1. Metagenomic analysis showed that Firmicutes and Bacteroides were the dominant bacteria in broiler and layer farms. Many aminoglycoside and tetracycline resistance genes were accumulated in these environments, and their absolute abundance was higher in broiler than in layer farms. A total of 526 MDRE were isolated with a similar distribution in both farms. The results of the K-B test showed that the resistance rate to seven antimicrobials including polymyxin B and meropenem in broiler poultry farms was significantly higher than that in layer poultry farms (P ≤ 0.05). PCR screening results revealed that the detection rates of mcr-1, aph(3’)Ia, aadA2, bla_oxa–1, bla_CTX–M, fosB, qnrD, sul1, tetA, and catA1 in broiler source MDRE were significantly higher than those in layers (P ≤0.05). A chimeric plasmid p20432-mcr which carried the novel integron In1866 was isolated from broiler source MDRE. The high frequency of conjugation (10^–1 to 10^–3) and a wide range of hosts made p20432-mcr likely to play an essential role in the high detection rate of mcr-1, aph(3’)-Ia, and aadA2 in broiler farms. These findings will help optimize disinfection and improve antimicrobial-resistant bacteria surveillance programs in poultry farms, especially broilers.

Improving genomically recoded Escherichia coli to produce proteins containing non-canonical amino acids

A genomically recoded Escherichia coli strain that lacks all amber codons and release factor 1 (C321.∆A) enables efficient genetic encoding of chemically diverse non-canonical amino acids (ncAAs) into proteins. While C321.∆A has opened new opportunities in chemical and synthetic biology, this strain has not been optimized for protein production, limiting its utility in widespread industrial and academic applications. To address this limitation, the construction of a series of genomically recoded organisms that are optimized for cellular protein production is described. It is demonstrated that the functional deactivation of nucleases (e.g., rne, endA) and proteases (e.g., lon) increases production of wild-type superfolder green fluorescent protein (sfGFP) and sfGFP containing two ncAAs up to ≈5-fold. Additionally, a genomic IPTG-inducible T7 RNA polymerase (T7RNAP) cassette into these strains is introduced. Using an optimized platform, the ability to introduce two identical N6 -(propargyloxycarbonyl)-L -Lysine residues site specifically into sfGFP with a 17-fold improvement in production relative to the parent strain is demonstrated. The authors envision that their library of organisms will provide the community with multiple options for increased expression of proteins with new and diverse chemistries.

Identification of Evolutionary Trajectories Associated with Antimicrobial Resistance Using Microfluidics

In vitro experimental evolution of pathogens to antibiotics is commonly used for the identification of clinical biomarkers associated with antibiotic resistance. Microdroplet emulsions allow exquisite control of spatial structure, species complexity, and selection microenvironments for such studies. We investigated the use of monodisperse microdroplets in experimental evolution. Using Escherichia coli adaptation to doxycycline, we examined how changes in environmental conditions such as droplet size, starting lambda value, selection strength, and incubation method affected evolutionary outcomes. We also examined the extent to which emulsions could reveal potentially new evolutionary trajectories and dynamics associated with antimicrobial resistance. Interestingly, we identified both expected and unexpected evolutionary trajectories including large-scale chromosomal rearrangements and amplification that were not observed in suspension culture methods. As microdroplet emulsions are well-suited for automation and provide exceptional control of conditions, they can provide a high-throughput approach for biomarker identification as well as preclinical evaluation of lead compounds.

Improving the Intensity of Integrated Expression for Microbial Production

Chromosomal integration of exogenous genes is preferred for industrially related fermentation, as plasmid-mediated fermentation leads to extra metabolic burden and genetic instability. Moreover, with the development and advancement of genome engineering and gene editing technologies, inserting genes into chromosomes has become more convenient; integration expression is extensively utilized in microorganisms for industrial bioproduction and expected to become the trend of recombinant protein expression. However, in actual research and application, it is important to enhance the expression of heterologous genes at the host genome level. Herein, we summarized the basic principles and characteristics of genomic integration; furthermore, we highlighted strategies to improve the expression of chromosomal integration of genes and pathways in host strains from three aspects, including chassis cell optimization, regulation of expression elements in gene expression cassettes, optimization of gene dose level and integration sites on chromosomes. Moreover, we reviewed and summarized the relevant studies on the application of integrated expression in the exploration of gene function and the various types of industrial microorganism production. Consequently, this review would serve as a reference for the better application of integrated expression.

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.

Altered expression of a quality control protease in E. coli reshapes the in vivo mutational landscape of a model enzyme

Protein mutational landscapes are shaped by the cellular environment, but key factors and their quantitative effects are often unknown. Here we show that Lon, a quality control protease naturally absent in common E. coli expression strains, drastically reshapes the mutational landscape of the metabolic enzyme dihydrofolate reductase (DHFR). Selection under conditions that resolve highly active mutants reveals that 23.3% of all single point mutations in DHFR are advantageous in the absence of Lon, but advantageous mutations are largely suppressed when Lon is reintroduced. Protein stability measurements demonstrate extensive activity-stability tradeoffs for the advantageous mutants and provide a mechanistic explanation for Lon’s widespread impact. Our findings suggest possibilities for tuning mutational landscapes by modulating the cellular environment, with implications for protein design and combatting antibiotic resistance.

Evolution of satellite plasmids can prolong the maintenance of newly acquired accessory genes in bacteria

Transmissible plasmids spread genes encoding antibiotic resistance and other traits to new bacterial species. Here we report that laboratory populations of Escherichia coli with a newly acquired IncQ plasmid often evolve 'satellite plasmids' with deletions of accessory genes and genes required for plasmid replication. Satellite plasmids are molecular parasites: their presence reduces the copy number of the full-length plasmid on which they rely for their continued replication. Cells with satellite plasmids gain an immediate fitness advantage from reducing burdensome expression of accessory genes. Yet, they maintain copies of these genes and the complete plasmid, which potentially enables them to benefit from and transmit the traits they encode in the future. Evolution of satellite plasmids is transient. Cells that entirely lose accessory gene function or plasmid mobility dominate in the long run. Satellite plasmids also evolve in Snodgrassella alvi colonizing the honey bee gut, suggesting that this mechanism may broadly contribute to the importance of IncQ plasmids as agents of bacterial gene transfer in nature.

New tools for recombinant protein production in Escherichia coli: A 5-year update

The production of proteins in sufficient amounts is key for their study or use as biotherapeutic agents. Escherichia coli is the host of choice for recombinant protein production given its fast growth, easy manipulation, and cost-effectiveness. As such, its protein production capabilities are continuously being improved. Also, the associated tools (such as plasmids and cultivation conditions) are subject of ongoing research to optimize product yield. In this work, we review the latest advances in recombinant protein production in E. coli.

Genome rearrangements in Escherichia coli during de novo acquisition of resistance to a single antibiotic or two antibiotics successively

Background

The ability of bacteria to acquire resistance to antibiotics relies to a large extent on their capacity for genome modification. Prokaryotic genomes are highly plastic and can utilize horizontal gene transfer, point mutations, and gene deletions or amplifications to realize genome expansion and rearrangements. The contribution of point mutations to de novo acquisition of antibiotic resistance is well-established. In this study, the internal genome rearrangement of Escherichia coli during to de novo acquisition of antibiotic resistance was investigated using whole-genome sequencing.

Results

Cells were made resistant to one of the four antibiotics and subsequently to one of the three remaining. This way the initial genetic rearrangements could be documented together with the effects of an altered genetic background on subsequent development of resistance. A DNA fragment including ampC was amplified by a factor sometimes exceeding 100 as a result of exposure to amoxicillin. Excision of prophage e14 was observed in many samples with a double exposure history, but not in cells exposed to a single antibiotic, indicating that the activation of the SOS stress response alone, normally the trigger for excision, was not sufficient to cause excision of prophage e14. Partial deletion of clpS and clpA occurred in strains exposed to enrofloxacin and tetracycline. Other deletions were observed in some strains, but not in replicates with the exact same exposure history. Various insertion sequence transpositions correlated with exposure to specific antibiotics.

Conclusions

Many of the genome rearrangements have not been reported before to occur during resistance development. The observed correlation between genome rearrangements and specific antibiotic pressure, as well as their presence in independent replicates indicates that these events do not occur randomly. Taken together, the observed genome rearrangements illustrate the plasticity of the E. coli genome when exposed to antibiotic stress.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5353-y) contains supplementary material, which is available to authorized users.

An Integrated Proteomic Approach Uncovers Novel Substrates and Functions of the Lon Protease in Escherichia coli

Controlling the cellular abundance and proper function of proteins by proteolysis is a universal process in all living organisms. In Escherichia coli, the ATP-dependent Lon protease is crucial for protein quality control and regulatory processes. To understand how diverse substrates are selected and degraded, unbiased global approaches are needed. We employed a quantitative Super-SILAC (stable isotope labeling with amino acids in cell culture) mass spectrometry approach and compared the proteomes of a lon mutant and a strain producing the protease to discover Lon-dependent physiological functions. To identify Lon substrates, we took advantage of a Lon trapping variant, which is able to translocate substrates but unable to degrade them. Lon-associated proteins were identified by label-free LC-MS/MS. The combination of both approaches revealed a total of 14 novel Lon substrates. Besides the identification of known pathways affected by Lon, for example, the superoxide stress response, our cumulative data suggests previously unrecognized fundamental functions of Lon in sulfur assimilation, nucleotide biosynthesis, amino acid and central energy metabolism.

Tetracycline Antibiotics and Resistance

Tetracyclines possess many properties considered ideal for antibiotic drugs, including activity against Gram-positive and -negative pathogens, proven clinical safety, acceptable tolerability, and the availability of intravenous (IV) and oral formulations for most members of the class. As with all antibiotic classes, the antimicrobial activities of tetracyclines are subject to both class-specific and intrinsic antibiotic-resistance mechanisms. Since the discovery of the first tetracyclines more than 60 years ago, ongoing optimization of the core scaffold has produced tetracyclines in clinical use and development that are capable of thwarting many of these resistance mechanisms. New chemistry approaches have enabled the creation of synthetic derivatives with improved in vitro potency and in vivo efficacy, ensuring that the full potential of the class can be explored for use against current and emerging multidrug-resistant (MDR) pathogens, including carbapenem-resistant Enterobacteriaceae, MDR Acinetobacter species, and Pseudomonas aeruginosa.

Comparative genomics and experimental evolution of Escherichia coli BL21(DE3) strains reveal the landscape of toxicity escape from membrane protein overproduction

Achieving sufficient yields of proteins in their functional form represents the first bottleneck in contemporary bioscience and biotechnology. To accomplish successful overexpression of membrane proteins in a workhorse organism such as E. coli, defined and rational optimization strategies based on an understanding of the genetic background of the toxicity-escape mechanism are desirable. To this end, we sequenced the genomes of E. coli C41(DE3) and its derivative C43(DE3), which were developed for membrane protein production. Comparative analysis of their genomes with those of their ancestral strain E. coli BL21(DE3) revealed various genetic changes in both strains. A series of E. coli variants that are able to tolerate transformation with or overexpression of membrane proteins were generated by in vitro evolution. Targeted sequencing of the evolved strains revealed the mutational hotspots among the acquired genetic changes. By these combinatorial approaches, we found non-synonymous changes in the lac repressor gene of the lac operon as well as nucleotide substitutions in the lacUV5 promoter of the DE3 region, by which the toxic effect to the host caused by overexpression of membrane proteins could be relieved. A mutation in lacI was demonstrated to be crucial for conferring tolerance to membrane protein overexpression.

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.

Horizontal gene transfer (HGT) is central to bacterial evolution. The outcome of an HGT event (fixation in a population, elimination, or separation as a subdominant clone) depends not only on the availability of a new gene but crucially on the fitness cost or benefit of the genomic incorporation of the foreign gene and its expression in recipient bacteria. Here we studied the fitness landscape for inter-species chromosomal replacement of an essential protein, dihydrofolate reductase (DHFR) encoded by the folA gene, by its orthologs from other mesophilic bacteria. We purified and biochemically characterized 33 out of 35 orthologous DHFRs and found that most of them are stable and more catalytically active than E. coli DHFR. However, the inter-species replacement of DHFR caused significant fitness loss for most transgenic strains due to low abundance of orthologous DHFRs in E. coli cytoplasm. Laboratory evolution resulted in an increase in orthologous DHFR abundance leading to a dramatic fitness improvement. Genomic and proteomic analyses of “naive” and evolved strains suggest a new function of protein homeostasis to discriminate between “self” and “non-self” proteins, thus creating fitness barriers to HGT.

Identification of riboflavin: revealing different metabolic characteristics between Escherichia coli BL21(DE3) and MG1655

There are many physiological differences between Escherichia coli B and K-12 strains, owing to their different origins. Deeper insight into the metabolic and regulative mechanisms of these strains will inform improved usage of these industrial workhorses. In the present study, we observed that BL21 fermentation broth gradually turned yellow during cultivation. By spectral analysis and liquid chromatography-mass spectrometry identification, we confirmed for the first time that the yellow substance accumulated in the fermentation broth is riboflavin. Comparing the enzyme sequences involved in riboflavin metabolism between BL21 and MG1655, we identified a site mutation on the 115 residue of bifunctional riboflavin kinase/FMN adenylyltransferase (RibF) in BL21. This His115Leu mutation was found to reduce enzyme activity to 55% of that of MG1655, which is probably one reason for riboflavin accumulation in BL21. Quantitative PCR analysis showed that genes of the entire branch of the riboflavin and FAD biosynthesis pathways in BL21 were up-regulated. Several physiological and metabolic characteristics of BL21 and MG1655 were found to be different, and may also be related to the riboflavin accumulation.

Biocontainment of genetically modified organisms by synthetic protein design

Genetically modified organisms (GMOs) are increasingly deployed at large scales and in open environments. Genetic biocontainment strategies are needed to prevent unintended proliferation of GMOs in natural ecosystems. Existing biocontainment methods are insufficient because they impose evolutionary pressure on the organism to eject the safeguard by spontaneous mutagenesis or horizontal gene transfer, or because they can be circumvented by environmentally available compounds. Here we computationally redesign essential enzymes in the first organism possessing an altered genetic code (Escherichia coli strain C321.ΔA) to confer metabolic dependence on non-standard amino acids for survival. The resulting GMOs cannot metabolically bypass their biocontainment mechanisms using known environmental compounds, and they exhibit unprecedented resistance to evolutionary escape through mutagenesis and horizontal gene transfer. This work provides a foundation for safer GMOs that are isolated from natural ecosystems by a reliance on synthetic metabolites.

The genetic basis of Escherichia coli pathoadaptation to macrophages

Antagonistic interactions are likely important driving forces of the evolutionary process underlying bacterial genome complexity and diversity. We hypothesized that the ability of evolved bacteria to escape specific components of host innate immunity, such as phagocytosis and killing by macrophages (MΦ), is a critical trait relevant in the acquisition of bacterial virulence. Here, we used a combination of experimental evolution, phenotypic characterization, genome sequencing and mathematical modeling to address how fast, and through how many adaptive steps, a commensal Escherichia coli (E. coli) acquire this virulence trait. We show that when maintained in vitro under the selective pressure of host MΦ commensal E. coli can evolve, in less than 500 generations, virulent clones that escape phagocytosis and MΦ killing in vitro, while increasing their pathogenicity in vivo, as assessed in mice. This pathoadaptive process is driven by a mechanism involving the insertion of a single transposable element into the promoter region of the E. coli yrfF gene. Moreover, transposition of the IS186 element into the promoter of Lon gene, encoding an ATP-dependent serine protease, is likely to accelerate this pathoadaptive process. Competition between clones carrying distinct beneficial mutations dominates the dynamics of the pathoadaptive process, as suggested from a mathematical model, which reproduces the observed experimental dynamics of E. coli evolution towards virulence. In conclusion, we reveal a molecular mechanism explaining how a specific component of host innate immunity can modulate microbial evolution towards pathogenicity.

Comparative multi-omics systems analysis of Escherichia coli strains B and K-12

BACKGROUND

Elucidation of a genotype-phenotype relationship is critical to understand an organism at the whole-system level. Here, we demonstrate that comparative analyses of multi-omics data combined with a computational modeling approach provide a framework for elucidating the phenotypic characteristics of organisms whose genomes are sequenced.

RESULTS

We present a comprehensive analysis of genome-wide measurements incorporating multifaceted holistic data - genome, transcriptome, proteome, and phenome - to determine the differences between Escherichia coli B and K-12 strains. A genome-scale metabolic network of E. coli B was reconstructed and used to identify genetic bases of the phenotypes unique to B compared with K-12 through in silico complementation testing. This systems analysis revealed that E. coli B is well-suited for production of recombinant proteins due to a greater capacity for amino acid biosynthesis, fewer proteases, and lack of flagella. Furthermore, E. coli B has an additional type II secretion system and a different cell wall and outer membrane composition predicted to be more favorable for protein secretion. In contrast, E. coli K-12 showed a higher expression of heat shock genes and was less susceptible to certain stress conditions.

CONCLUSIONS

This integrative systems approach provides a high-resolution system-wide view and insights into why two closely related strains of E. coli, B and K-12, manifest distinct phenotypes. Therefore, systematic understanding of cellular physiology and metabolism of the strains is essential not only to determine culture conditions but also to design recombinant hosts.

The quorum-sensing protein TraR of Agrobacterium tumefaciens is susceptible to intrinsic and TraM-mediated proteolytic instability

TraR of Agrobacterium tumefaciens is a LuxR-type transcription factor that regulates genes required for replication and conjugation of the tumour-inducing plasmid. TraR binds the pheromone 3-oxo-octanoylhomoserine lactone (OOHL) and requires this molecule for folding into a protease-resistant, soluble conformation. Even after binding to OOHL, TraR is degraded at readily detectable rates. Here we show that the N-terminal domain of TraR, which binds OOHL, is more resistant to degradation than the full length protein, suggesting that sites on the C-terminal DNA binding domain [TraR(170-234)] enhance protein turnover. A fusion between GFP and TraR(170-234) was poorly fluorescent, and truncations of this fusion protein allowed us to identify residues in this domain that contribute to protein degradation. TraR activity was previously shown to be inhibited by the antiactivator TraM. These proteins form 2:2 complexes that fail to bind DNA sequences. Here we show that TraM sharply decreased the accumulation of TraR in whole cells, indicating that TraM facilitates proteolysis of TraR. The TraM component of these complexes is spared from proteolysis, and could therefore act catalytically.

Using genomic sequencing for classical genetics in E. coli K12

We here develop computational methods to facilitate use of 454 whole genome shotgun sequencing to identify mutations in Escherichia coli K12. We had Roche sequence eight related strains derived as spontaneous mutants in a background without a whole genome sequence. They provided difference tables based on assembling each genome to reference strain E. coli MG1655 (NC_000913). Due to the evolutionary distance to MG1655, these contained a large number of both false negatives and positives. By manual analysis of the dataset, we detected all the known mutations (24 at nine locations) and identified and genetically confirmed new mutations necessary and sufficient for the phenotypes we had selected in four strains. We then had Roche assemble contigs de novo, which we further assembled to full-length pseudomolecules based on synteny with MG1655. This hybrid method facilitated detection of insertion mutations and allowed annotation from MG1655. After removing one genome with less than the optimal 20- to 30-fold sequence coverage, we identified 544 putative polymorphisms that included all of the known and selected mutations apart from insertions. Finally, we detected seven new mutations in a total of only 41 candidates by comparing single genomes to composite data for the remaining six and using a ranking system to penalize homopolymer sequencing and misassembly errors. An additional benefit of the analysis is a table of differences between MG1655 and a physiologically robust E. coli wild-type strain NCM3722. Both projects were greatly facilitated by use of comparative genomics tools in the CoGe software package (http://genomevolution.org/).

Construction and performance of heterologous polyketide-producing K-12- and B-derived Escherichia coli

AIMS

Escherichia coli has emerged as a viable heterologous host for the production of complex, polyketide natural compounds. In this study, polyketide biosynthesis was compared between different E. coli strains for the purpose of better understanding and improving heterologous production.

METHODS AND RESULTS

Both B and K-12 E. coli strains were genetically modified to support heterologous polyketide biosynthesis [specifically, 6-deoxyerythronolide B (6dEB)]. Polyketide production was analysed using a helper plasmid designed to overcome rare codon usage within E. coli. Each strain was analysed for recombinant protein production, precursor consumption, by-product production, and 6dEB biosynthesis. Of the strains tested for biosynthesis, 6dEB production was greatest for E. coli B strains. When comparing biosynthetic improvements as a function of mRNA stability vs codon bias, increased 6dEB titres were observed when additional rare codon tRNA molecules were provided.

CONCLUSIONS

Escherichia coli B strains and the use of tRNA supplementation led to improved 6dEB polyketide titres.

SIGNIFICANCE AND IMPACT OF THE STUDY

Given the medicinal potential and growing field of polyketide heterologous biosynthesis, the current study provides insight into host-specific genetic backgrounds and gene expression parameters aiding polyketide production through E. coli.

The Rut pathway for pyrimidine degradation: novel chemistry and toxicity problems

The Rut pathway is composed of seven proteins, all of which are required by Escherichia coli K-12 to grow on uracil as the sole nitrogen source. The RutA and RutB proteins are central: no spontaneous suppressors arise in strains lacking them. RutA works in conjunction with a flavin reductase (RutF or a substitute) to catalyze a novel reaction. It directly cleaves the uracil ring between N-3 and C-4 to yield ureidoacrylate, as established by both nuclear magnetic resonance (NMR) spectroscopy and mass spectrometry. Although ureidoacrylate appears to arise by hydrolysis, the requirements for the reaction and the incorporation of (18)O at C-4 from molecular oxygen indicate otherwise. Mass spectrometry revealed the presence of a small amount of product with the mass of ureidoacrylate peracid in reaction mixtures, and we infer that this is the direct product of RutA. In vitro RutB cleaves ureidoacrylate hydrolytically to release 2 mol of ammonium, malonic semialdehyde, and carbon dioxide. Presumably the direct products are aminoacrylate and carbamate, both of which hydrolyze spontaneously. Together with bioinformatic predictions and published crystal structures, genetic and physiological studies allow us to predict functions for RutC, -D, and -E. In vivo we postulate that RutB hydrolyzes the peracid of ureidoacrylate to yield the peracid of aminoacrylate. We speculate that RutC reduces aminoacrylate peracid to aminoacrylate and RutD increases the rate of spontaneous hydrolysis of aminoacrylate. The function of RutE appears to be the same as that of YdfG, which reduces malonic semialdehyde to 3-hydroxypropionic acid. RutG appears to be a uracil transporter.

Tuning a genetic switch: experimental evolution and natural variation of prophage induction

Genetic switches allow organisms to modulate their phenotype in response to environmental changes. Understanding the evolutionary processes by which switches are tuned is central to understanding how phenotypic variation is realized. Prophage induction by phage lambda is the classic example of a genetic switch and allows lambda to move between two different modes of transmission: as a lysogen it reproduces vertically as a component of the host genome; as a free phage it reproduces horizontally by infectious epidemic spread. We show that the lambda switch can respond rapidly to selection for alteration in sensitivity and threshold. Sequencing of candidate genes in the genetic circuitry underlying the switch revealed mutations of likely adaptive significance in some, but not all candidates, suggesting that the core genetic circuitry plays a limited role in the fine-tuning of the switch in vivo. The relative ease with which the switch could be tuned by selection was further indicated by extensive variation in sensitivity and threshold of its response function among wild lambdoid phages. Together, our findings emphasize the adaptive significance of a finely tuned switch and draw attention to the selective factors shaping prophage induction in natural phage populations.

Genome sequences of Escherichia coli B strains REL606 and BL21(DE3)

Escherichia coli K-12 and B have been the subjects of classical experiments from which much of our understanding of molecular genetics has emerged. We present here complete genome sequences of two E. coli B strains, REL606, used in a long-term evolution experiment, and BL21(DE3), widely used to express recombinant proteins. The two genomes differ in length by 72,304 bp and have 426 single base pair differences, a seemingly large difference for laboratory strains having a common ancestor within the last 67 years. Transpositions by IS1 and IS150 have occurred in both lineages. Integration of the DE3 prophage in BL21(DE3) apparently displaced a defective prophage in the lambda attachment site of B. As might have been anticipated from the many genetic and biochemical experiments comparing B and K-12 over the years, the B genomes are similar in size and organization to the genome of E. coli K-12 MG1655 and have >99% sequence identity over approximately 92% of their genomes. E. coli B and K-12 differ considerably in distribution of IS elements and in location and composition of larger mobile elements. An unexpected difference is the absence of a large cluster of flagella genes in B, due to a 41 kbp IS1-mediated deletion. Gene clusters that specify the LPS core, O antigen, and restriction enzymes differ substantially, presumably because of horizontal transfer. Comparative analysis of 32 independently isolated E. coli and Shigella genomes, both commensals and pathogenic strains, identifies a minimal set of genes in common plus many strain-specific genes that constitute a large E. coli pan-genome.

Understanding the differences between genome sequences of Escherichia coli B strains REL606 and BL21(DE3) and comparison of the E. coli B and K-12 genomes

Each difference between the genome sequences of Escherichia coli B strains REL606 and BL21(DE3) can be interpreted in light of known laboratory manipulations plus a gene conversion between ribosomal RNA operons. Two treatments with 1-methyl-3-nitro-1-nitrosoguanidine in the REL606 lineage produced at least 93 single-base-pair mutations ( approximately 90% GC-to-AT transitions) and 3 single-base-pair GC deletions. Two UV treatments in the BL21(DE3) lineage produced only 4 single-base-pair mutations but 16 large deletions. P1 transductions from K-12 into the two B lineages produced 317 single-base-pair differences and 9 insertions or deletions, reflecting differences between B DNA in BL21(DE3) and integrated restriction fragments of K-12 DNA inherited by REL606. Two sites showed selective enrichment of spontaneous mutations. No unselected spontaneous single-base-pair mutations were evident. The genome sequences revealed that a progenitor of REL606 had been misidentified, explaining initially perplexing differences. Limited sequencing of other B strains defined characteristic properties of B and allowed assembly of the inferred genome of the ancestral B of Delbrück and Luria. Comparison of the B and K-12 genomes shows that more than half of the 3793 proteins of their basic genomes are predicted to be identical, although approximately 310 appear to be functional in either B or K-12 but not in both. The ancestral basic genome appears to have had approximately 4039 coding sequences occupying approximately 4.0 Mbp. Repeated horizontal transfer from diverged Escherichia coli genomes and homologous recombination may explain the observed variable distribution of single-base-pair differences. Fifteen sites are occupied by phage-related elements, but only six by comparable elements at the same site. More than 50 sites are occupied by IS elements in both B and K, 16 in common, and likely founding IS elements are identified. A signature of widespread cryptic phage P4-type mobile elements was identified. Complex deletions (dense clusters of small deletions and substitutions) apparently removed nonessential genes from approximately 30 sites in the basic genomes.

Diversify or Die: Generation of Diversity in Response to Stress

When challenged with unfavorable conditions, microorganisms can develop a stress response that allows them to adapt to or survive in the new environment. A common feature of the numerous specific stress response pathways that have been described in a wide range of bacteria is that they are energy demanding and therefore often transient. In addition, stress responses may come too late or be insufficient to protect the cell or the population against very sudden or severe stresses. However, it seems that microorganisms can also enhance their chances of survival under stress by increasing the generation of diversity at the population level. This can be achieved either by creating genetic diversity by a variety of mechanisms involving for example constitutive or transient mutators and contingency loci, or by revealing phenotypic diversity that remained dormant due to a mechanism called genetic buffering. This review gives an overview of these emerging diversity-generating mechanisms, which seem to play an important role in the ability of microbial populations to overcome stress challenges.

Increased genome instability in Escherichia coli lon mutants: relation to emergence of multiple-antibiotic-resistant (Mar) mutants caused by insertion sequence elements and large tandem genomic amplifications

Thirteen spontaneous multiple-antibiotic-resistant (Mar) mutants of Escherichia coli AG100 were isolated on Luria-Bertani (LB) agar in the presence of tetracycline (4 microg/ml). The phenotype was linked to insertion sequence (IS) insertions in marR or acrR or unstable large tandem genomic amplifications which included acrAB and which were bordered by IS3 or IS5 sequences. Five different lon mutations, not related to the Mar phenotype, were also found in 12 of the 13 mutants. Under specific selective conditions, most drug-resistant mutants appearing late on the selective plates evolved from a subpopulation of AG100 with lon mutations. That the lon locus was involved in the evolution to low levels of multidrug resistance was supported by the following findings: (i) AG100 grown in LB broth had an important spontaneous subpopulation (about 3.7x10(-4)) of lon::IS186 mutants, (ii) new lon mutants appeared during the selection on antibiotic-containing agar plates, (iii) lon mutants could slowly grow in the presence of low amounts (about 2x MIC of the wild type) of chloramphenicol or tetracycline, and (iv) a lon mutation conferred a mutator phenotype which increased IS transposition and genome rearrangements. The association between lon mutations and mutations causing the Mar phenotype was dependent on the medium (LB versus MacConkey medium) and the antibiotic used for the selection. A previously reported unstable amplifiable high-level resistance observed after the prolonged growth of Mar mutants in a low concentration of tetracycline or chloramphenicol can be explained by genomic amplification.

Role for tandem duplication and lon protease in AcrAB-TolC- dependent multiple antibiotic resistance (Mar) in an Escherichia coli mutant without mutations in marRAB or acrRAB

A spontaneous mutant (M113) of Escherichia coli AG100 with an unstable multiple antibiotic resistance (Mar) phenotype was isolated in the presence of tetracycline. Two mutations were found: an insertion in the promoter of lon (lon3::IS186) that occurred first and a subsequent large tandem duplication, dupIS186, bearing the genes acrAB and extending from the lon3::IS186 to another IS186 present 149 kb away from lon. The decreased amount of Lon protease increased the amount of MarA by stabilization of the basal quantities of MarA produced, which in turn increased the amount of multidrug effux pump AcrAB-TolC. However, in a mutant carrying only a lon mutation, the overproduced pump mediated little, if any, increased multidrug resistance, indicating that the Lon protease was required for the function of the pump. This requirement was only partial since resistance was mediated when amounts of AcrAB in a lon mutant were further increased by a second mutation. In M113, amplification of acrAB on the duplication led to increased amounts of AcrAB and multidrug resistance. Spontaneous gene duplication represents a new mechanism for mediating multidrug resistance in E. coli through AcrAB-TolC.

Separate-component-stabilization system for protein and DNA production without the use of antibiotics

Plasmid instability is a significant concern in the industrial utilization of microorganisms for protein or DNA production. Here we report on the development of a new and highly effective stabilization system based on the use of the ccd antidote/poison genes. For the first time, we separated the antidote gene from the poison gene: localizing the former in the plasmid and integrating the latter in the bacterial chromosome. We show that this separate-component-stabilization (SCS) strategy: (i) allows for perfect stabilization without the use of antibiotics; (ii) increases three to five times the recombinant protein production levels; and (iii) does not require any specific modification of the protein production process or culture medium. We illustrate that point by using the classical T7 promotor (i.e., used in most expression systems). Finally, we demonstrate that the SCS system increases by five the yield in DNA production, a result especially important for the design and production of gene therapy constructs void of any antibiotic resistance gene.

SulA-dependent hypersensitivity to high pressure and hyperfilamentation after high-pressure treatment of Escherichia coli lon mutants

High-pressure treatment (>100 MPa) is known to induce several heat shock proteins as well as an SOS response in Escherichia coli. In the current work, we have investigated properties with respect to high-pressure treatment of mutants-deficient in Lon, a pressure-induced ATP-dependent protease that belongs to the heat shock regulon but that also has a link to the SOS regulon. We report that lon mutants show increased pressure sensitivity and exhibit hyperfilamentation during growth after high-pressure treatment. Both phenotypes could be entirely attributed to the action of the SOS protein SulA, a potent inhibitor of the cell division ring protein FtsZ and a specific target of the Lon protease, since they were suppressed by knock-out of SulA. Introduction of the lexA1 allele, which effectively blocks the entire SOS response, also suppressed the high pressure hypersensitivity of lon mutants, but not their UV hypersensitivity. These results indicate the existence of a SulA-dependent pathway of high-pressure-induced cell filamentation, and suggest involvement of the SOS response, and particularly of SulA, in high-pressure-mediated cell death in E. coli strains which are compromised in Lon function.

Directed evolution of single-chain Fv for cytoplasmic expression using the beta-galactosidase complementation assay results in proteins highly susceptible to protease degradation and aggregation

Background

Antibody fragments are molecules widely used for diagnosis and therapy. A large amount of protein is frequently required for such applications. New approaches using folding reporter enzymes have recently been proposed to increase soluble expression of foreign proteins in Escherichia coli. To date, these methods have only been used to screen for proteins with better folding properties but have never been used to select from a large library of mutants. In this paper we apply one of these methods to select mutations that increase the soluble expression of two antibody fragments in the cytoplasm of E. coli.

Results

We used the β-galactosidase α-complementation system to monitor and evolve two antibody fragments for high expression levels in E. coli cytoplasm. After four rounds of mutagenesis and selection from large library repertoires (>10⁷ clones), clones exhibiting high levels of β-galactosidase activity were isolated. These clones expressed a higher amount of soluble fusion protein than the wild type in the cytoplasm, particularly in a strain deficient in the cytoplasmic Lon protease. The increase in the soluble expression level of the unfused scFv was, however, much less pronounced, and the unfused proteins proved to be more aggregation prone than the wild type. In addition, the soluble expression levels were not correlated with the β-galactosidase activity present in the cells.

Conclusion

This is the first report of a selection for soluble protein expression using a fusion reporter method. Contrary to anticipated results, high enzymatic activity did not correlate with the soluble protein expression level. This was presumably due to free α-peptide released from the protein fusion by the host proteases. This means that the α-complementation assay does not sense the fusion expression level, as hypothesized, but rather the amount of free released α-peptide. Thus, the system does not select, in our case, for higher soluble protein expression level but rather for higher protease susceptibility of the fusion protein.

Detection and analysis of transpositionally active head-to-tail dimers in three additional Escherichia coli IS elements

This study demonstrates that Escherichia coli insertion elements IS3, IS150 and IS186 are able to form transpositionally active head-to-tail dimers which show similar structure and transpositional activity to the dimers of IS2, IS21 and IS30. These structures arise by joining of the left and right inverted repeats (IRs) of two elements. The resulting junction includes a spacer region (SR) of a few base pairs derived from the flanking sequence of one of the reacting IRs. Head-to-tail dimers of IS3, IS150 and IS186 are unstable due to their transpositional activity. They can be resolved in two ways that seem to form a general rule for those elements reported to form dimers. One way is a site-specific process (dimer dissolution) which is accompanied by the loss of one IS copy along with the SR. The other is 'classical' transposition where the joined ends integrate into the target DNA. In intramolecular transposition this often gives rise to deletion formation, whereas in intermolecular transposition it gives rise to replicon fusion. The results presented for IS3, IS150 and IS186 are in accordance with the IS dimer model, which is in turn consistent with models based on covalently closed minicircles.