Inheritance and organisation of the mitochondrial genome differ between two Saccharomyces yeasts

Petite-positive Saccharomyces yeasts can be roughly divided into the sensu stricto, including Saccharomyces cerevisiae, and sensu lato group, including Saccharomyces castellii; the latter was recently studied for transmission and the organisation of its mitochondrial genome. S. castellii mitochondrial molecules (mtDNA) carrying point mutations, which confer antibiotic resistance, behaved in genetic crosses as the corresponding point mutants of S. cerevisiae. While S. castellii generated spontaneous petite mutants in a similar way as S. cerevisiae, the petites exhibited a different inheritance pattern. In crosses with the wild type strains a majority of S. castellii petites was neutral, and the suppressivity in suppressive petites was never over 50%. The two yeasts also differ in organisation of their mtDNA molecules. The 25,753 bp sequence of S. castellii mtDNA was determined and the coding potential of both yeasts is similar. However, the S. castellii intergenic sequences are much shorter and do not contain sequences homologous to the S. cerevisiae biologically active intergenic sequences, as ori/rep/tra, which are responsible for the hyper-suppressive petite phenotype found in S. cerevisiae. The structure of one suppressive S. castellii mutant, CA38, was also determined. Apparently, a short direct intergenic repeat was involved in the generation of this petite mtDNA molecule.     

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.     

The intron of the mitochondrial 21S rRNA gene: distribution in different yeast species and sequence comparison between Kluyveromyces thermotolerans and Saccharomyces cerevisiae

We have screened numerous different yeast species for the presence of sequences homologous to the intron of the mitochondrial 21S rRNA gene of Saccharomyces cerevisiae (intron r1) and found them in all Kluyveromyces species, some of the Saccharomyces species and none of the other yeasts tested. We have determined the nucleotide sequence of the r1-intron in K. thermotolerans and compared it with that of S. cerevisiae. The two introns are inserted at the same position within the 21S rRNA gene. They contain homologous internal open reading frames (ORFs) initiated at the same AUG codon which can be aligned over their entire length. Several silent multi-substitutions indicate that these intronic ORFs represent selectively conserved functional genes. Other intron segments, on the contrary, reveal short blocks of extensive homology separated by non-homologous stretches and/or additions-deletions. Comparison of our two yeast r1-introns with equivalent introns of N. crassa and A. nidulans mitochondria reveals that introns with very similar RNA secondary structures can accommodate different types of ORFs.     

Diversity in organization and the origin of gene orders in the mitochondrial DNA molecules of the genus Saccharomyces

Sequencing of the Saccharomyces cerevisiae nuclear and mitochondrial genomes provided a new background for studies on the evolution of the genomes. In this study, mitochondrial genomes of a number of Saccharomyces yeasts were mapped by restriction enzyme analysis, the orders of the genes were determined, and two of the genes were sequenced. The genome organization, i.e., the size, presence of intergenic sequences, and gene order, as well as polymorphism within the coding regions, indicate that Saccharomyces mtDNA molecules are dynamic structures and have undergone numerous changes during their evolution. Since the separation and sexual isolation of different yeast lineages, the coding parts have been accumulating point mutations, presumably in a linear manner with the passage of time. However, the accumulation of other changes may not have been a simple function of time. Larger mtDNA molecules belonging to Saccharomyces sensu stricto yeasts have acquired extensive intergenic sequences, including guanosine-cytosine-rich clusters, and apparently have rearranged the gene order at higher rates than smaller mtDNAs belonging to the Saccharomyces sensu lato yeasts. While within the sensu stricto group transposition has been a predominant mechanism for the creation of novel gene orders, the sensu lato yeasts could have used both transposition- and inversion-based mechanisms.     

Order and orientation of genic sequences in circular mitochondrial DNA fromSaccharomyces exiguus: Implications for evolution of yeast mtDNAs

Mapping of sequences specifying the large and small ribosomal RNAs and six polypeptides in the circular 23.7 kbp mitochondrial DNA of Saccharomyces exiguus has shown that these genes have the same orientation and that a 5 gene cluster is common to this DNA and the 18.9 kbp mtDNA from Torulopsis glabrata. Included in the preserved region are juxtaposed sequences specifying ATPase subunits 6 and 9 which have the same order and orientation as analogous genes in the Escherichia coli unc operon. The above data, together with knowledge that these two sequences are dispersed in larger yeast mtDNAs, leads us to suggest that larger forms are derived from a smaller ancestral molecule that would have had some resemblance to the mtDNAs of S. exiguus and T. glabrata.     

After the duplication: gene loss and adaptation in Saccharomyces genomes

The ancient duplication of the Saccharomyces cerevisiae genome and subsequent massive loss of duplicated genes is apparent when it is compared to the genomes of related species that diverged before the duplication event. To learn more about the evolutionary effects of the duplication event, we compared the S. cerevisiae genome to other Saccharomyces genomes. We demonstrate that the whole genome duplication occurred before S. castellii diverged from S. cerevisiae. In addition to more accurately dating the duplication event, this finding allowed us to study the effects of the duplication on two separate lineages. Analyses of the duplication regions of the genomes indicate that most of the duplicated genes (approximately 85%) were lost before the speciation. Only a small amount of paralogous gene loss (4-6%) occurred after speciation. On the other hand, S. castellii appears to have lost several hundred genes that were not retained as duplicated paralogs. These losses could be related to genomic rearrangements that reduced the number of chromosomes from 16 to 9. In addition to S. castellii, other Saccharomyces sensu lato species likely diverged from S. cerevisiae after the duplication. A thorough analysis of these species will likely reveal other important outcomes of the whole genome duplication.     

A multispecies-based taxonomic microarray reveals interspecies hybridization and introgression in Saccharomyces cerevisiae

A multispecies-based taxonomic microarray targeting coding sequences of diverged orthologous genes in Saccharomyces cerevisiae, Saccharomyces paradoxus, Saccharomyces mikatae, Saccharomyces bayanus, Saccharomyces kudriavzevii, Naumovia castellii, Lachancea kluyveri and Candida glabrata was designed to allow identification of isolates of these species and their interspecies hybrids. Analysis of isolates of several Saccharomyces species and interspecies hybrids demonstrated the ability of the microarray to differentiate these yeasts on the basis of their specific hybridization patterns. Subsequent analysis of 183 supposed S. cerevisiae isolates of various ecological and geographical backgrounds revealed one misclassified S. bayanus or Saccharomyces uvarum isolate and four aneuploid interspecies hybrids, one between S. cerevisiae and S. bayanus and three between S. cerevisiae and S. kudriavzevii . Furthermore, this microarray design allowed the detection of multiple introgressed S. paradoxus DNA fragments in the genomes of three different S. cerevisiae isolates. These results show the power of multispecies-based microarrays as taxonomic tools for the identification of species and interspecies hybrids, and their ability to provide a more detailed characterization of interspecies hybrids and recombinants.     

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.     

Linear mitochondrial DNAs of yeasts: frequency of occurrence and general features.

In most yeast species, the mitochondrial DNA (mtDNA) has been reported to be a circular molecule. However, two cases of linear mtDNA with specific termini have previously been described. We examined the frequency of occurrence of linear forms of mtDNA among yeasts by pulsed-field gel electrophoresis. Among the 58 species from the genera Pichia and Williopsis that we examined, linear mtDNA was found with unexpectedly high frequency. Thirteen species contained a linear mtDNA, as confirmed by restriction mapping, and labeling, and electron microscopy. The mtDNAs from Pichia pijperi, Williopsis mrakii, and P. jadinii were studied in detail. In each case, the left and right terminal fragments shared homologous sequences. Between the terminal repeats, the order of mitochondrial genes was the same in all of the linear mtDNAs examined, despite a large variation of the genome size. This constancy of gene order is in contrast with the great variation of gene arrangement in circular mitochondrial genomes of yeasts. The coding sequences determined on several genes were highly homologous to those of the circular mtDNAs, suggesting that these two forms of mtDNA are not of distant origins.     

The mitochondrial genome of the fission yeast Schizosaccharomyces pombe. 2. Localization of genes by interspecific hybridization in strain ade7-50h- and cloning of the genome in small fragments

A series of 18 small overlapping restriction fragments has been cloned, covering the complete mitochondrial genome of Schizosaccharomyces pombe. By hybridizing mitochondrial gene probes from Saccharomyces cerevisiae and Neurospora crassa with restriction fragments of Schizosaccharomyces pombe mitochondrial DNA, the following homologous genes were localized on the mitochondrial genome of S. pombe: cob, cox1, cox2 and cox3, ATPase subunit 6 and 9 genes, the large rRNA gene and both types of open reading frames occurring in mitochondrial introns of various ascomycetes. The region of the genome, hybridizing with cob exon probes is separated by an intervening sequence of about 2500 bp, which is homologous with the first two introns of the cox1 gene in Saccharomyces cerevisiae (class II introns according to Michel et al. 1982). Similarly, in the cox1 homologous region, which covers about 4000 bp, two regions were detected hybridizing with class I intron probes, suggesting the existence of two cox1 introns in Schizosaccharomyces pombe. Hybridization with several specific exon probes with a determined order has revealed that cob, cox1, cox3 and the large rRNA gene are all transcribed from the same DNA strand. The low intensities of hybridization signals suggest a large evolutionary distance between Schizosaccharomyces pombe and Saccharomyces cerevisiae or Neurospora crassa mitochondrial genes. Considering the length of the mitochondrial DNA of Schizosaccharomyces pombe (about 19.4 kbp) and the expected length of the localized genes and intron sequences there is enough space left for encoding the expected set of tRNAs and the small rRNA gene. The existence of leader-, trailer-, ori- and spacer sequences or further unassigned reading frames is then restricted to a total length of about 3000 bp only.     

The RNA polymerase III-dependent family of genes in hemiascomycetes: comparative RNomics, decoding strategies, transcription and evolutionary implications

We present the first comprehensive analysis of RNA polymerase III (Pol III) transcribed genes in ten yeast genomes. This set includes all tRNA genes (tDNA) and genes coding for SNR6 (U6), SNR52, SCR1 and RPR1 RNA in the nine hemiascomycetes Saccharomyces cerevisiae, Saccharomyces castellii, Candida glabrata, Kluyveromyces waltii, Kluyveromyces lactis, Eremothecium gossypii, Debaryomyces hansenii, Candida albicans, Yarrowia lipolytica and the archiascomycete Schizosaccharomyces pombe. We systematically analysed sequence specificities of tRNA genes, polymorphism, variability of introns, gene redundancy and gene clustering. Analysis of decoding strategies showed that yeasts close to S.cerevisiae use bacterial decoding rules to read the Leu CUN and Arg CGN codons, in contrast to all other known Eukaryotes. In D.hansenii and C.albicans, we identified a novel tDNA-Leu (AAG), reading the Leu CUU/CUC/CUA codons with an unusual G at position 32. A systematic 'p-distance tree' using the 60 variable positions of the tRNA molecule revealed that most tDNAs cluster into amino acid-specific sub-trees, suggesting that, within hemiascomycetes, orthologous tDNAs are more closely related than paralogs. We finally determined the bipartite A- and B-box sequences recognized by TFIIIC. These minimal sequences are nearly conserved throughout hemiascomycetes and were satisfactorily retrieved at appropriate locations in other Pol III genes.     

Incipient mitochondrial evolution in yeasts. II. The complete sequence of the gene coding for cytochrome b in Saccharomyces douglasii reveals the presence of both new and conserved introns and discloses major differences in the fixation of mutations in evolution.

We have determined the complete sequence of the mitochondrial gene coding for cytochrome b in Saccharomyces douglasii. The gene is 6310 base-pairs long and is interrupted by four introns. The first one (1311 base-pairs) belongs to the group ID of secondary structure, contains a fragment open reading frame with a characteristic GIY ... YIG motif, is absent from Saccharomyces cerevisiae and is inserted in the same site in which introns 1 and 2 are inserted in Neurospora crassa and Podospora anserina, respectively. The next three S. douglasii introns are homologous to the first three introns of S. cerevisiae, are inserted at the same positions and display various degrees of similarity ranging from an almost complete identity (intron 2 and 4) to a moderate one (intron 3). We have compared secondary structures of intron RNAs, and nucleotide and amino acid sequences of cytochrome b exons and intron open reading frames in the two Saccharomyces species. The rules that govern fixation of mutations in exon and intron open reading frames are different: the relative proportion of mutations occurring in synonymous codons is low in some introns and high in exons. The overall frequency of mutations in cytochrome b exons is much smaller than in nuclear genes of yeasts, contrary to what has been found in vertebrates, where mitochondrial mutations are more frequent. The divergence of the cytochrome b gene is modular: various parts of the gene have changed with a different mode and tempo of evolution.     

Structural variations and optional introns in the mitochondrial DNAs of Neurospora strains isolated from nature

Mitochondrial DNAs from ten wild-type Neurospora crassa, Neurospora intermedia, and Neurospora sitophila strains collected from different geographical areas were screened for structural variations by restriction enzyme analysis. The different mtDNAs show much greater structural diversity, both within and among species, than had been apparent from previous studies of mtDNA from laboratory N. crassa strains. The mtDNAs range in size from 60 to 73 kb, and both the smallest and largest mtDNAs are found in N. crassa strains. In addition, four strains contain intramitochondrial plasmid DNAs that do not hybridize with the standard mtDNA. All of the mtDNA species have a basically similar organization. A 25-kb region that includes the rRNA genes and most tRNA genes shows very strong conservation of restriction sites in all strains. The 2.3-kb intron found in the large rRNA gene in standard N. crassa mtDNAs is present in all strains examined, including N. intermedia and N. sitophila strains. The size differences between the different mtDNAs are due to insertions or deletions that occur outside of the rRNA-tRNA region. Restriction enzyme and heteroduplex mapping suggest that four of these insertions are optional introns in the gene encoding cytochrome oxidase subunit I. Mitochondrial DNAs from different wild-type strains contain zero, one, three, or four of these introns.     

The complete mitochondrial DNA sequence of Mesostigma viride identifies this green alga as the earliest green plant divergence and predicts a highly compact mitochondrial genome in the ancestor of all green plants.

To gain insights into the nature of the mitochondrial genome in the common ancestor of all green plants, we have completely sequenced the mitochondrial DNA (mtDNA) of Mesostigma viride. This green alga belongs to a morphologically heterogeneous class (Prasinophyceae) that includes descendants of the earliest diverging green plants. Recent phylogenetic analyses of ribosomal RNAs (rRNAs) and concatenated proteins encoded by the chloroplast genome identified Mesostigma as a basal branch relative to the Streptophyta and the Chlorophyta, the two phyla that were previously thought to contain all extant green plants. The circular mitochondrial genome of Mesostigma resembles the mtDNAs of green algae occupying a basal position within the Chlorophyta in displaying a small size (42,424 bp) and a high gene density (86.6% coding sequences). It contains 65 genes that are conserved in other mtDNAs. Although none of these genes represents a novel coding sequence among green plant mtDNAs, four of them (rps1, sdh3, sdh4, and trnL[caa]) have not been reported previously in chlorophyte mtDNAs, and two others (rpl14 and trnI[gau]) have not been identified in the streptophyte mtDNAs examined so far (land-plant mtDNAs). Phylogenetic analyses of 19 concatenated mtDNA-encoded proteins favor the hypothesis that Mesostigma represents the earliest branch of green plant evolution. Four group I introns (two in rnl and two in cox1) and three group II introns (two in nad3 and one in cox2), two of which are trans-spliced at the RNA level, reside in Mesostigma mtDNA. The insertion sites of the three group II introns are unique to this mtDNA, suggesting that trans-splicing arose independently in the Mesostigma lineage and in the Streptophyta. The few structural features that can be regarded as ancestral in Mesostigma mtDNA predict that the common ancestor of all green plants had a compact mtDNA containing a minimum of 75 genes and perhaps two group I introns. Considering that the mitochondrial genome is much larger in size in land plants than in Mesostigma, we infer that mtDNA size began to increase dramatically in the Streptophyta either during the evolution of charophyte green algae or during the transition from charophytes to land plants.     

Characterization of the DNA binding features of Saccharomyces castellii Cdc13p

The principal function of Saccharomyces cerevisiae Cdc13p is to provide a loading platform to recruit complexes that provide end protection and telomere replication. We isolated the Saccharomyces castellii Cdc13p homolog (scasCdc13p) and characterized the in vitro DNA binding features of the purified recombinant scasCdc13p. The full-length scasCdc13p binds specifically to G-rich single-stranded telomeric DNA, and not to double-stranded DNA or the C-rich strand. Moreover, the minimal binding site for scasCdc13p is the octamer 5'-GTGTCTGG-3' of the S.castellii telomeric sequence. The scasCdc13p displayed a high affinity binding, where four individual nucleotide residues were found to be of most importance for the sequence specificity. Nonetheless, scasCdc13p binds the telomeric repeats from various other species, including the human. In spite of considerable divergence in telomere repeat length and sequence between these species, a conserved Cdc13p binding motif was detected. Among the budding yeasts this conserved Cdc13p binding site overlaps the Rap1p binding site. Together, these data implicate scasCdc13p as a telomere end-binding protein with a potential role in the regulation of telomere maintenance in vivo. Moreover, the results suggest that Rap1p and Cdc13p act together to preserve the conserved core present within the otherwise highly divergent btelomeric sequences among a wide variety of yeasts.     

Linear mitochondrial deoxyribonucleic acid from the yeast Hansenula mrakii.

The mitochondrial deoxyribonucleic acid (mtDNA) from a petite-negative yeast, Hansenula mrakii, was studied. A linear restriction map was constructed with 11 restriction enzymes. The linearity of the genome was confirmed by direct end labeling of the molecule, followed by restriction analysis. The molecular weight of the DNA was found to be 55,000 base pairs. This is the first linear mtDNA found in yeast species. Using specific gene probes obtained from Saccharomyces cerevisiae mtDNA, we have constructed a gene map of H. mrakii mtDNA. The arrangement of genes in this linear genome was very different from the circular mtDNA of other known yeasts.     

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.