The evolution of sex: empirical insights into the roles of epistasis and drift

Despite many years of theoretical and experimental work, the explanation for why sex is so common as a reproductive strategy continues to resist understanding. Recent empirical work has addressed key questions in this field, especially regarding rates of mutation accumulation in sexual and asexual organisms, and the roles of negative epistasis and drift as sources of adaptive constraint in asexually reproducing organisms. At the same time, new ideas about the evolution of sexual recombination are being tested, including intriguing suggestions of an important interplay between sex and genetic architecture, which indicate that sex and recombination could have affected their own evolution.     

Genetic Architecture Promotes the Evolution and Maintenance of Cooperation

When cooperation has a direct cost and an indirect benefit, a selfish behavior is more likely to be selected for than an altruistic one. Kin and group selection do provide evolutionary explanations for the stability of cooperation in nature, but we still lack the full understanding of the genomic mechanisms that can prevent cheater invasion. In our study we used Aevol, an agent-based, in silico genomic platform to evolve populations of digital organisms that compete, reproduce, and cooperate by secreting a public good for tens of thousands of generations. We found that cooperating individuals may share a phenotype, defined as the amount of public good produced, but have very different abilities to resist cheater invasion. To understand the underlying genetic differences between cooperator types, we performed bio-inspired genomics analyses of our digital organisms by recording and comparing the locations of metabolic and secretion genes, as well as the relevant promoters and terminators. Association between metabolic and secretion genes (promoter sharing, overlap via frame shift or sense-antisense encoding) was characteristic for populations with robust cooperation and was more likely to evolve when secretion was costly. In mutational analysis experiments, we demonstrated the potential evolutionary consequences of the genetic association by performing a large number of mutations and measuring their phenotypic and fitness effects. The non-cooperating mutants arising from the individuals with genetic association were more likely to have metabolic deleterious mutations that eventually lead to selection eliminating such mutants from the population due to the accompanying fitness decrease. Effectively, cooperation evolved to be protected and robust to mutations through entangled genetic architecture. Our results confirm the importance of second-order selection on evolutionary outcomes, uncover an important genetic mechanism for the evolution and maintenance of cooperation, and suggest promising methods for preventing gene loss in synthetically engineered organisms.     

Recombination and the evolution of coordinated phenotypic expression in a frequency-dependent game

A long standing question in evolutionary biology concerns the maintenance of adaptive combinations of traits in the presence of recombination. This problem may be solved if positive epistasis selects for reducing the rate of recombination between such traits, but this requires sufficiently strong epistasis. Here we use a model that we developed previously to analyze a frequency-dependent strategy game in asexual populations, to study how adaptive combinations of traits may be maintained in the presence of recombination when epistasis is too weak to select for genetic linkage. Previously, in the asexual case, our model demonstrated the evolution of adaptive associations between social foraging strategies and learning rules. We verify that these adaptive associations, which are represented by different two-locus haplotypes, can easily be broken by genetic recombination. We also confirm that a modifier allele that reduces the rate of recombination fails to evolve (due to weak epistasis). However, we find that under the same conditions of weak epistasis, there is an alternative mechanism that allows an association between traits to evolve. This is based on a genetic switch that responds to the presence of one social foraging allele by activating one of the two alternative learning alleles that are carried by all individuals. We suggest that such coordinated phenotypic expression by genetic switches offers a general and robust mechanism for the evolution of adaptive combinations of traits in the presence of recombination.     

The Evolution of Epistasis and the Advantage of Recombination in Populations of Bacteriophage T4

Experiments reported here test two hypotheses about the evolution of recombination: first, the Fisher-Muller concept that sexual organisms respond to selection more rapidly than do asexual ones, and second, that epistasis is more likely to evolve in the absence of recombination. Populations of bacteriophage T4 were selected by the drug proflavine in discrete generations and the change in mean population fitness was monitored. Three separate selection series yielded results supporting the Fisher-Muller hypothesis. The amount of epistasis evolved was measured by partitioning the T4 map into regions and comparing the sum of the proflavine resistances of each region with the resistance of the whole. Significantly more interactions were found in phage isolated from the populations with lower total recombination than in those from populations with higher recombination. The degree to which these experiments fit preconceived notions about natural selection suggests that microorganisms may be advantageously used in other population genetics experiments.     

A developmental perspective on the evolution of sexual size dimorphism of a moth

Males and females of almost all organisms exhibit sexual differences in body size, a phenomenon called sexual size dimorphism (SSD). How the sexes evolve to be different sizes, despite sharing the same genes that control growth and development, and hence a common genetic architecture, has remained elusive. Here, we show that the genetic architecture (heritabilities and genetic correlations) of the physiological mechanism that regulates size during the last stage of larval development of a moth, differs between the sexes, and thus probably facilitates, rather than hinders, the evolution of SSD. We further show that the endocrine system plays a critical role in generating SSD. Our results demonstrate that knowledge of the genetic architecture underlying the physiological process during development that ultimately produces SSD in adults can elucidate how males and females of organisms evolve to be of different sizes.     

RUNAWAY SEXUAL SELECTION LEADS TO GOOD GENES

Mate choice and sexual displays are widespread in nature, but their evolutionary benefits remain controversial. Theory predicts these traits can be favored by runaway sexual selection, in which preference and display reinforce one another due to genetic correlation; or by good genes benefits, in which mate choice is advantageous because extreme displays indicate a well‐adapted genotype. However, these hypotheses are not mutually exclusive, and the adaptive benefits underlying mate choice can themselves evolve. In particular, examining how and why sexual displays become indicators of good genes is challenging in natural systems. Here, we use experimental evolution in “digital organisms” to demonstrate the origins of condition‐dependent indicator displays following their spread due to a runaway process. Surprisingly, handicap‐like costs are not necessary for displays to become indicators of male viability. Instead, a pleiotropic genetic architecture underlies both displays and viability. Runaway sexual selection and good genes benefits should thus be viewed as interacting mechanisms that reinforce one another.     

Admixed sexual and facultatively asexual aphid lineages at mating sites

Cyclically parthenogenetic organisms may have facultative asexual counterparts. Such organisms, including aphids, are therefore interesting models for the study of ecological and genetic interactions between lineages differing in reproductive mode. Earlier studies on aphids have revealed major differences in the genetic outcomes of populations that are possibly resulting mostly either from sexual or from asexual reproduction. Besides, notable gene flow between sexual and asexual derivatives has been suspected, which could lead to the emergence of new asexual lineages. The present study examines the interplay between these lineages and is based on analyses of population structure of individuals that may contribute to the pool of sexual reproductive forms in the host alternating aphid Rhopalosiphum padi. Using a Bayesian assignment method, we first show that the sexual forms of R. padi on mating sites encompass two genetically distinct clusters of individuals in the western part of France. The first cluster included unique genotypes of sexual lineages, while the second cluster included facultatively asexual lineages in numerous copies, the reproductive mode of the two clusters being confirmed by reference clones. Sexual reproductive forms produced by sexual and facultatively asexual lineages are thus admixed at mating sites which gives a large opportunity for the two clusters to mate with each other. Nevertheless, this study also highlights, as previously demonstrated, that the two clusters retained high genetic differentiation. Possible explanations for the inferred limited genetic exchanges are advanced in the discussion, but further dedicated investigations are required to solve this paradox.     

Fitness-dependent mutation rates in finite populations

Mutation rate may be condition dependent, whereby individuals in poor condition, perhaps from high mutation load, have higher mutation rates than individuals in good condition. Agrawal (J. Evol. Biol.15, 2002, 1004) explored the basic properties of fitness-dependent mutation rate (FDMR) in infinite populations and reported some heuristic results for finite populations. The key parameter governing how infinite populations evolve under FDMR is the curvature (k) of the relationship between fitness and mutation rate. We extend Agrawal's analysis to finite populations and consider dominance and epistasis. In finite populations, the probability of long-term existence depends on k. In sexual populations, positive curvature leads to low equilibrium mutation rate, whereas negative curvature results in high mutation rate. In asexual populations, negative curvature results in rapid extinction via 'mutational meltdown', whereas positive curvature sometimes allows persistence. We speculate that fitness-dependent mutation rate may provide the conditions for genetic architecture to diverge between sexual and asexual taxa.     

Impact of epistasis and pleiotropy on evolutionary adaptation

Evolutionary adaptation is often likened to climbing a hill or peak. While this process is simple for fitness landscapes where mutations are independent, the interaction between mutations (epistasis) as well as mutations at loci that affect more than one trait (pleiotropy) are crucial in complex and realistic fitness landscapes. We investigate the impact of epistasis and pleiotropy on adaptive evolution by studying the evolution of a population of asexual haploid organisms (haplotypes) in a model of N interacting loci, where each locus interacts with K other loci. We use a quantitative measure of the magnitude of epistatic interactions between substitutions, and find that it is an increasing function of K. When haplotypes adapt at high mutation rates, more epistatic pairs of substitutions are observed on the line of descent than expected. The highest fitness is attained in landscapes with an intermediate amount of ruggedness that balance the higher fitness potential of interacting genes with their concomitant decreased evolvability. Our findings imply that the synergism between loci that interact epistatically is crucial for evolving genetic modules with high fitness, while too much ruggedness stalls the adaptive process.     

Why reproduce sexually?

There is reason to believe that intense selection, such that only a small minority at the top of the fitness distribution has any appreciable chance of survival, can sometimes give sexual reproduction an immediate (one-life-cycle) advantage over asexual. The advantage must be great enough to balance the 50% loss of genetic material in meiosis.One model shows the advantage to be frequency-dependent in life cycles in which there are several asexual generations and one sexual. The observed frequency of sexual reproduction in such a life cycle is explained as an evolutionary equilibrium by this model. In another model the optimum frequency of asexual reproduction drops to zero as fecundity and competition increase. This explains the exclusively sexual reproduction of such fecund organism as elms and oysters. Once lost, asexual reproduction may be difficult to evolve secondarily. This explains the presence of such exclusively sexual, low-fecundity organisms as the higher vertebrates.     

EVOLUTION OF TRADE-OFFS BETWEEN SEXUAL AND ASEXUAL PHASES AND THE ROLE OF REPRODUCTIVE PLASTICITY IN THE GENETIC ARCHITECTURE OF APHID LIFE HISTORIES

Life‐history theory postulates that evolution is constrained by trade‐offs (i.e., negative genetic correlations) among traits that contribute to fitness. However, in organisms with complex life cycles, trade‐offs may drastically differ between phases, putatively leading to different evolutionary trajectories. Here, we tested this possibility by examining changes in life‐history traits in an aphid species that alternates asexual and sexual reproduction in its life cycle. The quantitative genetics of reproductive and dispersal traits was studied in 23 lineages (genotypes) of the bird cherry‐oat aphid Rhopalosiphum padi, during both the sexual and asexual phases, which were induced experimentally under specific environmental conditions. We found large and significant heritabilities (broad‐sense) for all traits and several negative genetic correlations between traits (trade‐offs), which are related to reproduction (i.e., numbers of the various sexual or asexual morphs) or dispersal (i.e., numbers of winged or wingless morphs). These results suggest that R. padi exhibits lineage specialization both in reproductive and dispersal strategies. In addition, we found important differences in the structure of genetic variance–covariance matrices ( G ) between phases. These differences were due to two large, negative genetic correlations detected during the asexual phase only: (1) between fecundity and age at maturity and (2) between the production of wingless and winged parthenogenetic females. We propose that this differential expression in genetic architecture results from a reallocation scheme during the asexual phase, when sexual morphs are not produced. We also found significant G × E interaction and nonsignificant genetic correlations across phases, indicating that genotypes could respond independently to selection in each phase. Our results reveal a rather unique situation in which the same population and even the same genotypes express different genetic (co)variation under different environmental conditions, driven by optimal resource allocation criteria.     

Genetic variation of geminiviruses: comparison between sexual and asexual host plant populations

One of the most promising hypotheses for the evolution of sex is that sexual reproduction is advantageous because it increases the rate of adaptive evolution in response to parasites. To investigate this advantage of sex, we compared genetic variation of geminiviruses infecting sexual and asexual populations of Eupatorium (Asteraceae). The infection frequency was 37.5% in the sexual population and 87.8% in the asexual population. The lower infection frequency in the sexual population might be the result of higher genetic diversity of host plants. If geminiviruses have diverged to counter defence systems of genetically variable hosts, genetic diversity of viruses is expected to be higher in sexual host populations than in asexual host populations. To test this expectation, we used single-strand conformation polymorphism (SSCP) analysis to examine genetic diversity of the geminiviruses in a DNA region containing the open-reading frame (ORF) C4 gene, which is known to function as a host range determinant. As predicted, higher genetic diversity of viruses was observed in the sexual population: three SSCP types were found in the asexual population while six types were found in the sexual population. Sequencing of the polymerase chain reaction (PCR) products revealed further genetic diversity. Phylogenetic analysis of the sequences showed that the SSCP types belonged to four different clades. Several SSCP types from the same clade were found in the sexual population, whereas the asexual population included only one SSCP type from each clade. Amino acid replacements of ORF C4 are suggested to be accelerated in the sexual population. This evidence supports the hypothesis that sexual reproduction is advantageous as a defence against epidemic disease.     

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.     

Network models of TEM β-lactamase mutations coevolving under antibiotic selection show modular structure and anticipate evolutionary trajectories

Understanding how novel functions evolve (genetic adaptation) is a critical goal of evolutionary biology. Among asexual organisms, genetic adaptation involves multiple mutations that frequently interact in a non-linear fashion (epistasis). Non-linear interactions pose a formidable challenge for the computational prediction of mutation effects. Here we use the recent evolution of β-lactamase under antibiotic selection as a model for genetic adaptation. We build a network of coevolving residues (possible functional interactions), in which nodes are mutant residue positions and links represent two positions found mutated together in the same sequence. Most often these pairs occur in the setting of more complex mutants. Focusing on extended-spectrum resistant sequences, we use network-theoretical tools to identify triple mutant trajectories of likely special significance for adaptation. We extrapolate evolutionary paths (n = 3) that increase resistance and that are longer than the units used to build the network (n = 2). These paths consist of a limited number of residue positions and are enriched for known triple mutant combinations that increase cefotaxime resistance. We find that the pairs of residues used to build the network frequently decrease resistance compared to their corresponding singlets. This is a surprising result, given that their coevolution suggests a selective advantage. Thus, β-lactamase adaptation is highly epistatic. Our method can identify triplets that increase resistance despite the underlying rugged fitness landscape and has the unique ability to make predictions by placing each mutant residue position in its functional context. Our approach requires only sequence information, sufficient genetic diversity, and discrete selective pressures. Thus, it can be used to analyze recent evolutionary events, where coevolution analysis methods that use phylogeny or statistical coupling are not possible. Improving our ability to assess evolutionary trajectories will help predict the evolution of clinically relevant genes and aid in protein design.     

Genetic variation and asexual reproduction in the facultatively parthenogenetic cockroach Nauphoeta cinerea: implications for the evolution of sex

Asexual reproduction could offer up to a two‐fold fitness advantage over sexual reproduction, yet higher organisms usually reproduce sexually. Even in facultatively parthenogenetic species, where both sexual and asexual reproduction is sometimes possible, asexual reproduction is rare. Thus, the debate over the evolution of sex has focused on ecological and mutation‐elimination advantages of sex. An alternative explanation for the predominance of sex is that it is difficult for an organism to accomplish asexual reproduction once sexual reproduction has evolved. Difficulty in returning to asexuality could reflect developmental or genetic constraints. Here, we investigate the role of genetic factors in limiting asexual reproduction in Nauphoeta cinerea, an African cockroach with facultative parthenogenesis that nearly always reproduces sexually. We show that when N. cinerea females do reproduce asexually, offspring are genetically identical to their mothers. However, asexual reproduction is limited to a nonrandom subset of the genotypes in the population. Only females that have a high level of heterozygosity are capable of parthenogenetic reproduction and there is a strong familial influence on the ability to reproduce parthenogenetically. Although the mechanism by which genetic variation facilitates asexual reproduction is unknown, we suggest that heterosis may facilitate the switch from producing haploid meiotic eggs to diploid, essentially mitotic, eggs.     

Independently evolving species in asexual bdelloid rotifers

Asexuals are an important test case for theories of why species exist. If asexual clades displayed the same pattern of discrete variation as sexual clades, this would challenge the traditional view that sex is necessary for diversification into species. However, critical evidence has been lacking: all putative examples have involved organisms with recent or ongoing histories of recombination and have relied on visual interpretation of patterns of genetic and phenotypic variation rather than on formal tests of alternative evolutionary scenarios. Here we show that a classic asexual clade, the bdelloid rotifers, has diversified into distinct evolutionary species. Intensive sampling of the genus Rotaria reveals the presence of well-separated genetic clusters indicative of independent evolution. Moreover, combined genetic and morphological analyses reveal divergent selection in feeding morphology, indicative of niche divergence. Some of the morphologically coherent groups experiencing divergent selection contain several genetic clusters, in common with findings of cryptic species in sexual organisms. Our results show that the main causes of speciation in sexual organisms, population isolation and divergent selection, have the same qualitative effects in an asexual clade. The study also demonstrates how combined molecular and morphological analyses can shed new light on the evolutionary nature of species.     

A framework for detecting and characterizing genetic background-dependent imprinting effects

Genomic imprinting, where the effects of alleles depend on their parent-of-origin, can be an important component of the genetic architecture of complex traits. Although there has been a rapidly increasing number of studies of genetic architecture that have examined imprinting effects, none have examined whether imprinting effects depend on genetic background. Such effects are critical for the evolution of genomic imprinting because they allow the imprinting state of a locus to evolve as a function of genetic background. Here we develop a two-locus model of epistasis that includes epistatic interactions involving imprinting effects and apply this model to scan the mouse genome for loci that modulate the imprinting effects of quantitative trait loci (QTL). The inclusion of imprinting leads to nine orthogonal forms of epistasis, five of which do not appear in the usual two-locus decomposition of epistasis. Each form represents a change in the imprinting status of one locus across different classes of genotypes at the other locus. Our genome scan identified two different locus pairs that show complex patterns of epistasis, where the imprinting effect at one locus changes across genetic backgrounds at the other locus. Thus, our model provides a framework for the detection of genetic background-dependent imprinting effects that should provide insights into the background dependence and evolution of genomic imprinting. Our application of the model to a genome scan supports this assertion by identifying pairs of loci that show reciprocal changes in their imprinting status as the background provided by the other locus changes.     

The balance between sexual and asexual reproduction in plants living in variable environments

The balance between sexual and asexual propagule production is studied in an evolutionary model where plants produce the two kinds of propagules in genetically determined proportions. The male function of plants producing asexual propagules can be varied, and the sexual and asexual propagules carry different probabilities to turn into new reproductive individuals. These fitnesses may vary over years. The evolution of the population’s reproductive system is studied assuming modifier alleles with small effects. In this setting a balanced, mixed reproductive system can evolve, but only if the difference in fitness between the sexual and asexual propagules varies over years. When the two kinds of propagules are very similar to each other, as is often the case with sexual and asexual seed formation, evolution will tend towards a state dominated by the one or the other reproductive system.     

Understanding Evolutionary Potential in Virtual CPU Instruction Set Architectures

We investigate fundamental decisions in the design of instruction set architectures for linear genetic programs that are used as both model systems in evolutionary biology and underlying solution representations in evolutionary computation. We subjected digital organisms with each tested architecture to seven different computational environments designed to present a range of evolutionary challenges. Our goal was to engineer a general purpose architecture that would be effective under a broad range of evolutionary conditions. We evaluated six different types of architectural features for the virtual CPUs: (1) genetic flexibility: we allowed digital organisms to more precisely modify the function of genetic instructions, (2) memory: we provided an increased number of registers in the virtual CPUs, (3) decoupled sensors and actuators: we separated input and output operations to enable greater control over data flow. We also tested a variety of methods to regulate expression: (4) explicit labels that allow programs to dynamically refer to specific genome positions, (5) position-relative search instructions, and (6) multiple new flow control instructions, including conditionals and jumps. Each of these features also adds complication to the instruction set and risks slowing evolution due to epistatic interactions. Two features (multiple argument specification and separated I/O) demonstrated substantial improvements in the majority of test environments, along with versions of each of the remaining architecture modifications that show significant improvements in multiple environments. However, some tested modifications were detrimental, though most exhibit no systematic effects on evolutionary potential, highlighting the robustness of digital evolution. Combined, these observations enhance our understanding of how instruction architecture impacts evolutionary potential, enabling the creation of architectures that support more rapid evolution of complex solutions to a broad range of challenges.     

Sexual Development of Cellular Slime Molds

Macrocyst formation represents sexual development of cellular slime molds and begins with fusion between cells of compatible mating types. Homothallic and heterothallic strains as well as bisexual and asexual ones have been described. Macrocyst development requires certain environmental conditions such as darkness and excessive humidity. Sexual cell fusion has been analyzed at a molecular level in Dictyostelium discoideum , and several cell surface proteins related to it have been identified. Some of them are common to both mating types, while others are specific to one or other type. The involvement of cell-surface carbohydrates has also been suggested, though direct evidence for this is still lacking. Macrocyst formation is regulated by diffusible, pheromone like substances. Genetic studies on sexual development are scarce, probably because no suitable mutants have been available. However, several asexual mutants, as well as antibody and nucleotide probes, have recently been obtained, so mechanisms of sexual cell fusion may be understood in the near future. Considering the unique phylogenical position of cellular slime molds, analysis of sexual development in these organisms should contribute to the understanding of the mechanism and evolution of sexual reproduction systems.