Factors affecting the genetic load in Drosophila: synergistic epistasis and correlations among fitness components

The origination events of gametic sexual reproduction and anisogamy

The evolution of gametic sex (meiosis and fertilization) and subsequent transition from isogamy (fusion between two equal-sized gametes) to anisogamy (dimorphism into eggs and sperm, namely, females and males) is one of the largest enigmas of evolutionary biology. Meiosis entails genome-dilution cost and anisogamy entails male-production cost. Despite much progress has been made for the maintenance mechanisms of sex, its origination events under such “twofold cost of sex” are still unsolved. Here, we posit two hypothetical scenarios as follows: the “Seesaw Effect” hypothesizes that automictic selfing between isogametes effectively purged deleterious mutations from an organism’s lineage and simultaneously fixed the sex-controlling allele and all other loci (no genome-dilution cost raised). The high relatedness among homoeologous cell colonies led to multicellularization. The “inflated isogamy” hypothesizes that multicellularity increased the reproductive investment of both mates, resulting in excessively large isogametes. This redundancy induced cheating of one sex (evolving to male) to reduce gamete size. However, the other sex (evolving to female) allowed this cheat because her cost did not change. Therefore, anisogamy originated as a kind of commensalism but turned into beneficial for females because it solved the gamete limitation problem inherent to isogamy. Thus, smooth transition to anisogamy had been attained.

A Family of Fitness Landscapes Modeled through Gene Regulatory Networks

Fitness landscapes are a powerful metaphor for understanding the evolution of biological systems. These landscapes describe how genotypes are connected to each other through mutation and related through fitness. Empirical studies of fitness landscapes have increasingly revealed conserved topographical features across diverse taxa, e.g., the accessibility of genotypes and "ruggedness". As a result, theoretical studies are needed to investigate how evolution proceeds on fitness landscapes with such conserved features. Here, we develop and study a model of evolution on fitness landscapes using the lens of Gene Regulatory Networks (GRNs), where the regulatory products are computed from multiple genes and collectively treated as phenotypes. With the assumption that regulation is a binary process, we prove the existence of empirically observed, topographical features such as accessibility and connectivity. We further show that these results hold across arbitrary fitness functions and that a trade-off between accessibility and ruggedness need not exist. Then, using graph theory and a coarse-graining approach, we deduce a mesoscopic structure underlying GRN fitness landscapes where the information necessary to predict a population's evolutionary trajectory is retained with minimal complexity. Using this coarse-graining, we develop a bottom-up algorithm to construct such mesoscopic backbones, which does not require computing the genotype network and is therefore far more efficient than brute-force approaches. Altogether, this work provides mathematical results of high-dimensional fitness landscapes and a path toward connecting theory to empirical studies.

The target of selection matters: An established resistance-development-time negative genetic trade-off is not found when selecting on development time

Trade-offs are fundamental to evolutionary outcomes and play a central role in eco-evolutionary theory. They are often examined by experimentally selecting on one life-history trait and looking for negative correlations in other traits. For example, populations of the moth Plodia interpunctella selected to resist viral infection show a life-history cost with longer development times. However, we rarely examine whether the detection of such negative genetic correlations depends on the trait on which we select. Here, we examine a well-characterized negative genotypic trade-off between development time and resistance to viral infection in the moth Plodia interpunctella and test whether selection on a phenotype known to be a cost of resistance (longer development time) leads to the predicted correlated increase in resistance. If there is tight pleiotropic relationship between genes that determine development time and resistance underpinning this trade-off, we might expect increased resistance when we select on longer development time. However, we show that selecting for longer development time in this system selects for reduced resistance when compared to selection for shorter development time. This shows how phenotypes typically characterized by a trade-off can deviate from that trade-off relationship, and suggests little genetic linkage between the genes governing viral resistance and those that determine response to selection on the key life-history trait. Our results are important for both selection strategies in applied biological systems and for evolutionary modelling of host-parasite interactions.

Does operational sex ratio influence relative strength of purging selection in males versus females?

According to theory, sexual selection in males may efficiently purge mutation load of sexual populations, reducing or fully compensating 'the cost of males'. For this to occur, mutations not only need to be deleterious to both sexes, they also must affect males more than females. A frequently overlooked problem is that relative strength of selection on males versus females may vary between environments, with social conditions being particularly likely to affect selection in males and females differently. Here, we induced mutations in red flour beetles (Tribolium castaneum) and tested their effect in both sexes under three different operational sex ratios (1:2, 1:1 and 2:1). Induced mutations decreased fitness of both males and females, but their effect was not stronger in males. Surprisingly, operational sex ratio did not affect selection against deleterious mutations nor its relative strength in the sexes. Thus, our results show no support for the role of sexual selection in the evolutionary maintenance of sex.

Accelerated inbreeding depression suggests synergistic epistasis for deleterious mutations in Drosophila melanogaster

Epistasis may have important consequences for a number of issues in quantitative genetics and evolutionary biology. In particular, synergistic epistasis for deleterious alleles is relevant to the mutation load paradox and the evolution of sex and recombination. Some studies have shown evidence of synergistic epistasis for spontaneous or induced deleterious mutations appearing in mutation-accumulation experiments. However, many newly arising mutations may not actually be segregating in natural populations because of the erasing action of natural selection. A demonstration of synergistic epistasis for naturally segregating alleles can be achieved by means of inbreeding depression studies, as deleterious recessive allelic effects are exposed in inbred lines. Nevertheless, evidence of epistasis from these studies is scarce and controversial. In this paper, we report the results of two independent inbreeding experiments carried out with two different populations of Drosophila melanogaster. The results show a consistent accelerated inbreeding depression for fitness, suggesting synergistic epistasis among deleterious alleles. We also performed computer simulations assuming different possible models of epistasis and mutational parameters for fitness, finding some of them to be compatible with the results observed. Our results suggest that synergistic epistasis for deleterious mutations not only occurs among newly arisen spontaneous or induced mutations, but also among segregating alleles in natural populations.

DISSEQT-DIStribution-based modeling of SEQuence space Time dynamics

Rapidly evolving microbes are a challenge to model because of the volatile, complex, and dynamic nature of their populations. We developed the DISSEQT pipeline (DIStribution-based SEQuence space Time dynamics) for analyzing, visualizing, and predicting the evolution of heterogeneous biological populations in multidimensional genetic space, suited for population-based modeling of deep sequencing and high-throughput data. The pipeline is openly available on GitHub (https://github.com/rasmushenningsson/DISSEQT.jl, accessed 23 June 2019) and Synapse (https://www.synapse.org/#!Synapse: syn11425758, accessed 23 June 2019), covering the entire workflow from read alignment to visualization of results. Our pipeline is centered around robust dimension and model reduction algorithms for analysis of genotypic data with additional capabilities for including phenotypic features to explore dynamic genotype-phenotype maps. We illustrate its utility and capacity with examples from evolving RNA virus populations, which present one of the highest degrees of genetic heterogeneity within a given population found in nature. Using our pipeline, we empirically reconstruct the evolutionary trajectories of evolving populations in sequence space and genotype-phenotype fitness landscapes. We show that while sequence space is vastly multidimensional, the relevant genetic space of evolving microbial populations is of intrinsically low dimension. In addition, evolutionary trajectories of these populations can be faithfully monitored to identify the key minority genotypes contributing most to evolution. Finally, we show that empirical fitness landscapes, when reconstructed to include minority variants, can predict phenotype from genotype with high accuracy.

Do slower movers have lower reproductive success and higher mutation load?

Deleterious mutations occur frequently in eukaryotes, resulting in individuals carrying multiple alleles that decrease their fitness. At a population level, if unchecked, accumulation of this mutation load can ultimately lead to extinction. How selection counters the accumulation of mutation load, limiting declines in population fitness, is not well understood. Here, we use manipulative experiments in zebrafish (Danio rerio) to investigate the opportunities for selection on mutation load. Inducing high mutation load through mutagenesis, we applied one generation of within‐family selection on locomotor performance and characterized both the direct response to this selection and the indirect response of reproductive success. Offspring of slow swimming parents exhibited age‐dependent declines in swimming speed, whereas their cousins, with faster swimming parents, did not. This pattern mimics previously documented differences between high and low mutation load populations of zebrafish, suggesting that slow swimming siblings inherited (and transmitted) more mutations than their faster swimming siblings. Crosses among offspring of slow swimming fish had, on average, <75% of the reproductive success of crosses among offspring of fast swimming parents, or crosses of offspring of slow swimmers with offspring of fast swimmers. This evidence of mutationally correlated swimming speed and reproductive success reveals the potential for concordant selection on mutation load through different fitness components. There was no evidence that crosses within families (where parents potentially shared the same mutations inherited from their common ancestor) had lower reproductive success than crosses among families, suggesting that viability selection was not acting predominantly through lethal recessive homozygotes. Rather, patterns of reproductive success are suggestive of effects of mutation number per se on embryo viability. Overall, our results highlight the potential for early life mortality to remove deleterious mutations, and the need to account for this mortality when investigating the evolutionary dynamics of mutation load.

Genomic Tools and Selective Breeding in Molluscs

The production of most farmed molluscs, including mussels, oysters, scallops, abalone, and clams, is heavily dependent on natural seed from the plankton. Closing the lifecycle of species in hatcheries can secure independence from wild stocks and enables long-term genetic improvement of broodstock through selective breeding. Genomic techniques have the potential to revolutionize hatchery-based selective breeding by improving our understanding of the characteristics of mollusc genetics that can pose a challenge for intensive aquaculture and by providing a new suite of tools for genetic improvement. Here we review characteristics of the life history and genetics of molluscs including high fecundity, self-fertilization, high genetic diversity, genetic load, high incidence of deleterious mutations and segregation distortion, and critically assess their impact on the design and effectiveness of selective breeding strategies. A survey of the results of current breeding programs in the literature show that selective breeding with inbreeding control is likely the best strategy for genetic improvement of most molluscs, and on average growth rate can be improved by 10% per generation and disease resistance by 15% per generation across the major farmed species by implementing individual or family-based selection. Rapid advances in sequencing technology have resulted in a wealth of genomic resources for key species with the potential to greatly improve hatchery-based selective breeding of molluscs. In this review, we catalog the range of genomic resources currently available for molluscs of aquaculture interest and discuss the bottlenecks, including lack of high-quality reference genomes and the relatively high cost of genotyping, as well as opportunities for applying genomics-based selection.

Evolutionary constraints in fitness landscapes

In the last years, several genotypic fitness landscapes-combinations of a small number of mutations-have been experimentally resolved. To learn about the general properties of "real" fitness landscapes, it is key to characterize these experimental landscapes via simple measures of their structure, related to evolutionary features. Some of the most relevant measures are based on the selectively acessible paths and their properties. In this paper, we present some measures of evolutionary constraints based on (i) the similarity between accessible paths and (ii) the abundance and characteristics of "chains" of obligatory mutations, that are paths going through genotypes with a single fitter neighbor. These measures have a clear evolutionary interpretation. Furthermore, we show that chains are only weakly correlated to classical measures of epistasis. In fact, some of these measures of constraint are non-monotonic in the amount of epistatic interactions, but have instead a maximum for intermediate values. Finally, we show how these measures shed light on evolutionary constraints and predictability in experimentally resolved landscapes.

The Influence of Higher-Order Epistasis on Biological Fitness Landscape Topography

The effect of a mutation on the organism often depends on what other mutations are already present in its genome. Geneticists refer to such mutational interactions as epistasis. Pairwise epistatic effects have been recognized for over a century, and their evolutionary implications have received theoretical attention for nearly as long. However, pairwise epistatic interactions themselves can vary with genomic background. This is called higher-order epistasis, and its consequences for evolution are much less well understood. Here, we assess the influence that higher-order epistasis has on the topography of 16 published, biological fitness landscapes. We find that on average, their effects on fitness landscape declines with order, and suggest that notable exceptions to this trend may deserve experimental scrutiny. We conclude by highlighting opportunities for further theoretical and experimental work dissecting the influence that epistasis of all orders has on fitness landscape topography and on the efficiency of evolution by natural selection.

Polygenicity and Epistasis Underlie Fitness-Proximal Traits in the Caenorhabditis elegans Multiparental Experimental Evolution (CeMEE) Panel

Understanding the genetic basis of complex traits remains a major challenge in biology. Polygenicity, phenotypic plasticity, and epistasis contribute to phenotypic variance in ways that are rarely clear. This uncertainty can be problematic for estimating heritability, for predicting individual phenotypes from genomic data, and for parameterizing models of phenotypic evolution. Here, we report an advanced recombinant inbred line (RIL) quantitative trait locus mapping panel for the hermaphroditic nematode Caenorhabditis elegans, the C. elegans multiparental experimental evolution (CeMEE) panel. The CeMEE panel, comprising 507 RILs at present, was created by hybridization of 16 wild isolates, experimental evolution for 140-190 generations, and inbreeding by selfing for 13-16 generations. The panel contains 22% of single-nucleotide polymorphisms known to segregate in natural populations, and complements existing C. elegans mapping resources by providing fine resolution and high nucleotide diversity across > 95% of the genome. We apply it to study the genetic basis of two fitness components, fertility and hermaphrodite body size at time of reproduction, with high broad-sense heritability in the CeMEE. While simulations show that we should detect common alleles with additive effects as small as 5%, at gene-level resolution, the genetic architectures of these traits do not feature such alleles. We instead find that a significant fraction of trait variance, approaching 40% for fertility, can be explained by sign epistasis with main effects below the detection limit. In congruence, phenotype prediction from genomic similarity, while generally poor ([Formula: see text]), requires modeling epistasis for optimal accuracy, with most variance attributed to the rapidly evolving chromosome arms.

Sexual selection on spontaneous mutations strengthens the between-sex genetic correlation for fitness

A proposed benefit to sexual selection is that it promotes purging of deleterious mutations from populations. For this benefit to be realized, sexual selection, which is usually stronger on males, must purge mutations deleterious to both sexes. Here, we experimentally test the hypothesis that sexual selection on males purges deleterious mutations that affect both male and female fitness. We measured male and female fitness in two panels of spontaneous mutation-accumulation lines of the fly, Drosophila serrata, each established from a common ancestor. One panel of mutation accumulation lines limited both natural and sexual selection (LS lines), whereas the other panel limited natural selection, but allowed sexual selection to operate (SS lines). Although mutation accumulation caused a significant reduction in male and female fitness in both the LS and SS lines, sexual selection had no detectable effect on the extent of the fitness reduction. Similarly, despite evidence of mutational variance for fitness in males and females of both treatments, sexual selection had no significant impact on the amount of mutational genetic variance for fitness. However, sexual selection did reshape the between-sex correlation for fitness: significantly strengthening it in the SS lines. After 25 generations, the between-sex correlation for fitness was positive but considerably less than one in the LS lines, suggesting that, although most mutations had sexually concordant fitness effects, sex-limited, and/or sex-biased mutations contributed substantially to the mutational variance. In the SS lines this correlation was strong and could not be distinguished from unity. Individual-based simulations that mimick the experimental setup reveal two conditions that may drive our results: (1) a modest-to-large fraction of mutations have sex-limited (or highly sex-biased) fitness effects, and (2) the average fitness effect of sex-limited mutations is larger than the average fitness effect of mutations that affect both sexes similarly.

Two sides of the same coin: A population genetics perspective on lethal mutagenesis and mutational meltdown

The extinction of RNA virus populations upon application of a mutagenic drug is frequently referred to as evidence for the existence of an error threshold, above which the population cannot sustain the mutational load. To explain the extinction process after reaching this threshold, models of lethal mutagenesis have been proposed, in which extinction is described as a deterministic (and thus population size-independent) process. As a separate body of literature, the population genetics community has developed models of mutational meltdown, which focus on the stochastic (and thus population-size dependent) processes governing extinction. However, recent extensions of both models have blurred these boundaries. Here, we first clarify definitions in terms of assumptions, expectations, and relevant parameter spaces, and then assess similarities and differences. As concepts from both fields converge, we argue for a unified theoretical framework that is focused on the evolutionary processes at play, rather than dispute over terminology.

Role of epistasis on the fixation probability of a non-mutator in an adapted asexual population

The mutation rate of a well adapted population is prone to reduction so as to have a lower mutational load. We aim to understand the role of epistatic interactions between the fitness affecting mutations in this process. Using a multitype branching process, the fixation probability of a single non-mutator emerging in a large asexual mutator population is analytically calculated here. The mutator population undergoes deleterious mutations at constant, but at a much higher rate than that of the non-mutator. We find that antagonistic epistasis lowers the chances of mutation rate reduction, while synergistic epistasis enhances it. Below a critical value of epistasis, the fixation probability behaves non-monotonically with variation in the mutation rate of the background population. Moreover, the variation of this critical value of the epistasis parameter with the strength of the mutator is discussed in the appendix. For synergistic epistasis, when selection is varied, the fixation probability reduces overall, with damped oscillations.

Sexual selection expedites the evolution of pesticide resistance

The evolution of insecticide resistance by crop pests and disease vectors causes serious problems for agriculture and health. Sexual selection can accelerate or hinder adaptation to abiotic challenges in a variety of ways, but the effect of sexual selection on resistance evolution is little studied. Here, we examine this question using experimental evolution in the pest insect Tribolium castaneum. The experimental removal of sexual selection slowed the evolution of resistance in populations treated with pyrethroid pesticide, and also reduced the rate at which resistance was lost from pesticide-free populations. These results suggest that selection arising from variance in mating and fertilization success can augment natural selection on pesticide resistance, meaning that sexual selection should be considered when designing strategies to limit the evolution of pesticide resistance.

Epistasis and the Structure of Fitness Landscapes: Are Experimental Fitness Landscapes Compatible with Fisher's Geometric Model?

The fitness landscape defines the relationship between genotypes and fitness in a given environment and underlies fundamental quantities such as the distribution of selection coefficient and the magnitude and type of epistasis. A better understanding of variation in landscape structure across species and environments is thus necessary to understand and predict how populations will adapt. An increasing number of experiments investigate the properties of fitness landscapes by identifying mutations, constructing genotypes with combinations of these mutations, and measuring the fitness of these genotypes. Yet these empirical landscapes represent a very small sample of the vast space of all possible genotypes, and this sample is often biased by the protocol used to identify mutations. Here we develop a rigorous statistical framework based on Approximate Bayesian Computation to address these concerns and use this flexible framework to fit a broad class of phenotypic fitness models (including Fisher's model) to 26 empirical landscapes representing nine diverse biological systems. Despite uncertainty owing to the small size of most published empirical landscapes, the inferred landscapes have similar structure in similar biological systems. Surprisingly, goodness-of-fit tests reveal that this class of phenotypic models, which has been successful so far in interpreting experimental data, is a plausible in only three of nine biological systems. More precisely, although Fisher's model was able to explain several statistical properties of the landscapes-including the mean and SD of selection and epistasis coefficients-it was often unable to explain the full structure of fitness landscapes.

The effect of gene interactions on the long-term response to selection

The role of gene interactions in the evolutionary process has long been controversial. Although some argue that they are not of importance, because most variation is additive, others claim that their effect in the long term can be substantial. Here, we focus on the long-term effects of genetic interactions under directional selection assuming no mutation or dominance, and that epistasis is symmetrical overall. We ask by how much the mean of a complex trait can be increased by selection and analyze two extreme regimes, in which either drift or selection dominate the dynamics of allele frequencies. In both scenarios, epistatic interactions affect the long-term response to selection by modulating the additive genetic variance. When drift dominates, we extend Robertson's [Robertson A (1960)Proc R Soc Lond B Biol Sci153(951):234-249] argument to show that, for any form of epistasis, the total response of a haploid population is proportional to the initial total genotypic variance. In contrast, the total response of a diploid population is increased by epistasis, for a given initial genotypic variance. When selection dominates, we show that the total selection response can only be increased by epistasis when some initially deleterious alleles become favored as the genetic background changes. We find a simple approximation for this effect and show that, in this regime, it is the structure of the genotype-phenotype map that matters and not the variance components of the population.

Reduction in the cumulative effect of stress-induced inbreeding depression due to intragenerational purging in Drosophila melanogaster

Environmental stress generally exacerbates the harmful effects of inbreeding and it has been proposed that this could be exploited in purging deleterious alleles from threatened inbred populations. However, understanding what factors contribute to variability in the strength of inbreeding depression (ID) observed across adverse environmental conditions remains a challenge. Here, we examined how the nature and timing of stress affects ID and the potential for purging using inbred and outbred Drosophila melanogaster larvae exposed to biotic (larval competition, bacteria infection) and abiotic (ethanol, heat) stressors compared with unstressed controls. ID was measured during (larval survival) and after (male mating success) stress exposure. The level of stress imposed by each stressor was approximately equal, averaging a 42% reduction in outbred larval survival relative to controls. All stressors induced on average the same ID, causing a threefold increase in lethal equivalents for larval survival relative to controls. However, stress-induced ID in larval success was followed by a 30% reduction in ID in mating success of surviving males. We propose that this fitness recovery is due to 'intragenerational purging' whereby fitness correlations facilitate stress-induced purging that increases the average fitness of survivors in later life history stages. For biotic stressors, post-stress reductions in ID are consistent with intragenerational purging, whereas for abiotic stressors, there appeared to be an interaction between purging and stress-induced physiological damage. For all stressors, there was no net effect of stress on lifetime ID compared with unstressed controls, undermining the prediction that stress enhances the effectiveness of population-level purging across generations.

Estimating the mutation load in human genomes

Next-generation sequencing technology has facilitated the discovery of millions of genetic variants in human genomes. A sizeable fraction of these variants are predicted to be deleterious. Here, we review the pattern of deleterious alleles as ascertained in genome sequencing data sets and ask whether human populations differ in their predicted burden of deleterious alleles - a phenomenon known as mutation load. We discuss three demographic models that are predicted to affect mutation load and relate these models to the evidence (or the lack thereof) for variation in the efficacy of purifying selection in diverse human genomes. We also emphasize why accurate estimation of mutation load depends on assumptions regarding the distribution of dominance and selection coefficients - quantities that remain poorly characterized for current genomic data sets.

Mutation load under additive fitness effects

Under the traditional mutation load model based on multiplicative fitness effects, the load in a population is 1−e^−U, where U is the genomic deleterious mutation rate. Because this load becomes high under large U, synergistic epistasis has been proposed as one possible means of reducing the load. However, experiments on model organisms attempting to detect synergistic epistasis have often focused on a quadratic fitness model, with the resulting general conclusion being that epistasis is neither common nor strong. Here, I present a model of additive fitness effects and show that, unlike multiplicative effects, the equilibrium frequency of an allele under additivity is dependent on the average absolute fitness of the population. The additive model then results in a load of U/(U +1), which is much lower than 1−e^−U for large U. Numerical iterations demonstrate that this analytic derivation holds as a good approximation under biologically relevant values of selection coefficients and U. Additionally, regressions onto Drosophila mutation accumulation data suggest that the common method of inferring epistasis by detecting large quadratic terms from regressions is not always necessary, as the additive model fits the data well and results in synergistic epistasis. Furthermore, the additive model gives a much larger reduction in load than the quadratic model when predicted from the same data, indicating that it is important to consider this additive model in addition to the quadratic model when inferring epistasis from mutation accumulation data.

The effect of parasites on sex differences in selection

The life history strategies of males and females are often divergent, creating the potential for sex differences in selection. Deleterious mutations may be subject to stronger selection in males, owing to sexual selection, which can improve the mean fitness of females and reduce mutation load in sexual populations. However, sex differences in selection might also maintain sexually antagonistic genetic variation, creating a sexual conflict load. The overall impact of separate sexes on fitness is unclear, but the net effect is likely to be positive when there is a large sex difference in selection against deleterious mutations. Parasites can also have sex-specific effects on fitness, and there is evidence that parasites can intensify the fitness consequences of deleterious mutations. Using lines that accumulated mutations for over 60 generations, we studied the effect of the pathogenic bacterium Pseudomonas aeruginosa on sex differences in selection in the fruit fly Drosophila melanogaster. Pseudomonas infection increased the sex difference in selection, but may also have weakened the intersexual correlation for fitness. Our results suggest that parasites may increase the benefits of sexual selection.

Properties of selected mutations and genotypic landscapes under Fisher's geometric model

The fitness landscape-the mapping between genotypes and fitness-determines properties of the process of adaptation. Several small genotypic fitness landscapes have recently been built by selecting a handful of beneficial mutations and measuring fitness of all combinations of these mutations. Here, we generate several testable predictions for the properties of these small genotypic landscapes under Fisher's geometric model of adaptation. When the ancestral strain is far from the fitness optimum, we analytically compute the fitness effect of selected mutations and their epistatic interactions. Epistasis may be negative or positive on average depending on the distance of the ancestral genotype to the optimum and whether mutations were independently selected, or coselected in an adaptive walk. Simulations show that genotypic landscapes built from Fisher's model are very close to an additive landscape when the ancestral strain is far from the optimum. However, when it is close to the optimum, a large diversity of landscape with substantial roughness and sign epistasis emerged. Strikingly, small genotypic landscapes built from several replicate adaptive walks on the same underlying landscape were highly variable, suggesting that several realizations of small genotypic landscapes are needed to gain information about the underlying architecture of the fitness landscape.

What can we learn from fitness landscapes?

A combinatorially complete data set consists of studies of all possible combinations of a set of mutant sites in a gene or mutant alleles in a genome. Among the most robust conclusions from these studies is that epistasis between beneficial mutations often shows a pattern of diminishing returns, in which favorable mutations are less fit when combined than would be expected. Another robust inference is that the number of adaptive evolutionary paths is often limited to a relatively small fraction of the theoretical possibilities, owing largely to sign epistasis requiring evolutionary steps that would entail a decrease in fitness. Here we summarize these and other results while also examining issues that remain unresolved and future directions that seem promising.

Empirical fitness landscapes and the predictability of evolution

The genotype-fitness map (that is, the fitness landscape) is a key determinant of evolution, yet it has mostly been used as a superficial metaphor because we know little about its structure. This is now changing, as real fitness landscapes are being analysed by constructing genotypes with all possible combinations of small sets of mutations observed in phylogenies or in evolution experiments. In turn, these first glimpses of empirical fitness landscapes inspire theoretical analyses of the predictability of evolution. Here, we review these recent empirical and theoretical developments, identify methodological issues and organizing principles, and discuss possibilities to develop more realistic fitness landscape models.

Causes and consequences of genetic background effects illuminated by integrative genomic analysis

The phenotypic consequences of individual mutations are modulated by the wild-type genetic background in which they occur. Although such background dependence is widely observed, we do not know whether general patterns across species and traits exist or about the mechanisms underlying it. We also lack knowledge on how mutations interact with genetic background to influence gene expression and how this in turn mediates mutant phenotypes. Furthermore, how genetic background influences patterns of epistasis remains unclear. To investigate the genetic basis and genomic consequences of genetic background dependence of the scalloped(E3) allele on the Drosophila melanogaster wing, we generated multiple novel genome-level datasets from a mapping-by-introgression experiment and a tagged RNA gene expression dataset. In addition we used whole genome resequencing of the parental lines-two commonly used laboratory strains-to predict polymorphic transcription factor binding sites for SD. We integrated these data with previously published genomic datasets from expression microarrays and a modifier mutation screen. By searching for genes showing a congruent signal across multiple datasets, we were able to identify a robust set of candidate loci contributing to the background-dependent effects of mutations in sd. We also show that the majority of background-dependent modifiers previously reported are caused by higher-order epistasis, not quantitative noncomplementation. These findings provide a useful foundation for more detailed investigations of genetic background dependence in this system, and this approach is likely to prove useful in exploring the genetic basis of other traits as well.

The effects of sexual selection on life-history traits: an experimental study on guppies

Sexual selection is often prevented during captive breeding in order to maximize effective population size and retain genetic diversity. However, enforcing monogamy and thereby preventing sexual selection may affect population fitness either negatively by preventing the purging of deleterious mutations or positively by reducing sexual conflicts. To better understand the effect of sexual selection on the fitness of small populations, we compared components of female fitness and the expression of male secondary sexual characters in 19 experimental populations of guppies (Poecilia reticulata) maintained under polygamous or monogamous mating regimes over nine generations. In order to generate treatments that solely differed by their level of sexual selection, the middle-class neighbourhood breeding design was enforced in the monogamous populations, while in the polygamous populations, all females contributed similarly to the next generation with one male and one female offspring. This experimental design allowed potential sexual conflicts to increase in the polygamous populations because selection could not operate on adult-female traits. Clutch size and offspring survival showed a weak decline from generation to generation but did not differ among treatments. Offspring size, however, declined across generations, but more in monogamous than polygamous populations. By generation eight, orange- and black-spot areas were larger in males from the polygamous treatment, but these differences were not statistically significant. Overall, these results suggest that neither sexual conflict nor the purging of deleterious mutation had important effects on the fitness of our experimental populations. However, only few generations of enforced monogamy in a benign environment were sufficient to negatively affect offspring size, a trait potentially crucial for survival in the wild. Sexual selection may therefore, under certain circumstances, be beneficial over enforced monogamy during captive breeding.

Should evolutionary geneticists worry about higher-order epistasis?

Natural selection drives evolving populations up the fitness landscape, the projection from nucleotide sequence space to organismal reproductive success. While it has long been appreciated that topographic complexities on fitness landscapes can arise only as a consequence of epistatic interactions between mutations, evolutionary genetics has mainly focused on epistasis between pairs of mutations. Here we propose a generalization to the classical population genetic treatment of pairwise epistasis that yields expressions for epistasis among arbitrary subsets of mutations of all orders (pairwise, three-way, etc.). Our approach reveals substantial higher-order epistasis in almost every published fitness landscape. Furthermore we demonstrate that higher-order epistasis is critically important in two systems we know best. We conclude that higher-order epistasis deserves empirical and theoretical attention from evolutionary geneticists.

The conditional nature of genetic interactions: the consequences of wild-type backgrounds on mutational interactions in a genome-wide modifier screen

The phenotypic outcome of a mutation cannot be simply mapped onto the underlying DNA variant. Instead, the phenotype is a function of the allele, the genetic background in which it occurs and the environment where the mutational effects are expressed. While the influence of genetic background on the expressivity of individual mutations is recognized, its consequences on the interactions between genes, or the genetic network they form, is largely unknown. The description of genetic networks is essential for much of biology; yet if, and how, the topologies of such networks are influenced by background is unknown. Furthermore, a comprehensive examination of the background dependent nature of genetic interactions may lead to identification of novel modifiers of biological processes. Previous work in Drosophila melanogaster demonstrated that wild-type genetic background influences the effects of an allele of scalloped (sd), with respect to both its principal consequence on wing development and its interactions with a mutation in optomotor blind. In this study we address whether the background dependence of mutational interactions is a general property of genetic systems by performing a genome wide dominant modifier screen of the sd(E3) allele in two wild-type genetic backgrounds using molecularly defined deletions. We demonstrate that ~74% of all modifiers of the sd(E3) phenotype are background-dependent due in part to differential sensitivity to genetic perturbation. These background dependent interactions include some with qualitative differences in the phenotypic outcome, as well as instances of sign epistasis. This suggests that genetic interactions are often contingent on genetic background, with flexibility in genetic networks due to segregating variation in populations. Such background dependent effects can substantially alter conclusions about how genes influence biological processes, the potential for genetic screens in alternative wild-type backgrounds identifying new loci that contribute to trait expression, and the inferences of the topology of genetic networks.

Investigating the genetic load of an emblematic invasive species: the case of the invasive harlequin ladybird Harmonia axyridis

Introduction events can lead to admixture between genetically differentiated populations and bottlenecks in population size. These processes can alter the adaptive potential of invasive species by shaping genetic variation, but more importantly, they can also directly affect mean population fitness either increasing it or decreasing it. Which outcome is observed depends on the structure of the genetic load of the species. The ladybird Harmonia axyridis is a good example of invasive species where introduced populations have gone through admixture and bottleneck events. We used laboratory experiments to manipulate the relatedness among H. axyridis parental individuals to assess the possibility for heterosis or outbreeding depression in F₁ generation offspring for two traits related to fitness (lifetime performance and generation time). We found that inter-populations crosses had no major impact on the lifetime performance of the offspring produced by individuals from either native or invasive populations. Significant outbreeding depression was observed only for crosses between native populations for generation time. The absence of observed heterosis is indicative of a low occurrence of fixed deleterious mutations within both the native and invasive populations of H. axyridis. The observed deterioration of fitness in native inter-population crosses most likely results from genetic incompatibilities between native genomic backgrounds. We discuss the implications of these results for the structure of genetic load in H. axyridis in the light of the available information regarding the introduction history of this species.

Male-biased fitness effects of spontaneous mutations in Drosophila melanogaster

In populations with males and females, sexual selection may often represent a major component of overall selection. Sexual selection could act to eliminate deleterious alleles in concert with other forms of selection, thereby improving the fitness of sexual populations. Alternatively, the divergent reproductive strategies of the sexes could promote the maintenance of sexually antagonistic variation, causing sexual populations to be less fit. The net impact of sexual selection on fitness is not well understood, due in part to limited data on the sex-specific effects of spontaneous mutations on total fitness. Using a set of mutation accumulation lines of Drosophila melanogaster, we found that mutations were deleterious in both sexes and had larger effects on fitness in males than in females. This pattern is expected to reduce the mutation load of sexual females and promote the maintenance of sexual reproduction.

Sex enhances adaptation by unlinking beneficial from detrimental mutations in experimental yeast populations

BACKGROUND

The maintenance of sexuality is a classic problem in evolutionary biology because it is a less efficient mode of reproduction compared with asexuality; however, many organisms are sexual. Theoretical work suggests sex facilitates natural selection, and experimental data support this. However, there are fewer experimental studies that have attempted to determine the mechanisms underlying the advantage of sex. Two main classes of hypotheses have been proposed to explain its advantage: detrimental mutation clearance and beneficial mutation accumulation. Here we attempt to experimentally differentiate between these two classes by evolving Saccharomyces cerevisiae populations that differ only in their ability to undergo sex, and also manipulate mutation rate. We cannot manipulate the types of mutation that occur, but instead propagate populations in both stressful and permissive environments and assume that the extent of detrimental mutation clearance and beneficial mutation incorporation differs between them.

RESULTS

After 300 mitotic generations interspersed with 11 rounds of sex we found there was no change or difference in fitness between sexuals and asexuals propagated in the permissive environment, regardless of mutation rate. Sex conferred a greater extent of adaptation in the stressful environment, and wild-type and elevated mutation rate sexual populations adapted equivalently. However, the asexual populations with an elevated mutation rate appeared more retarded in their extent of adaptation compared to asexual wild-type populations.

CONCLUSIONS

Sex provided no advantage in the permissive environment where beneficial mutations were rare. We could not evaluate if sex functioned to clear detrimental mutations more effectively or not here as no additional fitness load was observed in the mutator populations. However, in the stressful environment, where detrimental mutations were likely of more consequence, and where beneficial mutations were apparent, sex provided an advantage. In the stressful environment asexuals were increasingly constrained in their extent of adaptation with increasing mutation rate. Sex appeared to facilitate adaptation not just by more rapidly combining beneficial mutations, but also by unlinking beneficial from detrimental mutations: sex allowed selection to operate on both types of mutations more effectively compared to asexual populations.

Exploring the complexity of the HIV-1 fitness landscape

Although fitness landscapes are central to evolutionary theory, so far no biologically realistic examples for large-scale fitness landscapes have been described. Most currently available biological examples are restricted to very few loci or alleles and therefore do not capture the high dimensionality characteristic of real fitness landscapes. Here we analyze large-scale fitness landscapes that are based on predictive models for in vitro replicative fitness of HIV-1. We find that these landscapes are characterized by large correlation lengths, considerable neutrality, and high ruggedness and that these properties depend only weakly on whether fitness is measured in the absence or presence of different antiretrovirals. Accordingly, adaptive processes on these landscapes depend sensitively on the initial conditions. While the relative extent to which mutations affect fitness on their own (main effects) or in combination with other mutations (epistasis) is a strong determinant of these properties, the fitness landscape of HIV-1 is considerably less rugged, less neutral, and more correlated than expected from the distribution of main effects and epistatic interactions alone. Overall this study confirms theoretical conjectures about the complexity of biological fitness landscapes and the importance of the high dimensionality of the genetic space in which adaptation takes place.

Evolutionary adaptation can be understood as populations moving uphill on landscapes, in which height corresponds to evolutionary fitness. Although such fitness landscapes are central to evolutionary theory, there is currently a lack of biologically realistic examples. Here we analyze large-scale fitness landscapes derived from in vitro fitness measurements of HIV-1. We find that these landscapes are very rugged and that, accordingly, adaptive processes on these landscapes depend sensitively on the initial conditions. Moreover, the landscapes contain large networks along which fitness changes only minimally. While the relative extent to which mutations affect fitness on their own or in combination with other mutations is a strong determinant of these properties, the fitness landscape of HIV-1 is considerably less rugged than expected from the individual and pair-wise effects of mutations. Overall this study confirms theoretical conjectures about the complexity of biological fitness landscapes and the importance of the high dimensionality of the genetic space in which adaptation takes place.

The causes of epistasis

Since Bateson's discovery that genes can suppress the phenotypic effects of other genes, gene interactions-called epistasis-have been the topic of a vast research effort. Systems and developmental biologists study epistasis to understand the genotype-phenotype map, whereas evolutionary biologists recognize the fundamental importance of epistasis for evolution. Depending on its form, epistasis may lead to divergence and speciation, provide evolutionary benefits to sex and affect the robustness and evolvability of organisms. That epistasis can itself be shaped by evolution has only recently been realized. Here, we review the empirical pattern of epistasis, and some of the factors that may affect the form and extent of epistasis. Based on their divergent consequences, we distinguish between interactions with or without mean effect, and those affecting the magnitude of fitness effects or their sign. Empirical work has begun to quantify epistasis in multiple dimensions in the context of metabolic and fitness landscape models. We discuss possible proximate causes (such as protein function and metabolic networks) and ultimate factors (including mutation, recombination, and the importance of natural selection and genetic drift). We conclude that, in general, pleiotropy is an important prerequisite for epistasis, and that epistasis may evolve as an adaptive or intrinsic consequence of changes in genetic robustness and evolvability.

Quantitative trait locus analysis of stage-specific inbreeding depression in the Pacific oyster Crassostrea gigas

Inbreeding depression and genetic load have been widely observed, but their genetic basis and effects on fitness during the life cycle remain poorly understood, especially for marine animals with high fecundity and high, early mortality (type-III survivorship). A high load of recessive mutations was previously inferred for the Pacific oyster Crassostrea gigas, from massive distortions of zygotic, marker segregation ratios in F(2) families. However, the number, genomic location, and stage-specific onset of mutations affecting viability have not been thoroughly investigated. Here, we again report massive distortions of microsatellite-marker segregation ratios in two F(2) hybrid families, but we now locate the causative deleterious mutations, using a quantitative trait locus (QTL) interval-mapping model, and we characterize their mode of gene action. We find 14-15 viability QTL (vQTL) in the two families. Genotypic frequencies at vQTL generally suggest selection against recessive or partially recessive alleles, supporting the dominance theory of inbreeding depression. No epistasis was detected among vQTL, so unlinked vQTL presumably have independent effects on survival. For the first time, we track segregation ratios of vQTL-linked markers through the life cycle, to determine their stage-specific expression. Almost all vQTL are absent in the earliest life stages examined, confirming zygotic viability selection; vQTL are predominantly expressed before the juvenile stage (90%), mostly at metamorphosis (50%). We estimate that, altogether, selection on vQTL caused 96% mortality in these families, accounting for nearly all of the actual mortality. Thus, genetic load causes substantial mortality in inbred Pacific oysters, particularly during metamorphosis, a critical developmental transition warranting further investigation.

The population genetics of X-autosome synthetic lethals and steriles

Epistatic interactions are widespread, and many of these interactions involve combinations of alleles at different loci that are deleterious when present in the same individual. The average genetic environment of sex-linked genes differs from that of autosomal genes, suggesting that the population genetics of interacting X-linked and autosomal alleles may be complex. Using both analytical theory and computer simulations, we analyzed the evolutionary trajectories and mutation-selection balance conditions for X-autosome synthetic lethals and steriles. Allele frequencies follow a set of fundamental trajectories, and incompatible alleles are able to segregate at much higher frequencies than single-locus expectations. Equilibria exist, and they can involve fixation of either autosomal or X-linked alleles. The exact equilibrium depends on whether synthetic alleles are dominant or recessive and whether fitness effects are seen in males, females, or both sexes. When single-locus fitness effects and synthetic incompatibilities are both present, population dynamics depend on the dominance of alleles and historical contingency (i.e., whether X-linked or autosomal mutations occur first). Recessive synthetic lethality can result in high-frequency X-linked alleles, and dominant synthetic lethality can result in high-frequency autosomal alleles. Many X-autosome incompatibilities in natural populations may be cryptic, appearing to be single-locus effects because one locus is fixed. We also discuss the implications of these findings with respect to standing genetic variation and the origins of Haldane's rule.

Dietary stress does not strengthen selection against single deleterious mutations in Drosophila melanogaster

Stress is generally thought to increase the strength of selection, although empirical results are mixed and general conclusions are difficult because data are limited. Here we compare the fitness effects of nine independent recessive mutations in Drosophila melanogaster in a high- and low-dietary-stress environment, estimating the strength of selection on these mutations arising from both a competitive measure of male reproductive success and productivity (female fecundity and the subsequent survival to adulthood of her offspring). The effect of stress on male reproductive success has not been addressed previously for individual loci and is of particular interest with respect to the alignment of natural and sexual selection. Our results do not support the hypothesis that stress increases the efficacy of selection arising from either fitness component. Results concerning the alignment of natural and sexual selection were mixed, although data are limited. In the low-stress environment, selection on mating success and productivity were concordant for five of nine mutations (four out of four when restricted to those with significant or near-significant productivity effects), whereas in the high-stress environment, selection aligned for seven of nine mutations (two out of two when restricted to those having significant productivity effects). General conclusions as to the effects of stress on the strength of selection and the alignment of natural and sexual selection await data from additional mutations, fitness components and stressors.

Reducing mutation load through sexual selection on males

Mutation load is a key parameter in evolutionary theories, but relatively little empirical information exists on the mutation load of populations, or the elimination of this load through selection. We manipulated the opportunity for sexual selection within a mutation accumulation divergence experiment to determine how sexual selection on males affected the accumulation of mutations contributing to sexual and nonsexual fitness. Sexual selection prevented the accumulation of mutations affecting male mating success, the target trait, as well as reducing mutation load on productivity, a nonsexual fitness component. Mutational correlations between mating success and productivity (estimated in the absence of sexual selection) were positive. Sexual selection significantly reduced these fitness component correlations. Male mating success significantly diverged between sexual selection treatments, consistent with the fixation of genetic differences. However, the rank of the treatments was not consistent across assays, indicating that the mutational effects on mating success were conditional on biotic and abiotic context. Our experiment suggests that greater insight into the genetic targets of natural and sexual selection can be gained by focusing on mutational rather than standing genetic variation, and on the behavior of trait variances rather than means.

Sex-specific effects of inbreeding in wild-caught Drosophila melanogaster under benign and stressful conditions

In animal populations, sib mating is often the primary source of inbreeding depression (ID). We used recently wild-caught Drosophila melanogaster to test whether such ID is amplified by environmental stress and, in males, by sexual selection. We also investigated whether increased ID because of stress (increased larval competition) persisted beyond the stressed stage and whether the effects of stress and sexual selection interacted. Sib mating resulted in substantial cumulative fitness losses (egg to adult reproduction) of 50% (benign) and 73% (stressed). Stress increased ID during the larval period (23% vs. 63%), but not during post-stress reproductive stages (36% vs. 31%), indicating larval stress may have purged some adult genetic load (although ID was uncorrelated across stages). Sexual selection exacerbated inbreeding depression, with inbred male offspring suffering a higher reproductive cost than females, independent of stress (57% vs. 14% benign, 49% vs. 11% stress).

Local adaptation does not always predict high mating success

The hypothesis that adaptation to local environments can increase mating success was tested using ten replicate lines of Drosophila melanogaster adapted either to 16 degrees C or to 25 degrees C. Competitive mating trials at both temperatures were performed with males taken from a pair of lines, one adapted to each temperature. There was no average increase in mating success for males adapted to the local environment. Although one pair of lines showed the expected pattern, another pair showed the reverse pattern. More data are needed on this hypothesis, preferably with lines that have more strongly adapted to local environments.

Evolutionary dynamics, epistatic interactions, and biological information

We investigate a definition of biological information that connects population genetics with the tools of information theory by focusing on the distribution of genotypes found in a population. Previous research has treated loci as non-interacting by making specific approximations in the calculation of information-theoretic quantities. We expand earlier mathematical forms to include epistasis, or interactions between mutations at all pairs of loci. Application of our improved measure of biological information to evolution on two-locus, two-allele fitness landscapes demonstrates that mutual information between loci reflects epistatic interaction of mutations. Finally, we consider four-locus, two-allele fitness landscapes with modular structure. As modular interactions are inherently epistatic, we demonstrate that our refined approximation provides insight into the underlying structure of these non-trivial fitness landscapes.