An Evolving Genetic Architecture Interacts with Hill-Robertson Interference to Determine the Benefit of Sex

Beyond the (geometric) mean: stochastic models undermine deterministic predictions of bet hedger evolution

Bet hedging is a ubiquitous strategy for risk reduction in the face of unpredictable environmental change where a lineage lowers its variance in fitness across environments at the expense of also lowering its arithmetic mean fitness. Classically, the benefit of bet hedging has been quantified using geometric mean fitness (GMF); bet hedging is expected to evolve if and only if it has a higher GMF than the wild-type. We build upon previous research on the effect of incorporating stochasticity in phenotypic distribution, environment, and reproduction to investigate the extent to which these sources of stochasticity will impact the evolution of real-world bet hedging traits. We utilize both individual-based simulations and Markov chain numerics to demonstrate that modeling stochasticity can alter the sign of selection for the bet hedger compared to deterministic predictions. We find that bet hedging can be deleterious at small population sizes and beneficial at larger population sizes. This non-monotonic dependence of the sign of selection on population size, known as sign inversion, exists across parameter space for both conservative and diversified bet hedgers. We apply our model to published data of bet hedging strategies to show that sign inversion exists for biologically relevant parameters in two study systems: Papaver dubium, an annual poppy with variable germination phenology, and Salmonella typhimurium, a pathogenic bacteria that exhibits antibiotic persistence. Taken together, our results suggest that GMF is not enough to predict when bet hedging is adaptive.

Genetic drift promotes and recombination hinders speciation on holey fitness landscapes

Dobzhansky and Muller proposed a general mechanism through which microevolution, the substitution of alleles within populations, can cause the evolution of reproductive isolation between populations and, therefore, macroevolution. As allopatric populations diverge, many combinations of alleles differing between them have not been tested by natural selection and may thus be incompatible. Such genetic incompatibilities often cause low fitness in hybrids between species. Furthermore, the number of incompatibilities grows with the genetic distance between diverging populations. However, what determines the rate and pattern of accumulation of incompatibilities remains unclear. We investigate this question by simulating evolution on holey fitness landscapes on which genetic incompatibilities can be identified unambiguously. We find that genetic incompatibilities accumulate more slowly among genetically robust populations and identify two determinants of the accumulation rate: recombination rate and population size. In large populations with abundant genetic variation, recombination selects for increased genetic robustness and, consequently, incompatibilities accumulate more slowly. In small populations, genetic drift interferes with this process and promotes the accumulation of genetic incompatibilities. Our results suggest a novel mechanism by which genetic drift promotes and recombination hinders speciation.

The advantage of sex: Reinserting fluctuating selection in the pluralist approach

The advantage of sex, and its fixation in some clades and species all over the eukaryote tree of life, is considered an evolutionary enigma, especially regarding its assumed two-fold cost. Several likely hypotheses have been proposed such as (1) a better response to the negative frequency-dependent selection imposed by the “Red Queen” hypothesis; (2) the competition between siblings induced by the Tangled Bank hypothesis; (3) the existence of genetic and of (4) ecological factors that can diminish the cost of sex to less than the standard assumed two-fold; and (5) a better maintenance of genetic diversity and its resulting phenotypic variation, providing a selective advantage in randomly fluctuating environments. While these hypotheses have mostly been studied separately, they can also act simultaneously. This was advocated by several studies which presented a pluralist point of view. Only three among the five causes cited above were considered yet in such a framework: the Red Queen hypothesis, the Tangled Bank and the genetic factors lowering the cost of sex. We thus simulated the evolution of a finite mutating population undergoing negative frequency-dependent selection on phenotypes and a two-fold (or less) cost of sexuality, experiencing randomly fluctuating selection along generations. The individuals inherited their reproductive modes, either clonal or sexual. We found that exclusive sexuality begins to fix in populations exposed to environmental variation that exceeds the width of one ecological niche (twice the standard deviation of a Gaussian response to environment). This threshold was lowered by increasing negative frequency-dependent selection and when reducing the two-fold cost of sex. It contributes advocating that the different processes involved in a short-term advantage of sex and recombination can act in combination to favor the fixation of sexual reproduction in populations.

Mitonuclear Mismatch is Associated With Increased Male Frequency, Outcrossing, and Male Sperm Size in Experimentally-Evolved C. elegans

We provide a partial test of the mitonuclear sex hypothesis with the first controlled study of how male frequencies and rates of outcrossing evolve in response to mitonuclear mismatch by allowing replicate lineages of C. elegans nematodes containing either mitochondrial or nuclear mutations of electron transport chain (ETC) genes to evolve under three sexual systems: facultatively outcrossing (wildtype), obligately selfing, and obligately outcrossing. Among facultatively outcrossing lines, we found evolution of increased male frequency in at least one replicate line of all four ETC mutant backgrounds tested—nuclear isp-1, mitochondrial cox-1 and ctb-1, and an isp-1 IV; ctb-1M mitonuclear double mutant—and confirmed for a single line set (cox-1) that increased male frequency also resulted in successful outcrossing. We previously found the same result for lines evolved from another nuclear ETC mutant, gas-1. For several lines in the current experiment, however, male frequency declined to wildtype levels (near 0%) in later generations. Male frequency did not change in lines evolved from a wildtype control strain. Additional phenotypic assays of lines evolved from the mitochondrial cox-1 mutant indicated that evolution of high male frequency was accompanied by evolution of increased male sperm size and mating success with tester females, but that it did not translate into increased mating success with coevolved hermaphrodites. Rather, hermaphrodites’ self-crossed reproductive fitness increased, consistent with sexually antagonistic coevolution. In accordance with evolutionary theory, males and sexual outcrossing may be most beneficial to populations evolving from a state of low ancestral fitness (gas-1, as previously reported) and less beneficial or deleterious to those evolving from a state of higher ancestral fitness (cox-1). In support of this idea, the obligately outcrossing fog-2 V; cox-1 M lines exhibited no fitness evolution compared to their ancestor, while facultatively outcrossing lines showed slight upward evolution of fitness, and all but one of the obligately selfing xol-1 X; cox-1 M lines evolved substantially increased fitness—even beyond wildtype levels. This work provides a foundation to directly test the effect of reproductive mode on the evolutionary dynamics of mitonuclear genomes, as well as whether compensatory mutations (nuclear or mitochondrial) can rescue populations from mitochondrial dysfunction.

Using phenotypic plasticity to understand the structure and evolution of the genotype-phenotype map

Deciphering the genotype-phenotype map necessitates relating variation at the genetic level to variation at the phenotypic level. This endeavour is inherently limited by the availability of standing genetic variation, the rate of spontaneous mutation to novo genetic variants, and possible biases associated with induced mutagenesis. An interesting alternative is to instead rely on the environment as a source of variation. Many phenotypic traits change plastically in response to the environment, and these changes are generally underlain by changes in gene expression. Relating gene expression plasticity to the phenotypic plasticity of more integrated organismal traits thus provides useful information about which genes influence the development and expression of which traits, even in the absence of genetic variation. We here appraise the prospects and limits of such an environment-for-gene substitution for investigating the genotype-phenotype map. We review models of gene regulatory networks, and discuss the different ways in which they can incorporate the environment to mechanistically model phenotypic plasticity and its evolution. We suggest that substantial progress can be made in deciphering this genotype-environment-phenotype map, by connecting theory on gene regulatory network to empirical patterns of gene co-expression, and by more explicitly relating gene expression to the expression and development of phenotypes, both theoretically and empirically.

Seasonal changes in recombination characteristics in a natural population of Drosophila melanogaster

Environmental seasonality is a potent evolutionary force, capable of maintaining polymorphism, promoting phenotypic plasticity and causing bet-hedging. In Drosophila, environmental seasonality has been reported to affect life-history traits, tolerance to abiotic stressors and immunity. Oscillations in frequencies of alleles underlying fitness-related traits were also documented alongside SNPs across the genome. Here, we test for seasonal changes in two recombination characteristics, crossover rate and crossover interference, in a natural D. melanogaster population from India using morphological markers of the three major chromosomes. We show that winter flies, collected after the dry season, have significantly higher desiccation tolerance than their autumn counterparts. This difference proved to hold also for hybrids with three independent marker stocks, suggesting its genetic rather than plastic nature. Significant between-season changes are documented for crossover rate (in 9 of 13 studied intervals) and crossover interference (in four of eight studied pairs of intervals); both single and double crossovers were usually more frequent in the winter cohort. The winter flies also display weaker plasticity of both recombination characteristics to desiccation. We ascribe the observed differences to indirect selection on recombination caused by directional selection on desiccation tolerance. Our findings suggest that changes in recombination characteristics can arise even after a short period of seasonal adaptation (~8-10 generations).

Divergence times and plastid phylogenomics within the intron-rich order Erythropeltales (Compsopogonophyceae, Rhodophyta)

The advent of high-throughput sequencing (HTS) has allowed for the use of large numbers of coding regions to produce robust phylogenies. These phylogenies have been used to highlight relationships at ancient diversifications (subphyla, class) and highlight the evolution of plastid genome structure. The Erythropeltales are an order in the Compsopogonophyceae, a group with unusual plastid genomes but with low taxon sampling. We use HTS to produce near complete plastid genomes of all genera, and multiple species within some genera, to produce robust phylogenies to investigate character evolution, dating of divergence in the group, and plastid organization, including intron patterns. Our results produce a fully supported phylogeny of the genera in the Erythropeltales and suggest that morphologies (upright versus crustose) have evolved multiple times. Our dated phylogeny also indicates that the order is very old (~800 Ma), with diversification occurring after the ice ages of the Cryogenian period (750-635 Ma). Plastid gene order is congruent with phylogenetic relationships and suggests that genome architecture does not change often. Our data also highlight the abundance of introns in the plastid genomes of this order. We also produce a nearly complete plastid genome of Tsunamia transpacifica (Stylonematophyceae) to add to the taxon sampling of genomes of this class. The use of plastid genomes clearly produces robust phylogenetic relationships that can be used to infer evolutionary events, and increased taxon sampling, especially in less well-known red algal groups, will provide additional insights into their evolution.

Selective Interference and the Evolution of Sex

Selection acts upon genes linked together on chromosomes. This physical connection reduces the efficiency by which selection can act because, in the absence of sex, alleles must rise and fall together in frequency with the genome in which they are found. This selective interference underlies such phenomena as clonal interference and Muller’s Ratchet and is broadly termed Hill-Robertson interference. In this review, I examine the potential for selective interference to account for the evolution and maintenance of sex, discussing the positive and negative evidence from both theoretical and empirical studies, and highlight the gaps that remain.

The evolutionary advantage of fitness-dependent recombination in diploids: A deterministic mutation-selection balance model

Recombination's omnipresence in nature is one of the most intriguing problems in evolutionary biology. The question of why recombination exhibits certain general features is no less interesting than that of why it exists at all. One such feature is recombination's fitness dependence (FD). The so far developed population genetics models have focused on the evolution of FD recombination mainly in haploids, although the empirical evidence for this phenomenon comes mostly from diploids. Using numerical analysis of modifier models for infinite panmictic populations, we show here that FD recombination can be evolutionarily advantageous in diploids subjected to purifying selection. We ascribe this advantage to the differential rate of disruption of lower- versus higher-fitness genotypes, which can be manifested in selected systems with at least three loci. We also show that if the modifier is linked to such selected system, it can additionally benefit from modifying this linkage in a fitness-dependent manner. The revealed evolutionary advantage of FD recombination appeared robust to crossover interference within the selected system, either positive or negative. Remarkably, FD recombination was often favored in situations where any constant nonzero recombination was evolutionarily disfavored, implying a relaxation of the rather strict constraints on major parameters (e.g., selection intensity and epistasis) required for the evolutionary advantage of nonzero recombination formulated by classical models.

Population structure promotes the evolution of costly sex in artificial gene networks

We build on previous observations that Hill-Robertson interference generates an advantage of sex that, in structured populations, can be large enough to explain the evolutionary maintenance of costly sex. We employed a gene network model that explicitly incorporates interactions between genes. Mutations in the gene networks have variable effects that depend on the genetic background in which they appear. Consequently, our simulations include two costs of sex-recombination and migration loads-that were missing from previous studies of the evolution of costly sex. Our results suggest a critical role for population structure that lies in its ability to align the long- and short-term advantages of sex. We show that the addition of population structure favored the evolution of sex by disproportionately decreasing the equilibrium mean fitness of asexual populations, primarily by increasing the strength of Muller's Ratchet. Population structure also increased the ability of the short-term advantage of sex to counter the primary limit to the evolution of sex in the gene network model-recombination load. On the other hand, highly structured populations experienced migration load in the form of Dobzhansky-Muller incompatibilities, decreasing the effective rate of migration between demes and, consequently, accelerating the accumulation of drift load in the sexual populations.

Sex and Mitonuclear Adaptation in Experimental Caenorhabditis elegans Populations

To reveal phenotypic and functional genomic patterns of mitonuclear adaptation, a laboratory adaptation study with Caenorhabditis elegans nematodes was conducted in which independently evolving lines were initiated from a low-fitness mitochondrial electron transport chain (ETC) mutant, gas-1 Following 60 generations of evolution in large population sizes with competition for food resources, two distinct classes of lines representing different degrees of adaptive response emerged: a low-fitness class that exhibited minimal or no improvement compared to the gas-1 mutant ancestor, and a high-fitness class containing lines that exhibited partial recovery of wild-type fitness. Many lines that achieved higher reproductive and competitive fitness levels were also noted to evolve high frequencies of males during the experiment, consistent with adaptation in these lines having been facilitated by outcrossing. Whole-genome sequencing and analysis revealed an enrichment of mutations in loci that occur in a gas-1-centric region of the C. elegans interactome and could be classified into a small number of functional genomic categories. A highly nonrandom pattern of mitochondrial DNA mutation was observed within high-fitness gas-1 lines, with parallel fixations of nonsynonymous base substitutions within genes encoding NADH dehydrogenase subunits I and VI. These mitochondrial gene products reside within ETC complex I alongside the nuclear-encoded GAS-1 protein, suggesting that rapid adaptation of select gas-1 recovery lines was driven by fixation of compensatory mitochondrial mutations.

Recombination drives the evolution of mutational robustness

Recombination can impose fitness costs as beneficial parental combinations of alleles are broken apart, a phenomenon known as recombination load. Computational models suggest that populations may evolve a reduced recombination load by reducing either the likelihood of recombination events (bring interacting loci in physical proximity) or the strength of interactions between loci (make loci more independent of one another). We review evidence for each of these possibilities and their consequences for the genotype-fitness relationship. In particular, we expect that reducing interaction strengths between loci will lead to genomes that are also robust to mutational perturbations, but reducing recombination rates alone will not. We note that both mechanisms most likely played a role in the evolution of extant populations, and that both can result in the frequently-observed pattern of physical linkage between interacting loci.

Synergy from reproductive division of labor and genetic complexity drive the evolution of sex

Computer experiments that mirror the evolutionary dynamics of sexual and asexual organisms as they occur in nature were used to test features proposed to explain the evolution of sexual recombination. Results show that this evolution is better described as a network of interactions between possible sexual forms, including diploidy, thelytoky, facultative sex, assortation, bisexuality, and division of labor between the sexes, rather than a simple transition from parthenogenesis to sexual recombination. Diploidy was shown to be fundamental for the evolution of sex; bisexual reproduction emerged only among anisogamic diploids with a synergistic division of reproductive labor; and facultative sex was more likely to evolve among haploids practicing assortative mating. Looking at the evolution of sex as a complex system through individual-based simulations explains better the diversity of sexual strategies known to exist in nature, compared to classical analytical models.

Sign of selection on mutation rate modifiers depends on population size

The influence of population size (N) on natural selection acting on alleles that affect fitness has been understood for almost a century. As N declines, genetic drift overwhelms selection and alleles with direct fitness effects are rendered neutral. Often, however, alleles experience so-called indirect selection, meaning they affect not the fitness of an individual but the fitness distribution of its offspring. Some of the best-studied examples of indirect selection include alleles that modify aspects of the genetic system such as recombination and mutation rates. Here, we use analytics, simulations, and experimental populations of Saccharomyces cerevisiae to examine the influence of N on indirect selection acting on alleles that increase the genomic mutation rate (mutators). Mutators experience indirect selection via genomic associations with beneficial and deleterious mutations they generate. We show that, as N declines, indirect selection driven by linked beneficial mutations is overpowered by drift before drift can neutralize the cost of the deleterious load. As a result, mutators transition from being favored by indirect selection in large populations to being disfavored as N declines. This surprising phenomenon of sign inversion in selective effect demonstrates that indirect selection on mutators exhibits a profound and qualitatively distinct dependence on N.

The role of recombination in evolutionary adaptation of Escherichia coli to a novel nutrient

The benefits and detriments of recombination for adaptive evolution have been studied both theoretically and experimentally, with conflicting predictions and observations. Most pertinent experiments examine recombination's effects in an unchanging environment and do not study its genomewide effects. Here, we evolved six replicate populations of either highly recombining R⁺ or lowly recombining R^- E. coli strains in a changing environment, by introducing the novel nutrients L-arabinose or indole into the environment. The experiment's ancestral strains are not viable on these nutrients, but 130 generations of adaptive evolution were sufficient to render them viable. Recombination conferred a more pronounced advantage to populations adapting to indole. To study the genomic changes associated with this advantage, we sequenced the genomes of 384 clones isolated from selected replicates at the end of the experiment. These genomes harbour complex changes that range from point mutations to large-scale DNA amplifications. Among several candidate adaptive mutations, those in the tryptophanase regulator tnaC stand out, because the tna operon in which it resides has a known role in indole metabolism. One of the highly recombining populations also shows a significant excess of large-scale segmental DNA amplifications that include the tna operon. This lineage also shows a unique and potentially adaptive combination of point mutations and DNA amplifications that may have originated independently from one another, to be joined later by recombination. Our data illustrate that the advantages of recombination for adaptive evolution strongly depend on the environment and that they can be associated with complex genomic changes.