Selection for recombination in small populations

Immunocytological analysis of meiotic recombination in two anole lizards (Squamata, Dactyloidae)

Although the evolutionary importance of meiotic recombination is not disputed, the significance of interspecies differences in the recombination rates and recombination landscapes remains under-appreciated. Recombination rates and distribution of chiasmata have been examined cytologically in many mammalian species, whereas data on other vertebrates are scarce. Immunolocalization of the protein of the synaptonemal complex (SYCP3), centromere proteins and the mismatch-repair protein MLH1 was used, which is associated with the most common type of recombination nodules, to analyze the pattern of meiotic recombination in the male of two species of iguanian lizards, Anolis carolinensis Voigt, 1832 and Deiroptyx coelestinus (Cope, 1862). These species are separated by a relatively long evolutionary history although they retain the ancestral iguanian karyotype. In both species similar and extremely uneven distributions of MLH1 foci along the macrochromosome bivalents were detected: approximately 90% of crossovers were located at the distal 20% of the chromosome arm length. Almost total suppression of recombination in the intermediate and proximal regions of the chromosome arms contradicts the hypothesis that "homogenous recombination" is responsible for the low variation in GC content across the anole genome. It also leads to strong linkage disequilibrium between the genes located in these regions, which may benefit conservation of co-adaptive gene arrays responsible for the ecological adaptations of the anoles.

Higher rates of sex evolve during adaptation to more complex environments

A leading hypothesis for the evolutionary maintenance of sexual reproduction proposes that sex is advantageous because it facilitates adaptation. Changes in the environment stimulate adaptation but not all changes are equivalent; a change may occur along one or multiple environmental dimensions. In two evolution experiments with the facultatively sexual rotifer Brachionus calyciflorus, we test how environmental complexity affects the evolution of sex by adapting replicate populations to various environments that differ from the original along one, two, or three environmental dimensions. Three different estimates of fitness (growth, lifetime reproduction, and population density) confirmed that populations adapted to their new environment. Growth measures revealed an intriguing cost of complex adaptations: populations that adapted to more complex environments lost greater amounts of fitness in the original environment. Furthermore, both experiments showed that B. calyciflorus became more sexual when adapting to a greater number of environmental dimensions. Common garden experiments confirmed that observed changes in sex were heritable. As environments in nature are inherently complex these findings help explain why sex is maintained in natural populations.

Extreme genomic erosion after recurrent demographic bottlenecks in the highly endangered Iberian lynx

BACKGROUND

Genomic studies of endangered species provide insights into their evolution and demographic history, reveal patterns of genomic erosion that might limit their viability, and offer tools for their effective conservation. The Iberian lynx (Lynx pardinus) is the most endangered felid and a unique example of a species on the brink of extinction.

RESULTS

We generate the first annotated draft of the Iberian lynx genome and carry out genome-based analyses of lynx demography, evolution, and population genetics. We identify a series of severe population bottlenecks in the history of the Iberian lynx that predate its known demographic decline during the 20th century and have greatly impacted its genome evolution. We observe drastically reduced rates of weak-to-strong substitutions associated with GC-biased gene conversion and increased rates of fixation of transposable elements. We also find multiple signatures of genetic erosion in the two remnant Iberian lynx populations, including a high frequency of potentially deleterious variants and substitutions, as well as the lowest genome-wide genetic diversity reported so far in any species.

CONCLUSIONS

The genomic features observed in the Iberian lynx genome may hamper short- and long-term viability through reduced fitness and adaptive potential. The knowledge and resources developed in this study will boost the research on felid evolution and conservation genomics and will benefit the ongoing conservation and management of this emblematic species.

The Nonstationary Dynamics of Fitness Distributions: Asexual Model with Epistasis and Standing Variation

Various models describe asexual evolution by mutation, selection, and drift. Some focus directly on fitness, typically modeling drift but ignoring or simplifying both epistasis and the distribution of mutation effects (traveling wave models). Others follow the dynamics of quantitative traits determining fitness (Fisher's geometric model), imposing a complex but fixed form of mutation effects and epistasis, and often ignoring drift. In all cases, predictions are typically obtained in high or low mutation rate limits and for long-term stationary regimes, thus losing information on transient behaviors and the effect of initial conditions. Here, we connect fitness-based and trait-based models into a single framework, and seek explicit solutions even away from stationarity. The expected fitness distribution is followed over time via its cumulant generating function, using a deterministic approximation that neglects drift. In several cases, explicit trajectories for the full fitness distribution are obtained for arbitrary mutation rates and standing variance. For nonepistatic mutations, especially with beneficial mutations, this approximation fails over the long term but captures the early dynamics, thus complementing stationary stochastic predictions. The approximation also handles several diminishing returns epistasis models (e.g., with an optimal genotype); it can be applied at and away from equilibrium. General results arise at equilibrium, where fitness distributions display a "phase transition" with mutation rate. Beyond this phase transition, in Fisher's geometric model, the full trajectory of fitness and trait distributions takes a simple form; robust to the details of the mutant phenotype distribution. Analytical arguments are explored regarding why and when the deterministic approximation applies.

Genomics of Natural Populations: How Differentially Expressed Genes Shape the Evolution of Chromosomal Inversions in Drosophila pseudoobscura

Chromosomal rearrangements can shape the structure of genetic variation in the genome directly through alteration of genes at breakpoints or indirectly by holding combinations of genetic variants together due to reduced recombination. The third chromosome of Drosophila pseudoobscura is a model system to test hypotheses about how rearrangements are established in populations because its third chromosome is polymorphic for >30 gene arrangements that were generated by a series of overlapping inversion mutations. Circumstantial evidence has suggested that these gene arrangements are selected. Despite the expected homogenizing effects of extensive gene flow, the frequencies of arrangements form gradients or clines in nature, which have been stable since the system was first described >80 years ago. Furthermore, multiple arrangements exist at appreciable frequencies across several ecological niches providing the opportunity for heterokaryotypes to form. In this study, we tested whether genes are differentially expressed among chromosome arrangements in first instar larvae, adult females and males. In addition, we asked whether transcriptional patterns in heterokaryotypes are dominant, semidominant, overdominant, or underdominant. We find evidence for a significant abundance of differentially expressed genes across the inverted regions of the third chromosome, including an enrichment of genes involved in sensory perception for males. We find the majority of loci show additivity in heterokaryotypes. Our results suggest that multiple genes have expression differences among arrangements that were either captured by the original inversion mutation or accumulated after it reached polymorphic frequencies, providing a potential source of genetic variation for selection to act upon. These data suggest that the inversions are favored because of their indirect effect of recombination suppression that has held different combinations of differentially expressed genes together in the various gene arrangement backgrounds.

Effect of manipulating recombination rates on response to selection in livestock breeding programs

Background

In this work, we performed simulations to explore the potential of manipulating recombination rates to increase response to selection in livestock breeding programs.

Methods

We carried out ten replicates of several scenarios that followed a common overall structure but differed in the average rate of recombination along the genome (expressed as the length of a chromosome in Morgan), the genetic architecture of the trait under selection, and the selection intensity under truncation selection (expressed as the proportion of males selected). Recombination rates were defined by simulating nine different chromosome lengths: 0.10, 0.25, 0.50, 1, 2, 5, 10, 15 and 20 Morgan, respectively. One Morgan was considered to be the typical chromosome length for current livestock species. The genetic architecture was defined by the number of quantitative trait variants (QTV) that affected the trait under selection. Either a large (10,000) or a small (1000 or 500) number of QTV was simulated. Finally, the proportions of males selected under truncation selection as sires for the next generation were equal to 1.2, 2.4, 5, or 10 %.

Results

Increasing recombination rate increased the overall response to selection and decreased the loss of genetic variance. The difference in cumulative response between low and high recombination rates increased over generations. At low recombination rates, cumulative response to selection tended to asymptote sooner and the genetic variance was completely eroded. If the trait under selection was affected by few QTV, differences between low and high recombination rates still existed, but the selection limit was reached at all rates of recombination.

Conclusions

Higher recombination rates can enhance the efficiency of breeding programs to turn genetic variation into response to selection. However, to increase response to selection significantly, the recombination rate would need to be increased 10- or 20-fold. The biological feasibility and consequences of such large increases in recombination rates are unknown.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0221-1) contains supplementary material, which is available to authorized users.

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

Sex is ubiquitous in the natural world, but the nature of its benefits remains controversial. Previous studies have suggested that a major advantage of sex is its ability to eliminate interference between selection on linked mutations, a phenomenon known as Hill-Robertson interference. However, those studies may have missed both important advantages and important disadvantages of sexual reproduction because they did not allow the distributions of mutational effects and interactions (i.e., the genetic architecture) to evolve. Here we investigate how Hill-Robertson interference interacts with an evolving genetic architecture to affect the evolutionary origin and maintenance of sex by simulating evolution in populations of artificial gene networks. We observed a long-term advantage of sex-equilibrium mean fitness of sexual populations exceeded that of asexual populations-that did not depend on population size. We also observed a short-term advantage of sex-sexual modifier mutations readily invaded asexual populations-that increased with population size, as was observed in previous studies. We show that the long- and short-term advantages of sex were both determined by differences between sexual and asexual populations in the evolutionary dynamics of two properties of the genetic architecture: the deleterious mutation rate ([Formula: see text]) and recombination load ([Formula: see text]). These differences resulted from a combination of selection to minimize [Formula: see text] which is experienced only by sexuals, and Hill-Robertson interference experienced primarily by asexuals. In contrast to the previous studies, in which Hill-Robertson interference had only a direct impact on the fitness advantages of sex, the impact of Hill-Robertson interference in our simulations was mediated additionally by an indirect impact on the efficiency with which selection acted to reduce [Formula: see text].

Conserved Genetic Architecture Underlying Individual Recombination Rate Variation in a Wild Population of Soay Sheep (Ovis aries)

Meiotic recombination breaks down linkage disequilibrium (LD) and forms new haplotypes, meaning that it is an important driver of diversity in eukaryotic genomes. Understanding the causes of variation in recombination rate is important in interpreting and predicting evolutionary phenomena and in understanding the potential of a population to respond to selection. However, despite attention in model systems, there remains little data on how recombination rate varies at the individual level in natural populations. Here we used extensive pedigree and high-density SNP information in a wild population of Soay sheep (Ovis aries) to investigate the genetic architecture of individual autosomal recombination rates. Individual rates were high relative to other mammal systems and were higher in males than in females (autosomal map lengths of 3748 and 2860 cM, respectively). The heritability of autosomal recombination rate was low but significant in both sexes (h² = 0.16 and 0.12 in females and males, respectively). In females, 46.7% of the heritable variation was explained by a subtelomeric region on chromosome 6; a genome-wide association study showed the strongest associations at locus RNF212, with further associations observed at a nearby ∼374-kb region of complete LD containing three additional candidate loci, CPLX1, GAK, and PCGF3. A second region on chromosome 7 containing REC8 and RNF212B explained 26.2% of the heritable variation in recombination rate in both sexes. Comparative analyses with 40 other sheep breeds showed that haplotypes associated with recombination rates are both old and globally distributed. Both regions have been implicated in rate variation in mice, cattle, and humans, suggesting a common genetic architecture of recombination rate variation in mammals.

Evolvability of an Optimal Recombination Rate

Evolution and maintenance of genetic recombination and its relation to the mutational process is a long-standing, fundamental problem in evolutionary biology that is linked to the general problem of evolution of evolvability. We explored a stochastic model of the evolution of recombination using additive fitness and infinite allele assumptions but no assumptions on the sign or magnitude of the epistasis and the distribution of mutation effects. In this model, fluctuating negative epistasis and predominantly deleterious mutations arise naturally as a consequence of the additive fitness and a reservoir from which new alleles arrive with a fixed distribution of fitness effects. Analysis of the model revealed a nonmonotonic effect of recombination intensity on fitness, with an optimal recombination rate value which maximized fitness in steady state. The optimal recombination rate depended on the mutation rate and was evolvable, that is, subject to selection. The predictions of the model were compatible with the observations on the dependence between genome rearrangement rate and gene flux in microbial genomes.

Experimental evolution of recombination and crossover interference in Drosophila caused by directional selection for stress-related traits

BACKGROUND

Population genetics predicts that tight linkage between new and/or pre-existing beneficial and deleterious alleles should decrease the efficiency of natural selection in finite populations. By decoupling beneficial and deleterious alleles and facilitating the combination of beneficial alleles, recombination accelerates the formation of high-fitness genotypes. This may impose indirect selection for increased recombination. Despite the progress in theoretical understanding, interplay between recombination and selection remains a controversial issue in evolutionary biology. Even less satisfactory is the situation with crossover interference, which is a deviation of double-crossover frequency in a pair of adjacent intervals from the product of recombination rates in the two intervals expected on the assumption of crossover independence. Here, we report substantial changes in recombination and interference in three long-term directional selection experiments with Drosophila melanogaster: for desiccation (~50 generations), hypoxia, and hyperoxia tolerance (>200 generations each).

RESULTS

For all three experiments, we found a high interval-specific increase of recombination frequencies in selection lines (up to 40-50% per interval) compared to the control lines. We also discovered a profound effect of selection on interference as expressed by an increased frequency of double crossovers in selection lines. Our results show that changes in interference are not necessarily coupled with increased recombination.

CONCLUSIONS

Our results support the theoretical predictions that adaptation to a new environment can promote evolution toward higher recombination. Moreover, this is the first evidence of selection for different recombination-unrelated traits potentially leading, not only to evolution toward increased crossover rates, but also to changes in crossover interference, one of the fundamental features of recombination.

Meiosis evolves: adaptation to external and internal environments

306 I. 306 II. 307 III. 312 IV. 317 V. 318 319 References 319 SUMMARY: Meiosis is essential for the fertility of most eukaryotes and its structures and progression are conserved across kingdoms. Yet many of its core proteins show evidence of rapid or adaptive evolution. What drives the evolution of meiosis proteins? How can constrained meiotic processes be modified in response to challenges without compromising their essential functions? In surveying the literature, we found evidence of two especially potent challenges to meiotic chromosome segregation that probably necessitate adaptive evolutionary responses: whole-genome duplication and abiotic environment, especially temperature. Evolutionary solutions to both kinds of challenge are likely to involve modification of homologous recombination and synapsis, probably via adjustments of core structural components important in meiosis I. Synthesizing these findings with broader patterns of meiosis gene evolution suggests that the structural components of meiosis coevolve as adaptive modules that may change in primary sequence and function while maintaining three-dimensional structures and protein interactions. The often sharp divergence of these genes among species probably reflects periodic modification of entire multiprotein complexes driven by genomic or environmental changes. We suggest that the pressures that cause meiosis to evolve to maintain fertility may cause pleiotropic alterations of global crossover rates. We highlight several important areas for future research.

Recombining without Hotspots: A Comprehensive Evolutionary Portrait of Recombination in Two Closely Related Species of Drosophila

Meiotic recombination rate varies across the genome within and between individuals, populations, and species in virtually all taxa studied. In almost every species, this variation takes the form of discrete recombination hotspots, determined in some mammals by a protein called PRDM9. Hotspots and their determinants have a profound effect on the genomic landscape, and share certain features that extend across the tree of life. Drosophila, in contrast, are anomalous in their absence of hotspots, PRDM9, and other species-specific differences in the determination of recombination. To better understand the evolution of meiosis and general patterns of recombination across diverse taxa, we present a truly comprehensive portrait of recombination across time, combining recently published cross-based contemporary recombination estimates from each of two sister species with newly obtained linkage-disequilibrium-based historic estimates of recombination from both of these species. Using Drosophila pseudoobscura and Drosophila miranda as a model system, we compare recombination rate between species at multiple scales, and we suggest that Drosophila replicate the pattern seen in human-chimpanzee in which recombination rate is conserved at broad scales. We also find evidence of a species-wide recombination modifier(s), resulting in both a present and historic genome-wide elevation of recombination rates in D. miranda, and identify broad scale effects on recombination from the presence of an inversion. Finally, we reveal an unprecedented view of the distribution of recombination in D. pseudoobscura, illustrating patterns of linked selection and where recombination is taking place. Overall, by combining these estimation approaches, we highlight key similarities and differences in recombination between Drosophila and other organisms.

The evolution of sex: A new hypothesis based on mitochondrial mutational erosion: Mitochondrial mutational erosion in ancestral eukaryotes would favor the evolution of sex, harnessing nuclear recombination to optimize compensatory nuclear coadaptation

The evolution of sex in eukaryotes represents a paradox, given the "twofold" fitness cost it incurs. We hypothesize that the mutational dynamics of the mitochondrial genome would have favored the evolution of sexual reproduction. Mitochondrial DNA (mtDNA) exhibits a high-mutation rate across most eukaryote taxa, and several lines of evidence suggest that this high rate is an ancestral character. This seems inexplicable given that mtDNA-encoded genes underlie the expression of life's most salient functions, including energy conversion. We propose that negative metabolic effects linked to mitochondrial mutation accumulation would have invoked selection for sexual recombination between divergent host nuclear genomes in early eukaryote lineages. This would provide a mechanism by which recombinant host genotypes could be rapidly shuffled and screened for the presence of compensatory modifiers that offset mtDNA-induced harm. Under this hypothesis, recombination provides the genetic variation necessary for compensatory nuclear coadaptation to keep pace with mitochondrial mutation accumulation.

Sexual reproduction with variable mating systems can resist asexuality in a rock-paper-scissors dynamics

While sex can be advantageous for a lineage in the long term, we still lack an explanation for its maintenance with the twofold cost per generation. Here we model an infinite diploid population where two autosomal loci determine, respectively, the reproductive mode, sexual versus asexual and the mating system, polygynous (costly sex) versus monogamous (assuming equal contribution of parents to offspring, i.e. non-costly sex). We show that alleles for costly sex can spread when non-costly sexual modes buffer the interaction between asexual and costly sexual strategies, even without twofold benefit of recombination with respect to asexuality. The three interacting strategies have intransitive fitness relationships leading to a rock-paper-scissors dynamics, so that alleles for costly sex cannot be eliminated by asexuals in most situations throughout the parameter space. Our results indicate that sexual lineages with variable mating systems can resist the invasion of asexuals and allow for long-term effects to accumulate, thus providing a solution to the persisting theoretical question of why sex was not displaced by asexuality along evolution.

The evolution of obligate sex: the roles of sexual selection and recombination

The evolution of sex is one of the greatest mysteries in evolutionary biology. An even greater mystery is the evolution of obligate sex, particularly when competing with facultative sex and not with complete asexuality. Here, we develop a stochastic simulation of an obligate allele invading a facultative population, where males are subject to sexual selection. We identify a range of parameters where sexual selection can contribute to the evolution of obligate sex: Especially when the cost of sex is low, mutation rate is high, and the facultative individuals do not reproduce sexually very often. The advantage of obligate sex becomes larger in the absence of recombination. Surprisingly, obligate sex can take over even when the population has a lower mean fitness as a result. We show that this is due to the high success of obligate males that can compensate the cost of sex.

Influences of dominance and evolution of sex in finite diploid populations

Most eukaryotes reproduce sexually. Although the benefits of sex in diploids mainly stem from recombination and segregation, the relative effects of recombination and segregation are relatively less known. In this study, we adopt an infinite loci model to illustrate how dominance coefficient of mutations affects the above-mentioned genetic events. However, we assume mutational effects to be independent and also ignore the effects of epistasis within loci. Our simulations show that with different levels of dominance, segregation and recombination may play different roles. In particular, recombination more commonly has a major impact on the evolution of sex when deleterious mutations are partially recessive. In contrast, when deleterious mutations are dominant, segregation becomes more important than recombination, a finding that is consistent with previous studies stating that segregation, rather than recombination, is more likely to drive the evolution of sex. Moreover, beneficial mutations alone remarkably increases the effects of recombination. We also note that populations favor sexual reproduction when deleterious mutations become more dominant or beneficial mutations become more recessive. Overall, these results illustrate that the existence of dominance is an important mechanism that affects the evolution of sex.

Mitonuclear Ecology

Eukaryotes were born of a chimeric union between two prokaryotes--the progenitors of the mitochondrial and nuclear genomes. Early in eukaryote evolution, most mitochondrial genes were lost or transferred to the nucleus, but a core set of genes that code exclusively for products associated with the electron transport system remained in the mitochondrion. The products of these mitochondrial genes work in intimate association with the products of nuclear genes to enable oxidative phosphorylation and core energy production. The need for coadaptation, the challenge of cotransmission, and the possibility of genomic conflict between mitochondrial and nuclear genes have profound consequences for the ecology and evolution of eukaryotic life. An emerging interdisciplinary field that I call "mitonuclear ecology" is reassessing core concepts in evolutionary ecology including sexual reproduction, two sexes, sexual selection, adaptation, and speciation in light of the interactions of mitochondrial and nuclear genomes.

Evolution in changing environments: modifiers of mutation, recombination, and migration

The production and maintenance of genetic and phenotypic diversity under temporally fluctuating selection and the signatures of environmental changes in the patterns of this variation have been important areas of focus in population genetics. On one hand, periods of constant selection pull the genetic makeup of populations toward local fitness optima. On the other, to cope with changes in the selection regime, populations may evolve mechanisms that create a diversity of genotypes. By tuning the rates at which variability is produced--such as the rates of recombination, mutation, or migration--populations may increase their long-term adaptability. Here we use theoretical models to gain insight into how the rates of these three evolutionary forces are shaped by fluctuating selection. We compare and contrast the evolution of recombination, mutation, and migration under similar patterns of environmental change and show that these three sources of phenotypic variation are surprisingly similar in their response to changing selection. We show that the shape, size, variance, and asymmetry of environmental fluctuation have different but predictable effects on evolutionary dynamics.

A phylogenetic test of the Red Queen Hypothesis: outcrossing and parasitism in the Nematode phylum

Sexual outcrossing is costly relative to selfing and asexuality, yet it is ubiquitous in nature, a paradox that has long puzzled evolutionary biologists. The Red Queen Hypothesis argues that outcrossing is maintained by antagonistic interactions between host and parasites. Most tests of this hypothesis focus on the maintenance of outcrossing in hosts. The Red Queen makes an additional prediction that parasitic taxa are more likely to be outcrossing than their free-living relatives. We test this prediction in the diverse Nematode phylum using phylogenetic comparative methods to evaluate trait correlations. In support of the Red Queen, we demonstrate a significant correlation between parasitism and outcrossing in this clade. We find that this correlation is driven by animal parasites, for which outcrossing is significantly enriched relative to both free-living and plant parasitic taxa. Finally, we test hypotheses for the evolutionary history underlying the correlation of outcrossing and animal parasitism. Our results demonstrate that selfing and asexuality are significantly less likely to arise on parasitic lineages than on free-living ones. The findings of this study are consistent with the Red Queen Hypothesis. Moreover, they suggest that the maintenance of genetic variation is an important factor in the persistence of parasitic lineages.

Strong artificial selection in domestic mammals did not result in an increased recombination rate

Recombination rates vary in intensity and location at the species, individual, sex and chromosome levels. Despite the fundamental biological importance of this process, the selective forces that operate to shape recombination rate and patterns are unclear. Domestication offers a unique opportunity to study the interplay between recombination and selection. In domesticates, intense selection for particular traits is imposed on small populations over many generations, resulting in organisms that differ, sometimes dramatically, in morphology and physiology from their wild ancestor. Although earlier studies suggested increased recombination rate in domesticates, a formal comparison of recombination rates between domestic mammals and their wild congeners was missing. In order to determine broad-scale recombination rate, we used immunolabeling detection of MLH1 foci as crossover markers in spermatocytes in three pairs of closely related wild and domestic species (dog and wolf, goat and ibex, and sheep and mouflon). In the three pairs, and contrary to previous suggestions, our data show that contemporary recombination rate is higher in the wild species. Subsequently, we inferred recombination breakpoints in sequence data for 16 genomic regions in dogs and wolves, each containing a locus associated with a dog phenotype potentially under selection during domestication. No difference in the number and distribution of recombination breakpoints was found between dogs and wolves. We conclude that our data indicate that strong directional selection did not result in changes in recombination in domestic mammals, and that both upper and lower bounds for crossover rates may be tightly regulated.

Effects of functionally asexual reproduction on quantitative genetic variation in the evening primroses (Oenothera, Onagraceae)

PREMISE OF THE STUDY

It has long been predicted that a loss of sexual reproduction leads to decreased heritable variation within populations and increased differentiation between populations. Despite an abundance of theory, there are few empirical tests of how sex affects genetic variation in phenotypic traits, especially for plants. Here we test whether repeated losses of two critical components of sex (recombination and segregation) in the evening primroses (Oenothera L., Onagraceae) affect quantitative genetic variation within and between populations.

METHODS

We sampled multiple genetic families from 3-5 populations from each of eight Oenothera species, which represented four independent transitions between sexual reproduction and a functionally asexual genetic system called "permanent translocation heterozygosity." We used quantitative genetics methods to partition genetic variation within and between populations for eight plant traits related to growth, leaf physiology, flowering, and resistance to herbivores.

KEY RESULTS

Heritability was, on average, 74% higher in sexual Oenothera populations than in functionally asexual populations, with plant growth rate, specific leaf area, and the percentage of leaf water content showing the strongest differences. By contrast, genetic differentiation among populations was 2.8× higher in functionally asexual vs. sexual Oenothera species. This difference was particularly strong for specific leaf area. Sexual populations tended to exhibit higher genetic correlations among traits, but this difference was weakly supported.

CONCLUSIONS

These results support the prediction that sexual reproduction maintains higher genetic variation within populations, which may facilitate adaptive evolution. We also found partial support for the prediction that a loss of sex leads to greater population differentiation, which may elevate speciation rates.

Measuring quantitative virulence in the wheat pathogen Zymoseptoria tritici using high-throughput automated image analysis

Zymoseptoria tritici, causal agent of Septoria tritici blotch on wheat, produces pycnidia in chlorotic and necrotic lesions on infected leaves. A high-throughput phenotyping method was developed based on automated digital image analysis that accurately measures the percentage of leaf area covered by lesions (PLACL) as well as pycnidia size and number. A seedling inoculation assay was conducted using 361 Z. tritici isolates originating from a controlled cross and two different winter wheat cultivars. Pycnidia size and density were found to be quantitative traits that showed a continuous distribution in the progeny. There was a weak correlation between pycnidia density and size (r = -0.27) and between pycnidia density and PLACL (r = 0.37). There were significant differences in PLACL and pycnidia density on resistant and susceptible cultivars. In all, >20% of the offspring exhibited significantly different pycnidia density on the two cultivars, consistent with host specialization. Automated image analysis provided greater accuracy and precision compared with traditional visual estimates of virulence. These results show that digital image analysis provides a powerful tool for measuring differences in quantitative virulence among strains of Z. tritici.

Discriminating micropathogen lineages and their reticulate evolution through graph theory-based network analysis: the case of Trypanosoma cruzi, the agent of Chagas disease

Micropathogens (viruses, bacteria, fungi, parasitic protozoa) share a common trait, which is partial clonality, with wide variance in the respective influence of clonality and sexual recombination on the dynamics and evolution of taxa. The discrimination of distinct lineages and the reconstruction of their phylogenetic history are key information to infer their biomedical properties. However, the phylogenetic picture is often clouded by occasional events of recombination across divergent lineages, limiting the relevance of classical phylogenetic analysis and dichotomic trees. We have applied a network analysis based on graph theory to illustrate the relationships among genotypes of Trypanosoma cruzi, the parasitic protozoan responsible for Chagas disease, to identify major lineages and to unravel their past history of divergence and possible recombination events. At the scale of T. cruzi subspecific diversity, graph theory-based networks applied to 22 isoenzyme loci (262 distinct Multi-Locus-Enzyme-Electrophoresis -MLEE) and 19 microsatellite loci (66 Multi-Locus-Genotypes -MLG) fully confirms the high clustering of genotypes into major lineages or "near-clades". The release of the dichotomic constraint associated with phylogenetic reconstruction usually applied to Multilocus data allows identifying putative hybrids and their parental lineages. Reticulate topology suggests a slightly different history for some of the main "near-clades", and a possibly more complex origin for the putative hybrids than hitherto proposed. Finally the sub-network of the near-clade T. cruzi I (28 MLG) shows a clustering subdivision into three differentiated lesser near-clades ("Russian doll pattern"), which confirms the hypothesis recently proposed by other investigators. The present study broadens and clarifies the hypotheses previously obtained from classical markers on the same sets of data, which demonstrates the added value of this approach. This underlines the potential of graph theory-based network analysis for describing the nature and relationships of major pathogens, thereby opening stimulating prospects to unravel the organization, dynamics and history of major micropathogen lineages.

Recombination accelerates adaptation on a large-scale empirical fitness landscape in HIV-1

Recombination has the potential to facilitate adaptation. In spite of the substantial body of theory on the impact of recombination on the evolutionary dynamics of adapting populations, empirical evidence to test these theories is still scarce. We examined the effect of recombination on adaptation on a large-scale empirical fitness landscape in HIV-1 based on in vitro fitness measurements. Our results indicate that recombination substantially increases the rate of adaptation under a wide range of parameter values for population size, mutation rate and recombination rate. The accelerating effect of recombination is stronger for intermediate mutation rates but increases in a monotonic way with the recombination rates and population sizes that we examined. We also found that both fitness effects of individual mutations and epistatic fitness interactions cause recombination to accelerate adaptation. The estimated epistasis in the adapting populations is significantly negative. Our results highlight the importance of recombination in the evolution of HIV-I.

One of the most challenging issues in evolutionary biology concerns the question of why most organisms exchange genetic material with each other, e.g. during sexual reproduction. Gene shuffling can create genetic diversity that facilitates adaptation to new environments, but theory shows that this effect is highly dependent on how different genes interact in determining the fitness of an organism. Using a large data set of fitness values based on HIV-1, we provide evidence that shuffling of genetic material indeed raises the level of genetic diversity, and as a result accelerates adaptation. Our results also propose genetic shuffling as a mechanism utilized by HIV to accelerate the evolution of multi-drug-resistant strains.

Domestication and the mitochondrial genome: comparing patterns and rates of molecular evolution in domesticated mammals and birds and their wild relatives

Studies of domesticated animals have led to the suggestion that domestication could have significant effects on patterns of molecular evolution. In particular, analyses of mitochondrial genome sequences from domestic dogs and yaks have yielded higher ratios of non-synonymous to synonymous substitutions in the domesticated lineages than in their wild relatives. These results are important because they imply that changes to selection or population size operating over a short timescale can cause significant changes to the patterns of mitochondrial molecular evolution. In this study, our aim is to test whether the impact on mitochondrial genome evolution is a general feature of domestication or whether it is specific to particular examples. We test whether domesticated mammals and birds have consistently different patterns of molecular evolution than their wild relatives for 16 phylogenetically independent comparisons of mitochondrial genome sequences. We find no consistent difference in branch lengths or dN/dS between domesticated and wild lineages. We also find no evidence that our failure to detect a consistent pattern is due to the short timescales involved or low genetic distance between domesticated lineages and their wild relatives. However, removing comparisons where the wild relative may also have undergone a bottleneck does reveal a pattern consistent with reduced effective population size in domesticated lineages. Our results suggest that, although some domesticated lineages may have undergone changes to selective regime or effective population size that could have affected mitochondrial evolution, it is not possible to generalize these patterns over all domesticated mammals and birds.

Why epistasis is important for selection and adaptation

Organisms are built from thousands of genes that interact in complex ways. Still, the mathematical theory of evolution is dominated by a gene-by-gene perspective in which genes are assumed to have the same effects regardless of genetic background. Gene interaction, or epistasis, plays a role in some theoretical developments such as the evolution of recombination, reproductive isolation, and canalization, but is strikingly missing from our standard accounts of phenotypic adaptation. This absence is most puzzling within the field of quantitative genetics, which, despite its polygenic perspective and elaborate statistical representation of epistasis, has not found a single important role for gene interaction in evolution. To the contrary, there is a widespread consensus that epistasis is evolutionary inert, and that all we need to know to predict evolutionary dynamics is the additive component of the genetic variance. This view may have roots in convenience, but also in theoretical results showing that the response to selection derived from epistatic variance components is not permanent and will decay when selection is relaxed. I show that these results are tied to a conceptual confusion, and are misleading as general statements about the significance of epistasis for the selection response and adaptation.

Relative effects of segregation and recombination on the evolution of sex in finite diploid populations

The mechanism of reproducing more viable offspring in response to selection is a major factor influencing the advantages of sex. In diploids, sexual reproduction combines genotype by recombination and segregation. Theoretical studies of sexual reproduction have investigated the advantage of recombination in haploids. However, the potential advantage of segregation in diploids is less studied. This study aimed to quantify the relative contribution of recombination and segregation to the evolution of sex in finite diploids by using multilocus simulations. The mean fitness of a sexually or asexually reproduced population was calculated to describe the long-term effects of sex. The evolutionary fate of a sex or recombination modifier was also monitored to investigate the short-term effects of sex. Two different scenarios of mutations were considered: (1) only deleterious mutations were present and (2) a combination of deleterious and beneficial mutations. Results showed that the combined effects of segregation and recombination strongly contributed to the evolution of sex in diploids. If deleterious mutations were only present, segregation efficiently slowed down the speed of Muller's ratchet. As the recombination level was increased, the accumulation of deleterious mutations was totally inhibited and recombination substantially contributed to the evolution of sex. The presence of beneficial mutations evidently increased the fixation rate of a recombination modifier. We also observed that the twofold cost of sex was easily to overcome in diploids if a sex modifier caused a moderate frequency of sex.

Current hypotheses for the evolution of sex and recombination

The evolution of sex is one of the most important and controversial problems in evolutionary biology. Although sex is almost universal in higher animals and plants, its inherent costs have made its maintenance difficult to explain. The most famous of these is the twofold cost of males, which can greatly reduce the fecundity of a sexual population, compared to a population of asexual females. Over the past century, multiple hypotheses, along with experimental evidence to support these, have been put forward to explain widespread costly sex. In this review, we outline some of the most prominent theories, along with the experimental and observational evidence supporting these. Historically, there have been 4 classes of theories: the ability of sex to fix multiple novel advantageous mutants (Fisher-Muller hypothesis); sex as a mechanism to stop the build-up of deleterious mutations in finite populations (Muller's ratchet); recombination creating novel genotypes that can resist infection by parasites (Red Queen hypothesis); and the ability of sex to purge bad genomes if deleterious mutations act synergistically (mutational deterministic hypothesis). Current theoretical and experimental evidence seems to favor the hypothesis that sex breaks down selection interference between new mutants, or it acts as a mechanism to shuffle genotypes in order to repel parasitic invasion. However, there is still a need to collect more data from natural populations and experimental studies, which can be used to test different hypotheses.

Experimental measures of pathogen competition and relative fitness

Competition among pathogen strains for limited host resources can have a profound effect on pathogen evolution. A better understanding of the principles and consequences of competition can be useful in designing more sustainable disease management strategies. The competitive ability and relative fitness of a pathogen strain are determined by its intrinsic biological properties, the resistance and heterogeneity of the corresponding host population, the population density and genetic relatedness of the competing strains, and the physical environment. Competitive ability can be inferred indirectly from fitness components, such as basic reproduction rate or transmission rate. However, pathogen strains that exhibit higher fitness components when they infect a host alone may not exhibit a competitive advantage when they co-infect the same host. The most comprehensive measures of competitive ability and relative fitness come from calculating selection coefficients in a mixed infection in a field setting. Mark-release-recapture experiments can be used to estimate fitness costs associated with unnecessary virulence and fungicide resistance.

Three decades of studies on chromosomal polymorphism of Drosophila willistoni and description of fifty different rearrangements

Drosophila willistoni (Insecta, Diptera) is considered a paradigm for evolutionary studies. Their chromosomes are characterized by multiple paracentric inversions that make it hard to identify and describe chromosomal polymorphisms. In the present report we attempted to systematize the description of all the 50 inversions found in the last three decades, since we have been studying the chromosomes of several individuals of 30 different populations, including the one used in the genome sequencing project (Gd-H4-1). We present the photographic register of 11 arrangements in the left arm of the X chromosome (XL), eight in the right arm (XR), 10 in the left arm of chromosome II (IIL), eight in its right arm (IIR) and 13 in chromosome III. This information also includes their breakpoints on the reference photomap. A clear geographic difference was detected in XL and XR, with different fixed arrangements depending on the origin of the population studied. Through the comparison of all X arrangements it was possible to infer the putative ancestral arrangements, i.e., those related to all the remaining arrangements through the small number of inversions that occurred in the past, which we will call XL-A and XR-A. In the autosomes (IIL/IIR and III), fixed inversions were detected, but most are segregating in different frequencies along the geographical distribution of the D. willistoni populations.

Antibody selection and amino acid reversions

Pathogens adapt to antibody surveillance through amino acid replacements in targeted protein regions, or epitopes, that interfere with antibody binding. However, such escape mutations may exact a fitness cost due to impaired protein function. Here, it is hypothesized that the recurring generation of specific neutralizing antibodies to an epitope region as it evolves in response to antibody selection will cause amino acid reversions by releasing early escape mutations from immune selection. The plausibility of this hypothesis was tested with stochastic simulation of adaptation at the molecular sequence level in finite populations. Under the conditions of strong selection and weak mutation, the rates of allele fixation and amino acid reversion increased with population size and selection coefficients. These rates decreased with population size, however, if mutation became strong, because clonal interference reduced the rate of adaptation. The model successfully predicts the rate of reversion per allele fixation for an important human immunodeficiency virus type 1 (HIV-1) antibody epitope region. Therefore, antibody selection may generate complex adaptive dynamics.

Selection and inertia in the evolution of holocentric chromosomes in sedges (Carex, Cyperaceae)

• Changes in chromosome number as a result of fission and fusion in holocentrics have direct and immediate effects on the recombination rate. We investigate the support for the classic hypothesis that environmental stability selects for increased recombination rates. • We employed a phylogenetic and cytogenetic data set from one of the most diverse angiosperm genera in the world, which has the largest nonpolyploid chromosome radiation (Carex, Cyperaceae; 2n = 12-124; 2100 spp.). We evaluated alternative Ornstein-Uhlenbeck models of chromosome number adaptation to the environment in an information-theoretic framework. • We found moderate support for a positive influence of lateral inflorescence unit size on chromosome number, which may be selected in a stable environment in which resources for reproductive investment are larger. We found weak support for a positive influence on chromosome number of water-saturated soils and among-month temperature constancy, which would be expected to be negatively select for pioneering species. Chromosome number showed a strong phylogenetic signal. • We argue that our finding of small but significant effects of life history and ecology is compatible with our original hypothesis regarding selection of optima in recombination rates: low recombination rate is optimal when inmediate fitness is required. By contrast, high recombination rate is optimal when stable environments allow for evolutionary innovation.

Limits to the rate of adaptive substitution in sexual populations

In large populations, many beneficial mutations may be simultaneously available and may compete with one another, slowing adaptation. By finding the probability of fixation of a favorable allele in a simple model of a haploid sexual population, we find limits to the rate of adaptive substitution, Λ, that depend on simple parameter combinations. When variance in fitness is low and linkage is loose, the baseline rate of substitution is Λ₀ = 2NU , where N is the population size, U is the rate of beneficial mutations per genome, and is their mean selective advantage. Heritable variance v in log fitness due to unlinked loci reduces Λ by e⁻⁴(v) under polygamy and e⁻⁸ (v) under monogamy. With a linear genetic map of length R Morgans, interference is yet stronger. We use a scaling argument to show that the density of adaptive substitutions depends on s, N, U, and R only through the baseline density: Λ/R = F (Λ₀/R). Under the approximation that the interference due to different sweeps adds up, we show that Λ/R ~(Λ₀/R) / (1 +2Λ₉/R) , implying that interference prevents the rate of adaptive substitution from exceeding one per centimorgan per 200 generations. Simulations and numerical calculations confirm the scaling argument and confirm the additive approximation for Λ₀/R ~ 1; for higher Λ₀/R , the rate of adaptation grows above R/2, but only very slowly. We also consider the effect of sweeps on neutral diversity and show that, while even occasional sweeps can greatly reduce neutral diversity, this effect saturates as sweeps become more common-diversity can be maintained even in populations experiencing very strong interference. Our results indicate that for some organisms the rate of adaptive substitution may be primarily recombination-limited, depending only weakly on the mutation supply and the strength of selection.

Collapse

Effect of epistasis and linkage on fixation probability in three-locus models: an ancestral recombination-selection graph approach

We study the probability of ultimate fixation of a single new mutant arising in an individual chosen at random at a locus linked to two other loci carrying previously arisen mutations. This is done using the Ancestral Recombination-Selection Graph (ARSG) in a finite population in the limit of a large population size, which is also known as the Ancestral Influence Graph (AIG). An analytical expansion of the fixation probability with respect to population-scaled recombination rates and selection intensities is obtained. The coefficients of the expansion are expressed in terms of the initial state of the population and the epistatic interactions among the selected loci. Under the assumption of weak selection at tightly linked loci, the sign of the leading term, which depends on the signs of epistasis and initial linkage disequilibrium, determines whether an increase in recombination rates increases the chance of ultimate fixation of the new mutant. If mutants are advantageous, this is the case when epistasis is positive or null and the initial linkage disequilibrium is negative, which is an expected state in a finite population under directional selection. Moreover, this is also the case for a neutral mutant modifier coding for higher recombination rates if the same conditions hold at the selected loci. Under the same conditions, deleterious mutants are disfavored for ultimate fixation and neutral modifiers for higher recombination rates still favored. The recombination rates between the modifier locus and the selected loci do not come into play in the leading terms of the approximation for the fixation probability, but they do in higher-order terms.

The evolution of sex is favoured during adaptation to new environments

Both theory and experiments have demonstrated that sex can facilitate adaptation, potentially yielding a group-level advantage to sex. However, it is unclear whether this process can help solve the more difficult problem of the maintenance of sex within populations. Using experimental populations of the facultatively sexual rotifer Brachionus calyciflorus, we show that rates of sex evolve to higher levels during adaptation but then decline as fitness plateaus. To assess the fitness consequences of genetic mixing, we directly compare the fitnesses of sexually and asexually derived genotypes that naturally occur in our experimental populations. Sexually derived genotypes are more fit than asexually derived genotypes when adaptive pressures are strong, but this pattern reverses as the pace of adaptation slows, matching the pattern of evolutionary change in the rate of sex. These fitness assays test the net effect of sex but cannot be used to disentangle whether selection on sex arises because highly sexual lineages become associated with different allele combinations or with different allele frequencies than less sexual lineages (i.e., "short-" or "long-term" effects, respectively). We infer which of these mechanisms provides an advantage to sex by performing additional manipulations to obtain fitness distributions of sexual and asexual progeny arrays from unbiased parents (rather than from naturally occurring, and thereby evolutionarily biased, parents). We find evidence that sex breaks down adaptive gene combinations, resulting in lower average fitness of sexual progeny (i.e., a short-term disadvantage to sex). As predicted by theory, the advantage to sex arises because sexually derived progeny are more variable in fitness, allowing for faster adaptation. This "long-term advantage" builds over multiple generations, eventually resulting in higher fitness of sexual types.

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.

Antagonistic experimental coevolution with a parasite increases host recombination frequency

Background

One of the big remaining challenges in evolutionary biology is to understand the evolution and maintenance of meiotic recombination. As recombination breaks down successful genotypes, it should be selected for only under very limited conditions. Yet, recombination is very common and phylogenetically widespread. The Red Queen Hypothesis is one of the most prominent hypotheses for the adaptive value of recombination and sexual reproduction. The Red Queen Hypothesis predicts an advantage of recombination for hosts that are coevolving with their parasites. We tested predictions of the hypothesis with experimental coevolution using the red flour beetle, Tribolium castaneum, and its microsporidian parasite, Nosema whitei.

Results

By measuring recombination directly in the individuals under selection, we found that recombination in the host population was increased after 11 generations of coevolution. Detailed insights into genotypic and phenotypic changes occurring during the coevolution experiment furthermore helped us to reconstruct the coevolutionary dynamics that were associated with this increase in recombination frequency. As coevolved lines maintained higher genetic diversity than control lines, and because there was no evidence for heterozygote advantage or for a plastic response of recombination to infection, the observed increase in recombination most likely represented an adaptive host response under Red Queen dynamics.

Conclusions

This study provides direct, experimental evidence for an increase in recombination frequency under host-parasite coevolution in an obligatory outcrossing species. Combined with earlier results, the Red Queen process is the most likely explanation for this observation.

The red queen coupled with directional selection favours the evolution of sex

Why sexual reproduction has evolved to be such a widespread mode of reproduction remains a major question in evolutionary biology. Although previous studies have shown that increased sex and recombination can evolve in the presence of host-parasite interactions (the 'Red Queen hypothesis' for sex), many of these studies have assumed that multiple loci mediate infection vs. resistance. Data suggest, however, that a major locus is typically involved in antigen presentation and recognition. Here, we explore a model where only one locus mediates host-parasite interactions, but a second locus is subject to directional selection. Even though the effects of these genes on fitness are independent, we show that increased rates of sex and recombination are favoured at a modifier gene that alters the rate of genetic mixing. This result occurs because of selective interference in finite populations (the 'Hill-Robertson effect'), which also favours sex. These results suggest that the Red Queen hypothesis may help to explain the evolution of sex by contributing a form of persistent selection, which interferes with directional selection at other loci and thereby favours sex and recombination.

Direct and indirect consequences of meiotic recombination: implications for genome evolution

There is considerable variation within eukaryotic genomes in the local rate of crossing over. Why is this and what effect does it have on genome evolution? On the genome scale, it is known that by shuffling alleles, recombination increases the efficacy of selection. By contrast, the extent to which differences in the recombination rate modulate the efficacy of selection between genomic regions is unclear. Recombination also has direct consequences on the origin and fate of mutations: biased gene conversion and other forms of meiotic drive promote the fixation of mutations in a similar way to selection, and recombination itself may be mutagenic. Consideration of both the direct and indirect effects of recombination is necessary to understand why its rate is so variable and for correct interpretation of patterns of genome evolution.

Recombination rate variation in closely related species

Despite their importance to successful meiosis and various evolutionary processes, meiotic recombination rates sometimes vary within species or between closely related species. For example, humans and chimpanzees share virtually no recombination hotspot locations in the surveyed portion of the genomes. However, conservation of recombination rates between closely related species has also been documented, raising an apparent contradiction. Here, we evaluate how and why conflicting patterns of recombination rate conservation and divergence may be observed, with particular emphasis on features that affect recombination, and the scale and method with which recombination is surveyed. Additionally, we review recent studies identifying features influencing fine-scale and broad-scale recombination patterns and informing how quickly recombination rates evolve, how changes in recombination impact selection and evolution in natural populations, and more broadly, which forces influence genome evolution.

Fixation probability in a two-locus model by the ancestral recombination-selection graph

We use the ancestral influence graph (AIG) for a two-locus, two-allele selection model in the limit of a large population size to obtain an analytic approximation for the probability of ultimate fixation of a single mutant allele A. We assume that this new mutant is introduced at a given locus into a finite population in which a previous mutant allele B is already segregating with a wild type at another linked locus. We deduce that the fixation probability increases as the recombination rate increases if allele A is either in positive epistatic interaction with B and allele B is beneficial or in no epistatic interaction with B and then allele A itself is beneficial. This holds at least as long as the recombination fraction and the selection intensity are small enough and the population size is large enough. In particular this confirms the Hill-Robertson effect, which predicts that recombination renders more likely the ultimate fixation of beneficial mutants at different loci in a population in the presence of random genetic drift even in the absence of epistasis. More importantly, we show that this is true from weak negative epistasis to positive epistasis, at least under weak selection. In the case of deleterious mutants, the fixation probability decreases as the recombination rate increases. This supports Muller's ratchet mechanism to explain the accumulation of deleterious mutants in a population lacking recombination.

Differential selection between the sexes and selection for sex

Anisogamy is known to generate an important cost for sexual reproduction (the famous "twofold cost of sex"). However, male-female differences may have other consequences on the evolution of sex, due to the fact that selective pressures may differ among the sexes. On the one hand, intralocus sexual conflict should favor asexual females, which can fix female-beneficial, male-detrimental alleles. On the other hand, it has been suggested repeatedly that sexual selection among males may help to purge the mutation load, providing an advantage to sexual females. However, no analytical model has computed the strength of selection acting on a modifier gene affecting the frequency of sexual reproduction when selection differs between the sexes. In this article, we analyze a two-locus model using two approaches: a quasi-linkage-equilibrium (QLE) analysis and a local stability analysis, whose predictions are verified using a multilocus simulation. We find that costly sex can be maintained when selection is stronger in males than in females, but acts in the same direction in both. Complete asexuality, however, evolves under any other form of selection. Finally, we discuss how experimental measurements of fitness variances and covariances between sexes could be used to determine the overall direction and strength on selection for sex arising from differences in selection between males and females.

Variation in human recombination rates and its genetic determinants

BACKGROUND

Despite the fundamental role of crossing-over in the pairing and segregation of chromosomes during human meiosis, the rates and placements of events vary markedly among individuals. Characterizing this variation and identifying its determinants are essential steps in our understanding of the human recombination process and its evolution.

STUDY DESIGN/RESULTS

Using three large sets of European-American pedigrees, we examined variation in five recombination phenotypes that capture distinct aspects of crossing-over patterns. We found that the mean recombination rate in males and females and the historical hotspot usage are significantly heritable and are uncorrelated with one another. We then conducted a genome-wide association study in order to identify loci that influence them. We replicated associations of RNF212 with the mean rate in males and in females as well as the association of Inversion 17q21.31 with the female mean rate. We also replicated the association of PRDM9 with historical hotspot usage, finding that it explains most of the genetic variance in this phenotype. In addition, we identified a set of new candidate regions for further validation.

SIGNIFICANCE

These findings suggest that variation at broad and fine scales is largely separable and that, beyond three known loci, there is no evidence for common variation with large effects on recombination phenotypes.

Evolution of the genomic recombination rate in murid rodents

Although very closely related species can differ in their fine-scale patterns of recombination hotspots, variation in the average genomic recombination rate among recently diverged taxa has rarely been surveyed. We measured recombination rates in eight species that collectively represent several temporal scales of divergence within a single rodent family, Muridae. We used a cytological approach that enables in situ visualization of crossovers at meiosis to quantify recombination rates in multiple males from each rodent group. We uncovered large differences in genomic recombination rate between rodent species, which were independent of karyotypic variation. The divergence in genomic recombination rate that we document is not proportional to DNA sequence divergence, suggesting that recombination has evolved at variable rates along the murid phylogeny. Additionally, we document significant variation in genomic recombination rate both within and between subspecies of house mice. Recombination rates estimated in F(1) hybrids reveal evidence for sex-linked loci contributing to the evolution of recombination in house mice. Our results provide one of the first detailed portraits of genomic-scale recombination rate variation within a single mammalian family and demonstrate that the low recombination rates in laboratory mice and rats reflect a more general reduction in recombination rate across murid rodents.

Genome structure and the benefit of sex

We examine the behavior of sexual and asexual populations in modular multipeaked fitness landscapes and show that sexuals can systematically reach different, higher fitness adaptive peaks than asexuals. Whereas asexuals must move against selection to escape local optima, sexuals reach higher fitness peaks reliably because they create specific genetic variants that "skip over" fitness valleys, moving from peak to peak in the fitness landscape. This occurs because recombination can supply combinations of mutations in functional composites or "modules," that may include individually deleterious mutations. Thus when a beneficial module is substituted for another less-fit module by sexual recombination it provides a genetic variant that would require either several specific simultaneous mutations in an asexual population or a sequence of individual mutations some of which would be selected against. This effect requires modular genomes, such that subsets of strongly epistatic mutations are tightly physically linked. We argue that such a structure is provided simply by virtue of the fact that genomes contain many genes each containing many strongly epistatic nucleotides. We briefly discuss the connections with "building blocks" in the evolutionary computation literature. We conclude that there are conditions in which sexuals can systematically evolve high-fitness genotypes that are essentially unevolvable for asexuals.

Abortive meiosis in the oogenesis of parthenogenetic Daphnia pulex

Most daphnid species adopt parthenogenesis and sexual reproduction differentially in response to varied environmental cues, resulting in the production of diploid progenies in both cases. Previous studies have reportedly suggested that daphnids produce their parthenogenetic eggs via apomixis; the nuclear division of mature oocytes should be an equational division similar to somatic mitosis. However, it seems premature to conclude that this has been unequivocally established in any daphnids. Therefore, the objective of our research was to precisely reveal the process and mechanism of parthenogenetic oogenesis and maintenance of diploidy in Daphnia pulex through histology, karyology, and immunohistochemistry. We found that, when a parthenogenetic egg entered the first meiosis, division was arrested in the early first anaphase. Then, two half-bivalents, which were dismembered from each bivalent, moved back to the equatorial plate and assembled to form a diploid equatorial plate. Finally, the sister chromatids were separated and moved to opposite poles in the same manner as the second meiotic division followed by the extrusion of one extremely small daughter cell (resembling a polar body). These results suggest that parthenogenetic D. pulex do not adopt typical apomixis. We hypothesize that D. pulex switches reproductive mode depending on whether the egg is fertilized or not.