Transcriptome analysis reveals multiple targets of erythritol-related transcription factor EUF1 in unconventional yeast Yarrowia Lipolytica.

BACKGROUND: Erythritol is a four-carbon polyol with an unclear role in metabolism of some unconventional yeasts. Its production has been linked to the osmotic stress response, but the mechanism of stress protection remains unclear. Additionally, erythritol can be used as a carbon source. In the yeast Yarrowia lipolytica, its assimilation is activated by the transcription factor Euf1. The study investigates whether this factor can link erythritol to other processes in the cell. RESULTS: The research was performed on two closely related strains of Y. lipolytica: MK1 and K1, where strain K1 has no functional Euf1. Cultures were carried out in erythritol-containing and erythritol-free media. Transcriptome analysis revealed the effect of Euf1 on the regulation of more than 150 genes. Some of these could be easily connected with different aspects of erythritol assimilation, such as: utilization pathway, a new potential isoform of transketolase, or polyol transporters. However, many of the upregulated genes have never been linked to metabolism of erythritol. The most prominent examples are the degradation pathway of branched-chain amino acids and the glyoxylate cycle. The high transcription of genes affected by Euf1 is still dependent on the erythritol concentration in the medium. Moreover, almost all up-regulated genes have an ATGCA motif in the promoter sequence. CONCLUSIONS: These findings may be particularly relevant given the increasing use of erythritol-induced promoters in genetic engineering of Y. lipolytica. Moreover, use of this yeast in biotechnological processes often takes place under osmotic stress conditions. Erythritol might be produce as a by-product, thus better understanding of its influence on cell metabolism could facilitate processes optimization.

A shuttle mutagenesis system for tagging genes in the yeast Yarrowia lipolytica.

A shuttle mutagenesis system was developed for the dimorphic yeast Yarrowia lipolytica. This system combines transposon insertions generated in Escherichia coli with the transformation of yeast with the Tn-mutagenized DNA. The mini-transposon mTn-3xHA/GFP, used in Saccharomyces cerevisiae for producing stable insertions, was adapted for use in the yeast Y. lipolytica. The mTnYl1 transposon (for mini-Tn of Y. lipolytica) confers resistance to tetracycline in E. coli. It also contains the Y. lipolytica URA3 gene for selection of yeast transformants, and the coding sequence for the S65T mutant form of GFP. The rare cutter endonuclease, I-SceI, restriction site, which enables identification of the chromosomal localization of mutagenized genes, was also incorporated. mTnYl1 was first tested on the ACO1 gene, which encodes an Acyl CoA oxidase isozyme. The mutagenesis system was further validated on a Y. lipolytica genomic DNA library constructed in a pHSS6 derivative vector. Mutants with a particular morphology or defective for alkane, fatty acids and oil degradation were obtained.

Genomic exploration of the hemiascomycetous yeasts: 5. Saccharomyces bayanus var. uvarum.

Saccharomyces bayanus var. uvarum investigated here is the species closest to Saccharomyces cerevisiae. Random sequence tags (RSTs) allowed us to identify homologues to 2789 open reading frames (ORFs) in S. cerevisiae, ORFs duplicated in S. uvarum but not in S. cerevisiae, centromeres, tRNAs, homologues of Ty1/2 and Ty4 retrotransposons, and a complete rDNA repeat. Only 13 RSTs seem to be homologous to sequences in other organisms but not in S. cerevisiae. As the synteny between the two species is very high, cases in which synteny is lost suggest special mechanisms of genome evolution. The corresponding RSTs revealed that S. uvarum can exist without any S. cerevisiae DNA introgression. Accession numbers are from AL397139 to AL402278 in the EMBL databank.

Genomic exploration of the hemiascomycetous yeasts: 6. Saccharomyces exiguus.

Random sequence tags were obtained from a genomic DNA library of Saccharomyces exiguus. The mitochondrial genome appeared to be at least 25.7 kb in size, with a different organization compared to Saccharomyces cerevisiae. An unusual putative 953 bp long terminal repeated element associated to Ty3 was found. A set of 1451 genes was identified homologous to S. cerevisiae open reading frames. Only five genes were identified outside the S. cerevisiae taxon, confirming that S. exiguus is phylogenetically closely related to S. cerevisiae. Unexpectedly, numerous duplicated genes were found whereas they are unique in S. cerevisiae. The sequences are deposited at EMBL under the accession numbers: AL407377-AL409955.

Genomic exploration of the hemiascomycetous yeasts: 7. Saccharomyces servazzii.

The genome of Saccharomyces servazzii was analyzed with 2570 random sequence tags totalling 2.3 Mb. BLASTX comparisons revealed a minimum of 1420 putative open reading frames with significant homology to Saccharomyces cerevisiae (58% aa identity on average), two with Schizosaccharomyces pombe and one with a human protein, confirming that S. servazzii is closely related to S. cerevisiae. About 25% of the S. servazzii genes were identified, assuming that the gene complement is identical in both yeasts. S. servazzii carries very few transposable elements related to Ty elements in S. cerevisiae. Most of the mitochondrial genes were identified in eight contigs altogether spanning 25 kb for a predicted size of 29 kb. A significant match with the Kluyveromyces lactis linear DNA plasmid pGKL-1 encoded RF4 killer protein suggests that a related plasmid exists in S. servazzii. The sequences have been deposited with EMBL under the accession numbers AL402279-AL404848.

Genomic exploration of the hemiascomycetous yeasts: 9. Saccharomyces kluyveri.

The genome of Saccharomyces kluyveri was explored through 2528 random sequence tags with an average length of 981 bp. The complete nuclear ribosomal DNA unit was found to be 8656 bp in length. Sequences homologous to retroelements of the gypsy and copia types were identified as well as numerous solo long terminal repeats. We identified at least 1406 genes homologous to Saccharomyces cerevisiae open reading frames, with on average 58.1% and 72.4% amino acid identity and similarity, respectively. In addition, by comparison with completely sequenced genomes and the SwissProt database, we found 27 novel S. kluyveri genes. Most of these genes belong to pathways or have functions absent from S. cerevisiae, such as the catabolic pathway of purines or pyrimidines, melibiose fermentation, sorbitol utilization, or degradation of pollutants. The sequences are deposited in EMBL under the accession numbers AL404849-AL407376.

Genomic exploration of the hemiascomycetous yeasts: 14. Debaryomyces hansenii var. hansenii.

By analyzing 2830 random sequence tags (RSTs), totalling 2.7 Mb, we explored the genome of the marine, osmo- and halotolerant yeast, Debaryomyces hansenii. A contig 29 kb in length harbors the entire mitochondrial genome. The genes encoding Cox1, Cox2, Cox3, Cob, Atp6, Atp8, Atp9, several subunits of the NADH dehydrogenase complex 1 and 11 tRNAs were unambiguously identified. An equivalent number of putative transposable elements compared to Saccharomyces cerevisiae were detected, the majority of which are more related to higher eukaryote copia elements. BLASTX comparisons of RSTs with databases revealed at least 1119 putative open reading frames with homology to S. cerevisiae and 49 to other genomes. Specific functions, including transport of metabolites, are clearly over-represented in D. hansenii compared to S. cerevisiae, consistent with the observed difference in physiology of the two species. The sequences have been deposited with EMBL under the accession numbers AL436045-AL438874.

Genomic exploration of the hemiascomycetous yeasts: 17. Yarrowia lipolytica.

A total of 4940 random sequence tags of the dimorphic yeast Yarrowia lipolytica, totalling 4.9 Mb, were analyzed. BLASTX comparisons revealed at least 1229 novel Y. lipolytica genes 1083 genes having homology with Saccharomyces cerevisiae genes and 146 with genes from various other genomes. This confirms the rapid sequence evolution assumed for Y. lipolytica. Functional analysis of newly discovered genes revealed that several enzymatic activities were increased compared to S. cerevisiae, in particular, transport activities, ion homeostasis, and various metabolism pathways. Most of the mitochondrial genes were identified in contigs spanning more than 47 kb. Matches to retrotransposons were observed, including a S. cerevisiae Ty3 and a LINE element. The sequences have been deposited with EMBL under the accession numbers AL409956-AL414895.

Genomic exploration of the hemiascomycetous yeasts: 3. Methods and strategies used for sequence analysis and annotation.

The primary analysis of the sequences for our Hemiascomycete random sequence tag (RST) project was performed using a combination of classical methods for sequence comparison and contig assembly, and of specifically written scripts and computer visualization routines. Comparisons were performed first against DNA and protein sequences from Saccharomyces cerevisiae, then against protein sequences from other completely sequenced organisms and, finally, against protein sequences from all other organisms. Blast alignments were individually inspected to help recognize genes within our random genomic sequences despite the fact that only parts of them were available. For each yeast species, validated alignments were used to infer the proper genetic code, to determine codon usage preferences and to calculate their degree of sequence divergence with S. cerevisiae. The quality of each genomic library was monitored from contig analysis of the DNA sequences. Annotated sequences were submitted to the EMBL database, and the general annotation tables produced served as a basis for our comparative description of the evolution, redundancy and function of the Hemiascomycete genomes described in other articles of this issue.

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'.

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

Since its completion more than 4 years ago, the sequence of Saccharomyces cerevisiae has been extensively used and studied. The original sequence has received a few corrections, and the identification of genes has been completed, thanks in particular to transcriptome analyses and to specialized studies on introns, tRNA genes, transposons or multigene families. In order to undertake the extensive comparative sequence analysis of this program, we have entirely revisited the S. cerevisiae sequence using the same criteria for all 16 chromosomes and taking into account publicly available annotations for genes and elements that cannot be predicted. Comparison with the other yeast species of this program indicates the existence of 50 novel genes in segments previously considered as 'intergenic' and suggests extensions for 26 of the previously annotated genes.

Genomic exploration of the hemiascomycetous yeasts: 18. Comparative analysis of chromosome maps and synteny with Saccharomyces cerevisiae.

We have analyzed the evolution of chromosome maps of Hemiascomycetes by comparing gene order and orientation of the 13 yeast species partially sequenced in this program with the genome map of Saccharomyces cerevisiae. From the analysis of nearly 8000 situations in which two distinct genes having homologs in S. cerevisiae could be identified on the sequenced inserts of another yeast species, we have quantified the loss of synteny, the frequency of single gene deletion and the occurrence of gene inversion. Traces of ancestral duplications in the genome of S. cerevisiae could be identified from the comparison with the other species that do not entirely coincide with those identified from the comparison of S. cerevisiae with itself. From such duplications and from the correlation observed between gene inversion and loss of synteny, a model is proposed for the molecular evolution of Hemiascomycetes. This model, which can possibly be extended to other eukaryotes, is based on the reiteration of events of duplication of chromosome segments, creating transient merodiploids that are subsequently resolved by single gene deletion events.

Genomic exploration of the hemiascomycetous yeasts: 19. Ascomycetes-specific genes.

Comparisons of the 6213 predicted Saccharomyces cerevisiae open reading frame (ORF) products with sequences from organisms of other biological phyla differentiate genes commonly conserved in evolution from 'maverick' genes which have no homologue in phyla other than the Ascomycetes. We show that a majority of the 'maverick' genes have homologues among other yeast species and thus define a set of 1892 genes that, from sequence comparisons, appear 'Ascomycetes-specific'. We estimate, retrospectively, that the S. cerevisiae genome contains 5651 actual protein-coding genes, 50 of which were identified for the first time in this work, and that the present public databases contain 612 predicted ORFs that are not real genes. Interestingly, the sequences of the 'Ascomycetes-specific' genes tend to diverge more rapidly in evolution than that of other genes. Half of the 'Ascomycetes-specific' genes are functionally characterized in S. cerevisiae, and a few functional categories are over-represented in them.

Genomic exploration of the hemiascomycetous yeasts: 20. Evolution of gene redundancy compared to Saccharomyces cerevisiae.

We have evaluated the degree of gene redundancy in the nuclear genomes of 13 hemiascomycetous yeast species. Saccharomyces cerevisiae singletons and gene families appear generally conserved in these species as singletons and families of similar size, respectively. Variations of the number of homologues with respect to that expected affect from 7 to less than 24% of each genome. Since S. cerevisiae homologues represent the majority of the genes identified in the genomes studied, the overall degree of gene redundancy seems conserved across all species. This is best explained by a dynamic equilibrium resulting from numerous events of gene duplication and deletion rather than by a massive duplication event occurring in some lineages and not in others.

Genomic exploration of the hemiascomycetous yeasts: 21. Comparative functional classification of genes.

We explored the biological diversity of hemiascomycetous yeasts using a set of 22000 newly identified genes in 13 species through BLASTX searches. Genes without clear homologue in Saccharomyces cerevisiae appeared to be conserved in several species, suggesting that they were recently lost by S. cerevisiae. They often identified well-known species-specific traits. Cases of gene acquisition through horizontal transfer appeared to occur very rarely if at all. All identified genes were ascribed to functional classes. Functional classes were differently represented among species. Species classification by functional clustering roughly paralleled rDNA phylogeny. Unequal distribution of rapidly evolving, ascomycete-specific, genes among species and functions was shown to contribute strongly to this clustering. A few cases of gene family amplification were documented, but no general correlation could be observed between functional differentiation of yeast species and variations of gene family sizes. Yeast biological diversity seems thus to result from limited species-specific gene losses or duplications, and for a large part from rapid evolution of genes and regulatory factors dedicated to specific functions.

Identification of group-I introns in the 28s rDNA of the entomopathogenic fungus Beauveria brongniartii.

The length of the 28s ribosomal DNA differs significantly between two strains (Bt102 and Bt114) of the entomopathogenic fungus Beauveria brongniartii. RFLP analysis on PCR products revealed the presence of three insertional elements of 350-450 bp in strain Bt114. One of the insertions has been cloned and sequenced and shown to possess all the characteristic sequences and secondary structures of a group-IC intron. Its length is 428 bp and it is devoid of any long open reading frame. The distribution of this intron elsewhere in the genome of Bt114, as well as in the chromosomal ribosomal DNA, was studied. It seems to be present as seven copies in different genes not corresponding to the mitochondrial DNA. The presence of the intron in other strains of B. brongniartii was examined by the hybridization method. Some of them seemed to possess introns with a similar core although others presented no homology with the cloned fragment.

28s rDNA group-I introns: a powerful tool for identifying strains of Beauveria brongniartii.

The nuclear ribosomal DNA of the entomopathogenic fungus Beauveria brongniartii is polymorphic in terms of both restriction site and length. Insertions of 350-450 bp long, identified as group-I introns, were detected in the 28s rDNA. A panel of 47 strains of B. brongniartii, two B. bassiana and one Metarhizium anisopliae of various geographical and biological origins were found to contain 14 variant forms of intron differing in size and restriction pattern, at four different positions. Twelve types of ribosomal large subunit were defined on the basis of variant distribution and compared with strain clustering based on internal transcribed spacers analysis. There was a correlation between the characteristic introns and isolates collected from the sugar cane pest Hoplochelus marginalis. Primers for polymerase chain reaction amplification were chosen from these variants, and used to develop a specific method for detecting strains pathogenic towards Hoplochelus.

Afut1, a retrotransposon-like element from Aspergillus fumigatus.

A repeated DNA sequence used for epidemiological studies of the human opportunistic pathogen Aspergillus fumigatus has been characterized. It is a retroelement of 6914 bp in length, bounded by long terminal repeats of 282 bp, with sequence and features characteristic of retroviruses and retrotransposons. A 5 bp duplication site was found at its borders. This element, designated Afut1, encodes amino acid sequences homologous to the reverse transcriptase, RNase H and endonuclease encoded by the pol genes of retroelements. Comparison of the peptidic sequences with other putative polypeptides of fungal LTR retrotransposons showed that Afut1 is a member of the gypsy group. This is the first report of a transposable element in A.fumigatus. Afut1 is a defective element: the putative coding domains contain multiple stop codons due exclusively to transitions from C:G to T:A.

Evolution of gene order in the genomes of two related yeast species.

Changes in gene order between the genomes of two related yeast species, Saccharomyces cerevisiae and Saccharomyces bayanus var. uvarum were studied. From the dataset of a previous low coverage sequencing of the S. bayanus var. uvarum genome, 35 different synteny breakpoints between neighboring genes and two cases of local gene inversion were characterized in detail. The number and the type of the chromosomal rearrangements that have led to these differences were identified. We show that evolution of gene order in the genomes of these two yeast species is driven mainly by gene duplication onto different chromosomes followed by differential loss of the repeated copies. In addition, local gene inversions also would result from a mechanism of gene duplication, but in an inverted orientation, followed by loss of the original copy. The identification of traces of anciently duplicated genes, called relics, show that the loss of duplicates is more frequently caused by the accumulation of numerous mutations in one of the two copies than by DNA deletion. Surprisingly, gross chromosomal rearrangements such as translocations have only a minor effect on gene order reshuffling as they account for <10% of the synteny breakpoints.

Homologous recombination and transposition generate chromosome I neopolymorphism during meiosis in Saccharomyces cerevisiae.

We have studied the meiotic segregation of a chromosome length polymorphism (CLP) in the yeast Saccharomyces cerevisiae. The neopolymorphism frequently observed within the smallest chromosomes (I, VI, III and IX) is not completely understood. We focused on the analysis of the structure of chromosome I in 88 segregants from a cross between YNN295 and FL100trp. Strain FL100trp is known to carry a reciprocal translocation between the left arm of chromosome III and the right arm of chromosome I. PCR and Southern hybridization analyses were performed and a method for the rapid detection of chromosome I rearrangements was developed. Seven chromosome I types were identified among the 88 segregants. We detected 22 recombination events between homologous chromosomes I and seven ectopic recombination events between FL100trp chromosome III and YNN295 chromosome I. These recombination events occurred in 20 of the 22 tetrads studied (91%). Nine tetrads (41%) showed two recombination events. This showed that homologous recombination involving polymorphic homologues or heterologous chromosomes is the main source of neopolymorphism. Only one of the seven chromosome I variants resulted from a transposition event rather than a recombination event. We demonstrated that a Tyl element had transposed within the translocated region of chromosome I, generating mutations in the 3' LTR, at the border between U5 and PBS.