Multicopy suppression underpins metabolic evolvability

Unveiling the critical role of K+ for xanthorhodopsin expression in E. coli

Xanthorhodopsin (XR), a retinal-binding 7-transmembrane protein isolated from the eubacterium Salinibacter ruber, utilizes two chromophores (retinal and salinixanthin (SAL)) as an outward proton pump and energy-donating carotenoid. However, research on XR has been impeded owing to limitations in achieving heterogeneous expression of stable forms and high production levels of both wild-type and mutants. We successfully expressed wild-type and mutant XRs in Escherichia coli in the presence of K+. Achieving XR expression requires significant K+ and a low inducer concentration. In particular, we highlight the significance of Ser-159 in helix E located near Gly-156 (a carotenoid-binding position) as a critical site for XR expression. Our findings indicate that replacing Ser-159 with a smaller amino acid, alanine, can enhance XR expression in a manner comparable to K+, implying that Ser-159 poses a steric hindrance for pigment formation in XR. In the presence of K+, the proton pumping and photocycle of the wild-type and mutants were characterized and compared; the wild-type result suggests similar properties to the first reported XR isolation from the S. ruber membrane fraction. We propose that the K+ gradient across the cell membrane of S. ruber serves to uphold the membrane potential of the organism and plays a role in the expression of proteins, such as XR, as demonstrated in our study. Our findings deepen the understanding of adaptive protein expression, particularly in halophilic organisms. We highlight salt selection as a promising strategy for improving protein yield and functionality.

Genomic factors shape carbon and nitrogen metabolic niche breadth across Saccharomycotina yeasts

Organisms exhibit extensive variation in ecological niche breadth, from very narrow (specialists) to very broad (generalists). Two general paradigms have been proposed to explain this variation: (i) trade-offs between performance efficiency and breadth and (ii) the joint influence of extrinsic (environmental) and intrinsic (genomic) factors. We assembled genomic, metabolic, and ecological data from nearly all known species of the ancient fungal subphylum Saccharomycotina (1154 yeast strains from 1051 species), grown in 24 different environmental conditions, to examine niche breadth evolution. We found that large differences in the breadth of carbon utilization traits between yeasts stem from intrinsic differences in genes encoding specific metabolic pathways, but we found limited evidence for trade-offs. These comprehensive data argue that intrinsic factors shape niche breadth variation in microbes.

Expanding the flexibility of base editing for high-throughput genetic screens in bacteria

Genome-wide screens have become powerful tools for elucidating genotype-to-phenotype relationships in bacteria. Of the varying techniques to achieve knockout and knockdown, CRISPR base editors are emerging as promising options. However, the limited number of available, efficient target sites hampers their use for high-throughput screening. Here, we make multiple advances to enable flexible base editing as part of high-throughput genetic screening in bacteria. We first co-opt the Streptococcus canis Cas9 that exhibits more flexible protospacer-adjacent motif recognition than the traditional Streptococcus pyogenes Cas9. We then expand beyond introducing premature stop codons by mutating start codons. Next, we derive guide design rules by applying machine learning to an essentiality screen conducted in Escherichia coli. Finally, we rescue poorly edited sites by combining base editing with Cas9-induced cleavage of unedited cells, thereby enriching for intended edits. The efficiency of this dual system was validated through a conditional essentiality screen based on growth in minimal media. Overall, expanding the scope of genome-wide knockout screens with base editors could further facilitate the investigation of new gene functions and interactions in bacteria.

A journey into the regulatory secrets of the de novo purine nucleotide biosynthesis

De novo purine nucleotide biosynthesis (DNPNB) consists of sequential reactions that are majorly conserved in living organisms. Several regulation events take place to maintain physiological concentrations of adenylate and guanylate nucleotides in cells and to fine-tune the production of purine nucleotides in response to changing cellular demands. Recent years have seen a renewed interest in the DNPNB enzymes, with some being highlighted as promising targets for therapeutic molecules. Herein, a review of two newly revealed modes of regulation of the DNPNB pathway has been carried out: i) the unprecedent allosteric regulation of one of the limiting enzymes of the pathway named inosine 5'-monophosphate dehydrogenase (IMPDH), and ii) the supramolecular assembly of DNPNB enzymes. Moreover, recent advances that revealed the therapeutic potential of DNPNB enzymes in bacteria could open the road for the pharmacological development of novel antibiotics.

Targeted Transcriptional Activation Using a CRISPR-Associated Transposon System

Synthetic perturbation of gene expression is central to our ability to reliably uncover genotype-phenotype relationships in microbes. Here, we present a novel transcription activation strategy that uses the Vibrio cholerae CRISPR-Associated Transposon (CAST) system to selectively insert promoter elements upstream of genes of interest. Through this strategy, we show robust activation of both recombinant and endogenous genes across the Escherichia coli chromosome. We then demonstrate the precise tuning of expression levels by exchanging the promoter elements being inserted. Finally, we demonstrate that CAST activation can be used to synthetically induce ampicillin-resistant phenotypes in E. coli.

The EcoCyc Database (2023)

EcoCyc is a bioinformatics database available online at EcoCyc.org that describes the genome and the biochemical machinery of Escherichia coli K-12 MG1655. The long-term goal of the project is to describe the complete molecular catalog of the E. coli cell, as well as the functions of each of its molecular parts, to facilitate a system-level understanding of E. coli. EcoCyc is an electronic reference source for E. coli biologists and for biologists who work with related microorganisms. The database includes information pages on each E. coli gene product, metabolite, reaction, operon, and metabolic pathway. The database also includes information on the regulation of gene expression, E. coli gene essentiality, and nutrient conditions that do or do not support the growth of E. coli. The website and downloadable software contain tools for the analysis of high-throughput data sets. In addition, a steady-state metabolic flux model is generated from each new version of EcoCyc and can be executed online. The model can predict metabolic flux rates, nutrient uptake rates, and growth rates for different gene knockouts and nutrient conditions. Data generated from a whole-cell model that is parameterized from the latest data on EcoCyc are also available. This review outlines the data content of EcoCyc and of the procedures by which this content is generated.

Mapping temperature-sensitive mutations at a genome scale to engineer growth switches in Escherichia coli

Temperature-sensitive (TS) mutants are a unique tool to perturb and engineer cellular systems. Here, we constructed a CRISPR library with 15,120 Escherichia coli mutants, each with a single amino acid change in one of 346 essential proteins. 1,269 of these mutants showed temperature-sensitive growth in a time-resolved competition assay. We reconstructed 94 TS mutants and measured their metabolism under growth arrest at 42°C using metabolomics. Metabolome changes were strong and mutant-specific, showing that metabolism of nongrowing E. coli is perturbation-dependent. For example, 24 TS mutants of metabolic enzymes overproduced the direct substrate metabolite due to a bottleneck in their associated pathway. A strain with TS homoserine kinase (ThrBF267D ) produced homoserine for 24 h, and production was tunable by temperature. Finally, we used a TS subunit of DNA polymerase III (DnaXL289Q ) to decouple growth from arginine overproduction in engineered E. coli. These results provide a strategy to identify TS mutants en masse and demonstrate their large potential to produce bacterial metabolites with nongrowing cells.

Genomic and ecological factors shaping specialism and generalism across an entire subphylum

Organisms exhibit extensive variation in ecological niche breadth, from very narrow (specialists) to very broad (generalists). Paradigms proposed to explain this variation either invoke trade-offs between performance efficiency and breadth or underlying intrinsic or extrinsic factors. We assembled genomic (1,154 yeast strains from 1,049 species), metabolic (quantitative measures of growth of 843 species in 24 conditions), and ecological (environmental ontology of 1,088 species) data from nearly all known species of the ancient fungal subphylum Saccharomycotina to examine niche breadth evolution. We found large interspecific differences in carbon breadth stem from intrinsic differences in genes encoding specific metabolic pathways but no evidence of trade-offs and a limited role of extrinsic ecological factors. These comprehensive data argue that intrinsic factors driving microbial niche breadth variation.

Cell Survival Enabled by Leakage of a Labile Metabolic Intermediate

Many metabolites are generated in one step of a biochemical pathway and consumed in a subsequent step. Such metabolic intermediates are often reactive molecules which, if allowed to freely diffuse in the intracellular milieu, could lead to undesirable side reactions and even become toxic to the cell. Therefore, metabolic intermediates are often protected as protein-bound species and directly transferred between enzyme active sites in multi-functional enzymes, multi-enzyme complexes, and metabolons. Sequestration of reactive metabolic intermediates thus contributes to metabolic efficiency. It is not known, however, whether this evolutionary adaptation can be relaxed in response to challenges to organismal survival. Here, we report evolutionary repair experiments on Escherichia coli cells in which an enzyme crucial for the biosynthesis of proline has been deleted. The deletion makes cells unable to grow in a culture medium lacking proline. Remarkably, however, cell growth is efficiently restored by many single mutations (12 at least) in the gene of glutamine synthetase. The mutations cause the leakage to the intracellular milieu of a highly reactive phosphorylated intermediate common to the biosynthetic pathways of glutamine and proline. This intermediate is generally assumed to exist only as a protein-bound species. Nevertheless, its diffusion upon mutation-induced leakage enables a new route to proline biosynthesis. Our results support that leakage of sequestered metabolic intermediates can readily occur and contribute to organismal adaptation in some scenarios. Enhanced availability of reactive molecules may enable the generation of new biochemical pathways and the potential of mutation-induced leakage in metabolic engineering is noted.

Pre-epidemic evolution of the MRSA USA300 clade and a molecular key for classification

Introduction

USA300 has remained the dominant community and healthcare associated methicillin-resistant Staphylococcus aureus (MRSA) clone in the United States and in northern South America for at least the past 20 years. In this time, it has experienced epidemic spread in both of these locations. However, its pre-epidemic evolutionary history and origins are incompletely understood. Large sequencing databases, such as NCBI, PATRIC, and Staphopia, contain clues to the early evolution of USA300 in the form of sequenced genomes of USA300 isolates that are representative of lineages that diverged prior to the establishment of the South American epidemic (SAE) clade and North American epidemic (NAE) clade. In addition, historical isolates collected prior to the emergence of epidemics can help reconstruct early events in the history of this lineage.

Methods

Here, we take advantage of the accrued, publicly available data, as well as two newly sequenced pre-epidemic historical isolates from 1996, and a very early diverging ACME-negative NAE genome, to understand the pre-epidemic evolution of USA300. We use database mining techniques to emphasize genomes similar to pre-epidemic isolates, with the goal of reconstructing the early molecular evolution of the USA300 lineage.

Results

Phylogenetic analysis with these genomes confirms that the NAE and SAE USA300 lineages diverged from a most recent common ancestor around 1970 with high confidence, and it also pinpoints the independent acquisition events of the of the ACME and COMER loci with greater precision than in previous studies. We provide evidence for a North American origin of the USA300 lineage and identify multiple introductions of USA300 into South and North America. Notably, we describe a third major USA300 clade (the pre-epidemic branching clade; PEB1) consisting of both MSSA and MRSA isolates circulating around the world that diverged from the USA300 lineage prior to the establishment of the South and North American epidemics. We present a detailed analysis of specific sequence characteristics of each of the major clades, and present diagnostic positions that can be used to classify new genomes.

How to Recruit a Promiscuous Enzyme to Serve a New Function

Promiscuous enzymes can be recruited to serve new functions when a genetic or environmental change makes catalysis of a novel reaction important for fitness or even survival. Subsequently, gene duplication and divergence can lead to evolution of an efficient and specialized new enzyme. Every organism likely has thousands of promiscuous enzyme activities that provide a vast reservoir of catalytic potential. However, much of this potential may not be accessible. We compiled kinetic parameters for promiscuous reactions catalyzed by 108 enzymes. The median value of k_cat/K_M is a very modest 31 M^-1 s^-1. Based upon the fluxes through metabolic pathways in E. coli, we estimate that many, if not most, promiscuous activities are too inefficient to impact fitness. However, mutations can elevate the level of an insufficient promiscuous activity by increasing enzyme expression, improving k_cat/K_M, or altering concentrations of the promiscuous and native substrates and allosteric regulators. Particularly in large bacterial populations, stochastic mutations may provide a viable pathway for recruitment of even inefficient promiscuous activities.

Heliorhodopsin binds and regulates glutamine synthetase activity

Photoreceptors are light-sensitive proteins found in various organisms that respond to light and relay signals into the cells. Heliorhodopsin, a retinal-binding membrane protein, has been recently discovered, however its function remains unknown. Herein, we investigated the relationship between Actinobacteria bacterium IMCC26103 heliorhodopsin (AbHeR) and an adjacent glutamine synthetase (AbGS) in the same operon. We demonstrate that AbHeR binds to AbGS and regulates AbGS activity. More specifically, the dissociation constant (K_d) value of the binding between AbHeR and AbGS is 6.06 μM. Moreover, the absence of positively charged residues within the intracellular loop of AbHeR impacted K_d value as they serve as critical binding sites for AbGS. We also confirm that AbHeR up-regulates the biosynthetic enzyme activity of AbGS both in vitro and in vivo in the presence of light. GS is a key enzyme involved in nitrogen assimilation that catalyzes the conversion of glutamate and ammonia to glutamine. Hence, the interaction between AbHeR and AbGS may be critical for nitrogen assimilation in Actinobacteria bacterium IMCC26103 as it survives in low-nutrient environments. Overall, the findings of our study describe, for the first time, to the best of our knowledge, a novel function of heliorhodopsin as a regulatory rhodopsin with the capacity to bind and regulate enzyme activity required for nitrogen assimilation.

A study of heliorhodopsin, an actinobacterial photoreceptor of unknown function, reveals that it interacts with glutamine synthetase, an enzyme involved in nitrogen assimilation, and regulates its activity in the presence of light, highlighting the diverse functions of rhodopsins in different organisms.

Bacterial genome reductions: Tools, applications, and challenges

Bacterial cells are widely used to produce value-added products due to their versatility, ease of manipulation, and the abundance of genome engineering tools. However, the efficiency of producing these desired biomolecules is often hindered by the cells’ own metabolism, genetic instability, and the toxicity of the product. To overcome these challenges, genome reductions have been performed, making strains with the potential of serving as chassis for downstream applications. Here we review the current technologies that enable the design and construction of such reduced-genome bacteria as well as the challenges that limit their assembly and applicability. While genomic reductions have shown improvement of many cellular characteristics, a major challenge still exists in constructing these cells efficiently and rapidly. Computational tools have been created in attempts at minimizing the time needed to design these organisms, but gaps still exist in modelling these reductions in silico. Genomic reductions are a promising avenue for improving the production of value-added products, constructing chassis cells, and for uncovering cellular function but are currently limited by their time-consuming construction methods. With improvements to and the creation of novel genome editing tools and in silico models, these approaches could be combined to expedite this process and create more streamlined and efficient cell factories.

Characterization of a novel β-alanine biosynthetic pathway consisting of promiscuous metabolic enzymes

Bacteria adapt to utilize the nutrients available in their environment through a sophisticated metabolic system composed of highly specialized enzymes. Although these enzymes can metabolize molecules other than those for which they evolved, their efficiency toward promiscuous substrates is considered too low to be of physiological relevance. Herein, we investigated the possibility that these promiscuous enzymes are actually efficient enough at metabolizing secondary substrates to modify the phenotype of the cell. For example, in the bacterium Acinetobacter baylyi ADP1 (ADP1), panD (coding for l-aspartate decarboxylase) encodes the only protein known to catalyze the synthesis of β-alanine, an obligate intermediate in CoA synthesis. However, we show that the ADP1 ΔpanD mutant could also form this molecule through an unknown metabolic pathway arising from promiscuous enzymes and grow as efficiently as the wildtype strain. Using metabolomic analyses, we identified 1,3-diaminopropane and 3-aminopropanal as intermediates in this novel pathway. We also conducted activity screening and enzyme kinetics to elucidate candidate enzymes involved in this pathway, including 2,4-diaminobutyrate aminotransferase (Dat) and 2,4-diaminobutyrate decarboxylase (Ddc) and validated this pathway in vivo by analyzing the phenotype of mutant bacterial strains. Finally, we experimentally demonstrate that this novel metabolic route is not restricted to ADP1. We propose that the occurrence of conserved genes in hundreds of genomes across many phyla suggests that this previously undescribed pathway is widespread in prokaryotes.

Gradients in gene essentiality reshape antibacterial research

Essential genes encode the processes that are necessary for life. Until recently, commonly applied binary classifications left no space between essential and non-essential genes. In this review, we frame bacterial gene essentiality in the context of genetic networks. We explore how the quantitative properties of gene essentiality are influenced by the nature of the encoded process, environmental conditions and genetic background, including a strain's distinct evolutionary history. The covered topics have important consequences for antibacterials, which inhibit essential processes. We argue that the quantitative properties of essentiality can thus be used to prioritize antibacterial cellular targets and desired spectrum of activity in specific infection settings. We summarize our points with a case study on the core essential genome of the cystic fibrosis pathobiome and highlight avenues for targeted antibacterial development.

Reduction of the Bacterial Genome by Transposon-Mediated Random Deletion

Genome reduction is an important strategy in synthetic biology for constructing functional chassis cells or minimal genomes. However, the limited knowledge of complex gene functions and interactions makes genome reduction by rational design encounter a bottleneck. Here, we present an iterative and random genome reduction method for Escherichia coli, named "transposon-mediated random deletion (TMRD)". TMRD generates random double-strand breaks (DSBs) in the genome by combining Tn5 transposition with the CRISPR/Cas9 system and allows genomic deletions of various sizes at random positions during DSB repair through the intracellular alternative end-joining mechanism. Using E. coli MG1655 as the original strain, a pool of cells with multiple random genomic deletions were obtained after five reduction cycles. The growth rates of the obtained cells were comparable to that of MG1655, while the electroporation efficiency increased by at least 2 magnitudes. TMRD can generate a small E. coli library carrying multiple and random genomic deletions while enriching the cells with environmental fitness in the population. TMRD has the potential to be widely applied in the construction of minimal genomes or chassis cells for metabolic engineering.

Second-Shell Amino Acid R266 Helps Determine N-Succinylamino Acid Racemase Reaction Specificity in Promiscuous N-Succinylamino Acid Racemase/o-Succinylbenzoate Synthase Enzymes

Catalytic promiscuity is the coincidental ability to catalyze nonbiological reactions in the same active site as the native biological reaction. Several lines of evidence show that catalytic promiscuity plays a role in the evolution of new enzyme functions. Thus, studying catalytic promiscuity can help identify structural features that predispose an enzyme to evolve new functions. This study identifies a potentially preadaptive residue in a promiscuous N-succinylamino acid racemase/o-succinylbenzoate synthase (NSAR/OSBS) enzyme from Amycolatopsis sp. T-1-60. This enzyme belongs to a branch of the OSBS family which includes many catalytically promiscuous NSAR/OSBS enzymes. R266 is conserved in all members of the NSAR/OSBS subfamily. However, the homologous position is usually hydrophobic in other OSBS subfamilies, whose enzymes lack NSAR activity. The second-shell amino acid R266 is close to the catalytic acid/base K263, but it does not contact the substrate, suggesting that R266 could affect the catalytic mechanism. Mutating R266 to glutamine in Amycolatopsis NSAR/OSBS profoundly reduces NSAR activity but moderately reduces OSBS activity. This is due to a 1000-fold decrease in the rate of proton exchange between the substrate and the general acid/base catalyst K263. This mutation is less deleterious for the OSBS reaction because K263 forms a cation-π interaction with the OSBS substrate and/or the intermediate, rather than acting as a general acid/base catalyst. Together, the data explain how R266 contributes to NSAR reaction specificity and was likely an essential preadaptation for the evolution of NSAR activity.

Building Endogenous Gene Connections through RNA Self-Assembly Controlled CRISPR/Cas9 Function

Construction of synthetic circuits that can artificially establish endogenous gene connections is essential to introduce new phenotypes for cellular behaviors. Given the diversity of endogenous genes, it lacks a general and easy-to-design toolbox to manipulate the genetic network. Here we present a type of self-assembly-induced RNA circuit that can directly build regulatory connections between endogenous genes. Inspired from the natural assembling process of guide RNA in the CRISPR/Cas9 complex, this design employs an independent trigger RNA strand to induce the formation of a ternary guide RNA assembly for functional control of CRISPR/Cas9. With this general principle, expressional regulations of endogenous genes can be controlled by totally independent endogenous small RNAs and mRNAs in E. coli via activatable CRISPR/Cas9 function. Moreover, the cellular phenotype of E. coli is successfully programmed with introduction of new gene connections. In addition, the functionality of this design is also verified in the mammalian system. This self-assembly-based RNA circuit exhibits a great flexibility and simplicity of design and provides a unique approach to build endogenous gene connections, which paves a broad way toward manipulation of cellular genetic networks.

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.

Analysis of metabolic network disruption in engineered microbial hosts due to enzyme promiscuity

Increasing understanding of metabolic and regulatory networks underlying microbial physiology has enabled creation of progressively more complex synthetic biological systems for biochemical, biomedical, agricultural, and environmental applications. However, despite best efforts, confounding phenotypes still emerge from unforeseen interplay between biological parts, and the design of robust and modular biological systems remains elusive. Such interactions are difficult to predict when designing synthetic systems and may manifest during experimental testing as inefficiencies that need to be overcome. Transforming organisms such as Escherichia coli into microbial factories is achieved via several engineering strategies, used individually or in combination, with the goal of maximizing the production of chosen target compounds. One technique relies on suppressing or overexpressing selected genes; another involves introducing heterologous enzymes into a microbial host. These modifications steer mass flux towards the set of desired metabolites but may create unexpected interactions. In this work, we develop a computational method, termed Metabolic Disruption Workflow (MDFlow), for discovering interactions and network disruptions arising from enzyme promiscuity – the ability of enzymes to act on a wide range of molecules that are structurally similar to their native substrates. We apply MDFlow to two experimentally verified cases where strains with essential genes knocked out are rescued by interactions resulting from overexpression of one or more other genes. We demonstrate how enzyme promiscuity may aid cells in adapting to disruptions of essential metabolic functions. We then apply MDFlow to predict and evaluate a number of putative promiscuous reactions that can interfere with two heterologous pathways designed for 3-hydroxypropionic acid (3-HP) production. Using MDFlow, we can identify putative enzyme promiscuity and the subsequent formation of unintended and undesirable byproducts that are not only disruptive to the host metabolism but also to the intended end-objective of high biosynthetic productivity and yield. As we demonstrate, MDFlow provides an innovative workflow to systematically identify incompatibilities between the native metabolism of the host and its engineered modifications due to enzyme promiscuity.

•

Engineering modifications to cellular hosts result in undesirable byproducts.

•

Metabolic Disruption: changes in engineered host due to enzyme promiscuity.

•

Metabolic Disruption Workflow (MDFlow) uncovers metabolic disruption.

•

MDFlow corroborates previously experimentally verified promiscuous interactions.

•

MDFlow compares disruption due to heterologous pathways targeting 3-HP production.

The impact of genetic diversity on gene essentiality within the Escherichia coli species

Bacteria from the same species can differ widely in their gene content. In Escherichia coli, the set of genes shared by all strains, known as the core genome, represents about half the number of genes present in any strain. Although recent advances in bacterial genomics have unravelled genes required for fitness in various experimental conditions, most studies have focused on single model strains. As a result, the impact of the species' genetic diversity on core processes of the bacterial cell remains largely under-investigated. Here, we have developed a CRISPR interference platform for high-throughput gene repression that is compatible with most E. coli isolates and closely related species. We have applied it to assess the importance of ~3,400 nearly ubiquitous genes in three growth conditions in 18 representative E. coli strains spanning most common phylogroups and lifestyles of the species. Our screens revealed extensive variations in gene essentiality between strains and conditions. Investigation of the genetic determinants for these variations highlighted the importance of epistatic interactions with mobile genetic elements. In particular, we have shown how prophage-encoded defence systems against phage infection can trigger the essentiality of persistent genes that are usually non-essential. This study provides broad insights into the evolvability of gene essentiality and argues for the importance of studying various isolates from the same species under diverse conditions.

L-Alanine Prototrophic Suppressors Emerge from L-Alanine Auxotroph through Stress-Induced Mutagenesis in Escherichia coli

In Escherichia coli, L-alanine is synthesized by three isozymes: YfbQ, YfdZ, and AvtA. When an E. coli L-alanine auxotrophic isogenic mutant lacking the three isozymes was grown on L-alanine-deficient minimal agar medium, L-alanine prototrophic mutants emerged considerably more frequently than by spontaneous mutation; the emergence frequency increased over time, and, in an L-alanine-supplemented minimal medium, correlated inversely with L-alanine concentration, indicating that the mutants were derived through stress-induced mutagenesis. Whole-genome analysis of 40 independent L-alanine prototrophic mutants identified 16 and 18 clones harboring point mutation(s) in pyruvate dehydrogenase complex and phosphotransacetylase-acetate kinase pathway, which respectively produce acetyl coenzyme A and acetate from pyruvate. When two point mutations identified in L-alanine prototrophic mutants, in pta (D656A) and aceE (G147D), were individually introduced into the original L-alanine auxotroph, the isogenic mutants exhibited almost identical growth recovery as the respective cognate mutants. Each original- and isogenic-clone pair carrying the pta or aceE mutation showed extremely low phosphotransacetylase or pyruvate dehydrogenase activity, respectively. Lastly, extracellularly-added pyruvate, which dose-dependently supported L-alanine auxotroph growth, relieved the L-alanine starvation stress, preventing the emergence of L-alanine prototrophic mutants. Thus, L-alanine starvation-provoked stress-induced mutagenesis in the L-alanine auxotroph could lead to intracellular pyruvate increase, which eventually induces L-alanine prototrophy.

Underground isoleucine biosynthesis pathways in E. coli

The promiscuous activities of enzymes provide fertile ground for the evolution of new metabolic pathways. Here, we systematically explore the ability of E. coli to harness underground metabolism to compensate for the deletion of an essential biosynthetic pathway. By deleting all threonine deaminases, we generated a strain in which isoleucine biosynthesis was interrupted at the level of 2-ketobutyrate. Incubation of this strain under aerobic conditions resulted in the emergence of a novel 2-ketobutyrate biosynthesis pathway based upon the promiscuous cleavage of O-succinyl-L-homoserine by cystathionine γ-synthase (MetB). Under anaerobic conditions, pyruvate formate-lyase enabled 2-ketobutyrate biosynthesis from propionyl-CoA and formate. Surprisingly, we found this anaerobic route to provide a substantial fraction of isoleucine in a wild-type strain when propionate is available in the medium. This study demonstrates the selective advantage underground metabolism offers, providing metabolic redundancy and flexibility which allow for the best use of environmental carbon sources.

Genomic and phenotypic evolution of Escherichia coli in a novel citrate-only resource environment

Evolutionary innovations allow populations to colonize new ecological niches. We previously reported that aerobic growth on citrate (Cit⁺) evolved in an Escherichia coli population during adaptation to a minimal glucose medium containing citrate (DM25). Cit⁺ variants can also grow in citrate-only medium (DM0), a novel environment for E. coli. To study adaptation to this niche, we founded two sets of Cit⁺ populations and evolved them for 2500 generations in DM0 or DM25. The evolved lineages acquired numerous parallel mutations, many mediated by transposable elements. Several also evolved amplifications of regions containing the maeA gene. Unexpectedly, some evolved populations and clones show apparent declines in fitness. We also found evidence of substantial cell death in Cit⁺ clones. Our results thus demonstrate rapid trait refinement and adaptation to the new citrate niche, while also suggesting a recalcitrant mismatch between E. coli physiology and growth on citrate.

Promiscuous activity of 3-isopropylmalate dehydrogenase produced at physiological level affords Escherichia coli growth on d-malate

Promiscuous activities of enzymes may serve as starting points for the evolution of new functions. However, most experimental examples of promiscuity affording an observable phenotype necessitate the artificial overexpression of the target enzyme. Here, we show that 3-isopropylmalate dehydrogenase (IPMDH), an enzyme involved in leucine biosynthesis, has a secondary activity on d-malate, which is sufficient for d-malate assimilation under physiological conditions where the enzyme is upregulated. In vitro, the turnover constant (k_cat ) of IPMDH for d-malate is about 30-fold lower than the k_cat for 3-isopropylmalate, yet sufficiently high to support the growth on d-malate. From an evolutionary perspective, our results highlight the possibility of phenotype emergence triggered by arbitrary changes in environmental conditions and prior to any mutational event.

Genome-Wide CRISPRi-Based Identification of Targets for Decoupling Growth from Production

Growth decoupling can be used to optimize microbial production of biobased compounds by inhibiting excess biomass formation and redirect carbon flux from growth to product formation. However, identifying suitable genetic targets through rational design is challenging. Here, we conduct a genome-wide CRISPRi screen to discover growth switches suitable for decoupling growth and production. Using an sgRNA library covering 12 238 loci in the Escherichia coli genome, we screen for targets that inhibit growth while allowing for continued protein production. In total, we identify 1332 sgRNAs that simultaneously decrease growth and maintain or increase accumulation of GFP. The top target sibB/ibsB shows more than 5-fold increase in GFP accumulation and 45% decrease in biomass formation. Overall, our genome-wide CRISPRi screen provides key targets for growth decoupling, and the approach can be applied to improve biobased production in other microorganisms.

In Vitro Selection of Peptides and Proteins-Advantages of mRNA Display

mRNA display is a robust in vitro selection technique that allows the selection of peptides and proteins with desired functions from libraries of trillions of variants. mRNA display relies upon a covalent linkage between a protein and its encoding mRNA molecule; the power of the technique stems from the stability of this link, and the large degree of control over experimental conditions afforded to the researcher. This article describes the major advantages that make mRNA display the method of choice among comparable in vivo and in vitro methods, including cell-surface display, phage display, and ribosomal display. We also describe innovative techniques that harness mRNA display for directed evolution, protein engineering, and drug discovery.

How enzyme promiscuity and horizontal gene transfer contribute to metabolic innovation

Promiscuity is the coincidental ability of an enzyme to catalyze its native reaction and additional reactions that are not biological functions in the same active site. Promiscuity plays a central role in enzyme evolution and is thus a useful property for protein and metabolic engineering. This review examines enzyme evolution holistically, beginning with evaluating biochemical support for four enzyme evolution models. As expected, there is strong biochemical support for the subfunctionalization and innovation-amplification-divergence models, in which promiscuity is a central feature. In many cases, however, enzyme evolution is more complex than the models indicate, suggesting much is yet to be learned about selective pressures on enzyme function. A complete understanding of enzyme evolution must also explain the ability of metabolic networks to integrate new enzyme activities. Hidden within metabolic networks are underground metabolic pathways constructed from promiscuous activities. We discuss efforts to determine the diversity and pervasiveness of underground metabolism. Remarkably, several studies have discovered that some metabolic defects can be repaired via multiple underground routes. In prokaryotes, metabolic innovation is driven by connecting enzymes acquired by horizontal gene transfer (HGT) into the metabolic network. Thus, we end the review by discussing how the combination of promiscuity and HGT contribute to evolution of metabolism in prokaryotes. Future studies investigating the contribution of promiscuity to enzyme and metabolic evolution will need to integrate deeper probes into the influence of evolution on protein biophysics, enzymology, and metabolism with more complex and realistic evolutionary models. ENZYMES: lactate dehydrogenase (EC 1.1.1.27), malate dehydrogenase (EC 1.1.1.37), OSBS (EC 4.2.1.113), HisA (EC 5.3.1.16), TrpF, PriA (EC 5.3.1.24), R-mandelonitrile lyase (EC 4.1.2.10), Maleylacetate reductase (EC 1.3.1.32).

Library Selection with a Randomized Repertoire of (βα)8-Barrel Enzymes Results in Unexpected Induction of Gene Expression

The potential of the frequently encountered (βα)₈-barrel fold to acquire new functions was tested by an approach combining random mutagenesis and selection in vivo. For this purpose, the genes encoding 52 different phosphate-binding (βα)₈-barrel proteins were subjected to error-prone PCR and cloned into an expression plasmid. The resulting mixed repertoire was used to transform different auxotrophic Escherichia coli strains, each lacking an enzyme with a phosphate-containing substrate. After plating of the different transformants on minimal medium, growth was observed only for two strains, lacking either the gene for the serine phosphatase SerB or the phosphoserine aminotransferase SerC. The same mutants of the E. coli genes nanE (encoding a putative N-acetylmannosamine-6-phosphate 2-epimerase) and pdxJ (encoding the pyridoxine 5'-phosphate synthase) were responsible for rescuing both ΔserB and ΔserC. Unexpectedly, the complementing NanE and PdxJ variants did not catalyze the SerB or SerC reactions in vitro. Instead, RT-qPCR, RNAseq, and transcriptome analysis showed that they rescue the deletions by enlisting the help of endogenous E. coli enzymes HisB and HisC through exclusive up-regulation of histidine operon transcription. While the promiscuous SerB activity of HisB is well-established, our data indicate that HisC is promiscuous for the SerC reaction, as well. The successful rescue of ΔserB and ΔserC through point mutations and recruitment of additional amino acids in NanE and PdxJ provides another example for the adaptability of the (βα)₈-barrel fold.

Systematic analysis reveals the prevalence and principles of bypassable gene essentiality

Gene essentiality is a variable phenotypic trait, but to what extent and how essential genes can become dispensable for viability remain unclear. Here, we investigate 'bypass of essentiality (BOE)' - an underexplored type of digenic genetic interaction that renders essential genes dispensable. Through analyzing essential genes on one of the six chromosome arms of the fission yeast Schizosaccharomyces pombe, we find that, remarkably, as many as 27% of them can be converted to non-essential genes by BOE interactions. Using this dataset we identify three principles of essentiality bypass: bypassable essential genes tend to have lower importance, tend to exhibit differential essentiality between species, and tend to act with other bypassable genes. In addition, we delineate mechanisms underlying bypassable essentiality, including the previously unappreciated mechanism of dormant redundancy between paralogs. The new insights gained on bypassable essentiality deepen our understanding of genotype-phenotype relationships and will facilitate drug development related to essential genes.

Genetic Code Evolution Investigated through the Synthesis and Characterisation of Proteins from Reduced-Alphabet Libraries

The universal genetic code of 20 amino acids is the product of evolution. It is believed that earlier versions of the code had fewer residues. Many theories for the order in which amino acids were integrated into the code have been proposed, considering factors ranging from prebiotic chemistry to codon capture. Several meta-analyses combined these theories to yield a feasible consensus chronology of the genetic code's evolution, but there is a dearth of experimental data to test the hypothesised order. We used combinatorial chemistry to synthesise libraries of random polypeptides that were based on different subsets of the 20 standard amino acids, thus representing different stages of a plausible history of the alphabet. Four libraries were comprised of the five, nine, and 16 most ancient amino acids, and all 20 extant residues for a direct side-by-side comparison. We characterised numerous variants from each library for their solubility and propensity to form secondary, tertiary or quaternary structures. Proteins from the two most ancient libraries were more likely to be soluble than those from the extant library. Several individual protein variants exhibited inducible protein folding and other traits typical of intrinsically disordered proteins. From these libraries, we can infer how primordial protein structure and function might have evolved with the genetic code.

Dual-barcoded shotgun expression library sequencing for high-throughput characterization of functional traits in bacteria

A major challenge in genomics is the knowledge gap between sequence and its encoded function. Gain-of-function methods based on gene overexpression are attractive avenues for phenotype-based functional screens, but are not easily applied in high-throughput across many experimental conditions. Here, we present Dual Barcoded Shotgun Expression Library Sequencing (Dub-seq), a method that uses random DNA barcodes to greatly increase experimental throughput. As a demonstration of this approach, we construct a Dub-seq library with Escherichia coli genomic DNA, performed 155 genome-wide fitness assays in 52 experimental conditions, and identified overexpression phenotypes for 813 genes. We show that Dub-seq data is reproducible, accurately recapitulates known biology, and identifies hundreds of novel gain-of-function phenotypes for E. coli genes, a subset of which we verified with assays of individual strains. Dub-seq provides complementary information to loss-of-function approaches and will facilitate rapid and systematic functional characterization of microbial genomes.

Gain of function methods based on gene overexpression are not easily applied to high-throughput screening across different experimental conditions. Here, the authors present Dub-seq, which separates overexpression library characterization from functional screening and uses random DNA barcodes to increase the experimental throughput.

Using a Chemical Genetic Screen to Enhance Our Understanding of the Antimicrobial Properties of Gallium against Escherichia coli

The diagnostic and therapeutic agent gallium offers multiple clinical and commercial uses including the treatment of cancer and the localization of tumors, among others. Further, this metal has been proven to be an effective antimicrobial agent against a number of microbes. Despite the latter, the fundamental mechanisms of gallium action have yet to be fully identified and understood. To further the development of this antimicrobial, it is imperative that we understand the mechanisms by which gallium interacts with cells. As a result, we screened the Escherichia coli Keio mutant collection as a means of identifying the genes that are implicated in prolonged gallium toxicity or resistance and mapped their biological processes to their respective cellular system. We discovered that the deletion of genes functioning in response to oxidative stress, DNA or iron–sulfur cluster repair, and nucleotide biosynthesis were sensitive to gallium, while Ga resistance comprised of genes involved in iron/siderophore import, amino acid biosynthesis and cell envelope maintenance. Altogether, our explanations of these findings offer further insight into the mechanisms of gallium toxicity and resistance in E. coli.

Balancing Specificity and Promiscuity in Enzyme Evolution: Multidimensional Activity Transitions in the Alkaline Phosphatase Superfamily

Highly proficient, promiscuous enzymes can be springboards for functional evolution, able to avoid loss of function during adaptation by their capacity to promote multiple reactions. We employ a systematic comparative study of structure, sequence, and substrate specificity to track the evolution of specificity and reactivity between promiscuous members of clades of the alkaline phosphatase (AP) superfamily. Construction of a phylogenetic tree of protein sequences maps out the likely transition zone between arylsulfatases (ASs) and phosphonate monoester hydrolases (PMHs). Kinetic analysis shows that all enzymes characterized have four chemically distinct phospho- and sulfoesterase activities, with rate accelerations ranging from 10¹¹- to 10¹⁷-fold for their primary and 10⁹- to 10¹²-fold for their promiscuous reactions, suggesting that catalytic promiscuity is widespread in the AP-superfamily. This functional characterization and crystallography reveal a novel class of ASs that is so similar in sequence to known PMHs that it had not been recognized as having diverged in function. Based on analysis of snapshots of catalytic promiscuity "in transition", we develop possible models that would allow functional evolution and determine scenarios for trade-off between multiple activities. For the new ASs, we observe largely invariant substrate specificity that would facilitate the transition from ASs to PMHs via trade-off-free molecular exaptation, that is, evolution without initial loss of primary activity and specificity toward the original substrate. This ability to bypass low activity generalists provides a molecular solution to avoid adaptive conflict.

Reframing gene essentiality in terms of adaptive flexibility

BACKGROUND

Essentiality assays are important tools commonly utilized for the discovery of gene functions. Growth/no growth screens of single gene knockout strain collections are also often utilized to test the predictive power of genome-scale models. False positive predictions occur when computational analysis predicts a gene to be non-essential, however experimental screens deem the gene to be essential. One explanation for this inconsistency is that the model contains the wrong information, possibly an incorrectly annotated alternative pathway or isozyme reaction. Inconsistencies could also be attributed to experimental limitations, such as growth tests with arbitrary time cut-offs. The focus of this study was to resolve such inconsistencies to better understand isozyme activities and gene essentiality.

RESULTS

In this study, we explored the definition of conditional essentiality from a phenotypic and genomic perspective. Gene-deletion strains associated with false positive predictions of gene essentiality on defined minimal medium for Escherichia coli were targeted for extended growth tests followed by population sequencing and transcriptome analysis. Of the twenty false positive strains available and confirmed from the Keio single gene knock-out collection, 11 strains were shown to grow with longer incubation periods making these actual true positives. These strains grew reproducibly with a diverse range of growth phenotypes. The lag phase observed for these strains ranged from less than one day to more than 7 days. It was found that 9 out of 11 of the false positive strains that grew acquired mutations in at least one replicate experiment and the types of mutations ranged from SNPs and small indels associated with regulatory or metabolic elements to large regions of genome duplication. Comparison of the detected adaptive mutations, modeling predictions of alternate pathways and isozymes, and transcriptome analysis of KO strains suggested agreement for the observed growth phenotype for 6 out of the 9 cases where mutations were observed.

CONCLUSIONS

Longer-term growth experiments followed by whole genome sequencing and transcriptome analysis can provide a better understanding of conditional gene essentiality and mechanisms of adaptation to such perturbations. Compensatory mutations are largely reproducible mechanisms and are in agreement with genome-scale modeling predictions to loss of function gene deletion events.

Structural and functional insight into serine hydroxymethyltransferase from Helicobacter pylori

Serine hydroxymethyltransferase (SHMT), encoded by the glyA gene, is a ubiquitous pyridoxal 5'-phosphate (PLP)-dependent enzyme that catalyzes the formation of glycine from serine. The thereby generated 5,10-methylene tetrahydrofolate (MTHF) is a major source of cellular one-carbon units and a key intermediate in thymidylate biosynthesis. While in virtually all eukaryotic and many bacterial systems thymidylate synthase ThyA, SHMT and dihydrofolate reductase (DHFR) are part of the thymidylate/folate cycle, the situation is different in organisms using flavin-dependent thymidylate synthase ThyX. Here the distinct catalytic reaction directly produces tetrahydrofolate (THF) and consequently in most ThyX-containing organisms, DHFR is absent. While the resulting influence on the folate metabolism of ThyX-containing bacteria is not fully understood, the presence of ThyX may provide growth benefits under conditions where the level of reduced folate derivatives is compromised. Interestingly, the third key enzyme implicated in generation of MTHF, serine hydroxymethyltransferase (SHMT), has a universal phylogenetic distribution, but remains understudied in ThyX-containg bacteria. To obtain functional insight into these ThyX-dependent thymidylate/folate cycles, we characterized the predicted SHMT from the ThyX-containing bacterium Helicobacter pylori. Serine hydroxymethyltransferase activity was confirmed by functional genetic complementation of a glyA-inactivated E. coli strain. A H. pylori ΔglyA strain was obtained, but exhibited markedly slowed growth and had lost the virulence factor CagA. Biochemical and spectroscopic evidence indicated formation of a characteristic enzyme-PLP-glycine-folate complex and revealed unexpectedly weak binding affinity of PLP. The three-dimensional structure of the H. pylori SHMT apoprotein was determined at 2.8Ǻ resolution, suggesting a structural basis for the low affinity of the enzyme for its cofactor. Stabilization of the proposed inactive configuration using small molecules has potential to provide a specific way for inhibiting HpSHMT.

Metabolic models and gene essentiality data reveal essential and conserved metabolism in prokaryotes

Essential metabolic reactions are shaping constituents of metabolic networks, enabling viable and distinct phenotypes across diverse life forms. Here we analyse and compare modelling predictions of essential metabolic functions with experimental data and thereby identify core metabolic pathways in prokaryotes. Simulations of 15 manually curated genome-scale metabolic models were integrated with 36 large-scale gene essentiality datasets encompassing a wide variety of species of bacteria and archaea. Conservation of metabolic genes was estimated by analysing 79 representative genomes from all the branches of the prokaryotic tree of life. We find that essentiality patterns reflect phylogenetic relations both for modelling and experimental data, which correlate highly at the pathway level. Genes that are essential for several species tend to be highly conserved as opposed to non-essential genes which may be conserved or not. The tRNA-charging module is highlighted as ancestral and with high centrality in the networks, followed closely by cofactor metabolism, pointing to an early information processing system supplied by organic cofactors. The results, which point to model improvements and also indicate faults in the experimental data, should be relevant to the study of centrality in metabolic networks and ancient metabolism but also to metabolic engineering with prokaryotes.

If we tried to list every known chemical reaction within an organism–human, plant or even bacteria–we would get quite a long and confusing read. But when this information is represented in so-called genome-scale metabolic networks, we have the means to access computationally each of those reactions and their interconnections. Some parts of the network have alternatives, while others are unique and therefore can be essential for growth. Here, we simulate growth and compare essential reactions and genes for the simplest type of unicellular species–prokaryotes–to understand which parts of their metabolism are universally essential and potentially ancestral. We show that similar patterns of essential reactions echo phylogenetic relationships (this makes sense, as the genome provides the building plan for the enzymes that perform those reactions). Our computational predictions correlate strongly with experimental essentiality data. Finally, we show that a crucial step of protein synthesis (tRNA charging) and the synthesis and transformation of small molecules that enzymes require (cofactors) are the most essential and conserved parts of metabolism in prokaryotes. Our results are a step further in understanding the biology and evolution of prokaryotes but can also be relevant in applied studies including metabolic engineering and antibiotic design.

The EcoCyc Database

EcoCyc is a bioinformatics database available at EcoCyc.org that describes the genome and the biochemical machinery of Escherichia coli K-12 MG1655. The long-term goal of the project is to describe the complete molecular catalog of the E. coli cell, as well as the functions of each of its molecular parts, to facilitate a system-level understanding of E. coli. EcoCyc is an electronic reference source for E. coli biologists and for biologists who work with related microorganisms. The database includes information pages on each E. coli gene product, metabolite, reaction, operon, and metabolic pathway. The database also includes information on E. coli gene essentiality and on nutrient conditions that do or do not support the growth of E. coli. The website and downloadable software contain tools for analysis of high-throughput data sets. In addition, a steady-state metabolic flux model is generated from each new version of EcoCyc and can be executed via EcoCyc.org. The model can predict metabolic flux rates, nutrient uptake rates, and growth rates for different gene knockouts and nutrient conditions. This review outlines the data content of EcoCyc and of the procedures by which this content is generated.

Mutational and transcriptional landscape of spontaneous gene duplications and deletions in Caenorhabditis elegans

Gene duplication and deletion are pivotal processes shaping the structural and functional repertoire of genomes, with implications for disease, adaptation, and evolution. We employed a mutation accumulation (MA) framework partnered with high-throughput genomics to assess the molecular and transcriptional characteristics of newly arisen gene copy-number variants (CNVs) in Caenorhabditis elegans populations subjected to varying intensity of selection. Here, we report a direct spontaneous genome-wide rate of gene duplication of 2.9 × 10^-5/gene per generation in C. elegans, the highest for any species to date. The rate of gene deletion is sixfold lower (5 × 10^-6/gene per generation). Deletions of highly expressed genes are particularly deleterious, given their paucity in even the N = 1 lines with minimal efficacy of selection. The increase in average transcript abundance of new duplicates arising under minimal selection is significantly greater than twofold compared with single copies of the same gene, suggesting that genes in segmental duplications are frequently overactive at inception. The average increase in transcriptional activity of gene duplicates is greater in the N = 1 MA lines than in MA lines with larger population bottlenecks. There is an inverse relationship between the ancestral transcription levels of new gene duplicates and population size, with duplicate copies of highly expressed genes less likely to accumulate in larger populations. Our results demonstrate a fitness cost of increased transcription following duplication, which results in purifying selection against new gene duplicates. However, on average, duplications also provide a significant increase in gene expression that can facilitate adaptation to novel environmental challenges.

Metabolite–Enzyme Coevolution: From Single Enzymes to Metabolic Pathways and Networks

How individual enzymes evolved is relatively well understood. However, individual enzymes rarely confer a physiological advantage on their own. Judging by its current state, the emergence of metabolism seemingly demanded the simultaneous emergence of many enzymes. Indeed, how multicomponent interlocked systems, like metabolic pathways, evolved is largely an open question. This complexity can be unlocked if we assume that survival of the fittest applies not only to genes and enzymes but also to the metabolites they produce. This review develops our current knowledge of enzyme evolution into a wider hypothesis of pathway and network evolution. We describe the current models for pathway evolution and offer an integrative metabolite–enzyme coevolution hypothesis. Our hypothesis addresses the origins of new metabolites and of new enzymes and the order of their recruitment. We aim to not only survey established knowledge but also present open questions and potential ways of addressing them.

Comparison of Alicyclobacillus acidocaldarius o-Succinylbenzoate Synthase to Its Promiscuous N-Succinylamino Acid Racemase/ o-Succinylbenzoate Synthase Relatives

Studying the evolution of catalytically promiscuous enzymes like those from the N-succinylamino acid racemase/ o-succinylbenzoate synthase (NSAR/OSBS) subfamily can reveal mechanisms by which new functions evolve. Some enzymes in this subfamily have only OSBS activity, while others catalyze OSBS and NSAR reactions. We characterized several NSAR/OSBS subfamily enzymes as a step toward determining the structural basis for evolving NSAR activity. Three enzymes were promiscuous, like most other characterized NSAR/OSBS subfamily enzymes. However, Alicyclobacillus acidocaldarius OSBS (AaOSBS) efficiently catalyzes OSBS activity but lacks detectable NSAR activity. Competitive inhibition and molecular modeling show that AaOSBS binds N-succinylphenylglycine with moderate affinity in a site that overlaps its normal substrate. On the basis of possible steric conflicts identified by molecular modeling and sequence conservation within the NSAR/OSBS subfamily, we identified one mutation, Y299I, that increased NSAR activity from undetectable to 1.2 × 10² M^-1 s^-1 without affecting OSBS activity. This mutation does not appear to affect binding affinity but instead affects k_cat, by reorienting the substrate or modifying conformational changes to allow both catalytic lysines to access the proton that is moved during the reaction. This is the first site known to affect reaction specificity in the NSAR/OSBS subfamily. However, this gain of activity was obliterated by a second mutation, M18F. Epistatic interference by M18F was unexpected because a phenylalanine at this position is important in another NSAR/OSBS enzyme. Together, modest NSAR activity of Y299I AaOSBS and epistasis between sites 18 and 299 indicate that additional sites influenced the evolution of NSAR reaction specificity in the NSAR/OSBS subfamily.

A bacteriophage enzyme induces bacterial metabolic perturbation that confers a novel promiscuous function

One key concept in the evolution of new functions is the ability of enzymes to perform promiscuous side-reactions that serve as a source of novelty that may become beneficial under certain conditions. Here, we identify a mechanism where a bacteriophage-encoded enzyme introduces novelty by inducing expression of a promiscuous bacterial enzyme. By screening for bacteriophage DNA that rescued an auxotrophic Escherichia coli mutant carrying a deletion of the ilvA gene, we show that bacteriophage-encoded S-adenosylmethionine (SAM) hydrolases reduce SAM levels. Through this perturbation of bacterial metabolism, expression of the promiscuous bacterial enzyme MetB is increased, which in turn complements the absence of IlvA. These results demonstrate how foreign DNA can increase the metabolic capacity of bacteria, not only by transfer of bona fide new genes, but also by bringing cryptic bacterial functions to light via perturbations of cellular physiology.

Factors driving metabolic diversity in the budding yeast subphylum

Background

Associations between traits are prevalent in nature, occurring across a diverse range of taxa and traits. Individual traits may co-evolve with one other, and these correlations can be driven by factors intrinsic or extrinsic to an organism. However, few studies, especially in microbes, have simultaneously investigated both across a broad taxonomic range. Here we quantify pairwise associations among 48 traits across 784 diverse yeast species of the ancient budding yeast subphylum Saccharomycotina, assessing the effects of phylogenetic history, genetics, and ecology.

Results

We find extensive negative (traits that tend to not occur together) and positive (traits that tend to co-occur) pairwise associations among traits, as well as between traits and environments. These associations can largely be explained by the biological properties of the traits, such as overlapping biochemical pathways. The isolation environments of the yeasts explain a minor but significant component of the variance, while phylogeny (the retention of ancestral traits in descendant species) plays an even more limited role. Positive correlations are pervasive among carbon utilization traits and track with chemical structures (e.g., glucosides and sugar alcohols) and metabolic pathways, suggesting a molecular basis for the presence of suites of traits. In several cases, characterized genes from model organisms suggest that enzyme promiscuity and overlapping biochemical pathways are likely mechanisms to explain these macroevolutionary trends. Interestingly, fermentation traits are negatively correlated with the utilization of pentose sugars, which are major components of the plant biomass degraded by fungi and present major bottlenecks to the production of cellulosic biofuels. Finally, we show that mammalian pathogenic and commensal yeasts have a suite of traits that includes growth at high temperature and, surprisingly, the utilization of a narrowed panel of carbon sources.

Conclusions

These results demonstrate how both intrinsic physiological factors and extrinsic ecological factors drive the distribution of traits present in diverse organisms across macroevolutionary timescales.

Electronic supplementary material

The online version of this article (10.1186/s12915-018-0498-3) contains supplementary material, which is available to authorized users.