Detection and Characterization of Fungal-Specific Proteins in Saccharomyces cerevisiae

In this study, I searched for fungal-specific proteins in the genome of the budding yeast Saccharomyces cerevisiae, inferred from a comparison of amino acid sequences. I used the GTOP (Genomes to Protein structures and functions) database of the DDBJ (DNA Data Bank of Japan), which consists of 21 genomes from Archaea, 203 genomes from Bacteria, and 50 genomes from Eucarya (including 18 fungal genomes). Among 5,874 proteins of S. cerevisiae, 1,551 have homologs only in Eucarya, and 504 of the 1,551 have homologs only in fungi. To find fungal-specific proteins, homologs of the homologs have been searched repeatedly. As a result, 132 of the 504 are characterized as fungal-specific proteins. The genes encoding the 132 fungal-specific proteins are not included in the list of essential genes for viability in the S. cerevisiae genome deletion project. Among the 132 proteins, 99 are S. cerevisiae-specific, and no protein that is distributed among 10 or more of the 18 fungal species exists. In addition, most of the fungal-specific proteins are very small and functionally unknown. My results show that the fungal-specific proteins have short evolutionary histories, suggesting that S. cerevisiae produces novel proteins and that ancestral fungi also produced small proteins most of which have disappeared or have been combined with other proteins during fungal evolution.     

Comparison of the Yeast Proteome to Other Fungal Genomes to Find Core Fungal Genes

The purpose of this research was to search for evolutionarily conserved fungal sequences to test the hypothesis that fungi have a set of core genes that are not found in other organisms, as these genes may indicate what makes fungi different from other organisms. By comparing 6355 predicted or known yeast (Saccharomyces cerevisiae) genes to the genomes of 13 other fungi using Standalone TBLASTN at an e-value <1E-5, a list of 3340 yeast genes was obtained with homologs present in at least 12 of 14 fungal genomes. By comparing these common fungal genes to complete genomes of animals (Fugu rubripes, Caenorhabditis elegans), plants (Arabidopsis thaliana, Oryza sativa), and bacteria (Agrobacterium tumefaciens, Xylella fastidiosa), a list of common fungal genes with homologs in these plants, animals, and bacteria was produced (938 genes), as well as a list of exclusively fungal genes without homologs in these other genomes (60 genes). To ensure that the 60 genes were exclusively fungal, these were compared using TBLASTN to the major sequence databases at GenBank: NR (nonredundant), EST (expressed sequence tags), GSS (genome survey sequences), and HTGS (unfinished high-throughput genome sequences). This resulted in 17 yeast genes with homologs in other fungal genomes, but without known homologs in other organisms. These 17 core, fungal genes were not found to differ from other yeast genes in GC content or codon usage patterns. More intensive study is required of these 17 genes and other common fungal genes to discover unique features of fungi compared to other organisms.Reviewing Editor: Prof. David Gottman     

Identification and analysis of cation channel homologues in human pathogenic fungi

Fungi are major causes of human, animal and plant disease. Human fungal infections can be fatal, but there are limited options for therapy, and resistance to commonly used anti-fungal drugs is widespread. The genomes of many fungi have recently been sequenced, allowing identification of proteins that may become targets for novel therapies. We examined the genomes of human fungal pathogens for genes encoding homologues of cation channels, which are prominent drug targets. Many of the fungal genomes examined contain genes encoding homologues of potassium (K(+)), calcium (Ca(2+)) and transient receptor potential (Trp) channels, but not sodium (Na(+)) channels or ligand-gated channels. Some fungal genomes contain multiple genes encoding homologues of K(+) and Trp channel subunits, and genes encoding novel homologues of voltage-gated K(v) channel subunits are found in Cryptococcus spp. Only a single gene encoding a homologue of a plasma membrane Ca(2+) channel was identified in the genome of each pathogenic fungus examined. These homologues are similar to the Cch1 Ca(2+) channel of Saccharomyces cerevisiae. The genomes of Aspergillus spp. and Cryptococcus spp., but not those of S. cerevisiae or the other pathogenic fungi examined, also encode homologues of the mitochondrial Ca(2+) uniporter (MCU). In contrast to humans, which express many K(+), Ca(2+) and Trp channels, the genomes of pathogenic fungi encode only very small numbers of K(+), Ca(2+) and Trp channel homologues. Furthermore, the sequences of fungal K(+), Ca(2+), Trp and MCU channels differ from those of human channels in regions that suggest differences in regulation and susceptibility to drugs.     

The machinery for cell polarity, cell morphogenesis, and the cytoskeleton in the Basidiomycete fungus Ustilago maydis-a survey of the genome sequence

Ustilago maydis, a Basidiomycete fungus that infects maize, exhibits two basic morphologies, a yeast-like and a filamentous form. The yeast-like cell is elongated, divides by budding, and the bud grows by tip extension. The filamentous form divides at the apical cell and grows by tip extension. The repertoire of morphologies is increased during interaction with its host, suggesting that plant signals play an important role in generation of additional morphologies. We have used Saccharomyces cerevisiae and Schizosaccharomyces pombe genes known to play a role in cell polarity and morphogenesis, and in the cytoskeleton as probes to survey the U. maydis genome. We have found that most of the yeast machinery is conserved in U. maydis, albeit the degree of similarity varies from strong to weak. The U. maydis genome contains the machinery for recognition and interpretation of the budding yeast axial and bipolar landmarks; however, genes coding for some of the landmark proteins are absent. Genes coding for cell polarity establishment, exocytosis, actin and microtubule organization, microtubule plus-end associated proteins, kinesins, and myosins are also present. Genes not present in S. cerevisiae and S. pombe include a homolog of mammalian Rac, a hybrid myosin-chitin synthase, and several kinesins that exhibit more similarity to their mammalian counterparts. We also used the U. maydis genes identified in this analysis to search other fungal and other eukaryotic genomes to identify the closest homologs. In most cases, not surprisingly, the closest homolog is among filamentous fungi, not the yeasts, and in some cases it is among mammals.     

ATG genes involved in non-selective autophagy are conserved from yeast to man, but the selective Cvt and pexophagy pathways also require organism-specific genes

ATG genes encode proteins that are required for macroautophagy, the Cvt pathway and/or pexophagy. Using the published Atg protein sequences, we have screened protein and DNA databases to identify putative functional homologs (orthologs) in 21 fungal species (yeast and filamentous fungi) of which the genome sequences were available. For comparison with Atg proteins in higher eukaryotes, also an analysis of Arabidopsis thaliana and Homo sapiens databases was included. This analysis demonstrated that Atg proteins required for non-selective macroautophagy are conserved from yeast to man, stressing the importance of this process in cell survival and viability. The A. thaliana and human genomes encode multiple proteins highly similar to specific fungal Atg proteins (paralogs), possibly representing cell type-specific isoforms. The Atg proteins specifically involved in the Cvt pathway and/or pexophagy showed poor conservation, and were generally not present in A. thaliana and man. Furthermore, Atg19, the receptor of Cvt cargo, was only detected in Saccharomyces cerevisiae. Nevertheless, Atg11, a protein that links receptor-bound cargo (peroxisomes, the Cvt complex) to the autophagic machinery was identified in all yeast species and filamentous fungi under study. This suggests that in fungi an organism-specific form of selective autophagy may occur, for which specialized Atg proteins have evolved.     

Trypsin-like proteins of the fungi as possible markers of pathogenicity

Sequences of peptidases with conserved motifs around the active site residues that are characteristic of trypsins (similar to trypsin peptidases, STP) were obtained from publicly-available fungal genomes and related databases. Among the 75 fungal genomes, 29 species of parasitic Ascomycota contained genes encoding STP and their homologs. Searches of non-redundant protein sequences, patented protein sequences, and expressed sequence tags resulted in another 18 STP sequences in 10 fungal species from Ascomycota, Basidiomycota, and Zygomycota. A comparison of fungi species containing STP sequences revealed that almost all are pathogens of plants, animals or fungi. A comparison of the primary structure of homologous proteins, including the residues responsible for substrate binding and specificity of the enzyme, revealed three groups of homologous sequences, all presumably from S1 family: trypsin-like peptidases, chymotrypsin-like peptidases and serine peptidases with unknown substrate specificity. Homologs that are presumably functionally inactive were predicted in all groups. The results in general support the hypothesis that the expression of trypsin-like peptidases in fungi represents a marker of fungal phytopathogenicity. A phylogenetic tree was constructed using peptidase and homolog amino acid sequences, demonstrating that all have noticeable differences and almost immediately deviate from the common root. Therefore, we conclude that the changes that occurred in STP of pathogenic fungi in the course of evolution represent specific adaptations to proteins of their respective hosts, and mutations in peptidase genes are important components of life-style changes and taxonomic divergence.     

Highly diverged homologs of Saccharomyces cerevisiae mitochondrial mRNA-specific translational activators have orthologous functions in other budding yeasts

Translation of mitochondrially coded mRNAs in Saccharomyces cerevisiae depends on membrane-bound mRNA-specific activator proteins, whose targets lie in the mRNA 5'-untranslated leaders (5'-UTLs). In at least some cases, the activators function to localize translation of hydrophobic proteins on the inner membrane and are rate limiting for gene expression. We searched unsuccessfully in divergent budding yeasts for orthologs of the COX2- and COX3-specific translational activator genes, PET111, PET54, PET122, and PET494, by direct complementation. However, by screening for complementation of mutations in genes adjacent to the PET genes in S. cerevisiae, we obtained chromosomal segments containing highly diverged homologs of PET111 and PET122 from Saccharomyces kluyveri and of PET111 from Kluyveromyces lactis. All three of these genes failed to function in S. cerevisiae. We also found that the 5'-UTLs of the COX2 and COX3 mRNAs of S. kluyveri and K. lactis have little similarity to each other or to those of S. cerevisiae. To determine whether the PET111 and PET122 homologs carry out orthologous functions, we deleted them from the S. kluyveri genome and deleted PET111 from the K. lactis genome. The pet111 mutations in both species prevented COX2 translation, and the S. kluyveri pet122 mutation prevented COX3 translation. Thus, while the sequences of these translational activator proteins and their 5'-UTL targets are highly diverged, their mRNA-specific functions are orthologous.     

Fungal apoptosis: function, genes and gene function

Cells of all living organisms are programmed to self-destruct under certain conditions. The most well known form of programmed cell death is apoptosis, which is essential for proper development in higher eukaryotes. In fungi, apoptotic-like cell death occurs naturally during aging and reproduction, and can be induced by environmental stresses and exposure to toxic metabolites. The core apoptotic machinery in fungi is similar to that in mammals, but the apoptotic network is less complex and of more ancient origin. Only some of the mammalian apoptosis-regulating proteins have fungal homologs, and the number of protein families is drastically reduced. Expression in fungi of animal proteins that do not have fungal homologs often affects apoptosis, suggesting functional conservation of these components despite the absence of protein-sequence similarity. Functional analysis of Saccharomyces cerevisiae apoptotic genes, and more recently of those in some filamentous species, has revealed partial conservation, along with substantial differences in function and mode of action between fungal and human proteins. It has been suggested that apoptotic proteins might be suitable targets for novel antifungal treatments. However, implementation of this approach requires a better understanding of fungal apoptotic networks and identification of the key proteins regulating apoptotic-like cell death in fungi.     

Conserved processes and lineage-specific proteins in fungal cell wall evolution

The cell wall is a defining organelle that differentiates fungi from its sister clades in the opisthokont superkingdom. With a sensitive technique to align low-complexity protein sequences, we have identified 187 cell wall-related proteins in Saccharomyces cerevisiae and determined the presence or absence of homologs in 17 other fungal genomes. There were both conserved and lineage-specific cell wall proteins, and the degree of conservation was strongly correlated with protein function. Some functional classes were poorly conserved and lineage specific: adhesins, structural wall glycoprotein components, and unannotated open reading frames. These proteins are primarily those that are constituents of the walls themselves. On the other hand, glycosyl hydrolases and transferases, proteases, lipases, proteins in the glycosyl phosphatidyl-inositol-protein synthesis pathway, and chaperones were strongly conserved. Many of these proteins are also conserved in other eukaryotes and are associated with wall synthesis in plants. This gene conservation, along with known similarities in wall architecture, implies that the basic architecture of fungal walls is ancestral to the divergence of the ascomycetes and basidiomycetes. The contrasting lineage specificity of wall resident proteins implies diversification. Therefore, fungal cell walls consist of rapidly diversifying proteins that are assembled by the products of an ancestral and conserved set of genes.     

Predicting the distribution,conservation, and functions of SNAREs and related proteins in fungi

Hyphal tip growth, the hallmark of the fungi, requires highly polarized and localized exocytosis, but how this requirement is met is unknown. Members of conserved protein families called SNAREs and Rabs mediate vesicle trafficking and fusion at virtually every step of the intracellular pathway in all examined eukaryotes. We have searched the available nearly complete fungal genomes, established the presence or absence of members of the SNARE and Rab families in these genomes, and predicted their evolutionary relationships to one another. Comparisons with the extensively studied Saccharomyces cerevisiae indicate that, in general, most of the members of these families (including those involved in mediating exocytosis) are conserved. The presence of exceptional SNAREs and Rabs in some fungi that are not conserved in S. cerevisiae may be indicative of specialized steps that occur in these fungi. The implications of these findings for current tip growth models are discussed.     

ATG Genes Involved in Non-Selective Autophagy are Conserved from Yeast to Man,but the Selective Cvt and Pexophagy Pathways also Require Organism-Specific Genes

ATG genes encode proteins that are required for macroautophagy, the Cvt pathway and/or pexophagy. Using the published Atg protein sequences, we have screened protein and DNA databases to identify putative functional homologs (orthologs) in 21 fungal species (yeast and filamentous fungi) of which the genome sequences were available. For comparison with Atg proteins in higher eukaryotes, also the genomes of Arabidopsis thaliana and Homo sapiens were included. This analysis demonstrated that Atg proteins required for non-selective macroautophagy are conserved from yeast to man, stressing the importance of this process in cell survival and viability. Remarkably, the A. thaliana and human genomes encode multiple proteins highly similar to specific Atg proteins (paralogs), the function of which is unknown. The Atg proteins specifically involved in the Cvt pathway and/or pexophagy showed poor conservation, and were generally not present in A. thaliana and man. Furthermore, the receptor of Cvt cargo, Atg19, was only detected in S. cerevisiae. Nevertheless, Atg11, a protein that links receptor-bound cargo (peroxisomes, Cvt bodies) to the autophagic machinery was identified in all yeast species and filamentous fungi under study. This suggests that in fungi an organism-specific form of selective autophagy may occur, for which specialized Atg proteins have evolved.     

Analysis of the pdx-1 (snz-1/sno-1) region of the Neurospora crassa genome: correlation of pyridoxine-requiring phenotypes with mutations in two structural genes

We report the analysis of a 36-kbp region of the Neurospora crassa genome, which contains homologs of two closely linked stationary phase genes, SNZ1 and SNO1, from Saccharomyces cerevisiae. Homologs of SNZ1 encode extremely highly conserved proteins that have been implicated in pyridoxine (vitamin B6) metabolism in the filamentous fungi Cercospora nicotianae and in Aspergillus nidulans. In N. crassa, SNZ and SNO homologs map to the region occupied by pdx-1 (pyridoxine requiring), a gene that has been known for several decades, but which was not sequenced previously. In this study, pyridoxine-requiring mutants of N. crassa were found to possess mutations that disrupt conserved regions in either the SNZ or SNO homolog. Previously, nearly all of these mutants were classified as pdx-1. However, one mutant with a disrupted SNO homolog was at one time designated pdx-2. It now appears appropriate to reserve the pdx-1 designation for the N. crassa SNZ homolog and pdx-2 for the SNO homolog. We further report annotation of the entire 36,030-bp region, which contains at least 12 protein coding genes, supporting a previous conclusion of high gene densities (12,000-13,000 total genes) for N. crassa. Among genes in this region other than SNZ and SNO homologs, there was no evidence of shared function. Four of the genes in this region appear to have been lost from the S. cerevisiae lineage.     

De novo origination of a new protein-coding gene in Saccharomyces cerevisiae

Origination of new genes is an important mechanism generating genetic novelties during the evolution of an organism. Processes of creating new genes using preexisting genes as the raw materials are well characterized, such as exon shuffling, gene duplication, retroposition, gene fusion, and fission. However, the process of how a new gene is de novo created from noncoding sequence is largely unknown. On the basis of genome comparison among yeast species, we have identified a new de novo protein-coding gene, BSC4 in Saccharomyces cerevisiae. The BSC4 gene has an open reading frame (ORF) encoding a 132-amino-acid-long peptide, while there is no homologous ORF in all the sequenced genomes of other fungal species, including its closely related species such as S. paradoxus and S. mikatae. The functional protein-coding feature of the BSC4 gene in S. cerevisiae is supported by population genetics, expression, proteomics, and synthetic lethal data. The evidence suggests that BSC4 may be involved in the DNA repair pathway during the stationary phase of S. cerevisiae and contribute to the robustness of S. cerevisiae, when shifted to a nutrient-poor environment. Because the corresponding noncoding sequences in S. paradoxus, S. mikatae, and S. bayanus also transcribe, we propose that a new de novo protein-coding gene may have evolved from a previously expressed noncoding sequence.     

Fungal Smn and Spf30 homologues are mainly present in filamentous fungi and genomes with many introns: implications for spinal muscular atrophy

Spinal muscular atrophy is an important rare genetic disease characterized by the loss of motor neurons, where the main gene responsible is smn1. Orthologous genes have only been characterized in a single fungal genome: Schizosaccharomyces pombe. We have searched for putative SMN orthologues in publically available fungal genomes, finding that they are predominately present in filamentous fungi. SMN binding partners and the SPF30 SMN paralogue, which are all involved in mRNA splicing, were found to be present in a similar but non-identical subset of fungal genomes. The Saccharomycces cerevisiae yeast genome contains neither smn1 orthologues nor paralogues and it has been suggested that this might be related to the low number of introns in this yeast. Here we have tested this hypothesis by looking at other fungal genomes. Significantly, we find that fungal genomes with high numbers of introns also possess an SMN orthologue or at least its paralogue, SPF30.     

Phylogenetic, genomic organization and expression analysis of hydrophobin genes in the ectomycorrhizal basidiomycete Laccaria bicolor

Hydrophobins are morphogenetic, small secreted hydrophobic fungal proteins produced in response to changing development and environmental conditions. These proteins are important in the interaction between certain fungi and their hosts. In mutualistic ectomycorrhizal fungi several hydrophobins form a subclass of mycorrhizal-induced small secreted proteins that are likely to be critical in the formation of the symbiotic interface with host root cells. In this study, two genomes of the ectomycorrhizal basidiomycete Laccaria bicolor strains S238N-H82 (from North America) and 81306 (from Europe) were surveyed to construct a comprehensive genome-wide inventory of hydrophobins and to explore their characteristics and roles during host colonization. The S238N-H82 L. bicolor hydrophobin gene family is composed of 12 genes while the 81306 strain encodes nine hydrophobins, all corresponding to class I hydrophobins. The three extra hydrophobin genes encoded by the S238N-H82 genome likely arose via gene duplication and are bordered by transposon rich regions. Expression profiles of the hydrophobin genes of L. bicolor varied greatly depending on life stage (e.g. free living mycelium vs. root colonization) and on the host root environment. We conclude from this study that the complex diversity and range of expression profiles of the Laccaria hydrophobin multi-gene family have likely been a selective advantage for this mutualist in colonizing a wide range of host plants.     

Fungal genome sequencing: basic biology to biotechnology

The genome sequences provide a first glimpse into the genomic basis of the biological diversity of filamentous fungi and yeast. The genome sequence of the budding yeast, Saccharomyces cerevisiae, with a small genome size, unicellular growth, and rich history of genetic and molecular analyses was a milestone of early genomics in the 1990s. The subsequent completion of fission yeast, Schizosaccharomyces pombe and genetic model, Neurospora crassa initiated a revolution in the genomics of the fungal kingdom. In due course of time, a substantial number of fungal genomes have been sequenced and publicly released, representing the widest sampling of genomes from any eukaryotic kingdom. An ambitious genome-sequencing program provides a wealth of data on metabolic diversity within the fungal kingdom, thereby enhancing research into medical science, agriculture science, ecology, bioremediation, bioenergy, and the biotechnology industry. Fungal genomics have higher potential to positively affect human health, environmental health, and the planet’s stored energy. With a significant increase in sequenced fungal genomes, the known diversity of genes encoding organic acids, antibiotics, enzymes, and their pathways has increased exponentially. Currently, over a hundred fungal genome sequences are publicly available; however, no inclusive review has been published. This review is an initiative to address the significance of the fungal genome-sequencing program and provides the road map for basic and applied research.     

Incipient mitochondrial evolution in yeasts. I. The physical map and gene order of Saccharomyces douglasii mitochondrial DNA discloses a translocation of a segment of 15,000 base-pairs and the presence of new introns in comparison with Saccharomyces cerevisiae

We have determined the physical and genetic map of the 73,000 base-pair mitochondrial genome of a novel yeast species Saccharomyces douglasii. Most of the protein and RNA-coding genes known to be present in the mitochondrial DNA of Saccharomyces cerevisiae have been identified and located on the S. douglasii mitochondrial genome. The nuclear genomes of the two species are thought to have diverged some 50 to 80 million years ago and their nucleo-mitochondrial hybrids are viable but respiratorily deficient. The mitochondrial genome of S. douglasii displays many interesting features in comparison with that of S. cerevisiae. The three mosaic genes present in both genomes are quite different with regard to their structure. The S. douglasii COXI gene has two new introns and is missing the five introns of the S. cerevisiae gene. The S. douglasii cytochrome b gene has one new intron and lacks two introns of the S. cerevisiae gene. Finally, the L-rRNA gene of S. douglasii, like that of S. cerevisiae, has one intron of which the structure is different. Another salient feature of the S. douglasii mitochondrial genome reported here is that the gene order is different in comparison with S. cerevisiae mitochondrial DNA. In particular, a segment of approximately 15,000 base-pairs including the genes coding for COXIII and S-rRNA has been translocated to a position between the genes coding for varl and L-rRNA.     

The mitochondrial genome of Arabidopsis is composed of both native and immigrant information

Plants contain large mitochondrial genomes, which are several times as complex as those in animals, fungi or algae. However, genome size is not correlated with information content. The mitochondrial genome (mtDNA) of Arabidopsis specifies only 58 genes in 367 kb, whereas the 184 kb mtDNA in the liverwort Marchantia polymorpha codes for 66 genes, and the 58 kb genome in the green alga Prototheca wickerhamii encodes 63 genes. In Arabidopsis’ mtDNA, genes for subunits of complex II, for several ribosomal proteins and for 16 tRNAs are missing, some of which have been transferred recently to the nuclear genome. Numerous integrated fragments originate from alien genomes, including 16 sequence stretches of plastid origin, 41 fragments of nuclear (retro)transposons and two fragments of fungal viruses. These immigrant sequences suggest that the large size of plant mitochondrial genomes is caused by secondary expansion as a result of integration and propagation, and is thus a derived trait established during the evolution of land plants.     

The mitochondrial genome of the arbuscular mycorrhizal fungus Gigaspora margarita reveals two unsuspected trans-splicing events of group I introns

? Arbuscular mycorrhizal fungi (AMF) are ubiquitous organisms that benefit ecosystems through the establishment of an association with the roots of most plants: the mycorrhizal symbiosis. Despite their ecological importance, however, these fungi have been poorly studied at the genome level. ? In this study, total DNA from the AMF Gigaspora margarita was subjected to a combination of 454 and Illumina sequencing, and the resulting reads were used to assemble its mitochondrial genome de novo. This genome was annotated and compared with those of other relatives to better comprehend the evolution of the AMF lineage. ? The mitochondrial genome of G. margarita is unique in many ways, exhibiting a large size (97 kbp) and elevated GC content (45%). This genome also harbors molecular events that were previously unknown to occur in fungal mitochondrial genomes, including trans-splicing of group I introns from two different genes coding for the first subunit of the cytochrome oxidase and for the small subunit of the rRNA. ? This study reports the second published genome from an AMF organelle, resulting in relevant DNA sequence information from this poorly studied fungal group, and providing new insights into the frequency, origin and evolution of trans-spliced group I introns found across the mitochondrial genomes of distantly related organisms.