Fluctuation domains in adaptive evolution

We derive an expression for the variation between parallel trajectories in phenotypic evolution, extending the well known result that predicts the mean evolutionary path in adaptive dynamics or quantitative genetics. We show how this expression gives rise to the notion of fluctuation domains-parts of the fitness landscape where the rate of evolution is very predictable (due to fluctuation dissipation) and parts where it is highly variable (due to fluctuation enhancement). These fluctuation domains are determined by the curvature of the fitness landscape. Regions of the fitness landscape with positive curvature, such as adaptive valleys or branching points, experience enhancement. Regions with negative curvature, such as adaptive peaks, experience dissipation. We explore these dynamics in the ecological scenarios of implicit and explicit competition for a limiting resource.     

Evidence for a multi‐peak adaptive landscape in the evolution of trophic morphology in terapontid fishes

Using the Australian marine‐freshwater terapontid fishes as a model system, we examined the role of dietary phenotypic optima in an adaptive macro‐evolutionary landscape. Comparative modelling relying on both a priori and data‐driven identification of selective regimes suggested multi‐peak models as best describing much of the dietary phenotypic landscape of terapontids. Both approaches identified common phenotypic optima for different lineages of marine and freshwater herbivores, and minimal differentiation between carnivores and omnivores, irrespective of their phylogenetic relationships, as the model best describing morphological evolution. Significant correlations also existed between these phenotypic axes and proportions of non‐animal dietary items in species’ diets. While simulation results provided evidence for a multi‐peak adaptive landscape in the evolution of trophic morphology in terapontids, they could not rule out chance convergence in these adaptive peaks. However, they do provide scope for identifying areas for more detailed, functionally specific study of phenotypic convergence in herbivorous terapontid trophic habits. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 623–634.     

Fitness trade-offs in the evolution of dihydrofolate reductase and drug resistance in Plasmodium falciparum

Background

Patterns of emerging drug resistance reflect the underlying adaptive landscapes for specific drugs. In Plasmodium falciparum, the parasite that causes the most serious form of malaria, antifolate drugs inhibit the function of essential enzymes in the folate pathway. However, a handful of mutations in the gene coding for one such enzyme, dihydrofolate reductase, confer drug resistance. Understanding how evolution proceeds from drug susceptibility to drug resistance is critical if new antifolate treatments are to have sustained usefulness.

Methodology/Principal Findings

We use a transgenic yeast expression system to build on previous studies that described the adaptive landscape for the antifolate drug pyrimethamine, and we describe the most likely evolutionary trajectories for the evolution of drug resistance to the antifolate chlorcycloguanil. We find that the adaptive landscape for chlorcycloguanil is multi-peaked, not all highly resistant alleles are equally accessible by evolution, and there are both commonalities and differences in adaptive landscapes for chlorcycloguanil and pyrimethamine.

Conclusions/Significance

Our findings suggest that cross-resistance between drugs targeting the same enzyme reflect the fitness landscapes associated with each particular drug and the position of the genotype on both landscapes. The possible public health implications of these findings are discussed.     

On Simplified Global Nonlinear Function for Fitness Landscape: A Case Study of Inverse Protein Folding

The construction of fitness landscape has broad implication in understanding molecular evolution, cellular epigenetic state, and protein structures. We studied the problem of constructing fitness landscape of inverse protein folding or protein design, with the aim to generate amino acid sequences that would fold into an a priori determined structural fold which would enable engineering novel or enhanced biochemistry. For this task, an effective fitness function should allow identification of correct sequences that would fold into the desired structure. In this study, we showed that nonlinear fitness function for protein design can be constructed using a rectangular kernel with a basis set of proteins and decoys chosen a priori. The full landscape for a large number of protein folds can be captured using only 480 native proteins and 3,200 non-protein decoys via a finite Newton method. A blind test of a simplified version of fitness function for sequence design was carried out to discriminate simultaneously 428 native sequences not homologous to any training proteins from 11 million challenging protein-like decoys. This simplified function correctly classified 408 native sequences (20 misclassifications, 95% correct rate), which outperforms several other statistical linear scoring function and optimized linear function. Our results further suggested that for the task of global sequence design of 428 selected proteins, the search space of protein shape and sequence can be effectively parametrized with just about 3,680 carefully chosen basis set of proteins and decoys, and we showed in addition that the overall landscape is not overly sensitive to the specific choice of this set. Our results can be generalized to construct other types of fitness landscape.     

Selection and evolution of bacteriophages in cellstat

Objectives of this work were as follows: 1. to establish a laboratory experimental system utilizable in a biophysical approach to molecular evolution: and 2. to provide real world parameters to theories of molecular evolution, especially to eigen's theory of quasi-species.Secretion type bacteriophage fd of E. coli, closely related phages and artificial chimera phages of fd, and a virulent phage Qβ of E. coli were cultured continuously in a specially designed fermenter called a “cellstat”. A phage is cultured in a flow of host bacterial cells. Due to its high dilution rate, the mutant cell could not be selected in the cellstat. It was therefore recognized that the cellstat is suitable for study of the selection and evolution process of a bacteriophage under well-defined environmental conditions without interference from host cell mutations.Population dynamics of bacteriophages of various types in the cellstat were studied theoretically by computer simulation and experimentally. A genetically invariable pure population of phage behaves like an open non-linear chemical reaction system. An invariable mixed population shows a selection process, while a variable population generates an evolution process.Kinetic constants describing the dynamics were determined by curve fitting between the theoretical and the experimental curve obtained from competition experiments and from biological relaxation experiments. One of the most important kinetic parameters thus obtained was the selection coefficient, and its dependence on the base sequence of phage DNA. We drew a local landscape of the selection coefficient near the fd sequence on the base sequence space. From this landscape we were able to confirm the importance of slightly deleterious mutants in molecular evolution. We also confirmed the possibility of developing an evolutionary molecular engineering using a cellstat as an evolution reactor and fd phage as a working replicon.Novelties of this work were as follows: 1. the first stable continuous culture of a bacteriophage was achieved with a cellstat; 2. a local landscape of selection coefficient near the fd sequence on the sequence space was the first experimental drawing of such a map; 3. a biological relaxation method was realized to measure kinetic constants of a biological kinetic process, or molecular evolution; and 4. a practical engineering process of evolutionary molecular engineering was proposed.     

Biomolecular information gained through in vitro evolution

An in vitro evolution is a simplified Darwinian evolution in well-controlled surroundings. This evolution process can be modeled as a hill-climbing or adaptive walk on a fitness landscape in sequence space. The evolving molecular system gains at least two kinds of information originating from the converged sequences and the fitness increment of the evolving biopolymer as the adaptive walker. These two represent two aspects of the biomolecular information, its extent and its content, respectively. Here, we review studies related to formulation of the “content” and “extent” of biomolecular information. The two aspects are interconnected through physicochemical properties of the biopolymer, contrary to the case of conventional information, which seems to be independent of matter. The interconnection was analyzed based on the analogy between the evolution process and thermodynamics. The linear combination of the two by a temperature-like fluctuation factor resulted in a free-energy-like monotonically increasing function during the evolution process.     

Molecular evolution of rbcL in three gymnosperm families: identifying adaptive and coevolutionary patterns

Background

The chloroplast-localized ribulose-1, 5-biphosphate carboxylase/oxygenase (Rubisco), the primary enzyme responsible for autotrophy, is instrumental in the continual adaptation of plants to variations in the concentrations of CO₂. The large subunit (LSU) of Rubisco is encoded by the chloroplast rbcL gene. Although adaptive processes have been previously identified at this gene, characterizing the relationships between the mutational dynamics at the protein level may yield clues on the biological meaning of such adaptive processes. The role of such coevolutionary dynamics in the continual fine-tuning of RbcL remains obscure.

Results

We used the timescale and phylogenetic analyses to investigate and search for processes of adaptive evolution in rbcL gene in three gymnosperm families, namely Podocarpaceae, Taxaceae and Cephalotaxaceae. To understand the relationships between regions identified as having evolved under adaptive evolution, we performed coevolutionary analyses using the software CAPS. Importantly, adaptive processes were identified at amino acid sites located on the contact regions among the Rubisco subunits and on the interface between Rubisco and its activase. Adaptive amino acid replacements at these regions may have optimized the holoenzyme activity. This hypothesis was pinpointed by evidence originated from our analysis of coevolution that supported the correlated evolution between Rubisco and its activase. Interestingly, the correlated adaptive processes between both these proteins have paralleled the geological variation history of the concentration of atmospheric CO₂.

Conclusions

The gene rbcL has experienced bursts of adaptations in response to the changing concentration of CO₂ in the atmosphere. These adaptations have emerged as a result of a continuous dynamic of mutations, many of which may have involved innovation of functional Rubisco features. Analysis of the protein structure and the functional implications of such mutations put forward the conclusion that this evolutionary scenario has been possible through a complex interplay between adaptive mutations, often structurally destabilizing, and compensatory mutations. Our results unearth patterns of evolution that have likely optimized the Rubisco activity and uncover mutational dynamics useful in the molecular engineering of enzymatic activities.

Reviewers

This article was reviewed by Prof. Christian Blouin (nominated by Dr W Ford Doolittle), Dr Endre Barta (nominated by Dr Sandor Pongor), and Dr Nicolas Galtier.     

Mainstreaming Caenorhabditis elegans in experimental evolution

Experimental evolution provides a powerful manipulative tool for probing evolutionary process and mechanism. As this approach to hypothesis testing has taken purchase in biology, so too has the number of experimental systems that use it, each with its own unique strengths and weaknesses. The depth of biological knowledge about Caenorhabditis nematodes, combined with their laboratory tractability, positions them well for exploiting experimental evolution in animal systems to understand deep questions in evolution and ecology, as well as in molecular genetics and systems biology. To date, Caenorhabditis elegans and related species have proved themselves in experimental evolution studies of the process of mutation, host–pathogen coevolution, mating system evolution and life-history theory. Yet these organisms are not broadly recognized for their utility for evolution experiments and remain underexploited. Here, we outline this experimental evolution work undertaken so far in Caenorhabditis, detail simple methodological tricks that can be exploited and identify research areas that are ripe for future discovery.     

Altered retinoic acid signalling underpins dentition evolution

Small variations in signalling pathways have been linked to phenotypic diversity and speciation. In vertebrates, teeth represent a reservoir of adaptive morphological structures that are prone to evolutionary change. Cyprinid fish display an impressive diversity in tooth number, but the signals that generate such diversity are unknown. Here, we show that retinoic acid (RA) availability influences tooth number size in Cyprinids. Heterozygous adult zebrafish heterozygous for the cyp26b1 mutant that encodes an enzyme able to degrade RA possess an extra tooth in the ventral row. Expression analysis of pharyngeal mesenchyme markers such as dlx2a and lhx6 shows lateral, anterior and dorsal expansion of these markers in RA-treated embryos, whereas the expression of the dental epithelium markers dlx2b and dlx3b is unchanged. Our analysis suggests that changes in RA signalling play an important role in the diversification of teeth in Cyprinids. Our work illustrates that through subtle changes in the expression of rate-limiting enzymes, the RA pathway is an active player of tooth evolution in fish.     

Adaptive Landscape by Environment Interactions Dictate Evolutionary Dynamics in Models of Drug Resistance

The adaptive landscape analogy has found practical use in recent years, as many have explored how their understanding can inform therapeutic strategies that subvert the evolution of drug resistance. A major barrier to applications of these concepts is a lack of detail concerning how the environment affects adaptive landscape topography, and consequently, the outcome of drug treatment. Here we combine empirical data, evolutionary theory, and computer simulations towards dissecting adaptive landscape by environment interactions for the evolution of drug resistance in two dimensions—drug concentration and drug type. We do so by studying the resistance mediated by Plasmodium falciparum dihydrofolate reductase (DHFR) to two related inhibitors—pyrimethamine and cycloguanil—across a breadth of drug concentrations. We first examine whether the adaptive landscapes for the two drugs are consistent with common definitions of cross-resistance. We then reconstruct all accessible pathways across the landscape, observing how their structure changes with drug environment. We offer a mechanism for non-linearity in the topography of accessible pathways by calculating of the interaction between mutation effects and drug environment, which reveals rampant patterns of epistasis. We then simulate evolution in several different drug environments to observe how these individual mutation effects (and patterns of epistasis) influence paths taken at evolutionary “forks in the road” that dictate adaptive dynamics in silico. In doing so, we reveal how classic metrics like the IC₅₀ and minimal inhibitory concentration (MIC) are dubious proxies for understanding how evolution will occur across drug environments. We also consider how the findings reveal ambiguities in the cross-resistance concept, as subtle differences in adaptive landscape topography between otherwise equivalent drugs can drive drastically different evolutionary outcomes. Summarizing, we discuss the results with regards to their basic contribution to the study of empirical adaptive landscapes, and in terms of how they inform new models for the evolution of drug resistance.     

Environmental correlates of phenotypic evolution in ecologically diverse Liolaemus lizards

Evolutionary correlations between phenotypic and environmental traits characterize adaptive radiations. However, the lizard genus Liolaemus, one of the most ecologically diverse terrestrial vertebrate radiations on earth, has so far shown limited or mixed evidence of adaptive diversification in phenotype. Restricted use of comprehensive environmental data, incomplete taxonomic representation and not considering phylogenetic uncertainty may have led to contradictory evidence. We compiled a 26‐taxon dataset for the Liolaemus gracilis species group, representing much of the ecological diversity represented within Liolaemus and used environmental data to characterize how environments occupied by species'' relate to phenotypic evolution. Our analyses, explicitly accounting for phylogenetic uncertainty, suggest diversification in phenotypic traits toward the present, with body shape evolution rapidly evolving in this group. Body shape evolution correlates with the occupation of different structural habitats indicated by vegetation axes suggesting species have adapted for maximal locomotory performance in these habitats. Our results also imply that the effects of phylogenetic uncertainty and model misspecification may be more extensive on univariate, relative to multivariate analyses of evolutionary correlations, which is an important consideration in analyzing data from rapidly radiating adaptive radiations.     

Epigenetic and phenotypic variability in populations of Schistosoma mansoni– a possible kick‐off for adaptive host/parasite evolution

Epigenetics, the science of heritable but modifiable information, is now a well‐accepted component of many research fields. Nevertheless, epigenetics has not yet found broad appreciation in one of the most exciting fields of biology: the comprehension of evolution. This is surprising, since the reason for the existence of this alternative information‐transmitting system lies certainly in the evolutionary advantage it provides. Theoretical considerations support a model in which epigenetic mechanisms allow for increasing phenotypic variability and permit populations to explore the adaptive landscape without modifications of the genotype. The data presented here support the view that modulating the epigenotype of the human bloodfluke Schistosoma mansoni by treatment of larvae with histone deacetylase inhibitor leads indeed to an increase of phenotypic variability. It is therefore conceivable that environmentally induced changes in the epigenotype release new phenotypes on which selection can act and that this process is the first step in adaptive evolution.     

The effect of demographic population factors and individual biological parameters on the rate of the neutral molecular evolution

The dependence of the rate of neutral molecular evolution on both biological parameters and demographic population factors has been investigated. Real genetic data and an individual-based model of the population dynamics were used for the study. The first part of the study deals with tracing the neutral molecular evolution occurring in the model concurrently with adaptive speciation. The second part concerns the effect of individual biological parameters and demographic population factors of members of the Order Testudines (Turtles) on the relative rate of evolution of the mitochondrial CytB gene. It has been shown that demographic population factors and individual biological parameters affect the rate of the neutral molecular evolution.     

The economics of managing evolution

Humans are altering biological systems at unprecedented rates, and these alterations often have longer-term evolutionary impacts. Most obvious is the spread of resistance to pesticides and antibiotics. There are a wide variety of management strategies available to slow this evolution, and there are many reasons for using them. In this paper, we focus on the economic aspects of evolution management and ask: When is it economically beneficial for an individual decision-maker to invest in evolution management? We derive a simple dimensionless inequality showing that it is cost-effective to manage evolution when the percentage increase in the effective life span of the biological resource that management generates is larger than the percentage increase in annual profit that could be obtained by not managing evolution. We show how this inequality can be used to determine optimal investment choices for single decision-makers, to determine Nash equilibrium investment choices for multiple interacting decision-makers, and to examine how these equilibrium choices respond to regulatory interventions aimed at stimulating investment in evolution management. Our results are illustrated with examples involving Bacillus thuringiensis (Bt) crops and antibiotic use in fish farming.

Humans are altering biological systems at unprecedented rates and these alterations often have longer-term evolutionary impacts, such as the spread of resistance to pesticides and antibiotics. In this study, a simple mathematical criterion is derived determining when it is economically beneficial to invest in strategies for controlling and managing evolutionary change.     

Forecasting of phenotypic and genetic outcomes of experimental evolution in Pseudomonas protegens

Experimental evolution with microbes is often highly repeatable under identical conditions, suggesting the possibility to predict short-term evolution. However, it is not clear to what degree evolutionary forecasts can be extended to related species in non-identical environments, which would allow testing of general predictive models and fundamental biological assumptions. To develop an extended model system for evolutionary forecasting, we used previous data and models of the genotype-to-phenotype map from the wrinkly spreader system in Pseudomonas fluorescens SBW25 to make predictions of evolutionary outcomes on different biological levels for Pseudomonas protegens Pf-5. In addition to sequence divergence (78% amino acid and 81% nucleotide identity) for the genes targeted by mutations, these species also differ in the inability of Pf-5 to make cellulose, which is the main structural basis for the adaptive phenotype in SBW25. The experimental conditions were changed compared to the SBW25 system to test if forecasts were extendable to a non-identical environment. Forty-three mutants with increased ability to colonize the air-liquid interface were isolated, and the majority had reduced motility and was partly dependent on the Pel exopolysaccharide as a structural component. Most (38/43) mutations are expected to disrupt negative regulation of the same three diguanylate cyclases as in SBW25, with a smaller number of mutations in promoter regions, including an uncharacterized polysaccharide synthase operon. A mathematical model developed for SBW25 predicted the order of the three main pathways and the genes targeted by mutations, but differences in fitness between mutants and mutational biases also appear to influence outcomes. Mutated regions in proteins could be predicted in most cases (16/22), but parallelism at the nucleotide level was low and mutational hot spot sites were not conserved. This study demonstrates the potential of short-term evolutionary forecasting in experimental populations and provides testable predictions for evolutionary outcomes in other Pseudomonas species.     

The speciation of conger eel galectins by rapid adaptive evolution

Many cases of accelerated evolution driven by positive Darwinian selection are identified in the genes of venomous and reproductive proteins. This evolutional phenomenon might have important consequences in their gene-products' functions, such as multiple specific toxins for quick immobilization of the prey and the establishment of barriers to fertilization that might lead to speciation, and in the molecular evolution of novel genes. Recently, we analyzed the molecular evolution of two galectins isolated from the skin mucus of conger eel (Conger myriaster), named congerins I and II, by cDNA cloning and X-ray structural analysis, and we found that they have evolved in the rapid adaptive manner to emergence of a new structure including strand-swapping and a unique new ligand-binding site. In this review article we summarize and discuss the molecular evolution, especially the rapid adaptive evolution, and the structure-function relationships of conger eel galectins. Published in 2004.     

Adaptive landscape and functional diversity of Neotropical cichlids: implications for the ecology and evolution of Cichlinae (Cichlidae; Cichliformes)

Morphological, lineage and ecological diversity can vary substantially even among closely related lineages. Factors that influence morphological diversification, especially in functionally relevant traits, can help to explain the modern distribution of disparity across phylogenies and communities. Multivariate axes of feeding functional morphology from 75 species of Neotropical cichlid and a stepwise‐AIC algorithm were used to estimate the adaptive landscape of functional morphospace in Cichlinae. Adaptive landscape complexity and convergence, as well as the functional diversity of Cichlinae, were compared with expectations under null evolutionary models. Neotropical cichlid feeding function varied primarily between traits associated with ram feeding vs. suction feeding/biting and secondarily with oral jaw muscle size and pharyngeal crushing capacity. The number of changes in selective regimes and the amount of convergence between lineages was higher than expected under a null model of evolution, but convergence was not higher than expected under a similarly complex adaptive landscape. Functional disparity was compatible with an adaptive landscape model, whereas the distribution of evolutionary change through morphospace corresponded with a process of evolution towards a single adaptive peak. The continentally distributed Neotropical cichlids have evolved relatively rapidly towards a number of adaptive peaks in functional trait space. Selection in Cichlinae functional morphospace is more complex than expected under null evolutionary models. The complexity of selective constraints in feeding morphology has likely been a significant contributor to the diversity of feeding ecology in this clade.     

Glucosinolate structures in evolution

By 2000, around 106 natural glucosinolates (GSLs) were probably documented. In the past decade, 26 additional natural GSL structures have been elucidated and documented. Hence, the total number of documented GSLs from nature by 2011 can be estimated to around 132. A considerable number of additional suggested structures are concluded not to be sufficiently documented. In many cases, NMR spectroscopy would have provided the missing structural information. Of the GSLs documented in the past decade, several are of previously unexpected structures and occur at considerable levels. Most originate from just four species: Barbarea vulgaris, Arabidopsis thaliana, Eruca sativa and Isatis tinctoria. Acyl derivatives of known GSLs comprised 15 of the 26 newly documented structures, while the remaining exhibited new substitution patterns or chain length, or contained a mercapto group or related thio-functionality.GSL identification methods are reviewed, and the importance of using authentic references and structure-sensitive detection methods such as MS and NMR is stressed, especially when species with relatively unknown chemistry are analyzed. An example of qualitative GSL analysis is presented with experimental details (group separation and HPLC of both intact and desulfated GSLs, detection and structure determination by UV, MS, NMR and susceptibility to myrosinase) with emphasis on the use of NMR for structure elucidation of even minor GSLs and GSL hydrolysis products. The example includes identification of a novel GSL, (R)-2-hydroxy-2-(3-hydroxyphenyl)ethylglucosinolate.Recent investigations of GSL evolution, based on investigations of species with well established phylogeny, are reviewed. From the relatively few such investigations, it is already clear that GSL profiles are regularly subject to evolution. This result is compatible with natural selection for specific GSL side chains. The probable existence of structure-specific GSL catabolism in intact plants suggests that biochemical evolution of GSLs has more complex implications than the mere liberation of a different hydrolysis product upon tissue disruption.