A synthetic differentiation circuit in Escherichia coli for suppressing mutant takeover.

Differentiation is crucial for multicellularity. However, it is inherently susceptible to mutant cells that fail to differentiate. These mutants outcompete normal cells by excessive self-renewal. It remains unclear what mechanisms can resist such mutant expansion. Here, we demonstrate a solution by engineering a synthetic differentiation circuit in Escherichia coli that selects against these mutants via a biphasic fitness strategy. The circuit provides tunable production of synthetic analogs of stem, progenitor, and differentiated cells. It resists mutations by coupling differentiation to the production of an essential enzyme, thereby disadvantaging non-differentiating mutants. The circuit selected for and maintained a positive differentiation rate in long-term evolution. Surprisingly, this rate remained constant across vast changes in growth conditions. We found that transit-amplifying cells (fast-growing progenitors) underlie this environmental robustness. Our results provide insight into the stability of differentiation and demonstrate a powerful method for engineering evolutionarily stable multicellular consortia.

Dominance and multi-locus interaction.

Dominance is usually considered a constant value that describes the relative difference in fitness or phenotype between heterozygotes and the average of homozygotes at a focal polymorphic locus. However, the observed dominance can vary with the genetic background of the focal locus. Here, alleles at other loci modify the observed phenotype through position effects or dominance modifiers that are sometimes associated with pathogen resistance, lineage, sex, or mating type. Theoretical models have illustrated how variable dominance appears in the context of multi-locus interaction (epistasis). Here, we review empirical evidence for variable dominance and how the observed patterns may be captured by proposed epistatic models. We highlight how integrating epistasis and dominance is crucial for comprehensively understanding adaptation and speciation.

Cryo-EM structure and functional landscape of an RNA polymerase ribozyme.

The emergence of an RNA replicase capable of self-replication is considered an important stage in the origin of life. RNA polymerase ribozymes (PR) - including a variant that uses trinucleotide triphosphates (triplets) as substrates - have been created by in vitro evolution and are the closest functional analogues of the replicase, but the structural basis for their function is poorly understood. Here we use single-particle cryogenic electron microscopy (cryo-EM) and high-throughput mutation analysis to obtain the structure of a triplet polymerase ribozyme (TPR) apoenzyme and map its functional landscape. The cryo-EM structure at 5-Å resolution reveals the TPR as an RNA heterodimer comprising a catalytic subunit and a noncatalytic, auxiliary subunit, resembling the shape of a left hand with thumb and fingers at a 70° angle. The two subunits are connected by two distinct kissing-loop (KL) interactions that are essential for polymerase function. Our combined structural and functional data suggest a model for templated RNA synthesis by the TPR holoenzyme, whereby heterodimer formation and KL interactions preorganize the TPR for optimal primer-template duplex binding, triplet substrate discrimination, and templated RNA synthesis. These results provide a better understanding of TPR structure and function and should aid the engineering of more efficient PRs.

Biophysical principles predict fitness of SARS-CoV-2 variants.

SARS-CoV-2 employs its spike protein's receptor binding domain (RBD) to enter host cells. The RBD is constantly subjected to immune responses, while requiring efficient binding to host cell receptors for successful infection. However, our understanding of how RBD's biophysical properties contribute to SARS-CoV-2's epidemiological fitness remains largely incomplete. Through a comprehensive approach, comprising large-scale sequence analysis of SARS-CoV-2 variants and the identification of a fitness function based on binding thermodynamics, we unravel the relationship between the biophysical properties of RBD variants and their contribution to viral fitness. We developed a biophysical model that uses statistical mechanics to map the molecular phenotype space, characterized by dissociation constants of RBD to ACE2, LY-CoV016, LY-CoV555, REGN10987, and S309, onto an epistatic fitness landscape. We validate our findings through experimentally measured and machine learning (ML) estimated binding affinities, coupled with infectivity data derived from population-level sequencing. Our analysis reveals that this model effectively predicts the fitness of novel RBD variants and can account for the epistatic interactions among mutations, including explaining the later reversal of Q493R. Our study sheds light on the impact of specific mutations on viral fitness and delivers a tool for predicting the future epidemiological trajectory of previously unseen or emerging low-frequency variants. These insights offer not only greater understanding of viral evolution but also potentially aid in guiding public health decisions in the battle against COVID-19 and future pandemics.

In vivo functional phenotypes from a computational epistatic model of evolution.

Computational models of evolution are valuable for understanding the dynamics of sequence variation, to infer phylogenetic relationships or potential evolutionary pathways and for biomedical and industrial applications. Despite these benefits, few have validated their propensities to generate outputs with in vivo functionality, which would enhance their value as accurate and interpretable evolutionary algorithms. We demonstrate the power of epistasis inferred from natural protein families to evolve sequence variants in an algorithm we developed called sequence evolution with epistatic contributions (SEEC). Utilizing the Hamiltonian of the joint probability of sequences in the family as fitness metric, we sampled and experimentally tested for in vivo [Formula: see text]-lactamase activity in Escherichia coli TEM-1 variants. These evolved proteins can have dozens of mutations dispersed across the structure while preserving sites essential for both catalysis and interactions. Remarkably, these variants retain family-like functionality while being more active than their wild-type predecessor. We found that depending on the inference method used to generate the epistatic constraints, different parameters simulate diverse selection strengths. Under weaker selection, local Hamiltonian fluctuations reliably predict relative changes to variant fitness, recapitulating neutral evolution. SEEC has the potential to explore the dynamics of neofunctionalization, characterize viral fitness landscapes, and facilitate vaccine development.

Genotypic-phenotypic landscape computation based on first principle and deep learning.

The relationship between genotype and fitness is fundamental to evolution, but quantitatively mapping genotypes to fitness has remained challenging. We propose the Phenotypic-Embedding theorem (P-E theorem) that bridges genotype-phenotype through an encoder-decoder deep learning framework. Inspired by this, we proposed a more general first principle for correlating genotype-phenotype, and the P-E theorem provides a computable basis for the application of first principle. As an application example of the P-E theorem, we developed the Co-attention based Transformer model to bridge Genotype and Fitness model, a Transformer-based pre-train foundation model with downstream supervised fine-tuning that can accurately simulate the neutral evolution of viruses and predict immune escape mutations. Accordingly, following the calculation path of the P-E theorem, we accurately obtained the basic reproduction number (${R}_0$) of SARS-CoV-2 from first principles, quantitatively linked immune escape to viral fitness and plotted the genotype-fitness landscape. The theoretical system we established provides a general and interpretable method to construct genotype-phenotype landscapes, providing a new paradigm for studying theoretical and computational biology.

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.

Master regulators of biological systems in higher dimensions.

A longstanding goal of biology is to identify the key genes and species that critically impact evolution, ecology, and health. Network analysis has revealed keystone species that regulate ecosystems and master regulators that regulate cellular genetic networks. Yet these studies have focused on pairwise biological interactions, which can be affected by the context of genetic background and other species present, generating higher-order interactions. The important regulators of higher-order interactions are unstudied. To address this, we applied a high-dimensional geometry approach that quantifies epistasis in a fitness landscape to ask how individual genes and species influence the interactions in the rest of the biological network. We then generated and also reanalyzed 5-dimensional datasets (two genetic, two microbiome). We identified key genes (e.g., the rbs locus and pykF) and species (e.g., Lactobacilli) that control the interactions of many other genes and species. These higher-order master regulators can induce or suppress evolutionary and ecological diversification by controlling the topography of the fitness landscape. Thus, we provide a method and mathematical justification for exploration of biological networks in higher dimensions.

Fluctuating selection maintains distinct species phenotypes in an ecological community in the wild.

Species' phenotypic characteristics often remain unchanged over long stretches of geological time. Stabilizing selection-in which fitness is highest for intermediate phenotypes and lowest for the extremes-has been widely invoked as responsible for this pattern. At the community level, such stabilizing selection acting individually on co-occurring species is expected to produce a rugged fitness landscape on which different species occupy distinct fitness peaks. However, even with an explosion of microevolutionary field studies over the past four decades, evidence for persistent stabilizing selection driving long-term stasis is lacking. Nonetheless, biologists continue to invoke stabilizing selection as a major factor explaining macroevolutionary patterns. Here, by directly measuring natural selection in the wild, we identified a complex community-wide fitness surface in which four Anolis lizard species each occupy a distinct fitness peak close to their mean phenotype. The presence of local fitness optima within species, and fitness valleys between species, presents a barrier to adaptive evolutionary change and acts to maintain species differences through time. However, instead of continuously operating stabilizing selection, we found that species were maintained on these peaks by the combination of many independent periods among which selection fluctuated in form, strength, direction, or existence and in which stabilizing selection rarely occurred. Our results suggest that lack of substantial phenotypic evolutionary change through time may be the result of selection, but not persistent stabilizing selection as classically envisioned.

Increased Positive Selection in Highly Recombining Genes Does not Necessarily Reflect an Evolutionary Advantage of Recombination.

It is commonly thought that the long-term advantage of meiotic recombination is to dissipate genetic linkage, allowing natural selection to act independently on different loci. It is thus theoretically expected that genes with higher recombination rates evolve under more effective selection. On the other hand, recombination is often associated with GC-biased gene conversion (gBGC), which theoretically interferes with selection by promoting the fixation of deleterious GC alleles. To test these predictions, several studies assessed whether selection was more effective in highly recombining genes (due to dissipation of genetic linkage) or less effective (due to gBGC), assuming a fixed distribution of fitness effects (DFE) for all genes. In this study, I directly derive the DFE from a gene's evolutionary history (shaped by mutation, selection, drift, and gBGC) under empirical fitness landscapes. I show that genes that have experienced high levels of gBGC are less fit and thus have more opportunities for beneficial mutations. Only a small decrease in the genome-wide intensity of gBGC leads to the fixation of these beneficial mutations, particularly in highly recombining genes. This results in increased positive selection in highly recombining genes that is not caused by more effective selection. Additionally, I show that the death of a recombination hotspot can lead to a higher dN/dS than its birth, but with substitution patterns biased towards AT, and only at selected positions. This shows that controlling for a substitution bias towards GC is therefore not sufficient to rule out the contribution of gBGC to signatures of accelerated evolution. Finally, although gBGC does not affect the fixation probability of GC-conservative mutations, I show that by altering the DFE, gBGC can also significantly affect nonsynonymous GC-conservative substitution patterns.

Discovering pathways through ribozyme fitness landscapes using information theoretic quantification of epistasis.

The identification of catalytic RNAs is typically achieved through primarily experimental means. However, only a small fraction of sequence space can be analyzed even with high-throughput techniques. Methods to extrapolate from a limited data set to predict additional ribozyme sequences, particularly in a human-interpretable fashion, could be useful both for designing new functional RNAs and for generating greater understanding about a ribozyme fitness landscape. Using information theory, we express the effects of epistasis (i.e., deviations from additivity) on a ribozyme. This representation was incorporated into a simple model of the epistatic fitness landscape, which identified potentially exploitable combinations of mutations. We used this model to theoretically predict mutants of high activity for a self-aminoacylating ribozyme, identifying potentially active triple and quadruple mutants beyond the experimental data set of single and double mutants. The predictions were validated experimentally, with nine out of nine sequences being accurately predicted to have high activity. This set of sequences included mutants that form a previously unknown evolutionary "bridge" between two ribozyme families that share a common motif. Individual steps in the method could be examined, understood, and guided by a human, combining interpretability and performance in a simple model to predict ribozyme sequences by extrapolation.

Fitness and Functional Landscapes of the E. coli RNase III Gene rnc.

How protein properties such as protein activity and protein essentiality affect the distribution of fitness effects (DFE) of mutations are important questions in protein evolution. Deep mutational scanning studies typically measure the effects of a comprehensive set of mutations on either protein activity or fitness. Our understanding of the underpinnings of the DFE would be enhanced by a comprehensive study of both for the same gene. Here, we compared the fitness effects and in vivo protein activity effects of ∼4,500 missense mutations in the E. coli rnc gene. This gene encodes RNase III, a global regulator enzyme that cleaves diverse RNA substrates including precursor ribosomal RNA and various mRNAs including its own 5' untranslated region (5'UTR). We find that RNase III's ability to cleave dsRNA is the most important determinant of the fitness effects of rnc mutations. The DFE of RNase III was bimodal, with mutations centered around neutral and deleterious effects, consistent with previously reported DFE's of enzymes with a singular physiological role. Fitness was buffered to small effects on RNase III activity. The enzyme's RNase III domain, which contains the RNase III signature motif and all active site residues, was more sensitive to mutation than its dsRNA binding domain, which is responsible for recognition and binding to dsRNA. Differential effects on fitness and functional scores for mutations at highly conserved residues G97, G99, and F188 suggest that these positions may be important for RNase III cleavage specificity.

Genes Vary Greatly in Their Propensity for Collateral Fitness Effects of Mutations.

Mutations can have deleterious fitness effects when they decrease protein specific activity or decrease active protein abundance. Mutations will also be deleterious when they cause misfolding or misinteractions that are toxic to the cell (i.e., independent of whether the mutations affect specific activity and abundance). The extent to which protein evolution is shaped by these and other collateral fitness effects is unclear in part because little is known of their frequency and magnitude. Using deep mutational scanning (DMS), we previously found at least 42% of missense mutations in the TEM-1 ß-lactamase antibiotic resistance gene cause deleterious collateral fitness effects. Here, we used DMS to comprehensively determine the collateral fitness effects of missense mutations in three genes encoding the antibiotic resistance proteins New Delhi metallo-ß-lactamase (NDM-1), chloramphenicol acetyltransferase I (CAT-I), and 2″-aminoglycoside nucleotidyltransferase (AadB). AadB (20%), CAT-I (0.9%), and NDM-1 (0.2%) were less susceptible to deleterious collateral fitness effects than TEM-1 (42%) indicating that genes have different propensities for these effects. As was observed with TEM-1, all the studied deleterious aadB mutants increased aggregation. However, aggregation did not correlate with collateral fitness effects for many of the deleterious mutants of CAT-I and NDM-1. Select deleterious mutants caused unexpected phenotypes to emerge. The introduction of internal start codons in CAT-1 caused loss of the episome and a mutation in aadB made its cognate antibiotic essential for growth. Our study illustrates how the complexity of the cell provides a rich environment for collateral fitness effects and new phenotypes to emerge.

Sampling Strategies for Experimentally Mapping Molecular Fitness Landscapes Using High-Throughput Methods.

Empirical studies of genotype-phenotype-fitness maps of proteins are fundamental to understanding the evolutionary process, in elucidating the space of possible genotypes accessible through mutations in a landscape of phenotypes and fitness effects. Yet, comprehensively mapping molecular fitness landscapes remains challenging since all possible combinations of amino acid substitutions for even a few protein sites are encoded by an enormous genotype space. High-throughput mapping of genotype space can be achieved using large-scale screening experiments known as multiplexed assays of variant effect (MAVEs). However, to accommodate such multi-mutational studies, the size of MAVEs has grown to the point where a priori determination of sampling requirements is needed. To address this problem, we propose calculations and simulation methods to approximate minimum sampling requirements for multi-mutational MAVEs, which we combine with a new library construction protocol to experimentally validate our approximation approaches. Analysis of our simulated data reveals how sampling trajectories differ between simulations of nucleotide versus amino acid variants and among mutagenesis schemes. For this, we show quantitatively that marginal gains in sampling efficiency demand increasingly greater sampling effort when sampling for nucleotide sequences over their encoded amino acid equivalents. We present a new library construction protocol that efficiently maximizes sequence variation, and demonstrate using ultradeep sequencing that the library encodes virtually all possible combinations of mutations within the experimental design. Insights learned from our analyses together with the methodological advances reported herein are immediately applicable toward pooled experimental screens of arbitrary design, enabling further assay upscaling and expanded testing of genotype space.

A Long-Lived Alpine Perennial Advances Flowering under Warmer Conditions but Not Enough to Maintain Reproductive Success.

AbstractAssessing whether phenological shifts in response to climate change confer a fitness advantage requires investigating the relationships among phenology, fitness, and environmental drivers of selection. Despite widely documented advancements in phenology with warming climate, we lack empirical estimates of how selection on phenology varies in response to continuous climate drivers or how phenological shifts in response to warming conditions affect fitness. We leverage an unusual long-term dataset with repeated, individual measurements of phenology and reproduction in a long-lived alpine plant. We analyze phenotypic plasticity in flowering phenology in relation to two climate drivers, snowmelt timing and growing degree days (GDDs). Plants flower earlier with increased GDDs and earlier snowmelt, and directional selection also favors earlier flowering under these conditions. However, reproduction still declines with warming and early snowmelt, even when flowering is early. Furthermore, the steepness of this reproductive decline increases dramatically with warming conditions, resulting in very little fruit production regardless of flowering time once GDDs exceed approximately 225 degree days or snowmelt occurs before May 15. Even though advancing phenology confers a fitness advantage relative to stasis, these shifts are insufficient to maintain reproduction under warming, highlighting limits to the potential benefits of phenological plasticity under climate change.

Mutations in yeast are deleterious on average regardless of the degree of adaptation to the testing environment.

The role of spontaneous mutations in evolution depends on the distribution of their effects on fitness. Despite a general consensus that new mutations are deleterious on average, a handful of mutation accumulation experiments in diverse organisms instead suggest that beneficial and deleterious mutations can have comparable fitness impacts, i.e. the product of their respective rates and effects can be roughly equal. We currently lack a general framework for predicting when such a pattern will occur. One idea is that beneficial mutations will be more evident in genotypes that are not well adapted to the testing environment. We tested this prediction experimentally in the laboratory yeast Saccharomyces cerevisiae by allowing nine replicate populations to adapt to novel environments with complex sets of stressors. After >1000 asexual generations interspersed with 41 rounds of sexual reproduction, we assessed the mean effect of induced mutations on yeast growth in both the environment to which they had been adapting and the alternative novel environment. The mutations were deleterious on average, with the severity depending on the testing environment. However, we found no evidence that the adaptive match between genotype and environment is predictive of mutational fitness effects.

Elevated [CO2] and temperature augment gas exchange and shift the fitness landscape in a montane forb.

Climate change is simultaneously increasing carbon dioxide concentrations ([CO2]) and temperature. These factors could interact to influence plant physiology and performance. Alternatively, increased [CO2] may offset costs associated with elevated temperatures. Furthermore, the interaction between elevated temperature and [CO2] may differentially affect populations from along an elevational gradient and disrupt local adaptation. We conducted a multifactorial growth chamber experiment to examine the interactive effects of temperature and [CO2] on fitness and ecophysiology of diverse accessions of Boechera stricta (Brassicaceae) sourced from a broad elevational gradient in Colorado. We tested whether increased [CO2] would enhance photosynthesis across accessions, and whether warmer conditions would depress the fitness of high-elevation accessions owing to steep reductions in temperature with increasing elevation in this system. Elevational clines in [CO2] are not as evident, making it challenging to predict how locally adapted ecotypes will respond to elevated [CO2]. This experiment revealed that elevated [CO2] increased photosynthesis and intrinsic water use efficiency across all accessions. However, these instantaneous responses to treatments did not translate to changes in fitness. Instead, increased temperatures reduced the probability of reproduction for all accessions. Elevated [CO2] and increased temperatures interacted to shift the adaptive landscape, favoring lower elevation accessions for the probability of survival and fecundity. Our results suggest that elevated temperatures and [CO2] associated with climate change could have severe negative consequences, especially for high-elevation populations.

Adenovirus-vectored vaccine containing multidimensionally conserved parts of the HIV proteome is immunogenic in rhesus macaques.

An effective vaccine that can protect against HIV infection does not exist. A major reason why a vaccine is not available is the high mutability of the virus, which enables it to evolve mutations that can evade human immune responses. This challenge is exacerbated by the ability of the virus to evolve compensatory mutations that can partially restore the fitness cost of immune-evading mutations. Based on the fitness landscapes of HIV proteins that account for the effects of coupled mutations, we designed a single long peptide immunogen comprising parts of the HIV proteome wherein mutations are likely to be deleterious regardless of the sequence of the rest of the viral protein. This immunogen was then stably expressed in adenovirus vectors that are currently in clinical development. Macaques immunized with these vaccine constructs exhibited T-cell responses that were comparable in magnitude to animals immunized with adenovirus vectors with whole HIV protein inserts. Moreover, the T-cell responses in immunized macaques strongly targeted regions contained in our immunogen. These results suggest that further studies aimed toward using our vaccine construct for HIV prophylaxis and cure are warranted.

Adaptive evolution of hybrid bacteria by horizontal gene transfer.

Horizontal gene transfer (HGT) is an important factor in bacterial evolution that can act across species boundaries. Yet, we know little about rate and genomic targets of cross-lineage gene transfer and about its effects on the recipient organism's physiology and fitness. Here, we address these questions in a parallel evolution experiment with two Bacillus subtilis lineages of 7% sequence divergence. We observe rapid evolution of hybrid organisms: gene transfer swaps ∼12% of the core genome in just 200 generations, and 60% of core genes are replaced in at least one population. By genomics, transcriptomics, fitness assays, and statistical modeling, we show that transfer generates adaptive evolution and functional alterations in hybrids. Specifically, our experiments reveal a strong, repeatable fitness increase of evolved populations in the stationary growth phase. By genomic analysis of the transfer statistics across replicate populations, we infer that selection on HGT has a broad genetic basis: 40% of the observed transfers are adaptive. At the level of functional gene networks, we find signatures of negative, positive, and epistatic selection, consistent with hybrid incompatibilities and adaptive evolution of network functions. Our results suggest that gene transfer navigates a complex cross-lineage fitness landscape, bridging epistatic barriers along multiple high-fitness paths.

The ecological causes of functional distinctiveness in communities.

Recent work has shown that evaluating functional trait distinctiveness, the average trait distance of a species to other species in a community offers promising insights into biodiversity dynamics and ecosystem functioning. However, the ecological mechanisms underlying the emergence and persistence of functionally distinct species are poorly understood. Here, we address the issue by considering a heterogeneous fitness landscape whereby functional dimensions encompass peaks representing trait combinations yielding positive population growth rates in a community. We identify four ecological cases contributing to the emergence and persistence of functionally distinct species. First, environmental heterogeneity or alternative phenotypic designs can drive positive population growth of functionally distinct species. Second, sink populations with negative population growth can deviate from local fitness peaks and be functionally distinct. Third, species found at the margin of the fitness landscape can persist but be functionally distinct. Fourth, biotic interactions (positive or negative) can dynamically alter the fitness landscape. We offer examples of these four cases and guidelines to distinguish between them. In addition to these deterministic processes, we explore how stochastic dispersal limitation can yield functional distinctiveness. Our framework offers a novel perspective on the relationship between fitness landscape heterogeneity and the functional composition of ecological assemblages.