Flux, toxicity, and expression costs generate complex genetic interactions in a metabolic pathway

Proteome allocation and the evolution of metabolic cross-feeding

In a common instance of metabolic cross-feeding (MCF), an organism incompletely metabolizes nutrients and releases metabolites that are used by another to produce energy or building blocks. Why would the former waste edible food, and why does this preferentially occur at specific locations in a metabolic pathway have challenged evolutionary theory for decades. To address these questions, we combine adaptive dynamics with an explicit model of cell metabolism, including enzyme-driven catalysis of metabolic reactions and the cellular constraints acting on the proteome that may incur a cost to expressing all enzymes along a pathway. After pointing out that cells should in principle prioritize upstream reactions when metabolites are restrained inside the cell, we show that the occurrence of permeability-driven MCF is rare and requires that an intermediate metabolite be extremely diffusive. Indeed, only at very high levels of membrane permeability (consistent with those of acetate and glycerol, for instance) and under distinctive sets of parameters should the population diversify and MCF evolve. These results help understand the origins of simple microbial communities, such as those that readily evolve in short-term evolutionary experiments, and may later be extended to investigate how evolution has progressively built up today's extremely diverse ecosystems.

Mutational biases favor complexity increases in protein interaction networks after gene duplication

Biological systems can gain complexity over time. While some of these transitions are likely driven by natural selection, the extent to which they occur without providing an adaptive benefit is unknown. At the molecular level, one example is heteromeric complexes replacing homomeric ones following gene duplication. Here, we build a biophysical model and simulate the evolution of homodimers and heterodimers following gene duplication using distributions of mutational effects inferred from available protein structures. We keep the specific activity of each dimer identical, so their concentrations drift neutrally without new functions. We show that for more than 60% of tested dimer structures, the relative concentration of the heteromer increases over time due to mutational biases that favor the heterodimer. However, allowing mutational effects on synthesis rates and differences in the specific activity of homo- and heterodimers can limit or reverse the observed bias toward heterodimers. Our results show that the accumulation of more complex protein quaternary structures is likely under neutral evolution, and that natural selection would be needed to reverse this tendency.

Transcriptome Analysis Reveals the Involvement of Mitophagy and Peroxisome in the Resistance to QoIs in Corynespora cassiicola

Quinone outside inhibitor fungicides (QoIs) are crucial fungicides for controlling plant diseases, but resistance, mainly caused by G143A, has been widely reported with the high and widespread use of QoIs. However, two phenotypes of Corynespora casiicola (RI and RII) with the same G143A showed significantly different resistance to QoIs in our previous study, which did not match the reported mechanisms. Therefore, transcriptome analysis of RI and RII strains after trifloxystrobin treatment was used to explore the new resistance mechanism in this study. The results show that 332 differentially expressed genes (DEGs) were significantly up-regulated and 448 DEGs were significantly down-regulated. The results of GO and KEGG enrichment showed that DEGs were most enriched in ribosomes, while also having enrichment in peroxide, endocytosis, the lysosome, autophagy, and mitophagy. In particular, mitophagy and peroxisome have been reported in medicine as the main mechanisms of reactive oxygen species (ROS) scavenging, while the lysosome and endocytosis are an important organelle and physiological process, respectively, that assist mitophagy. The oxidative stress experiments showed that the oxidative stress resistance of the RII strains was significantly higher than that of the RI strains: specifically, it was more than 1.8-fold higher at a concentration of 0.12% H2O2. This indicates that there is indeed a significant difference in the scavenging capacity of ROS between the two phenotypic strains. Therefore, we suggest that QoIs' action caused a high production of ROS, and that scavenging mechanisms such as mitophagy and peroxisomes functioned in RII strains to prevent oxidative stress, whereas RI strains were less capable of resisting oxidative stress, resulting in different resistance to QoIs. In this study, it was first revealed that mitophagy and peroxisome mechanisms available for ROS scavenging are involved in the resistance of pathogens to fungicides.

Deep mutational scanning reveals the molecular determinants of RNA polymerase-mediated adaptation and tradeoffs

RNA polymerase (RNAP) is emblematic of complex biological systems that control multiple traits involving trade-offs such as growth versus maintenance. Laboratory evolution has revealed that mutations in RNAP subunits, including RpoB, are frequently selected. However, we lack a systems view of how mutations alter the RNAP molecular functions to promote adaptation. We, therefore, measured the fitness of thousands of mutations within a region of rpoB under multiple conditions and genetic backgrounds, to find that adaptive mutations cluster in two modules. Mutations in one module favor growth over maintenance through a partial loss of an interaction associated with faster elongation. Mutations in the other favor maintenance over growth through a destabilized RNAP-DNA complex. The two molecular handles capture the versatile RNAP-mediated adaptations. Combining both interaction losses simultaneously improved maintenance and growth, challenging the idea that growth-maintenance tradeoff resorts only from limited resources, and revealing how compensatory evolution operates within RNAP.

The genotype-phenotype relationship and evolutionary genetics in the light of the Metabolic Control Analysis

Metabolic control analysis has long been used as a systemic model of the genotype-phenotype (GP) relationship. By considering kinetic parameters and enzyme concentrations as reflecting the genotype level and metabolic fluxes or pools as phenotypes related to fitness, MCA has given a biological basis to the relationship between these two levels. The non-linear and concave relationship between enzymes and fluxes can account for common genetic effects that reductionist approaches have been powerless to explain, such as the dominance of active alleles over less active alleles, the various types of epistasis and heterosis, and reveals the structural links between these genetic effects. The summation property of the flux control coefficients accounts for the L-shaped distribution of Quantitative Trait Locus (QTL) effects, irrespective of other possible causes. Metabolic models of response to selection results in evolutionary scenarios that are markedly different from those derived from the classical infinitesimal model of quantitative genetics. In particular, evolution towards selective neutrality appears to be a consequence of the diminishing return of the flux-enzyme relationship. In this paper, we survey the historical and recent achievements of MCA in genetics, quantitative genetics and evolution, focusing on epistasis and the evolution of flux in relation to enzyme concentrations.

A null model for the distribution of fitness effects of mutations

The distribution of fitness effects (DFE) of new mutations is key to our understanding of many evolutionary processes. Theoreticians have developed several models to help understand the patterns seen in empirical DFEs. Many such models reproduce the broad patterns seen in empirical DFEs but these models often rely on structural assumptions that cannot be tested empirically. Here, we investigate how much of the underlying "microscopic" biological processes involved in the mapping of new mutations to fitness can be inferred from "macroscopic" observations of the DFE. We develop a null model by generating random genotype-to-fitness maps and show that the null DFE is that with the largest possible information entropy. We further show that, subject to one simple constraint, this null DFE is a Gompertz distribution. Finally, we illustrate how the predictions of this null DFE match empirically measured DFEs from several datasets, as well as DFEs simulated from Fisher's geometric model. This suggests that a match between models and empirical data is often not a very strong indication of the mechanisms underlying the mapping of mutation to fitness.

Environment-dependent epistasis increases phenotypic diversity in gene regulatory networks

Mutations to gene regulatory networks can be maladaptive or a source of evolutionary novelty. Epistasis confounds our understanding of how mutations affect the expression patterns of gene regulatory networks, a challenge exacerbated by the dependence of epistasis on the environment. We used the toolkit of synthetic biology to systematically assay the effects of pairwise and triplet combinations of mutant genotypes on the expression pattern of a gene regulatory network expressed in Escherichia coli that interprets an inducer gradient across a spatial domain. We uncovered a preponderance of epistasis that can switch in magnitude and sign across the inducer gradient to produce a greater diversity of expression pattern phenotypes than would be possible in the absence of such environment-dependent epistasis. We discuss our findings in the context of the evolution of hybrid incompatibilities and evolutionary novelties.

Epistasis decreases with increasing antibiotic pressure but not temperature

Predicting mutational effects is essential for the control of antibiotic resistance (ABR). Predictions are difficult when there are strong genotype-by-environment (G × E), gene-by-gene (G × G or epistatic) or gene-by-gene-by-environment (G × G × E) interactions. We quantified G × G × E effects in Escherichia coli across environmental gradients. We created intergenic fitness landscapes using gene knock-outs and single-nucleotide ABR mutations previously identified to vary in the extent of G × E effects in our environments of interest. Then, we measured competitive fitness across a complete combinatorial set of temperature and antibiotic dosage gradients. In this way, we assessed the predictability of 15 fitness landscapes across 12 different but related environments. We found G × G interactions and rugged fitness landscapes in the absence of antibiotic, but as antibiotic concentration increased, the fitness effects of ABR genotypes quickly overshadowed those of gene knock-outs, and the landscapes became smoother. Our work reiterates that some single mutants, like those conferring resistance or susceptibility to antibiotics, have consistent effects across genetic backgrounds in stressful environments. Thus, although epistasis may reduce the predictability of evolution in benign environments, evolution may be more predictable in adverse environments. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Designing artificial pathways for improving chemical production

Metabolic engineering exploits manipulation of catalytic and regulatory elements to improve a specific function of the host cell, often the synthesis of interesting chemicals. Although naturally occurring pathways are significant resources for metabolic engineering, these pathways are frequently inefficient and suffer from a series of inherent drawbacks. Designing artificial pathways in a rational manner provides a promising alternative for chemicals production. However, the entry barrier of designing artificial pathway is relatively high, which requires researchers a comprehensive and deep understanding of physical, chemical and biological principles. On the other hand, the designed artificial pathways frequently suffer from low efficiencies, which impair their further applications in host cells. Here, we illustrate the concept and basic workflow of retrobiosynthesis in designing artificial pathways, as well as the most currently used methods including the knowledge- and computer-based approaches. Then, we discuss how to obtain desired enzymes for novel biochemistries, and how to trim the initially designed artificial pathways for further improving their functionalities. Finally, we summarize the current applications of artificial pathways from feedstocks utilization to various products synthesis, as well as our future perspectives on designing artificial pathways.

L’impact des mutations neutres sur l’évolvabilité et l’évolution des génomes

Les mutations bénéfiques à forts effets sont rares et les mutations délétères sont éliminées par la sélection naturelle. La majorité des mutations qui s’accumulent dans les génomes ont donc des effets sélectifs très faibles, voire nuls ; elles sont alors appelées mutations neutres. Au cours des deux dernières décennies, il a été montré que les mutations, même en l’absence d’effet sur la valeur sélective des organismes, affectent leur évolvabilité, en donnant accès à de nouveaux phénotypes par le biais de mutations apparaissant ultérieurement, et qui n’auraient pas été disponibles autrement. En plus de cet effet, de nombreuses mutations neutres – indépendamment de leurs effets sélectifs – peuvent affecter la mutabilité de séquences d’ADN voisines, et moduler l’efficacité de la recombinaison homologue. De telles mutations ne modifient pas le spectre des phénotypes accessibles, mais plutôt la vitesse à laquelle de nouveaux phénotypes seront produits, un processus qui a des conséquences à long terme mais aussi potentiellement à court terme, en lien avec l’émergence de cancers.

On the incongruence of genotype-phenotype and fitness landscapes

The mapping from genotype to phenotype to fitness typically involves multiple nonlinearities that can transform the effects of mutations. For example, mutations may contribute additively to a phenotype, but their effects on fitness may combine non-additively because selection favors a low or intermediate value of that phenotype. This can cause incongruence between the topographical properties of a fitness landscape and its underlying genotype-phenotype landscape. Yet, genotype-phenotype landscapes are often used as a proxy for fitness landscapes to study the dynamics and predictability of evolution. Here, we use theoretical models and empirical data on transcription factor-DNA interactions to systematically study the incongruence of genotype-phenotype and fitness landscapes when selection favors a low or intermediate phenotypic value. Using the theoretical models, we prove a number of fundamental results. For example, selection for low or intermediate phenotypic values does not change simple sign epistasis into reciprocal sign epistasis, implying that genotype-phenotype landscapes with only simple sign epistasis motifs will always give rise to single-peaked fitness landscapes under such selection. More broadly, we show that such selection tends to create fitness landscapes that are more rugged than the underlying genotype-phenotype landscape, but this increased ruggedness typically does not frustrate adaptive evolution because the local adaptive peaks in the fitness landscape tend to be nearly as tall as the global peak. Many of these results carry forward to the empirical genotype-phenotype landscapes, which may help to explain why low- and intermediate-affinity transcription factor-DNA interactions are so prevalent in eukaryotic gene regulation.

How do mutations change phenotypic traits and organismal fitness? This question is often addressed in the context of a classic metaphor of evolutionary theory—the fitness landscape. A fitness landscape is akin to a physical landscape, in which genotypes define spatial coordinates, and fitness defines the elevation of each coordinate. Evolution then acts like a hill-climbing process, in which populations ascend fitness peaks as a consequence of mutation and selection. It is becoming increasingly common to construct such landscapes using experimental data from high-throughput sequencing technologies and phenotypic assays, in systems such as macromolecules and gene regulatory circuits. Although these landscapes are typically defined by molecular phenotypes, and are therefore more appropriately referred to as genotype-phenotype landscapes, they are often used to study evolutionary dynamics. This requires the assumption that the molecular phenotype is a reasonable proxy for fitness, which need not be the case. For example, selection may favor a low or intermediate phenotypic value, causing incongruence between a fitness landscape and its underlying genotype-phenotype landscape. Here, we study such incongruence using a diversity of theoretical models and experimental data from gene regulatory systems. We regularly find incongruence, in that fitness landscapes tend to comprise more peaks than their underlying genotype-phenotype landscapes. However, using evolutionary simulations, we show that this increased ruggedness need not impede adaptation.

Environmental selection and epistasis in an empirical phenotype-environment-fitness landscape

Fitness landscapes, mappings of genotype/phenotype to their effects on fitness, are invaluable concepts in evolutionary biochemistry. Although widely discussed, measurements of phenotype-fitness landscapes in proteins remain scarce. Here, we quantify all single mutational effects on fitness and phenotype (EC₅₀) of VIM-2 β-lactamase across a 64-fold range of ampicillin concentrations. We then construct a phenotype-fitness landscape that takes variations in environmental selection pressure into account. We found that a simple, empirical landscape accurately models the ~39,000 mutational data points, suggesting that the evolution of VIM-2 can be predicted on the basis of the selection environment. Our landscape provides new quantitative knowledge on the evolution of the β-lactamases and proteins in general, particularly their evolutionary dynamics under subinhibitory antibiotic concentrations, as well as the mechanisms and environmental dependence of non-specific epistasis.

Emergence and propagation of epistasis in metabolic networks

Epistasis is often used to probe functional relationships between genes, and it plays an important role in evolution. However, we lack theory to understand how functional relationships at the molecular level translate into epistasis at the level of whole-organism phenotypes, such as fitness. Here, I derive two rules for how epistasis between mutations with small effects propagates from lower- to higher-level phenotypes in a hierarchical metabolic network with first-order kinetics and how such epistasis depends on topology. Most importantly, weak epistasis at a lower level may be distorted as it propagates to higher levels. Computational analyses show that epistasis in more realistic models likely follows similar, albeit more complex, patterns. These results suggest that pairwise inter-gene epistasis should be common, and it should generically depend on the genetic background and environment. Furthermore, the epistasis coefficients measured for high-level phenotypes may not be sufficient to fully infer the underlying functional relationships.

Fitness variation across subtle environmental perturbations reveals local modularity and global pleiotropy of adaptation

Building a genotype-phenotype-fitness map of adaptation is a central goal in evolutionary biology. It is difficult even when adaptive mutations are known because it is hard to enumerate which phenotypes make these mutations adaptive. We address this problem by first quantifying how the fitness of hundreds of adaptive yeast mutants responds to subtle environmental shifts. We then model the number of phenotypes these mutations collectively influence by decomposing these patterns of fitness variation. We find that a small number of inferred phenotypes can predict fitness of the adaptive mutations near their original glucose-limited evolution condition. Importantly, inferred phenotypes that matter little to fitness at or near the evolution condition can matter strongly in distant environments. This suggests that adaptive mutations are locally modular - affecting a small number of phenotypes that matter to fitness in the environment where they evolved - yet globally pleiotropic - affecting additional phenotypes that may reduce or improve fitness in new environments.

Epistatic evidence for gender-dependant slow neurotransmission signalling in substance use disorders: PPP1R12B versus PPP1R1B

BACKGROUND

Slow neurotransmission including DARPP-32 signalling is implicated in substance use disorders (SUDs) by experimental systems but not yet in the human aetiology. PPP1R12B, encoding another protein in the DARPP-32 family, hasn't been studied in the brain.

METHODS

Brain-regional gene activity was assessed in three different animal models of SUDs for mRNA level alterations. Genetic associations were assessed by meta-analysis of pre-existing dbGaP GWAS datasets for main effects and epistasis with known genetic risks, followed by cell type-specific pathway delineation. Parkinson's disease (PD) was included as a dopamine-related disease control for SUDs.

FINDINGS

In animal models of SUDs, environmentally-altered PPP1R12B expression sex-dependently involves motivation-related brain regions. In humans with polysubstance abuse, meta-analysis of pre-existing datasets revealed that PPP1R12B and PPP1R1B, although expressed in dopamine vs. dopamine-recipient neurons, exerted similar interactions with known genetic risks such as ACTR1B and DRD2 in men but with ADH1B, HGFAC and DRD3 in women. These interactions reached genome-wide significances (P_meta<10^-20) for SUDs but not for PD (disease selectivity: P = 4.8 × 10^-142, OR = 6.7 for PPP1R12B; P = 8.0 × 10^-8, OR = 2.1 for PPP1R1B). CADM2 was the common risk in the molecular signalling regardless of gender and cell type.

INTERPRETATION

Gender-dependant slow neurotransmission may convey both genetic and environmental vulnerabilities selectively to SUDs.

FUNDING

Grants from National Institute on Drug Abuse (NIDA) and National Institute on Alcohol Abuse and Alcoholism (NIAAA) of U.S.A. and National Natural Science Foundation of China (NSFC).