An integrated approach to characterize genetic interaction networks in yeast metabolism

Although experimental and theoretical efforts have been applied to globally map genetic interactions, we still do not understand how gene-gene interactions arise from the operation of biomolecular networks. To bridge the gap between empirical and computational studies, we i, quantitatively measured genetic interactions between ～185,000 metabolic gene pairs in Saccharomyces cerevisiae, ii, superposed the data on a detailed systems biology model of metabolism and iii, introduced a machine-learning method to reconcile empirical interaction data with model predictions. We systematically investigated the relative impacts of functional modularity and metabolic flux coupling on the distribution of negative and positive genetic interactions. We also provide a mechanistic explanation for the link between the degree of genetic interaction, pleiotropy and gene dispensability. Last, we show the feasibility of automated metabolic model refinement by correcting misannotations in NAD biosynthesis and confirming them by in vivo experiments.     

Distributions of epistasis in microbes fit predictions from a fitness landscape model

How do the fitness effects of several mutations combine? Despite its simplicity, this question is central to the understanding of multilocus evolution. Epistasis (the interaction between alleles at different loci), especially epistasis for fitness traits such as reproduction and survival, influences evolutionary predictions "almost whenever multilocus genetics matters". Yet very few models have sought to predict epistasis, and none has been empirically tested. Here we show that the distribution of epistasis can be predicted from the distribution of single mutation effects, based on a simple fitness landscape model. We show that this prediction closely matches the empirical measures of epistasis that have been obtained for Escherichia coli and the RNA virus vesicular stomatitis virus. Our results suggest that a simple fitness landscape model may be sufficient to quantitatively capture the complex nature of gene interactions. This model may offer a simple and widely applicable alternative to complex metabolic network models, in particular for making evolutionary predictions.     

Pervasive joint influence of epistasis and plasticity on mutational effects in Escherichia coli

The effects of mutations on phenotype and fitness may depend on the environment (phenotypic plasticity), other mutations (genetic epistasis) or both. Here we examine the fitness effects of 18 random insertion mutations in E. coli in two resource environments and five genetic backgrounds. We tested each mutation for plasticity and epistasis by comparing its fitness effects across these ecological and genetic contexts. Some mutations had no measurable effect in any of these contexts. None of the mutations had effects on phenotypic plasticity that were independent of genetic background. However, half the mutations had epistatic interactions such that their effects differed among genetic backgrounds, usually in an environment-dependent manner. Also, the pattern of mutational effects across backgrounds indicated that epistasis had been shaped primarily by unique events in the evolutionary history of a population rather than by repeatable events associated with shared environmental history.     

Negative epistasis between the malaria-protective effects of alpha+-thalassemia and the sickle cell trait

The hemoglobinopathies, disorders of hemoglobin structure and production, protect against death from malaria. In sub-Saharan Africa, two such conditions occur at particularly high frequencies: presence of the structural variant hemoglobin S and alpha(+)-thalassemia, a condition characterized by reduced production of the normal alpha-globin component of hemoglobin. Individually, each is protective against severe Plasmodium falciparum malaria, but little is known about their malaria-protective effects when inherited in combination. We investigated this question by studying a population on the coast of Kenya and found that the protection afforded by each condition inherited alone was lost when the two conditions were inherited together, to such a degree that the incidence of both uncomplicated and severe P. falciparum malaria was close to baseline in children heterozygous with respect to the mutation underlying the hemoglobin S variant and homozygous with respect to the mutation underlying alpha(+)-thalassemia. Negative epistasis could explain the failure of alpha(+)-thalassemia to reach fixation in any population in sub-Saharan Africa.     

Construction of a mouse yeast artificial chromosome library in a recombination-deficient strain of yeast.

We have constructed a new generation yeast artificial chromosome (YAC) library from female C57BL/10 mice in a recombination-deficient strain of Saccharomyces cerevisiae carrying a mutation in the RAD52 gene. The YAC library contains 41,568 clones with an average insert size of 240 kilobases, representing a greater than threefold coverage of the mouse genome. Currently, the library can be screened by polymerase chain reaction and we have isolated positive clones at a number of loci in the mouse genome. This rad52 library should enable a long-term assessment of the effect of one of the yeast recombination pathway genes on both, genome-wide YAC clone stability and the frequency of chimaeric clones.     

Robustness against mutations in genetic networks of yeast

There are two principal mechanisms that are responsible for the ability of an organism's physiological and developmental processes to compensate for mutations. In the first, genes have overlapping functions, and loss-of-function mutations in one gene will have little phenotypic effect if there are one or more additional genes with similar functions. The second mechanism has its origin in interactions between genes with unrelated functions, and has been documented in metabolic and regulatory gene networks. Here I analyse, on a genome-wide scale, which of these mechanisms of robustness against mutations is more prevalent. I used functional genomics data from the yeast Saccharomyces cerevisiae to test hypotheses related to the following: if gene duplications are mostly responsible for robustness, then a correlation is expected between the similarity of two duplicated genes and the effect of mutations in one of these genes. My results demonstrate that interactions among unrelated genes are the major cause of robustness against mutations. This type of robustness is probably an evolved response of genetic networks to stabilizing selection.     

Systematic screen for human disease genes in yeast

High similarity between yeast and human mitochondria allows functional genomic study of Saccharomyces cerevisiae to be used to identify human genes involved in disease. So far, 102 heritable disorders have been attributed to defects in a quarter of the known nuclear-encoded mitochondrial proteins in humans. Many mitochondrial diseases remain unexplained, however, in part because only 40-60% of the presumed 700-1,000 proteins involved in mitochondrial function and biogenesis have been identified. Here we apply a systematic functional screen using the pre-existing whole-genome pool of yeast deletion mutants to identify mitochondrial proteins. Three million measurements of strain fitness identified 466 genes whose deletions impaired mitochondrial respiration, of which 265 were new. Our approach gave higher selection than other systematic approaches, including fivefold greater selection than gene expression analysis. To apply these advantages to human disorders involving mitochondria, human orthologs were identified and linked to heritable diseases using genomic map positions.     

Population genomic analysis of outcrossing and recombination in yeast

The budding yeast Saccharomyces cerevisiae has been used by humans for millennia to make wine, beer and bread. More recently, it became a key model organism for studies of eukaryotic biology and for genomic analysis. However, relatively little is known about the natural lifestyle and population genetics of yeast. One major question is whether genetically diverse yeast strains mate and recombine in the wild. We developed a method to infer the evolutionary history of a species from genome sequences of multiple individuals and applied it to whole-genome sequence data from three strains of Saccharomyces cerevisiae and the sister species Saccharomyces paradoxus. We observed a pattern of sequence variation among yeast strains in which ancestral recombination events lead to a mosaic of segments with shared genealogy. Based on sequence divergence and the inferred median size of shared segments (approximately 2,000 bp), we estimated that although any two strains have undergone approximately 16 million cell divisions since their last common ancestor, only 314 outcrossing events have occurred during this time (roughly one every 50,000 divisions). Local correlations in polymorphism rates indicate that linkage disequilibrium in yeast should extend over kilobases. Our results provide the initial foundation for population studies of association between genotype and phenotype in S. cerevisiae.     

Predicting phenotypic variation in yeast from individual genome sequences

A central challenge in genetics is to predict phenotypic variation from individual genome sequences. Here we construct and evaluate phenotypic predictions for 19 strains of Saccharomyces cerevisiae. We use conservation-based methods to predict the impact of protein-coding variation within genes on protein function. We then rank strains using a prediction score that measures the total sum of function-altering changes in different sets of genes reported to influence over 100 phenotypes in genome-wide loss-of-function screens. We evaluate our predictions by comparing them with the observed growth rate and efficiency of 15 strains tested across 20 conditions in quantitative experiments. The median predictive performance, as measured by ROC AUC, was 0.76, and predictions were more accurate when the genes reported to influence a trait were highly connected in a functional gene network.     

Quantitative trait loci mapped to single-nucleotide resolution in yeast

Identifying the genetic variation underlying quantitative trait loci remains problematic. Consequently, our molecular understanding of genetically complex, quantitative traits is limited. To address this issue directly, we mapped three quantitative trait loci that control yeast sporulation efficiency to single-nucleotide resolution in a noncoding regulatory region (RME1) and to two missense mutations (TAO3 and MKT1). For each quantitative trait locus, the responsible polymorphism is rare among a diverse set of 13 yeast strains, suggestive of genetic heterogeneity in the control of yeast sporulation. Additionally, under optimal conditions, we reconstituted approximately 92% of the sporulation efficiency difference between the two genetically distinct parents by engineering three nucleotide changes in the appropriate yeast genome. Our results provide the highest resolution to date of the molecular basis of a quantitative trait, showing that the interaction of a few genetic variants can have a profound phenotypic effect.     

Resurrecting ancestral alcohol dehydrogenases from yeast

Modern yeast living in fleshy fruits rapidly convert sugars into bulk ethanol through pyruvate. Pyruvate loses carbon dioxide to produce acetaldehyde, which is reduced by alcohol dehydrogenase 1 (Adh1) to ethanol, which accumulates. Yeast later consumes the accumulated ethanol, exploiting Adh2, an Adh1 homolog differing by 24 (of 348) amino acids. As many microorganisms cannot grow in ethanol, accumulated ethanol may help yeast defend resources in the fruit. We report here the resurrection of the last common ancestor of Adh1 and Adh2, called Adh(A). The kinetic behavior of Adh(A) suggests that the ancestor was optimized to make (not consume) ethanol. This is consistent with the hypothesis that before the Adh1-Adh2 duplication, yeast did not accumulate ethanol for later consumption but rather used Adh(A) to recycle NADH generated in the glycolytic pathway. Silent nucleotide dating suggests that the Adh1-Adh2 duplication occurred near the time of duplication of several other proteins involved in the accumulation of ethanol, possibly in the Cretaceous age when fleshy fruits arose. These results help to connect the chemical behavior of these enzymes through systems analysis to a time of global ecosystem change, a small but useful step towards a planetary systems biology.     

Epistatic buffering of fitness loss in yeast double deletion strains

Interactions between deleterious mutations have been insufficiently studied, despite the fact that their strength and direction are critical for understanding the evolution of genetic recombination and the buildup of mutational load in populations. We compiled a list of 758 yeast gene deletions causing growth defects (from the Munich Information Center for Protein Sequences database and ref. 7). Using BY4741 and BY4742 single-deletion strains, we carried out 639 random crosses and assayed growth curves of the resulting progeny. We show that the maximum growth rate averaged over strains lacking deletions and those with double deletions is higher than that of strains with single deletions, indicating a positive epistatic effect. This tendency is shared by genes belonging to a variety of functional classes. Based on our data and former theoretical work, we suggest that epistasis is likely to diminish the negative effects of mutations when the ability to produce biomass at high rates contributes significantly to fitness.