Adaptive Evolution: Can we read the future from a tree?

  1. Michael Lässig  Is a corresponding author
  2. Marta Łuksza
  1. University of Cologne, Germany
  2. Institute for Advanced Study, United States

Darwinian evolution is a dynamical principle that connects the past and the future. According to this principle, fitness differences between the individuals in a population are an important driving force of evolution. Biologists have long used fitness effects to explain observed evolutionary changes. For good reasons, however, they have been hesitant to make predictions about the future of a species. Given the bewildering complexity of what is possible in evolution, attempting to say what will happen in a specific instance may appear futile. Moreover, we cannot test any predictions, because we have not seen the evolutionary past and will not see the future.

Recently, however, evolutionary biology is gaining predictive power in an increasing number of systems, which include viruses, bacteria and populations of cancer cells. In these systems, high mutation rates make evolution happen in front of our eyes. Every year, for example, the human influenza virus replaces 2% of the amino acids in the protein domains that interact with the immune system of its host. Using modern genome sequencing, we can now monitor the genetic history of entire populations and reconstruct their genealogical trees. Such trees show how the individuals of today's populations are connected to their evolutionary ancestors. Now, in eLife, Richard Neher, Colin Russell and Boris Shraiman investigate how much these trees can tell us about the future of a population (Neher et al., 2014).

Inferring evolutionary patterns from genealogical trees has a long history. Geneticists use probabilistic methods to map mutations onto specific tree branches (Figure 1A). Counting how often these mutations appear in different lineages tells us which fitness effects are predominant in a population (McDonald and Kreitman, 1991; Strelkowa and Lässig, 2012). From the statistics of the genealogical tree itself, epidemiologists infer the growth rate of pathogen populations and use that information to predict the future course of an epidemic (Figure 1B, Stadler, 2010). Neher, Russell and Shraiman—who are at the Max Planck Institute for Developmental Biology, the University of Cambridge, and the University of California at Santa Barbara, respectively—extend this genealogy-based inference to genetic changes within a population (Figure 1C). This required developing new ways to extract information from genealogical trees: predictions must now be made for clades of genetically similar individuals, so we need a model that captures growth rate differences between different clades within one genealogical tree.

Fitness inference from genealogical trees.

Lineages in these trees connect the individuals in a population sample and their evolutionary ancestors, which are the nodes of the tree. Evolutionarily successful lineages have descendants in the far future and are marked by thick lines; all other lineages are lost in the evolutionary process. (A) The relative numbers of mutations in successful and in lost lineages measure the predominant fitness effects in a population (orange dots: amino acid changes, blue dots: synonymous mutations). (B) The global statistics of nodes and branches measures the absolute rate of exponential population growth (indicated by the shaded area). (C) The local statistics of nodes and branches measures growth rate differences between clades. Neher and colleagues use this information to predict clade evolution.

To meet this challenge, Neher and colleagues build on a formalism that is rooted in statistical physics and has become a major new development in population genetics (Tsimring et al., 1996; Rouzine et al., 2003; Desai and Fisher, 2007). The basic idea is simple. Given that fitness differences within a population are carried by genetic mutations, we can imagine splitting each mutation and its fitness effect into ever-smaller pieces. This leads to a model in which the overall fitness variation of a population is made up of many small-effect mutations. By the law of large numbers, the fitness distribution then becomes bell-shaped. Such distributions are called travelling fitness waves (Tsimring et al., 1996). In a given lineage, the accumulation of many small fitness effects follows a diffusive random walk. This picture applies to fast adaptive processes in asexual populations where the expansion of a successful clade is fuelled by multiple beneficial mutations—for example, when viruses evolve to escape their hosts' immune defences.

Neher and colleagues link their fitness wave model to simpler heuristic measures of growth, which can easily be used to analyse data from a large genealogical tree. Specifically, they look at the local tree ‘volume’ λ(τ), which sums all tree branches in the vicinity of a given node with a discounting scale τ. This quantity provides a (nonlinear) measure of how fast the number of individuals grows around that node. For example, in a subtree growing exponentially with rate r, the volume λ(τ) equals simply τ/(1–τr). By interpreting this growth rate as fitness, Neher and colleagues obtain a measure of fitness differences between clades. A substantial fraction of the local tree volume is generated by small-effect mutations ‘hitch-hiking’ in successful clades (for example synonymous mutations, which do not change a protein). This explains why the local tree volume is closely related to fitness measures used in previous prediction schemes (Łuksza and Lässig, 2014).

The key strength of this method is that it uses only the information contained in a genealogical tree. Thus, it can be applied in cases where we do not know which functions undergo adaptive evolution or where in the genome they are encoded. This feature is also important for interpreting the results: genealogy-based inference reveals growth rate differences within a population sample, but it remains agnostic about their cause. In the fitness wave model, adaptive evolution is that cause, but the demographic structure of a population or variations in sampling density may generate a similar signal in tree data.

Neher and colleagues apply their method to predict the evolution of the human influenza virus A/H3N2. This is a challenging problem: one year in advance, we need to forecast the prevalent clades circulating in a given winter season. Despite the simplicity of their method, Neher and colleagues predict the ancestor sequence of next year's clades with remarkable accuracy for the majority of northern winters between 1995 and 2013.

We do not yet know in detail how the genetic evolution of the influenza virus is related to its interactions with the human immune system. These ‘antigenic’ properties determine how effective influenza vaccines are. They depend on a smaller number of mutations, some of which have individually large effects (Koel et al., 2013). Thus, prediction schemes geared towards antigenic properties must go beyond examining the overall sequence genealogies and weigh mutations by their antigenic effect (Bedford et al., 2014; Łuksza and Lässig, 2014).

Altogether, as Neher and colleagues show, current predictions reach about halfway between random picks and optimal predictions. This poses big conceptual and practical questions: How much can future methods improve on that score? And where does the inherent unpredictability of evolution start? Prediction is the ultimate test of any dynamical principle. Quantitative evolutionary science is being put to that test now.

References

    1. Rouzine IM
    2. Wakeley J
    3. Coffin JM
    (2003) The solitary wave of asexual evolution
    Proceedings of the National Academy of Sciences of USA 100:587–592.
    https://doi.org/10.1073/pnas.242719299

Article and author information

Author details

  1. Michael Lässig

    Institute for Theoretical Physics, University of Cologne, Cologne, Germany
    For correspondence
    mlaessig@uni-koeln.de
    Competing interests
    The authors declare that no competing interests exist.
  2. Marta Łuksza

    Institute for Advanced Study, Princeton, United States
    Competing interests
    The authors declare that no competing interests exist.

Publication history

  1. Version of Record published:

Copyright

© 2014, Lässig and Łuksza

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 1,874
    views
  • 151
    downloads
  • 2
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Michael Lässig
  2. Marta Łuksza
(2014)
Adaptive Evolution: Can we read the future from a tree?
eLife 3:e05060.
https://doi.org/10.7554/eLife.05060
  1. Further reading

Further reading

    1. Evolutionary Biology
    Ljiljana Mihajlovic, Bharat Ravi Iyengar ... Yolanda Schaerli
    Research Article

    Gene duplication drives evolution by providing raw material for proteins with novel functions. An influential hypothesis by Ohno (1970) posits that gene duplication helps genes tolerate new mutations and thus facilitates the evolution of new phenotypes. Competing hypotheses argue that deleterious mutations will usually inactivate gene duplicates too rapidly for Ohno’s hypothesis to work. We experimentally tested Ohno’s hypothesis by evolving one or exactly two copies of a gene encoding a fluorescent protein in Escherichia coli through several rounds of mutation and selection. We analyzed the genotypic and phenotypic evolutionary dynamics of the evolving populations through high-throughput DNA sequencing, biochemical assays, and engineering of selected variants. In support of Ohno’s hypothesis, populations carrying two gene copies displayed higher mutational robustness than those carrying a single gene copy. Consequently, the double-copy populations experienced relaxed purifying selection, evolved higher phenotypic and genetic diversity, carried more mutations and accumulated combinations of key beneficial mutations earlier. However, their phenotypic evolution was not accelerated, possibly because one gene copy rapidly became inactivated by deleterious mutations. Our work provides an experimental platform to test models of evolution by gene duplication, and it supports alternatives to Ohno’s hypothesis that point to the importance of gene dosage.

    1. Cell Biology
    2. Evolutionary Biology
    Paul Richard J Yulo, Nicolas Desprat ... Heather L Hendrickson
    Research Article

    Maintenance of rod-shape in bacterial cells depends on the actin-like protein MreB. Deletion of mreB from Pseudomonas fluorescens SBW25 results in viable spherical cells of variable volume and reduced fitness. Using a combination of time-resolved microscopy and biochemical assay of peptidoglycan synthesis, we show that reduced fitness is a consequence of perturbed cell size homeostasis that arises primarily from differential growth of daughter cells. A 1000-generation selection experiment resulted in rapid restoration of fitness with derived cells retaining spherical shape. Mutations in the peptidoglycan synthesis protein Pbp1A were identified as the main route for evolutionary rescue with genetic reconstructions demonstrating causality. Compensatory pbp1A mutations that targeted transpeptidase activity enhanced homogeneity of cell wall synthesis on lateral surfaces and restored cell size homeostasis. Mechanistic explanations require enhanced understanding of why deletion of mreB causes heterogeneity in cell wall synthesis. We conclude by presenting two testable hypotheses, one of which posits that heterogeneity stems from non-functional cell wall synthesis machinery, while the second posits that the machinery is functional, albeit stalled. Overall, our data provide support for the second hypothesis and draw attention to the importance of balance between transpeptidase and glycosyltransferase functions of peptidoglycan building enzymes for cell shape determination.