A genome-informed higher rank classification of the biotechnologically important fungal subphylum Saccharomycotina.

The subphylum Saccharomycotina is a lineage in the fungal phylum Ascomycota that exhibits levels of genomic diversity similar to those of plants and animals. The Saccharomycotina consist of more than 1 200 known species currently divided into 16 families, one order, and one class. Species in this subphylum are ecologically and metabolically diverse and include important opportunistic human pathogens, as well as species important in biotechnological applications. Many traits of biotechnological interest are found in closely related species and often restricted to single phylogenetic clades. However, the biotechnological potential of most yeast species remains unexplored. Although the subphylum Saccharomycotina has much higher rates of genome sequence evolution than its sister subphylum, Pezizomycotina, it contains only one class compared to the 16 classes in Pezizomycotina. The third subphylum of Ascomycota, the Taphrinomycotina, consists of six classes and has approximately 10 times fewer species than the Saccharomycotina. These data indicate that the current classification of all these yeasts into a single class and a single order is an underappreciation of their diversity. Our previous genome-scale phylogenetic analyses showed that the Saccharomycotina contains 12 major and robustly supported phylogenetic clades; seven of these are current families (Lipomycetaceae, Trigonopsidaceae, Alloascoideaceae, Pichiaceae, Phaffomycetaceae, Saccharomycodaceae, and Saccharomycetaceae), one comprises two current families (Dipodascaceae and Trichomonascaceae), one represents the genus Sporopachydermia, and three represent lineages that differ in their translation of the CUG codon (CUG-Ala, CUG-Ser1, and CUG-Ser2). Using these analyses in combination with relative evolutionary divergence and genome content analyses, we propose an updated classification for the Saccharomycotina, including seven classes and 12 orders that can be diagnosed by genome content. This updated classification is consistent with the high levels of genomic diversity within this subphylum and is necessary to make the higher rank classification of the Saccharomycotina more comparable to that of other fungi, as well as to communicate efficiently on lineages that are not yet formally named. Taxonomic novelties: New classes: Alloascoideomycetes M. Groenew., Hittinger, Opulente & A. Rokas, Dipodascomycetes M. Groenew., Hittinger, Opulente & A. Rokas, Lipomycetes M. Groenew., Hittinger, Opulente, A. Rokas, Pichiomycetes M. Groenew., Hittinger, Opulente & A. Rokas, Sporopachydermiomycetes M. Groenew., Hittinger, Opulente & A. Rokas, Trigonopsidomycetes M. Groenew., Hittinger, Opulente & A. Rokas. New orders: Alloascoideomycetes: Alloascoideales M. Groenew., Hittinger, Opulente & A. Rokas; Dipodascomycetes: Dipodascales M. Groenew., Hittinger, Opulente & A. Rokas; Lipomycetes: Lipomycetales M. Groenew., Hittinger, Opulente & A. Rokas; Pichiomycetes: Alaninales M. Groenew., Hittinger, Opulente & A. Rokas, Pichiales M. Groenew., Hittinger, Opulente & A. Rokas, Serinales M. Groenew., Hittinger, Opulente & A. Rokas; Saccharomycetes: Phaffomycetales M. Groenew., Hittinger, Opulente & A. Rokas, Saccharomycodales M. Groenew., Hittinger, Opulente & A. Rokas; Sporopachydermiomycetes: Sporopachydermiales M. Groenew., Hittinger, Opulente & A. Rokas; Trigonopsidomycetes: Trigonopsidales M. Groenew., Hittinger, Opulente & A. Rokas. New families: Alaninales: Pachysolenaceae M. Groenew., Hittinger, Opulente & A. Rokas; Pichiales: Pichiaceae M. Groenew., Hittinger, Opulente & A. Rokas; Sporopachydermiales: Sporopachydermiaceae M. Groenew., Hittinger, Opulente & A. Rokas. Citation: Groenewald M, Hittinger CT, Bensch K, Opulente DA, Shen X-X, Li Y, Liu C, LaBella AL, Zhou X, Limtong S, Jindamorakot S, Gonçalves P, Robert V, Wolfe KH, Rosa CA, Boekhout T, Cadez N, Péter G, Sampaio JP, Lachance M-A, Yurkov AM, Daniel H-M, Takashima M, Boundy-Mills K, Libkind D, Aoki K, Sugita T, Rokas A (2023). A genome-informed higher rank classification of the biotechnologically important fungal subphylum Saccharomycotina. Studies in Mycology 105: 1-22. doi: 10.3114/sim.2023.105.01 This study is dedicated to the memory of Cletus P. Kurtzman (1938-2017), a pioneer of yeast taxonomy.

Yesterday's polyploids and the mystery of diploidization.

Thirty years after Susumu Ohno proposed that vertebrate genomes are degenerate polyploids, the extent to which genome duplication contributed to the evolution of the vertebrate genome, if at all, is still uncertain. Sequence-level studies on model organisms whose genomes show clearer evidence of ancient polyploidy are invaluable because they indicate what the evolutionary products of genome duplication can look like. The greatest mystery is the molecular basis of diploidization, the evolutionary process by which a polyploid genome turns into a diploid one.

Many parallel losses of infA from chloroplast DNA during angiosperm evolution with multiple independent transfers to the nucleus.

We used DNA sequencing and gel blot surveys to assess the integrity of the chloroplast gene infA, which codes for translation initiation factor 1, in >300 diverse angiosperms. Whereas most angiosperms appear to contain an intact chloroplast infA gene, the gene has repeatedly become defunct in approximately 24 separate lineages of angiosperms, including almost all rosid species. In four species in which chloroplast infA is defunct, transferred and expressed copies of the gene were found in the nucleus, complete with putative chloroplast transit peptide sequences. The transit peptide sequences of the nuclear infA genes from soybean and Arabidopsis were shown to be functional by their ability to target green fluorescent protein to chloroplasts in vivo. Phylogenetic analysis of infA sequences and assessment of transit peptide homology indicate that the four nuclear infA genes are probably derived from four independent gene transfers from chloroplast to nuclear DNA during angiosperm evolution. Considering this and the many separate losses of infA from chloroplast DNA, the gene has probably been transferred many more times, making infA by far the most mobile chloroplast gene known in plants.

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.

Relationship of codon bias to mRNA concentration and protein length in Saccharomyces cerevisiae.

In 1982, Ikemura reported a strikingly unequal usage of different synonymous codons, in five Saccharomyces cerevisiae nuclear genes having high protein levels. To study this trend in detail, we examined data from three independent studies that used oligonucleotide arrays or SAGE to estimate mRNA concentrations for nearly all genes in the genome. Correlation coefficients were calculated for the relationship of mRNA concentration to four commonly used measures of synonymous codon usage bias: the codon adaptation index (CAI), the codon bias index (CBI), the frequency of optimal codons (F(op)), and the effective number of codons (N(c)). mRNA concentration was best approximated as an exponential function of each of these four measures. Of the four, the CAI was the most strongly correlated with mRNA concentration (r(s)=0.62+/-0.01, n=2525, p<10(-17)). When we controlled for CAI, mRNA concentration and protein length were negatively correlated (partial r(s)=-0.23+/-0.01, n=4765, p<10(-17)). This may result from selection to reduce the size of abundant proteins to minimize transcriptional and translational costs. When we controlled for mRNA concentration, protein length and CAI were positively correlated (partial r(s)=0.16+/-0.01, n=4765, p<10(-17)). This may reflect more effective selection in longer genes against missense errors during translation. The correlation coefficients between the mRNA levels of individual genes, as measured by different investigators and methods, were low, in the range r(s)=0.39-0.68.

Estimation of synteny conservation and genome compaction between pufferfish (Fugu) and human.

BACKGROUND: Knowledge of the amount of gene order and synteny conservation between two species gives insights to the extent and mechanisms of divergence. The vertebrate Fugu rubripes (pufferfish) has a small genome with little repetitive sequence which makes it attractive as a model genome. Genome compaction and synteny conservation between human and Fugu were studied using data from public databases. METHODS: Intron length and map positions of human and Fugu orthologues were compared to analyse relative genome compaction and synteny conservation respectively. The divergence of these two genomes by genome rearrangement was simulated and the results were compared to the real data. RESULTS: Analysis of 199 introns in 22 orthologous genes showed an eight-fold average size reduction in Fugu, consistent with the ratio of total genome sizes. There was no consistent pattern relating the size reduction in individual introns or genes to gene base composition in either species. For genes that are neighbours in Fugu (genes from the same cosmid or GenBank entry), 40-50% have conserved synteny with a human chromosome. This figure may be underestimated by as much as two-fold, due to problems caused by incomplete human genome sequence data and the existence of dispersed gene families. Some genes that are neighbours in Fugu have human orthologues that are several megabases and tens of genes apart. This is probably caused by small inversions or other intrachromosomal rearrangements. CONCLUSIONS: Comparison of observed data to computer simulations suggests that 4000-16 000 chromosomal rearrangements have occurred since Fugu and human shared a common ancestor, implying a faster rate of rearrangement than seen in human/mouse comparisons.

Evidence for a high frequency of simultaneous double-nucleotide substitutions.

Point mutations are generally assumed to involve changes of single nucleotides. Nevertheless, the nature and known mechanisms of mutation do not exclude the possibility that several adjacent nucleotides may change simultaneously in a single mutational event. Two independent approaches are used here to estimate the frequency of simultaneous double-nucleotide substitutions. The first examines switches between TCN and AGY (where N is any nucleotide and Y is a pyrimidine) codons encoding absolutely conserved serine residues in a number of proteins from diverse organisms. The second reveals double-nucleotide substitutions in primate noncoding sequences. These two complementary approaches provide similar high estimates for the rate of doublet substitutions, on the order of 0.1 per site per billion years.

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.

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.

G+C content variation along and among Saccharomyces cerevisiae chromosomes.

Past analyses of the genome of the yeast Saccharomyces cerevisiae have revealed substantial regional variation in G+C content. Important questions remain, though, as to the origin, nature, significance, and generality of this variation. We conducted an extensive analysis of the yeast genome to try to answer these questions. Our results indicate that open reading frames (ORFs) with similar G+C contents at silent codon positions are significantly clustered on chromosomes. This clustering can be explained by very short range correlations of silent-site G+C contents at neighboring ORFs. ORFs of high silent-site G+C content are disproportionately concentrated on shorter chromosomes, which causes a negative relationship between chromosome length and G+C content. Contrary to previous reports, there is no correlation between gene density and silent-site G+C content in yeast. Chromosome III is atypical in many regards, and possible reasons for this are discussed.

Proteinases and associated genes of parasitic helminths.

Many parasites have deployed proteinases to accomplish some of the tasks imposed by a parasitic life style, including tissue penetration, digestion of host tissue for nutrition and evasion of host immune responses. Information on proteinases from trematodes, cestodes and nematode parasites is reviewed, concentrating on those worms of major medical and economical importance. Their biochemical characterization is discussed, along with their putative biological roles and, where available, their associated genes. For example, proteinases expressed by the various stages of the schistosome life-cycle, in particular the well-characterized cercarial elastase which is involved in the penetration of the host skin and the variety of proteinases, such as cathepsin B (Sm31), cathepsin L1, cathepsin L2, cathepsin D, cathepsin C and legumain (Sm32), which are believed to be involved in the catabolism of host haemoglobin. The various endo- and exoproteinases of Fasciola hepatica, the causative agent of liver fluke disease, are reviewed, and recent reports of how these enzymes have been successfully employed in cocktail vaccines are discussed. The various proteinases of cestodes and of the diverse superfamilies of parasitic nematodes are detailed, with special attention being given to those parasites for which most is known, including species of Taenia, Echinococcus, Spirometra, Necator, Acylostoma and Haemonchus. By far the largest number of papers in the literature and entries to the sequence data bases dealing with proteinases of parasitic helminths report on enzymes belonging to the papain superfamily of cysteine proteinases. Accordingly, the final section of the review is devoted to a phylogenetic analysis of this superfamily using over 150 published sequences. This analysis shows that the papain superfamily can be divided into two major branches. Branch A contains the cathepin Bs, the cathepsin Cs and a novel family termed cathepsin Xs, while Branch B contains the cruzipains, cathepsin Ls, papain-like and aleurain/cathepsin H-like proteinases. The relationships of the helminth proteinases, and similar proteinases from protozoan parasites and other organisms, within these groups are discussed.

Gene duplication and gene conversion in the Caenorhabditis elegans genome.

A comprehensive analysis of duplication and gene conversion for 7394 Caenorhabditis elegans genes (about half the expected total for the genome) is presented. Of the genes examined, 40% are involved in duplicated gene pairs. Intrachromosomal or cis gene duplications occur approximately two times more often than expected. In general the closer the members of duplicated gene pairs are, the more likely it is that gene orientation is conserved. Gene conversion events are detectable between only 2% of the duplicated pairs. Even given the excesses of cis duplications, there is an excess of gene conversion events between cis duplicated pairs on every chromosome except the X chromosome. The relative rates of cis and trans gene conversion and the negative correlation between conversion frequency and DNA sequence divergence for unconverted regions of converted pairs are consistent with previous experimental studies in yeast. Three recent, regional duplications, each spanning three genes are described. All three have already undergone substantial deletions spanning hundreds of base pairs. The relative rates of duplication and deletion may contribute to the compactness of the C. elegans genome.

Proteome composition and codon usage in spirochaetes: species-specific and DNA strand-specific mutational biases.

The genomes of the spirochaetes Borrelia burgdorferi and Treponema pallidum show strong strand-specific skews in nucleotide composition, with the leading strand in replication being richer in G and T than the lagging strand in both species. This mutation bias results in codon usage and amino acid composition patterns that are significantly different between genes encoded on the two strands, in both species. There are also substantial differences between the species, with T.pallidum having a much higher G+C content than B. burgdorferi. These changes in amino acid and codon compositions represent neutral sequence change that has been caused by strong strand- and species-specific mutation pressures. Genes that have been relocated between the leading and lagging strands since B. burgdorferi and T.pallidum diverged from a common ancestor now show codon and amino acid compositions typical of their current locations. There is no evidence that translational selection operates on codon usage in highly expressed genes in these species, and the primary influence on codon usage is whether a gene is transcribed in the same direction as replication, or opposite to it. The dnaA gene in both species has codon usage patterns distinctive of a lagging strand gene, indicating that the origin of replication lies downstream of this gene, possibly within dnaN. Our findings strongly suggest that gene-finding algorithms that ignore variability within the genome may be flawed.

Base composition skews, replication orientation, and gene orientation in 12 prokaryote genomes.

Variation in GC content, GC skew and AT skew along genomic regions was examined at third codon positions in completely sequenced prokaryotes. Eight out of nine eubacteria studied show GC and AT skews that change sign at the origin of replication. The leading strand in DNA replication is G-T rich at codon position 3 in six eubacteria, but C-T rich in two Mycoplasma species. In M. genitalium the AT and GC skews are symmetrical around the origin and terminus of replication, whereas its GC content variation has been shown to have a centre of symmetry elsewhere in the genome. Borrelia burgdorferi and Treponema pallidum show extraordinary extents of base composition skew correlated with direction of DNA replication. Base composition skews measured at third codon positions probably reflect mutational biases, whereas those measured over all bases in a sequence (or at codon positions 1 and 2) can be strongly affected by protein considerations due to the tendency in some bacteria for genes to be transcribed in the same direction that they are replicated. Consequently in some species the direction of skew for total genomic DNA is opposite to that for codon position 3.

A subset of conserved tRNA genes in plastid DNA of nongreen plants.

The plastid genome of the nonphotosynthetic parasitic plant Epifagus virginiana contains only 17 of the 30 tRNA genes normally found in angiosperm plastid DNA. Although this is insufficient for translation, the genome is functional, so import of cytosolic tRNAs into plastids has been suggested. This raises the question of whether the tRNA genes that remain in E. virginiana plastid DNA are active or have just fortuitously escaped deletion. We report the sequences of 20 plastid tRNA loci from Orobanche minor, which shares a nonphotosynthetic ancestor with E. virginiana. The two species have 9 intact tRNA genes in common, the others being defunct in one or both species. The intron-containing trnLUAA gene is absent from E. virginiana, but it is intact, transcribed, and spliced in O. minor. The shared intact genes are better conserved than intergenic sequences, which indicates that these genes are being maintained by natural selection and, therefore, must be functional. For the most part, the tRNA species conserved in nonphotosynthetic plastids are also those that have never been found to be imported in plant mitochondria, which suggests that the same rules may govern tRNA import in the two organelles. A small photosynthesis gene, psbI, is still intact in O. minor, and computer simulations show that some small nonessential genes have an appreciable chance of escaping deletion.

Extent of genomic rearrangement after genome duplication in yeast.

Whole-genome duplication approximately 10(8) years ago was proposed as an explanation for the many duplicated chromosomal regions in Saccharomyces cerevisiae. Here we have used computer simulations and analytic methods to estimate some parameters describing the evolution of the yeast genome after this duplication event. Computer simulation of a model in which 8% of the original genes were retained in duplicate after genome duplication, and 70-100 reciprocal translocations occurred between chromosomes, produced arrangements of duplicated chromosomal regions very similar to the map of real duplications in yeast. An analytical method produced an independent estimate of 84 map disruptions. These results imply that many smaller duplicated chromosomal regions exist in the yeast genome in addition to the 55 originally reported. We also examined the possibility of determining the original order of chromosomal blocks in the ancestral unduplicated genome, but this cannot be done without information from one or more additional species. If the genome sequence of one other species (such as Kluyveromyces lactis) were known it should be possible to identify 150-200 paired regions covering the whole yeast genome and to reconstruct approximately two-thirds of the original order of blocks of genes in yeast. Rates of interchromosome translocation in yeast and mammals appear similar despite their very different rates of homologous recombination per kilobase.

Evolution of gene order and chromosome number in Saccharomyces, Kluyveromyces and related fungi.

The extent to which the order of genes along chromosomes is conserved between Saccharomyces cerevisiae and related species was studied by analysing data from DNA sequence database. As expected, the extent of gene order conservation decreases with increasing evolutionary distance. About 59% of adjacent gene pairs in Kluyveromyces lactis or K. marxianus are also adjacent in S. cerevisiae, and a further 16% of Kluyveromyces neighbours can be explained in terms of the inferred ancestral gene order in Saccharomyces prior to the occurrence of an ancient whole-genome duplication. Only 13% of Candida albicans linkages, and no Schizosaccharomyces pombe linkages, are conserved. Analysis of gene order arrangements, chromosome numbers, and ribosomal RNA sequences suggests that genome duplication occurred before the divergence of the four species in Saccharomyces sensu stricto (all of which have 16 chromosomes), but after this lineage had diverged from Saccharomyces kluyveri and the Kluyveromyces lactislmarxianus species assemblage.

Eukaryote genome duplication - where's the evidence?

Several eukaryotes, including maize, yeast and Xenopus, are degenerate polyploids formed by relatively recent whole-genome duplications. Ohno's conjecture that more ancient genome duplications occurred in an ancestor of vertebrates is probably at least partly true but the present shortage of gene sequence and map information from vertebrates makes it difficult to either prove or disprove this hypothesis. Candidate paralogous segments in mammalian genomes have been identified but the lack of statistical rigour means that many of the proposals in the literature are probably artefacts.

Molecular evidence for an ancient duplication of the entire yeast genome.

Gene duplication is an important source of evolutionary novelty. Most duplications are of just a single gene, but Ohno proposed that whole-genome duplication (polyploidy) is an important evolutionary mechanism. Many duplicate genes have been found in Saccharomyces cerevisiae, and these often seem to be phenotypically redundant. Here we show that the arrangement of duplicated genes in the S. cerevisiae genome is consistent with Ohno's hypothesis. We propose a model in which this species is a degenerate tetraploid resulting from a whole-genome duplication that occurred after the divergence of Saccharomyces from Kluyveromyces. Only a small fraction of the genes were subsequently retained in duplicate (most were deleted), and gene order was rearranged by many reciprocal translocations between chromosomes. Protein pairs derived from this duplication event make up 13% of all yeast proteins, and include pairs of transcription factors, protein kinases, myosins, cyclins and pheromones. Tetraploidy may have facilitated the evolution of anaerobic fermentation in Saccharomyces.

Accelerated evolution of sites undergoing mRNA editing in plant mitochondria and chloroplasts.

The selective constraints influencing mRNA editing in plant organelles are largely unknown. To investigate these, we compared patterns of editing between monocot and dicot mitochondrial mRNA. On average, 24% of sites that are edited form C to U in one species have been substituted during evolution by a genomic T in the other: this is four times the rate of evolution seen at nonedited synonymously variable C residues. A similar, but weaker trend (not statistically significant) is seen at sites edited in chloroplast mRNA. The elevated substitution rate does not appear to be a consequence of a higher mutability of the trinucleotide motif (T-C-purine) associated with editing. nor to be a result of reverse transcription from mature mRNA. Selection to replace the genomic C with a T may account for the accelerated evolution, either due to elimination of inefficient transcripts and protein products or as a consequence of the prior loss of components outside the edit site which are necessary for editing; the latter hypothesis is supported by the frequent loss of editing without genomic mutation at third codon positions. Whatever the cause, the rapid rate of evolution indicates that editing confers little selective advantage at most sites.