Genomic exploration of the hemiascomycetous yeasts: 4. The genome of Saccharomyces cerevisiae revisited

Expression of virus-like particles (VLPs) of foot-and-mouth disease virus (FMDV) using Saccharomyces cerevisiae

We engineered Saccharomyces cerevisiae to express structural proteins of foot-and-mouth disease virus (FMDV) and produce virus-like particles (VLPs). The gene, which encodes four structural capsid proteins (VP0 (VP4 and VP2), VP3, and VP1), followed by a translational "ribosomal skipping" sequence consisting of 2A and protease 3C, was codon-optimized and chemically synthesized. The cloned gene was used to transform S. cerevisiae 2805 strain. Western blot analysis revealed that the polyprotein consisting of VP0, VP3, and VP1 was processed into the discrete capsid proteins. Western blot analysis of 3C confirmed the presence of discrete 3C protein, suggesting that the 2A sequence functioned as a "ribosomal skipping" signal in the yeast for an internal re-initiation of 3C translation from a monocistronic transcript, thereby indicating polyprotein processing by the discrete 3C protease. Moreover, a band corresponding to only VP2, which was known to be non-enzymatically processed from VP0 to both VP4 and VP2 during viral assembly, further validated the assembly of processed capsid proteins into VLPs. Electron microscopy showed the presence of the characteristic icosahedral VLPs. Our results clearly demonstrate that S. cerevisiae processes the viral structural polyprotein using a viral 3C protease and the resulting viral capsid subunits are assembled into virion particles. KEY POINTS: • Ribosomal skipping by self-cleaving FMDV peptide in S. cerevisiae. • Proteolytic processing of a structural polyprotein from a monocistronic transcript. • Assembly of the processed viral capsid proteins into a virus-like particle.

Metabolic engineering of Saccharomyces cerevisiae for the biosynthesis of a fungal pigment from the phytopathogenic fungus Cladosporium phlei

BACKGROUND

Cladosporium phlei is a phytopathogenic fungus that produces a pigment called phleichrome. This fungal perylenequinone plays an important role in the production of a photosensitizer that is a necessary component of photodynamic therapy. We applied synthetic biology to produce phleichrome using Saccharomyces cerevisiae.

RESULTS

The gene Cppks1, which encodes a non-reducing polyketide synthase (NR-PKS) responsible for the biosynthesis of phleichrome in C. phlei, was cloned into a yeast episomal vector and used to transform S. cerevisiae. In addition, a gene encoding a phosphopantetheinyl transferase (PPTase) of Aspergillus nidulans was cloned into a yeast integrative vector and also introduced into S. cerevisiae for the enzymatic activation of the protein product of Cppks1. Co-transformed yeasts were screened on a leucine/uracil-deficient selective medium and the presence of both integrative as well as episomal recombinant plasmids in the yeast were confirmed by colony PCR. The episomal vector for Cppks1 expression was so dramatically unstable during cultivation that most cells lost their episomal vector rapidly in nonselective media. This loss was also observed to a less degree in selective media. This data strongly suggests that the presence of the Cppks1 gene exerts a significant detrimental effect on the growth of transformed yeast cells and that selection pressure is required to maintain the Cppks1-expressing vector. The co-transformants on the selective medium showed the distinctive changes in pigmentation after a period of prolonged cultivation at 20 °C and 25 °C, but not at 30 °C. Furthermore, thin layer chromatography (TLC) revealed the presence of a spot corresponding with the purified phleichrome in the extract from the cells of the co-transformants. Liquid chromatography (LC/MS/MS) verified that the newly expressed pigment was indeed phleichrome.

CONCLUSION

Our results indicate that metabolic engineering by multiple gene expression is possible and capable of producing fungal pigment phleichrome in S. cerevisiae. This result adds to our understanding of the characteristics of fungal PKS genes, which exhibit complex structures and diverse biological activities.

A vast evolutionarily transient translatome contributes to phenotype and fitness

Translation is the process by which ribosomes synthesize proteins. Ribosome profiling recently revealed that many short sequences previously thought to be noncoding are pervasively translated. To identify protein-coding genes in this noncanonical translatome, we combine an integrative framework for extremely sensitive ribosome profiling analysis, iRibo, with high-powered selection inferences tailored for short sequences. We construct a reference translatome for Saccharomyces cerevisiae comprising 5,400 canonical and almost 19,000 noncanonical translated elements. Only 14 noncanonical elements were evolving under detectable purifying selection. A representative subset of translated elements lacking signatures of selection demonstrated involvement in processes including DNA repair, stress response, and post-transcriptional regulation. Our results suggest that most translated elements are not conserved protein-coding genes and contribute to genotype-phenotype relationships through fast-evolving molecular mechanisms.

Deciphering the Origin, Evolution, and Physiological Function of the Subtelomeric Aryl-Alcohol Dehydrogenase Gene Family in the Yeast Saccharomyces cerevisiae

Homology searches indicate that Saccharomyces cerevisiae strain BY4741 contains seven redundant genes that encode putative aryl-alcohol dehydrogenases (AAD). Yeast AAD genes are located in subtelomeric regions of different chromosomes, and their functional role(s) remain enigmatic. Here, we show that two of these genes, AAD4 and AAD14, encode functional enzymes that reduce aliphatic and aryl-aldehydes concomitant with the oxidation of cofactor NADPH, and that Aad4p and Aad14p exhibit different substrate preference patterns. Other yeast AAD genes are undergoing pseudogenization. The 5′ sequence of AAD15 has been deleted from the genome. Repair of an AAD3 missense mutation at the catalytically essential Tyr⁷³ residue did not result in a functional enzyme. However, ancestral-state reconstruction by fusing Aad6 with Aad16 and by N-terminal repair of Aad10 restores NADPH-dependent aryl-alcohol dehydrogenase activities. Phylogenetic analysis indicates that AAD genes are narrowly distributed in wood-saprophyte fungi and in yeast that occupy lignocellulosic niches. Because yeast AAD genes exhibit activity on veratraldehyde, cinnamaldehyde, and vanillin, they could serve to detoxify aryl-aldehydes released during lignin degradation. However, none of these compounds induce yeast AAD gene expression, and Aad activities do not relieve aryl-aldehyde growth inhibition. Our data suggest an ancestral role for AAD genes in lignin degradation that is degenerating as a result of yeast's domestication and use in brewing, baking, and other industrial applications.

IMPORTANCE Functional characterization of hypothetical genes remains one of the chief tasks of the postgenomic era. Although the first Saccharomyces cerevisiae genome sequence was published over 20 years ago, 22% of its estimated 6,603 open reading frames (ORFs) remain unverified. One outstanding example of this category of genes is the enigmatic seven-member AAD family. Here, we demonstrate that proteins encoded by two members of this family exhibit aliphatic and aryl-aldehyde reductase activity, and further that such activity can be recovered from pseudogenized AAD genes via ancestral-state reconstruction. The phylogeny of yeast AAD genes suggests that these proteins may have played an important ancestral role in detoxifying aromatic aldehydes in ligninolytic fungi. However, in yeast adapted to niches rich in sugars, AAD genes become subject to mutational erosion. Our findings shed new light on the selective pressures and molecular mechanisms by which genes undergo pseudogenization.

Small proteins: untapped area of potential biological importance

Polypeptides containing ≤100 amino acid residues (AAs) are generally considered to be small proteins (SPs). Many studies have shown that some SPs are involved in important biological processes, including cell signaling, metabolism, and growth. SP generally has a simple domain and has an advantage to be used as model system to overcome folding speed limits in protein folding simulation and drug design. But SPs were once thought to be trivial molecules in biological processes compared to large proteins. Because of the constraints of experimental methods and bioinformatics analysis, many genome projects have used a length threshold of 100 amino acid residues to minimize erroneous predictions and SPs are relatively under-represented in earlier studies. The general protein discovery methods have potential problems to predict and validate SPs, and very few effective tools and algorithms were developed specially for SPs identification. In this review, we mainly consider the diverse strategies applied to SPs prediction and discuss the challenge for differentiate SP coding genes from artifacts. We also summarize current large-scale discovery of SPs in species at the genome level. In addition, we present an overview of SPs with regard to biological significance, structural application, and evolution characterization in an effort to gain insight into the significance of SPs.

Re-annotation of protein-coding genes in the genome of saccharomyces cerevisiae based on support vector machines

The annotation of the well-studied organism, Saccharomyces cerevisiae, has been improving over the past decade while there are unresolved debates over the amount of biologically significant open reading frames (ORFs) in yeast genome. We revisited the total count of protein-coding genes in S. cerevisiae S288c genome using a theoretical approach by combining the Support Vector Machine (SVM) method with six widely used measurements of sequence statistical features. The accuracy of our method is over 99.5% in 10-fold cross-validation. Based on the annotation data in Saccharomyces Genome Database (SGD), we studied the coding capacity of all 1744 ORFs which lack experimental results and suggested that the overall number of chromosomal ORFs encoding proteins in yeast should be 6091 by removing 488 spurious ORFs. The importance of the present work lies in at least two aspects. First, cross-validation and retrospective examination showed the fidelity of our method in recognizing ORFs that likely encode proteins. Second, we have provided a web service that can be accessed at http://cobi.uestc.edu.cn/services/yeast/, which enables the prediction of protein-coding ORFs of the genus Saccharomyces with a high accuracy.

Hundreds of putatively functional small open reading frames in Drosophila

Background

The relationship between DNA sequence and encoded information is still an unsolved puzzle. The number of protein-coding genes in higher eukaryotes identified by genome projects is lower than was expected, while a considerable amount of putatively non-coding transcription has been detected. Functional small open reading frames (smORFs) are known to exist in several organisms. However, coding sequence detection methods are biased against detecting such very short open reading frames. Thus, a substantial number of non-canonical coding regions encoding short peptides might await characterization.

Results

Using bio-informatics methods, we have searched for smORFs of less than 100 amino acids in the putatively non-coding euchromatic DNA of Drosophila melanogaster, and initially identified nearly 600,000 of them. We have studied the pattern of conservation of these smORFs as coding entities between D. melanogaster and Drosophila pseudoobscura, their presence in syntenic and in transcribed regions of the genome, and their ratio of conservative versus non-conservative nucleotide changes. For negative controls, we compared the results with those obtained using random short sequences, while a positive control was provided by smORFs validated by proteomics data.

Conclusions

The combination of these analyses led us to postulate the existence of at least 401 functional smORFs in Drosophila, with the possibility that as many as 4,561 such functional smORFs may exist.

Characterization of a new multigene family encoding isomaltases in the yeast Saccharomyces cerevisiae, the IMA family

It has been known for a long time that the yeast Saccharomyces cerevisiae can assimilate alpha-methylglucopyranoside and isomaltose. We here report the identification of 5 genes (YGR287c, YIL172c, YJL216c, YJL221c and YOL157c), which, similar to the SUCx, MALx, or HXTx multigene families, are located in the subtelomeric regions of different chromosomes. They share high nucleotide sequence identities between themselves (66-100%) and with the MALx2 genes (63-74%). Comparison of their amino acid sequences underlined a substitution of threonine by valine in region II, one of the four highly conserved regions of the alpha-glucosidase family. This change was previously shown to be sufficient to discriminate alpha-1,4- to alpha-1,6-glucosidase activity in YGR287c (Yamamoto, K., Nakayama, A., Yamamoto, Y., and Tabata, S. (2004) Eur. J. Biochem. 271, 3414-3420). We showed that each of these five genes encodes a protein with alpha-glucosidase activity on isomaltose, and we therefore renamed these genes IMA1 to IMA5 for IsoMAltase. Our results also illustrated that sequence polymorphisms among this family led to interesting variability of gene expression patterns and of catalytic efficiencies on different substrates, which altogether should account for the absence of functional redundancy for growth on isomaltose. Indeed, deletion studies revealed that IMA1/YGR287c encodes the major isomaltase and that growth on isomaltose required the presence of AGT1, which encodes an alpha-glucoside transporter. Expressions of IMA1 and IMA5/YJL216c were strongly induced by maltose, isomaltose, and alpha-methylglucopyranoside, in accordance with their regulation by the Malx3p-transcription system. The physiological relevance of this IMAx multigene family in S. cerevisiae is discussed.

ORF organization and gene recognition in the yeast genome

Some rules on gene recognition and ORF organization in the Saccharomyces cerevisiae genome are demonstrated by statistical analyses of sequence data. This study includes: (a) The random frame rule-that the six reading frames W1, W2, W3, C1, C2 and C3 in the double-stranded genome are randomly occupied by ORFs (related phenomena on ORF overlapping are also discussed). (b) The inhomogeneity rule-coding and non-coding ORFs differ in inhomogeneity of base composition in the three codon positions. By use of the inhomogeneity index (IHI), one can make a distinction between coding (IHI > 14) and non-coding (IHI < or = 14) ORFs at 95% accuracy. We find that 'spurious' ORFs (with IHI < or = 14) are distributed mainly in three classes of ORFs, namely, those with 'similarity to unknown proteins', those with 'no similarity', or 'questionable ORFs'. The total number of spurious ORFs (which are unlikely to be regarded as coding ORFs) is estimated to be 470. (c) The evaluation of ORF length distribution shows that below 200 amino acids the occurrence of ATG initiator ORFs is close to random.

Functional analyses of PAU genes in Saccharomyces cerevisiae

PAU genes constitute the largest gene family in Saccharomyces cerevisiae, with 24 members mostly located in the subtelomeric regions of chromosomes. Little information is available about PAU genes, other than expression data for some members. In this study, we systematically compared the sequences of all 24 members, examined the expression of PAU3, PAU5, DAN2, PAU17 and PAU20 in response to stresses, and investigated the stability of all Pau proteins. The chromosomal localization, synteny and sequence analyses revealed that PAU genes could have been amplified by segmental and retroposition duplication through mechanisms of chromosomal end translocation and Ty-associated recombination. The coding sequences diverged through nucleotide substitution and insertion/deletion of one to four codons, thus causing changes in amino acids, truncation or extension of Pau proteins. Pairwise comparison of non-coding regions revealed little homology in flanking sequences of some members. All 24 PAU promoters contain a TATA box, and 22 PAU promoters contain at least one copy of the anaerobic response element and the aerobic repression motif. Differential expression was observed among PAU3, PAU5, PAU17, PAU20 and DAN2 in response to stress, with PAU5 having the highest capacity to be induced by anaerobic conditions, low temperature and wine fermentations. Furthermore, Pau proteins with 124 aa were less stable than those with 120 or 122 aa. Our results indicate that duplicated PAU genes have been evolving, and the individual Pau proteins might possess specific roles for the adaptation of S. cerevisiae to certain environmental stresses.

Why are there still over 1000 uncharacterized yeast genes?

The yeast genetics community has embraced genomic biology, and there is a general understanding that obtaining a full encyclopedia of functions of the approximately 6000 genes is a worthwhile goal. The yeast literature comprises over 40,000 research papers, and the number of yeast researchers exceeds the number of genes. There are mutated and tagged alleles for virtually every gene, and hundreds of high-throughput data sets and computational analyses have been described. Why, then, are there >1000 genes still listed as uncharacterized on the Saccharomyces Genome Database, 10 years after sequencing the genome of this powerful model organism? Examination of the currently uncharacterized gene set suggests that while some are small or newly discovered, the vast majority were evident from the initial genome sequence. Most are present in multiple genomics data sets, which may provide clues to function. In addition, roughly half contain recognizable protein domains, and many of these suggest specific metabolic activities. Notably, the uncharacterized gene set is highly enriched for genes whose only homologs are in other fungi. Achieving a full catalog of yeast gene functions may require a greater focus on the life of yeast outside the laboratory.

Mutations in the SF1-U2AF59-U2AF23 Complex Cause Exon Skipping in Schizosaccharomyces pombe

To identify genes involved in the mechanism to ensure ordered 5' to 3' exon joining in constitutively spliced pre-mRNAs, we screened for mutants that cause exon skipping in the fission yeast Schizosaccharomyces pombe using a reporter plasmid, which contains the ura4+ gene with the nda3 intron 1-exon 2-intron 2 sequence. The reporter plasmid was designed to produce the functional ura4+ mRNA, when the central nda3 exon is skipped during the splicing reaction. We mutagenized cells harboring the plasmid by UV irradiation and isolated 34 ura+ mutants that grew on minimal medium. Of those, eight mutants were found to be temperature sensitive (ts) for growth. Complementation analyses revealed that the ts mutants belong to three distinct complementation groups named ods (ordered splicing) 1, 2, and 3. RT-PCR analyses showed that products of exon skipping were actually generated in the ods mutants. We cloned the genes responsible for the ods mutations, and found that ods1+, ods2+, and ods3+ encode splicing factors Prp2p/U2AF59, U2AF23, and SF1, respectively, which form a SF1-U2AF59-U2AF23 complex involved in recognition of the branch-point and 3' splice site sequences in a pre-mRNA. We also showed that mutations in the SF1-U2AF59-U2AF23 binding sequences in the reporter plasmid result in exon skipping in wild-type S. pombe cells. In addition, drugs that decrease the rate of transcription elongation were found to suppress the exon skipping in the ods mutants. These results suggest that co-transcriptional recognition of a nascent pre-mRNA by the SF1-U2AF59-U2AF23 complex is essential for ordered exon joining in constitutive splicing in S. pombe.

Saccharomyces cerevisiae S288C genome annotation: a working hypothesis

The S. cerevisiae genome is the most well-characterized eukaryotic genome and one of the simplest in terms of identifying open reading frames (ORFs), yet its primary annotation has been updated continually in the decade since its initial release in 1996 (Goffeau et al., 1996). The Saccharomyces Genome Database (SGD; www.yeastgenome.org) (Hirschman et al., 2006), the community-designated repository for this reference genome, strives to ensure that the S. cerevisiae annotation is as accurate and useful as possible. At SGD, the S. cerevisiae genome sequence and annotation are treated as a working hypothesis, which must be repeatedly tested and refined. In this paper, in celebration of the tenth anniversary of the completion of the S. cerevisiae genome sequence, we discuss the ways in which the S. cerevisiae sequence and annotation have changed, consider the multiple sources of experimental and comparative data on which these changes are based, and describe our methods for evaluating, incorporating and documenting these new data.

Functional genomics of genes with small open reading frames (sORFs) in S. cerevisiae

Genes with small open reading frames (sORFs; <100 amino acids) represent an untapped source of important biology. sORFs largely escaped analysis because they were difficult to predict computationally and less likely to be targeted by genetic screens. Thus, the substantial number of sORFs and their potential importance have only recently become clear. To investigate sORF function, we undertook the first functional studies of sORFs in any system, using the model eukaryote Saccharomyces cerevisiae. Based on independent experimental approaches and computational analyses, evidence exists for 299 sORFs in the S. cerevisiae genome, representing approximately 5% of the annotated ORFs. We determined that a similar percentage of sORFs are annotated in other eukaryotes, including humans, and 184 of the S. cerevisiae sORFs exhibit similarity with ORFs in other organisms. To investigate sORF function, we constructed a collection of gene-deletion mutants of 140 newly identified sORFs, each of which contains a strain-specific "molecular barcode," bringing the total number of sORF deletion strains to 247. Phenotypic analyses of the new gene-deletion strains identified 22 sORFs required for haploid growth, growth at high temperature, growth in the presence of a nonfermentable carbon source, or growth in the presence of DNA damage and replication-arrest agents. We provide a collection of sORF deletion strains that can be integrated into the existing deletion collection as a resource for the yeast community for elucidating gene function. Moreover, our analyses of the S. cerevisiae sORFs establish that sORFs are conserved across eukaryotes and have important biological functions.

A novel system of genetic transformation allows multiple integrations of a desired gene in Saccharomyces cerevisiae chromosomes

Increasing industrial competitiveness and productivity demand that recombinant yeast strains, used in many different processes, be constantly adapted and/or genetically improved to suit changing requirements. Among yeasts, Saccharomyces cerevisiae is the best-studied organism, and the most frequently employed yeast in industrial processes. In the present study, laboratory strains and industrial S. cerevisiae strains were stably transformed with a novel vector containing the glucoamylase cDNA of Aspergillus awamori flanked by delta-sequences (deltaGlucodelta), and lacking a positive selection marker. Co-transformation with known plasmids allowed selection by auxotrophic complementation of the leu2 mutation and/or geneticin resistance (G418). In all cases, several copies of the deltaGlucodelta vector were inserted into the genome of the yeast cell without selective pressure, showing 100% stability after 80 generations. Transformation frequency of the new vector was similar for S. cerevisiae laboratory strains and industrial wild-type S. cerevisiae strains. This novel genetic transformation system is versatile and suitable to introduce several stable copies of a desired expression cassette into the genome of different S. cerevisiae yeast strains.

Genome Snapshot: a new resource at the Saccharomyces Genome Database (SGD) presenting an overview of the Saccharomyces cerevisiae genome

Sequencing and annotation of the entire Saccharomyces cerevisiae genome has made it possible to gain a genome-wide perspective on yeast genes and gene products. To make this information available on an ongoing basis, the Saccharomyces Genome Database (SGD) () has created the Genome Snapshot (). The Genome Snapshot summarizes the current state of knowledge about the genes and chromosomal features of S.cerevisiae. The information is organized into two categories: (i) number of each type of chromosomal feature annotated in the genome and (ii) number and distribution of genes annotated to Gene Ontology terms. Detailed lists are accessible through SGD's Advanced Search tool (), and all the data presented on this page are available from the SGD ftp site ().

Yeast evolution and comparative genomics

Recent sequencing efforts and experiments have advanced our understanding of genome evolution in yeasts, particularly the Saccharomyces yeasts. The ancestral genome of the Saccharomyces sensu stricto complex has been subject to both whole-genome duplication, followed by massive sequence loss and divergence, and segmental duplication. In addition the subtelomeric regions are subject to further duplications and rearrangements via ectopic exchanges. Translocations and other gross chromosomal rearrangements that break down syntenic relationships occur; however, they do not appear to be a driving force of speciation. Analysis of single genomes has been fruitful for hypothesis generation such as the whole-genome duplication, but comparative genomics between close and more distant species has proven to be a powerful tool in testing these hypotheses as well as elucidating evolutionary processes acting on the genome. Future work on population genomics and experimental evolution will keep yeast at the forefront of studies in genome evolution.

Biochemical and genetic analysis of the yeast proteome with a movable ORF collection

Functional analysis of the proteome is an essential part of genomic research. To facilitate different proteomic approaches, a MORF (moveable ORF) library of 5854 yeast expression plasmids was constructed, each expressing a sequence-verified ORF as a C-terminal ORF fusion protein, under regulated control. Analysis of 5573 MORFs demonstrates that nearly all verified ORFs are expressed, suggests the authenticity of 48 ORFs characterized as dubious, and implicates specific processes including cytoskeletal organization and transcriptional control in growth inhibition caused by overexpression. Global analysis of glycosylated proteins identifies 109 new confirmed N-linked and 345 candidate glycoproteins, nearly doubling the known yeast glycome.

Comparative genomics of microsporidia

Microsporidia have been known for some time to possess among the smallest genomes of any eukaryote. There is now a completely sequenced microsporidian genome, as well as several other large-scale sequencing efforts, so the nature of these genomes is becoming apparent. This paper reviews some of the characteristics of microsporidian genomes in general, and some of the recent discoveries made through comparative genomic analyses. In general, microsporidian genomes are both reduced and compacted. Reduction takes place through gene loss, which is understandable in obligate intracellular parasites that rely on their host for many metabolites. Compaction is a more complex process, and is as yet not fully understood. It is clear from genomes surveyed thus far that the remaining genes are tightly packed and that there is little non-coding sequence, resulting in some extraordinary arrangements, including overlapping genes. Compaction also seems to affect certain aspects of genome evolution, like the frequency of rearrangements. The force behind this compaction is not known, and is especially interesting in light of the fact that surveys of genomes that are significantly different in size yield similar complements of protein-coding genes. There are some interesting exceptions, including catalase, photolyase and some mitochondrial proteins, but the rarity of these raises an interesting question as to what accounts for the significant differences seen in the genome sizes among microsporidia.

Fungal BLAST and Model Organism BLASTP Best Hits: new comparison resources at the Saccharomyces Genome Database (SGD)

The Saccharomyces Genome Database (SGD; http://www.yeastgenome.org/) is a scientific database of gene, protein and genomic information for the yeast Saccharomyces cerevisiae. SGD has recently developed two new resources that facilitate nucleotide and protein sequence comparisons between S.cerevisiae and other organisms. The Fungal BLAST tool provides directed searches against all fungal nucleotide and protein sequences available from GenBank, divided into categories according to organism, status of completeness and annotation, and source. The Model Organism BLASTP Best Hits resource displays, for each S.cerevisiae protein, the single most similar protein from several model organisms and presents links to the database pages of those proteins, facilitating access to curated information about potential orthologs of yeast proteins.

A highly conserved binding site in vesicle-associated membrane protein-associated protein (VAP) for the FFAT motif of lipid-binding proteins

A variety of lipid-binding proteins contain a recently described motif, designated FFAT (two phenylalanines in an acidic tract), which binds to vesicle-associated-membrane protein-associated protein (VAP). VAP is a conserved integral membrane protein of the endoplasmic reticulum that contains at its amino terminus a domain related to the major sperm protein of nematode worms. Here we have studied the FFAT-VAP interaction in Saccharomyces cerevisiae, where the VAP homologue Scs2 regulates phospholipid metabolism via an interaction with the FFAT motif of Opi1. By introducing mutations at random into Scs2, we found that mutations that abrogated binding to FFAT were clustered in the most highly conserved region. Using site-directed mutagenesis, we identified several critical residues, including two lysines widely separated in the primary sequence. By examining all other conserved basic residues, we identified a third residue that was moderately important for binding FFAT. Modeling VAP on the known structure of major sperm protein showed that the critical residues form a patch on a positively charged face of the protein. In vivo functional studies of SCS22, a second SCS2-like gene in S. cerevisiae, showed that SCS2 was the dominant gene in the regulation of Opi1, with a minor contribution from SCS22. We then established that reduction in the affinity of Scs2 mutants for FFAT correlated well with loss of function, indicating the importance of these residues for binding FFAT motifs. Finally, we found that human VAP-A could substitute for Scs2 but that it functioned poorly, suggesting that other factors modulate the binding of Scs2 to proteins with FFAT motifs.

Yos1p is a novel subunit of the Yip1p-Yif1p complex and is required for transport between the endoplasmic reticulum and the Golgi complex

Yeast Yip1p is a member of a conserved family of transmembrane proteins that interact with Rab GTPases. Previous studies also have indicated a role for Yip1p in the biogenesis of endoplasmic reticulum (ER)-derived COPII transport vesicles. In this report, we describe the identification and characterization of the uncharacterized open reading frame YER074W-A as a novel multicopy suppressor of the thermosensitive yip1-4 strain. We have termed this gene Yip One Suppressor 1 (YOS1). Yos1p is essential for growth and for function of the secretory pathway; depletion or inactivation of Yos1p blocks transport between the ER and the Golgi complex. YOS1 encodes an integral membrane protein of 87 amino acids that is conserved in eukaryotes. Yos1p localizes to ER and Golgi membranes and is efficiently packaged into ER-derived COPII transport vesicles. Yos1p associates with Yip1p and Yif1p, indicating Yos1p is a novel subunit of the Yip1p-Yif1p complex.

Charting gene regulatory networks: strategies, challenges and perspectives

One of the foremost challenges in the post-genomic era will be to chart the gene regulatory networks of cells, including aspects such as genome annotation, identification of cis-regulatory elements and transcription factors, information on protein-DNA and protein-protein interactions, and data mining and integration. Some of these broad sets of data have already been assembled for building networks of gene regulation. Even though these datasets are still far from comprehensive, and the approach faces many important and difficult challenges, some strategies have begun to make connections between disparate regulatory events and to foster new hypotheses. In this article we review several different genomics and proteomics technologies, and present bioinformatics methods for exploring these data in order to make novel discoveries.

Large-scale exploration of growth inhibition caused by overexpression of genomic fragments in Saccharomyces cerevisiae

We have screened the genome of Saccharomyces cerevisiae for fragments that confer a growth-retardation phenotype when overexpressed in a multicopy plasmid with a tetracycline-regulatable (Tet-off) promoter. We selected 714 such fragments with a mean size of 700 base-pairs out of around 84,000 clones tested. These include 493 in-frame open reading frame fragments corresponding to 454 distinct genes (of which 91 are of unknown function), and 162 out-of-frame, antisense and intergenic genomic fragments, representing the largest collection of toxic inserts published so far in yeast.

Methods in comparative genomics: genome correspondence, gene identification and regulatory motif discovery

In Kellis et al. (2003), we reported the genome sequences of S. paradoxus, S. mikatae, and S. bayanus and compared these three yeast species to their close relative, S. cerevisiae. Genomewide comparative analysis allowed the identification of functionally important sequences, both coding and noncoding. In this companion paper we describe the mathematical and algorithmic results underpinning the analysis of these genomes. (1) We present methods for the automatic determination of genome correspondence. The algorithms enabled the automatic identification of orthologs for more than 90% of genes and intergenic regions across the four species despite the large number of duplicated genes in the yeast genome. The remaining ambiguities in the gene correspondence revealed recent gene family expansions in regions of rapid genomic change. (2) We present methods for the identification of protein-coding genes based on their patterns of nucleotide conservation across related species. We observed the pressure to conserve the reading frame of functional proteins and developed a test for gene identification with high sensitivity and specificity. We used this test to revisit the genome of S. cerevisiae, reducing the overall gene count by 500 genes (10% of previously annotated genes) and refining the gene structure of hundreds of genes. (3) We present novel methods for the systematic de novo identification of regulatory motifs. The methods do not rely on previous knowledge of gene function and in that way differ from the current literature on computational motif discovery. Based on genomewide conservation patterns of known motifs, we developed three conservation criteria that we used to discover novel motifs. We used an enumeration approach to select strongly conserved motif cores, which we extended and collapsed into a small number of candidate regulatory motifs. These include most previously known regulatory motifs as well as several noteworthy novel motifs. The majority of discovered motifs are enriched in functionally related genes, allowing us to infer a candidate function for novel motifs. Our results demonstrate the power of comparative genomics to further our understanding of any species. Our methods are validated by the extensive experimental knowledge in yeast and will be invaluable in the study of complex genomes like that of the human.

Yeast genome sequencing: the power of comparative genomics

For decades, unicellular yeasts have been general models to help understand the eukaryotic cell and also our own biology. Recently, over a dozen yeast genomes have been sequenced, providing the basis to resolve several complex biological questions. Analysis of the novel sequence data has shown that the minimum number of genes from each species that need to be compared to produce a reliable phylogeny is about 20. Yeast has also become an attractive model to study speciation in eukaryotes, especially to understand molecular mechanisms behind the establishment of reproductive isolation. Comparison of closely related species helps in gene annotation and to answer how many genes there really are within the genomes. Analysis of non-coding regions among closely related species has provided an example of how to determine novel gene regulatory sequences, which were previously difficult to analyse because they are short and degenerate and occupy different positions. Comparative genomics helps to understand the origin of yeasts and points out crucial molecular events in yeast evolutionary history, such as whole-genome duplication and horizontal gene transfer(s). In addition, the accumulating sequence data provide the background to use more yeast species in model studies, to combat pathogens and for efficient manipulation of industrial strains.

Gene relics in the genome of the yeast Saccharomyces cerevisiae

There is increasing evidence that DNA duplication is a common and ongoing process that plays a major role in molecular evolution of genomes and that a large fraction of the duplicated gene copies becomes non-functional by accumulation of deleterious mutations. In order to describe this phenomenon, we systematically searched the 6404 intergenic regions (IRs) of the genome of Saccharomyces cerevisiae for traces of coding sequences presenting degenerated but still recognizable sequence similarity with active open reading frames (5823 annotated ORFs). We detected a total of 124 anciently coding regions, or "gene relics", showing similarity to a total of 149 distinct active ORFs. This set of relics shows a continuum of sequence degeneration from those whose sequence is slightly altered compared to the functional ORF (classically defined as pseudogenes), to those that contains so many deleterious mutations, as to reach the limit of recognition. Gene relics are more concentrated in the subtelomeric regions of the chromosomes, reflecting the high plasticity of these regions. The presence of relics also revealed ancestral duplication events of chromosomal segments that were previously undetected. Some of these segments are intermingled with the more easily recognizable ancestral blocks of duplication, indicating successive duplication events. We present a compilation of all the data available, leading to a total of 278 pseudogenes in the genome of S. cerevisiae.

Differential evolution of the Saccharomyces cerevisiae DUP240 paralogs and implication of recombination in phylogeny

Multigene families are observed in all genomes sequenced so far and are the reflection of key evolutionary mechanisms. The DUP240 family, identified in Saccharomyces cerevisiae strain S288C, is composed of 10 paralogs: seven are organized as two tandem repeats and three are solo ORFs. To investigate the evolution of the three solo paralogs, YAR023c, YCR007c and YHL044w, we performed a comparative analysis between 15 S.cerevisiae strains. These three ORFs are present in all strains and the conservation of synteny indicates that they are not frequently involved in chromosomal reshaping, in contrast to the DUP240 ORFs organized in tandem repeats. Our analysis of nucleotide and amino acid variations indicates that YAR023c and YHL044w fix mutations more easily than YCR007c, although they all belong to the same multigene family. This comparative analysis was also conducted with five arbitrarily chosen Ascomycetes-specific genes and five arbitrarily chosen common genes (genes that have a homolog in at least one non-Ascomycetes organism). Ascomycetes-specific genes appear to be diverging faster than common genes in the S.cerevisiae species, a situation that was previously described between different yeast species. Our results point to the strong contribution, during DNA sequence evolution, of allelic recombination besides nucleotide substitution.

A "polyORFomic" analysis of prokaryote genomes using disabled-homology filtering reveals conserved but undiscovered short ORFs

Prokaryote gene annotation is complicated by large numbers of short open reading frames (ORFs) that arise naturally from genetic code design. Historically, many hypothetical ORFs have been annotated as genes in microbes, usually with an arbitrary length threshold (e.g. greater than 100 codons). Given the use of such thresholds, what is the extent of genuine undiscovered short genes in the current sampling of prokaryote genomes? To assess rigorously the potential under-annotation of short ORFs with homology, we exhaustively compared the polyORFome--all possible ORFs in 64 prokaryotes (53 bacteria and 11 archaea) plus budding yeast--to itself and to all known proteins. The novelty of our analysis is that, firstly, sequence comparisons to/between both annotated and un-annotated ORFs are considered, and secondly a two-step disabled-homology filter is applied to set aside putative pseudogenes and spurious ORFs. We find that un-annotated homologous short ORFs (uhORFs) correspond to a small but non-negligible fraction of the annotated prokaryote proteomes (0.5-3.8%, depending on selection criteria). Moreover, the disabled-homology filter indicates that about a third of uhORFs correspond to putative pseudogenes or spurious ORFs. Our analysis shows that the use of annotation length thresholds is unnecessary, as there are manageable numbers of short ORF homologies conserved (without disablements) across microbial genomes. Data on uhORFs are available from http://pseudogene.org/polyo

Candida tropicalis expresses two mitochondrial 2-enoyl thioester reductases that are able to form both homodimers and heterodimers

Here we report on the cloning of a Candida tropicalis gene, ETR2, that is closely related to ETR1. Both genes encode enzymatically active 2-enoyl thioester reductases involved in mitochondrial synthesis of fatty acids (fatty acid synthesis type II) and respiratory competence. The 5'- and 3'-flanking (coding) regions of ETR2 and ETR1 are about 90% (97%) identical, indicating that the genes have evolved via gene duplication. The gene products differ in three amino acid residues: Ile67 (Val), Ala92 (Thr), and Lys251 (Arg) in Etr2p (Etr1p). Quantitative PCR analysis and reverse transcriptase-PCR indicated that both genes were expressed about equally in fermenting and ETR1 predominantly respiring yeast cells. Like the situation with ETR1, expression of ETR2 in respiration-deficient Saccharomyces cerevisiae mutant cells devoid of Ybr026p/Etr1p was able to restore growth on glycerol. Triclosan that is used as an antibacterial agent against fatty acid synthesis type II 2-enoyl thioester reductases inhibited growth of FabI overexpressing mutant yeast cells but was not able to inhibit respiratory growth of the ETR2- or ETR1-complemented mutant yeast cells. Resolving of crystal structures obtained via Etr2p and Etr1p co-crystallization indicated that all possible dimer variants occur in the same asymmetric unit, suggesting that similar dimer formation also takes place in vivo.

A novel design of whole-genome microarray probes for Saccharomyces cerevisiae which minimizes cross-hybridization

Background

Numerous DNA microarray hybridization experiments have been performed in yeast over the last years using either synthetic oligonucleotides or PCR-amplified coding sequences as probes. The design and quality of the microarray probes are of critical importance for hybridization experiments as well as subsequent analysis of the data.

Results

We present here a novel design of Saccharomyces cerevisiae microarrays based on a refined annotation of the genome and with the aim of reducing cross-hybridization between related sequences. An effort was made to design probes of similar lengths, preferably located in the 3'-end of reading frames. The sequence of each gene was compared against the entire yeast genome and optimal sub-segments giving no predicted cross-hybridization were selected. A total of 5660 novel probes (more than 97% of the yeast genes) were designed. For the remaining 143 genes, cross-hybridization was unavoidable. Using a set of 18 deletant strains, we have experimentally validated our cross-hybridization procedure. Sensitivity, reproducibility and dynamic range of these new microarrays have been measured. Based on this experience, we have written a novel program to design long oligonucleotides for microarray hybridizations of complete genome sequences.

Conclusions

A validated procedure to predict cross-hybridization in microarray probe design was defined in this work. Subsequently, a novel Saccharomyces cerevisiae microarray (which minimizes cross-hybridization) was designed and constructed. Arrays are available at Eurogentec S. A. Finally, we propose a novel design program, OliD, which allows automatic oligonucleotide design for microarrays. The OliD program is available from authors.

Perturbation of transcription elongation influences the fidelity of internal exon inclusion in Saccharomyces cerevisiae

Unknown mechanisms exist to ensure that exons are not skipped during biogenesis of mRNA. Studies have connected transcription elongation with regulated alternative exon inclusion. To determine whether the relative rates of transcription elongation and spliceosome assembly might play a general role in enforcing constitutive exon inclusion, we measured exon skipping for a natural two-intron gene in which the internal exon is constitutively included in the mRNA. Mutations in this gene that subtly reduce recognition of the intron 1 branchpoint cause exon skipping, indicating that rapid recognition of the first intron is important for enforcing exon inclusion. To test the role of transcription elongation, we treated cells to increase or decrease the rate of transcription elongation. Consistent with the "first come, first served" model, we found that exon skipping in vivo is inhibited when transcription is slowed by RNAP II mutants or when cells are treated with inhibitors of elongation. Expression of the elongation factor TFIIS stimulates exon skipping, and this effect is eliminated when lac repressor is targeted to DNA encoding the second intron. A mutation in U2 snRNA promotes exon skipping, presumably because a delay in recognition of the first intron allows elongating RNA polymerase to transcribe the downstream intron. This indicates that the relative rates of elongation and splicing are tuned so that the fidelity of exon inclusion is enhanced. These findings support a general role for kinetic coordination of transcription elongation and splicing during the transcription-dependent control of splicing.

Reinvestigation of the Saccharomyces cerevisiae genome annotation by comparison to the genome of a related fungus: Ashbya gossypii

BACKGROUND

The recently sequenced genome of the filamentous fungus Ashbya gossypii revealed remarkable similarities to that of the budding yeast Saccharomyces cerevisiae both at the level of homology and synteny (conservation of gene order). Thus, it became possible to reinvestigate the S. cerevisiae genome in the syntenic regions leading to an improved annotation.

RESULTS

We have identified 23 novel S. cerevisiae open reading frames (ORFs) as syntenic homologs of A. gossypii genes; for all but one, homologs are present in other eukaryotes including humans. Other comparisons identified 13 overlooked introns and suggested 69 potential sequence corrections resulting in ORF extensions or ORF fusions with improved homology to the syntenic A. gossypii homologs. Of the proposed corrections, 25 were tested and confirmed by resequencing. In addition, homologs of nearly 1,000 S. cerevisiae ORFs, presently annotated as hypothetical, were found in A. gossypii at syntenic positions and can therefore be considered as authentic genes. Finally, we suggest that over 400 S. cerevisiae ORFs that overlap other ORFs in S. cerevisiae and for which no homolog can be detected in A. gossypii should be regarded as spurious.

CONCLUSIONS

Although, the S. cerevisiae genome is rightly considered as one of the most accurately sequenced and annotated eukaryotic genomes, we have shown that it still benefits substantially from comparison to the completed sequence and syntenic gene map of A. gossypii, an evolutionarily related fungus. This type of approach will strongly support the annotation of more complex genomes such as the human and murine genomes.

Sequencing and comparison of yeast species to identify genes and regulatory elements

Identifying the functional elements encoded in a genome is one of the principal challenges in modern biology. Comparative genomics should offer a powerful, general approach. Here, we present a comparative analysis of the yeast Saccharomyces cerevisiae based on high-quality draft sequences of three related species (S. paradoxus, S. mikatae and S. bayanus). We first aligned the genomes and characterized their evolution, defining the regions and mechanisms of change. We then developed methods for direct identification of genes and regulatory motifs. The gene analysis yielded a major revision to the yeast gene catalogue, affecting approximately 15% of all genes and reducing the total count by about 500 genes. The motif analysis automatically identified 72 genome-wide elements, including most known regulatory motifs and numerous new motifs. We inferred a putative function for most of these motifs, and provided insights into their combinatorial interactions. The results have implications for genome analysis of diverse organisms, including the human.

C. elegans ORFeome version 1.1: experimental verification of the genome annotation and resource for proteome-scale protein expression

To verify the genome annotation and to create a resource to functionally characterize the proteome, we attempted to Gateway-clone all predicted protein-encoding open reading frames (ORFs), or the 'ORFeome,' of Caenorhabditis elegans. We successfully cloned approximately 12,000 ORFs (ORFeome 1.1), of which roughly 4,000 correspond to genes that are untouched by any cDNA or expressed-sequence tag (EST). More than 50% of predicted genes needed corrections in their intron-exon structures. Notably, approximately 11,000 C. elegans proteins can now be expressed under many conditions and characterized using various high-throughput strategies, including large-scale interactome mapping. We suggest that similar ORFeome projects will be valuable for other organisms, including humans.

Schizosaccharomyces pombe essential genes: a pilot study

After completion of the Schizosaccharomyces pombe genome sequence, we have carried out a pilot gene deletion project to assess the feasibility of a genome-wide deletion project and to estimate the percentage of essential genes. Using a PCR-based gene deletion procedure, we investigated 100 genes within a 253-kb region of chromosome II. Eight of nine genes located within a region of 18 kb could not be deleted, suggesting that systematic deletion of all fission yeast genes may be difficult to achieve using this PCR approach. The percentage of essential genes was found to be 17.5%. Further deletion of selected S. pombe genes revealed that whether a gene is essential or not is correlated with the timing of its appearance on the tree of life and its conservation within all branches of the tree. None of the investigated ancient genes in fission yeast that have been lost in the Saccharomyces cerevisiae lineage are essential. In agreement with S. cerevisiae and Caenorhabditis elegans genome analyses, our data suggest that natural selection has preferentially kept the genes required for vital functions. We propose that many of the essential eukaryotic genes appeared with the first eukaryotic cell and have remained conserved in all species.

Systematic discovery of new genes in the Saccharomyces cerevisiae genome

We used genome-wide comparative analysis of predicted protein sequences to identify many novel small genes, named smORFs for small open reading frames, within the budding yeast genome. Further analysis of 117 of these new genes showed that 84 are transcribed. We extended our analysis of one smORF conserved from yeast to human. This investigation provides an updated and comprehensive annotation of the yeast genome, validates additional concepts in the study of genomes in silico, and increases the expected numbers of coding sequences in a genome with the corresponding impact on future functional genomics and proteomics studies.

Phospho-regulation of kinetochore-microtubule attachments by the Aurora kinase Ipl1p

The Aurora kinase Ipl1p plays a crucial role in regulating kinetochore-microtubule attachments in budding yeast, but the underlying basis for this regulation is not known. To identify Ipl1p targets, we first purified 28 kinetochore proteins from yeast protein extracts. These studies identified five previously uncharacterized kinetochore proteins and defined two additional kinetochore subcomplexes. We then used mass spectrometry to identify 18 phosphorylation sites in 7 of these 28 proteins. Ten of these phosphorylation sites are targeted directly by Ipl1p, allowing us to identify a consensus phosphorylation site for an Aurora kinase. Our systematic mutational analysis of the Ipl1p phosphorylation sites demonstrated that the essential microtubule binding protein Dam1p is a key Ipl1p target for regulating kinetochore-microtubule attachments in vivo.

Functional analysis of yeast gene families involved in metabolism of vitamins B1and B6

In order to clarify their physiological functions, we have undertaken a characterization of the three-membered gene families SNZ1-3 and SNO1-3. In media lacking vitamin B(6), SNZ1 and SNO1 were both required for growth in certain conditions, but neither SNZ2, SNZ3, SNO2 nor SNO3 were required. Copies 2 and 3 of the gene products have, in spite of their extremely close sequence similarity, slightly different functions in the cell. We have also found that copies 2 and 3 are activated by the lack of thiamine and that the Snz proteins physically interact with the thiamine biosynthesis Thi5 protein family. Whereas copy 1 is required for conditions in which B(6) is essential for growth, copies 2 and 3 seem more related with B(1) biosynthesis during the exponential phase.

Genome-wide nuclear morphology screen identifies novel genes involved in nuclear architecture and gene-silencing in Saccharomyces cerevisiae

Organisation of the cell nucleus is crucial for the regulation of gene expression but little is known about how nuclei are structured. To address this issue, we designed a genomic screen to identify factors involved in nuclear architecture in Saccharomyces cerevisiae. This screen is based on microscopic monitoring of nuclear pore complexes and nucleolar proteins fused with the green fluorescent protein in a collection of approximately 400 individual deletion mutants. Among the 12 genes identified by this screen, most affect both the nuclear envelope and the nucleolar morphology. Corresponding gene products are localised preferentially to the nucleus or close to the nuclear periphery. Interestingly, these nuclear morphology alterations were associated with chromatin-silencing defects. These genes provide a molecular context to explore the functional link between nuclear architecture and gene silencing.

Parallel identification of new genes in Saccharomyces cerevisiae

Short open reading frames (ORFs) occur frequently in primary genome sequence. Distinguishing bona fide small genes from the tens of thousands of short ORFs is one of the most challenging aspects of genome annotation. Direct experimental evidence is often required. Here we use a combination of expression profiling and mass spectrometry to verify the independent transcription of 138 and the translation of 50 previously nonannotated genes in the Saccharomyces cerevisiae genome. Through combined evidence, we propose the addition of 62 new genes to the genome and provide experimental support for the inclusion of 10 previously identified genes.

Functional analysis of the Saccharomyces cerevisiae DUP240 multigene family reveals membrane-associated proteins that are not essential for cell viability

The DUP240 gene family of Saccharomyces cerevisiae is composed of 10 members. They encode proteins of about 240 amino acids which contain two predicted transmembrane domains. Database searches identified only one homologue in the closely related species Saccharomyces bayanus, indicating that the DUP240 genes encode proteins specific to Saccharomyces sensu stricto. The short-flanking homology PCR gene-replacement strategy with a variety of selective markers for replacements, and classical genetic methods, were used to generate strains deleted for all 10 DUP240 genes. All of the knock-out strains were viable and had similar growth kinetics to the wild-type. Two-hybrid screens, hSos1p fusions and GFP fusions were carried out; the results indicated that the Dup240 proteins are membrane associated, and that some of them are concentrated around the plasma membrane.

Mnd1p: an evolutionarily conserved protein required for meiotic recombination

We used a functional genomics approach to identify a gene required for meiotic recombination, YGL183c or MND1. MND1 was spliced in meiotic cells, extending the annotated YGL183c ORF N terminus by 45 aa. Saccharomyces cerevisiae mnd1-1 mutants, in which the majority of the MND1 coding sequence was removed, arrested before the first meiotic division with a phenotype reminiscent of dmc1 mutants. Physical and genetic analysis showed that these cells initiated recombination, but did not form heteroduplex DNA or double Holliday junctions, suggesting that Mnd1p is involved in strand invasion. Orthologs of MND1 were identified in protists, several yeasts, plants, and mammals, suggesting that its function has been conserved throughout evolution.

How many protein-coding genes are there in the Saccharomyces cerevisiae genome?

We have compared the results of estimations of the total number of protein-coding genes in the Saccharomyces cerevisiae genome, which have been obtained by many laboratories since the yeast genome sequence was published in 1996. We propose that there are 5300-5400 genes in the genome. This makes the first estimation of the number of intronless ORFs longer than 100 codons, based on the features of the set of genes with phenotypes known in 1997 to be correct. This estimation assumed that the set of the first 2300 genes with known phenotypes was representative for the whole set of protein-coding genes in the genome. The same method used in this paper for the approximation of the total number of protein-coding sequences among more than 40 000 ORFs longer than 20 codons gives a result that is only slightly higher. This suggests that there are still some non-coding ORFs in the databases and a few dozen small ORFs, not yet annotated, which probably code for proteins.

A question of size: the eukaryotic proteome and the problems in defining it

We discuss the problems in defining the extent of the proteomes for completely sequenced eukaryotic organisms (i.e. the total number of protein-coding sequences), focusing on yeast, worm, fly and human. (i) Six years after completion of its genome sequence, the true size of the yeast proteome is still not defined. New small genes are still being discovered, and a large number of existing annotations are being called into question, with these questionable ORFs (qORFs) comprising up to one-fifth of the 'current' proteome. We discuss these in the context of an ideal genome-annotation strategy that considers the proteome as a rigorously defined subset of all possible coding sequences ('the orfome'). (ii) Despite the greater apparent complexity of the fly (more cells, more complex physiology, longer lifespan), the nematode worm appears to have more genes. To explain this, we compare the annotated proteomes of worm and fly, relating to both genome-annotation and genome evolution issues. (iii) The unexpectedly small size of the gene complement estimated for the complete human genome provoked much public debate about the nature of biological complexity. However, in the first instance, for the human genome, the relationship between gene number and proteome size is far from simple. We survey the current estimates for the numbers of human genes and, from this, we estimate a range for the size of the human proteome. The determination of this is substantially hampered by the unknown extent of the cohort of pseudogenes ('dead' genes), in combination with the prevalence of alternative splicing. (Further information relating to yeast is available at http://genecensus.org/yeast/orfome)

A small reservoir of disabled ORFs in the yeast genome and its implications for the dynamics of proteome evolution

We surveyed the sequenced Saccharomyces cerevisiae genome (strain S288C) comprehensively for open reading frames (ORFs) that could encode full-length proteins but contain obvious mid-sequence disablements (frameshifts or premature stop codons). These pseudogenic features are termed disabled ORFs (dORFs). Using homology to annotated yeast ORFs and non-yeast proteins plus a simple region extension procedure, we have found 183 dORFs. Combined with the 38 existing annotations for potential dORFs, we have a total pool of up to 221 dORFs, corresponding to less than approximately 3% of the proteome. Additionally, we found 20 pairs of annotated ORFs for yeast that could be merged into a single ORF (termed a mORF) by read-through of the intervening stop codon, and may comprise a complete ORF in other yeast strains. Focussing on a core pool of 98 dORFs with a verifying protein homology, we find that most dORFs are substantially decayed, with approximately 90% having two or more disablements, and approximately 60% having four or more. dORFs are much more yeast-proteome specific than live yeast genes (having about half the chance that they are related to a non-yeast protein). They show a dramatically increased density at the telomeres of chromosomes, relative to genes. A microarray study shows that some dORFs are expressed even though they carry multiple disablements, and thus may be more resistant to nonsense-mediated decay. Many of the dORFs may be involved in responding to environmental stresses, as the largest functional groups include growth inhibition, flocculation, and the SRP/TIP1 family. Our results have important implications for proteome evolution. The characteristics of the dORF population suggest the sorts of genes that are likely to fall in and out of usage (and vary in copy number) in a strain-specific way and highlight the role of subtelomeric regions in engendering this diversity. Our results also have important implications for the effects of the [PSI+] prion. The dORFs disabled by only a single stop and the mORFs (together totalling 35) provide an estimate for the extent of the sequence population that can be resurrected readily through the demonstrated ability of the [PSI+] prion to cause nonsense-codon read-through. Also, the dORFs and mORFs that we find have properties (e.g. growth inhibition, flocculation, vanadate resistance, stress response) that are potentially related to the ability of [PSI+] to engender substantial phenotypic variation in yeast strains under different environmental conditions. (See genecensus.org/pseudogene for further information.)

The genome sequence of Schizosaccharomyces pombe

We have sequenced and annotated the genome of fission yeast (Schizosaccharomyces pombe), which contains the smallest number of protein-coding genes yet recorded for a eukaryote: 4,824. The centromeres are between 35 and 110 kilobases (kb) and contain related repeats including a highly conserved 1.8-kb element. Regions upstream of genes are longer than in budding yeast (Saccharomyces cerevisiae), possibly reflecting more-extended control regions. Some 43% of the genes contain introns, of which there are 4,730. Fifty genes have significant similarity with human disease genes; half of these are cancer related. We identify highly conserved genes important for eukaryotic cell organization including those required for the cytoskeleton, compartmentation, cell-cycle control, proteolysis, protein phosphorylation and RNA splicing. These genes may have originated with the appearance of eukaryotic life. Few similarly conserved genes that are important for multicellular organization were identified, suggesting that the transition from prokaryotes to eukaryotes required more new genes than did the transition from unicellular to multicellular organization.