The high-throughput yeast deletion fitness data and the theories of dominance

Schrödinger's yeast: the challenge of using transformation to compare fitness among Saccharomyces cerevisiae that differ in ploidy or zygosity

How the number of genome copies modifies the effect of random mutations remains poorly known. In yeast, researchers have investigated these effects for knock-out or other large-effect mutations, but have not accounted for differences at the mating-type locus. We set out to compare fitness differences among strains that differ in ploidy and/or zygosity using a panel of spontaneously arising mutations acquired in haploid yeast from a previous study. To ensure no genetic differences, even at the mating-type locus, we embarked on a series of transformations, which first sterilized and then temporarily introduced plasmid-borne mating types. Despite these attempts to equalize the haplotypes, fitness variation introduced during transformation swamped the differences among the original mutation-accumulation lines. While colony size looked normal, we observed a bi-modality in the maximum growth rate of our transformed yeast and determined that many of the slow growing lines were respiratory deficient ("petite"). Not previously reported, we found that yeast that were TID1/RDH54 knockouts were less likely to become petite. Even for lines with the same petite status, however, we found no correlation in fitness between the two replicate transformations performed. These results pose a challenge for any study using transformation to measure the fitness effect of genetic differences among strains. By attempting to hold haplotypes constant, we introduced more mutations that overwhelmed our ability to measure fitness differences between the genetic states. In this study, we transformed over one hundred different lines of yeast, using two independent transformations, and found that this common laboratory procedure can cause large changes to the microbe studied. Our study provides a cautionary tale of the need to use multiple transformants in fitness assays.

Evolution of the spinach sex-linked region within a rarely recombining pericentromeric region

Sex chromosomes have evolved independently in many different plant lineages. Here, we describe reference genomes for spinach (Spinacia oleracea) X and Y haplotypes by sequencing homozygous XX females and YY males. The long arm of 185-Mb chromosome 4 carries a 13-Mb X-linked region (XLR) and 24.1-Mb Y-linked region (YLR), of which 10 Mb is Y specific. We describe evidence that this reflects insertions of autosomal sequences creating a "Y duplication region" or "YDR" whose presence probably directly reduces genetic recombination in the immediately flanking regions, although both the X and Y sex-linked regions are within a large pericentromeric region of chromosome 4 that recombines rarely in meiosis of both sexes. Sequence divergence estimates using synonymous sites indicate that YDR genes started diverging from their likely autosomal progenitors about 3 MYA, around the time when the flanking YLR stopped recombining with the XLR. These flanking regions have a higher density of repetitive sequences in the YY than the XX assembly and include slightly more pseudogenes compared with the XLR, and the YLR has lost about 11% of the ancestral genes, suggesting some degeneration. Insertion of a male-determining factor would have caused Y linkage across the entire pericentromeric region, creating physically small, highly recombining, terminal pseudoautosomal regions. These findings provide a broader understanding of the origin of sex chromosomes in spinach.

Fitness dependence preserves selection for recombination across diverse mixed mating strategies

Meiotic recombination and the factors affecting its rate and fate in nature have inspired many studies in theoretical evolutionary biology. Classical theoretical models have inferred that recombination can be favored under a rather restricted parameter range. Thus, the ubiquity of recombination in nature remains an open question. However, these models assumed constant recombination with an equal rate across all individuals within the population, whereas empirical evidence suggests that recombination may display certain sensitivity to ecological stressors and/or genotype fitness. Models assuming condition-dependent recombination show that such a strategy can often be favored over constant recombination. Moreover, in our recent model with panmictic populations subjected to purifying selection, fitness-dependent recombination was quite often favored even when any constant recombination was rejected. By using numerical modeling, we test whether such a 'recombination-rescuing potential' of fitness dependence holds also beyond panmixia, given the recognized effect of mating strategy on the evolution of recombination. We show that deviations from panmixia generally increase the recombination-rescuing potential of fitness dependence, with the strongest effect under intermediate selfing or high clonality. We find that under partial clonality, the evolutionary advantage of fitness-dependent recombination is determined mostly by selection against heterozygotes and additive-by-additive epistasis, while under partial selfing, additive-by-dominance epistasis is also a driver.

When and how do sex-linked regions become sex chromosomes?

The attention given to heteromorphism and genetic degeneration of "classical sex chromosomes" (Y chromosomes in XY systems, and the W in ZW systems that were studied first and are best described) has perhaps created the impression that the absence of recombination between sex chromosomes is inevitable. I here argue that continued recombination is often to be expected, that absence of recombination is surprising and demands further study, and that the involvement of selection in reduced recombination is not yet well understood. Despite a long history of investigations of sex chromosome pairs, there is a need for more quantitative approaches to studying sex-linked regions. I describe a scheme to help understand the relationships between different properties of sex-linked regions. Specifically, I focus on their sizes (differentiating between small regions and extensive fully sex-linked ones), the times when they evolved, and their differentiation, and review studies using DNA sequencing in nonmodel organisms that are providing information about the processes causing these properties.

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

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

Strong dominance of functional alleles over gene deletions in both intensely growing and deeply starved yeast cells

Previous studies with diploid yeast have shown that the deletion of one allele at a single locus typically has little impact on fitness under conditions promoting fast growth. Here, we confirm and quantify this finding. The strong dominance of functional over nonfunctional alleles is predicted by the metabolic control theory which assumes that the cell is a system of metabolic fluxes and that the total metabolic rate is equivalent to fitness. To test whether these requirements are critical, we tested dominance under conditions of long-term starvation when metabolism is low and thus the metabolic activities of proteins are likely inadequate or imbalanced. More fundamentally, the central assumption of the model, that high metabolic rate translates into high fitness, appears implausible. Contrary to these conjectures, we found that the mean rate of survival of starving heterozygotes was affected only slightly more than was the mean rate of growth under good conditions. Under none of the two treatments the central prediction of the model, that fitness of heterozygous strains is higher for the enzymatic proteins than for nonenzymatic ones, was confirmed. Our data add to growing uncertainty whether the metabolic control theory is sufficient to explain the remarkable ubiquity of strong genetic dominance.

Intrinsic differences between males and females determine sex-specific consequences of inbreeding

BACKGROUND

Inbreeding increases homozygosity and exposes deleterious recessive alleles, generally decreasing the fitness of inbred individuals. Interestingly, males and females are usually affected differently by inbreeding, though the more vulnerable sex depends on the species and trait measured.

RESULTS

We used the soil-dwelling nematode Caenorhabditis remanei to examine sex-specific inbreeding depression across nine lineages, five levels of inbreeding, and hundreds of thousands of progeny. Female nematodes consistently suffered greater fitness losses than their male counterparts, especially at high levels of inbreeding.

CONCLUSIONS

These results suggest that females experience stronger selection on genes contributing to reproductive traits. Inbreeding depression in males may be further reduced by sex chromosome hemizygosity, which affects the dominance of some mutations, as well as by the absence of sexual selection. Determining the relative contributions of sex-specific expression, genes on the sex chromosomes, and the environment they are filtered through-including opportunities for sexual selection-may explain the frequent though inconsistent records of sex differences in inbreeding depression, along with their implications for conservation and the evolution of mating systems.

Dominance from the perspective of gene-gene and gene-chemical interactions

In this study, we used genetic interaction (GI) and gene-chemical interaction (GCI) data to compare mutations with different dominance phenotypes. Our analysis focused primarily on Saccharomyces cerevisiae, where haploinsufficient genes (HI; genes with dominant loss-of-function mutations) were found to be participating in gene expression processes, namely, the translation and regulation of gene transcription. Non-ribosomal HI genes (mainly regulators of gene transcription) were found to have more GIs and GCIs than haplosufficient (HS) genes. Several properties seem to lead to the enrichment of interactions, most notably, the following: importance, pleiotropy, gene expression level and gene expression variation. Importantly, after these properties were appropriately considered in the analysis, the correlation between dominance and GI/GCI degrees was still observed. Strikingly, for the GCIs of heterozygous strains, haploinsufficiency was the only property significantly correlated with the number of GCIs. We found ribosomal HI genes to be depleted in GIs/GCIs. This finding can be explained by their high variation in gene expression under different genetic backgrounds and environmental conditions. We observed the same distributions of GIs among non-ribosomal HI, ribosomal HI and HS genes in three other species: Schizosaccharomyces pombe, Drosophila melanogaster and Homo sapiens. One potentially interesting exception was the lack of significant differences in the degree of GIs between non-ribosomal HI and HS genes in Schizosaccharomyces pombe.

Enhancer Runaway and the Evolution of Diploid Gene Expression

Evidence is mounting that the evolution of gene expression plays a major role in adaptation and speciation. Understanding the evolution of gene regulatory regions is indeed an essential step in linking genotypes and phenotypes and in understanding the molecular mechanisms underlying evolutionary change. The common view is that expression traits (protein folding, expression timing, tissue localization and concentration) are under natural selection at the individual level. Here, we use a theoretical approach to show that, in addition, in diploid organisms, enhancer strength (i.e., the ability of enhancers to activate transcription) may increase in a runaway process due to competition for expression between homologous enhancer alleles. These alleles may be viewed as self-promoting genetic elements, as they spread without conferring a benefit at the individual level. They gain a selective advantage by getting associated to better genetic backgrounds: deleterious mutations are more efficiently purged when linked to stronger enhancers. This process, which has been entirely overlooked so far, may help understand the observed overrepresentation of cis-acting regulatory changes in between-species phenotypic differences, and sheds a new light on investigating the contribution of gene expression evolution to adaptation.

With the advent of new sequencing technologies, the evolution of gene expression regulation is becoming a subject of intensive research. In this paper, we report an entirely new phenomenon acting on the evolution of gene regulatory sequences. We show that in a small genomic region around genes there is a selection pressure to increase expression, such that stronger enhancers are favored. This leads to an open-ended escalation of enhancer strength. This outcome is not a particular case and we expect it to occur for all genes in nearly all eukaryotic diploid organisms. We also show that this escalation is not stopped by stabilizing selection on expression profiles. Indeed, regulators may coevolve to maintain optimal phenotypes despite the enhancer strength escalation. This widespread phenomenon can significantly shift our understanding of gene regulatory regions and opens a wide array of possible tests.

Effects of Interference Between Selected Loci on the Mutation Load, Inbreeding Depression, and Heterosis

A classical prediction from single-locus models is that inbreeding increases the efficiency of selection against partially recessive deleterious alleles (purging), thereby decreasing the mutation load and level of inbreeding depression. However, previous multilocus simulation studies found that increasing the rate of self-fertilization of individuals may not lead to purging and argued that selective interference among loci causes this effect. In this article, I derive simple analytical approximations for the mutation load and inbreeding depression, taking into account the effects of interference between pairs of loci. I consider two classical scenarios of nonrandomly mating populations: a single population undergoing partial selfing and a subdivided population with limited dispersal. In the first case, correlations in homozygosity between loci tend to reduce mean fitness and increase inbreeding depression. These effects are stronger when deleterious alleles are more recessive, but only weakly depend on the strength of selection against deleterious alleles and on recombination rates. In subdivided populations, interference increases inbreeding depression within demes, but decreases heterosis between demes. Comparisons with multilocus, individual-based simulations show that these analytical approximations are accurate as long as the effects of interference stay moderate, but fail for high deleterious mutation rates and low dominance coefficients of deleterious alleles.

Shifting fitness landscapes in response to altered environments

The role of adaptation in molecular evolution has been contentious for decades. Here, we shed light on the adaptive potential in Saccharomyces cerevisiae by presenting systematic fitness measurements for all possible point mutations in a region of Hsp90 under four environmental conditions. Under elevated salinity, we observe numerous beneficial mutations with growth advantages up to 7% relative to the wild type. All of these beneficial mutations were observed to be associated with high costs of adaptation. We thus demonstrate that an essential protein can harbor adaptive potential upon an environmental challenge, and report a remarkable fit of the data to a version of Fisher's geometric model that focuses on the fitness trade-offs between mutations in different environments.

Predicting complex phenotype-genotype interactions to enable yeast engineering: Saccharomyces cerevisiae as a model organism and a cell factory

There is an increasing use of systems biology approaches in both "red" and "white" biotechnology in order to enable medical, medicinal, and industrial applications. The intricate links between genotype and phenotype may be explained through the use of the tools developed in systems biology, synthetic biology, and evolutionary engineering. Biomedical and biotechnological research are among the fields that could benefit most from the elucidation of this complex relationship. Researchers have studied fitness extensively to explain the phenotypic impacts of genetic variations. This elaborate network of dependencies and relationships so revealed are further complicated by the influence of environmental effects that present major challenges to our achieving an understanding of the cellular mechanisms leading to healthy or diseased phenotypes or optimized production yields. An improved comprehension of complex genotype-phenotype interactions and their accurate prediction should enable us to more effectively engineer yeast as a cell factory and to use it as a living model of human or pathogen cells in intelligent screens for new drugs. This review presents different methods and approaches undertaken toward improving our understanding and prediction of the growth phenotype of the yeast Saccharomyces cerevisiae as both a model and a production organism.