A simple rule for predicting function of microbial communities.

Microbial communities perform many important functions, such as carbon sequestration, decomposition, pathogen resistance, etc., but quantitatively predicting functions of new communities remains a major challenge. In this issue of Cell, Diaz-Colunga et al. report a new simple statistical regularity that enables such predictions.

Best Practices in Designing, Sequencing, and Identifying Random DNA Barcodes.

Random DNA barcodes are a versatile tool for tracking cell lineages, with applications ranging from development to cancer to evolution. Here, we review and critically evaluate barcode designs as well as methods of barcode sequencing and initial processing of barcode data. We first demonstrate how various barcode design decisions affect data quality and propose a new design that balances all considerations that we are currently aware of. We then discuss various options for the preparation of barcode sequencing libraries, including inline indices and Unique Molecular Identifiers (UMIs). Finally, we test the performance of several established and new bioinformatic pipelines for the extraction of barcodes from raw sequencing reads and for error correction. We find that both alignment and regular expression-based approaches work well for barcode extraction, and that error-correction pipelines designed specifically for barcode data are superior to generic ones. Overall, this review will help researchers to approach their barcoding experiments in a deliberate and systematic way.

Evolutionary stalling and a limit on the power of natural selection to improve a cellular module.

Cells consist of molecular modules which perform vital biological functions. Cellular modules are key units of adaptive evolution because organismal fitness depends on their performance. Theory shows that in rapidly evolving populations, such as those of many microbes, adaptation is driven primarily by common beneficial mutations with large effects, while other mutations behave as if they are effectively neutral. As a consequence, if a module can be improved only by rare and/or weak beneficial mutations, its adaptive evolution would stall. However, such evolutionary stalling has not been empirically demonstrated, and it is unclear to what extent stalling may limit the power of natural selection to improve modules. Here we empirically characterize how natural selection improves the translation machinery (TM), an essential cellular module. We experimentally evolved populations of Escherichia coli with genetically perturbed TMs for 1,000 generations. Populations with severe TM defects initially adapted via mutations in the TM, but TM adaptation stalled within about 300 generations. We estimate that the genetic load in our populations incurred by residual TM defects ranges from 0.5 to 19%. Finally, we found evidence that both epistasis and the depletion of the pool of beneficial mutations contributed to evolutionary stalling. Our results suggest that cellular modules may not be fully optimized by natural selection despite the availability of adaptive mutations.

Modular epistasis and the compensatory evolution of gene deletion mutants.

Screens for epistatic interactions have long been used to characterize functional relationships corresponding to protein complexes, metabolic pathways, and other functional modules. Although epistasis between adaptive mutations is also common in laboratory evolution experiments, the functional basis for these interactions is less well characterized. Here, we quantify the extent to which gene function (as determined by a genome-wide screen for epistasis among deletion mutants) influences the rate and genetic basis of compensatory adaptation in a set of 37 gene deletion mutants nested within 16 functional modules. We find that functional module has predictive power: mutants with deletions in the same module tend to adapt more similarly, on average, than those with deletions in different modules. At the same time, initial fitness also plays a role: independent of the specific functional modules involved, adaptive mutations tend to be systematically more beneficial in less-fit genetic backgrounds, consistent with a general pattern of diminishing returns epistasis. We measured epistatic interactions between initial gene deletion mutations and the mutations that accumulate during compensatory adaptation and found a general trend towards positive epistasis (i.e. mutations tend to be more beneficial in the background in which they arose). In two functional modules, epistatic interactions between the initial gene deletions and the mutations in their descendant lines caused evolutionary entrenchment, indicating an intimate functional relationship. Our results suggest that genotypes with similar epistatic interactions with gene deletion mutations will also have similar epistatic interactions with adaptive mutations, meaning that genome scale maps of epistasis between gene deletion mutations can be predictive of evolutionary dynamics.

Coordinated Evolution of Influenza A Surface Proteins.

The surface proteins hemagglutinin (HA) and neuraminidase (NA) of human influenza A virus evolve under selection pressures to escape adaptive immune responses and antiviral drug treatments. In addition to these external selection pressures, some mutations in HA are known to affect the adaptive landscape of NA, and vice versa, because these two proteins are physiologically interlinked. However, the extent to which evolution of one protein affects the evolution of the other one is unknown. Here we develop a novel phylogenetic method for detecting the signatures of such genetic interactions between mutations in different genes - that is, inter-gene epistasis. Using this method, we show that influenza surface proteins evolve in a coordinated way, with mutations in HA affecting subsequent spread of mutations in NA and vice versa, at many sites. Of particular interest is our finding that the oseltamivir-resistance mutations in NA in subtype H1N1 were likely facilitated by prior mutations in HA. Our results illustrate that the adaptive landscape of a viral protein is remarkably sensitive to its genomic context and, more generally, that the evolution of any single protein must be understood within the context of the entire evolving genome.

Loss and recovery of genetic diversity in adapting populations of HIV.

The evolution of drug resistance in HIV occurs by the fixation of specific, well-known, drug-resistance mutations, but the underlying population genetic processes are not well understood. By analyzing within-patient longitudinal sequence data, we make four observations that shed a light on the underlying processes and allow us to infer the short-term effective population size of the viral population in a patient. Our first observation is that the evolution of drug resistance usually occurs by the fixation of one drug-resistance mutation at a time, as opposed to several changes simultaneously. Second, we find that these fixation events are accompanied by a reduction in genetic diversity in the region surrounding the fixed drug-resistance mutation, due to the hitchhiking effect. Third, we observe that the fixation of drug-resistance mutations involves both hard and soft selective sweeps. In a hard sweep, a resistance mutation arises in a single viral particle and drives all linked mutations with it when it spreads in the viral population, which dramatically reduces genetic diversity. On the other hand, in a soft sweep, a resistance mutation occurs multiple times on different genetic backgrounds, and the reduction of diversity is weak. Using the frequency of occurrence of hard and soft sweeps we estimate the effective population size of HIV to be 1.5 x 10(5) (95% confidence interval [0.8 x 10(5),4.8 x 10(5)]). This number is much lower than the actual number of infected cells, but much larger than previous population size estimates based on synonymous diversity. We propose several explanations for the observed discrepancies. Finally, our fourth observation is that genetic diversity at non-synonymous sites recovers to its pre-fixation value within 18 months, whereas diversity at synonymous sites remains depressed after this time period. These results improve our understanding of HIV evolution and have potential implications for treatment strategies.

Prevalence of epistasis in the evolution of influenza A surface proteins.

The surface proteins of human influenza A viruses experience positive selection to escape both human immunity and, more recently, antiviral drug treatments. In bacteria and viruses, immune-escape and drug-resistant phenotypes often appear through a combination of several mutations that have epistatic effects on pathogen fitness. However, the extent and structure of epistasis in influenza viral proteins have not been systematically investigated. Here, we develop a novel statistical method to detect positive epistasis between pairs of sites in a protein, based on the observed temporal patterns of sequence evolution. The method rests on the simple idea that a substitution at one site should rapidly follow a substitution at another site if the sites are positively epistatic. We apply this method to the surface proteins hemagglutinin and neuraminidase of influenza A virus subtypes H3N2 and H1N1. Compared to a non-epistatic null distribution, we detect substantial amounts of epistasis and determine the identities of putatively epistatic pairs of sites. In particular, using sequence data alone, our method identifies epistatic interactions between specific sites in neuraminidase that have recently been demonstrated, in vitro, to confer resistance to the drug oseltamivir; these epistatic interactions are responsible for widespread drug resistance among H1N1 viruses circulating today. This experimental validation demonstrates the predictive power of our method to identify epistatic sites of importance for viral adaptation and public health. We conclude that epistasis plays a large role in shaping the molecular evolution of influenza viruses. In particular, sites with , which would normally not be identified as positively selected, can facilitate viral adaptation through epistatic interactions with their partner sites. The knowledge of specific interactions among sites in influenza proteins may help us to predict the course of antigenic evolution and, consequently, to select more appropriate vaccines and drugs.

Environment-independent distribution of mutational effects emerges from microscopic epistasis.

Predicting how new mutations alter phenotypes is difficult because mutational effects vary across genotypes and environments. Recently discovered global epistasis, in which the fitness effects of mutations scale with the fitness of the background genotype, can improve predictions, but how the environment modulates this scaling is unknown. We measured the fitness effects of ~100 insertion mutations in 42 strains of Saccharomyces cerevisiae in six laboratory environments and found that the global epistasis scaling is nearly invariant across environments. Instead, the environment tunes one global parameter, the background fitness at which most mutations switch sign. As a consequence, the distribution of mutational effects is predictable across genotypes and environments. Our results suggest that the effective dimensionality of genotype-to-phenotype maps across environments is surprisingly low.

Young proteins experience more variable selection pressures than old proteins.

It is well known that young proteins tend to experience weaker purifying selection and evolve more quickly than old proteins. Here, we show that, in addition, young proteins tend to experience more variable selection pressures over time than old proteins. We demonstrate this pattern in three independent taxonomic groups: yeast, Drosophila, and mammals. The increased variability of selection pressures on young proteins is highly significant even after controlling for the fact that young proteins are typically shorter and experience weaker purifying selection than old proteins. The majority of our results are consistent with the hypothesis that the function of a young gene tends to change over time more readily than that of an old gene. At the same time, our results may be caused in part by young genes that serve constant functions over time, but nevertheless appear to evolve under changing selection pressures due to depletion of adaptive mutations. In either case, our results imply that the evolution of a protein-coding sequence is partly determined by its age and origin, and not only by the phenotypic properties of the encoded protein. We discuss, via specific examples, the consequences of these findings for understanding of the sources of evolutionary novelty.

Evolutionary repeatability of emergent properties of ecological communities.

Most species belong to ecological communities where their interactions give rise to emergent community-level properties, such as diversity and productivity. Understanding and predicting how these properties change over time has been a major goal in ecology, with important practical implications for sustainability and human health. Less attention has been paid to the fact that community-level properties can also change because member species evolve. Yet, our ability to predict long-term eco-evolutionary dynamics hinges on how repeatably community-level properties change as a result of species evolution. Here, we review studies of evolution of both natural and experimental communities and make the case that community-level properties at least sometimes evolve repeatably. We discuss challenges faced in investigations of evolutionary repeatability. In particular, only a handful of studies enable us to quantify repeatability. We argue that quantifying repeatability at the community level is critical for approaching what we see as three major open questions in the field: (i) Is the observed degree of repeatability surprising? (ii) How is evolutionary repeatability at the community level related to repeatability at the level of traits of member species? (iii) What factors affect repeatability? We outline some theoretical and empirical approaches to addressing these questions. Advances in these directions will not only enrich our basic understanding of evolution and ecology but will also help us predict eco-evolutionary dynamics. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Environment-independent distribution of mutational effects emerges from microscopic epistasis.

Predicting how new mutations alter phenotypes is difficult because mutational effects vary across genotypes and environments. Recently discovered global epistasis, where the fitness effects of mutations scale with the fitness of the background genotype, can improve predictions, but how the environment modulates this scaling is unknown. We measured the fitness effects of ~100 insertion mutations in 42 strains of Saccharomyces cerevisiae in six laboratory environments and found that the global-epistasis scaling is nearly invariant across environments. Instead, the environment tunes one global parameter, the background fitness at which most mutations switch sign. As a consequence, the distribution of mutational effects is remarkably predictable across genotypes and environments. Our results suggest that the effective dimensionality of genotype-to-phenotype maps across environments is surprisingly low.

Mutualism-enhancing mutations dominate early adaptation in a two-species microbial community.

Species interactions drive evolution while evolution shapes these interactions. The resulting eco-evolutionary dynamics and their repeatability depend on how adaptive mutations available to community members affect fitness and ecologically relevant traits. However, the diversity of adaptive mutations is not well characterized, and we do not know how this diversity is affected by the ecological milieu. Here we use barcode lineage tracking to address this question in a community of yeast Saccharomyces cerevisiae and alga Chlamydomonas reinhardtii that have a net commensal relationship that results from a balance between competitive and mutualistic interactions. We find that yeast has access to many adaptive mutations with diverse ecological consequences, in particular those that increase and reduce the yields of both species. The presence of the alga does not change which mutations are adaptive in yeast (that is, there is no fitness trade-off for yeast between growing alone or with alga), but rather shifts selection to favour yeast mutants that increase the yields of both species and make the mutualism stronger. Thus, in the presence of the alga, adaptative mutations contending for fixation in yeast are more likely to enhance the mutualism, even though cooperativity is not directly favoured by natural selection in our system. Our results demonstrate that ecological interactions not only alter the trajectory of evolution but also dictate its repeatability; in particular, weak mutualisms can repeatably evolve to become stronger.

Upper Bound on the Mutational Burden Imposed by a CRISPR-Cas9 Gene-Drive Element.

CRISPR-Cas9 gene drives (CCGDs) are powerful tools for genetic control of wild populations, useful for eradication of disease vectors, conservation of endangered species and other applications. However, Cas9 alone and in a complex with gRNA can cause double-stranded DNA breaks at off-target sites, which could increase the mutational load and lead to loss of heterozygosity (LOH). These undesired effects raise potential concerns about the long-term evolutionary safety of CCGDs, but the magnitude of these effects is unknown. To estimate how the presence of a CCGD or a Cas9 alone in the genome affects the rates of LOH events and de novo mutations, we carried out a mutation accumulation experiment in yeast Saccharomyces cerevisiae. Despite its substantial statistical power, our experiment revealed no detectable effect of CCGD or Cas9 alone on the genome-wide rates of mutations or LOH events, suggesting that these rates are affected by less than 30%. Nevertheless, we found that Cas9 caused a slight but significant shift towards more interstitial and fewer terminal LOH events, and the CCGD caused a significant difference in the distribution of LOH events on Chromosome V. Taken together, our results show that these genetic elements impose a weak and likely localized additional mutational burden in the yeast model. Although the mutagenic effects of CCGDs need to be further evaluated in other systems, our results suggest that the effect of CCGDs on off-target mutation rates and genetic diversity may be acceptable.

Multidrug therapy and evolution of antibiotic resistance: when order matters.

The evolution of drug resistance among pathogenic bacteria has led public health workers to rely increasingly on multidrug therapy to treat infections. Here, we compare the efficacy of combination therapy (i.e., using two antibiotics simultaneously) and sequential therapy (i.e., switching two antibiotics) in minimizing the evolution of multidrug resistance. Using in vitro experiments, we show that the sequential use of two antibiotics against Pseudomonas aeruginosa can slow down the evolution of multiple-drug resistance when the two antibiotics are used in a specific order. A simple population dynamics model reveals that using an antibiotic associated with high costs of resistance first minimizes the chance of multidrug resistance evolution during sequential therapy under limited mutation supply rate. As well as presenting a novel approach to multidrug therapy, this work shows that costs of resistance not only influences the persistence of antibiotic-resistant bacteria but also plays an important role in the emergence of resistance.

The dynamics of adaptation on correlated fitness landscapes.

Evolutionary theory predicts that a population in a new environment will accumulate adaptive substitutions, but precisely how they accumulate is poorly understood. The dynamics of adaptation depend on the underlying fitness landscape. Virtually nothing is known about fitness landscapes in nature, and few methods allow us to infer the landscape from empirical data. With a view toward this inference problem, we have developed a theory that, in the weak-mutation limit, predicts how a population's mean fitness and the number of accumulated substitutions are expected to increase over time, depending on the underlying fitness landscape. We find that fitness and substitution trajectories depend not on the full distribution of fitness effects of available mutations but rather on the expected fixation probability and the expected fitness increment of mutations. We introduce a scheme that classifies landscapes in terms of the qualitative evolutionary dynamics they produce. We show that linear substitution trajectories, long considered the hallmark of neutral evolution, can arise even when mutations are strongly selected. Our results provide a basis for understanding the dynamics of adaptation and for inferring properties of an organism's fitness landscape from temporal data. Applying these methods to data from a long-term experiment, we infer the sign and strength of epistasis among beneficial mutations in the Escherichia coli genome.

Isolation and characterization of adenoviruses persistently shed from the gastrointestinal tract of non-human primates.

Adenoviruses are important human pathogens that have been developed as vectors for gene therapies and genetic vaccines. Previous studies indicated that human infections with adenoviruses are self-limiting in immunocompetent hosts with evidence of some persistence in adenoid tissue. We sought to better understand the natural history of adenovirus infections in various non-human primates and discovered that healthy populations of great apes (chimpanzees, bonobos, gorillas, and orangutans) and macaques shed substantial quantities of infectious adenoviruses in stool. Shedding in stools from asymptomatic humans was found to be much less frequent, comparable to frequencies reported before. We purified and fully sequenced 30 novel adenoviruses from apes and 3 novel adenoviruses from macaques. Analyses of the new ape adenovirus sequences (as well as the 4 chimpanzee adenovirus sequences we have previously reported) together with 22 complete adenovirus genomes available from GenBank revealed that (a) the ape adenoviruses could clearly be classified into species corresponding to human adenovirus species B, C, and E, (b) there was evidence for intraspecies recombination between adenoviruses, and (c) the high degree of phylogenetic relatedness of adenoviruses across their various primate hosts provided evidence for cross species transmission events to have occurred in the natural history of B and E viruses. The high degree of asymptomatic shedding of live adenovirus in non-human primates and evidence for zoonotic transmissions warrants caution for primate handling and housing. Furthermore, the presence of persistent and/or latent adenovirus infections in the gut should be considered in the design and interpretation of human and non-human primate studies with adenovirus vectors.

The dynamics of social networks among female Asian elephants.

BACKGROUND: Patterns in the association of individuals can shed light on the underlying conditions and processes that shape societies. Here we characterize patterns of association in a population of wild Asian Elephants at Uda Walawe National Park in Sri Lanka. We observed 286 individually-identified adult female elephants over 20 months and examined their social dynamics at three levels of organization: pairs of individuals (dyads), small sets of direct companions (ego-networks), and the population level (complete networks). RESULTS: Corroborating previous studies of this and other Asian elephant populations, we find that the sizes of elephant groups observed in the field on any particular day are typically small and that rates of association are low. In contrast to earlier studies, our longitudinal observations reveal that individuals form larger social units that can be remarkably stable across years while associations among such units change across seasons. Association rates tend to peak in dry seasons as opposed to wet seasons, with some cyclicity at the level of dyads. In addition, we find that individuals vary substantially in their fidelity to companions. At the ego-network level, we find that despite these fluctuations, individuals associate with a pool of long-term companions. At the population level, social networks do not exhibit any clear seasonal structure or hierarchical stratification. CONCLUSIONS: This detailed longitudinal study reveals different social dynamics at different levels of organization. Taken together, these results demonstrate that low association rates, seemingly small group sizes, and fission-fusion grouping behavior mask hidden stability in the extensive and fluid social affiliations in this population of Asian elephants.

The population genetics of dN/dS.

Evolutionary pressures on proteins are often quantified by the ratio of substitution rates at non-synonymous and synonymous sites. The dN/dS ratio was originally developed for application to distantly diverged sequences, the differences among which represent substitutions that have fixed along independent lineages. Nevertheless, the dN/dS measure is often applied to sequences sampled from a single population, the differences among which represent segregating polymorphisms. Here, we study the expected dN/dS ratio for samples drawn from a single population under selection, and we find that in this context, dN/dS is relatively insensitive to the selection coefficient. Moreover, the hallmark signature of positive selection over divergent lineages, dN/dS>1, is violated within a population. For population samples, the relationship between selection and dN/dS does not follow a monotonic function, and so it may be impossible to infer selection pressures from dN/dS. These results have significant implications for the interpretation of dN/dS measurements among population-genetic samples.

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.