Random exploration of the Kluyveromyces lactis genome and comparison with that of Saccharomyces cerevisiae

Deep functional analysis of synII, a 770-kilobase synthetic yeast chromosome

Here, we report the successful design, construction, and characterization of a 770-kilobase synthetic yeast chromosome II (synII). Our study incorporates characterization at multiple levels-including phenomics, transcriptomics, proteomics, chromosome segregation, and replication analysis-to provide a thorough and comprehensive analysis of a synthetic chromosome. Our Trans-Omics analyses reveal a modest but potentially relevant pervasive up-regulation of translational machinery observed in synII, mainly caused by the deletion of 13 transfer RNAs. By both complementation assays and SCRaMbLE (synthetic chromosome rearrangement and modification by loxP-mediated evolution), we targeted and debugged the origin of a growth defect at 37°C in glycerol medium, which is related to misregulation of the high-osmolarity glycerol response. Despite the subtle differences, the synII strain shows highly consistent biological processes comparable to the native strain.

Genome-wide metabolic (re-) annotation of Kluyveromyces lactis

BACKGROUND

Even before having its genome sequence published in 2004, Kluyveromyces lactis had long been considered a model organism for studies in genetics and physiology. Research on Kluyveromyces lactis is quite advanced and this yeast species is one of the few with which it is possible to perform formal genetic analysis. Nevertheless, until now, no complete metabolic functional annotation has been performed to the proteins encoded in the Kluyveromyces lactis genome.

RESULTS

In this work, a new metabolic genome-wide functional re-annotation of the proteins encoded in the Kluyveromyces lactis genome was performed, resulting in the annotation of 1759 genes with metabolic functions, and the development of a methodology supported by merlin (software developed in-house). The new annotation includes novelties, such as the assignment of transporter superfamily numbers to genes identified as transporter proteins. Thus, the genes annotated with metabolic functions could be exclusively enzymatic (1410 genes), transporter proteins encoding genes (301 genes) or have both metabolic activities (48 genes). The new annotation produced by this work largely surpassed the Kluyveromyces lactis currently available annotations. A comparison with KEGG's annotation revealed a match with 844 (~90%) of the genes annotated by KEGG, while adding 850 new gene annotations. Moreover, there are 32 genes with annotations different from KEGG.

CONCLUSIONS

The methodology developed throughout this work can be used to re-annotate any yeast or, with a little tweak of the reference organism, the proteins encoded in any sequenced genome. The new annotation provided by this study offers basic knowledge which might be useful for the scientific community working on this model yeast, because new functions have been identified for the so-called metabolic genes. Furthermore, it served as the basis for the reconstruction of a compartmentalized, genome-scale metabolic model of Kluyveromyces lactis, which is currently being finished.

Comparative study on synteny between yeasts and vertebrates

We studied synteny conservation between 18 yeast species and 13 vertebrate species in order to provide a comparative analysis of the chromosomal plasticity in these 2 phyla. By computing the regions of conserved synteny between all pairwise combinations of species within each group, we show that in vertebrates, the number of conserved synteny blocks exponentially increases along with the divergence between orthologous protein and that concomitantly; the number of genes per block exponentially decreases. The same trends are found in yeasts but only when the mean protein divergence between orthologs remains below 36%. When the average protein divergence exceeds this threshold, the total number of recognizable synteny blocks gradually decreases due to the repeated accumulation of rearrangements. We also show that rearrangement rates are on average 3-fold higher in vertebrates than in yeasts, and are estimated to be of 2 rearrangements/Myr. However, the genome sizes being on average 200 times larger in vertebrates than in yeasts, the normalized rates of chromosome rearrangements (per Mb) are about 50-fold higher in yeast than in vertebrate genomes.

Genome-wide analysis of Kluyveromyces lactis in wild-type and rag2 mutant strains

The use of heterologous DNA arrays from Saccharomyces cerevisiae has been tested and revealed as a suitable tool to compare the transcriptomes of S. cerevisiae and Kluyveromyces lactis, two yeasts with notable differences in their respirofermentative metabolism. The arrays have also been applied to study the changes in the K. lactis transcriptome owing to mutation in the RAG2 gene coding for the glycolytic enzyme phosphoglucose isomerase. Comparison of the rag2 mutant growing in 2% glucose versus 2% fructose has been used as a model to elucidate the importance of transcriptional regulation of metabolic routes, which may be used to reoxidize the NADPH produced in the pentose phosphate pathway. At this transcriptional level, routes related to the oxidative stress response become an interesting alternative for NADPH use.

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.

The role of adaptor protein Ste50-dependent regulation of the MAPKKK Ste11 in multiple signalling pathways of yeast

In Saccharomyces cerevisiae, Ste50 functions in cell signalling between the activated G protein and the mitogen-activated protein kinase (MAPK) kinase kinase (MAPKKK) Ste11. ScSte50 is an essential component of three MAPK-mediated signalling pathways, which control the mating response, invasive/filamentous growth and osmotolerance (HOG pathway), respectively. ScSte50 signalling may also contribute to cell wall integrity in vegetative cells. The protein contains a sterile alpha motif (SAM) and a putative Ras-associated domain (RAD), which are essential for signal transduction. Ste50 and Ste11 interact constitutively via their SAM regions. Ste50 interacts weakly and probably transiently with the pheromone receptor-bound heterotrimeric G protein G(alpha beta gamma), and with the small G proteins Cdc42, Ras1 and Ras2. It is specifically the RAD region of Ste50 that mediates the interactions with Cdc42 and Ras. Homologues of ScSTE50 are also found in other fungi, like S. kluyveri, Hansenula polymorpha, Candida albicans and Neurospora crassa. In this review, the role of Ste50 as an adaptor that links the G protein-associated Cdc42-Ste20 kinase complex to the effector kinase Ste11 and thus modulates signal transduction, especially in the pheromone-response pathway of S. cerevisiae, is discussed.

Krh1p and Krh2p act downstream of the Gpa2p G(alpha) subunit to negatively regulate haploid invasive growth

The yeast G(alpha) subunit Gpa2p and its coupled receptor Gpr1p function in a signaling pathway that is required for the transition to pseudohyphal and invasive growth. A two-hybrid screen using a constitutively active allele of GPA2 identified the KRH1 gene as encoding a potential binding partner of Gpa2p. Strains containing deletions of KRH1 and its homolog KRH2 were hyper-invasive and displayed a high level of expression of FLO11, a gene involved in pseudohyphal and invasive growth. Therefore, KRH1 and KRH2 encode negative regulators of the invasive growth pathway. Cells containing krh1Delta krh2Delta mutations also displayed increased sensitivity to heat shock and decreased sporulation efficiency, indicating that Krh1p and Krh2p regulate multiple processes controlled by the cAMP/PKA pathway. The krh1Delta krh2Delta mutations suppressed the effect of a gpa2Delta mutation on FLO11 expression and eliminated the effect of a constitutively active GPA2 allele on induction of FLO11 and heat shock sensitivity, suggesting that Krh1p and Krh2p act downstream of Gpa2p. The Sch9p kinase was not required for the signal generated by deletion of KRH1 and KRH2; however, the cAMP-dependent kinase Tpk2p was required for generation of this signal. These results support a model in which activation of Gpa2p relieves the inhibition exerted by Krh1p and Krh2p on components of the cAMP/PKA signaling pathway.

Trehalose-mediated inhibition of the plasma membrane H+-ATPase from Kluyveromyces lactis: dependence on viscosity and temperature

The effect of increasing trehalose concentrations on the kinetics of the plasma membrane H+-ATPase from Kluyveromyces lactis was studied at different temperatures. At 20 degrees C, increasing concentrations of trehalose (0.2 to 0.8 M) decreased V(max) and increased S(0.5) (substrate concentration when initial velocity equals 0.5 V(max)), mainly at high trehalose concentrations (0.6 to 0.8 M). The quotient V(max)/S(0.5) decreased from 5.76 micromol of ATP mg of protein(-1) x min(-1) x mM(-1) in the absence of trehalose to 1.63 micromol of ATP mg of protein(-1) x min(-1) x mM(-1) in the presence of 0.8 M trehalose. The decrease in V(max) was linearly dependent on solution viscosity (eta), suggesting that inhibition was due to hindering of protein domain diffusional motion during catalysis and in accordance with Kramer's theory for reactions in solution. In this regard, two other viscosity-increasing agents, sucrose and glycerol, behaved similarly, exhibiting the same viscosity-enzyme inhibition correlation predicted. In the absence of trehalose, increasing the temperature up to 40 degrees C resulted in an exponential increase in V(max) and a decrease in enzyme cooperativity (n), while S(0.5) was not modified. As temperature increased, the effect of trehalose on V(max) decreased to become negligible at 40 degrees C, in good correlation with the temperature-mediated decrease in viscosity. The trehalose-mediated increase in S(0.5) was similar at all temperatures tested, and thus, trehalose effects on V(max)/S(0.5) were always observed. Trehalose increased the activation energy for ATP hydrolysis. Trehalose-mediated inhibition of enzymes may explain why yeast rapidly hydrolyzes trehalose when exiting heat shock.

Gene order evolution and paleopolyploidy in hemiascomycete yeasts

The wealth of comparative genomics data from yeast species allows the molecular evolution of these eukaryotes to be studied in great detail. We used "proximity plots" to visually compare chromosomal gene order information from 14 hemiascomycetes, including the recent Génolevures survey, to Saccharomyces cerevisiae. Contrary to the original reports, we find that the Génolevures data strongly support the hypothesis that S. cerevisiae is a degenerate polyploid. Using gene order information alone, 70% of the S. cerevisiae genome can be mapped into "sister" regions that tile together with almost no overlap. This map confirms and extends the map of sister regions that we constructed previously by using duplicated genes, an independent source of information. Combining gene order and gene duplication data assigns essentially the whole genome into sister regions, the largest gap being only 36 genes long. The 16 centromere regions of S. cerevisiae form eight pairs, indicating that an ancestor with eight chromosomes underwent complete doubling; alternatives such as segmental duplications can be ruled out. Gene arrangements in Kluyveromyces lactis and four other species agree quantitatively with what would be expected if they diverged from S. cerevisiae before its polyploidization. In contrast, Saccharomyces exiguus, Saccharomyces servazzii, and Candida glabrata show higher levels of gene adjacency conservation, and more cases of imperfect conservation, suggesting that they split from the S. cerevisiae lineage after polyploidization. This finding is confirmed by sequences around the C. glabrata TRP1 and IPP1 loci, which show that it contains sister regions derived from the same duplication event as that of S. cerevisiae.

A genetic map of Blumeria graminis based on functional genes, avirulence genes, and molecular markers

A genetic map of the powdery mildew fungus, Blumeria graminis f. sp. hordei, an obligate biotrophic pathogen of barley, is presented. The linkage analysis was conducted on 81 segregating haploid progeny isolates from a cross between 2 isolates differing in seven avirulence genes. A total of 359 loci were mapped, comprising 182 amplified fragment length polymorphism markers, 168 restriction fragment length polymorphism markers including 42 LTR-retrotransposon loci and 99 expressed sequence tags (ESTs), all the seven avirulence genes, and a marker closely linked to the mating type gene. The markers are distributed over 34 linkage groups covering a total of 2114 cM. Five avirulence genes were found to be linked and mapped in clusters of three and two, and two were unlinked. The Avr(a6) gene was found to be closely linked to markers suitable for a map-based cloning approach. A linkage between ESTs allowed us to demonstrate examples of synteny between genes in B. graminis and Neurospora crassa.

Fungi and humans: closer than you think

The budding yeast, Saccharomyces cerevisiae, has long been used as a model system to study the functions of human genes. Now that the genome sequences from several other fungal species are nearly complete, we can characterize the genetic diversity in the fungal kingdom at the genomic level. This diversity means that the number of human genes with homologues in the fungal kingdom is double that with homologues in S. cerevisiae only. Therefore, functional studies of human genes in the fungal model systems should look beyond S. cerevisiae.

The sequence of a 15 769 bp segment of Pichia anomala identifies the SEC61 and FBP1 genes and five new open reading frames

We have determined the sequence of a 15 769 bp DNA segment of Pichia anomala. The sequence contains seven complete open reading frames (ORFs) longer than 100 amino acids and a putative tRNA gene. Two of the ORFs code for the well-characterized genes SEC61 (which codes for the core subunit of the ER translocation complex) and FBP1 (encoding fructose-1,6-bisphosphatase). A gene coding for a protein similar to S. cerevisiae YDL054c was found between the two genes. These three genes show a different organization (intermingled triples) in three yeast species: Saccharomyces cerevisiae, Candida albicans and P. anomala. Two out of the four remaining ORFs show weak homology with different proteins from other species and the other two show non-significant similarity with previously sequenced genes. The nucleotide sequence has been submitted to the EMBL database under Accession No. AJ306295.

Surveying Saccharomyces genomes to identify functional elements by comparative DNA sequence analysis

Comparative sequence analysis has facilitated the discovery of protein coding genes and important functional sequences within proteins, but has been less useful for identifying functional sequence elements in nonprotein-coding DNA because the relatively rapid rate of change of nonprotein-coding sequences and the relative simplicity of non-coding regulatory sequence elements necessitates the comparison of sequences of relatively closely related species. We tested the use of comparative DNA sequence analysis to aid identification of promoter regulatory elements, nonprotein-coding RNA genes, and small protein-coding genes by surveying random DNA sequences of several Saccharomyces yeast species, with the goal of learning which species are best suited for comparisons with S. cerevisiae. We also determined the DNA sequence of a few specific promoters and RNA genes of several Saccharomyces species to determine the degree of conservation of known functional elements within the genome. Our results lead us to conclude that comparative DNA sequence analysis will enable identification of functionally conserved elements within the yeast genome, and suggest a path for obtaining this information.

Regions of microsynteny in Magnaporthe grisea and Neurospora crassa

A bacterial artificial chromosome (BAC) clone containing 110,467 bp of genomic DNA from Magnaporthe grisea was sequenced, annotated, and compared to the genomes of Neurospora crassa, Candida albicans, and Saccharomyces cerevisiae. Twenty-six open reading frames (ORFs), involved in multiple biochemical pathways, were identified in the BAC sequence. A region of 53 kb, containing 18 of the 26 ORFs, was found to be syntenic to a portion of the N. crassa genome. Subregions of complete colinearity as well as interrupted colinearity were present. No synteny was evident with either C. albicans or S. cerevisiae. The identification of syntenic regions containing highly conserved genes across two genera that have been evolutionarily separated for approximately 200 million years elicits many biological questions as to the function and identity of these genes.

Gene order in a 10 275 bp fragment of Yarrowia lipolytica, including adjacent YlURA5 and YlSEC65 genes conserved in four yeast species

We have determined the sequence of a 10275 bp DNA segment of Yarrowia lipolytica located on chromosome VI. The sequence contains six complete open reading frames (ORFs) longer than 100 amino acids and two more partial ORFs at both ends. Two of the ORFs encode for the well-characterized genes YlURA5 (orotate phosphoribosyltransferase) and YlSEC65 (encoding a subunit of the signal recognition particle). These two genes show an identical organization-located on opposite strands and in opposite orientations-in four yeast species: Saccharomyces cerevisiae, Kluyveromyces lactis, Candida albicans and Y. lipolytica. One ORF and the two partial ORFs code for putative proteins showing significant homology with proteins from other organisms. YlVI-108w (partial) and YlVI-103w show 39% and 54% identity, respectively, with YDR430c and YHR088w from S. cerevisiae. YlVI-102c (partial) shows significant homology with a matrix protein, lustrin A from Haliotis rufescens, and with the PGRS subfamily (Gly-rich proteins) of Mycobacterium tuberculosis. The three remaining ORFs show weak or non-significant homology with previously sequenced genes. The nucleotide sequence has been submitted to the EMBL database under Accession No. AI006754.

Uncoupling of 3'-phosphatase and 5'-kinase functions in budding yeast. Characterization of Saccharomyces cerevisiae DNA 3'-phosphatase (TPP1)

Polynucleotide kinase is a bifunctional enzyme containing both DNA 3'-phosphatase and 5'-kinase activities seemingly suited to the coupled repair of single-strand nicks in which the phosphate has remained with the 3'-base. We show that the yeast Saccharomyces cerevisiae is able to repair transformed dephosphorylated linear plasmids by non-homologous end joining with considerable efficiency independently of the end-processing polymerase Pol4p. Homology searches and biochemical assays did not reveal a 5'-kinase that would account for this repair, however. Instead, open reading frame YMR156C (here named TPP1) is shown to encode only a polynucleotide kinase-type 3'-phosphatase. Tpp1p bears extensive similarity to the ancient L-2-halo-acid dehalogenase and DDDD phosphohydrolase superfamilies, but is specific for double-stranded DNA. It is present at high levels in cell extracts in a functional form and so does not represent a pseudogene. Moreover, the phosphatase-only nature of this gene is shared by Saccharomyces mikatae YMR156C and Arabidopsis thaliana K15M2.3. Repair of 3'-phosphate and 5'-hydroxyl lesions is thus uncoupled in budding yeast as compared with metazoans. Repair of transformed dephosphorylated plasmids, and 5'-hydroxyl blocking lesions more generally, likely proceeds by a cycle of base removal and resynthesis.

Cloning and sequence of the LYS2 homologue gene from the osmotolerant yeast Pichia sorbitophila

We have isolated the Pichia sorbitophila LYS2 (PsLYS2) gene by complementation of a lys2 Saccharomyces cerevisiae mutant. The sequenced DNA fragment contains a putative ORF of 4197 bp and the deduced translation product shares a global identity of 66% and 58% to the Lys2 protein homologues of Candida albicans and S. cerevisiae, respectively. Analysis of PsLYS2 sequence suggests that, similarly to S. cerevisiae LYS2, it codes for a polypeptide having two separate enzymatic activities which reside in different domains of the protein, including an adenylate domain, an acyl-carrier site and a short-chain reductase domain. Several GCN4- and NIT2-binding motifs have been matched in the promotor sequence of PsLYS2. In addition, upstream of the sequenced PsLYS2 sequence, we have found the 3'-terminal half of a gene of same orientation encoding a RAD16-like protein, a genomic organization similar to that of C. albicans.

Genomic exploration of the hemiascomycetous yeasts: 8. Zygosaccharomyces rouxii

This paper reports the genomic analysis of strain CBS732 of Zygosaccharomyces rouxii, a homothallic diploid yeast. We explored the sequences of 4934 random sequencing tags of about 1 kb in size and compared them to the Saccharomyces cerevisiae gene products. Approximately 2250 nuclear genes, 57 tRNAs, the rDNA locus, the endogenous pSR1 plasmid and 15 mitochondrial genes were identified. According to 18S and 25S rRNA cladograms and to synteny analysis, Z. rouxii could be placed among the S. cerevisiae sensu lato yeasts.

Genomic exploration of the hemiascomycetous yeasts: 13. Pichia angusta

As part of a comparative genomics project on 13 hemiascomycetous yeasts, the Pichia angusta type strain was studied using a partial random sequencing strategy. With coverage of 0.5 genome equivalents, about 2500 novel protein-coding genes were identified, probably corresponding to more than half of the P. angusta protein-coding genes, 6% of which do not have homologs in Saccharomyces cerevisiae. Some of them contain one or two introns, on average three times shorter than those in S. cerevisiae. We also identified 28 tRNA genes, a few retrotransposons similar to Ty5 of S. cerevisiae, solo long terminal repeats, the whole ribosomal DNA cluster, and segments of mitochondrial DNA. The P. angusta sequences were deposited in EMBL under the accession numbers AL430961 to AL436044.

Genomic exploration of the hemiascomycetous yeasts: 15. Pichia sorbitophila

This paper reports the genomic analysis of the strain CBS7064 of Pichia sorbitophila, a homothallic diploid yeast. We sequenced 4829 random sequence tags of about 1 kb and compared them to the Saccharomyces cerevisiae gene products. Approximately 1300 nuclear genes, 22 tRNAs, the rDNA locus, and six mitochondrial genes have been identified. The analysis of the rDNA genes has permitted to classify this organism close to the Candida species. Accession numbers from AL414896 to AL419724 at EMBL databank.

Genomic exploration of the hemiascomycetous yeasts: 1. A set of yeast species for molecular evolution studies

The identification of molecular evolutionary mechanisms in eukaryotes is approached by a comparative genomics study of a homogeneous group of species classified as Hemiascomycetes. This group includes Saccharomyces cerevisiae, the first eukaryotic genome entirely sequenced, back in 1996. A random sequencing analysis has been performed on 13 different species sharing a small genome size and a low frequency of introns. Detailed information is provided in the 20 following papers. Additional tables available on websites describe the ca. 20000 newly identified genes. This wealth of data, so far unique among eukaryotes, allowed us to examine the conservation of chromosome maps, to identify the 'yeast-specific' genes, and to review the distribution of gene families into functional classes. This project conducted by a network of seven French laboratories has been designated 'Génolevures'.

Prevalence of small inversions in yeast gene order evolution

Gene order evolution in two eukaryotes was studied by comparing the Saccharomyces cerevisiae genome sequence to extensive new data from whole-genome shotgun and cosmid sequencing of Candida albicans. Gene order is substantially different between these two yeasts, with only 9% of gene pairs that are adjacent in one species being conserved as adjacent in the other. Inversion of small segments of DNA, less than 10 genes long, has been a major cause of rearrangement, which means that even where a pair of genes has been conserved as adjacent, the transcriptional orientations of the two genes relative to one another are often different. We estimate that about 1,100 single-gene inversions have occurred since the divergence between these species. Other genes that are adjacent in one species are in the same neighborhood in the other, but their precise arrangement has been disrupted, probably by multiple successive multigene inversions. We estimate that gene adjacencies have been broken as frequently by local rearrangements as by chromosomal translocations or long-distance transpositions. A bias toward small inversions has been suggested by other studies on animals and plants and may be general among eukaryotes.

Yeast chromosomes have been significantly reshaped during their evolutionary history

The structure of the first eukaryotic genome, belonging to Saccharomyces cerevisiae, has been deduced; however, very little is known about its origin. In order to trace events that led to the current state of the Saccharomyces nuclear genomes, random fragments of genomic DNA from three yeasts were sequenced and compared to the S. cerevisiae database sequence. Whereas, S. cerevisiae and Saccharomyces bayanus show perfect synteny, a significant portion of the analysed fragments from Saccharomyces servazzii and Saccharomyces kluyveri show a different arrangement of genes when compared to S. cerevisiae. When the sequenced fragments were probed to the corresponding karyotype, a group of genes present on a single chromosome of S. servazzii and S. kluyveri had homologues scattered on several S. cerevisiae chromosomes. Apparently, extensive reorganisation of the chromosomes has taken place during evolution of the Saccharomyces yeasts. In addition, while one gross duplication could have taken place, at least a few genes have been duplicated independently at different time-points in the evolution.

Organization of specific genomic regions of Zygosaccharomyces rouxii and Pichia sorbitophila: comparison with Saccharomyces cerevisiae

The genomes of Zygosaccharomyces rouxii and Pichia sorbitophila were partially explored. The genome of Z. rouxii CBS 732 consists of seven chromosomes with an approximate size of 1.0-2.75 Mb, 12.8 Mb in total. Five of the chromosomes were labelled with specific probes. Three Z. rouxii genomic DNA fragments were sequenced; all 10 ORFs found were without introns and they have homologues in S. cerevisiae. Gene order comparison revealed that the organization is partially conserved in both species. The genome of P. sorbitophila CBS 7064 consists of seven chromosomes with an approximate size of 1.0-2.9 Mb, 13.9 Mb in total. Three of the chromosomes were labelled with specific probes. The sequencing of a 5.2 kb genomic DNA fragment revealed three ORFs, but no conservation of their organization was found, although all of them have their respective homologues in S. cerevisiae. According to our results, the presence of two overlapping ORFs in S. cerevisiae (YJL107c-YJL108c) could be interpreted as the result of a frameshift mutation.

Prm1p, a pheromone-regulated multispanning membrane protein, facilitates plasma membrane fusion during yeast mating

Cell fusion occurs throughout development, from fertilization to organogenesis. The molecular mechanisms driving plasma membrane fusion in these processes remain unknown. While yeast mating offers an excellent model system in which to study cell fusion, all genes previously shown to regulate the process act at or before cell wall breakdown; i.e., well before the two plasma membranes have come in contact. Using a new strategy in which genomic data is used to predict which genes may possess a given function, we identified PRM1, a gene that is selectively expressed during mating and that encodes a multispanning transmembrane protein. Prm1p localizes to sites of cell-cell contact where fusion occurs. In matings between Deltaprm1 mutants, a large fraction of cells initiate zygote formation and degrade the cell wall separating mating partners but then fail to fuse. Electron microscopic analysis reveals that the two plasma membranes in these mating pairs are tightly apposed, remaining separated only by a uniform gap of approximately 8 nm. Thus, the phenotype of Deltaprm1 mutants defines a new step in the mating reaction in which membranes are juxtaposed, possibly through a defined adherence junction, yet remain unfused. This phenotype suggests a role for Prm1p in plasma membrane fusion.

Kluyveromyces marxianus exhibits an ancestral Saccharomyces cerevisiae genome organization downstream of ADH2

In Saccharomyces cerevisiae, the alcohol dehydrogenase genes ADH1 and ADH5 are part of a duplicated block of genome, thought to originate from a genome-wide duplication posterior to the divergence from the Kluyveromyces lineage. We report here the characterization of Kluyveromyces marxianus ADH2 and the five genes found in its immediate downstream region, MRPS9, YOL087C, RPB5, RIB7 and SPP381. The order of these six genes reflects the structure of the ancestral S. cerevisiae genome before the duplication that formed the blocks including ADH1 on chromosome XV and ADH5 on chromosome II, indicating these ADH genes share a direct ancestor. On the one hand, the two genes found immediately downstream of KmADH2 are located, for the first, downstream ADH5 and, for the second, downstream ADH1 in S. cerevisiae. On the other hand, the order of the paralogs included in the blocks of ADH1 and ADH5 in S. cerevisiae suggests that two of them have been inverted within one block after its formation, and that inversion is confirmed by the gene order observed in K. marxianus.

Genetics and molecular physiology of the yeast Kluyveromyces lactis

With the recent development of powerful molecular genetic tools, Kluyveromyces lactis has become an excellent alternative yeast model organism for studying the relationships between genetics and physiology. In particular, comparative yeast research has been providing insights into the strikingly different physiological strategies that are reflected by dominance of respiration over fermentation in K. lactis versus Saccharomyces cerevisiae. Other than S. cerevisiae, whose physiology is exceptionally affected by the so-called glucose effect, K. lactis is adapted to aerobiosis and its respiratory system does not underlie glucose repression. As a consequence, K. lactis has been successfully established in biomass-directed industrial applications and large-scale expression of biotechnically relevant gene products. In addition, K. lactis maintains species-specific phenomena such as the "DNA-killer system, " analyses of which are promising to extend our knowledge about microbial competition and the fundamentals of plasmid biology.

Comparative genetic and physiological studies of the MAP kinase Mpk1p from Kluyveromyces lactis and Saccharomyces cerevisiae

MAP kinases are essential components of signal transduction pathways in yeasts and higher eukaryotes. Here, we report on the isolation of the gene encoding the MAP kinase KlMpk1p by complementation of the respective Saccharomyces cerevisiae deletion mutant with a genomic library from Kluyveromyces lactis. Sequencing revealed the presence of an open reading frame capable of encoding a protein of 520 amino acid residues with a deduced molecular mass of 59.726 Da. The deduced protein sequence displayed a high degree of similarity to known MAP kinases from yeast to man, with an overall identity of 70 % to ScMpk1p. One-hybrid analysis demonstrated the presence of a cryptic transcriptional activation domain in the C-terminal part of the protein. Deletion of this sequence in ScMpk1p resulted in a reduced MAP kinase activity (measured by an indirect assay), an increased sensitivity towards caffeine and an increased resistance against Calcofluor white. Complete deletion mutants of Klmpk1 display an osmo-remedial phenotype on rich medium, but are capable of growth in the absence of osmotic stabilization on synthetic medium. As Scmpk1 deletion mutants, they are sensitive to cell surface destabilizing agents such as Calcofluor white and SDS, and growth is inhibited in the presence of 5 mM caffeine. Overexpression of KlMPK1 did not produce a growth defect in S. cerevisiae or in K. lactis.

The conserved genome organisation of non-falciparum malaria species: the need to know more

The current knowledge on genomes of non-falciparum malaria species and the potential of model malaria parasites for functional analyses are reviewed and compared with those of the most pathogenic human parasite, Plasmodium falciparum. There are remarkable similarities in overall genome composition among the different species at the level of chromosome organisation and chromosome number, conserved order of individual genes, and even conserved functions of specific gene domains and regulatory control elements. With the initiative taken to sequence the genome of P. falciparum, a wealth of information is already becoming available to the scientific community. In order to exploit the biological information content of a complete genome sequence, simple storage of the bulk of sequence data will be inadequate. The requirement for functional analyses to determine the biological role of the open reading frames is commonly accepted and knowledge of the genomes of the animal model malaria species will facilitate these analyses. Detailed comparative genome information and sequencing of additional Plasmodium genomes will provide a deeper insight into the evolutionary history of the species, the biology of the parasite, and its interactions with the mammalian host and mosquito vector. Therefore, an extended and integrated approach will enhance our knowledge of malaria and will ultimately lead to a more rational approach that identifies and evaluates new targets for anti-malarial drug and vaccine development.

Large-scale comparison of fungal sequence information: mechanisms of innovation in Neurospora crassa and gene loss in Saccharomyces cerevisiae

We report a large-scale comparison of sequence data from the filamentous fungus Neurospora crassa with the complete genome sequence of Saccharomyces cerevisiae. N. crassa is considerably more morphologically and developmentally complex than S. cerevisiae. We found that N. crassa has a much higher proportion of "orphan" genes than S. cerevisiae, suggesting that its morphological complexity reflects the acquisition or maintenance of novel genes, consistent with its larger genome. Our results also indicate the loss of specific genes from S. cerevisiae. Surprisingly, some of the genes lost from S. cerevisiae are involved in basic cellular processes, including translation and ion (especially calcium) homeostasis. Horizontal gene transfer from prokaryotes appears to have played a relatively modest role in the evolution of the N. crassa genome. Differences in the overall rate of molecular evolution between N. crassa and S. cerevisiae were not detected. Our results indicate that the current public sequence databases have fairly complete samples of gene families with ancient conserved regions, suggesting that further sequencing will not substantially change the proportion of genes with homologs among distantly related groups. Models of the evolution of fungal genomes compatible with these results, and their functional implications, are discussed.

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.

Yeast genome evolution in the post-genome era

The Saccharomyces cerevisiae genome sequence, augmented by new data on gene expression and function, continues to yield new findings about eukaryote genome evolution. Analysis of the duplicate gene pairs formed by whole-genome duplication indicates that selection for increased levels of gene expression was a significant factor determining which genes were retained as duplicates and which were returned to a single-copy state, possibly in addition to selection for novel gene functions. Proteome comparisons between worm and yeast show that genes for core metabolic processes are shared among eukaryotes and unchanging in function, while comparisons between different yeast species identify 'orphan' genes as the most rapidly evolving fraction of the proteome. Natural hybridisation among yeast species is frequent, but its long-term evolutionary significance is unknown.

Updated map of duplicated regions in the yeast genome

We have updated the map of duplicated chromosomal segments in the Saccharomyces cerevisiae genome originally published by Wolfe and Shields in 1997 (Nature 387, 708-713). The new analysis is based on the more sensitive Smith Waterman search method instead of BLAST. The parameters used to identify duplicated chromosomal regions were optimized such as to maximize the amount of the genome placed into paired regions, under the assumption that the hypothesis that the entire genome was duplicated in a single event is correct. The core of the new map, with 52 pairs of regions containing three or more duplicated genes, is largely unchanged from our original map. 39 tRNA gene pairs and one snRNA pair have been added. To find additional pairs of genes that may have been formed by whole genome duplication, we searched through the parts of the genome that are not covered by this core map, looking for putative duplicated chromosomal regions containing only two duplicate genes instead of three, or having lower-scoring gene pairs. This approach identified a further 32 candidate paired regions, bringing the total number of protein-coding genes on the duplication map to 905 (16% of the proteome). The updated map suggests that a second copy of the ribosomal DNA array has been deleted from chromosome IV.