Resolving the paradox of sex and recombination

Linkage Mapping Reveals Strong Chiasma Interference in Sockeye Salmon: Implications for Interpreting Genomic Data

Meiotic recombination is fundamental for generating new genetic variation and for securing proper disjunction. Further, recombination plays an essential role during the rediploidization process of polyploid-origin genomes because crossovers between pairs of homeologous chromosomes retain duplicated regions. A better understanding of how recombination affects genome evolution is crucial for interpreting genomic data; unfortunately, current knowledge mainly originates from a few model species. Salmonid fishes provide a valuable system for studying the effects of recombination in nonmodel species. Salmonid females generally produce thousands of embryos, providing large families for conducting inheritance studies. Further, salmonid genomes are currently rediploidizing after a whole genome duplication and can serve as models for studying the role of homeologous crossovers on genome evolution. Here, we present a detailed interrogation of recombination patterns in sockeye salmon (Oncorhynchus nerka). First, we use RAD sequencing of haploid and diploid gynogenetic families to construct a dense linkage map that includes paralogous loci and location of centromeres. We find a nonrandom distribution of paralogs that mainly cluster in extended regions distally located on 11 different chromosomes, consistent with ongoing homeologous recombination in these regions. We also estimate the strength of interference across each chromosome; results reveal strong interference and crossovers are mostly limited to one per arm. Interference was further shown to continue across centromeres, but metacentric chromosomes generally had at least one crossover on each arm. We discuss the relevance of these findings for both mapping and population genomic studies.

Contribution of recombination to the evolutionary history of HIV

PURPOSE OF REVIEW

An improved understanding of how recombination affects the evolutionary history of HIV is crucial to understand its current and future evolution. The present review aims to disentangle the manifold effects of recombination on HIV by discussing its effects on the evolutionary history and the adaptive potential of HIV in the context of concepts from evolutionary genetics and genomics.

RECENT FINDINGS

The increasing occurrence of secondary contacts between divergent subtype populations (during coinfection) results in increased observations of recombinants worldwide. Recombination is heterogeneous along the HIV genome. Consequences of recombination of HIV evolution are, in combination with other demographic processes, expected to either homogenize the genetic composition of HIV populations (homogenization) or provide the potential for novel adaptations (diversification). New methods in population genomics allow deep characterization of recombinant genome (the segment composition and origin) and their evolutionary trajectories.

SUMMARY

HIV recombinants increase worldwide and invade geographical regions where pure subtypes were previously predominant. This trend is expected to continue in the future, as ease to travel worldwide increases opportunities for recombination between divergent HIV strains. While the effects of recombination in HIV are much researched, more effort is required to characterize current HIV recombinant composition and dynamics. This can be achieved with new population genetic and genomic methods.

The status of supergenes in the 21st century: recombination suppression in Batesian mimicry and sex chromosomes and other complex adaptations

I review theoretical models for the evolution of supergenes in the cases of Batesian mimicry in butterflies, distylous plants and sex chromosomes. For each of these systems, I outline the genetic evidence that led to the proposal that they involve multiple genes that interact during ‘complex adaptations’, and at which the mutations involved are not unconditionally advantageous, but show advantages that trade‐off against some disadvantages. I describe recent molecular genetic studies of these systems and questions they raise about the evolution of suppressed recombination. Nonrecombining regions of sex chromosomes have long been known, but it is not yet fully understood why recombination suppression repeatedly evolved in systems in distantly related taxa, but does not always evolve. Recent studies of distylous plants are tending to support the existence of recombination‐suppressed genome regions, which may include modest numbers of genes and resemble recently evolved sex‐linked regions. For Batesian mimicry, however, molecular genetic work in two butterfly species suggests a new supergene scenario, with a single gene mutating to produce initial adaptive phenotypes, perhaps followed by modifiers specifically refining and perfecting the new phenotype.

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.

From molecules to mating: Rapid evolution and biochemical studies of reproductive proteins

Sexual reproduction and the exchange of genetic information are essential biological processes for species across all branches of the tree of life. Over the last four decades, biochemists have continued to identify many of the factors that facilitate reproduction, but the molecular mechanisms that mediate this process continue to elude us. However, a recurring observation in this research has been the rapid evolution of reproductive proteins. In animals, the competing interests of males and females often result in arms race dynamics between pairs of interacting proteins. This phenomenon has been observed in all stages of reproduction, including pheromones, seminal fluid components, and gamete recognition proteins. In this article, we review how the integration of evolutionary theory with biochemical experiments can be used to study interacting reproductive proteins. Examples are included from both model and non-model organisms, and recent studies are highlighted for their use of state-of-the-art genomic and proteomic techniques.

SIGNIFICANCE

Despite decades of research, our understanding of the molecular mechanisms that mediate fertilization remain poorly characterized. To date, molecular evolutionary studies on both model and non-model organisms have provided some of the best inferences to elucidating the molecular underpinnings of animal reproduction. This review article details how biochemical and evolutionary experiments have jointly enhanced the field for 40 years, and how recent work using high-throughput genomic and proteomic techniques have shed additional insights into this crucial biological process.

Selective Advantages of a Parasexual Cycle for the Yeast Candida albicans

The yeast Candida albicans can mate. However, in the natural environment mating may generate progeny (fusants) fitter than clonal lineages too rarely to render mating biologically significant: C. albicans has never been observed to mate in its natural environment, the human host, and the population structure of the species is largely clonal. It seems incapable of meiosis, and most isolates are diploid and carry both mating-type-like (MTL) locus alleles, preventing mating. Only chromosome loss or localized loss of heterozygosity can generate mating-competent cells, and recombination of parental alleles is limited. To determine if mating is a biologically significant process, we investigated if mating is under selection. The ratio of nonsynonymous to synonymous mutations in mating genes and the frequency of mutations abolishing mating indicated that mating is under selection. The MTL locus is located on chromosome 5, and when we induced chromosome 5 loss in 10 clinical isolates, most of the resulting MTL-homozygotes could mate with each other, producing fusants. In laboratory culture, a novel environment favoring novel genotypes, some fusants grew faster than their parents, in which loss of heterozygosity had reduced growth rates, and also faster than their MTL-heterozygous ancestors-albeit often only after serial propagation. In a small number of experiments in which co-inoculation of an oral colonization model with MTL-homozygotes yielded small numbers of fusants, their numbers declined over time relative to those of the parents. Overall, our results indicate that mating generates genotypes superior to existing MTL-heterozygotes often enough to be under selection.

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.

Genetics of dioecy and causal sex chromosomes in plants

Dioecy (separate male and female individuals) ensures outcrossing and is more prevalent in animals than in plants. Although it is common in bryophytes and gymnosperms, only 5% of angiosperms are dioecious. In dioecious higher plants, flowers borne on male and female individuals are, respectively deficient in functional gynoecium and androecium. Dioecy is inherited via three sex chromosome systems: XX/XY, XX/X0 and WZ/ZZ, such that XX or WZ is female and XY, X0 or ZZ are males. The XX/XY system generates the rarer XX/X0 and WZ/ZZ systems. An autosome pair begets XY chromosomes. A recessive loss-of-androecium mutation (ana) creates X chromosome and a dominant gynoecium-suppressing (GYS) mutation creates Y chromosome. The ana/ANA and gys/GYS loci are in the sex-determining region (SDR) of the XY pair. Accumulation of inversions, deleterious mutations and repeat elements, especially transposons, in the SDR of Y suppresses recombination between X and Y in SDR, making Y labile and increasingly degenerate and heteromorphic from X. Continued recombination between X and Y in their pseudoautosomal region located at the ends of chromosomal arms allows survival of the degenerated Y and of the species. Dioecy is presumably a component of the evolutionary cycle for the origin of new species. Inbred hermaphrodite species assume dioecy. Later they suffer degenerate-Y-led population regression. Cross-hybridization between such extinguishing species and heterologous species, followed by genome duplication of segregants from hybrids, give rise to new species.

Female sterility associated with increased clonal propagation suggests a unique combination of androdioecy and asexual reproduction in populations of Cardamine amara (Brassicaceae)

BACKGROUND AND AIMS

The coexistence of hermaphrodites and female-sterile individuals, or androdioecy, has been documented in only a handful of plants and animals. This study reports its existence in the plant species Cardamine amara (Brassicaceae), in which female-sterile individuals have shorter pistils than seed-producing hermaphrodites.

METHODS

Morphological analysis, in situ manual pollination, microsatellite genotyping and differential gene expression analysis using Arabidopsis microarrays were used to delimit variation between female-sterile individuals and hermaphrodites.

KEY RESULTS

Female sterility in C. amara appears to be caused by disrupted ovule development. It was associated with a 2.4- to 2.9-fold increase in clonal propagation. This made the pollen number of female-sterile genets more than double that of hermaphrodite genets, which fulfils a condition of co-existence predicted by simple androdioecy theories. When female-sterile individuals were observed in wild androdioecious populations, their ramet frequencies ranged from 5 to 54 %; however, their genet frequencies ranged from 11 to 29 %, which is consistent with the theoretically predicted upper limit of 50 %.

CONCLUSIONS

The results suggest that a combination of sexual reproduction and increased asexual proliferation by female-sterile individuals probably explains the invasion and maintenance of female sterility in otherwise hermaphroditic populations. To our knowledge, this is the first report of the coexistence of female sterility and hermaphrodites in the Brassicaceae.

Fluorescent in situ hybridization shows DIPLOSPOROUS located on one of the NOR chromosomes in apomictic dandelions (Taraxacum) in the absence of a large hemizygous chromosomal region

Apomixis in dandelions (Taraxacum: Asteraceae) is encoded by two unlinked dominant loci and a third yet undefined genetic factor: diplosporous omission of meiosis (DIPLOSPOROUS, DIP), parthenogenetic embryo development (PARTHENOGENESIS, PAR), and autonomous endosperm formation, respectively. In this study, we determined the chromosomal position of the DIP locus in Taraxacum by using fluorescent in situ hybridization (FISH) with bacterial artificial chromosomes (BACs) that genetically map within 1.2-0.2 cM of DIP. The BACs showed dispersed fluorescent signals, except for S4-BAC 83 that displayed strong unique signals as well. Under stringent blocking of repeats by C0t-DNA fragments, only a few fluorescent foci restricted to defined chromosome regions remained, including one on the nucleolus organizer region (NOR) chromosomes that contains the 45S rDNAs. FISH with S4-BAC 83 alone and optimal blocking showed discrete foci in the middle of the long arm of one of the NOR chromosomes only in triploid and tetraploid diplosporous dandelions, while signals in sexual diploids were lacking. This agrees with the genetic model of a single dose, dominant DIP allele, absent in sexuals. The length of the DIP region is estimated to cover a region of 1-10 Mb. FISH in various accessions of Taraxacum and the apomictic sister species Chondrilla juncea, confirmed the chromosomal position of DIP within Taraxacum but not outside the genus. Our results endorse that, compared to other model apomictic species, expressing either diplospory or apospory, the genome of Taraxacum shows a more similar and less diverged chromosome structure at the DIP locus. The different levels of allele sequence divergence at apomeiosis loci may reflect different terms of asexual reproduction. The association of apomeiosis loci with repetitiveness, dispersed repeats, and retrotransposons commonly observed in apomictic species may imply a functional role of these shared features in apomictic reproduction, as is discussed.

Contrasting reproductive strategies of triploid hybrid males in vertebrate mating systems

The scarcity of parthenogenetic vertebrates is often attributed to their 'inferior' mode of clonal reproduction, which restricts them to self-reproduce their own genotype lineage and leaves little evolutionary potential with regard to speciation and evolution of sexual reproduction. Here, we show that for some taxa, such uniformity does not hold. Using hybridogenetic water frogs (Pelophylax esculentus) as a model system, we demonstrate that triploid hybrid males from two geographic regions exhibit very different reproductive modes. With an integrative data set combining field studies, crossing experiments, flow cytometry and microsatellite analyses, we found that triploid hybrids from Central Europe are rare, occur in male sex only and form diploid gametes of a single clonal lineage. In contrast, triploid hybrids from north-western Europe are widespread, occur in both sexes and produce recombined haploid gametes. These differences translate into contrasting reproductive roles between regions. In Central Europe, triploid hybrid males sexually parasitize diploid hybrids and just perpetuate their own genotype--which is the usual pattern in parthenogens. In north-western Europe, on the other hand, the triploid males are gamete donors for diploid hybrids, thereby stabilizing the mixed 2n-3n hybrid populations. By demonstrating these contrasting roles in male reproduction, we draw attention to a new significant evolutionary potential for animals with nonsexual reproduction, namely reproductive plasticity.

Recombination, chromosome number and eusociality in the Hymenoptera

Extraordinarily high rates of recombination have been observed in some eusocial species. The most popular explanation is that increased recombination increases genetic variation among workers, which in turn increases colony performance, for example by increasing parasite resistance. However, support for the generality of higher recombination rates among eusocial organisms remains weak, due to low sample size and a lack of phylogenetic independence of observations. Recombination rate, although difficult to measure directly, is correlated with chromosome number. As predicted, several authors have noted that chromosome numbers are higher among the eusocial species of Hymenoptera (ants, bees and wasps). Here, we present a formal comparative analysis of karyotype data from 1567 species of Hymenoptera. Contrary to earlier studies, we find no evidence for an absolute difference between chromosome number in eusocial and solitary species of Hymenoptera. However, we find support for an increased rate of chromosome number change in eusocial taxa. We show that among eusocial taxa colony size is able to explain some of the variation in chromosome number: intermediate-sized colonies have more chromosomes than those that are either very small or very large. However, we were unable to detect effects of a number of other colony characteristics predicted to affect recombination rate - including colony relatedness and caste number. Taken together, our results support the view that a eusocial lifestyle has led to variable selection pressure for increased recombination rates, but that identifying the factors contributing to this variable selection will require further theoretical and empirical effort.

Unisexual reproduction

Sexual reproduction is ubiquitous throughout the eukaryotic kingdom, but the capacity of pathogenic fungi to undergo sexual reproduction has been a matter of intense debate. Pathogenic fungi maintained a complement of conserved meiotic genes but the populations appeared to be clonally derived. This debate was resolved first with the discovery of an extant sexual cycle and then unisexual reproduction. Unisexual reproduction is a distinct form of homothallism that dispenses with the requirement for an opposite mating type. Pathogenic and nonpathogenic fungi previously thought to be asexual are able to undergo robust unisexual reproduction. We review here recent advances in our understanding of the genetic and molecular basis of unisexual reproduction throughout fungi and the impact of unisex on the ecology and genomic evolution of fungal species.

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.

Termite queens close the sperm gates of eggs to switch from sexual to asexual reproduction

Males and females are in conflict over genetic transmission in the evolution of parthenogenesis, because it enhances female reproductive output but deprives the males' genetic contribution. For males, any trait that coerces females into sexual reproduction should increase their fitness. However, in the termite Reticulitermes speratus, queens produce their replacements (neotenic queens) parthenogenetically while using normal sexual reproduction to produce other colony members. Here, we show that termite queens produce parthenogenetic offspring in the presence of kings by closing the micropyles (sperm gates; i.e., openings for sperm entry) of their eggs. Our field survey showed that termite eggs show large variation in numbers of micropyles, with some having none. Microsatellite analysis showed that embryos of micropyleless eggs develop parthenogenetically, whereas those of eggs with micropyles are fertilized and develop sexually. Surveys of eggs among queens of different age groups showed that queens begin to lay micropyleless eggs when they are older and thus, need to produce their replacements parthenogenetically. In addition, we found clear seasonality in new neotenic queen differentiation and micropyleless egg production. This micropyle-dependent parthenogenesis is the first identification, to our knowledge, of the mechanism through which females control egg fertilization over time in diploid animals, implying a novel route of the evolution of parthenogenesis in favor of female interests without interference from males.

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.

White cells facilitate opposite- and same-sex mating of opaque cells in Candida albicans

Modes of sexual reproduction in eukaryotic organisms are extremely diverse. The human fungal pathogen Candida albicans undergoes a phenotypic switch from the white to the opaque phase in order to become mating-competent. In this study, we report that functionally- and morphologically-differentiated white and opaque cells show a coordinated behavior during mating. Although white cells are mating-incompetent, they can produce sexual pheromones when treated with pheromones of the opposite mating type or by physically interacting with opaque cells of the opposite mating type. In a co-culture system, pheromones released by white cells induce opaque cells to form mating projections, and facilitate both opposite- and same-sex mating of opaque cells. Deletion of genes encoding the pheromone precursor proteins and inactivation of the pheromone response signaling pathway (Ste2-MAPK-Cph1) impair the promoting role of white cells (MTLa) in the sexual mating of opaque cells. White and opaque cells communicate via a paracrine pheromone signaling system, creating an environment conducive to sexual mating. This coordination between the two different cell types may be a trade-off strategy between sexual and asexual lifestyles in C. albicans.

In eukaryotic organisms, cells often undergo differentiation into distinct cell types in order to fulfill specialized roles. To achieve a certain function, different cell types may behave coordinately to complete a task that they may otherwise be incapable of completing independently. The human fungal pathogen Candida albicans can exist as two functionally and morphologically distinct cell types: white and opaque. The white cell type is thought to be the default state and may be the majority cell population in nature. However, only the minority opaque cells are mating-competent. In this study, we report that white and opaque cells show a coordinated behavior in the process of mating. When in the presence of opaque cells with an opposite mating type, white cells release sexual pheromones, and thus create an environment conducive for both opposite- and same-sex mating of opaque cells. The two cell types communicate via a paracrine pheromone signaling system. We propose that this communal coordination between white and opaque cells may not only support the fungus to be a successful commensal and pathogen in the host, but may also increase the fitness of the fungus during evolution over time.

Effective population size does not predict codon usage bias in mammals

Synonymous codons are not used at equal frequency throughout the genome, a phenomenon termed codon usage bias (CUB). It is often assumed that interspecific variation in the intensity of CUB is related to species differences in effective population sizes (N_e), with selection on CUB operating less efficiently in species with small N_e. Here, we specifically ask whether variation in N_e predicts differences in CUB in mammals and report two main findings. First, across 41 mammalian genomes, CUB was not correlated with two indirect proxies of N_e (body mass and generation time), even though there was statistically significant evidence of selection shaping CUB across all species. Interestingly, autosomal genes showed higher codon usage bias compared to X-linked genes, and high-recombination genes showed higher codon usage bias compared to low recombination genes, suggesting intraspecific variation in N_e predicts variation in CUB. Second, across six mammalian species with genetic estimates of N_e (human, chimpanzee, rabbit, and three mouse species: Mus musculus, M. domesticus, and M. castaneus), N_e and CUB were weakly and inconsistently correlated. At least in mammals, interspecific divergence in N_e does not strongly predict variation in CUB. One hypothesis is that each species responds to a unique distribution of selection coefficients, confounding any straightforward link between N_e and CUB.

Multidimensional epistasis and the transitory advantage of sex

Identifying and quantifying the benefits of sex and recombination is a long-standing problem in evolutionary theory. In particular, contradictory claims have been made about the existence of a benefit of recombination on high dimensional fitness landscapes in the presence of sign epistasis. Here we present a comparative numerical study of sexual and asexual evolutionary dynamics of haploids on tunably rugged model landscapes under strong selection, paying special attention to the temporal development of the evolutionary advantage of recombination and the link between population diversity and the rate of adaptation. We show that the adaptive advantage of recombination on static rugged landscapes is strictly transitory. At early times, an advantage of recombination arises through the possibility to combine individually occurring beneficial mutations, but this effect is reversed at longer times by the much more efficient trapping of recombining populations at local fitness peaks. These findings are explained by means of well-established results for a setup with only two loci. In accordance with the Red Queen hypothesis the transitory advantage can be prolonged indefinitely in fluctuating environments, and it is maximal when the environment fluctuates on the same time scale on which trapping at local optima typically occurs.

The relations between recombination rate and patterns of molecular variation and evolution in Drosophila

Genetic recombination affects levels of variability and the efficacy of selection because natural selection acting at one site affects evolutionary processes at linked sites. The variation in local recombination rates across the Drosophila genome provides excellent material for testing hypotheses concerning the evolutionary consequences of recombination. The current state of knowledge from studies of Drosophila genomics and population genetics is reviewed here. Selection at linked sites has influenced the relations between recombination rates and patterns of molecular variation and evolution, such that higher rates of recombination are associated with both higher levels of variability and a greater efficacy of selection. It seems likely that background selection against deleterious mutations is a major factor contributing to these patterns in genome regions in which crossing over is rare or absent, whereas selective sweeps of positively selected mutations probably play an important role in regions with crossing over.

Sex in the Mucoralean fungi

Sexual development is extant in virtually all eukaryotic species, including throughout the kingdom Fungi. Positioned within the opisthokonts along with metazoans, fungi serve as model systems to elucidate the genetics and impact of sexual development. Basal fungal lineages such as the Mucoralean fungi provide a unique basis to study sexual reproduction, in which common ancestral traits found in both animal and fungal lineages may be conserved. This review discusses the sexual development, sex loci, and evolution of the sex locus in the Mucoralean fungi, which sheds light on our understanding of the evolution and functions of sex.

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.

The effects of recombination on phenotypic exploration and robustness in evolution

Recombination is a commonly used genetic operator in artificial and computational evolutionary systems. It has been empirically shown to be essential for evolutionary processes. However, little has been done to analyze the effects of recombination on quantitative genotypic and phenotypic properties. The majority of studies only consider mutation, mainly due to the more serious consequences of recombination in reorganizing entire genomes. Here we adopt methods from evolutionary biology to analyze a simple, yet representative, genetic programming method, linear genetic programming. We demonstrate that recombination has less disruptive effects on phenotype than mutation, that it accelerates novel phenotypic exploration, and that it particularly promotes robust phenotypes and evolves genotypic robustness and synergistic epistasis. Our results corroborate an explanation for the prevalence of recombination in complex living organisms, and helps elucidate a better understanding of the evolutionary mechanisms involved in the design of complex artificial evolutionary systems and intelligent algorithms.

A high-density linkage map enables a second-generation collared flycatcher genome assembly and reveals the patterns of avian recombination rate variation and chromosomal evolution

Detailed linkage and recombination rate maps are necessary to use the full potential of genome sequencing and population genomic analyses. We used a custom collared flycatcher 50 K SNP array to develop a high-density linkage map with 37 262 markers assigned to 34 linkage groups in 33 autosomes and the Z chromosome. The best-order map contained 4215 markers, with a total distance of 3132 cm and a mean genetic distance between markers of 0.12 cm. Facilitated by the array being designed to include markers from most scaffolds, we obtained a second-generation assembly of the flycatcher genome that approaches full chromosome sequences (N50 super-scaffold size 20.2 Mb and with 1.042 Gb (of 1.116 Gb) anchored to and mostly ordered and oriented along chromosomes). We found that flycatcher and zebra finch chromosomes are entirely syntenic but that inversions at mean rates of 1.5–2.0 event (6.6–7.5 Mb) per My have changed the organization within chromosomes, rates high enough for inversions to potentially have been involved with many speciation events during avian evolution. The mean recombination rate was 3.1 cm/Mb and correlated closely with chromosome size, from 2 cm/Mb for chromosomes >100 Mb to >10 cm/Mb for chromosomes <10 Mb. This size dependence seemed entirely due to an obligate recombination event per chromosome; if 50 cm was subtracted from the genetic lengths of chromosomes, the rate per physical unit DNA was constant across chromosomes. Flycatcher recombination rate showed similar variation along chromosomes as chicken but lacked the large interior recombination deserts characteristic of zebra finch chromosomes.

Effects of new mutations on fitness: insights from models and data

The rates and properties of new mutations affecting fitness have implications for a number of outstanding questions in evolutionary biology. Obtaining estimates of mutation rates and effects has historically been challenging, and little theory has been available for predicting the distribution of fitness effects (DFE); however, there have been recent advances on both fronts. Extreme-value theory predicts the DFE of beneficial mutations in well-adapted populations, while phenotypic fitness landscape models make predictions for the DFE of all mutations as a function of the initial level of adaptation and the strength of stabilizing selection on traits underlying fitness. Direct experimental evidence confirms predictions on the DFE of beneficial mutations and favors distributions that are roughly exponential but bounded on the right. A growing number of studies infer the DFE using genomic patterns of polymorphism and divergence, recovering a wide range of DFE. Future work should be aimed at identifying factors driving the observed variation in the DFE. We emphasize the need for further theory explicitly incorporating the effects of partial pleiotropy and heterogeneity in the environment on the expected DFE.

Role of sex and migration in adaptation to sink environments

Understanding the effects of sex and migration on adaptation to novel environments remains a key problem in evolutionary biology. Using a single-cell alga Chlamydomonas reinhardtii, we investigated how sex and migration affected rates of evolutionary rescue in a sink environment, and subsequent changes in fitness following evolutionary rescue. We show that sex and migration affect both the rate of evolutionary rescue and subsequent adaptation. However, their combined effects change as the populations adapt to a sink habitat. Both sex and migration independently increased rates of evolutionary rescue, but the effect of sex on subsequent fitness improvements, following initial rescue, changed with migration, as sex was beneficial in the absence of migration but constraining adaptation when combined with migration. These results suggest that sex and migration are beneficial during the initial stages of adaptation, but can become detrimental as the population adapts to its environment.

High-resolution sex-specific linkage maps of the mouse reveal polarized distribution of crossovers in male germline

Since the publication of the first comprehensive linkage map for the laboratory mouse, the architecture of recombination as a basic biological process has become amenable to investigation in mammalian model organisms. Here we take advantage of high-density genotyping and the unique pedigree structure of the incipient Collaborative Cross to investigate the roles of sex and genetic background in mammalian recombination. Our results confirm the observation that map length is longer when measured through female meiosis than through male meiosis, but we find that this difference is modified by genotype at loci on both the X chromosome and the autosomes. In addition, we report a striking concentration of crossovers in the distal ends of autosomes in male meiosis that is absent in female meiosis. The presence of this pattern in both single- and double-recombinant chromosomes, combined with the absence of a corresponding asymmetry in the distribution of double-strand breaks, indicates a regulated sequence of events specific to male meiosis that is anchored by chromosome ends. This pattern is consistent with the timing of chromosome pairing and evolutionary constraints on male recombination. Finally, we identify large regions of reduced crossover frequency that together encompass 5% of the genome. Many of these "cold regions" are enriched for segmental duplications, suggesting an inverse local correlation between recombination rate and mutation rate for large copy number variants.

Comparative analysis of the gonadal transcriptomes of the all-female species Poecilia formosa and its maternal ancestor Poecilia mexicana

Background

The Amazon molly, Poecilia formosa (Teleostei: Poeciliinae) is an unisexual, all-female species. It evolved through the hybridisation of two closely related sexual species and exhibits clonal reproduction by sperm dependent parthenogenesis (or gynogenesis) where the sperm of a parental species is only used to activate embryogenesis of the apomictic, diploid eggs but does not contribute genetic material to the offspring.

Here we provide and describe the first de novo assembled transcriptome of the Amazon molly in comparison with its maternal ancestor, the Atlantic molly Poecilia mexicana. The transcriptome data were produced through sequencing of single end libraries (100 bp) with the Illumina sequencing technique.

Results

83,504,382 reads for the Amazon molly and 81,625,840 for the Atlantic molly were assembled into 127,283 and 78,961 contigs for the Amazon molly and the Atlantic molly, respectively. 63% resp. 57% of the contigs could be annotated with gene ontology terms after sequence similarity comparisons. Furthermore, we were able to identify genes normally involved in reproduction and especially in meiosis also in the transcriptome dataset of the apomictic reproducing Amazon molly.

Conclusions

We assembled and annotated the transcriptome of a non-model organism, the Amazon molly, without a reference genome (de novo). The obtained dataset is a fundamental resource for future research in functional and expression analysis. Also, the presence of 30 meiosis-specific genes within a species where no meiosis is known to take place is remarkable and raises new questions for future research.

Impact of the competition between mating types on the cultivation of Tuber melanosporum: Romeo and Juliet and the matter of space and time

Major breakthroughs in our understanding of the life cycles of the symbiotic ascomycetes belonging to the genus Tuber have occurred over the last several years. A number of Tuber species produce edible fruiting bodies, known as truffles, that are marketed worldwide. A better understanding of the basic biological characteristics of Tuber spp. is likely to have tremendous practical relevance for their cultivation. Tuber melanosporum produces the most valuable black truffles and its genome has been recently sequenced. This species is now serving as a model for studying the biology of truffles. Here, we review recent progress in the understanding of sexual reproduction modalities in T. melanosporum. The practical relevance of these findings is outlined. In particular, the discoveries that T. melanosporum is heterothallic and that strains of different mating types compete to persist on the roots of host plants suggest that the spatial and temporal distributional patterns of strains of different mating types are key determinants of truffle fructification. The spatial segregation of the two mating types in areas where T. melanosporum occurs likely limits truffle production. Thus, host plant inoculation techniques and agronomic practices that might be pursued to manage T. melanosporum orchards with a balanced presence of the two mating partners are described.

Replicated high-density genetic maps of two great tit populations reveal fine-scale genomic departures from sex-equal recombination rates

Linking variation in quantitative traits to variation in the genome is an important, but challenging task in the study of life-history evolution. Linkage maps provide a valuable tool for the unravelling of such trait-gene associations. Moreover, they give insight into recombination landscapes and between-species karyotype evolution. Here we used genotype data, generated from a 10k single-nucleotide polymorphism (SNP) chip, of over 2000 individuals to produce high-density linkage maps of the great tit (Parus major), a passerine bird that serves as a model species for ecological and evolutionary questions. We created independent maps from two distinct populations: a captive F2-cross from The Netherlands (NL) and a wild population from the United Kingdom (UK). The two maps contained 6554 SNPs in 32 linkage groups, spanning 2010 cM and 1917 cM for the NL and UK populations, respectively, and were similar in size and marker order. Subtle levels of heterochiasmy within and between chromosomes were remarkably consistent between the populations, suggesting that the local departures from sex-equal recombination rates have evolved. This key and surprising result would have been impossible to detect if only one population was mapped. A comparison with zebra finch Taeniopygia guttata, chicken Gallus gallus and the green anole lizard Anolis carolinensis genomes provided further insight into the evolution of avian karyotypes.

Horizontal gene transfer can rescue prokaryotes from Muller's ratchet: benefit of DNA from dead cells and population subdivision

Horizontal gene transfer (HGT) is a major factor in the evolution of prokaryotes. An intriguing question is whether HGT is maintained during evolution of prokaryotes owing to its adaptive value or is a byproduct of selection driven by other factors such as consumption of extracellular DNA (eDNA) as a nutrient. One hypothesis posits that HGT can restore genes inactivated by mutations and thereby prevent stochastic, irreversible deterioration of genomes in finite populations known as Muller’s ratchet. To examine this hypothesis, we developed a population genetic model of prokaryotes undergoing HGT via homologous recombination. Analysis of this model indicates that HGT can prevent the operation of Muller’s ratchet even when the source of transferred genes is eDNA that comes from dead cells and on average carries more deleterious mutations than the DNA of recipient live cells. Moreover, if HGT is sufficiently frequent and eDNA diffusion sufficiently rapid, a subdivided population is shown to be more resistant to Muller’s ratchet than an undivided population of an equal overall size. Thus, to maintain genomic information in the face of Muller’s ratchet, it is more advantageous to partition individuals into multiple subpopulations and let them “cross-reference” each other’s genetic information through HGT than to collect all individuals in one population and thereby maximize the efficacy of natural selection. Taken together, the results suggest that HGT could be an important condition for the long-term maintenance of genomic information in prokaryotes through the prevention of Muller’s ratchet.

Diversity of Prdm9 zinc finger array in wild mice unravels new facets of the evolutionary turnover of this coding minisatellite

In humans and mice, meiotic recombination events cluster into narrow hotspots whose genomic positions are defined by the PRDM9 protein via its DNA binding domain constituted of an array of zinc fingers (ZnFs). High polymorphism and rapid divergence of the Prdm9 gene ZnF domain appear to involve positive selection at DNA-recognition amino-acid positions, but the nature of the underlying evolutionary pressures remains a puzzle. Here we explore the variability of the Prdm9 ZnF array in wild mice, and uncovered a high allelic diversity of both ZnF copy number and identity with the caracterization of 113 alleles. We analyze features of the diversity of ZnF identity which is mostly due to non-synonymous changes at codons -1, 3 and 6 of each ZnF, corresponding to amino-acids involved in DNA binding. Using methods adapted to the minisatellite structure of the ZnF array, we infer a phylogenetic tree of these alleles. We find the sister species Mus spicilegus and M. macedonicus as well as the three house mouse (Mus musculus) subspecies to be polyphyletic. However some sublineages have expanded independently in Mus musculus musculus and M. m. domesticus, the latter further showing phylogeographic substructure. Compared to random genomic regions and non-coding minisatellites, none of these patterns appears exceptional. In silico prediction of DNA binding sites for each allele, overlap of their alignments to the genome and relative coverage of the different families of interspersed repeated elements suggest a large diversity between PRDM9 variants with a potential for highly divergent distributions of recombination events in the genome with little correlation to evolutionary distance. By compiling PRDM9 ZnF protein sequences in Primates, Muridae and Equids, we find different diversity patterns among the three amino-acids most critical for the DNA-recognition function, suggesting different diversification timescales.

Interactions between insect herbivores and plant mating systems

Self-pollination is common in plants, and limited seed and pollen dispersal can create localized inbreeding even within outcrossing plants. Consequently, insects regularly encounter inbred plants in nature. Because inbreeding results in elevated homozygosity, greater expression of recessive alleles, and subsequent phenotypic changes in inbred plants, inbreeding may alter plant-insect interactions. Recent research has found that plant inbreeding alters resistance and tolerance to herbivores, alters the attraction and susceptibility of plants to insects that vector plant pathogens, and alters visitation rates of insect pollinators. These results suggest that interactions with insects can increase or decrease inbreeding depression (the loss of fitness due to self-fertilization) and subsequently alter the evolution of selfing within plant populations. Future work needs to focus on the mechanisms underlying genetic variation in the effects of inbreeding on plant-insect interactions and the consequences of altered plant-insect interactions on the evolution of plant defense and plant mating systems.

The genomic landscape of meiotic crossovers and gene conversions in Arabidopsis thaliana

Knowledge of the exact distribution of meiotic crossovers (COs) and gene conversions (GCs) is essential for understanding many aspects of population genetics and evolution, from haplotype structure and long-distance genetic linkage to the generation of new allelic variants of genes. To this end, we resequenced the four products of 13 meiotic tetrads along with 10 doubled haploids derived from Arabidopsis thaliana hybrids. GC detection through short reads has previously been confounded by genomic rearrangements. Rigid filtering for misaligned reads allowed GC identification at high accuracy and revealed an ∼80-kb transposition, which undergoes copy-number changes mediated by meiotic recombination. Non-crossover associated GCs were extremely rare most likely due to their short average length of ∼25-50 bp, which is significantly shorter than the length of CO-associated GCs. Overall, recombination preferentially targeted non-methylated nucleosome-free regions at gene promoters, which showed significant enrichment of two sequence motifs. DOI: http://dx.doi.org/10.7554/eLife.01426.001.

Genetic interaction networks: better understand to better predict

A genetic interaction (GI) between two genes generally indicates that the phenotype of a double mutant differs from what is expected from each individual mutant. In the last decade, genome scale studies of quantitative GIs were completed using mainly synthetic genetic array technology and RNA interference in yeast and Caenorhabditis elegans. These studies raised questions regarding the functional interpretation of GIs, the relationship of genetic and molecular interaction networks, the usefulness of GI networks to infer gene function and co-functionality, the evolutionary conservation of GI, etc. While GIs have been used for decades to dissect signaling pathways in genetic models, their functional interpretations are still not trivial. The existence of a GI between two genes does not necessarily imply that these two genes code for interacting proteins or that the two genes are even expressed in the same cell. In fact, a GI only implies that the two genes share a functional relationship. These two genes may be involved in the same biological process or pathway; or they may also be involved in compensatory pathways with unrelated apparent function. Considering the powerful opportunity to better understand gene function, genetic relationship, robustness and evolution, provided by a genome-wide mapping of GIs, several in silico approaches have been employed to predict GIs in unicellular and multicellular organisms. Most of these methods used weighted data integration. In this article, we will review the later knowledge acquired on GI networks in metazoans by looking more closely into their relationship with pathways, biological processes and molecular complexes but also into their modularity and organization. We will also review the different in silico methods developed to predict GIs and will discuss how the knowledge acquired on GI networks can be used to design predictive tools with higher performances.

The evolution of resistance against good and bad infections

Opportunities for genetic exchange are abundant between bacteria and foreign genetic elements (FGEs) such as conjugative plasmids, transposable elements and bacteriophages. The genetic novelty that may arise from these forms of genetic exchange is potentially beneficial to bacterial hosts, but there are also potential costs, which may be considerable in the case of phage infection. Some bacterial resistance mechanisms target both beneficial and deleterious forms of genetic exchange. Using a general epidemiological model, we explored under which conditions such resistance mechanisms may evolve. We considered a population of hosts that may be infected by FGEs that either confer a benefit or are deleterious to host fitness, and we analysed the epidemiological and evolutionary outcomes of resistance evolving under different cost/benefit scenarios. We show that the degree of co-infection between these two types of infection is particularly important in determining the evolutionarily stable level of host resistance. We explore these results using the example of CRISPR-Cas, a form of bacterial immunity that targets a variety of FGEs, and we show the potential role of bacteriophage infection in selecting for resistance mechanisms that in turn limit the acquisition of plasmid-borne antibiotic resistance. Finally, beyond microbes, we discuss how endosymbiotic associations may have shaped the evolution of host immune responses to pathogens.

Should evolutionary geneticists worry about higher-order epistasis?

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

Mating type gene homologues and putative sex pheromone-sensing pathway in arbuscular mycorrhizal fungi, a presumably asexual plant root symbiont

The fungal kingdom displays a fascinating diversity of sex-determination systems. Recent advances in genomics provide insights into the molecular mechanisms of sex, mating type determination, and evolution of sexual reproduction in many fungal species in both ancient and modern phylogenetic lineages. All major fungal groups have evolved sexual differentiation and recombination pathways. However, sexuality is unknown in arbuscular mycorrhizal fungi (AMF) of the phylum Glomeromycota, an ecologically vital group of obligate plant root symbionts. AMF are commonly considered an ancient asexual lineage dating back to the Ordovician, approximately 460 M years ago. In this study, we used genomic and transcriptomic surveys of several AMF species to demonstrate the presence of conserved putative sex pheromone-sensing mitogen-activated protein (MAP) kinases, comparable to those described in Ascomycota and Basidiomycota. We also find genes for high mobility group (HMG) transcription factors, homologous to SexM and SexP genes in the Mucorales. The SexM genes show a remarkable sequence diversity among multiple copies in the genome, while only a single SexP sequence was detected in some isolates of Rhizophagus irregularis. In the Mucorales and Microsporidia, the sexM gene is flanked by genes for a triosephosphate transporter (TPT) and a RNA helicase, but we find no evidence for synteny in the vicinity of the Sex locus in AMF. Nonetheless, our results, together with previous observations on meiotic machinery, suggest that AMF could undergo a complete sexual reproduction cycle.

Evolutionary advantage via common action of recombination and neutrality

We investigate evolution models with recombination and neutrality. We consider the Crow-Kimura (parallel) mutation-selection model with the neutral fitness landscape, in which there is a central peak with high fitness A, and some of 1-point mutants have the same high fitness A, while the fitness of other sequences is 0. We find that the effect of recombination and neutrality depends on the concrete version of both neutrality and recombination. We consider three versions of neutrality: (a) all the nearest neighbor sequences of the peak sequence have the same high fitness A; (b) all the l-point mutations in a piece of genome of length l≥1 are neutral; (c) the neutral sequences are randomly distributed among the nearest neighbors of the peak sequences. We also consider three versions of recombination: (I) the simple horizontal gene transfer (HGT) of one nucleotide; (II) the exchange of a piece of genome of length l, HGT-l; (III) two-point crossover recombination (2CR). For the case of (a), the 2CR gives a rather strong contribution to the mean fitness, much stronger than that of HGT for a large genome length L. For the random distribution of neutral sequences there is a critical degree of neutrality ν(c), and for μ<μ(c) and (μ(c)-μ) is not large, the 2CR suppresses the mean fitness while HGT increases it; for ν much larger than ν(c), the 2CR and HGT-l increase the mean fitness larger than that of the HGT. We also consider the recombination in the case of smooth fitness landscapes. The recombination gives some advantage in the evolutionary dynamics, where recombination distinguishes clearly the mean-field-like evolutionary factors from the fluctuation-like ones. By contrast, mutations affect the mean-field-like and fluctuation-like factors similarly. Consequently, recombination can accelerate the non-mean-field (fluctuation) type dynamics without considerably affecting the mean-field-like factors.

The evolutionary role of males in C. elegans

Although the nematode worm Caenorhabditis elegans reproduces primarily as a self-fertilizing hermaphrodite, males are maintained in natural populations at low frequency. In this commentary, I discuss the evolutionary forces that maintain males and the role males might play in this mating system.

An alignment-free test for recombination

MOTIVATION

Why recombination? is one of the central questions in biology. This has led to a host of methods for quantifying recombination from sequence data. These methods are usually based on aligned DNA sequences. Here, we propose an efficient alignment-free alternative.

RESULTS

Our method is based on the distribution of match lengths, which we look up using enhanced suffix arrays. By eliminating the alignment step, the test becomes fast enough for application to whole bacterial genomes. Using simulations we show that our test has similar power as established tests when applied to long pairs of sequences. When applied to 58 genomes of Escherichia coli, we pick up the strongest recombination signal from a 125 kb horizontal gene transfer engineered 20 years ago.

AVAILABILITY AND IMPLEMENTATION

We have implemented our method in the command-line program rush. Its C sources and documentation are available under the GNU General Public License from http://guanine.evolbio.mpg.de/rush/.

Competing islands limit the rate of adaptation in structured populations

Beneficial mutations can co-occur when population structure slows down adaptation. Here, we consider the process of adaptation in asexual populations distributed over several locations ("islands"). New beneficial mutations arise at constant rate ub, and each mutation has the same selective advantage s>0. We assume that populations evolve within islands according to the successional mutations regime of Desai and Fisher (2007), that is, the time to local fixation of a mutation is short compared to the expected waiting time until the next mutation occurs. To study the rate of adaptation, we introduce an approximate model, the successional mutations (SM) model, which can be simulated efficiently and yields accurate results for a wide range of parameters. In the SM model, mutations fix instantly within islands, and migrants can take over the destination island if they are fitter than the residents. For the special case of a population distributed equally across two islands with population size N, we approximate the model further for small and large migration rates in comparison to the mutation rate. These approximations lead to explicit formulas for the rate of adaptation which fit the original model for a large range of parameter values. For the d island case we provide some heuristics on how to extend the explicit formulas and check these with computer simulations. We conclude that the SM model is a good approximation of the adaptation process in a structured population, at least if mutation or migration is limited.

Unisexual and heterosexual meiotic reproduction generate aneuploidy and phenotypic diversity de novo in the yeast Cryptococcus neoformans

Aneuploidy is known to be deleterious and underlies several common human diseases, including cancer and genetic disorders such as trisomy 21 in Down's syndrome. In contrast, aneuploidy can also be advantageous and in fungi confers antifungal drug resistance and enables rapid adaptive evolution. We report here that sexual reproduction generates phenotypic and genotypic diversity in the human pathogenic yeast Cryptococcus neoformans, which is globally distributed and commonly infects individuals with compromised immunity, such as HIV/AIDS patients, causing life-threatening meningoencephalitis. C. neoformans has a defined a-α opposite sexual cycle; however, >99% of isolates are of the α mating type. Interestingly, α cells can undergo α-α unisexual reproduction, even involving genotypically identical cells. A central question is: Why would cells mate with themselves given that sex is costly and typically serves to admix preexisting genetic diversity from genetically divergent parents? In this study, we demonstrate that α-α unisexual reproduction frequently generates phenotypic diversity, and the majority of these variant progeny are aneuploid. Aneuploidy is responsible for the observed phenotypic changes, as chromosome loss restoring euploidy results in a wild-type phenotype. Other genetic changes, including diploidization, chromosome length polymorphisms, SNPs, and indels, were also generated. Phenotypic/genotypic changes were not observed following asexual mitotic reproduction. Aneuploidy was also detected in progeny from a-α opposite-sex congenic mating; thus, both homothallic and heterothallic sexual reproduction can generate phenotypic diversity de novo. Our study suggests that the ability to undergo unisexual reproduction may be an evolutionary strategy for eukaryotic microbial pathogens, enabling de novo genotypic and phenotypic plasticity and facilitating rapid adaptation to novel environments.