Brevia

Genotype to Phenotype: A Complex Problem

See allHide authors and affiliations

Science  23 Apr 2010:
Vol. 328, Issue 5977, pp. 469
DOI: 10.1126/science.1189015

Abstract

We generated a high-resolution whole-genome sequence and individually deleted 5100 genes in Σ1278b, a Saccharomyces cerevisiae strain closely related to reference strain S288c. Similar to the variation between human individuals, Σ1278b and S288c average 3.2 single-nucleotide polymorphisms per kilobase. A genome-wide comparison of deletion mutant phenotypes identified a subset of genes that were conditionally essential by strain, including 44 essential genes unique to Σ1278b and 13 unique to S288c. Genetic analysis indicates the conditional phenotype was most often governed by complex genetic interactions, depending on multiple background-specific modifiers. Our comprehensive analysis suggests that the presence of a complex set of modifiers will often underlie the phenotypic differences between individuals.

Rapid genome-sequencing methods coupled with whole-genome transcription profiling suggests that it may be possible to predict phenotype from a genotype. Human genetic association studies of common single-nucleotide polymorphisms (SNPs) explain only a fraction of phenotypic variation among individuals (1). This may be due to rare SNPs (1), structural (2) and epigenetic variants, or multiple alleles with additive effects or synergistic genetic interactions associated with complex combinations of genetic variation (3).

To address the genotype-to-phenotype problem, we developed a simple comparative model for the budding yeast Saccharomyces cerevisiae that enables a comprehensive assessment of the genetic mechanisms leading to different phenotypes for the same mutation in two different genetic backgrounds. The strain Σ1278b mates and forms viable meiotic progeny with the reference strain, S288c, and the divergence between the two strains is roughly equal to the divergence between the genomes of two humans (4, 5).

We sequenced and assembled the 12-Mb Σ1278b genome, annotating 6923 open reading frames (ORFs) and RNAs, of which 6848 have orthologs within S288c (5). The order of genes between the strains was the same (except in the highly variable subtelomeric regions), and the sequence of 46% of the Σ1278b ORFs was identical to those in S288c. Differences between the strains were largely due to small insertions and deletions or SNPs, with an average SNP density of 3.2 per kilobase.

We deleted ~5100 genes within Σ1278b (5) to systematically compare identical deletion mutants (6). In particular, we identified “conditional essentials,” those genes required uniquely for viability in either strain (fig. S3 and table S1). We scored colonies as dead or alive and surveyed all vital pathways for individual-specific genetic interactions. We expected such conditional essential genes to be rare because the genomes of Σ1278b and S288c are nearly identical.

Although 894 genes were essential in both S288c and Σ1278b, 44 genes were essential only in Σ1278b and 13 genes were essential only in S288c (Fig. 1A). The conditional subsets included genes of various functions; however, the Σ1278b subset was enriched for genes involved in mRNA metabolic process, whereas the S288c set was enriched for genes annotated to SRP-dependent cotranslational targeting. These biological biases suggest that these phenotypes result from genetic interactions associated with an individual genotype.

Fig. 1

(A) Most S288c essential genes are also essential in Σ1278b (94%); however, ~5% are essential only in the Σ1278b genetic background, whereas ~1% are essential only in S288c. (B) Conditional essential genes in Σ1278b are the consequence of complex genetics. χ2 tests indicated the number of modifiers associated with conditional essentiality (5).

Hybrid strain crosses and tetrad analysis focusing on 18 mutants that were lethal in Σ1278b with wild-type levels of fitness when deleted in S288c were used to investigate conditional essentiality. We mated viable haploid S288c deletion mutants to wild-type Σ1278b and analyzed the hybrid diploid progeny by tetrad analysis. The number of viable meiotic progeny carrying the deletion allele is related to the number of unlinked background-specific modifiers (5) that contribute to the genetic interaction. In all 18 cases, the conditional phenotype was associated with numerous modifier genes that differ between strains. The simplest cases, SKI7 and BEM1, are likely due to a genetic interaction with at least two or more modifiers, but all other cases were more complex (Fig. 1B). Thus, our analysis showed that conditional essentiality is almost always a consequence of complex genetic interactions involving multiple modifiers associated with strain-specific genetic variation rather than classic digenic synthetic lethality (6, 7).

Our genome-wide survey of conditionally essential genes demonstrates that in most cases a complex set of background-specific modifiers influence a mutation whose phenotype differs between individuals. These results raise the possibility that similar complex modifiers may largely explain the difficulty in identifying the genetic basis for individual phenotypes. The potential for genetic interactions to control individual phenotypes becomes even more important if different combinations of alleles can lead to the same physiological state. The ability to identify these conditional essential phenotypes in yeast provides a framework to unravel the fundamental principles of genetic networks resulting from natural variation, including those that underlie human disease.

Supporting Online Material

www.sciencemag.org/cgi/content/full/328/5977/469/DC1

Materials and Methods

Figs. S1 and S2

Tables S1 to S3

References

  • * Present address: Molecular, Cellular, and Developmental Biology, 347 UCB, University of Colorado, Boulder, CO 80309, USA.

  • These authors contributed equally to this work.

References and Notes

  1. Materials and methods are available as supporting information on Science Online.
  2. Supported by NIH grants DK076284 (R.D.D.), GM035010 (G.R.F.), GM069676 (D.K.G.), P01 NS055923 (D.K.G.), ACS PF-09-072-01-MBC (D.A.B.); NSF grant DEB-0614588; a Senior Scholar Award in Aging from the Ellison Medical Foundation (P.C.P.); the Natural Sciences and Engineering Research Council (NSERC) of Canada; and Howard Hughes Medical Institute (HHMI) (C.B.). Sequence and assembly data has been deposited at DNA Data Bank of Japan/European Molecular Biology Laboratory/GenBank, project accession ACVY00000000.
View Abstract

Navigate This Article