Determinants of divergent adaptation and Dobzhansky-Muller interaction in experimental yeast populations

Reticulate evolution and rapid development of reproductive barriers upon secondary contact in a forest fungus

Reproductive barriers between sister species of the mushroom-forming fungi tend to be stronger in sympatry, leading to speculation on whether they are being reinforced by selection against hybrids. We have used population genomic analyses together with in vitro crosses of a global sample of the wood decay fungus Trichaptum abietinum to investigate reproductive barriers within this species complex and the processes that have shaped them. Our phylogeographic analyses show that T. abietinum is delimited into six major genetic groups: one in Asia, two in Europe, and three in North America. The groups present in Europe are interfertile and admixed, whereas our crosses show that the North American groups are reproductively isolated. In Asia, a more complex pattern appears, with partial intersterility between subgroups that likely originated independently and more recently than the reproductive barriers in North America. We found pre-mating barriers in T. abietinum to be moderately correlated with genomic divergence, whereas mean growth reduction of the mated hybrids showed a strong correlation with increasing genomic divergence. Genome-wide association analyses identified candidate genes with programmed cell death annotations, which are known to be involved in intersterility in distantly related fungi, although their link here remains unproven. Our demographic modeling and phylogenetic network analyses fit a scenario where reproductive barriers in Trichaptum abietinum could have been reinforced upon secondary contact between groups that diverged in allopatry during the Pleistocene glacial cycles. Our combination of experimental and genomic approaches demonstrates how T. abietinum is a tractable system for studying speciation mechanisms.

Escherichia coli adaptation under prolonged resource exhaustion is characterized by extreme parallelism and frequent historical contingency

Like many other non-sporulating bacterial species, Escherichia coli is able to survive prolonged periods of resource exhaustion, by entering a state of growth called long-term stationary phase (LTSP). In July 2015, we initiated a set of evolutionary experiments aimed at characterizing the dynamics of E. coli adaptation under LTSP. In these experiments populations of E. coli were allowed to initially grow on fresh rich media, but were not provided with any new external growth resources since their establishment. Utilizing whole genome sequencing data obtained for hundreds of clones sampled at 12 time points spanning the first six years of these experiments, we reveal several novel aspects of the dynamics of adaptation. First, we show that E. coli continuously adapts genetically, up to six years under resource exhaustion, through the highly convergent accumulation of mutations. We further show that upon entry into LTSP, long-lasting lineages are established. This lineage structure is in itself convergent, with similar lineages arising across independently evolving populations. The high parallelism with which adaptations occur under LTSP, combined with the LTSP populations' lineage structure, enable us to screen for pairs of loci displaying a significant association in the occurrence of mutations, suggestive of a historical contingency. We find that such associations are highly frequent and that a third of convergently mutated loci are involved in at least one such association. Combined our results demonstrate that LTSP adaptation is characterized by remarkably high parallelism and frequent historical contingency.

Rapid evolution of pre-zygotic reproductive barriers in allopatric populations

IMPORTANCE

A population diversifies into two or more species-such a process is known as speciation. In sexually reproducing microorganisms, which barriers arise first-pre-mating or post-mating? In this work, we quantify the relative strengths of these barriers and demonstrate that pre-mating barriers arise first in allopatrically evolving populations of yeast, Saccharomyces cerevisiae. These defects arise because of the altered kinetics of mating of the participating groups. Thus, our work provides an understanding of how adaptive changes can lead to diversification among microbial populations.

Enforcement of Postzygotic Species Boundaries in the Fungal Kingdom

Understanding the molecular basis of speciation is a primary goal in evolutionary biology. The formation of the postzygotic reproductive isolation that causes hybrid dysfunction, thereby reducing gene flow between diverging populations, is crucial for speciation. Using various advanced approaches, including chromosome replacement, hybrid introgression and transcriptomics, population genomics, and experimental evolution, scientists have revealed multiple mechanisms involved in postzygotic barriers in the fungal kingdom. These results illuminate both unique and general features of fungal speciation. Our review summarizes experiments on fungi exploring how Dobzhansky-Muller incompatibility, killer meiotic drive, chromosome rearrangements, and antirecombination contribute to postzygotic reproductive isolation. We also discuss possible evolutionary forces underlying different reproductive isolation mechanisms and the potential roles of the evolutionary arms race under the Red Queen hypothesis and epigenetic divergence in speciation.

Evolution and molecular bases of reproductive isolation

The most challenging problem in speciation research is disentangling the relative strength and order in which different reproductive barriers evolve. Here, we review recent developments in the study of reproductive isolation in yeasts. With over a thousand genome-sequenced isolates readily available for testing the viability, sterility, and fitness of both intraspecies and interspecies hybrid crosses, Saccharomyces yeasts are an ideal model to study such fundamental questions. Our survey demonstrates that, while chromosomal-level mutations are widespread at the intraspecific level, anti-recombination-driven chromosome missegregation is the primary reproductive barrier between species. Finally, despite their strength, all of these postzygotic barriers can be resolved through the asexual life history of hybrids.

Interpreting phylogenetic conflict: Hybridization in the most speciose genus of lichen-forming fungi

While advances in sequencing technologies have been invaluable for understanding evolutionary relationships, increasingly large genomic data sets may result in conflicting evolutionary signals that are often caused by biological processes, including hybridization. Hybridization has been detected in a variety of organisms, influencing evolutionary processes such as generating reproductive barriers and mixing standing genetic variation. Here, we investigate the potential role of hybridization in the diversification of the most speciose genus of lichen-forming fungi, Xanthoparmelia. As Xanthoparmelia is projected to have gone through recent, rapid diversification, this genus is particularly suitable for investigating and interpreting the origins of phylogenomic conflict. Focusing on a clade of Xanthoparmelia largely restricted to the Holarctic region, we used a genome skimming approach to generate 962 single-copy gene regions representing over 2 Mbp of the mycobiont genome. From this genome-scale dataset, we inferred evolutionary relationships using both concatenation and coalescent-based species tree approaches. We also used three independent tests for hybridization. Although different species tree reconstruction methods recovered largely consistent and well-supported trees, there was widespread incongruence among individual gene trees. Despite challenges in differentiating hybridization from ILS in situations of recent rapid radiations, our genome-wide analyses detected multiple potential hybridization events in the Holarctic clade, suggesting one possible source of trait variability in this hyperdiverse genus. This study highlights the value in using a pluralistic approach for characterizing genome-scale conflict, even in groups with well-resolved phylogenies, while highlighting current challenges in detecting the specific impacts of hybridization.

Interspecific hybridization as a driver of fungal evolution and adaptation

Cross-species gene transfer is often associated with bacteria, which have evolved several mechanisms that facilitate horizontal DNA exchange. However, the increased availability of whole-genome sequences has revealed that fungal species also exchange DNA, leading to intertwined lineages, blurred species boundaries or even novel species. In contrast to prokaryotes, fungal DNA exchange originates from interspecific hybridization, where two genomes are merged into a single, often highly unstable, polyploid genome that evolves rapidly into stabler derivatives. The resulting hybrids can display novel combinations of genetic and phenotypic variation that enhance fitness and allow colonization of new niches. Interspecific hybridization led to the emergence of important pathogens of humans and plants (for example, various Candida and 'powdery mildew' species, respectively) and industrially important yeasts, such as Saccharomyces hybrids that are important in the production of cold-fermented lagers or cold-cellared Belgian ales. In this Review, we discuss the genetic processes and evolutionary implications of fungal interspecific hybridization and highlight some of the best-studied examples. In addition, we explain how hybrids can be used to study molecular mechanisms underlying evolution, adaptation and speciation, and serve as a route towards development of new variants for industrial applications.

Defining and Disrupting Species Boundaries in Saccharomyces

The genus Saccharomyces is an evolutionary paradox. On the one hand, it is composed of at least eight clearly phylogenetically delineated species; these species are reproductively isolated from each other, and hybrids usually cannot complete their sexual life cycles. On the other hand, Saccharomyces species have a long evolutionary history of hybridization, which has phenotypic consequences for adaptation and domestication. A variety of cellular, ecological, and evolutionary mechanisms are responsible for this partial reproductive isolation among Saccharomyces species. These mechanisms have caused the evolution of diverse Saccharomyces species and hybrids, which occupy a variety of wild and domesticated habitats. In this article, we introduce readers to the mechanisms isolating Saccharomyces species, the circumstances in which reproductive isolation mechanisms are effective and ineffective, and the evolutionary consequences of partial reproductive isolation. We discuss both the evolutionary history of the genus Saccharomyces and the human history of taxonomists and biologists struggling with species concepts in this fascinating genus.

Environment-dependent costs and benefits of recombination in independently evolved populations of Escherichia coli

Understanding of the causes by which reproductive isolation arises remains limited. We examine the role of adaptation in driving reproductive isolation among 12 Escherichia coli populations evolved in two different environments. We found that, regardless of whether parents were selected in the same or different environments, the average fitness of recombinants was lower than the expected, consistent with a prevailing influence of incompatibility between independently accumulated mutations. Exceptions to this pattern occurred among recombinants of some parents evolved in different environments. These recombinants were less fit than expected in the selective environment of one parent, but more fit than expected in the selective environment of the other parent. Our results indicate that both parallel and divergent adaptation can quickly lead to intrinsic genetic barriers contributing to the initial stages of speciation and show that these barriers can be complex, for example, depending on the environment in which recombinant offspring are tested.

Mitochondria and the Origin of Species: Bridging Genetic and Ecological Perspectives on Speciation Processes

Mitochondria have been known to be involved in speciation through the generation of Dobzhansky-Muller incompatibilities, where functionally neutral co-evolution between mitochondrial and nuclear genomes can cause dysfunction when alleles are recombined in hybrids. We propose that adaptive mitochondrial divergence between populations can not only produce intrinsic (Dobzhansky-Muller) incompatibilities, but could also contribute to reproductive isolation through natural and sexual selection against migrants, post-mating prezygotic isolation, as well as by causing extrinsic reductions in hybrid fitness. We describe how these reproductive isolating barriers can potentially arise through adaptive divergence of mitochondrial function in the absence of mito-nuclear coevolution, a departure from more established views. While a role for mitochondria in the speciation process appears promising, we also highlight critical gaps of knowledge: (1) many systems with a potential for mitochondrially-mediated reproductive isolation lack crucial evidence directly linking reproductive isolation and mitochondrial function; (2) it often remains to be seen if mitochondrial barriers are a driver or a consequence of reproductive isolation; (3) the presence of substantial gene flow in the presence of mito-nuclear incompatibilities raises questions whether such incompatibilities are strong enough to drive speciation to completion; and (4) it remains to be tested how mitochondrial effects on reproductive isolation compare when multiple mechanisms of reproductive isolation coincide. We hope this perspective and the proposed research plans help to inform future studies of mitochondrial adaptation in a manner that links genotypic changes to phenotypic adaptations, fitness, and reproductive isolation in natural systems, helping to clarify the importance of mitochondria in the formation and maintenance of biological diversity.

Adaptation by Loss of Heterozygosity in Saccharomyces cerevisiae Clones Under Divergent Selection

Loss of heterozygosity (LOH) is observed during vegetative growth and reproduction of diploid genotypes through mitotic crossovers, aneuploidy caused by nondisjunction, and gene conversion. We aimed to test the role that LOH plays during adaptation of two highly heterozygous Saccharomyces cerevisiae genotypes to multiple environments over a short time span in the laboratory. We hypothesized that adaptation would be observed through parallel LOH events across replicate populations. Using genome resequencing of 70 clones, we found that LOH was widespread with 5.2 LOH events per clone after ∼500 generations. The most common mode of LOH was gene conversion (51%) followed by crossing over consistent with either break-induced replication or double Holliday junction resolution. There was no evidence that LOH involved nondisjunction of whole chromosomes. We observed parallel LOH in both an environment-specific and environment-independent manner. LOH largely involved recombining existing variation between the parental genotypes, but also was observed after de novo, presumably beneficial, mutations occurred in the presence of canavanine, a toxic analog of arginine. One highly parallel LOH event involved the ENA salt efflux pump locus on chromosome IV, which showed repeated LOH to the allele from the European parent, an allele originally derived by introgression from S. paradoxus Using CRISPR-engineered LOH we showed that the fitness advantage provided by this single LOH event was 27%. Overall, we found extensive evidence that LOH could be adaptive and is likely to be a greater source of initial variation than de novo mutation for rapid evolution of diploid genotypes.

The Past and Future of Experimental Speciation

Speciation is the result of evolutionary processes that generate barriers to gene flow between populations, facilitating reproductive isolation. Speciation is typically studied via theoretical models and snapshot tests in natural populations. Experimental speciation enables real-time direct tests of speciation theory and has been long touted as a critical complement to other approaches. We argue that, despite its promise to elucidate the evolution of reproductive isolation, experimental speciation has been underutilised and lags behind other contributions to speciation research. We review recent experiments and outline a framework for how experimental speciation can be implemented to address current outstanding questions that are otherwise challenging to answer. Greater uptake of this approach is necessary to rapidly advance understanding of speciation.

The evolution of hybrid fitness during speciation

The evolution of postzygotic reproductive isolation is an important component of speciation. But before isolation is complete there is sometimes a phase of heterosis in which hybrid fitness exceeds that of the two parental species. The genetics and evolution of heterosis and postzygotic isolation have typically been studied in isolation, precluding the development of a unified theory of speciation. Here, we develop a model that incorporates both positive and negative gene interactions, and accounts for the evolution of both heterosis and postzygotic isolation. We parameterize the model with recent data on the fitness effects of 10,000 mutations in yeast, singly and in pairwise epistatic combinations. The model makes novel predictions about the types of interactions that contribute to declining hybrid fitness. We reproduce patterns familiar from earlier models of speciation (e.g. Haldane's Rule and Darwin's Corollary) and identify new mechanisms that may underlie these patterns. Our approach provides a general framework for integrating experimental data from gene interaction networks into speciation theory and makes new predictions about the genetic mechanisms of speciation.

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

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

Contingency and determinism in evolution: Replaying life’s tape

Historical processes display some degree of “contingency,” meaning their outcomes are sensitive to seemingly inconsequential events that can fundamentally change the future. Contingency is what makes historical outcomes unpredictable. Unlike many other natural phenomena, evolution is a historical process. Evolutionary change is often driven by the deterministic force of natural selection, but natural selection works upon variation that arises unpredictably through time by random mutation, and even beneficial mutations can be lost by chance through genetic drift. Moreover, evolution has taken place within a planetary environment with a particular history of its own. This tension between determinism and contingency makes evolutionary biology a kind of hybrid between science and history. While philosophers of science examine the nuances of contingency, biologists have performed many empirical studies of evolutionary repeatability and contingency. Here, we review the experimental and comparative evidence from these studies. Replicate populations in evolutionary “replay” experiments often show parallel changes, especially in overall performance, although idiosyncratic outcomes show that the particulars of a lineage’s history can affect which of several evolutionary paths is taken. Comparative biologists have found many notable examples of convergent adaptation to similar conditions, but quantification of how frequently such convergence occurs is difficult. On balance, the evidence indicates that evolution tends to be surprisingly repeatable among closely related lineages, but disparate outcomes become more likely as the footprint of history grows deeper. Ongoing research on the structure of adaptive landscapes is providing additional insight into the interplay of fate and chance in the evolutionary process.

Modeling the contributions of chromosome segregation errors and aneuploidy to Saccharomyces hybrid sterility

Errors in meiosis can be important postzygotic barriers between different species. In Saccharomyces hybrids, chromosomal missegregation during meiosis I produces gametes with missing or extra chromosomes. Gametes with missing chromosomes are inviable, but we do not understand how extra chromosomes (disomies) influence hybrid gamete inviability. We designed a model predicting rates of missegregation in interspecific hybrid meioses assuming several different mechanisms of disomy tolerance, and compared predictions from the model with observations of sterility in hybrids between Saccharomyces yeast species. Sterility observations were consistent with the hypothesis that chromosomal missegregation causes hybrid sterility, and the model indicated that missegregation probabilities of 13-50% per chromosome can cause observed values of 90-99% hybrid sterility regardless of how cells tolerate disomies. Missing chromosomes in gametes are responsible for most infertility, but disomies may kill as many as 11% of the gametes produced by hybrids between Saccharomyces cerevisiae and Saccharomyces paradoxus. Copyright © 2017 John Wiley & Sons, Ltd.

The adaptive landscape of wildtype and glycosylation-deficient populations of the industrial yeast Pichia pastoris

Background

The effects of long-term environmental adaptation and the implications of major cellular malfunctions are still poorly understood for non-model but biotechnologically relevant species. In this study we performed a large-scale laboratory evolution experiment with 48 populations of the yeast Pichia pastoris in order to establish a general adaptive landscape upon long-term selection in several glucose-based growth environments. As a model for a cellular malfunction the implications of OCH1 mannosyltransferase knockout-mediated glycosylation-deficiency were analyzed.

Results

In-depth growth profiling of evolved populations revealed several instances of genotype-dependent growth trade-off/cross-benefit correlations in non-evolutionary growth conditions. On the genome level a high degree of mutational convergence was observed among independent populations. Environment-dependent mutational hotspots were related to osmotic stress-, Rim - and cAMP signaling pathways. In agreement with the observed growth phenotypes, our data also suggest diverging compensatory mutations in glycosylation-deficient populations. High osmolarity glycerol (HOG) pathway loss-of-functions mutations, including genes such as SSK2 and SSK4, represented a major adaptive strategy during environmental adaptation. However, genotype-specific HOG-related mutations were predominantly observed in opposing environmental conditions. Surprisingly, such mutations emerged during salt stress adaptation in OCH1 knockout populations and led to growth trade-offs in non-adaptive conditions that were distinct from wildtype HOG-mutants. Further environment-dependent mutations were identified for a hitherto uncharacterized species-specific Gal4-like transcriptional regulator involved in environmental sensing.

Conclusion

We show that metabolic constraints such as glycosylation-deficiency can contribute to evolution on the molecular level, even in non-diverging growth environments. Our dataset suggests universal adaptive mechanisms involving cellular stress response and cAMP/PKA signaling but also the existence of highly species-specific strategies involving unique transcriptional regulators, improving our biological understanding of distinct Ascomycetes species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3952-7) contains supplementary material, which is available to authorized users.

Intrinsic incompatibilities evolving as a by-product of divergent ecological selection: Considering them in empirical studies on divergence with gene flow

The possibility of intrinsic barriers to gene flow is often neglected in empirical research on local adaptation and speciation with gene flow, for example when interpreting patterns observed in genome scans. However, we draw attention to the fact that, even with gene flow, divergent ecological selection may generate intrinsic barriers involving both ecologically selected and other interacting loci. Mechanistically, the link between the two types of barriers may be generated by genes that have multiple functions (i.e., pleiotropy), and/or by gene interaction networks. Because most genes function in complex networks, and their evolution is not independent of other genes, changes evolving in response to ecological selection can generate intrinsic barriers as a by-product. A crucial question is to what extent such by-product barriers contribute to divergence and speciation-that is whether they stably reduce gene flow. We discuss under which conditions by-product barriers may increase isolation. However, we also highlight that, depending on the conditions (e.g., the amount of gene flow and the strength of selection acting on the intrinsic vs. the ecological barrier component), the intrinsic incompatibility may actually destabilize barriers to gene flow. In practice, intrinsic barriers generated as a by-product of divergent ecological selection may generate peaks in genome scans that cannot easily be interpreted. We argue that empirical studies on divergence with gene flow should consider the possibility of both ecological and intrinsic barriers. Future progress will likely come from work combining population genomic studies, experiments quantifying fitness and molecular studies on protein function and interactions.

Genomics of Compensatory Adaptation in Experimental Populations of Aspergillus nidulans

Knowledge of the number and nature of genetic changes responsible for adaptation is essential for understanding and predicting evolutionary trajectories. Here, we study the genomic basis of compensatory adaptation to the fitness cost of fungicide resistance in experimentally evolved strains of the filamentous fungus Aspergillus nidulans The original selection experiment tracked the fitness recovery of lines founded by an ancestral strain that was resistant to fludioxonil, but paid a fitness cost in the absence of the fungicide. We obtained whole-genome sequence data for the ancestral A. nidulans strain and eight experimentally evolved strains. We find that fludioxonil resistance in the ancestor was likely conferred by a mutation in histidine kinase nikA, part of the two-component signal transduction system of the high-osmolarity glycerol (HOG) stress response pathway. To compensate for the pleiotropic negative effects of the resistance mutation, the subsequent fitness gains observed in the evolved lines were likely caused by secondary modification of HOG pathway activity. Candidate genes for the compensatory fitness increases were significantly overrepresented by stress response functions, and some were specifically associated with the HOG pathway itself. Parallel evolution at the gene level was rare among evolved lines. There was a positive relationship between the predicted number of adaptive steps, estimated from fitness data, and the number of genomic mutations, determined by whole-genome sequencing. However, the number of genomic mutations was, on average, 8.45 times greater than the number of adaptive steps inferred from fitness data. This research expands our understanding of the genetic basis of adaptation in multicellular eukaryotes and lays out a framework for future work on the genomics of compensatory adaptation in A. nidulans.

Widespread Genetic Incompatibilities between First-Step Mutations during Parallel Adaptation of Saccharomyces cerevisiae to a Common Environment

Independently evolving populations may adapt to similar selection pressures via different genetic changes. The interactions between such changes, such as in a hybrid individual, can inform us about what course adaptation may follow and allow us to determine whether gene flow would be facilitated or hampered following secondary contact. We used Saccharomyces cerevisiae to measure the genetic interactions between first-step mutations that independently evolved in the same biosynthetic pathway following exposure to the fungicide nystatin. We found that genetic interactions are prevalent and predominantly negative, with the majority of mutations causing lower growth when combined in a double mutant than when alone as a single mutant (sign epistasis). The prevalence of sign epistasis is surprising given the small number of mutations tested and runs counter to expectations for mutations arising in a single biosynthetic pathway in the face of a simple selective pressure. Furthermore, in one third of pairwise interactions, the double mutant grew less well than either single mutant (reciprocal sign epistasis). The observation of reciprocal sign epistasis among these first adaptive mutations arising in the same genetic background indicates that partial postzygotic reproductive isolation could evolve rapidly between populations under similar selective pressures, even with only a single genetic change in each. The nature of the epistatic relationships was sensitive, however, to the level of drug stress in the assay conditions, as many double mutants became fitter than the single mutants at higher concentrations of nystatin. We discuss the implications of these results both for our understanding of epistatic interactions among beneficial mutations in the same biochemical pathway and for speciation.

Spiraling Complexity: A Test of the Snowball Effect in a Computational Model of RNA Folding

Genetic incompatibilities can emerge as a byproduct of genetic divergence. According to Dobzhansky and Muller, an allele that fixes in one population may be incompatible with an allele at a different locus in another population when the two alleles are brought together in hybrids. Orr showed that the number of Dobzhansky-Muller incompatibilities (DMIs) should accumulate faster than linearly-i.e., snowball-as two lineages diverge. Several studies have attempted to test the snowball effect using data from natural populations. One limitation of these studies is that they have focused on predictions of the Orr model, but not on its underlying assumptions. Here, we use a computational model of RNA folding to test both predictions and assumptions of the Orr model. Two populations are allowed to evolve in allopatry on a holey fitness landscape. We find that the number of inviable introgressions (an indicator for the number of DMIs) snowballs, but does so more slowly than expected. We show that this pattern is explained, in part, by the fact that DMIs can disappear after they have arisen, contrary to the assumptions of the Orr model. This occurs because DMIs become progressively more complex (i.e., involve alleles at more loci) as a result of later substitutions. We also find that most DMIs involve >2 loci, i.e., they are complex. Reproductive isolation does not snowball because DMIs do not act independently of each other. We conclude that the RNA model supports the central prediction of the Orr model that the number of DMIs snowballs, but challenges other predictions, as well as some of its underlying assumptions.

Fitness Trade-Offs Lead to Suppressor Tolerance in Yeast

Genetic variation among individuals within a population provides the raw material for phenotypic diversity upon which natural selection operates. Some given variants can act on multiple standing genomic variations simultaneously and release previously inaccessible phenotypes, leading to increased adaptive potential upon challenging environments. Previously, we identified such a variant related to a tRNA nonsense suppressor in yeast. When introduced into other genetic backgrounds, the suppressor led to an increased population phenotypic variance on various culture conditions, conferring background and environment specific selective advantages. Nonetheless, most isolates are intolerant to the suppressor on rich media due to a severe fitness cost. Here, we found that the tolerance to suppressor is related to a surprising level of fitness outburst, showing a trade-off effect to accommodate the cost of carrying the suppressor. To dissect the genetic basis of such trade-offs, we crossed strains with contrasting tolerance levels on rich media, and analyzed the fitness distribution patterns in the offspring. Combining quantitative tetrad analysis and bulk segregant analysis, we identified two genes, namely MKT1 and RGA1, involved in suppressor tolerance. We showed that alleles from the tolerant parent for both genes conferred a significant gain of fitness, which increased the suppressor tolerance. Our results present a detailed dissection of suppressor tolerance in yeast and provide insights into the molecular basis of trade-offs between fitness and evolutionary potential.

The Hidden Complexity of Mendelian Traits across Natural Yeast Populations

Mendelian traits are considered to be at the lower end of the complexity spectrum of heritable phenotypes. However, more than a century after the rediscovery of Mendel's law, the global landscape of monogenic variants, as well as their effects and inheritance patterns within natural populations, is still not well understood. Using the yeast Saccharomyces cerevisiae, we performed a species-wide survey of Mendelian traits across a large population of isolates. We generated offspring from 41 unique parental pairs and analyzed 1,105 cross/trait combinations. We found that 8.9% of the cases were Mendelian. Further tracing of causal variants revealed background-specific expressivity and modified inheritances, gradually transitioning from Mendelian to complex traits in 30% of the cases. In fact, when taking into account the natural population diversity, the hidden complexity of traits could be substantial, confounding phenotypic predictability even for simple Mendelian traits.

A Genetic Incompatibility Accelerates Adaptation in Yeast

During mismatch repair (MMR) MSH proteins bind to mismatches that form as the result of DNA replication errors and recruit MLH factors such as Mlh1-Pms1 to initiate excision and repair steps. Previously, we identified a negative epistatic interaction involving naturally occurring polymorphisms in the MLH1 and PMS1 genes of baker’s yeast. Here we hypothesize that a mutagenic state resulting from this negative epistatic interaction increases the likelihood of obtaining beneficial mutations that can promote adaptation to stress conditions. We tested this by stressing yeast strains bearing mutagenic (incompatible) and non-mutagenic (compatible) mismatch repair genotypes. Our data show that incompatible populations adapted more rapidly and without an apparent fitness cost to high salt stress. The fitness advantage of incompatible populations was rapid but disappeared over time. The fitness gains in both compatible and incompatible strains were due primarily to mutations in PMR1 that appeared earlier in incompatible evolving populations. These data demonstrate a rapid and reversible role (by mating) for genetic incompatibilities in accelerating adaptation in eukaryotes. They also provide an approach to link experimental studies to observational population genomics.

In nature, bacterial populations with high mutation rates can adapt faster to new environments by acquiring beneficial mutations. However, such populations also accumulate harmful mutations that reduce their fitness. We show that the model eukaryote baker’s yeast can use a similar mutator strategy to adapt to new environments. The mutator state that we observed resulted from an incompatibility involving two genes, MLH1 and PMS1, that work together to remove DNA replication errors through a spellchecking mismatch repair mechanism. This incompatibility can occur through mating between baker’s yeast from different genetic backgrounds, yielding mutator offspring containing an MLH1-PMS1 combination not present in either parent. Interestingly, these offspring adapted more rapidly to stress, compared to the parental strains, and did so without an overall loss in fitness. DNA sequencing analyses of baker’s yeast strains from across the globe support the presence of incompatible hybrid yeast strains in nature. These observations provide a powerful model to understand how the segregation of defects in DNA mismatch repair can serve as an effective strategy to enable eukaryotes to adapt to changing environments.

Temporal expression profiling identifies pathways mediating effect of causal variant on phenotype

Even with identification of multiple causal genetic variants for common human diseases, understanding the molecular processes mediating the causal variants’ effect on the disease remains a challenge. This understanding is crucial for the development of therapeutic strategies to prevent and treat disease. While static profiling of gene expression is primarily used to get insights into the biological bases of diseases, it makes differentiating the causative from the correlative effects difficult, as the dynamics of the underlying biological processes are not monitored. Using yeast as a model, we studied genome-wide gene expression dynamics in the presence of a causal variant as the sole genetic determinant, and performed allele-specific functional validation to delineate the causal effects of the genetic variant on the phenotype. Here, we characterized the precise genetic effects of a functional MKT1 allelic variant in sporulation efficiency variation. A mathematical model describing meiotic landmark events and conditional activation of MKT1 expression during sporulation specified an early meiotic role of this variant. By analyzing the early meiotic genome-wide transcriptional response, we demonstrate an MKT1-dependent role of novel modulators, namely, RTG1/3, regulators of mitochondrial retrograde signaling, and DAL82, regulator of nitrogen starvation, in additively effecting sporulation efficiency. In the presence of functional MKT1 allele, better respiration during early sporulation was observed, which was dependent on the mitochondrial retrograde regulator, RTG3. Furthermore, our approach showed that MKT1 contributes to sporulation independent of Puf3, an RNA-binding protein that steady-state transcription profiling studies have suggested to mediate MKT1-pleiotropic effects during mitotic growth. These results uncover interesting regulatory links between meiosis and mitochondrial retrograde signaling. In this study, we highlight the advantage of analyzing allele-specific transcriptional dynamics of mediating genes. Applications in higher eukaryotes can be valuable for inferring causal molecular pathways underlying complex dynamic processes, such as development, physiology and disease progression.

The causal path from a genetic variant to a complex phenotype such as disease progression is often not known. Studying gene expression variation is one approach to identify the mediating genes, however, it is difficult to distinguish causative from correlative genes. This becomes a challenge especially when studying developmental and physiological traits, since they involve dynamic processes contributing to the variation and only single static expression profiling is performed. As a proof of concept, we addressed this challenge here in yeast, by studying genome-wide gene expression in the presence of the causative polymorphism of MKT1 as the sole genetic variant, during the time phase when it contributes to sporulation efficiency variation. Our analysis during early sporulation identified mitochondrial retrograde signaling and nitrogen starvation as novel regulators, acting additively to regulate sporulation efficiency. Furthermore, we showed that PUF3, a known interactor of MKT1 had an independent role in sporulation. Our results highlight the role of differential mitochondrial signaling for efficient meiosis, providing insights into the factors regulating infertility. In addition, our study has implications for characterizing the molecular effects of causal genetic variants on dynamic biological processes during development and disease progression.

Rapid selective sweep of pre-existing polymorphisms and slow fixation of new mutations in experimental evolution of Desulfovibrio vulgaris

To investigate the genetic basis of microbial evolutionary adaptation to salt (NaCl) stress, populations of Desulfovibrio vulgaris Hildenborough (DvH), a sulfate-reducing bacterium important for the biogeochemical cycling of sulfur, carbon and nitrogen, and potentially the bioremediation of toxic heavy metals and radionuclides, were propagated under salt stress or non-stress conditions for 1200 generations. Whole-genome sequencing revealed 11 mutations in salt stress-evolved clone ES9-11 and 14 mutations in non-stress-evolved clone EC3-10. Whole-population sequencing data suggested the rapid selective sweep of the pre-existing polymorphisms under salt stress within the first 100 generations and the slow fixation of new mutations. Population genotyping data demonstrated that the rapid selective sweep of pre-existing polymorphisms was common in salt stress-evolved populations. In contrast, the selection of pre-existing polymorphisms was largely random in EC populations. Consistently, at 100 generations, stress-evolved population ES9 showed improved salt tolerance, namely increased growth rate (2.0-fold), higher biomass yield (1.8-fold) and shorter lag phase (0.7-fold) under higher salinity conditions. The beneficial nature of several mutations was confirmed by site-directed mutagenesis. All four tested mutations contributed to the shortened lag phases under higher salinity condition. In particular, compared with the salt tolerance improvement in ES9-11, a mutation in a histidine kinase protein gene lytS contributed 27% of the growth rate increase and 23% of the biomass yield increase while a mutation in hypothetical gene DVU2472 contributed 24% of the biomass yield increase. Our results suggested that a few beneficial mutations could lead to dramatic improvements in salt tolerance.

Post-zygotic sterility and cytonuclear compatibility limits in S. cerevisiae xenomitochondrial cybrids

Nucleo-mitochondrial interactions, particularly those determining the primary divergence of biological species, can be studied by means of xenomitochondrial cybrids, which are cells where the original mitochondria are substituted by their counterparts from related species. Saccharomyces cerevisiae cybrids are prepared simply by the mating of the ρ(0) strain with impaired karyogamy and germinating spores from other Saccharomyces species and fall into three categories. Cybrids with compatible mitochondrial DNA (mtDNA) from Saccharomyces paradoxus CBS 432 and Saccharomyces cariocanus CBS 7994 are metabolically and genetically similar to cybrids containing mtDNA from various S. cerevisiae. Cybrids with mtDNA from other S. paradoxus strains, S. cariocanus, Saccharomyces kudriavzevii, and Saccharomyces mikatae require a period of adaptation to establish efficient oxidative phosphorylation. They exhibit a temperature-sensitive phenotype, slower growth rate on a non-fermentable carbon source and a long lag phase after the shift from glucose. Their decreased respiration capacity and reduced cytochrome aa3 content is associated with the inefficient splicing of cox1I3β, the intron found in all Saccharomyces species but not in S. cerevisiae. The splicing defect is compensated in cybrids by nuclear gain-of-function and can be alternatively suppressed by overexpression of MRP13 gene for mitochondrial ribosomal protein or the MRS2, MRS3, and MRS4 genes involved in intron splicing. S. cerevisiae with Saccharomyces bayanus mtDNA is unable to respire and the growth on ethanol-glycerol can be restored only after mating to some mit (-) strains. The nucleo-mitochondrial compatibility limit of S. cerevisiae and other Saccharomyces was set between S. kudriavzevii and S. bayanus at the divergence from S. cerevisiae about 15 MYA. The MRS1-cox1 S. cerevisiae/S. paradoxus cytonuclear Dobzhansky-Muller pair has a neglible impact on the separation of species since its imperfection is compensated for by gain-of-function mutation.

Estimation of epistatic variance components and heritability in founder populations and crosses

Genetic association studies have explained only a small proportion of the estimated heritability of complex traits, leaving the remaining heritability “missing.” Genetic interactions have been proposed as an explanation for this, because they lead to overestimates of the heritability and are hard to detect. Whether this explanation is true depends on the proportion of variance attributable to genetic interactions, which is difficult to measure in outbred populations. Founder populations exhibit a greater range of kinship than outbred populations, which helps in fitting the epistatic variance. We extend classic theory to founder populations, giving the covariance between individuals due to epistasis of any order. We recover the classic theory as a limit, and we derive a recently proposed estimator of the narrow sense heritability as a corollary. We extend the variance decomposition to include dominance. We show in simulations that it would be possible to estimate the variance from pairwise interactions with samples of a few thousand from strongly bottlenecked human founder populations, and we provide an analytical approximation of the standard error. Applying these methods to 46 traits measured in a yeast (Saccharomyces cerevisiae) cross, we estimate that pairwise interactions explain 10% of the phenotypic variance on average and that third- and higher-order interactions explain 14% of the phenotypic variance on average. We search for third-order interactions, discovering an interaction that is shared between two traits. Our methods will be relevant to future studies of epistatic variance in founder populations and crosses.

Widespread host-dependent hybrid unfitness in the pea aphid species complex

Linking adaptive divergence to hybrid unfitness is necessary to understand the ecological factors contributing to reproductive isolation and speciation. To date, this link has been demonstrated in few model systems, most of which encompass ecotypes that occupy relatively early stages in the speciation process. Here we extend these studies by assessing how host-plant adaptation conditions hybrid fitness in the pea aphid, Acyrthosiphon pisum. We made crosses between and within five pea aphid biotypes adapted to different host plants and representing various stages of divergence within the complex. Performance of F1 hybrids and nonhybrids was assessed on a "universal" host that is favorable to all pea aphid biotypes in laboratory conditions. Although hybrids performed equally well as nonhybrids on the universal host, their performance was much lower than nonhybrids on the natural hosts of their parental populations. Hence, hybrids, rather than being intrinsically deficient, are maladapted to their parents' hosts. Interestingly, the impact of this maladaptation was stronger in certain hybrids from crosses involving the most divergent biotype, suggesting that host-dependent postzygotic isolation has continued to evolve late in divergence. Even though host-independent deficiencies are not excluded, hybrid maladaptation to parental hosts supports the hypothesis of ecological speciation in this complex.

The underlying structure of adaptation under strong selection in 12 experimental yeast populations

The aims of this study were to determine (i) whether adaptation under strong selection occurred through mutations in a narrow target of one or a few nucleotide sites or a broad target of numerous sites and (ii) whether the programs of adaptation previously observed from three experimental populations were unique or shared among populations that underwent parallel evolution. We used archived population samples from a previous study, representing 500 generations of experimental evolution in 12 populations under strong selection, 6 populations in a high-salt environment and 6 populations in a low-glucose environment. Each set of six populations included four with sexual reproduction and two with exclusively asexual reproduction. Populations were sampled as resequenced genomes of 115 individuals and as bulk samples from which frequencies of mutant alleles were estimated. In a high-salt environment, a broad target of 11 mutations within the proton exporter, PMA1, was observed among the six populations, in addition to expansions of the ENA gene cluster. This pattern was shared among populations that underwent parallel evolution. In a low-glucose environment, two programs of adaptation were observed. The originally observed pattern of mutation in MDS3/MKT1 in population M8 was a narrow target of a single nucleotide, unique to this population. Among the other five populations, the three mutations were shared in a broad target, sensing/signaling genes RAS1 and RAS2. RAS1/RAS2 mutations were not observed in the high-salt populations; PMA1 mutations were observed only in a high-salt environment.

Chromosomal rearrangements as a major mechanism in the onset of reproductive isolation in Saccharomyces cerevisiae

Understanding the molecular basis of how reproductive isolation evolves between individuals from the same species offers valuable insight into patterns of genetic differentiation as well as the onset of speciation [1, 2]. The yeast Saccharomyces cerevisiae constitutes an ideal model partly due to its vast ecological range, high level of genetic diversity [3-6], and laboratory-amendable sexual reproduction. Between S. cerevisiae and its sibling species in the Saccharomyces sensu stricto complex, reproductive isolation acts postzygotically and could be attributed to chromosomal rearrangements [7], cytonuclear incompatibility [8, 9], and antirecombination [10, 11], although the implication of these mechanisms at the incipient stage of speciation remains unclear due to further divergence in the nascent species. Recently, several studies assessed the onset of intraspecific reproductive isolation in S. cerevisiae by evaluating the effect of the mismatch repair system [12-14] or by fostering incipient speciation using the same initial genetic background [15-18]. Nevertheless, the overall genetic diversity within this species was largely overlooked, and no systematic evaluation has been performed. Here, we carried out the first species-wide survey for postzygotic reproductive isolation in S. cerevisiae. We crossed 60 natural isolates sampled from diverse niches with the reference strain S288c and identified 16 cases of reproductive isolation with reduced offspring viabilities ranging from 44% to 86%. Using different mapping strategies, we identified reciprocal translocations in a large fraction of all isolates surveyed, indicating that large-scale chromosomal rearrangements might play a major role in the onset of reproductive isolation in this species.

The reproducibility of adaptation in the light of experimental evolution with whole genome sequencing

A key question in evolutionary biology is the reproducibility of adaptation. This question can now be quantitatively analyzed using experimental evolution coupled to whole genome sequencing (WGS). With complete sequence data, one can assess convergence among replicate populations. In turn, convergence reflects the action of natural selection and also the breadth of the field of possible adaptive solutions. That is, it provides insight into how many genetic solutions or adaptive paths may lead to adaptation in a given environment. Convergence is both a property of an adaptive landscape and, reciprocally, a tool to study that landscape. In this chapter we present the links between convergence and the properties of adaptive landscapes with respect to two types of microbial experimental evolution. The first tries to reconstruct a full adaptive landscape using a handful of carefully identified mutations (the reductionist approach), while the second uses WGS of replicate experiments to infer properties of the adaptive landscape. Reductionist approaches have highlighted the importance of epistasis in shaping the adaptive landscape, but have also uncovered a wide diversity of landscape architectures. The WGS approach has uncovered a very high diversity of beneficial mutations that affect a limited set of genes or functions and also suggests some shortcomings of the reductionist approach. We conclude that convergence may be better defined at an integrated level, such as the genic level or even at a phenotypic level, and that integrated mechanistic models derived from systems biology may offer an interesting perspective for the analysis of convergence at all levels.

Genomics and the Evolution of Phenotypic Traits

Evolutionary genetics has entered an unprecedented era of discovery, catalyzed in large part by the development of technologies that provide information about genome sequence and function. An important benefit is the ability to move beyond a handful of model organisms in lab settings to identify the genetic basis for evolutionarily interesting traits in many organisms in natural settings. Other benefits are the abilities to identify causal mutations and validate their phenotypic consequences more readily and in many more species. Genomic technologies have reinvigorated interest in some of the most fundamental and persistent questions in evolutionary genetics, revealed previously unsuspected evolutionary phenomena, and opened the door to a wide range of new questions.

Evolution, selection and isolation: a genomic view of speciation in fungal plant pathogens

895 I. 895 II. 896 III. 898 IV. 900 V. 902 VI. 904 VII. 905 905 References 905 SUMMARY: Speciation of fungal plant pathogens has been associated with host jumps, host domestication, clonal divergence, and hybridization. Although we have substantial insight into the speciation histories of several important plant pathogens, we still know very little about the underlying genetics of reproductive isolation. Studies in Saccharomyces cerevisiae, Neurospora crassa, and nonfungal model systems illustrate that reproductive barriers can evolve by different mechanisms, including genetic incompatibilities between neutral and adaptive substitutions, reinforcement selection, and chromosomal rearrangements. Advances in genome sequencing and sequence analyses provide a new framework to identify those traits that have driven the divergence of populations or caused reproductive isolation between species of fungal plant pathogens. These traits can be recognized based on signatures of strong divergent selection between species or through the association of allelic combination conferring hybrid inferiority. Comparative genome analyses also provide information about the contribution of genome rearrangements to speciation. This is particularly relevant for species of fungal pathogens with extreme levels of genomic rearrangements and within-species genome plasticity.

Epistatic interactions determine the mutational pathways and coexistence of lineages in clonal Escherichia coli populations

Understanding how diversity emerges in a single niche is not fully understood. Rugged fitness landscapes and epistasis between beneficial mutations could explain coexistence among emerging lineages. To provide an experimental test of this notion, we investigated epistasis among four pleiotropic mutations in rpoS, mglD, malT, and hfq present in two coexisting lineages that repeatedly fixed in experimental populations of Escherichia coli. The mutations were transferred into the ancestral background individually or in combination of double or triple alleles. The combined competitive fitness of two or three beneficial mutations from the same lineage was consistently lower than the sum of the competitive fitness of single mutants--a clear indication of negative epistasis within lineages. We also found sign epistasis (i.e., the combined fitness of two beneficial mutations lower than the ancestor), not only from two different lineages (i.e., hfq and rpoS) but also from the same lineage (i.e., mglD and malT). The sign epistasis between loci of different lineages indeed indicated a rugged fitness landscape, providing an epistatic explanation for the coexistence of distinct rpoS and hfq lineages in evolving populations. The negative and sign epistasis between beneficial mutations within the same lineage can further explain the order of mutation acquisition.

Genetic and genomic architecture of the evolution of resistance to antifungal drug combinations

The evolution of drug resistance in fungal pathogens compromises the efficacy of the limited number of antifungal drugs. Drug combinations have emerged as a powerful strategy to enhance antifungal efficacy and abrogate drug resistance, but the impact on the evolution of drug resistance remains largely unexplored. Targeting the molecular chaperone Hsp90 or its downstream effector, the protein phosphatase calcineurin, abrogates resistance to the most widely deployed antifungals, the azoles, which inhibit ergosterol biosynthesis. Here, we evolved experimental populations of the model yeast Saccharomyces cerevisiae and the leading human fungal pathogen Candida albicans with azole and an inhibitor of Hsp90, geldanamycin, or calcineurin, FK506. To recapitulate a clinical context where Hsp90 or calcineurin inhibitors could be utilized in combination with azoles to render resistant pathogens responsive to treatment, the evolution experiment was initiated with strains that are resistant to azoles in a manner that depends on Hsp90 and calcineurin. Of the 290 lineages initiated, most went extinct, yet 14 evolved resistance to the drug combination. Drug target mutations that conferred resistance to geldanamycin or FK506 were identified and validated in five evolved lineages. Whole-genome sequencing identified mutations in a gene encoding a transcriptional activator of drug efflux pumps, PDR1, and a gene encoding a transcriptional repressor of ergosterol biosynthesis genes, MOT3, that transformed azole resistance of two lineages from dependent on calcineurin to independent of this regulator. Resistance also arose by mutation that truncated the catalytic subunit of calcineurin, and by mutation in LCB1, encoding a sphingolipid biosynthetic enzyme. Genome analysis revealed extensive aneuploidy in four of the C. albicans lineages. Thus, we identify molecular determinants of the transition of azole resistance from calcineurin dependence to independence and establish multiple mechanisms by which resistance to drug combinations evolves, providing a foundation for predicting and preventing the evolution of drug resistance.

Fungal infections are a leading cause of mortality worldwide and are difficult to treat due to the limited number of antifungal drugs, whose effectiveness is compromised by the emergence of drug resistance. A powerful strategy to combat drug resistance is combination therapy. Inhibiting the molecular chaperone Hsp90 or its downstream effector calcineurin cripples fungal stress responses and abrogates drug resistance. Here we provide the first analysis of the genetic and genomic changes that underpin the evolution of resistance to antifungal drug combinations in the leading human fungal pathogen, Candida albicans, and model yeast, Saccharomyces cerevisiae. We evolved experimental populations with combinations of inhibitors of Hsp90 or calcineurin and the most widely used antifungal in the clinic, the azoles, which inhibit ergosterol biosynthesis. We harnessed whole-genome sequencing to identify diverse resistance mutations among the 14 lineages that evolved resistance to the drug combination. These included mutations in genes encoding the drug targets, a transcriptional regulator of multidrug transporters, a transcriptional repressor of ergosterol biosynthesis enzymes, and a regulator of sphingolipid biosynthesis. We also identified extensive aneuploidies in several C. albicans lineages. Our study reveals multiple mechanisms by which resistance to drug combination can evolve, suggesting new strategies to combat drug resistance.

The Loci of repeated evolution: a catalog of genetic hotspots of phenotypic variation

What is the nature of the genetic changes underlying phenotypic evolution? We have catalogued 1008 alleles described in the literature that cause phenotypic differences among animals, plants, and yeasts. Surprisingly, evolution of similar traits in distinct lineages often involves mutations in the same gene ("gene reuse"). This compilation yields three important qualitative implications about repeated evolution. First, the apparent evolution of similar traits by gene reuse can be traced back to two alternatives, either several independent causative mutations or a single original mutational event followed by sorting processes. Second, hotspots of evolution-defined as the repeated occurrence of de novo mutations at orthologous loci and causing similar phenotypic variation-are omnipresent in the literature with more than 100 examples covering various levels of analysis, including numerous gain-of-function events. Finally, several alleles of large effect have been shown to result from the aggregation of multiple small-effect mutations at the same hotspot locus, thus reconciling micromutationist theories of adaptation with the empirical observation of large-effect variants. Although data heterogeneity and experimental biases prevented us from extracting quantitative trends, our synthesis highlights the existence of genetic paths of least resistance leading to viable evolutionary change.

Indirect evolution of hybrid lethality due to linkage with selected locus in Mimulus guttatus

Most species are superbly and intricately adapted to the environments in which they live. Adaptive evolution by natural selection is the primary force shaping biological diversity. Differences between closely related species in ecologically selected characters such as habitat preference, reproductive timing, courtship behavior, or pollinator attraction may prevent interbreeding in nature, causing reproductive isolation. But does ecological adaptation cause reproductive incompatibilities such as hybrid sterility or lethality? Although several genes causing hybrid incompatibilities have been identified, there is intense debate over whether the genes that contribute to ecological adaptations also cause hybrid incompatibilities. Thirty years ago, a genetic study of local adaptation to copper mine soils in the wildflower Mimulus guttatus identified a locus that appeared to cause copper tolerance and hybrid lethality in crosses to other populations. But do copper tolerance and hybrid lethality have the same molecular genetic basis? Here we show, using high-resolution genome mapping, that copper tolerance and hybrid lethality are not caused by the same gene but are in fact separately controlled by two tightly linked loci. We further show that selection on the copper tolerance locus indirectly caused the hybrid incompatibility allele to go to high frequency in the copper mine population because of hitchhiking. Our results provide a new twist on Darwin's original supposition that hybrid incompatibilities evolve as an incidental by-product of ordinary adaptation to the environment.

Genetic background affects epistatic interactions between two beneficial mutations

The phenotypic effect of mutations can depend on their genetic background, a phenomenon known as epistasis. Many experimental studies have found that epistasis is pervasive, and some indicate that it may follow a general pattern dependent on the fitness effect of the interacting mutations. These studies have, however, typically examined the effect of interactions between a small number of focal mutations in a single genetic background. Here, we extend this approach by considering how the interaction between two beneficial mutations that were isolated from a population of laboratory evolved Escherichia coli changes when they are added to divergent natural isolate strains of E. coli. We find that interactions between the focal mutations and the different genetic backgrounds are common. Moreover, the pair-wise interaction between the focal mutations also depended on their genetic background, being more negative in backgrounds with higher absolute fitness. Together, our results indicate the presence of interactions between focal mutations, but also caution that these interactions depend quantitatively on the wider genetic background.

The tempo and modes of evolution of reproductive isolation in fungi

Reproductive isolation is an essential ingredient of speciation, and much has been learned in recent years about the evolution of reproductive isolation and the genetics of reproductive barriers in animals and plants. Fungi have been neglected on these aspects, despite being tractable model eukaryotes. Here, we used a model fitting approach to look at the importance of different barriers to gene flow to explain the decrease of reproductive compatibility with genetic distance in fungi. We found support for the occurrence of reinforcement in the presyngamy compatibility among basidiomycetes. In contrast, no evidence for reinforcement was detected in ascomycetes, concurring with the idea that host/habitat adaptation in this group can pleiotropically cause reproductive isolation. We found no evidence of a snowballing accumulation of postsyngamic reproductive incompatibilities in either ascomycetes or the complex of anther smut fungi. Together with previous studies, our results suggest that ecologically based barriers to gene flow and karyotypic differences may have an important role in hybrid inviability and sterility in fungi. Interestingly, hybrid sterility appeared to evolve faster than hybrid inviability in fungi.

Large-scale candidate gene scan reveals the role of chemoreceptor genes in host plant specialization and speciation in the pea aphid

Understanding the drivers of speciation is critical to interpreting patterns of biodiversity. The identification of the genetic changes underlying adaptation and reproductive isolation is necessary to link barriers to gene flow to the causal origins of divergence. Here, we present a novel approach to the genetics of speciation, which should complement the commonly used approaches of quantitative trait locus mapping and genome-wide scans for selection. We present a large-scale candidate gene approach by means of sequence capture, applied to identifying the genetic changes underlying reproductive isolation in the pea aphid, a model system for the study of ecological speciation. Targeted resequencing enabled us to scale up the candidate gene approach, specifically testing for the role of chemosensory gene families in host plant specialization. Screening for the signature of divergence under selection at 172 candidate and noncandidate loci, we revealed a handful of loci that show high levels of differentiation among host races, which almost all correspond to odorant and gustatory receptor genes. This study offers the first indication that some chemoreceptor genes, often tightly linked together in the genome, could play a key role in local adaptation and reproductive isolation in the pea aphid and potentially other phytophagous insects. Our approach opens a new route toward the functional genomics of ecological speciation.

A stochastic model for the development of Bateson-Dobzhansky-Muller incompatibilities that incorporates protein interaction networks

Speciation is characterized by the development of reproductive isolating barriers between diverging groups. Intrinsic post-zygotic barriers of the type envisioned by Bateson, Dobzhansky, and Muller are deleterious epistatic interactions among loci that reduce hybrid fitness, leading to reproductive isolation. The first formal population genetic model of the development of these barriers was published by Orr in 1995, and here we develop a more general model of this process by incorporating finite protein-protein interaction networks, which reduce the probability of deleterious interactions in vivo. Our model shows that the development of deleterious interactions is limited by the density of the protein-protein interaction network. We have confirmed our analytical predictions of the number of possible interactions given the number of allele substitutions by using simulations on the Saccharomyces cerevisiae protein-protein interaction network. These results allow us to define the rate at which deleterious interactions are expected to form, and hence the speciation rate, for any protein-protein interaction network.

Evolutionary insight from whole-genome sequencing of experimentally evolved microbes

Experimental evolution (EE) combined with whole-genome sequencing (WGS) has become a compelling approach to study the fundamental mechanisms and processes that drive evolution. Most EE-WGS studies published to date have used microbes, owing to their ease of propagation and manipulation in the laboratory and relatively small genome sizes. These experiments are particularly suited to answer long-standing questions such as: How many mutations underlie adaptive evolution, and how are they distributed across the genome and through time? Are there general rules or principles governing which genes contribute to adaptation, and are certain kinds of genes more likely to be targets than others? How common is epistasis among adaptive mutations, and what does this reveal about the variety of genetic routes to adaptation? How common is parallel evolution, where the same mutations evolve repeatedly and independently in response to similar selective pressures? Here, we summarize the significant findings of this body of work, identify important emerging trends and propose promising directions for future research. We also outline an example of a computational pipeline for use in EE-WGS studies, based on freely available bioinformatics tools.

Systematic detection of epistatic interactions based on allele pair frequencies

Epistatic genetic interactions are key for understanding the genetic contribution to complex traits. Epistasis is always defined with respect to some trait such as growth rate or fitness. Whereas most existing epistasis screens explicitly test for a trait, it is also possible to implicitly test for fitness traits by searching for the over- or under-representation of allele pairs in a given population. Such analysis of imbalanced allele pair frequencies of distant loci has not been exploited yet on a genome-wide scale, mostly due to statistical difficulties such as the multiple testing problem. We propose a new approach called Imbalanced Allele Pair frequencies (ImAP) for inferring epistatic interactions that is exclusively based on DNA sequence information. Our approach is based on genome-wide SNP data sampled from a population with known family structure. We make use of genotype information of parent-child trios and inspect 3×3 contingency tables for detecting pairs of alleles from different genomic positions that are over- or under-represented in the population. We also developed a simulation setup which mimics the pedigree structure by simultaneously assuming independence of the markers. When applied to mouse SNP data, our method detected 168 imbalanced allele pairs, which is substantially more than in simulations assuming no interactions. We could validate a significant number of the interactions with external data, and we found that interacting loci are enriched for genes involved in developmental processes.