High-throughput screens have begun to reveal the protein interaction network that underpins most cellular functions in the yeast Saccharomyces cerevisiae. How the organization of this network affects the evolution of the proteins that compose it is a fundamental question in molecular evolution. We show that the connectivity of well-conserved proteins in the network is negatively correlated with their rate of evolution. Proteins with more interactors evolve more slowly not because they are more important to the organism, but because a greater proportion of the protein is directly involved in its function. At sites important for interaction between proteins, evolutionary changes may occur largely by coevolution, in which substitutions in one protein result in selection pressure for reciprocal changes in interacting partners. We confirm one predicted outcome of this process-namely, that interacting proteins evolve at similar rates.A protein's rate of evolution is thought to depend both on its dispensability to the organism and on the proportion of potential amino acid changes that are compatible with proper protein function (1). We recently analyzed functional genomic data (2) in conjunction with genomic comparisons (3) to confirm and further characterize the relation between protein dispensability and evolutionary rate (4). Here we apply a similar approach to investigate how protein function constrains evolution. Early studies of the structure and function of individual proteins suggested that, because molecular interactions require precisely specified structures, they impose constraints on sequence evolution (5, 6). Recent advances in the rapid detection of protein-protein interactions (7-9), as well as in the sequencing of complete genomes, allow us to expand the scale on which the evolutionary effects of molecular interactions are investigated and shift from a focus on individual proteins to a broad survey of the proteome and characterization of the general relation between protein interaction and evolution.We compiled a list of 3541 interactions between 2445 different yeast proteins (10). To estimate the evolutionary rates of these proteins, we compared putatively orthologous sequences between Saccharomyces cerevisiae and the nematode Caenorhabditis elegans (11). A subclass of putative orthologs, which we called "well-conserved orthologs," exhibited Ͼ50% amino acid identity over aligned regions; 1531 sequence pairs met our criteria for putative orthologs, and 309 of these were in the well-conserved class. For each pair of orthologs, we estimated the evolutionary distance (K) that separates the two sequences, where K is defined as the number of substitutions per amino acid site that have taken place since the fungi-animal split (12). There were 164 yeast proteins for which we had both an estimate of the number of interactors and a well-conserved ortholog in the nematode. Among these proteins, there is a significant negative correlation between each protein's number of interactors I and protein evolutionary rate...

Deleting a gene in an organism often has little phenotypic effect, owing to two mechanisms of compensation. The first is the existence of duplicate genes: that is, the loss of function in one copy can be compensated by the other copy or copies. The second mechanism of compensation stems from alternative metabolic pathways, regulatory networks, and so on. The relative importance of the two mechanisms has not been investigated except for a limited study, which suggested that the role of duplicate genes in compensation is negligible. The availability of fitness data for a nearly complete set of single-gene-deletion mutants of the Saccharomyces cerevisiae genome has enabled us to carry out a genome-wide evaluation of the role of duplicate genes in genetic robustness against null mutations. Here we show that there is a significantly higher probability of functional compensation for a duplicate gene than for a singleton, a high correlation between the frequency of compensation and the sequence similarity of two duplicates, and a higher probability of a severe fitness effect when the duplicate copy that is more highly expressed is deleted. We estimate that in S. cerevisiae at least a quarter of those gene deletions that have no phenotype are compensated by duplicate genes.

High mutation rate in mammalian mitochondrial DNA generates a highly divergent pool of alleles even within species that have dispersed and expanded in size recently. Phylogenetic analysis of 277 human mitochondrial genomes revealed a significant (P , 0.01) excess of rRNA and nonsynonymous base substitutions among hotspots of recurrent mutation. Most hotspots involved transitions from guanine to adenine that, with thymine-to-cytosine transitions, illustrate the asymmetric bias in codon usage at synonymous sites on the heavy-strand DNA. The mitochondrion-encoded tRNA Thr varied significantly more than any other tRNA gene. Threonine and valine codons were involved in 259 of the 414 amino acid replacements observed. The ratio of nonsynonymous changes from and to threonine and valine differed significantly (P ¼ 0.003) between populations with neutral (22/58) and populations with significantly negative Tajima's D values (70/76), independent of their geographic location. In contrast to a recent suggestion that the excess of nonsilent mutations is characteristic of Arctic populations, implying their role in cold adaptation, we demonstrate that the surplus of nonsynonymous mutations is a general feature of the young branches of the phylogenetic tree, affecting also those that are found only in Africa. We introduce a new calibration method of the mutation rate of synonymous transitions to estimate the coalescent times of mtDNA haplogroups.

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.

Wolfram syndrome is an autosomal recessive disorder characterized by juvenile diabetes mellitus, diabetes insipidus, optic atrophy and a number of neurological symptoms including deafness, ataxia and peripheral neuropathy. Mitochondrial DNA deletions have been described in a few patients and a locus has been mapped to 4p16 by linkage analysis. Susceptibility to psychiatric illness is reported to be high in affected individuals and increased in heterozygous carriers in Wolfram syndrome families. We screened four candidate genes in a refined critical linkage interval covered by an unfinished genomic sequence of 600 kb. One of these genes, subsequently named wolframin, codes for a predicted transmembrane protein which was expressed in various tissues, including brain and pancreas, and carried loss-of-function mutations in both alleles in Wolfram syndrome patients.

X-linked congenital stationary night blindness (XLCSNB) is characterized by impaired scotopic vision with associated ocular symptoms such as myopia, hyperopia, nystagmus and reduced visual acuity. Genetic mapping in families with XLCSNB revealed two different loci on the proximal short arm of the X chromosome. These two genetic subtypes can be distinguished on the basis of electroretinogram (ERG) responses and psychophysical testing as a complete (CSNB1) and an incomplete (CSNB2) form. The CSNB1 locus has been mapped to a 5-cM linkage interval in Xp11.4 (refs 2,5-7). Here we construct and analyse a contig between the markers DXS993 and DXS228, leading to the identification of a new gene mutated in CSNB1 patients. It is partially deleted in 3 families and mutation analysis in a further 21 families detected another 13 different mutations. This gene, designated NYX, encodes a protein of 481 amino acids (nyctalopin) and is expressed at low levels in tissues including retina, brain, testis and muscle. The predicted polypeptide is a glycosylphosphatidylinositol (GPI)-anchored extracellular protein with 11 typical and 2 cysteine-rich, leucine-rich repeats (LRRs). This motif is important for protein-protein interactions and members of the LRR superfamily are involved in cell adhesion and axon guidance. Future functional analysis of nyctalopin might therefore give insight into the fine-regulation of cell-cell contacts in the retina.

In this study yeast mitochondria were used as a model system to apply, evaluate, and integrate different genomic approaches to define the proteins of an organelle. Liquid chromatography mass spectrometry applied to purified mitochondria identified 546 proteins. By expression analysis and comparison to other proteome studies, we demonstrate that the proteomic approach identifies primarily highly abundant proteins. By expanding our evaluation to other types of genomic approaches, including systematic deletion phenotype screening, expression profiling, subcellular localization studies, protein interaction analyses, and computational predictions, we show that an integration of approaches moves beyond the limitations of any single approach. We report the success of each approach by benchmarking it against a reference set of known mitochondrial proteins, and predict approximately 700 proteins associated with the mitochondrial organelle from the integration of 22 datasets. We show that a combination of complementary approaches like deletion phenotype screening and mass spectrometry can identify over 75% of the known mitochondrial proteome. These findings have implications for choosing optimal genome-wide approaches for the study of other cellular systems, including organelles and pathways in various species. Furthermore, our systematic identification of genes involved in mitochondrial function and biogenesis in yeast expands the candidate genes available for mapping Mendelian and complex mitochondrial disorders in humans.

The MitoP2 database () integrates information on mitochondrial proteins, their molecular functions and associated diseases. The central database features are manually annotated reference proteins localized or functionally associated with mitochondria supplied for yeast, human and mouse. MitoP2 enables (i) the identification of putative orthologous proteins between these species to study evolutionarily conserved functions and pathways; (ii) the integration of data from systematic genome-wide studies such as proteomics and deletion phenotype screening; (iii) the prediction of novel mitochondrial proteins using data integration and the assignment of evidence scores; and (iv) systematic searches that aim to find the genes that underlie common and rare mitochondrial diseases. The data and analysis files are referenced to data sources in PubMed and other online databases and can be easily downloaded. MitoP2 users can explore the relationship between mitochondrial dysfunctions and disease and utilize this information to conduct systems biology approaches on mitochondria.