Abstract

Budding yeasts inhabit a range of environments by exploiting various metabolic traits. The genetic bases for these traits are mostly unknown, preventing their addition or removal in a chassis organism for metabolic engineering. Insight into the evolution of orthologs, paralogs and xenologs in the yeast pan-genome can help bridge these genotypes; however, existing phylogenomic databases do not span diverse yeasts, and sometimes cannot distinguish between these homologs. To help understand the molecular evolution of these traits in yeasts, we created Analyzing Yeasts by Reconstructing Ancestry of Homologs (AYbRAH), an open-source database of predicted and manually curated ortholog groups for 33 diverse fungi and yeasts in Dikarya, spanning 600 million years of evolution. OrthoMCL and OrthoDB were used to cluster protein sequence into ortholog and homolog groups, respectively; MAFFT and PhyML reconstructed the phylogeny of all homolog groups. Ortholog assignments for enzymes and small metabolite transporters were compared to their phylogenetic reconstruction, and curated to resolve any discrepancies. Information on homolog and ortholog groups can be viewed in the AYbRAH web portal (https://lmse.github.io/aybrah/), including functional annotations, predictions for mitochondrial localization and transmembrane domains, literature references and phylogenetic reconstructions. Ortholog assignments in AYbRAH were compared to HOGENOM, KEGG Orthology, OMA, eggNOG and PANTHER. PANTHER and OMA had the most congruent ortholog groups with AYbRAH, while the other phylogenomic databases had greater amounts of under-clustering, over-clustering or no ortholog annotations for proteins. Future plans are discussed for AYbRAH, and recommendations are made for other research communities seeking to create curated ortholog databases.

Introduction

Yeasts are unicellular fungi that exploit diverse habitats on every continent, including the gut of wood boring beetles, insect frass, tree exudate, rotting wood, rotting cactus tissue, soil, brine solutions and fermenting juice (1). The most widely studied yeasts are true budding yeasts, which span roughly 400 million years of evolution in the subphylum Saccharomycotina (2), and possess a broad range of traits important to metabolic engineering. These include citrate and lipid accumulation in Yarrowia (3) and Lipomyces (4), thermotolerance in multiple lineages (5, 6), acid tolerance in Pichia (7) and Zygosaccharomyces (8), methanol utilization in Komagataella (9), osmotolerance in Debaryomyces (10), xylose to ethanol fermentation in multiple yeast lineages (11–13), alternative nuclear codon assignments (14), glucose and acetic acid co-consumption in Zygosaccharomyces (15) and aerobic ethanol production (the Crabtree effect) in multiple lineages (16–19). The complete genetic bases of these traits are mostly unknown, preventing their addition or removal in a chassis organism for biotechnology.

The distinction between orthologs, paralogs, ohnologs and xenologs plays an important role in bridging the genotype–phenotype gap across the tree of life (20). Briefly, orthologs are genes that arise from speciation and typically have a conserved function; paralogs and ohnologs emerge from locus and whole genome duplications, respectively, and may have a novel function; xenologs derive from horizontal gene transfer between organisms and do not necessarily have conserved function (21, 22). Knowledge of these types of genes has played an important role in deciphering Saccharomyces cerevisiae’s physiology. For example, the Adh2p paralog in S. cerevisiae consumes ethanol and evolved from an ancient Adh1p duplication whose kinetics favored ethanol production (23); the Saccharomycetaceae Whole Genome Duplication led to the MPC2 and MPC3 ohnologs in the Saccharomyces genus, which encode the fermentative and respirative subunits of the mitochondrial pyruvate carrier (24), respectively; the URA1 xenolog from Lactobacillales enables uracil to be synthesized anaerobically in most Saccharomycetaceae yeasts (25). These examples demonstrate how understanding the origin of genes has narrowed the genotype–phenotype gap for fermentation in Saccharomycetaceae.

Many genomics studies have focused on the Saccharomycetaceae family, and to a lesser extent the CTG clade (26), but more can be learned about yeast metabolism by studying its evolution over a longer time horizon, especially with yeasts having deeper phylogeny (27). If we could study the metabolism of the mother of all budding yeasts, which we refer to as the Proto-Yeast, we could track the gains and losses of orthologs and function in all of her descendants to bridge various genotype–phenotype gaps. Proto-Yeast has evolved from her original state, making this direct study impossible, but we can reconstruct her metabolism through her living descendants. In recent yeasts, dozens of yeasts with deep phylogeny have been sequenced (28), paving the way for greater insight into the evolution of metabolism in yeasts beyond Saccharomycetaceae.

Ortholog databases are critical to facilitating comparative genomics studies and inferring protein function. Most of these databases are constructed using graph-based methods that rely on sequence similarity, while fewer databases use tree-based methods (29). Existing ortholog databases do not span diverse yeasts (Figure 1), and sometimes cannot distinguish between orthologs and paralogs (Tables S1 and S2). In addition to these databases, orthologs are identified on an ad hoc basis with OrthoMCL for comparative genomics studies (30, 31), or with the reciprocal best hit (RBH) method for genome-scale network reconstructions (GENREs) (32); these ortholog assignments often lack transparency or traceability, and therefore cannot be scrutinized or continuously improved by research communities. To solve these outlined problems, and ultimately improve our understanding of budding yeast physiology, we present Analyzing Yeasts by Reconstructing Ancestry of Homologs (AYbRAH; Figure 2). AYbRAH, derived from the Hebrew name Abra, mother of many, is an open-source database of predicted and manually curated orthologs, their function and their origin. The initial AYbRAH database was constructed using OrthoMCL and OrthoDB. PhyML was used to reconstruct the phylogeny of each homolog group. AYbRAH ortholog assignments for enzymes and small metabolite transporters were compared against their phylogenetic reconstruction and curated to resolve any discrepancies. We discuss the information available in the AYbRAH web portal (https://lmse.github.io/aybrah/), issues that arose from reviewing the accuracy of ortholog predictions, compare AYbRAH to established phylogenomic databases, discuss the benefits of open-source ortholog databases, future directions for AYbRAH, and offer recommendations to research communities looking to develop ortholog databases for other taxa.

Ortholog database coverage for fungal and yeast genomes in AYbRAH, YGOB, CGOB, PANTHER, HOGENOM, KO, OMA and eggNOG. Ortholog assignments based on the manual curation of sequence similarity and synteny are shown in green columns, tree-based methods in red columns, graph-based methods in blue columns and a hybrid graph and tree-based method in the purple column. Many ortholog databases are well represented in Saccharomycetaceae and the CTG clade, which had their genomes sequenced during the 2000s (26). AYbRAH has ortholog assignments for species in Pichiaceae, Phaffomycetaceae and several incertae sedis families, which are not well represented in other ortholog databases, as these yeasts were recently sequenced (28). The well established phylogenomic databases span other yeast species not shown in this phylogeny, but they mostly belong to Saccharomycetaceae or the CTG clade.
Figure 1

Ortholog database coverage for fungal and yeast genomes in AYbRAH, YGOB, CGOB, PANTHER, HOGENOM, KO, OMA and eggNOG. Ortholog assignments based on the manual curation of sequence similarity and synteny are shown in green columns, tree-based methods in red columns, graph-based methods in blue columns and a hybrid graph and tree-based method in the purple column. Many ortholog databases are well represented in Saccharomycetaceae and the CTG clade, which had their genomes sequenced during the 2000s (26). AYbRAH has ortholog assignments for species in Pichiaceae, Phaffomycetaceae and several incertae sedis families, which are not well represented in other ortholog databases, as these yeasts were recently sequenced (28). The well established phylogenomic databases span other yeast species not shown in this phylogeny, but they mostly belong to Saccharomycetaceae or the CTG clade.

Table 1

Fungal and yeast strain genomes in AYbRAH. Protein sequences were downloaded from UniProt or MycoCosm. Species were assigned to monophyletic or paraphyletic groups based on divergence time with S. cerevisiae

SpeciesStrainGroupDatabaseReference
Rhodotorula graminisWP1Saccharomycotina outgroupMycoCosm(73)
Saitoella complicataNRRL Y-17804MycoCosm(28)
Schizosaccharomyces pombe972h-UniProt(74)
Aspergillus nigerCBS 513.88UniProt(75)
Neurospora crassaCBS708.71UniProt(76)
Trichoderma reeseiQM6aUniProt(77)
Lipomyces starkeyiNRRL Y-11557basal SaccharomycotinaMycoCosm(28)
Yarrowia lipolyticaCLIB 122UniProt(78)
Blastobotrys adeninivoransLS3MycoCosm(79)
Nadsonia fulvescens var. elongataDSM 6959MycoCosm(28)
Ascoidea rubescensNRRL Y17699MycoCosm(28)
Pachysolen tannophilusNRRL Y-2460PichiaceaeMycoCosm(28)
Komagataella phaffiiGS115UniProt(80)
Kuraishia capsulataCBS 1993UniProt(81)
Ogataea arabinofermentansNRRL YB-2248MycoCosm(28)
Ogataea parapolymorphaNRRL Y-7560UniProt(83)
Dekkera bruxellensisCBS 2499MycoCosm(82)
Pichia membranifaciensNRRL Y-2026MycoCosm(28)
Pichia kudriavzeviiSD108UniProt(84)
Babjeviella inositovoraNRRL Y-12698CTG cladeMycoCosm(28)
Metschnikowia bicuspidataNRRL YB-4993MycoCosm(28)
Meyerozyma guilliermondiiCBS 566UniProt(85)
Debaryomyces hanseniiCBS 767UniProt(78)
Scheffersomyces stipitisCBS 6054UniProt(86)
Spathaspora passalidarumNRRL Y-27907UniProt(30)
Wickerhamomyces anomalusNRRL Y-366-8Phaffomycetaceae & SaccharomycodaceaeMycoCosm(28)
Cyberlindnera jadiniiNRRL Y-1542MycoCosm(28)
Hanseniaspora valbyensisNRRL Y-1626MycoCosm(28)
Kluyveromyces lactisCBS 2359SaccharomycetaceaeUniProt(78)
Lachancea thermotoleransCBS 6340UniProt(87)
Zygosaccharomyces rouxiiCBS 732UniProt(87)
Saccharomyces cerevisiaeS288CUniProt(88)
Vanderwaltozyma polysporaDSM 70294UniProt(89)
SpeciesStrainGroupDatabaseReference
Rhodotorula graminisWP1Saccharomycotina outgroupMycoCosm(73)
Saitoella complicataNRRL Y-17804MycoCosm(28)
Schizosaccharomyces pombe972h-UniProt(74)
Aspergillus nigerCBS 513.88UniProt(75)
Neurospora crassaCBS708.71UniProt(76)
Trichoderma reeseiQM6aUniProt(77)
Lipomyces starkeyiNRRL Y-11557basal SaccharomycotinaMycoCosm(28)
Yarrowia lipolyticaCLIB 122UniProt(78)
Blastobotrys adeninivoransLS3MycoCosm(79)
Nadsonia fulvescens var. elongataDSM 6959MycoCosm(28)
Ascoidea rubescensNRRL Y17699MycoCosm(28)
Pachysolen tannophilusNRRL Y-2460PichiaceaeMycoCosm(28)
Komagataella phaffiiGS115UniProt(80)
Kuraishia capsulataCBS 1993UniProt(81)
Ogataea arabinofermentansNRRL YB-2248MycoCosm(28)
Ogataea parapolymorphaNRRL Y-7560UniProt(83)
Dekkera bruxellensisCBS 2499MycoCosm(82)
Pichia membranifaciensNRRL Y-2026MycoCosm(28)
Pichia kudriavzeviiSD108UniProt(84)
Babjeviella inositovoraNRRL Y-12698CTG cladeMycoCosm(28)
Metschnikowia bicuspidataNRRL YB-4993MycoCosm(28)
Meyerozyma guilliermondiiCBS 566UniProt(85)
Debaryomyces hanseniiCBS 767UniProt(78)
Scheffersomyces stipitisCBS 6054UniProt(86)
Spathaspora passalidarumNRRL Y-27907UniProt(30)
Wickerhamomyces anomalusNRRL Y-366-8Phaffomycetaceae & SaccharomycodaceaeMycoCosm(28)
Cyberlindnera jadiniiNRRL Y-1542MycoCosm(28)
Hanseniaspora valbyensisNRRL Y-1626MycoCosm(28)
Kluyveromyces lactisCBS 2359SaccharomycetaceaeUniProt(78)
Lachancea thermotoleransCBS 6340UniProt(87)
Zygosaccharomyces rouxiiCBS 732UniProt(87)
Saccharomyces cerevisiaeS288CUniProt(88)
Vanderwaltozyma polysporaDSM 70294UniProt(89)
Table 1

Fungal and yeast strain genomes in AYbRAH. Protein sequences were downloaded from UniProt or MycoCosm. Species were assigned to monophyletic or paraphyletic groups based on divergence time with S. cerevisiae

SpeciesStrainGroupDatabaseReference
Rhodotorula graminisWP1Saccharomycotina outgroupMycoCosm(73)
Saitoella complicataNRRL Y-17804MycoCosm(28)
Schizosaccharomyces pombe972h-UniProt(74)
Aspergillus nigerCBS 513.88UniProt(75)
Neurospora crassaCBS708.71UniProt(76)
Trichoderma reeseiQM6aUniProt(77)
Lipomyces starkeyiNRRL Y-11557basal SaccharomycotinaMycoCosm(28)
Yarrowia lipolyticaCLIB 122UniProt(78)
Blastobotrys adeninivoransLS3MycoCosm(79)
Nadsonia fulvescens var. elongataDSM 6959MycoCosm(28)
Ascoidea rubescensNRRL Y17699MycoCosm(28)
Pachysolen tannophilusNRRL Y-2460PichiaceaeMycoCosm(28)
Komagataella phaffiiGS115UniProt(80)
Kuraishia capsulataCBS 1993UniProt(81)
Ogataea arabinofermentansNRRL YB-2248MycoCosm(28)
Ogataea parapolymorphaNRRL Y-7560UniProt(83)
Dekkera bruxellensisCBS 2499MycoCosm(82)
Pichia membranifaciensNRRL Y-2026MycoCosm(28)
Pichia kudriavzeviiSD108UniProt(84)
Babjeviella inositovoraNRRL Y-12698CTG cladeMycoCosm(28)
Metschnikowia bicuspidataNRRL YB-4993MycoCosm(28)
Meyerozyma guilliermondiiCBS 566UniProt(85)
Debaryomyces hanseniiCBS 767UniProt(78)
Scheffersomyces stipitisCBS 6054UniProt(86)
Spathaspora passalidarumNRRL Y-27907UniProt(30)
Wickerhamomyces anomalusNRRL Y-366-8Phaffomycetaceae & SaccharomycodaceaeMycoCosm(28)
Cyberlindnera jadiniiNRRL Y-1542MycoCosm(28)
Hanseniaspora valbyensisNRRL Y-1626MycoCosm(28)
Kluyveromyces lactisCBS 2359SaccharomycetaceaeUniProt(78)
Lachancea thermotoleransCBS 6340UniProt(87)
Zygosaccharomyces rouxiiCBS 732UniProt(87)
Saccharomyces cerevisiaeS288CUniProt(88)
Vanderwaltozyma polysporaDSM 70294UniProt(89)
SpeciesStrainGroupDatabaseReference
Rhodotorula graminisWP1Saccharomycotina outgroupMycoCosm(73)
Saitoella complicataNRRL Y-17804MycoCosm(28)
Schizosaccharomyces pombe972h-UniProt(74)
Aspergillus nigerCBS 513.88UniProt(75)
Neurospora crassaCBS708.71UniProt(76)
Trichoderma reeseiQM6aUniProt(77)
Lipomyces starkeyiNRRL Y-11557basal SaccharomycotinaMycoCosm(28)
Yarrowia lipolyticaCLIB 122UniProt(78)
Blastobotrys adeninivoransLS3MycoCosm(79)
Nadsonia fulvescens var. elongataDSM 6959MycoCosm(28)
Ascoidea rubescensNRRL Y17699MycoCosm(28)
Pachysolen tannophilusNRRL Y-2460PichiaceaeMycoCosm(28)
Komagataella phaffiiGS115UniProt(80)
Kuraishia capsulataCBS 1993UniProt(81)
Ogataea arabinofermentansNRRL YB-2248MycoCosm(28)
Ogataea parapolymorphaNRRL Y-7560UniProt(83)
Dekkera bruxellensisCBS 2499MycoCosm(82)
Pichia membranifaciensNRRL Y-2026MycoCosm(28)
Pichia kudriavzeviiSD108UniProt(84)
Babjeviella inositovoraNRRL Y-12698CTG cladeMycoCosm(28)
Metschnikowia bicuspidataNRRL YB-4993MycoCosm(28)
Meyerozyma guilliermondiiCBS 566UniProt(85)
Debaryomyces hanseniiCBS 767UniProt(78)
Scheffersomyces stipitisCBS 6054UniProt(86)
Spathaspora passalidarumNRRL Y-27907UniProt(30)
Wickerhamomyces anomalusNRRL Y-366-8Phaffomycetaceae & SaccharomycodaceaeMycoCosm(28)
Cyberlindnera jadiniiNRRL Y-1542MycoCosm(28)
Hanseniaspora valbyensisNRRL Y-1626MycoCosm(28)
Kluyveromyces lactisCBS 2359SaccharomycetaceaeUniProt(78)
Lachancea thermotoleransCBS 6340UniProt(87)
Zygosaccharomyces rouxiiCBS 732UniProt(87)
Saccharomyces cerevisiaeS288CUniProt(88)
Vanderwaltozyma polysporaDSM 70294UniProt(89)

AYbRAH workflow for ortholog curation. A total of 33 fungal and yeast proteomes were downloaded from UniProt and MycoCosm. BLASTP computed the sequence similarity between all proteins. OrthoMCL clustered the proteins into putative Fungal Ortholog Groups (FOGs) using the BLASTP results. FOGs were clustered into HOmolog Groups (HOGs) using Fungi-level homolog assignments from OrthoDB. Multiple sequence alignments for each HOG were obtained with MAFFT, and 100 bootstrap phylogenetic trees were reconstructed with PhyML. The consensus phylogenetic trees for enzymes and transporters were reviewed and curated to differentiate between orthologs, paralogs, ohnologs and xenologs.
Figure 2

AYbRAH workflow for ortholog curation. A total of 33 fungal and yeast proteomes were downloaded from UniProt and MycoCosm. BLASTP computed the sequence similarity between all proteins. OrthoMCL clustered the proteins into putative Fungal Ortholog Groups (FOGs) using the BLASTP results. FOGs were clustered into HOmolog Groups (HOGs) using Fungi-level homolog assignments from OrthoDB. Multiple sequence alignments for each HOG were obtained with MAFFT, and 100 bootstrap phylogenetic trees were reconstructed with PhyML. The consensus phylogenetic trees for enzymes and transporters were reviewed and curated to differentiate between orthologs, paralogs, ohnologs and xenologs.

Annotation features of a sample phylogenetic tree in AYbRAH. Square and circle leaves indicate protein sequences in Basidiomycota or Ascomycota, respectively. Leaf nodes are colored based on taxonomic groups. Circle leaves are used for proteins with no paralogs in the same species, whereas sphere leaves are used to designate proteins with paralogs in the same species. Vertical bold lines indicate species-lineage expansions, which are sometimes called in-paralogs or co-orthologs (61). Horizontal bold lines designate S. cerevisiae proteins, which is the most widely studied eukaryote. Dashed lines indicate the most anciently diverged protein sequence in the ortholog group. Ortholog groups can be identified by color groups to help the visual inspection of ortholog assignments. The leaf names include a three-letter species code and a sequence accession. Internal nodes are labeled with the bootstrap values from phylogenetic reconstruction with PhyML.
Figure 3

Annotation features of a sample phylogenetic tree in AYbRAH. Square and circle leaves indicate protein sequences in Basidiomycota or Ascomycota, respectively. Leaf nodes are colored based on taxonomic groups. Circle leaves are used for proteins with no paralogs in the same species, whereas sphere leaves are used to designate proteins with paralogs in the same species. Vertical bold lines indicate species-lineage expansions, which are sometimes called in-paralogs or co-orthologs (61). Horizontal bold lines designate S. cerevisiae proteins, which is the most widely studied eukaryote. Dashed lines indicate the most anciently diverged protein sequence in the ortholog group. Ortholog groups can be identified by color groups to help the visual inspection of ortholog assignments. The leaf names include a three-letter species code and a sequence accession. Internal nodes are labeled with the bootstrap values from phylogenetic reconstruction with PhyML.

Localization predictions for internal NADH dehydrogenase (NDI1_YEAST) in AYbRAH. (A) Histogram plots are shown for mitochondrial localization predictions of Ndi1p orthologs Ndi1p predicted by Predotar, TargetP and MitoProt. (B) Transmembrane domain predictions computed for orthologous proteins by the Phobius web server.
Figure 4

Localization predictions for internal NADH dehydrogenase (NDI1_YEAST) in AYbRAH. (A) Histogram plots are shown for mitochondrial localization predictions of Ndi1p orthologs Ndi1p predicted by Predotar, TargetP and MitoProt. (B) Transmembrane domain predictions computed for orthologous proteins by the Phobius web server.

Methods

Initial construction of AYbRAH

AYbRAH was created by combining several algorithms and databases in a pipeline (Figure 2). A total of 212 836 protein sequences from 33 organisms (Table 1) in Dikarya were downloaded from UniProt (33) and MycoCosm (34). OrthoMCL (35) clustered protein sequences into putative Fungal Ortholog Groups (FOGs); default parameters were used for BLASTP and OrthoMCL. The FOGs from OrthoMCL were coalesced into HOmolog Groups (HOGs) using Fungi-level homolog group assignments from OrthoDB v8 (36).

AYbRAH curation

Multiple sequence alignments were obtained for each HOG with MAFFT v7.245 (37) using a gap and extension penalty of 1.5. A total of 100 bootstrap trees were reconstructed for each HOG with PhyML v3.2.0 (38), optimized for tree topology and branch length. Consensus phylogenetic trees were generated for each HOG with SumTrees from DendroPy v4.1.0 (39), and trees were rendered with ETE v3 (40). The phylogenetic reconstruction for enzymes and metabolite transporters were reviewed when OrthoMCL failed to differentiate between orthologs and paralogs, caused by over-clustering (Figure 5), or when orthologous proteins were dispersed into multiple ortholog groups, caused by under-clustering (Figure 6). Orthologs were identified by visual inspection of the phylogenetic trees or with a custom ETE 3-based script (40).

Example of over-clustering by OrthoMCL with the hexokinase family and its curation in AYbRAH. A gene duplication of HXK2 in Pichiaceae led to the HXK3 paralog. HXK2 was subsequently lost in Ogataea parapolymorpha but maintained in Komagataella phaffii. OrthoMCL was unable to differentiate between the Hxk2p and Hxk3p orthologs. Both ortholog groups are also assigned to the same Fungi-level ortholog group in OrthoDB.
Figure 5

Example of over-clustering by OrthoMCL with the hexokinase family and its curation in AYbRAH. A gene duplication of HXK2 in Pichiaceae led to the HXK3 paralog. HXK2 was subsequently lost in Ogataea parapolymorpha but maintained in Komagataella phaffii. OrthoMCL was unable to differentiate between the Hxk2p and Hxk3p orthologs. Both ortholog groups are also assigned to the same Fungi-level ortholog group in OrthoDB.

Example of under-clustering by OrthoMCL in the FLO8 ortholog group and its curation in AYbRAH. OrthoMCL dispersed the Flo8p proteins into multiple ortholog groups due to the low sequence similarity between the proteins. The proteins were merged into one ortholog group.
Figure 6

Example of under-clustering by OrthoMCL in the FLO8 ortholog group and its curation in AYbRAH. OrthoMCL dispersed the Flo8p proteins into multiple ortholog groups due to the low sequence similarity between the proteins. The proteins were merged into one ortholog group.

Annotating additional proteins

Additional steps were required to assign proteins to ortholog groups because OrthoMCL did not cluster all related proteins to ortholog groups, or because whole genome protein annotations were incomplete. First, proteins in OrthoDB homolog groups were added to new FOGs if they were not assigned to any FOG by OrthoMCL. Next, each organism had its genome nucleotide sequence queried by a protein sequence of the species closest relative for each FOG using TBLASTN (expect threshold of 1e-20). Annotated proteins were then queried against the TBLASTN hits to determine which proteins were annotated but not assigned to a FOG by OrthoMCL (misidentified) and which proteins were unannotated despite a match in its nucleotide sequence (unidentified). Proteins identified via TBLASTN with a sequence length <75% of the mean FOG sequence length were discarded from the candidate list. The remaining proteins were assigned to a HOG by its best hit via BLASTP, and to a FOG with pplacer (41) via the MAFFT add alignment option. The following examples highlight how misidentified and unidentified protein annotations were resolved in AYbRAH, respectively. First, Cybja1_169606 (A0A1E4RV95), which encodes NADP-dependent isocitrate dehydrogenase in Cyberlindnera jadinii, was not assigned to any ortholog group by OrthoMCL despite its high sequence similarity to other proteins. It was added to FOG00618 by pplacer (41) with a likelihood weight ratio of 1. Second, no 60S ribosomal protein L6 (FOG00006) was present in Meyerozyma guilliermondii’s protein annotation; it was identified by TBLASTN, annotated as mgu_AYbRAH_00173, and added to FOG00006 by pplacer with a 0.79 likelihood weight ratio (41).

Comparison of ortholog groups

AYbRAH ortholog assignments were compared to OMA (42), PANTHER (43), HOGENOM (44), eggNOG (45) and KEGG Orthology (46). Phylogenomic annotations were downloaded from UniProt. Ortholog groups were assessed as congruent, over-clustered, under-clustered, over and under-clustered or no ortholog assignment relative to AYbRAH. AYbRAH ortholog groups were only compared with a database if an ortholog group in AYbRAH had proteins from species present in the other ortholog database. For example, FOG19691 consists of proteins from Ascoidea rubescens, Pachysolen tannophilus, Kuraishia capsulata, Ogataea parapolymorpha, Dekkera bruxellensis, Pichia kudriavzevii, Pichia membranifaciens, Babjeviella inositovora, Wickerhamomyces anomalus and C. jadinii. None of the phylogenomic databases have ortholog assignments for these organisms, and therefore cannot be compared with AYbRAH. Evolview v2 (47) was used to map ortholog databases coverage onto the yeast species tree.

Subcellular localization prediction

Subcellular localization predictions for all proteins in the pan-genome were computed with MitoProt II (48), Predotar (49) and TargetP (50). The Phobius web server (51) was used to predict transmembrane domains for all proteins.

Literature references

Literature references for characterized proteins were assigned to FOGs in AYbRAH. Additional references were obtained from paperBLAST (52), UniProt (33), Saccharomyces Genome Database (53), PomBase (54), Candida Genome Database (55) and Aspergillus Genome Database (56).

AYbRAH overview

AYbRAH v0.1 and v0.2.3 database statistics are summarized in Table 2. In total, there are 214 498 protein sequences in the pan-genome for 33 yeasts and fungi; Pezizomycotina fungi were included in the database as an outgroup because they have genes that were present in Proto-Yeast’s ancestor, but subsequently lost. AYbRAH has 187 555 proteins (87% of the pan-proteome) that were assigned to 22 538 FOGs and 18 202 HOGs. Ortholog assignments are available in an Excel spreadsheet, a tab-separated file, orthoXML (57) and a JSON format.

Table 2

AYbRAH ortholog database statistics before and after curation. The initial ortholog assignments were obtained with OrthoMCL and OrthoDB. Additional proteins were annotated using TBLASTN. Ortholog groups for enzymes and small metabolite transporters were manually curated by visual inspection of homolog phylogeny and by identifying ortholog groups with an ETE 3 script (40). Ortholog groups were modified by adding unannotated proteins to existing groups via pplacer (41) or by collapsing multiple ortholog groups into a single ortholog group if there were no gene duplications in the homolog group (under-clustering)

AYbRAH
v0.1v0.2.3
Proteins212 551214 498
Proteins in AYbRAH169 118 (79%)187 555 (87%)
Fungal ortholog groups14 24922 538
Homolog groups018 202
Manually curated ortholog groups0625
Electronically modified ortholog groups03760
AYbRAH
v0.1v0.2.3
Proteins212 551214 498
Proteins in AYbRAH169 118 (79%)187 555 (87%)
Fungal ortholog groups14 24922 538
Homolog groups018 202
Manually curated ortholog groups0625
Electronically modified ortholog groups03760
Table 2

AYbRAH ortholog database statistics before and after curation. The initial ortholog assignments were obtained with OrthoMCL and OrthoDB. Additional proteins were annotated using TBLASTN. Ortholog groups for enzymes and small metabolite transporters were manually curated by visual inspection of homolog phylogeny and by identifying ortholog groups with an ETE 3 script (40). Ortholog groups were modified by adding unannotated proteins to existing groups via pplacer (41) or by collapsing multiple ortholog groups into a single ortholog group if there were no gene duplications in the homolog group (under-clustering)

AYbRAH
v0.1v0.2.3
Proteins212 551214 498
Proteins in AYbRAH169 118 (79%)187 555 (87%)
Fungal ortholog groups14 24922 538
Homolog groups018 202
Manually curated ortholog groups0625
Electronically modified ortholog groups03760
AYbRAH
v0.1v0.2.3
Proteins212 551214 498
Proteins in AYbRAH169 118 (79%)187 555 (87%)
Fungal ortholog groups14 24922 538
Homolog groups018 202
Manually curated ortholog groups0625
Electronically modified ortholog groups03760

The AYbRAH web portal

AYbRAH has a web page for each HOG with information on gene names, descriptions, gene origin (paralog, ohnolog and xenolog), literature references, localization predictions and phylogenetic reconstruction. A sample webpage for the acetyl-CoA synthetase can be seen in Supplementary Information. Protein families can be searched by FOG (FOG00404) or HOG (HOG00229) identification codes, gene names (ACS1), ordered locus (YAL054C), UniProt entry names (ACS1_YEAST) or protein accession codes from UniProt (Q01574), NCBI RefSeq (NP_009347.1) or EMBL (CAA47054.1).

A sample phylogenetic tree rendered by ETE v3 (40) and descriptions of its annotation features is shown in Figure 3 for the acetyl-CoA synthetase family (HOG00229). The initial ortholog assignments by OrthoMCL did not distinguish between the ACS1 (FOG00404) and ACS2 (FOG00405) paralogs. From this phylogeny, we can see that ACS2 arose from a duplication from ACS1, because the basal species (Rhodotorula graminis, Schizosaccharomyces pombe, Pezizomycotina fungi) do not have ACS2, and the ACS2 subtree has high bootstrap support (79%). Therefore, ACS1 is the parent ortholog group to ACS2. This multi-level hierarchical relationship for ortholog groups was adopted in AYbRAH and was recently recommended by (58); current ortholog databases and Clusters of Orthologous Groups (COGs) collections treat these ortholog groups as equal or siblings. Discrepancies in ortholog assignments can be identified by comparing bootstrap support values for subtrees and ortholog assignments, as was done with ACS1 and ACS2. Issues may be reported on GitHub or pull requests can be initiated for large changes to ortholog groups.

Snapshots for mitochondrial localization and transmembrane domain predictions are shown in Figures 4A and B for internal alternative NADH dehydrogenase, encoded by NDI1 (FOG00846). Reviewing localization predictions for orthologous proteins with multiple algorithms enables researchers to make prudent decisions about protein localization, rather than relying on one method for one protein sequence. For example, Cybja1_131289 encodes internal alternative NADH dehydrogenase, yet its mitochondrial localization probability is 0.0019 with MitoProt II; all other mitochondrial predictions for Ndi1p orthologs are greater than 0.80 with MitoProt II. A review of the upstream nucleotide sequence of Cybja1_131289 indicates additional start codons that were not included in the protein annotation. MitoProt II predicts a mitochondrial localization probability of 0.5191 for the full protein sequence, which is more consistent with its orthologs.

AYbRAH curation

OrthoMCL and OrthoDB are less computationally intensive than phylogenetic-based methods, but they are not always accurate (59). Curation was required to resolve incorrect ortholog assignments due to over-clustering and under-clustering.

Over-clustering by OrthoMCL

Over-clustering has been described in past studies (60), which occurs when graph-based methods create ortholog groups that do not distinguish between orthologs and paralogs. Over-clustering by OrthoMCL was common in gene families with many duplications or high sequence similarities, such as the aldehyde dehydrogenase (HOG00216) and the major facilitator superfamily (HOG01031); adjusting parameters for BLASTP and OrthoMCL did not help differentiate between orthologs and paralogs in HOG00216 and neither did adding more proteomes to the OrthoMCL pipeline (results not shown). Figure 5 illustrates an example of over-clustering with a subset of the hexokinase family (HOG00193). In this phylogenetic reconstruction, one hexokinase gene was present in the ancestral yeast species, but a gene duplication in Pichiaceae led to the HXK3 paralog; the HXK2 ortholog is subsequently not maintained in O. parapolymorpha’s genome. OrthoMCL assigned the HXK3 paralog to the same ortholog group as HXK2. The RBH method, commonly used for ortholog identification (62), would have also falsely identified O. parapolymorpha’s HXK3 as orthologous to S. cerevisiae’s HXK2. This example highlights how the greediness of graph-based methods can misidentify orthologs, which has been shown for yeast ohnologs (59), and how incorrect ortholog assignments can be made with pairwise comparisons. Paralogs were identified from over-clustered ortholog groups by finding nodes with high bootstrap support in the consensus phylogenetic trees for homologs and migrating the proteins to new ortholog groups; in some cases orthologs were identified by reviewing the sequence alignment of homologs.

Under-clustering by OrthoMCL

Under-clustering occurs when orthologous proteins are assigned to multiple ortholog groups. OrthoMCL was more prone to under-clustering for short protein sequences and proteins with low sequence similarity, such as subunits in the electron transport chain complexes and Flo8p. Figure 6 demonstrates under-clustering with a subset of the Flo8p family that was incorrectly assigned to multiple ortholog groups by OrthoMCL. Under-clustering was mostly resolved via a Python script that coalesced proteins into a new ortholog group when multiple FOGs were present in a HOG yet no organism had any gene duplications.

Comparison of AYbRAH to other ortholog identification methods

BLASTP scoring metrics

BLASTP is used as the basis for many ortholog predictions, including graph-based methods (29) and RBH (62). The distribution of percent identity, log(bit score) and −log(expect value) for proteins identified as orthologs to S. cerevisiae in AYbRAH are shown in Figure 7. Taxonomic groups include the Saccharomycotina outgroup, basal Saccharomycotina, Pichiaceae, CTG clade, Phaffomycetaceae and Saccharomycodaceae and Saccharomycetaceae (Table 1). The approximate divergence time with S. cerevisiae is 400–600 million years with the Saccharomycotina outgroup, 200–400 million years with the basal Saccharomycotina yeasts, 200 million years with Pichiaceae and CTG clades, 100–200 million years with Phaffomycetaceae and Saccharomycodaceae and 0–100 million years with Saccharomycetaceae. The distributions of percent identity, log(bit score), and -log(expect value) for proteins with 100–400 million years of divergence with S. cerevisiae are similar; however, the distributions skew differently for percent identity and -log(expect value) for the Saccharomycotina outgroup (400 million years of divergence) and Saccharomycetaceae (100 million years of divergence). Distributions for percent identity, log(bit score) and -log(expect value) for each species in AYbRAH are shown in Figures S1, S2 and S3. These results highlight the need to use phylogenetic methods and hidden Markov models to identify orthologs over long evolutionary timescales (43), but also enable orthologs to be identified by synteny and sequence similarity over smaller evolutionary time ranges (63, 64).

Distribution of BLASTP percent identities, logarithm of bit scores and negative logarithm of expect values for proteins orthologous to S. cerevisiae. The bottom half of orthologous proteins in the Saccharomycotina outgroup and Saccharomycetaceae has a percent identities of <40% and 58%, respectively; the bottom half of the expect-value ranges is >1e-60 and 1e-125 for the same groups. The wide and skewed distribution in the Saccharomycotina outgroup highlights the difficulty in making pairwise ortholog predictions for proteins with >400 millions of divergence in Dikarya fungi with BLASTP results; however, orthologs can be easily identified in the Saccharomycetaceae family because of their high sequence similarities and low expect values.
Figure 7

Distribution of BLASTP percent identities, logarithm of bit scores and negative logarithm of expect values for proteins orthologous to S. cerevisiae. The bottom half of orthologous proteins in the Saccharomycotina outgroup and Saccharomycetaceae has a percent identities of <40% and 58%, respectively; the bottom half of the expect-value ranges is >1e-60 and 1e-125 for the same groups. The wide and skewed distribution in the Saccharomycotina outgroup highlights the difficulty in making pairwise ortholog predictions for proteins with >400 millions of divergence in Dikarya fungi with BLASTP results; however, orthologs can be easily identified in the Saccharomycetaceae family because of their high sequence similarities and low expect values.

Comparison of AYbRAH to well-established phylogenomic databases

Ortholog assignments in AYbRAH were compared with OMA, PANTHER, HOGENOM, eggNOG and KO (Table 3). OMA and PANTHER have the highest number of congruous ortholog groups with AYbRAH. Interestingly, PANTHER tends to over-cluster protein sequences into ortholog groups, while OMA tends to under-cluster. HOGENOM, eggNOG and KO have a high fraction of proteins not assigned to any ortholog groups, indicating that AYbRAH is able to identify more ortholog groups with OrthoMCL and OrthoDB.

Table 3

Comparison of ortholog assignments between AYbRAH and well-established phylogenomic databases. OMA and PANTHER are the most congruous with AYbRAH. Bold numbers indicate the greatest source of incongruency with AYbRAH. OMA and PANTHER are predicted to have more under-clustered and over-clustered groups relative to AYbRAH, respectively. HOGENOM, eggNOG and KO have a large number of proteins with no ortholog assignment

Ortholog DatabaseFOGs comparedCongruent groupsOver-clustered groupsUnder-clustered groupsOver and under-clustered groupsNo ortholog group assignment
OMA850559%5%19%3%14%
PANTHER701458%29%1%4%8%
HOGENOM939350%14%11%1%24%
eggNOG782748%10%4%1%37%
KO902722%16%0%0%62%
Ortholog DatabaseFOGs comparedCongruent groupsOver-clustered groupsUnder-clustered groupsOver and under-clustered groupsNo ortholog group assignment
OMA850559%5%19%3%14%
PANTHER701458%29%1%4%8%
HOGENOM939350%14%11%1%24%
eggNOG782748%10%4%1%37%
KO902722%16%0%0%62%
Table 3

Comparison of ortholog assignments between AYbRAH and well-established phylogenomic databases. OMA and PANTHER are the most congruous with AYbRAH. Bold numbers indicate the greatest source of incongruency with AYbRAH. OMA and PANTHER are predicted to have more under-clustered and over-clustered groups relative to AYbRAH, respectively. HOGENOM, eggNOG and KO have a large number of proteins with no ortholog assignment

Ortholog DatabaseFOGs comparedCongruent groupsOver-clustered groupsUnder-clustered groupsOver and under-clustered groupsNo ortholog group assignment
OMA850559%5%19%3%14%
PANTHER701458%29%1%4%8%
HOGENOM939350%14%11%1%24%
eggNOG782748%10%4%1%37%
KO902722%16%0%0%62%
Ortholog DatabaseFOGs comparedCongruent groupsOver-clustered groupsUnder-clustered groupsOver and under-clustered groupsNo ortholog group assignment
OMA850559%5%19%3%14%
PANTHER701458%29%1%4%8%
HOGENOM939350%14%11%1%24%
eggNOG782748%10%4%1%37%
KO902722%16%0%0%62%

Ten ortholog groups were randomly selected from the over-clustered groups in PANTHER and under-clustered groups in OMA to determine the source of the incongruency. It was found that 3 of the ten over-clustered ortholog groups in PANTHER were correctly annotated in AYbRAH, 1 ortholog group was correctly identified in PANTHER but under-clustered in AYbRAH, 1 ortholog group was not correctly identified in either database and 5 ortholog groups required further curation since the phylogenies are ambiguous. All ten ortholog groups from OMA were under-clustered, suggesting a systematic bias to not cluster proteins with lower sequence similarity; i.e., proteins identified as orthologous in AYbRAH were separated into two or more ortholog groups in OMA. Therefore, the PANTHER database is most closely aligned with AYbRAH. All other databases appear to be more prone to over-clustering or not have any annotation.

Orthology is inherently defined by phylogeny (65, 66). Clustering-based methods are well suited to cluster proteins into homolog groups, but it is not clear how these methods can properly identify orthologous proteins with one-dimensional sequence similarity alone, or identify xenologs without knowledge of a species tree. In our experience adding more diverse proteomes to OrthoMCL did not improve differentiation between orthologs and paralogs. PANTHER had a higher accuracy than other phylogenomic databases in our comparison with AYbRAH, despite PANTHER having fewer proteomes in its pan-genome. This is likely an outcome of its phylogenetic reconstruction of PANTHER families and its continued curation for two decades. Therefore, future methods should consider mapping new proteomes to existing databases, such as eggNOG-mapper (67) and TreeGrafter (68), rather than recomputing ortholog assignments, but also have a component of community curation.

Applications of a curated ortholog database

Ortholog databases offer additional benefits beyond simply identifying orthologous proteins. These databases can be used to identify gene targets for functional characterization to functional genome annotation to streamlining GENRE; Galperin et al. (58) recently outlined some of the benefits and challenges to ortholog databases for microbial genomics. First, a curated ortholog database can serve as a repository for orthologs that have been screened and orthologs that require screening (69). Rather than characterizing all the orthologs in a handful of model organisms, research communities can broaden their efforts to understand the orthologs that do not exist in model organisms and the set of orthologs that do not have a conserved function with orthologs in model organisms. Second, a curated ortholog database can be used to improve and simplify genome annotation (69). Genes from newly sequenced organisms can be mapped to curated ortholog groups rather than using protein sequences from ortholog databases as queries in TBLASTN searches (70). New ortholog groups can be created for de novo genes or genes from recent duplications. Pulling annotations from a curated ortholog database has the advantage of unifying the names and descriptions of genes between organisms, as has been proposed for ribosomal subunits (71), and can reduce the number of genes that are misannotated or annotated as conserved hypothetical proteins. Finally, a curated ortholog database can be used to improve the quality and quantity of GENREs. GENREs inherently require a great deal of curation to identify orthologous proteins and their function, which is often not transparent. Refocusing this effort to curate ortholog groups and their function in open-source knowledgebase for pan-genomes can allow for improvements to be pushed to all GENREs, and for GENREs to be compiled for any taxonomic level, from kingdom to strain.

Future plans for AYbRAH

Integration with PANTHER

OrthoDB was chosen to cluster ortholog groups in AYbRAH into homolog groups because it spans more taxa than other phylogenomic databases and has ortholog assignments for different taxonomic ranks; however, it is less specific than PANTHER, despite the latter only having a few fungal proteome annotations. Future updates to AYbRAH will migrate the AYbRAH homolog group backbone from OrthoDB to PANTHER, and add the remaining fungi in PANTHER to increase its phylogenomic span. These include other fungal model organisms, fungi and yeasts having pathogenicity to humans or plants or fungi and yeasts occupying the following important taxonomic ranks: Batrachochytrium dendrobatidis, Cryptococcus neoformans, Puccinia graminis, Ustilago maydis, Emericella nidulans, Neosartorya fumigata, Phaeosphaeria nodorum, Sclerotinia sclerotiorum, Candida albicans and Eremothecium gossypii.

Reconciling AYbRAH with YGOB and CGOB

The Yeast Gene Order Browser (YGOB) (63) and Candida Gene Order Browser (CGOB) (72) are the gold standard for ortholog databases in yeast genomics and were created using sequence similarity and synteny. YGOB and CGOB span roughly 112 and 239 million years of evolution, respectively, while AYbRAH spans 600 million years of evolution (2). Although AYbRAH has a broader pan-genomic coverage, YGOB and CGOB are expected to have better paralog and ohnolog assignments than AYbRAH because of its use of synteny. Future versions of AYbRAH will be reconciled with YGOB and CGOB.

Coordinate-based protein annotations

It has been noted that genome protein annotations sometimes contain inaccuracies (72). For example, the protein translation Cybja1_131289 does not include its full N-terminal sequence. Another surprising shortfall of some genome annotations are genes that do not have any annotation. Spathaspora passildarium’s genome encodes have two PHO3 homologs in tandem, but only one protein is currently annotated. AYbRAH will adopt the genomic coordinate-based system used in YGOB and CGOB (72) to improve protein annotations.

Conclusion

In conclusion, we developed AYbRAH as an open-source ortholog database for yeasts and fungi because existing phylogenomic databases do not span diverse yeasts and sometimes cannot distinguish between orthologs, paralogs and xenologs. Manual curation was required for gene families with high sequence similarity, often arising from recent gene duplications, and with gene families with low sequence similarity. Curated ortholog databases can be implemented for other taxa to improve their genome annotations using PANTHER and other tree-based methods.

Abbreviations

(AYbRAH) Analyzing Yeasts by Reconstructing Ancestry of Homologs; (CGOB) Candida Gene Order Browser; (COG) Clusters of Orthologous Groups; (FOG) Fungal Ortholog Group; (GENRE) Genome-scale Network REconstruction; (HOG) Homolog Group; (YGOB) Yeast Gene Order Browser; (RBH) Reciprocal Best Hit

Availability of data

AYbRAH database files and additional files, such as phylogenetic trees and sequence alignments, can be found at https://github.com/LMSE/aybrah.

Acknowledgements

The authors gratefully acknowledge Prof. Belinda Chang and Ryan Schott for their advice with the phylogenetic analysis and Dean Robson for his help implementing the search function in the AYbRAH web portal.

Funding

NSERC Bioconversion Network, Industrial Biocatalysis Network, Genome Canada, Ontario Ministry of Research and Innovation, and NSERC CREATE M3 (to K.C.).

Conflict of interest. None declared.

Database URL: https://lmse.github.io/aybrah/

References

1

Kurtzman
,
C.
,
Fell
,
J.W.
and
Boekhout
,
T.
(
2011
)
The Yeasts: A Taxonomic Study
.
Elsevier
.
Burlington, MA, USA
.

2

Hedges
,
S.B.
,
Marin
,
J.
,
Suleski
,
M.
et al.  (
2015
)
Tree of life reveals clock-like speciation and diversification
.
Mol. Biol. Evol.
32
,
835
845
.

3

Aiba
,
S.
and
Matsuoka
,
M.
(
1979
)
Identification of metabolic model: citrate production from glucose by Candida lipolytica
.
Biotechnol. Bioeng.
,
21
,
1373
1386
.

4

Boulton
,
C.A.
and
Ratledge
,
C.
(
1983
)
Use of transition studies in continuous cultures of Lipomyces starkeyi, an oleaginous yeast, to investigate the physiology of lipid accumulation
.
Microbiology
,
129
,
2871
2876
.

5

Banat
,
I.
,
Singh
,
D.
and
Marchant
,
R.
(
1996
)
The use of a thermotolerant fermentative Kluyveromyces marxianus IMB3 yeast strain for ethanol production
.
Acta Biotechnologica
,
16
,
215
223
.

6

Ryabova
,
O.B.
,
Chmil
,
O.M.
and
Sibirny
,
A.A.
(
2003
)
Xylose and cellobiose fermentation to ethanol by the thermotolerant methylotrophic yeast Hansenula polymorpha
.
FEMS Yeast Res.
,
4
,
157
164
.

7

Rush
,
B.J.
and
Fosmer
,
A.M.
(
2013
)
Methods for succinate production
.
US Patent App.
14
/374
,
464
.

8

Lindberg
,
L.
,
Santos
,
A.X.
,
Riezman
,
H.
et al.  (
2013
)
Lipidomic profiling of Saccharomyces cerevisiae and Zygosaccharomyces bailii reveals critical changes in lipid composition in response to acetic acid stress
.
PLoS One
,
8
,
e73936
.

9

Hall
,
C.
,
Brachat
,
S.
and
Dietrich
,
F.S.
(
2005
)
Contribution of horizontal gene transfer to the evolution of Saccharomyces cerevisiae
.
Eukaryot. Cell
,
4
,
1102
1115
.

10

Larsson
,
C.
and
Gustafsson
,
L.
(
1993
)
The role of physiological state in osmotolerance of the salt-tolerant yeast Debaryomyces hansenii
.
Can. J. Microbiol.
,
39
,
603
609
.

11

Schneider
,
H.
,
Wang
,
P.
,
Chan
,
Y.
et al.  (
1981
)
Conversion of D-xylose into ethanol by the yeast Pachysolen tannophilus
.
Biotechnol. Lett.
,
3
,
89
92
.

12

Slininger
,
P.
,
Bothast
,
R.
,
van Cauwenberge
,
J.
et al.  (
1982
)
Conversion of D-xylose to ethanol by the yeast Pachysolen tannophilus
.
Biotechnol. Bioeng.
,
24
,
371
384
.

13

Toivola
,
A.
,
Yarrow
,
D.
,
Van Den Bosch
,
E.
et al.  (
1984
)
Alcoholic fermentation of D-xylose by yeasts
.
Appl. Environ. Microbiol.
,
47
,
1221
1223
.

14

Mühlhausen
,
S.
,
Findeisen
,
P.
,
Plessmann
,
U.
et al.  (
2016
)
A novel nuclear genetic code alteration in yeasts and the evolution of codon reassignment in eukaryotes
.
Genome Res.
,
26
,
945
955
.

15

Sousa
,
M.J.
,
Rodrigues
,
F.
,
Coôrte-Real
,
M.
et al.  (
1998
)
Mechanisms underlying the transport and intracellular metabolism of acetic acid in the presence of glucose in the yeast Zygosaccharomyces bailii
.
Microbiology
,
144
,
665
670
.

16

de Deken
,
R
. (
1966
)
The Crabtree effect: a regulatory system in yeast
.
Microbiology
,
44
,
149
156
.

17

van Urk
,
H.
,
Voll
,
W.L.
,
Scheffers
,
W.A.
et al.  (
1990
)
Transient-state analysis of metabolic fluxes in Crabtree-positive and Crabtree-negative yeasts
.
Appl. Environ. Microbiol.
,
56
,
281
287
.

18

Blank
,
L.M.
,
Lehmbeck
,
F.
and
Sauer
,
U.
(
2005
)
Metabolic-flux and network analysis in fourteen hemiascomycetous yeasts
.
FEMS Yeast Res.
,
5
,
545
558
.

19

Christen
,
S.
and
Sauer
,
U.
(
2011
)
Intracellular characterization of aerobic glucose metabolism in seven yeast species by 13C flux analysis and metabolomics
.
FEMS Yeast Res.
,
11
,
263
272
.

20

Koonin
,
E.V.
(
2001
)
An apology for orthologs-or brave new memes
.
Genome Biol.
,
2
,
1005
1
.

21

Jensen
,
R.A.
(
2001
)
Orthologs and paralogs - we need to get it right
.
Genome Biol.
,
2
,
1002
1
.

22

Koonin
,
E.V.
(
2005
)
Orthologs, paralogs, and evolutionary genomics
.
Annu. Rev. Genet.
,
39
,
309
338
.

23

Thomson
,
J.M.
,
Gaucher
,
E.A.
,
Burgan
,
M.F.
et al.  (
2005
)
Resurrecting ancestral alcohol dehydrogenases from yeast
.
Nat. Genet.
,
37
,
630
635
.

24

Bender
,
T.
,
Pena
,
G.
and
Martinou
,
J.-C.
(
2015
)
Regulation of mitochondrial pyruvate uptake by alternative pyruvate carrier complexes
.
EMBO J.
,
34
,
911
924
.

25

Hall
,
C.
,
Brachat
,
S.
and
Dietrich
,
F.S.
(
2005
)
Contribution of horizontal gene transfer to the evolution of Saccharomyces cerevisiae
.
Eukaryot. Cell
,
4
,
1102
1115
.

26

Dujon
,
B.
(
2010
)
Yeast evolutionary genomics
.
Nat. Rev. Genet.
,
11
,
512
524
.

27

Gaucher
,
E.A.
,
Kratzer
,
J.T.
and
Randall
,
R.N.
(
2010
)
Deep phylogeny–how a tree can help characterize early life on Earth
.
Cold Spring Harb. Perspect. Biol.
,
2
,
a002238
.

28

Riley
,
R.
,
Haridas
,
S.
,
Wolfe
,
K.H.
et al.  (
2016
)
Comparative genomics of biotechnologically important yeasts
.
Proc. Natl. Acad. Sci.
,
113
,
9882
9887
.

29

Kuzniar
,
A.
,
van Ham
,
R.C.
,
Pongor
,
S.
et al.  (
2008
)
The quest for orthologs: finding the corresponding gene across genomes
.
Trends Genet.
,
24
,
539
551
.

30

Wohlbach
,
D.J.
,
Kuo
,
A.
,
Sato
,
T.K.
et al.  (
2011
)
Comparative genomics of xylose-fermenting fungi for enhanced biofuel production
.
Proc. Natl. Acad. Sci.
,
108
,
13212
13217
.

31

Papini
,
M.
,
Nookaew
,
I.
,
Uhlén
,
M.
et al.  (
2012
)
Scheffersomyces stipitis
:
a comparative systems biology study with the Crabtree positive yeast Saccharomyces cerevisiae
.
Microb. Cell Fact.
,
11
,
1
.

32

Caspeta
,
L.
,
Shoaie
,
S.
,
Agren
,
R.
et al.  (
2012
)
Genome-scale metabolic reconstructions of Pichia stipitis and Pichia pastoris and in silico evaluation of their potentials
.
BMC Syst. Biol.
,
6
,
1
.

33

Uniprot Consortium
(
2014
)
UniProt: a hub for protein information
.
Nucleic Acids Res.
,
43
,
D204
D212
.

34

Grigoriev
,
I.V.
,
Nikitin
,
R.
,
Haridas
,
S.
et al.  (
2013
)
MycoCosm portal: gearing up for 1000 fungal genomes
.
Nucleic Acids Res.
,
42
,
699
704
.

35

Fischer
,
S.
,
Brunk
,
B.P.
,
Chen
,
F.
et al.  (
2011
)
Using OrthoMCL to assign proteins to OrthoMCL-DB groups or to cluster proteomes into new ortholog groups
.
Curr. Protoc. Bioinformatics
,
6
,
1
9
.

36

Kriventseva
,
E.V.
,
Tegenfeldt
,
F.
,
Petty
,
T.J.
et al.  (
2015
)
OrthoDB v8: update of the hierarchical catalog of orthologs and the underlying free software
.
Nucleic Acids Res.
,
43
,
D250
D256
.

37

Katoh
,
K.
and
Standley
,
D.M.
(
2013
)
MAFFT multiple sequence alignment software version 7: improvements in performance and usability
.
Mol. Biol. Evol.
,
30
,
772
780
.

38

Guindon
,
S.
,
Dufayard
,
J.-F.
,
Lefort
,
V.
et al.  (
2010
)
New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0
.
Syst. Biol.
,
59
,
307
321
.

39

Sukumaran
,
J.
and
Holder
,
M.T.
(
2010
)
DendroPy: a python library for phylogenetic computing
.
Bioinformatics
,
26
,
1569
1571
.

40

Huerta-Cepas
,
J.
,
Serra
,
F.
and
Bork
,
P.
(
2016a
)
ETE 3: reconstruction, analysis, and visualization of phylogenomic data
.
Mol. Biol. Evol.
,
33
,
1635
1638
.

41

Matsen
,
F.
,
Kodner
,
R.
and
Armbrust
,
E.
(
2010
)
pplacer: linear time maximum-likelihood and Bayesian phylogenetic placement of sequences onto a fixed reference tree
.
BMC Bioinformatics
,
11
,
538
.

42

Altenhoff
,
A.M.
,
Škunca
,
N.
,
Glover
,
N.
et al.  (
2015
)
The OMA orthology database in 2015: function predictions, better plant support, synteny view and other improvements
.
Nucleic Acids Res.
,
43
,
D240
D249
.

43

Mi
,
H.
,
Huang
,
X.
,
Muruganujan
,
A.
et al.  (
2016
)
PANTHER version 11: expanded annotation data from Gene Ontology and Reactome pathways, and data analysis tool enhancements
.
Nucleic Acids Res.
,
45
,
D183
D189
.

44

Penel
,
S.
,
Arigon
,
A.-M.
,
Dufayard
,
J.-F.
et al.  (
2009
)
Databases of homologous gene families for comparative genomics
.
BMC Bioinformatics
,
10
,
S3
.

45

Huerta-Cepas
,
J.
,
Szklarczyk
,
D.
,
Forslund
,
K.
et al.  (
2016b
)
eggNOG 4.5: a hierarchical orthology framework with improved functional annotations for eukaryotic, prokaryotic and viral sequences
.
Nucleic Acids Res.
,
44
,
D286
D293
.

46

Mao
,
X.
,
Cai
,
T.
,
Olyarchuk
,
J.G.
et al.  (
2005
)
Automated genome annotation and pathway identification using the KEGG orthology (KO) as a controlled vocabulary
.
Bioinformatics
,
21
,
3787
3793
.

47

He
,
Z.
,
Zhang
,
H.
,
Gao
,
S.
et al.  (
2016
)
Evolview v2: an online visualization and management tool for customized and annotated phylogenetic trees
.
Nucleic Acids Res.
,
W236
W241
.

48

Claros
,
M.G.
and
Vincens
,
P.
(
1996
)
Computational method to predict mitochondrially imported proteins and their targeting sequences
.
Eur. J. Biochem.
,
241
,
779
786
.

49

Small
,
I.
,
Peeters
,
N.
,
Legeai
,
F.
et al.  (
2004
)
Predotar: a tool for rapidly screening proteomes for N-terminal targeting sequences
.
Proteomics
,
4
,
1581
1590
.

50

Emanuelsson
,
O.
,
Brunak
,
S.
,
von Heijne
,
G.
et al.  (
2007
)
Locating proteins in the cell using TargetP, SignalP and related tools
.
Nat. Protoc.
,
2
,
953
971
.

51

Käll
,
L.
,
Krogh
,
A.
and
Sonnhammer
,
E.L.
(
2007
)
Advantages of combined transmembrane topology and signal peptide prediction—the Phobius web server
.
Nucleic Acids Res.
,
35
,
W429
W432
.

52

Price
,
M.N.
and
Arkin
,
A.P.
(
2017
)
PaperBLAST: text mining papers for information about homologs
.
mSystems
,
2
,
e00039
e00017
.

53

Cherry
,
J.M.
,
Hong
,
E.L.
,
Amundsen
,
C.
et al.  (
2011
)
Saccharomyces
Genome Database: the genomics resource of budding yeast
.
Nucleic Acids Res.
,
40
,
D700
D705
.

54

McDowall
,
M.D.
,
Harris
,
M.A.
,
Lock
,
A.
et al.  (
2014
)
PomBase 2015: updates to the fission yeast database
.
Nucleic Acids Res.
,
43
,
D656
D661
.

55

Inglis
,
D.O.
,
Arnaud
,
M.B.
,
Binkley
,
J.
et al.  (
2011
)
The Candida genome database incorporates multiple Candida species: multispecies search and analysis tools with curated gene and protein information for Candida albicans and Candida glabrata
.
Nucleic Acids Res.
,
40
,
D667
D674
.

56

Cerqueira
,
G.C.
,
Arnaud
,
M.B.
,
Inglis
,
D.O.
et al.  (
2013
)
The Aspergillus Genome Database: multispecies curation and incorporation of RNA-Seq data to improve structural gene annotations
.
Nucleic Acids Res.
,
42
,
D705
D710
.

57

Schmitt
,
T.
,
Messina
,
D.N.
,
Schreiber
,
F.
et al.  (
2011
)
Letter to the editor: SeqXML and Orthoxml_ standards for sequence and orthology information
.
Brief. Bioinform.
,
12
,
485
488
.

58

Galperin
,
M.Y.
,
Kristensen
,
D.M.
,
Makarova
,
K.S.
et al.  (
2017
)
Microbial genome analysis: the COG approach
.
Brief. Bioinform.
, https://academic.oup.com/bib/advance-article-abstract/doi/10.1093/bib/bbx117/4158183.

59

Salichos
,
L.
and
Rokas
,
A.
(
2011
)
Evaluating ortholog prediction algorithms in a yeast model clade
.
PLoS One
,
6
,
e18755
.

60

Jothi
,
R.
,
Zotenko
,
E.
,
Tasneem
,
A.
et al.  (
2006
)
COCO-CL: hierarchical clustering of homology relations based on evolutionary correlations
.
Bioinformatics
,
22
,
779
788
.

61

Remm
,
M.
,
Storm
,
C.E.
and
Sonnhammer
,
E.L.
(
2001
)
Automatic clustering of orthologs and in-paralogs from pairwise species comparisons
.
J. Mol. Biol.
,
314
,
1041
1052
.

62

Moreno-Hagelsieb
,
G.
and
Latimer
,
K.
(
2008
)
Choosing BLAST options for better detection of orthologs as reciprocal best hits
.
Bioinformatics
,
24
,
319
324
.

63

Byrne
,
K.P.
and
Wolfe
,
K.H.
(
2005
)
The Yeast Gene Order Browser: combining curated homology and syntenic context reveals gene fate in polyploid species
.
Genome Res.
,
15
,
1456
1461
.

64

Scannell
,
D.R.
,
Byrne
,
K.P.
,
Gordon
,
J.L.
et al.  (
2006
)
Multiple rounds of speciation associated with reciprocal gene loss in polyploid yeasts
.
Nature
,
440
,
341
.

65

Fitch
,
W.M.
(
1970
)
Distinguishing homologous from analogous proteins
.
Syst. Zool.
,
19
,
99
113
.

66

Gabaldón
,
T.
(
2008
)
Large-scale assignment of orthology: back to phylogenetics?
Genome Biol.
,
9
,
235
.

67

Huerta-Cepas
,
J.
,
Forslund
,
K.
,
Pedro Coelho
,
L.
et al.  (
2017
)
Fast genome-wide functional annotation through orthology assignment by eggNOG-Mapper
.
Mol. Biol. Evol.
,
34
,
2115
2122
.

68

Tang
,
H.
,
Finn
,
R.D.
and
Thomas
,
P.D.
(
2018
)
TreeGrafter: phylogenetic tree-based annotation of proteins with Gene Ontology terms and other annotations
.
Bioinformatics
,
bty625
.

69

Galperin
,
M.Y.
and
Koonin
,
E.V.
(
2004
)
‘Conserved hypothetical’ proteins: prioritization of targets for experimental study
.
Nucleic Acids Res.
,
32
,
5452
5463
.

70

Proux-Wéra
,
E.
,
Armisén
,
D.
,
Byrne
,
K.P.
et al.  (
2012
)
A pipeline for automated annotation of yeast genome sequences by a conserved-synteny approach
.
BMC Bioinformatics
,
13
,
237
.

71

Ban
,
N.
,
Beckmann
,
R.
,
Cate
,
J.H.
et al.  (
2014
)
A new system for naming ribosomal proteins
.
Curr. Opin. Struct. Biol.
,
24
,
165
169
.

72

Maguire
,
S.L.
,
OhÉigeartaigh
,
S.S.
,
Byrne
,
K.P.
et al.  (
2013
)
Comparative genome analysis and gene finding in Candida species using CGOB
.
Mol. Biol. Evol.
,
30
,
1281
1291
.

73

Firrincieli
,
A.
,
Otillar
,
R.
,
Salamov
,
A.
et al.  (
2015
)
Genome sequence of the plant growth promoting endophytic yeast Rhodotorula graminis WP1
.
Front. Microbiol.
,
6
,
978
.

74

Wood
,
V.
,
Gwilliam
,
R.
,
Rajandream
,
M.-A.
et al.  (
2002
)
The genome sequence of Schizosaccharomyces pombe
.
Nature
,
415
,
871
.

75

Pel
,
H.J.
,
de Winde
,
J.H.
,
Archer
,
D.B.
et al.  (
2007
)
Genome sequencing and analysis of the versatile cell factory Aspergillus niger CBS 513.88
.
Nat. Biotechnol.
,
25
,
221
.

76

Galagan
,
J.E.
,
Calvo
,
S.E.
,
Borkovich
,
K.A.
et al.  (
2003
)
The genome sequence of the filamentous fungus Neurospora crassa
.
Nature
,
422
,
859
.

77

Martinez
,
D.
,
Berka
,
R.M.
,
Henrissat
,
B.
et al.  (
2008
)
Genome sequencing and analysis of the biomass-degrading fungus Trichoderma reesei (syn. Hypocrea jecorina)
.
Nat. Biotechnol.
,
26
,
553
.

78

Dujon
,
B.
,
Sherman
,
D.
,
Fischer
,
G.
et al.  (
2004
)
Genome evolution in yeasts
.
Nature
,
430
,
35
.

79

Kunze
,
G.
,
Gaillardin
,
C.
,
Czernicka
,
M.
et al.  (
2014
)
The complete genome of Blastobotrys (Arxula) adeninivorans LS3-a yeast of biotechnological interest
.
Biotechnol. Biofuels
,
7
,
66
.

80

De Schutter
,
K.
,
Lin
,
Y.-C.
,
Tiels
,
P.
et al.  (
2009
)
Genome sequence of the recombinant protein production host Pichia pastoris
.
Nat. Biotechnol.
,
27
,
561
.

81

Morales
,
L.
,
Noel
,
B.
,
Porcel
,
B.
et al.  (
2013
)
Complete DNA sequence of Kuraishia capsulata illustrates novel genomic features among budding yeasts (Saccharomycotina)
.
Genome Biol. Evol.
,
5
,
2524
2539
.

82

Piškur
,
J.
,
Ling
,
Z.
,
Marcet-Houben
,
M.
et al.  (
2012
)
The genome of wine yeast Dekkera bruxellensis provides a tool to explore its food-related properties
.
Int. J. Food Microbiol.
,
157
,
202
209
.

83

Ravin
,
N.V.
,
Eldarov
,
M.A.
,
Kadnikov
,
V.V.
et al.  (
2013
)
Genome sequence and analysis of methylotrophic yeast Hansenula polymorpha DL1
.
BMC Genomics
,
14
,
837
.

84

Xiao
,
H.
,
Shao
,
Z.
,
Jiang
,
Y.
et al.  (
2014
)
Exploiting Issatchenkia orientalis SD108 for succinic acid production
.
Microb. Cell Fact.
,
13
,
121
.

85

Butler
,
G.
,
Rasmussen
,
M.D.
,
Lin
,
M.F.
et al.  (
2009
)
Evolution of pathogenicity and sexual reproduction in eight Candida genomes
.
Nature
,
459
,
657
.

86

Jeffries
,
T.W.
,
Grigoriev
,
I.V.
,
Grimwood
,
J.
et al.  (
2007
)
Genome sequence of the lignocellulose-bioconverting and xylose-fermenting yeast Pichia stipitis
.
Nat. Biotechnol.
,
25
,
319
326
.

87

Souciet
,
J.-L.
,
Dujon
,
B.
,
Gaillardin
,
C.
et al.  (
2009
)
Comparative genomics of protoploid Saccharomycetaceae
.
Genome Res.
19
,
1696
1709
.

88

Goffeau
,
A.
,
Barrell
,
B.G.
,
Bussey
,
H.
et al.  (
1996
)
Life with 6000 genes
.
Science
,
274
,
546
567
.

89

Scannell
,
D.R.
,
Frank
,
A.C.
,
Conant
,
G.C.
et al.  (
2007
)
Independent sorting-out of thousands of duplicated gene pairs in two yeast species descended from a whole-genome duplication
.
Proc. Natl. Acad. Sci.
,
104
,
8397
8402
.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.