Genome Biology and Evolution Advance Access published July 22, 2014 doi:10.1093/gbe/evu152 Molecular Phylogeny of Sequenced Saccharomycetes Reveals Polyphyly of the Alternative Yeast Codon Usage Stefanie Mühlhausen1, Martin Kollmar1§ 1 Group Systems Biology of Motor Proteins, Department of NMR-based Structural § Corresponding author, Tel.: +49-551-2012260, Fax: +49-551-2012202 Email addresses: SM: [email protected] MK: [email protected] © The Author(s) 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. 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. -1- Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 Biology, Max-Planck-Institute for Biophysical Chemistry, Göttingen, Germany Abstract The universal genetic code defines the translation of nucleotide triplets, called codons, into amino acids. In many Saccharomycetes a unique alteration of this code affects the translation of the CUG codon, which is normally translated as leucine. Most of the species encoding CUG alternatively as serine belong to the Candida genus and were grouped into a so-called CTG clade. However, the “Candida genus” is not a monophyletic group and several Candida species are known to use the standard CUG ancestor of the Candida, or to several branches independently leading to a polyphyletic alternative yeast codon usage (AYCU). In order to resolve the monophyly or polyphyly of the AYCU, we performed a phylogenomics analysis of 26 motor and cytoskeletal proteins from 60 sequenced yeast species. By investigating the CUG codon positions with respect to sequence conservation at the respective alignment positions we were able to unambiguously assign the standard code or AYCU. Quantitative analysis of the highly conserved leucine and serine alignment positions showed, that 61.1% and 17% of the CUG codons coding for leucine and serine, respectively, are at highly conserved positions, while only 0.6% and 2.3% of the CUG codons, respectively, are at positions conserved in the respective other amino acid. Plotting the codon usage onto the phylogenetic tree revealed the polyphyly of the AYCU with Pachysolen tannophilus and the CTG clade branching independently within a time span of 30 to 100 Mya. Keywords Genetic code, codon reassignment, codon usage, evolution, Candida -2- Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 translation. The codon identity could have been changed in a single branch, the Introduction The standard genetic code has been altered in many organisms. In eukaryotes, natural alterations have been identified in mitochondria, in which the universal stop codon UGA is translated as tryptophane (Yokobori et al. 2001; Watanabe & Yokobori 2011), in ciliates, in which the function of the stop codons UAA and UGA was changed to code for glutamine (Tourancheau et al. 1995), and in some Candida yeasts codon CUG translation was changed to serine. In taxonomic terms, the definition of the genus Candida is not very specific. Candida species comprise a group of about 850 organisms (Robert et al. 2013), which can be distantly related as for example Candida glabrata, Ca. albicans, and Ca. caseinolytica. Although it has been claimed that almost all Candida species, with the exceptions of Ca. glabrata and Ca. krusei, belong to a single clade (Butler et al. 2009), which is commonly referred to as “Candida clade”, many analyses showed that Candida species are spread over the entire Saccharomycetes clade (Diezmann et al. 2004; Kurtzman & Suzuki 2010; Kurtzman et al. 2011; Kurtzman 2011). In addition, different names have been given to the same yeast species in the past depending on the reproductive stage, in which they had been identified. For example, the names Meyerozyma guilliermondii (teleomorph), Ca. guilliermondii (anamorph), Yamadazyma guilliermondii and Pichia guilliermondii are used for the same holomorph (whole fungus, including anamorph and teleomorph). -3- Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 (Ohama et al. 1993; Sugita & Nakase 1999; Santos et al. 2011), in which the leucine Because Ca. albicans is the most prominent representative of the genus Candida and was the first shown to use the alternative CUG encoding, the term “Candida clade” is often used synonymous to “CTG clade” and most studies about CUG encoding concentrate on Ca. albicans and closely related Candida species. Therefore, depending on which node is used for defining the “Candida clade”, several species of the branch containing Ca. albicans might decode CUG as leucine and many species outside the “Candida clade” are also named Candida. Similarly, there are many nonCandida species decoding CUG as serine. In order to resolve this ambiguity, one of single branch or to several branches independently. The latter scenario would be in agreement with the proposed timing of the code alteration (~ 270 Ma) and the split of the Saccharomyces and Candida clades about 170 Ma (Massey et al. 2003; Santos et al. 2011) implying that the CUG codon was highly ambiguous for approximately 100 Myr. Species that branched within this time span could have adopted either of the two decoding possibilities. Most phylogenetic analyses of yeast species used 28s rRNA, 18s rRNA, actin and translation elongation factor-1α for reconstructing single or combined trees, which is a very successful approach when analyzing hundreds of species (Kurtzman 2011). However, higher accuracy for tree topologies is obtained in phylogenomic studies that use dozens to thousands of concatenated homologs in tree computations (Fitzpatrick et al. 2006). For the latter approach it is important to distinguish orthologs from paralogs. Only a few yeast species have been included in phylogenomic studies so far. In addition to assigning the CUG codon usage by phylogenetic comparisons it has been suggested to base the translation of codon CUG on the presence of Q9 as the -4- Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 the major questions is whether the change of identity of the CUG codon happened to a major ubiquinone (Sugita & Nakase 1999). Ubiquinone (Coenzyme Q, CoQ) is an electron carrier in the respiratory chain and involved in oxidative stress resistance (Szkopińska 2000; Turunen et al. 2004). Most enzymes of the ubiquinone biosynthesis pathway are conserved between Saccharomyces cerevisiae and Ca. albicans. However, the Ca. albicans ubiquinone has a nonaprenyl side chain (Q9) whereas the budding yeast has six isoprene units (Q6). Although a great variety of CoQ types are present in the genus Candida (Schauer & Hanschke 1999), the species analyzed at that time, which were known to decode CUG as serine, lacked galactose Based on these ideas, an extensive study revealed 89 yeast strains to possess coenzyme Q9 (except for Cl. lusitaniae) and to miss galactose in the cell (Sugita & Nakase 1999). However, in vitro translation assays showed 11 of these species to use codon CUG as leucine. The incongruence in many instances of the presence of Q9, absence of galactose, and molecular phylogenetic data has also been reported by others (Kurtzman & Robnett 1997; Diezmann et al. 2004; Kurtzman & Suzuki 2010). Sixty yeast species (not counting strains and data under embargo) have already been sequenced and assembled (http://www.diark.org, visited May 20, 2014). Only half of these genomes, which include both yeasts of the “CTG clade” and of the Saccharomyces clade, have been analyzed in detail yet. In addition, sequencing efforts increasingly focus on genomes of rarely studied species, whose phylogenetic relationships are not clear. For most of these genomes, gene predictions were even done and made available using both the Standard Code and the Alternative Yeast Nuclear Code (AYNC) resulting in mixed datasets. In order to determine whether there was a time span with a CUG codon usage ambiguity or whether the codon usage -5- Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 in the cells and possessed Q9, with the exception of Clavispora lusitaniae (Q8). alteration could be attributed to a single event, we performed a phylogenomic analysis of highly validated protein sequences of 26 members of the actin, actin-related, tubulin, myosin, kinesin, dynein, and actin-capping protein families of the sequenced yeast species. CUG codon translation has been assigned by analyzing the respective amino acids of the sequences in the context of the sequence alignments. Materials and Methods Some fungal actin-related proteins and myosins have been extracted from previously published datasets (Odronitz & Kollmar 2007; Hammesfahr & Kollmar 2012). The sequences were updated based on newer genome assemblies if necessary and corrected for CUG usage. The data for the other proteins and the species not included in the published datasets have essentially been obtained as described (Odronitz & Kollmar 2007). Shortly, the corresponding gene regions have been identified in TBLASTN searches starting with the respective protein sequences of homologs of Saccharomyces cerevisiae. The respective genomic regions were submitted to AUGUSTUS (Stanke & Morgenstern 2005) to obtain gene predictions. However, feature sets are only available for a few species of the Saccharomycetes clade. Therefore, all hits were subsequently manually analyzed at the genomic DNA level. When necessary, gene predictions were corrected by comparison with the homologs already included in the multiple sequence alignments. Especially, the short exons at the 5’ ends of the tubulin and actin genes were not always identified in automatic gene predictions. EST data to help in the annotation are only available for a few species in public databases. -6- Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 Identification and Annotation of the Proteins In the last years, genome-sequencing efforts have been extended from sequencing species from new branches to sequencing closely related organisms. Here, these species include for example Saccharomyces and Eremothecium species. Protein sequences from these closely related species have been obtained by using the crossspecies functionality of WebScipio (Odronitz et al. 2008; Hatje et al. 2011). Nevertheless, for all these genomes TBLASTN searches have been performed. With this strategy, we sought to ensure that we would not miss more divergent protein family homologs, which might have been derived by species-specific inventions or therefore these cross-species gene reconstructions only worked for a few species. All sequence-related data (protein names, corresponding species, sequences, and gene structure reconstructions) and references to genome sequencing centers are available at CyMoBase (http://www.cymobase.org, (Odronitz & Kollmar 2006)). A list of the analyzed species, their abbreviations as used in the alignments and trees, as well as references and accession numbers is also available as supplementary file S1, Supplementary Material online. WebScipio (Hatje et al. 2011; Odronitz et al. 2008) was used to reconstruct the gene structure of each sequence. Throughout this study we use the teleomorph names of all species. Generating the Multiple Sequence Alignment The protein sequences were added to the already existing multiple sequence alignments (Odronitz & Kollmar 2007; Hammesfahr & Kollmar 2012). In detail, we first aligned every newly predicted sequence to its supposed closest relative using ClustalW (Thompson et al. 2002) and added it then to the multiple sequence alignment. During the subsequent sequence validation process, we manually adjusted -7- Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 duplications. However, the Saccharomycetes belong to the fast-evolving species and the obtained alignment by removing wrongly predicted sequence regions and filling gaps. Still, in those sequences derived from low-coverage genomes many gaps remained. To maintain the integrity of exons preceded or followed by gaps, gaps reflecting missing parts of the genomes were added to the multiple sequence alignment. Sequence homologs of the same protein family (e.g., α-tubulin, β-tubulin, and γ-tubulin) were kept in one single protein family alignment. The actin/actinrelated protein alignment therefore contains 11 sub-families, the CapZ alignment 2, the tubulin alignment 3, the myosin alignment 3, and the kinesin alignment 6. Only sequence alignments can be obtained from CyMoBase. Computing Sequence Conservation The residue conservation at alignment positions was calculated with the conservation code toolbox as implemented by (Capra & Singh 2007). Conservation was estimated with the property-entropy method, an entropy measurement refined with respect to chemical properties of amino acids. Scores were calculated with conservation of adjacent amino acids incorporated (window size 3) and not (window size 0). Except for window size and scoring method, standard parameters were used. Unexpectedly, although window sizes were applied, the rest of the alignment still seemed to have an influence, albeit marginal, on the conservation scores. In addition, amino acids given as “X” are replaced by hyphens “-“ by the software, which denote gap positions in the alignment. Gap positions, however, indicate a biological relevant deletion of this position in the sequence and thus have a strong influence on the conservation score. Therefore, conservation estimates were performed on small sections of the concatenated alignment including Schizosaccharomyectes sequences. For each serine -8- Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 single dynein heavy chain homologs have been found in the analyzed species. The and leucine position, the alignment was reduced to this position and 15 adjacent positions in each direction. Sequences with CUG codons at the respective leucine or serine position were removed from the alignment sections. Computing and Visualising Phylogenetic Trees For calculating phylogenetic trees, the alignments of the CapZ, actin and actin-related, tubulin, kinesin, and the myosin protein families were split into sub-family alignments and the respective alignments were concatenated (supplementary file S2, Schizosaccharomycetes (Sc. pombe, Sc. octosporus, Sc. japonicus, and Sc. cryophilus) were taken. Gblocks v.0.91b (Talavera & Castresana 2007) was run with standard parameters, and half gap positions allowed to reduce the concatenated alignment from 28,202 amino acids to 9,788 amino acids in 230 blocks. The phylogenetic trees were generated using four different methods: Neighbor-Joining (NJ), maximum-likelihood (ML), Bayesian inference and split networks. 1) ClustalW v.2.0.10 (Thompson et al. 2002) was used to calculate unrooted trees with the NJ method. For each dataset, bootstrapping with 1,000 replicates was performed. 2) ML analysis with estimated proportion of invariable sites and bootstrapping (1.000 replicates) was performed using RAxML v.7.3.1 (Stamatakis et al. 2008). To this end, ProtTest v.3.2 was used first to determine the most appropriate of the available 112 amino acid substitution models (Darriba et al. 2011). Within ProtTest, the tree topology was calculated with the BioNJ algorithm and both the branch lengths and the model of protein evolution were optimized simultaneously. The Akaike Information Criterion with a modification to control for small sample size (AICc, with alignment length representing sample size) identified the WAG model with gamma model of rate -9- Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 Supplementary Material online). As outgroup, sequences from heterogeneity and site-specific evolutionary rates to be the best. 3) The ML analysis was repeated on an additional dataset generated from the concatenated alignment. In this dataset, amino acids encoded by CUG codons were substituted by “X”. Gblocks was invoked as described above, with standard parameters and half gap positions allowed, to reduce the alignment to 8,903 positions in 297 blocks. The phylogenetic tree was generated with RAxML under the WAG+Γ+F model and 1,000 replicates were performed. 4) Posterior probabilities were generated using MrBayes v3.2.1. (Ronquist & Huelsenbeck 2003). Using the mixed amino-acid option, two performed. MrBayes used the WAG model for all protein alignments. The trees were calculated based on a reduced dataset generated by Gblocks using standard parameters allowing no gap positions. This dataset includes 5,038 positions in 131 blocks. Trees were sampled every 1,000th generation and the first 25% of the trees were discarded as “burn-in” before generating a consensus tree. 5) An unrooted phylogenetic split network was generated with SplitsTree v. 4.1.3.1 (Huson & Bryant 2006). The Neighbor-Net method as implemented in SplitsTree was used to identify alternative splits. Phylogenetic trees and networks were visualized with FigTree v.1.3.1 (Rambaut and Drummond) and SplitsTree v.4.13.1, respectively, and are available as supplementary fig. S1, Supplementary Material online. Estimating the divergence times The divergence times of species using AYCU were estimated with a penalizedlikelihood approach as implemented in treePL (Smith & O’Meara 2012). The underlying tree topology was generated with RAxML under the WAG+Γ+F model on the alignment with CUG codons assigned. The splits between Sa. cerevisiae and Ca. - 10 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 independent runs with 100,000 generations, 4 chains, and a random starting tree were albicans, and Sa. cerevisiae and Sc. pombe (Heckman et al. 2001; Douzery et al. 2004) were constrained both individually and combined. Results Phylogeny of Sequenced Yeast Species In order to obtain a reliable phylogeny of the sequenced yeast species we choose a small-scale phylogenomics approach. The more data (i.e. orthologs) included in the important parameter for the reliability of the computed trees is the quality of the underlying sequence data. Here, we manually determined 26 motor and cytoskeletal proteins for each of the 60 sequenced yeast species. These 26 proteins belong to six major protein families, the actin/actin-related proteins, the tubulins, the myosins, the kinesins, the dynein heavy chain proteins, and the actin capping proteins of the CapZ complex. We included all proteins of these six families in the analysis except some kinesins that are unique to certain species. Many of the actin and tubulin genes contain very short exons of one to three residues that are not included in datasets of automatic gene predictions. By manually inspecting the genomic DNA sequences we could completely reconstruct all genes as long as the genomic sequences did not contain gaps. The concatenated sequences consist of on average 19,300 residues per species amounting to 28,202 alignment positions. Gblocks (Talavera & Castresana 2007) was used with less stringent and more stringent parameters to reduce the alignment to 9,788 and 5,038 positions, respectively. NJ and ML methods were used to construct trees from the extended alignment, and the Bayesian approach to compute a tree from the shorter alignment (supplementary fig. S1, Supplementary Material - 11 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 alignment the better and more robust the phylogenetic trees should become. Another online). To evaluate the influence of misassigned CUG codons on the tree topology all CUG codons were translated with “X” and the tree reconstruction repeated. The topology of this unbiased tree is identical to that based on the alignment with assigned CUG codons. As outgroup we included four Schizosaccharomyces species. The constructed trees have almost the same topology with the major exceptions being the placing of Kluyveromyces in the NJ tree, in which they group to Lachancea and Eremothecium. To highlight alternative phylogenetic relationships between the species we also constructed a phylogenetic network using the Neighbour-Net method trees show that Yarrowia lipolytica is the earliest diverging Saccharomycetes. The Saccharomycetaceae form a highly supported group together with the early separating Pfaffomycetaceae (supplementary fig. S1, Supplementary Material online). In our analysis, Pichia kudriavzevii and Ogataea parapolymorpha branch together with Dekkera bruxellensis, Kuraishia capsulata, and Komagataella pastoris. This branch is more closely related to the Ca. albicans species group than to the Saccharomycetaceae (supplementary fig. S1, Supplementary Material online), but the species encode CUG as leucine and not as serine. The species, for which the AYCU has been determined previously (e.g., Ca. albicans, Ca. dubliniensis and Ca. parapsilosis), are part of a highly supported group in all phylogenetic trees (fig. 1 and supplementary fig. S1, Supplementary Material online). Therefore, we refer to this group from here on as CTG clade. The CTG clade consists of several “Candida” species including Ca. albicans and Ca. tenuis, and also many others such as species from Metschnikowiaceae and Debaryomycetaceae. In all trees, Pachysolen groups outside the CTG clade as sister to the Pichiaceae/Ogataea (supplementary fig. S1, Supplementary Material online). In the phylogenetic network, Pachysolen is placed at - 12 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 (fig. 1). In agreement with other studies (Kurtzman & Robnett 2013b, 2013a), the the origin of the CTG clade but there are many connections supporting alternative topologies. Interestingly, species of the Candida genus are spread all over the tree grouping to the Saccharomycetes and to various subgroups of the CTG clade (fig. 1). Conservation of Amino Acids at CUG Codon Positions in Class I Myosins In order to determine, whether CUG codons are coding for conserved amino acids (leucine or serine) we performed a thorough analysis of the alignments of all protein subfamilies. As representative example a part of the alignment of the class I myosins supplementary fig. S2, Supplementary Material online). This analysis provided several surprising findings: 1) The class I myosin alignment contains 33 positions with 100% leucine conservation but only 11 positions that show 100% serine conservation. 2) Except for Wickerhamomyces anomalus and Candida tenuis (no CUG codons in their class I myosin genes), and Cl. lusitaniae, Ca. tropicalis, and Pachysolen (two to three CUG codons in their class I myosin genes but not at conserved positions), all Saccharomycetaceae, Pfaffomycetaceae, Pichiaceae, and Ogataea species, and D. bruxellensis, K. capsulata, K. pastoris, and Y. lipolytica, have at least one CUG codon at one of the 100% conserved leucine positions. 3) Except for the sequences of the Debaryomyces and Spathaspora species investigated here, whose CUG positions are in almost all cases identical, none of the other closely related species (e.g., Saccharomyces or Pichiaceae) have conserved CUG positions (fig. 2). This is especially apparent for the six analyzed Saccharomyces species that have duplicated class I myosins but do not have even a single conserved CUG position between the two paralogs. In addition, CUG positions are not conserved within any of the orthologs. There is one position in the alignment (position 322), at - 13 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 is presented in figure 2 (the full class I myosin alignment is available in which ten of the species of the CTG clade have a CUG-encoded serine whereas serines, alanines and glycines are found at that position in the other species (fig. 2). Conservation of Amino Acids at CUG Codon Positions in General The analysis of the complete data with respect to amino acid similarity at alignment positions, at which at least one leucine/serine is present, showed that positions with leucines are stronger conserved than positions with serines although in total serines have been observed at as many alignment positions as leucines (fig. 3A and conservation of each position alone (window size = 0) and within its closest environment (window size = 3). The analysis of CUG codons present at leucine positions showed that CUG codons are not evenly distributed to all positions with leucines but enriched at highly conserved alignment positions (fig. 3A). In contrast, there are far less CUG codons at serine positions, and these are in general present at less conserved positions compared with the leucine encoding CUG codons. However, although small there are a considerable number of CUG codons at highly conserved and thus discriminative serine positions (fig. 3A). In order to determine the amino acid conservation, it is extremely important to resolve protein sub-family relationships. For example, the kinesin-1 sub-family members have a highly conserved serine at the same position where the kinesin-5 proteins have a highly conserved leucine (fig. 3B). Assignment of the Codon Usage Scheme by Amino Acid Conservation Patterns In order to determine the CUG codon translation we investigated whether the CUG codons within the analyzed sequences of each species are at highly conserved leucine - 14 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 supplementary fig. S3, Supplementary Material online). We determined the or serine positions in the alignment. Therefore, we first determined the percentage of CUG codons at extremely conserved leucine and serine positions (conservation score of ≥90%; fig. 4). The analysis resulted in three groups of species: The first group, including Sa. cerevisiae, had 2.8-26.4% (on average 14.3%) of the CUG codons at highly conserved leucine positions. The second group, including Ca. albicans, had CUG codons at conserved serine positions although the total numbers were considerably lower than those for the leucine positions (0.9 to 7.5%; on average 3.4%), Pachysolen tannophilus did not have any CUG codons at either conserved a conservation score of more than or equal to 80% (fig. 4). In the same group of species as at the conservation score of more than or equal to 90%, 8.4-44.4% (on average 24.5%) of the CUG codons are at leucine positions. At that conservation level, one of the CUG codons of Yarrowia is found at a conserved serine position. At the same conservation score (≥80%), all other species have CUG codons at highly conserved serine positions (2.2-11.8%; on average 6.1%). In five species, one to two of the CUG codons are at conserved leucine positions (Candida metapsilosis, Ca. tenuis, Cl. lusitaniae, Spathaspora arborariae and Millerozyma). Thus, the species clearly separate into two groups, one with many CUG codons almost exclusively at highly conserved leucine positions, and the other with CUG codons at serine positions (fig. 4). At a modest sequence conservation level (conservation score ≥50%), 5591 (61.1%) of the CUG codons of the first group are at leucine positions (the maximum is 70.2% of the CUG codons in Saccharomyces kudriavzevii). At this conservation level 54 CUG codons (0.6%) are at modestly conserved serine positions (conservation score ≥50%). At the same conservation level, 17% of the CUG codons of the second group (in total 332 codons) are at serine positions, whereas 44 CUG codons (2.4 %) - 15 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 position. A similar situation was found when analyzing the CUG codon distribution at are at leucine positions with the same conservation level. The same pattern can be observed when using scores calculated column wise instead of window wise (window size 0; supplementary fig. S4, Supplementary Material online). CUG Codon Positions Shared between Species It was proposed that all CUG codons had been removed during the time of codon ambiguity and later reintroduced within the Saccharomycetaceae and Candida branches (Massey et al. 2003). If this model was true, species of the positions within but not across these two branches. In order to evaluate CUG codon position conservation we determined the number of CUG codons shared between every two species (fig. 5). This CUG codon position analysis is independent from CUG codon assignment. For comparison, we analyzed the conservation of the CUG codon positions at all alignment positions (fig. 5, lower triangle) and at positions with a modest conservation level of more than or equal to 50% (fig 5, upper triangle). Independent of whether sequence conservation at alignment positions is required the yeast species group into two distinct classes. Species of the first group containing the Saccharomycetaceae, Pfaffomycetaceae, Pichiaceae, and Ogataea species, and Dekkera bruxellensis, Komagataella pastoris, Kuraishia capsulate, and Yarrowia lipolytica share considerably more CUG positions than those from the second group formed by the CTG-clade species. When all CUG positions are considered, there are a few positions shared between the two groups. However, if only CUG positions at a modest conservation level of more than or equal to 50% are evaluated, the two groups do not share any CUG positions. The few CUG positions shared at all alignment positions between the two species groups reflect CUG positions in disordered and - 16 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 Saccharomycetaceae and Candida would be expected to have shared CUG codon non-conserved protein regions (see, e.g., the alignment of the class-1 myosin tail regions in supplementary fig. S2, Supplementary Material online). Interestingly, Pachysolen, does not share any CUG positions with any other species (fig. 5). Mapping the Assigned Codon Usage onto the Species Tree In order to determine the mono- or polyphyly of the AYCU assignment, we mapped the assigned codon usage onto the reconstructed species tree (fig. 6). This mapping supports the polyphyly of the AYCU assignment because the CTG clade and Saccharomycetales. The CUG codon position conservation is also in accordance with the species tree and the codon usage assignment. The CTG clade species group together in the species tree and share CUG codon positions implying that at least some of them have been introduced in the last common ancestor of the CTG clade. Pachysolen diverged separate from the CTG clade and therefore does not share any CUG codon positions with species from the CTG clade. Timeline for the Separation of the CTG Clade In order to estimate the time span, in which the CUG codon usage might have been ambiguous, we applied published branching time estimations to our species trees. However, the divergence time estimations for the splits between Archiascomycetes and other Ascomycota and within the Saccharomycetes differ by more than 700 Myr although extensive molecular and fossil record data have been used in respective studies ((Heckman et al. 2001; Hedges et al. 2004; Douzery et al. 2004; Berbee & Taylor 2007), for more references see (Hedges et al. 2006)). Therefore, we choose some of the most extreme estimates for the calibration (table 1). The divergence time - 17 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 Pachysolen branched at different points during the evolution of the of the CTG clade was calibrated on the basis of the topology and branch lengths of the ML tree under penalized-likelihood constraints as implemented in treePL (Smith & O’Meara 2012). The splits of Sc. pombe and Sa. cerevisiae, and of Sa. cerevisiae and Ca. albicans were set to 1,144 and 841 Ma (Heckman et al. 2001), and 587 and 235 Ma (Douzery et al. 2004), respectively. Under these two sets of constraints, the separation of the CTG clade was estimated to 440 and 135 Ma, respectively (table 1), implying that the reassignment of the CUG codon in the CTG clade happened at most 400 and 100 Myr, respectively, after the split of the Saccharomycetaceae and the obtained for the emergence of Pachysolen indicating that ancient yeast species emerging in the time of 220 to 190 Ma (740 to 640 Ma, respectively) could have independently reassigned the CUG codon translation. This time span is the minimum time of CUG codon ambiguity. The maximum time span depends on the last split before the first appearance of the AYCU (split of the Saccharomycetaceae and the branch containing the CTG clade, Pichiaceae, and Ogataea; fig. 6) and both the first split after the last assignment of the AYCU (split of the branch containing Pichiaceae and Ogataea from Pachysolen) and the divergence of the CTG clade (135 and 440 Ma, respectively). The maximum time span of codon ambiguity is thus estimated to be 100 Myr (400 Myr, respectively). Discussion The AYCU has already been assigned to ten of the sequenced yeast genomes, Ca. albicans (Jones et al. 2004; Butler et al. 2009), Ca. dubliniensis (Jackson et al. 2009), Ca. parapsilosis (Butler et al. 2009), Ca. tropicalis (Butler et al. 2009), Cl. lusitaniae - 18 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 CTG clade (fig. 6). Similarly, a date of at least 190 Ma (640 Ma, respectively) was (Butler et al. 2009), Debaryomyces hansenii (Dujon et al. 2004; Butler et al. 2009), Lodderomyces elongisporus (Butler et al. 2009), Mey. guilliermondii (Butler et al. 2009), Scheffersomyces stipites (Jeffries et al. 2007), and Spathaspora passalidarum (Wohlbach et al. 2011) and the standard codon usage was shown for several of the Saccharomyces species. Mass spectrometry and Edman sequencing of SMKT demonstrated that Millerozyma farinosa also uses the AYCU (Suzuki et al. 2002). An in vitro translation study showed 78 Candida species to use the AYCU including the above listed Candida strains (Sugita & Nakase 1999), and in addition Ca. recently (Wohlbach et al. 2011). Nevertheless, the gene prediction datasets for Ca. tanzawaensis and Ca. tenuis available at the respective sequencing centers include predictions both with AYCU and with standard codon usage. However, for highly resolved phylogenetic studies, protein ortholog assignments and functional studies, naming only a few applications, the accuracy of the protein sequences is important. The best and most accurate way to determine a species codon usage would be the experimental analysis of some proteins known to contain CUG codons by, for example, mass spectrometry. However, as long as proteins cannot be enriched this is extremely difficult to conduct, and the most prevalent endogenous proteins might not contain CUG codons. For example, our analysis showed that actin genes from only a few yeasts (11 out of 60 analyzed) contain CUG codons. Other approaches only provide indirect evidence. As discussed, the presence or absence of cellular components such as galactose or Q9 does not correlate with CUG codon usage. Here, we hypothesized, that the codon usage can unambiguously be assigned by analyzing three types of data: A highly resolved phylogenetic tree, sequence conservation at CUG coding positions, and conservation of CUG codon positions across species. To - 19 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 tanzawaensis and Ca. tenuis, for which the genome sequence had been assembled obtain significant numbers of sequences we were restricted to analyse sequenced genomes. Up to May 2014, genomes of 60 yeast species have been sequenced, assembled and released. We did not include genomes in the analysis, which were sequenced with low-coverage, and only included one strain per species (e.g., 73 different Sa. cerevisiae strains have been sequenced, of which only the reference strain was included in our analysis). To avoid bias by misassigning CUG codons we excluded all CUG positions from the alignment and the sequence conservation calculations. proteins, which are key-components of eukaryotic cells. We decided to only use proteins in the analysis that are at least present in all unikonts (the eukaryotic superkingdom including Fungi, Metazoa, Apusozoa and Amoebozoa) to ensure that we do not have extensive missing data. The topology of the species tree was almost identical independent of the reconstruction method used (NJ, ML and Bayesian; fig. 6). Although this particular set of species has never been analyzed before, the branching of most sub-branches and individual species is comparable to those found in previously published yeast trees (see references in Introduction section). The species known to use either the standard code (Saccharomycetaceae) or the AYCU (see above) grouped into two distinct clades. It seemed obvious to assume, but had to be shown, that all other species within these two groups, for which the codon usage has not yet been confirmed, use the same codon usage as the known species. In terms of a mono- or polyphyletic origin of the AYCU it would be necessary to resolve the codon usage of the species of the other major branches (Phaffomycetaceae clade and the - 20 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 For the reconstruction of the yeast tree we concatenated 26 cytoskeletal and motor branch containing the Pichiaceae and Ogataea) and the species that diverged before the combined Saccharomycetaceae/CTG clades such as Y. lipolytica. To determine the most probable CUG codon usage for each species we analyzed the amino acid conservation at the respective CUG codon positions in the protein sequence alignments. Amino acid positions highly conserved over hundreds of million years (divergence time of Schizosaccharomyces and Saccharomycetes) are expected to represent structurally and functionally important sites in proteins. A small residues in the other species, seems very unlikely. Similarly, large hydrophobic residues at positions conserved in small and polar amino acids should be strongly disfavoured. Our data show that the assignment of the standard codon usage is facilitated by the enrichment of CUGs at highly conserved leucine positions. In this respect only determining the class I myosin of a species with unknown CUG codon usage and comparing this sequence to the available data (fig. 2 and supplementary fig. S2, Supplementary Material online) would in most cases be enough to find out whether the respective species uses the standard or the alternative yeast CUG codon translation. If CUG codons can be mapped to conserved leucine positions, it will be highly probable that the respective species translates CUG as leucine. The AYCU will be assigned to a species if CUG codons are mapped to highly conserved serine positions or if CUG codons are not present at highly conserved leucine positions (proof by contradiction). In addition to the class I myosins, several proteins could be identified and compared until a conserved serine position was found. The latter approach would lead to an unambiguous assignment and was used in our analysis here. It has been proposed that CUG codons were not reintroduced randomly in CTG - 21 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 polar residue at a position, which is characterized by conserved large hydrophobic clade genomes but by avoiding structurally sensitive sites (Rocha et al. 2011). This is in agreement with the different biochemical properties of leucine and serine and the associated consequences for protein folding and solubility. Leucine positions are in general stronger conserved than serine positions, leucines are usually located in the core of protein domains, and mutations at leucine positions are therefore restricted. Compared with serine, the introduction of leucine encoding CUG codons is possible by synonymous mutations from CUN codons or by transition from thymine from the UUG codon. Transitions are more favourable and generated at higher frequencies than conservation of the CUG codons within the yeast species using the Standard Code. In contrast, mutating an existing serine codon to a CUG codon would require at least two mutations. Mutations from other codons than leucine or serine codons would either require the less likely transversions (e.g., GUG to CUG) or very unusual amino acid mutations (e.g., from proline to serine). Therefore, it is not surprising to find less CUG codons in CTG clade species, and those codons at less conserved positions. Overall, we could demonstrate, that especially the CUG codons coding for leucine, and also a substantial number of the CUG codons in species using the AYCU, are at conserved alignment positions (fig. 4). Although seven of the CUG codons coding for leucine and serine, respectively, are at positions supporting the respective opposite codon translation based on an amino acid conservation level of more than or equal to 80%, the majority of the CUG codon positions supports either the standard or the AYCU. This is in accordance with protein structure and in vivo studies showing that single point mutations and even partial reversion of the CUG identity are tolerated in Ca. albicans although large numbers of mutations cause profound cellular changes such as morphological variation and altered gene expression (Cutfield et al. 2000; - 22 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 transversions. Such silent substitutions seem to be very common given the low Gomes et al. 2007; Miranda et al. 2007; Rocha et al. 2011; Miranda et al. 2013; Bezerra et al. 2013). Because we excluded all CUG codon positions from the analysis, our assignments are independent from previous knowledge. Our data unambiguously supports the assignment of the standard codon usage to the Saccharomycetaceae, and the AYCU to the eight species, which were known to translate CUG as serine. From these data we conclude, that the CUG codon assignment can unambiguously be resolved by analyzing a considerable amount of sequences and that our CUG codon translation assignment is robust. Therefore, we assigned the AYCU to the CTG clade assignment of the CUG codon usage, we acknowledge that the final judgement should come from experiments. When plotting the assigned AYCU onto the species in the phylogenetic trees it became evident that the CUG codon usage must have been ambiguous for a long time in the ancestry of the Saccharomycetes leading to polyphyly of the AYCU (fig. 6). This means that ancestral species that emerged during this time span could have adopted either the standard or the AYC usage. In the case of a long time span of codon ambiguity, the original CUG codons in the ancestor of the Saccharomycetes were most likely erased and later reintroduced in the emerging branches either as leucine codons or as serine codons. In the less likely and mutually exclusive scenario, the CUG ambiguity appeared only once in the ancestor of the CTG clade. The combination of the reconstructed species trees with our CUG codon assignment strongly supports the first scenario of a polyphyletic origin of the AYCU. - 23 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 species and Pachysolen. Although our data strongly support an unambiguous Based on the phylogeny of the analyzed species Pachysolen does not group to the CTG clade. In addition, Pachysolen does not have any CUG codons at highly conserved leucine positions, but at conserved serine positions. The phylogenetic grouping of Pachysolen has not been analyzed in detail so far, and not at all in the context of a taxonomic sampling as used in our study. The evolutionary grouping and CUG codon usage were also not analyzed in the study of the Pachysolen draft genome sequence (Liu et al. 2012). A previous study based on 28s rDNA and 18s rDNA only resolved the Saccharomycetes but not the grouping of Pachysolen and none of the addition, support for the branchings was low (Suh et al. 2006). In contrast, we used a phylogenomics approach for the tree reconstruction resulting in strongly supported branching of Pachysolen. The separate branching of Pachysolen is in agreement with the previously proposed model of a long time span of codon usage ambiguity (Sugita & Nakase 1999). If the CUG codon usage was ambiguous for a long time, it would be natural to assume that species diverged within this time span, and that those species could have adapted either of the usages. Supposed, that the codon usage was ambiguous for a long time, it seems highly non-parsimonious, that either species did not branch within this time span or that all species except the ancestor of the CTG clade adapted the standard codon usage. A monophyletic group of the species using the AYCU is also highly unlikely given the missing conservation of CUG codon positions between the CTG clade species and Pachysolen. To date the codon reassignments we dated the ML tree based on various separation time estimates for the splits of Sc. pombe and Sa. cerevisiae, and Sa. cerevisiae and Ca. albicans (fig. 6 and table 1). Accordingly, there - 24 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 other clades shown to be monophyletic in our trees (including the CTG clade). In is a maximum time span of codon ambiguity of 100-400 Myr ranging from the first possible appearance of the serine tRNASerCAG (separation of the Saccharomycetes and the branch containing the CTG clade, the Pichiacea and Ogataea) to the split of the most basal species of the CTG clade (Ca. tenuis). Assuming an unambiguous usage of the CUG codon in Pachysolen and the CTG clade, the minimal time span of CUG codon ambiguity was 30 to 100 million years (fig. 6 and table 1). This is well in agreement with a previous study, which suggested an ambiguous usage of the CUG codon for 100 million years (Sugita & Nakase 1999). However, the presented time assume that further sequenced species will change the time estimates in two directions (this might already happen by analyzing the yeast species whose data is still under embargo). The inclusion of more basal CTG clade species will decrease the current codon ambiguity time. The identification of further species using the AYCU, which do not group to the CTG clade or to Pachysolen, instead might increase both the codon ambiguity and codon reassignment time ranges. Such species include yeasts, which might even branch before the split of the Saccharomycetaceae and the CTG clade/Pichiacea/Ogataea, as well as additional species diverging in the sister branch of the CTG clade, and species separating early in the Saccharomycetaceae branch. Acknowledgments We would like to thank Prof. Christian Griesinger for his continuous generous support, and Prof. Toni Gabaldón and Dr. Qinhong Wang for their permission to use the Ca.metapsilosis CANME and Ca.maltosa Xu316 genome assembly data prior to publication. This project has been funded by the Deutsche Forschungsgemeinschaft - 25 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 estimates are based on the species analyzed in this study. It is well reasonable to DFG Grant KO 2251/13-1 to M.K. and a Synaptic Systems fellowship to S.M., and was partly supported by the Göttingen Graduate School of Neurosciences and Molecular Biosciences (DFG Grants GSC 226/1 and GSC 226/2). Authors contribution SM and MK assembled and annotated genes, analyzed the data and wrote the manuscript. Berbee ML, Taylor JW. 2007. Rhynie chert: a window into a lost world of complex plant-fungus interactions. New Phytol. 174:475–479. doi: 10.1111/j.14698137.2007.02080.x. Bezerra AR et al. 2013. Reversion of a fungal genetic code alteration links proteome instability with genomic and phenotypic diversification. Proc. Natl. Acad. Sci. U.S.A. 110:11079–11084. doi: 10.1073/pnas.1302094110. Butler G et al. 2009. Evolution of pathogenicity and sexual reproduction in eight Candida genomes. Nature. 459:657–662. doi: 10.1038/nature08064. Capra JA, Singh M. 2007. Predicting functionally important residues from sequence conservation. Bioinformatics. 23:1875–1882. doi: 10.1093/bioinformatics/btm270. Crooks GE, Hon G, Chandonia JM, Brenner SE. 2004. WebLogo: a sequence logo generator. Genome Res. 14:1188 – 1190. Cutfield JF, Sullivan PA, Cutfield SM. 2000. Minor structural consequences of alternative CUG codon usage (Ser for Leu) in Candida albicans exoglucanase. Protein Eng. 13:735–738. Darriba D, Taboada GL, Doallo R, Posada D. 2011. ProtTest 3: fast selection of bestfit models of protein evolution. Bioinformatics. 27:1164–1165. doi: 10.1093/bioinformatics/btr088. Diezmann S, Cox CJ, Schönian G, Vilgalys RJ, Mitchell TG. 2004. Phylogeny and evolution of medical species of Candida and related taxa: a multigenic analysis. J. Clin. Microbiol. 42:5624–5635. doi: 10.1128/JCM.42.12.5624-5635.2004. - 26 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 References Douzery EJP, Snell EA, Bapteste E, Delsuc F, Philippe H. 2004. The timing of eukaryotic evolution: Does a relaxed molecular clock reconcile proteins and fossils? Proc Natl Acad Sci U S A. 101:15386–15391. doi: 10.1073/pnas.0403984101. Dujon B et al. 2004. Genome evolution in yeasts. Nature. 430:35–44. doi: 10.1038/nature02579. Fitzpatrick DA, Logue ME, Stajich JE, Butler G. 2006. A fungal phylogeny based on 42 complete genomes derived from supertree and combined gene analysis. BMC Evol. Biol. 6:99. doi: 10.1186/1471-2148-6-99. Gomes AC et al. 2007. A genetic code alteration generates a proteome of high diversity in the human pathogen Candida albicans. Genome Biol. 8:R206. doi: 10.1186/gb-2007-8-10-r206. Hatje K et al. 2011. Cross-species protein sequence and gene structure prediction with fine-tuned Webscipio 2.0 and Scipio. BMC Res. Notes. 4:265. doi: 10.1186/17560500-4-265. Heckman DS et al. 2001. Molecular Evidence for the Early Colonization of Land by Fungi and Plants. Science. 293:1129–1133. doi: 10.1126/science.1061457. Hedges SB, Blair JE, Venturi ML, Shoe JL. 2004. A molecular timescale of eukaryote evolution and the rise of complex multicellular life. BMC Evolutionary Biology. 4:2. doi: 10.1186/1471-2148-4-2. Hedges SB, Dudley J, Kumar S. 2006. TimeTree: a public knowledge-base of divergence times among organisms. Bioinformatics. 22:2971–2972. doi: 10.1093/bioinformatics/btl505. Huson DH, Bryant D. 2006. Application of phylogenetic networks in evolutionary studies. Mol. Biol. Evol. 23:254–267. doi: 10.1093/molbev/msj030. Jackson AP et al. 2009. Comparative genomics of the fungal pathogens Candida dubliniensis and Candida albicans. Genome Res. 19:2231–2244. doi: 10.1101/gr.097501.109. Jeffries TW et al. 2007. Genome sequence of the lignocellulose-bioconverting and xylose-fermenting yeast Pichia stipitis. Nat. Biotechnol. 25:319–326. doi: 10.1038/nbt1290. Jones T et al. 2004. The diploid genome sequence of Candida albicans. Proc. Natl. Acad. Sci. U.S.A. 101:7329–7334. doi: 10.1073/pnas.0401648101. Kurtzman C, Fell JW, Boekhout T. 2011. The Yeasts: A Taxonomic Study. 5th edition. London, Burlington, San Diego: Elsevier. - 27 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 Hammesfahr B, Kollmar M. 2012. Evolution of the eukaryotic dynactin complex, the activator of cytoplasmic dynein. BMC Evol. Biol. 12:95. doi: 10.1186/1471-2148-1295. Kurtzman CP. 2011. Phylogeny of the ascomycetous yeasts and the renaming of Pichia anomala to Wickerhamomyces anomalus. Antonie Van Leeuwenhoek. 99:13– 23. doi: 10.1007/s10482-010-9505-6. Kurtzman CP, Robnett CJ. 2013a. Alloascoidea hylecoeti gen. nov., comb. nov., Alloascoidea africana comb. nov., Ascoidea tarda sp. nov., and Nadsonia starkeyihenricii comb. nov., new members of the Saccharomycotina (Ascomycota). FEMS Yeast Res. 13:423–432. doi: 10.1111/1567-1364.12044. Kurtzman CP, Robnett CJ. 1997. Identification of clinically important ascomycetous yeasts based on nucleotide divergence in the 5’ end of the large-subunit (26S) ribosomal DNA gene. J. Clin. Microbiol. 35:1216–1223. Kurtzman CP, Robnett CJ. 2013b. Relationships among genera of the Saccharomycotina (Ascomycota) from multigene phylogenetic analysis of type species. FEMS Yeast Res. 13:23–33. doi: 10.1111/1567-1364.12006. Liu X, Kaas RS, Jensen PR, Workman M. 2012. Draft genome sequence of the yeast Pachysolen tannophilus CBS 4044/NRRL Y-2460. Eukaryotic Cell. 11:827. doi: 10.1128/EC.00114-12. Massey SE et al. 2003. Comparative evolutionary genomics unveils the molecular mechanism of reassignment of the CTG codon in Candida spp. Genome Res. 13:544– 557. doi: 10.1101/gr.811003. Miranda I et al. 2007. A genetic code alteration is a phenotype diversity generator in the human pathogen Candida albicans. PLoS ONE. 2:e996. doi: 10.1371/journal.pone.0000996. Miranda I et al. 2013. Candida albicans CUG mistranslation is a mechanism to create cell surface variation. MBio. 4: pii e00285-13. doi: 10.1128/mBio.00285-13. Odronitz F, Kollmar M. 2007. Drawing the tree of eukaryotic life based on the analysis of 2,269 manually annotated myosins from 328 species. Genome Biol. 8:R196. doi: 10.1186/gb-2007-8-9-r196. Odronitz F, Kollmar M. 2006. Pfarao: a web application for protein family analysis customized for cytoskeletal and motor proteins (CyMoBase). BMC Genomics. 7:300. doi: 10.1186/1471-2164-7-300. Odronitz F, Pillmann H, Keller O, Waack S, Kollmar M. 2008. WebScipio: an online tool for the determination of gene structures using protein sequences. BMC Genomics. 9:422. doi: 10.1186/1471-2164-9-422. Ohama T et al. 1993. Non-universal decoding of the leucine codon CUG in several Candida species. Nucleic Acids Res. 21:4039–4045. - 28 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 Kurtzman CP, Suzuki M. 2010. Phylogenetic analysis of ascomycete yeasts that form coenzyme Q-9 and the proposal of the new genera Babjeviella, Meyerozyma, Millerozyma, Priceomyces, and Scheffersomyces. Mycoscience. 51:2–14. doi: 10.1007/s10267-009-0011-5. Rambaut A, Drummond A. 2009. FigTree v1.3.1. http://tree.bio.ed.ac.uk/software/figtree/, visited May 9, 2014. Robert V et al. 2013. MycoBank gearing up for new horizons. IMA Fungus. 4:371– 379. doi: 10.5598/imafungus.2013.04.02.16. Rocha R, Pereira PJB, Santos MAS, Macedo-Ribeiro S. 2011. Unveiling the structural basis for translational ambiguity tolerance in a human fungal pathogen. PNAS. 108:14091–14096. doi: 10.1073/pnas.1102835108. Ronquist F, Huelsenbeck JP. 2003. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 19:1572–1574. Santos MAS, Gomes AC, Santos MC, Carreto LC, Moura GR. 2011. The genetic code of the fungal CTG clade. C. R. Biol. 334:607–611. doi: 10.1016/j.crvi.2011.05.008. Smith SA, O’Meara BC. 2012. treePL: divergence time estimation using penalized likelihood for large phylogenies. Bioinformatics. 28:2689–2690. doi: 10.1093/bioinformatics/bts492. Stamatakis A, Hoover P, Rougemont J. 2008. A rapid bootstrap algorithm for the RAxML Web servers. Syst. Biol. 57:758–771. doi: 10.1080/10635150802429642. Stanke M, Morgenstern B. 2005. AUGUSTUS: a web server for gene prediction in eukaryotes that allows user-defined constraints. Nucleic Acids Res. 33:W465–467. doi: 10.1093/nar/gki458. Sugita T, Nakase T. 1999. Non-universal usage of the leucine CUG codon and the molecular phylogeny of the genus Candida. Syst. Appl. Microbiol. 22:79–86. doi: 10.1016/S0723-2020(99)80030-7. Suh S-O, Blackwell M, Kurtzman CP, Lachance M-A. 2006. Phylogenetics of Saccharomycetales, the ascomycete yeasts. Mycologia. 98:1006–1017. Suzuki C, Kashiwagi T, Hirayama K. 2002. Alternative CUG codon usage (Ser for Leu) in Pichia farinosa and the effect of a mutated killer gene in Saccharomyces cerevisiae. Protein Eng. 15:251–255. Szkopińska A. 2000. Ubiquinone. Biosynthesis of quinone ring and its isoprenoid side chain. Intracellular localization. Acta Biochim. Pol. 47:469–480. Talavera G, Castresana J. 2007. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst Biol. 56:564–577. doi: 10.1080/10635150701472164. Thompson JD, Gibson TJ, Higgins DG. 2002. Multiple sequence alignment using ClustalW and ClustalX. Curr Protoc Bioinformatics. Chapter 2:Unit 2.3. doi: 10.1002/0471250953.bi0203s00. - 29 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 Schauer F, Hanschke R. 1999. [Taxonomy and ecology of the genus Candida]. Mycoses. 42 Suppl 1:12–21. Tourancheau AB, Tsao N, Klobutcher LA, Pearlman RE, Adoutte A. 1995. Genetic code deviations in the ciliates: evidence for multiple and independent events. EMBO J. 14:3262–3267. Turunen M, Olsson J, Dallner G. 2004. Metabolism and function of coenzyme Q. Biochim. Biophys. Acta. 1660:171–199. Watanabe K, Yokobori S-I. 2011. tRNA Modification and Genetic Code Variations in Animal Mitochondria. J Nucleic Acids. 2011:623095. doi: 10.4061/2011/623095. Wohlbach DJ et al. 2011. Comparative genomics of xylose-fermenting fungi for enhanced biofuel production. Proc. Natl. Acad. Sci. U.S.A. 108:13212–13217. doi: 10.1073/pnas.1103039108. Figure Legends Figure 1 Phylogenetic relationship between Saccharomycetes The unrooted phylogenetic network was generated using the Neighbor-Net method as implemented in SplitsTree 4.1.3.1. The Schizosaccharomycetes species were included as outgroup. The network strongly supports the Saccharomycetaceae and the CTG clade (highlighted in orange). The grouping of Pachysolen tannophilus is not unambiguously resolved. Species of the genus "Candida" are highlighted in blue (teleomorph names) and green (if anamorphs of the species are called "Candida") showing the paraphyly (or mis-assingment) of this genus. Orange and purple dots mark species, for which alternative or standard codon usage has already been shown elsewhere. - 30 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 Yokobori S, Suzuki T, Watanabe K. 2001. Genetic code variations in mitochondria: tRNA as a major determinant of genetic code plasticity. J. Mol. Evol. 53:314–326. doi: 10.1007/s002390010221. Figure 2 Sequence alignment of the yeast class I myosins highlighting leucines and serines encoded by CUG The protein sequence alignment represents part of the class I myosin alignment (for the complete alignment see fig. S2, Supplementary Material online). Numbers on the left denote the residue numbers of the first amino acids of the sequences in this section of the alignment. All CUG positions occurring in the aligned class I myosin genes are highlighted. We assigned the most probable translation scheme to each species and translated the CUG codons accordingly. Blue and green boxes indicate Figure 3 Conservation of serine and leucine positions A) The charts show the amino acid conservation at all alignment positions of the Gblocks reduced concatenated alignment of the 26 cytoskeletal and motor proteins, at which at least one leucine (upper charts) or one serine (lower charts) is present. The sequence conservation has been determined based on the Property-Entropy divergence, as described in (Capra & Singh 2007). With a window size of 0, each column is scored independently (left row), while the surrounding 3 columns are also taken into account with a window size of 3 (right row). Blue bars represent the number of alignment positions with a conservation score for leucine and serine residues, respectively, within the given half-bounded intervals. Red bars denote the number of alignment positions with respect to conservation, at which at least one CUG codon is present independent of its translation. Green bars give the total numbers of CUG codons at the respective alignment positions. B) The weblogos (Crooks et al. 2004) show the sequence conservation of two kinesin sub-families, kinesin-1 and kinesin-5, within the family-defining motor domain around the highly - 31 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 CUG codons coding for leucine and serine, respectively. conserved switch II and α-helix α4 motifs. At the position within α4 marked by a grey bar kinesin-1 sequences contain a highly conserved serine while kinesin-5 sequences contain a highly conserved leucine indicating the need to resolve subfamily relationships when determining CUG codon usage by sequence conservation. Figure 4 CUG codon usage at conserved leucine and serine positions The graph presents the CUG codon usage with respect to alignment position alignment positions with conservation scores of ≥90%, ≥80%, and ≥50% and a window size of 3, respectively, in the concatenated alignment of the 26 motor and cytoskeletal proteins. Red and blue colours denote the percentages of CUG codons present at alignment positions enriched in serines and leucines, respectively. For comparison, we plotted the percentages of CUG codons at positions of the assigned codon translation to the left (% CUG codons at leucine positions for species using the standard code and % CUG codons at serine positions for species using the AYCU) and the percentages of CUG codons present at alignment positions enriched in the respective other amino acid to the right. When considering only highly conserved alignment positions (≥90% conserved) the CUG codon translation assignment is unambiguous. Species using the AYCU are highlighted in bold. Figure 5 CUG codon positions conserved between species The heatmap represents the number of CUG codon positions shared between every two species. The upper triangle shows the number of shared CUG codon positions at those positions in the concatenated alignment of the full-length sequences, which have - 32 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 conservation. For each species we determined the percentage of CUG codons at a conservation score of at least 50%, the lower triangle the number of shared CUG codons at all alignment positions. The diagonal represents the total number of CUG codons in the respective species. The number of CUG codons is coloured on a logarithmic scale. Species encoding CUG as serine are typed in bold. Figure 6 Timeline of the CUG codon reassignment The tree presents the maximum-likelihood topology generated under the Γ + WAGF model in RAxML showing branch lengths for the concatenated alignments of 26 bootstrap replicates), MrBayes (posterior probabilities) and ClustalW trees (1,000 replicates) is indicated (for more details see fig. S1, Supplementary Material online). Species using the AYCU are highlighted in bold. With the splits of Sc. pombe and Sa. cerevisiae and Sa. cerevisiae and Ca. albicans set to 587 mya and 235 mya, respectively, the divergence of the CTG clade was estimated to 190 mya using treePL (Smith & O’Meara 2012). The scale bar denotes amino acid substitutions per site. The width and colour of the branches to extant species represent the total number of CUG codons in the respective concatenated sequences. Tables Table 1 Estimation of the time span of the CUG codon reassignment Sc. pombe-Sa. cerevisiae Sa. cerevisiae-Ca. albicans CUG codon ambiguity CUG codon reassignment ( 1,850 841a 520-841 740-855) 1,144a 730 415-730 605-695 - 33 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 cytoskeletal and motor proteins. Support for major branches of the RAxML (1,000 1,144a 841a 440-841 640-740 485 235b 135-235 195-225 587b 360 205-360 295-340 587b 235b 135-235 190-220 Estimated times in million years for the split between Sa. cerevisiae and Ca. albicans, Sa. cerevisiae and Sc. pombe and the reassignment of the CUG codon were calculated with a Penalized-likelihood program for the phylogenetic tree generated with RAxML a Time constrains for the splits between Sa. cerevisiae and Ca. albicans, and Sa. cerevisiae and Sc. pombe as reported by (Heckman et al. 2001) and b(Douzery et al. 2004), respectively, were included both individually and combined. For each CUG codon reassignment time computation, constrained time estimates are typed in bold. The CUG codon reassignment time estimates are given as differences between the time estimates for the divergence of P. tannophilus from the Pichiaceae and Ogataea species (first number), and the first occurrence of the tRNASerCAG (second number; separation of the CTG clade from the branch containing P. tannophilus, the Pichiaceae and Ogataea species). The CUG codon ambiguity time estimates, however, are given as differences between the time estimates for the latest possible fixation of the tRNASerCAG (first number; separation of Ca. tenuis, the most basal species in the CTG clade) and the earliest possible invention of the tRNASerCAG (second number; separation of the Saccharomycetaceae and the CTG clade). The time estimate for the divergence of Sc. pombe and Sa. cerevisiae (1,850 mya) based on the constraint of 841 mya for the split of Sa. cerevisiae and Ca. albicans predates the emergence of eukaryotes and should be ignored. - 34 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 using the WAG+Γ+F model. Supporting Information Figure S1 Phylogenetic trees This file contains the phylogenetic trees generated with ClustalW, RAxML and MrBayes based on the concatenated alignments of 26 motor and cytoskeletal proteins. Bootstrap support values and posterior probabilities are reported in absolute values dated by constraining the time for the splits between Sa. cerevisiae and Ca. albicans, Sa. cerevisiae and Sc. pombe is included. The scale bar denotes million years. Moreover, a RAxML tree was generated on a modified alignment, in which CTG encoded amino acids were replaced by “X”. In all trees, major branches as well as the “CTG-clade” are highlighted with colours. Figure S2 Alignment of the yeast class I myosins with all CUG positions highlighted. All CUG positions occurring in the aligned class I myosin genes are highlighted. Blue and green boxes indicate CUG codons coding for leucine and serine, respectively. - 35 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 (ClustalW) and relative values (RAxML and MrBayes). In addition, the RAxML tree Figure S3 Conservation of serine and leucine positions Opposed to fig. 3, alignment positions were scored on all sequences of our dataset, including inferred CUG translations. Scores were calculated with window sizes of 0 (charts on the left) and 3 (charts on the right). Figure S4 CUG codon usage at conserved leucine and serine positions Alignment positions were scored without sequences containing CUG codons. Opposed to fig. 4, the window size for score computations was set to 0. This table lists the teleomorph name (used as reference) of every analyzed species, names of anamorphs, alternatively used scientific names and the species abbreviations as used throughout this analysis. For every species the number of CapZ, actin and actin-related, tubulin, kinesin, myosin, and dynein heavy chain genes is given. Genome analyses, if published, accession numbers and sequencing centres are referenced. File S2 Concatenated alignment This zip-compressed file contains the concatenated multiple sequence alignment of the 26 cytoskeletal and motor proteins. File S3 Conservation of serine and leucine positions This file contains the raw data to fig. 3 and fig. S3. - 36 - Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 File S1 List of analyzed species File S4 CUG codon usage at conserved leucine and serine positions This file contains the raw data to fig. 4 and fig. S4. Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 - 37 - ru ab ck erh Wi om rha rap isisp ora Tetra nd am cke Tet rli pisis pora Vande ne yc mo es my ces blat tae ra in an phaffi i CBS i far S B C 67 ii C 7 BS C 09 88 sN rrii R NR RL RL Y-3 Y-1 66 031 628 0.1 4 4417 Yarrowia lipolytica CLIB122 olyspo ra DSM delbrueckii C ii n e s an en r. h s a n m v ii zy s ha n e o s 4 r e 98 4 li le myc s hanC 23 1 YL M ryo yce TC R NR a om M s i b 20 u 7 y e 2 n r 4 e a D eb TCC da t A i e D d a Can itani s u l ra 77 2 o a p l o s i c i fruct Clav a i w o k ni Metsch 6260 C C T A ii d n o guillierm Meyerozyma a BR alu zyma p Torulaspora ii N om cife CBS rwalto jad a s no 4 6 70 BS 1146 70294 Schizosaccharomyces osporus Schizoct osacch cryopharomyces ilus Schi zosa ccha p omb romyc Sch e es izos japaocchar nicuomy s ces Zygosaccharomyces bailii ISA1307 ii CBS732 x u o r s e c y m o Zygosacchar ns is N ia R CB nia a R f L ric nag Y- S 4 a Nak a n 12 21 nis aC ase 63 h om ii C BS yce BS 25 0 s ba 8 17 797 cilli spo rus Cand CBS ida g 772 labra 0 ta CB Nakase S138 omyce s delph ens Candida niva riensis CBS 9 is CBS 2170 983 Candida bracarensis CBS 10154 ne re te llii ai hs za c Kaz ach sta Ka um Na tan ym ac as ad ym oz ov oz um ov Na nia spo Klu ra v Klu yv ine ae yv erom T02 er om yce /19 AF yc s m arx es ian lac us tis va NR r. m RL arx Yian 11 us 40 K 10 -2 54 CD 2 ii U 86 m 18 C ha TC er ii A ick ar w 95 stu es 08 15 ae yc C1 #72 m es TC PG ro yc ii A BV ve om yp eD s s uy er aria go Kl yv bal m Klu ciu cym 651 the ium -12 mo hec RL Y t R N i Ere mo 6340 ver CBS Ere kluy rans cea otole han therm a 2644 Lac ance i NCYC Lach ancea walti Lach ellii CBS 4332 cc ha rom yce nse Candida cast Sa Ha Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 H-6 a l o c i r arbo -17217c s e c y rom RL Y e S288815 02 R N Saccha oxus evisia IFO 1 O 18 3-6C d a r er tae i IF 62 c es pa s c y e a evi us c m k y o i r a h m avz an rom s a Sacc e h c yc dri bay Sac u m ro s k ces a e h c c y y m Sac m o ar aro h h cc cc a Sa S CT C1 ccha romyces Pha ffo ia dr my ce tac ea e ku Wi be e Sc hiz osa ia ch Pi Cy Candida orthopsilosis Co 90-125 Candida met apsilosis CAN ME Candid Cand a albican s WO-1 Can ida dub linie Ca dida tr nsis nd CD3 o pica i 6 d a l S i s p m Sp at alto MYAh a 340 Sc t a sa 4 Xu he has spor 3 p ffe 16 or a ar a p bo rs om rar as s yc ali iae U es da FM ru st Gm ip 19 i ti NR .1A sC RL BS Y27 60 90 54 7 er kk De Pic hia ce ae psilosis Candida para us gispor es el4o2n39 romyc Lodde NRRL YB- 60 Y-24 RRL ilus N noph 93 n ta n S 19 ysole a CB Pach 5 ulat ap s S11 -1 ia c sG ori DL aish Kur ast ha la p orp ael lym gat o 9 ma ap 49 r Ko pa S2 ea CB is ata 12 ns Og ii M lle ev xe vz O a tae a g CT Gc la d 75 55 Reported standard codon usage tac eae Reported alternative codon usage CnbMyo1A CdnMyo1A CacaMyo1A NabMyo1A NdMyo1A CnbMyo1B CglMyo1A Sab_aMyo1B Sc_cMyo1B SaaMyo1B SakMyo1B SmiMyo1B Sap_aMyo1B NacMyo1 KaaMyo1 VpMyo1 TtpMyo1A TtpMyo1B TtbMyo1A TtbMyo1B ErgMyo1 ErcMyo1 HsvMyo1 KlMyo1 KmmMyo1 KlwMyo1 KaMyo1 KnMyo1 LwMyo1 LatMyo1 Lak_aMyo1 Sab_aMyo1A Sc_cMyo1A SaaMyo1A SakMyo1A SmiMyo1A Sap_aMyo1A ZrMyo1 TodMyo1 CglMyo1B WaMyo1Alpha WicMyo1 CyjMyo1 YlMyo1 CdnMyo1B Kop_bMyo1 NdMyo1B KcMyo1 OgpMyo1 DebMyo1 PiuMyo1 NabMyo1B ZbMyo1 Ca_bMyo1 StaMyo1 ShpMyo1 CllMyo1 MefMyo1 CnmMyo1 DhhMyo1 DehMyo1 CatMyo1 LoeMyo1 MrgMyo1 ShsMyo1 CameMyo1 MifMyo1Alpha CadMyo1 CaoMyo1 CapMyo1 Ct_aMyo1 PtaMyo1 CacaMyo1B NadMyo1 278 278 279 278 278 278 279 278 278 278 278 278 278 278 278 278 278 278 282 280 277 278 277 276 276 276 277 278 278 278 282 278 277 278 268 278 278 278 278 278 278 278 276 279 278 279 277 280 284 283 284 278 278 287 285 284 286 285 285 285 285 285 292 283 281 286 284 288 286 286 284 282 281 278 310 320 330 340 350 360 370 380 390 400 ....|....|....|....|....|....|....|....|....|....|....|....|....|....|....|....|....|....|....|....| KAMQIIGLSQDEQDQIFRMLAAILWIGNISFIENEEGNAQVRDTSVTDFVAYLLQID-SPSLTKSLVERIVETNHGMKRGSVYHVPLNIVQATAVRDALA KAMQVIGLSQEEQDQIFRMLASILWIGNISFIENEEGNAQVRDTSVTDFVAYLLQID-SPSLTKALVERIVETNHGMKRGSVYHVPLNIVQATAVRDALA KAARTIGLSQEEQDQIFRQLAAILWIGNLTFVENEEGNAQVRDTSVSDFIAYLLQVE-SELLIKCIVERTVETNHGMKRGSVYHVPLNIVQANAVRDALA KAMQTIGLAQEEQDNIFRMLAAILWIGNISFIENEEGNAQIRDTSVTDFVAYLLQIN-SDILITAITERTVETNHGMKRGSVYHIPMNIVQATAVRDALA KAMQIIGLSQDEQDQIFRMLASILWIGNISFIENEEGNAQVRDTSVTDFVAYLLQVD-SPSLTKALVERIVETSHGMKRGSVYHVPLNIVQATAVRDALA KAMEIIGLSKEEQDQIFRMLAVILWIGNISFIENEEGNAQIRDTSVTAFVAYLLQVQ-EELLIKSLIERIIETGHGAKRGSTYHSPLNIVQATAVRDALA KAMQVIGLAQEEQDQIFRMLAAILWIGNISFIENEEGNAQVRDTSVTDFVAYLLQVD-SQSLIKALVERIVETNHGSRRGSVYHVPLNIVQATAVRDALA KAMKVIGLGQEEQDQIFRMLAAILWIGNISFIENEEGNAQVRDTSVTDFVAYLLQID-GQLLVKSLVERIMETNHGMKRGSVYHVPLNIVQADAVRDALA KAMRVIGLGQEEQDQIFRMLAAILWIGNVSFIENEEGNAQVRDTSVTDFVAYLLQID-SQLLIKSLVERIMETNHGMKRGSVYHVPLNIVQADAVRDALA KAMRVIGLEQDEQDQIFRMLAAILWIGNVSFIENEEGNAQVRDTSVTDFVAYLLQID-SQLLLKSLVERIMETSHGMKRGSVYHVPLNVVQADAVRDALA KAMRVIGLGQEEQDQIFRMLAAILWIGNVSFTENEEGNSQVRDTSVTDFVAYLLQID-SQLLIKSLVERIMETNHGMKRGSVYHVPLNIVQADAVRDALA KAMRVIGLGQEEQDQIFRMLAAILWIGNVSFIENEEGNAQVRDTSVTDFVAYLLQID-SQLLIKSLVERIMETNHGMKRGSVYHVPLNVVQADAVRDALA KAMRVIGLGQEEQDQIFRMLAAILWIGNVSFIENEEGNAQVRDTSVTDFVAYLLQID-SQLLTKSLVERIMETNHGMKRGSVYHVPLNIVQADAVRDALA KAMEVIGLHQEEQDQIFRMLAAVLWIGNVSFVENEEGNAQVRDTSVTDFVAYLLQID-APLLIKSLVERVMETNHGMRRGSVYHVPLNIVQATAVRDALA KAMQVIGLTQQEQDEIFRMLSAILWIGNITFIENEEGNAQVRDTSVTDFVAYLLQVD-SQLLIKSIVERTMETNHGMRRGSIYHVPLNIVQATAVKDALA KAMQVIGLSQDEQDQIFRMLSSILWVGNISFVENEEGNAQVRDTSVTDFVAYLLQID-SQVLMKALVERTMETSHGMRRGSVYHVPLNIVQATAVRDALA KAMQIIGLSQDEQDQIFRMLASILWIGNISFVENEEGNAQVRDTSVTDFVAYLLQID-ASILIKGLVERTMETSHGMRRGSVYHVPLNIVQANAAKDALA KAMQIIGLSQDEQDQIFRMLAFILWIGNISFVENEEGNAQVRDTSVTDFVAYLLQID-ASILIKGLVERTMETSHGMRRGSVYHVPLNIVQANAAEDALA KAMEIIGLDQNEQDQIFRMLAAILWIGNITFEENDEGNAQVRDTSVTDFVAYLLEVD-SPLLIKSLVERIMETSHGSQRGSVYHVPLNITQATAVRDALA KAMQVIGLSQDEQDQIFRMLSAILWIGNITFVENDEGNAQVADSSVTDFVAYLLQVD-AGVLVKSLVERIMETNHGMKRGSVYHVPLNRVQATAVRDALA EAMNVIGLSQAEQDEIFRLLSAILWIGNVTFMEDDEGNAKIADTSITDFVAYLLQVD-AGLLVKSLVERTIETTHGMRRGSIYNVPLNIVQATAVRDALA KAMGVIGLSQAEQDEIFRMLAAILWIGNITFAENDEGNAQVGDTTVTDFLAYLLQVD-PDLLLKSLVERVIETNHGMKRGSIYNVPLNIVQATAVRDALA KAMNTIGVTQQEQDEIFRMMAAILWIGNVTFVENDEGNAQVKDTSVTDFVAYLLQVD-SQILIKCLVERIMETNHGMKRGSVYHVPLNRVQATAVRDALA KAMQVIGLSQEEQDQIFRMLSAILWIGNVTFVENNEGNAEVRDTSVTDFVAYLMQVD-SGLLIKCLVERVMETGHGSRRGSVYHVPLNVVQATAVRDALA KAMDVIGLSQDEQDQIFRMLSAILWIGNISFIENDEGNSQVRDTSVTDFVAYLLQVD-ASLLIKCLVERVMETGHGARRGSVYHVPLNVVQANAVRDALA KAMEIIGLSQDEQDQIFRMLSAILWIGNISFVENDEGNSEVRDTSVTDFVAYLMQVD-SSLLVKCLVERVMETGHGSRRGSVYHVPLNVVQATAVKDALA KAMQVIGLSQDEQDQIFRMLAAILWIGNISFVENEEGNAQIRDTSVTDFVAYLLQVD-SNLLIKSLVERVMETGHGSKRGSIYHVPLNIVQASAVKDALA KAMQIIGLSQEEQDQIFRVLSSILWIGNITFAENEEGSAAVRDTSVTDFVAYLLQID-AALLTQSLVERTMETNHGMRRGSIYHVPLNVVQATAVRDALA NAMSVIGITQHEQDEVFRFLAAILWIGNISFTENEEGNAQVRDTSVTDFVAYLLQVD-SNLLVQALVERVMETNHGMKRGSIYHVPLNVVQATAVKDALA KAMQVIGITQQEQDELFRLLSAILWIGNISFTENEEGNAQVCDTSVTDFVAYLLQVD-AHFLTQALVERVMETNHGMKRGSVYHVPLNVVQATAVKDALA KAMQVIGLSQEEQDQIFRLLATILWIGNISFTEDEEGNAQVRDTSVTDFVAYLLQVD-SQLLIKCLVERVMETSHGMRRGSVYHVPLNIVQATAVKDALA EAMRTIGLVQEEQDQIFRMLAAILWIGNISFIENEEGNAQVRDTSVTDFVAYLLQVD-APLLIKCLVERIMQTSHGMKRGSVYHVPLNAVQATAARDALA EAMRTIGLVQEEQDQIFRMLAAILWIGNISFIENEEGNAQVGDTSVTDFVAYLLQVD-ASLLVKCLVERIMQTSHGMKRGSVYHVPLNPVQATAVRDALA EAMKTIGLVQEEQDQIFRMLAAILWIGNISFVENEEGNAQVRDTSVTDFVAYLLQVD-ASLLIKCLVERIMQTSHGMKRGSVYHVPLNAVQATAVKDALA ------GLGQEDQYQFFRMLAAIMWMANISFIENEEANAQVRDTSVTDFAAYILQVD-APFLIKCMVDRIMQTILGMKRGSVYHVPLNTVQATAVRDALA EAMRTIGLAQEEQDQIFRMLAAILWVGNISFIENEEGNAQVRDTSVTDFVAYLLQVD-ASLLIKCLVERIMQTSHGMKRGSVYHVPLNPVQATAVRDALA EAMRTIGLAQEEQDQIFRMLAAILWIGNISFIENEEGNAQIRDTSVTDFVAYLLQVD-ASLLMKCLVERIMQTSHGMKRGSVYHVPLNPVQATAVRDALA AAMRVIGLSQEEQDEIFRVLAAILWTGNITFMENDEGNAQVRDTSVTDFVAYLLQVD-SQLLVKSLVERIMETNHGMRRGSVYHVPLTTVHATAVKDALA KAMQVIGLSQEEQDQIFRLLAAILWIGNITFVEDEEGNAKIRDTSVTDFVAYLLQVE-SELLIKSLVERVLETNHGMRRGSVYHVPLNAVQATAVKDALA NAMRTIGLTKSEQDQIFRALAAILWIGNISFVENEAGNAEIRDKSVTTFVAYLLEVQ-EELLIKALIERIIETTHGAKRGSTYHSPLNIIQATAVRDALA NAMQTVGITQPEQDQIFRVLAAILWIGNITFVENEEGNAQVRDTSVTDFVAYLLQVD-SETLNKAVTERVMETSHGMRRGSVYNVPLNITQATAVRDALA SAMQTIGVTQDEQDQIFRILAAILWIGNITFMENDEGNSQVKDTSVTDFVAYLLQVD-SEILIKSITQRVMETSHGMKRGSVYEVPLNITQAIAVKDALA KAMQIIGIDQEQQDQIFRILAAILWIGNITFVEDEQGNSTVRDSSVTKFVAYLISVD-EQTLINSLITRVMTTMN-----ETYDIPLNPTQATAVKNALA AAMNLIGLTQAEQDNLFKLLAAILWIGNMSFVEDKDGNAAIADVSVPNFVAYLLEVD-AESVVKAVTQRIMETSRGGRRGSVYEVALNIAQATSVRDALA NAMNIIGLSKEEQDQIFRMLAAILWIGNISFIEDEQGNAQIRDTSVTTFVAYLLQVQ-EELLIKSLVERIIETGHGAKRGSTYHSPLNIVQATAVRDALA KAMQTIGLTQDEQDQVFRILAAILWIGNISFVENEDGNAQIRDSSVTDFVAYLLQVN-AQILTNSLIERIVETSHGSKRGSIYHVPLNITQATAVRDALS RAMDIIGLAKEEQDQIFRLLAAILWIGNISFVENEQGNAQIRDTSVTTFVAYLLQVQ-EDMLIKSLVERVIETSHGAKRGSTYHSPLNIVQATAVRDALA KAMRTIGMTQEEQDELFSLLAAVLWIGNVSFAEDAEGNSTIRDTSVTDFVAYLLQVD-SQILCKSLTERTLETNHGMKRGSIYHVPLNMTQATAVRDALA ASMTTIGLTQEEQDQVFRVLAAILWIGNITFVEDAEGNAQVRDTGVTDFVAYLLQVN-SEVLVKSIIERTMETSHGMRRGSVYHVPLNPVQATASRDALA KAMDVIGITDEERDQIFRILAGILWIGNITFTEDEEGNAAIADTSVTDFVAYLLQVD-AQTLCKSVVERTIQTFHGMRRGSIYHSPLNPVQATASRDALA KAMDTIGLSSEERDHVFRLLSAILWIGNVSFIEDEEGNSKIRDESVTNFVAYLLQVD-ATLLCKSLTERTMETSHGMKRGSIYHVPLNTTQATAVKDALA RAMETIGISQEEQDNIFRLLAAILWTGNISFVENEEGNAQIRDTSVTDFVAYLLQVD-SQMLIGSIIERIMETSHGSRRGSIYNVPLNIVQATAVKDALA KAMQVIGLTQEEQDEIFRVLAAILWIGNITFGENDEGNAEVRDTSVTDFVAYLLQVD-SQLLIKSLVERIMETNHGMRRGSVYHVPLNTVQATAVRDALA NAMKIIGLTQQEQDNIFRMLASILWIGNISFVEDENGNAAIRDDSVTNFAAYLLDVN-PEILKKAIIEKTIETSHGMRRGSTYHSPLNIVQATAVRDALA NAMNIIGLTQDEQDNIFKMLAAILWIGNISFVEDESGNAAIRDETVTQFVAYLLDVPSAEIMKKAIIERTIETTHGSRRGSTYHVPLNIVQATSVRDALA NAMNIIGLTQEEQDNIFRMLAAILWVGNISFVEDESGNAAIRDETVTQFVAYLLDVS-PEIMKKAIIERTIETTHGSRRGSTYHVPLNIVQATSVRDALA RAMEIIGLSQAEQDNIFRMLASILWIGNISFVEDESGNATIRDAGVTNFVAYLLEVS-PEILQKAIVERIIETSHGMRRGSTYHVPLNIVQATAVRDALA NAMNVIGLAQVEQDNIFRMLAAILWIGNISFVEDENGNAAVRDEGVTNFVAYLLEVN-AEILKKSIIERIIETSHGMRRGSTYHVPLNITQATAVRDALA NAMQVIGLTQQEQDSIFRMLASILWIGNISFVEDEHGNAAIRDESVTNFVGYLLDVS-GEIVKKSIIERTIETSFGSRRGSTYHSPLNIVQATAVRDALA RAMQVIGLSQEEQDNIFRMLASILWIGNISFVEDEHGNAAIRDESVTAFVAYLLDVN-AETLKTSLIQRVMQTSHGMRRGSTYHVPLNIVQATSVRDALA RAMQVIGLSQDEQDNIFRMLASILWVGNVSFVEDDNGNAAVRDESVTAFIAYLLDVD-AETLKTSLIQRVMQTSHGMRRGSTYHVPLNIVQATSVRDALA NAMKVIGITPQEQDHIFRMLAAILWVGNISYVEDEHGNAAIRDESVINFVAYLLETD-AESVKKSITEKIVQTSHGMRRGSTYHSPLNIVQATAVRDALA NAMNIIGLSQAEQDNIFRILASILWTGNISFVEDESGNAAIRDDSVTNFVAYLLDVN-AEILKKAITERTIETSHGMKRGSTYHVPLNIVQATAVRDALA QAMNTIGLSKAEQDNIFRSLASILWIGNISFVENEDGNAAIRDDTVTTFVAYLLEVD-ANVLKKSILERVIETSHGMRRGSTYHVPLNIVQATASRDALA SAMKIIGLTELEQNNIFRMIASILWIGNVSFVEDESGNAAIRDDSVTQFVAYLLEVN-PEILKKAIVERVIETTHGMRRGSTYHVPLNIVQATSVRDALA NAMNVIGLTQDEQDNIFRMLASILWIGNISFVEDESGNAAIRDDSVTNFVAYLLDVN-PEILKKAIIERTIETSHGMKRGSTYHVPLNIVQATAVRDALA KAMQIIGLTQDEQDNIFRVLASILWIGNISFVEDENGNAIPRDESVTNFVAYLLNVE-PVILKTAFVQRVMQTSHGMRRGSTYHVPLNIVQATAVRDALA NAMKIIGLTQQEQDNIFRMLAVILWIGNISFIEDENGNATIRDDSVTNFVAYLLDVN-SEILKKAIIERTIETSHGMRRGSTYHSPLNIVQATAVRDALA NAMNVIGLSQEEQDNIFRMLASILWIGNISFVEDESGNAAIRDDSVTNFVAYLLDVN-PEILKKAIIERTIETSHGMKRGSTYHVPLNIVQATAVRDALA NAMNVIGLSQDEQDNIFRMLASILWIGNISFIEDESGNAAIRDDSVTNFVAYLLDAN-PDILKKAIIERTIETSHGMKRGSTYHVPLNIVQATAVRDALA NAMQIIGLSQDEQDQIFRMLASILWVGNISFVEDENGNAAIRDESVTNFVAYLLGVN-TEILKKAIIERTIETSHGMRRGSTYHSPLNIVQATAVRDALA NAMSIIGLSQEEQDQIFRILSAILWIGNVLFVENEDANSAIGDQGVINYIAYLLQVD-AEQLAKSLTERIMETSHGMKRGSIYHVPLNLTQAYAVRDALA NAMKTIGLSQTEQDHIFRLLAAILWIGNISFIENEEGSAQIRDTSVTDFVAYLLEVD-SSLLITSIVERILETSHGTKRGSIYHTPLNIVQATAVRDGLA KAMQTIGLAQEEQDQIFRMLAAILWIGNISFIENEEGNAQVRDTSVTDFVAYLLQVD-AAVLIKALVERIMETSHGMRRGSVYHVPLNIVQATAVRDALA A window size 0 window size 3 900 2500 800 2000 leucine residues 700 occurence occurence 600 1500 1000 500 400 300 200 500 0 0 0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 5 0 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 conservation score [%] conservation score [%] 450 350 400 300 350 occurence 300 200 150 250 200 150 100 100 50 0 5 0 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 5 0 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 conservation score [%] B 4 kinesin-1 bits 3 2 I K L I N A S R M E N A QV R H Y E Q V Y 4 kinesin-5 bits 3 2 1 T K R Q R V I DG R A S F N D A S P KN L C D H I S MYF I F V M VG L FVR I L Y I V K Y F E 0 N N G R I SQ TVTT I 0 N LVDLAGSEK KTG A G LE AK IN SLS LG VIN LTD S S Q D N A R R I R T A N K NLVDLAGSE SG A E G IN SL LG VIN L K S LF 1 S VT I K R T SN A I DR SQ T S V M K L V L Q I G KL L ST L Q T Q N K AS TST YD I R M RA V N A N SA GK S K AS H I L N F H A C E QT FAR VMDV P R conservation score [%] S L M SS PE NC AQ RK K LQ A A G N N D SS GV Y E H N L I K N NTMVST Q V 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 0 50 T A F LQL S M T I H S I K Y T R I N FA S W I A G N F S V D H L H G K T G S A K E Q S Y A N I E S T N K Q K SIV V HMK T E R DA V F Q M E L G Q S Q I I S S T I D N L M G I K A A F G H R V KS N V K N S N X AS V KS S Q N I K N V F S L T MV G T E KN S C I I N M H G R C N Q G L NAP S QT T AS N A I Q 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 occurence 250 switch II Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 100 α4 C serine residues # alignment positions with a score in the respective range # alignment positions with a score in the respective range and with at least 1 CUG codon # CUG codons at alignment positions with a score in the respective range [%] CUG codons 0 10 20 30 40 50 60 70 80 90 100 0 10 20 30 40 50 60 70 80 90 100 Yarrowia lipolytica CLIB122 Wickerhamomyces anomalus NRRL Y-366 Wickerhamomyces ciferrii NRRL Y-1031 Cyberlindnera jadinii NBRC 0988 Kluyveromyces aestuarii ATCC 18862 Kluyveromyces wickerhamii UCD 54-210 Kluyveromyces marxianus var. marxianus KCTC 17555 Kluyveromyces lactis NRRL Y-1140 Lachancea thermotolerans CBS 6340 Lachancea waltii NCYC 2644 Eremothecium gossypii ATCC 10895 Eremothecium cymbalariae DBVPG#7215 Lachancea kluyveri NRRL Y-12651 Tetrapisispora phaffii CBS 4417 Vanderwaltozyma polyspora DSM 70294 Tetrapisispora blattae CBS 6284 Torulaspora delbrueckii CBS 1146 Zygosaccharomyces bailii ISA1307 Zygosaccharomyces rouxii CBS732 Candida castellii CBS 4332 Nakaseomyces bacillisporus CBS 7720 Candida glabrata CBS138 Nakaseomyces delphensis CBS 2170 Candida nivariensis CBS 9983 Candida bracarensis CBS 10154 Saccharomyces bayanus 623-6C Saccharomyces kudriavzevii IFO 1802 Saccharomyces mikatae IFO 1815 Saccharomyces cerevisiae S288c Saccharomyces paradoxus NRRL Y-17217 Saccharomyces arboricola H-6 Kazachstania africana CBS 2517 Kazachstania naganishii CBS 8797 Naumovozyma dairenensis CBS 421 Naumovozyma castellii NRRL Y-12630 Hanseniaspora vineae T02/19AF Komagataella pastoris GS115 Pichia kudriavzevii M12 Dekkera bruxellensis CBS 2499 Ogataea parapolymorpha DL-1 Kuraishia capsulata CBS 1993 Pachysolen tannophilus NRRL Y-2460 Candida tenuis NRRL Y-1498 Spathaspora passalidarum NRRL Y-27907 Spathaspora arborariae UFMG-19.1A Candida metapsilosis CANME Candida orthopsilosis Co 90-125 Candida parapsilosis Lodderomyces elongisporus NRRL YB-4239 Candida tropicalis MYA-3404 Candida maltosa Xu316 Candida albicans WO-1 Candida dubliniensis CD36 Scheffersomyces stipitis CBS 6054 Millerozyma farinosa CBS7064 Debaryomyces hansenii var. hansenii MTCC 234 Debaryomyces hansenii CBS767 Metschnikowia fructicola 277 Clavispora lusitaniae ATCC 42720 Meyerozyma guilliermondii ATCC 6260 [%] CUG codons Serine, min. 90% conservation score Leucine, min. 90% conservation score Serine, min. 80% conservation score Leucine, min. 80% conservation score Serine, min. 50% conservation sore Leucine, min. 50% conservation sore Yarrowia lipolytica CLIB122 Wickerhamomyces anomalus NRRL Y-366 Wickerhamomyces ciferrii NRRL Y-1031 Cyberlindnera jadinii NBRC 0988 Kluyveromyces aestuarii ATCC 18862 Kluyveromyces wickerhamii UCD 54-210 Kluyveromyces marxianus var. marxianus KCTC 17555 Kluyveromyces lactis NRRL Y-1140 Lachancea thermotolerans CBS 6340 Lachancea waltii NCYC 2644 Eremothecium gossypii ATCC 10895 Eremothecium cymbalariae DBVPG#7215 Lachancea kluyveri NRRL Y-12651 Tetrapisispora phaffii CBS 4417 Vanderwaltozyma polyspora DSM 70294 Tetrapisispora blattae CBS 6284 Torulaspora delbrueckii CBS 1146 Zygosaccharomyces bailii ISA1307 Zygosaccharomyces rouxii CBS732 Candida castellii CBS 4332 Nakaseomyces bacillisporus CBS 7720 Candida glabrata CBS138 Nakaseomyces delphensis CBS 2170 Candida nivariensis CBS 9983 Candida bracarensis CBS 10154 Saccharomyces bayanus 623-6C Saccharomyces kudriavzevii IFO 1802 Saccharomyces mikatae IFO 1815 Saccharomyces cerevisiae S288c Saccharomyces paradoxus NRRL Y-17217 Saccharomyces arboricola H-6 Kazachstania africana CBS 2517 Kazachstania naganishii CBS 8797 Naumovozyma dairenensis CBS 421 Naumovozyma castellii NRRL Y-12630 Hanseniaspora vineae T02/19AF Komagataella pastoris GS115 Pichia kudriavzevii M12 Dekkera bruxellensis CBS 2499 Ogataea parapolymorpha DL-1 Kuraishia capsulata CBS 1993 Pachysolen tannophilus NRRL Y-2460 Candida tenuis NRRL Y-1498 Spathaspora passalidarum NRRL Y-27907 Spathaspora arborariae UFMG-19.1A Candida metapsilosis CANME Candida orthopsilosis Co 90-125 Candida parapsilosis Lodderomyces elongisporus NRRL YB-4239 Candida tropicalis MYA-3404 Candida maltosa Xu316 Candida albicans WO-1 Candida dubliniensis CD36 Scheffersomyces stipitis CBS 6054 Millerozyma farinosa CBS7064 Debaryomyces hansenii var. hansenii MTCC 234 Debaryomyces hansenii CBS767 Metschnikowia fructicola 277 Clavispora lusitaniae ATCC 42720 Meyerozyma guilliermondii ATCC 6260 # CUGs 770 # CUGs at positions with a conservation score of ≥50% Yarrowia lipolytica CLIB122 Wickerhamomyces anomalus NRRL Y-366 Wickerhamomyces ciferrii NRRL Y-1031 Cyberlindnera jadinii NBRC 0988 Kluyveromyces aestuarii ATCC 18862 Kluyveromyces wickerhamii UCD 54-210 Kluyveromyces marxianus var. marxianus KCTC 17555 Kluyveromyces lactis NRRL Y-1140 Lachancea thermotolerans CBS 6340 Lachancea waltii NCYC 2644 Eremothecium gossypii ATCC 10895 Eremothecium cymbalariae DBVPG#7215 Lachancea kluyveri NRRL Y-12651 Tetrapisispora phaffii CBS 4417 Vanderwaltozyma polyspora DSM 70294 Tetrapisispora blattae CBS 6284 Torulaspora delbrueckii CBS 1146 Zygosaccharomyces bailii ISA1307 Zygosaccharomyces rouxii CBS732 Candida castellii CBS 4332 Nakaseomyces bacillisporus CBS 7720 Candida glabrata CBS138 Nakaseomyces delphensis CBS 2170 Candida nivariensis CBS 9983 Candida bracarensis CBS 10154 Saccharomyces bayanus 623-6C Saccharomyces kudriavzevii IFO 1802 Saccharomyces mikatae IFO 1815 Saccharomyces cerevisiae S288c Saccharomyces paradoxus NRRL Y-17217 Saccharomyces arboricola H-6 Kazachstania africana CBS 2517 Kazachstania naganishii CBS 8797 Naumovozyma dairenensis CBS 421 Naumovozyma castellii NRRL Y-12630 Hanseniaspora vineae T02/19AF Komagataella pastoris GS115 Pichia kudriavzevii M12 Dekkera bruxellensis CBS 2499 Ogataea parapolymorpha DL-1 Kuraishia capsulata CBS 1993 Pachysolen tannophilus NRRL Y-2460 Candida tenuis NRRL Y-1498 Spathaspora passalidarum NRRL Y-27907 Spathaspora arborariae UFMG-19.1A Candida metapsilosis CANME Candida orthopsilosis Co 90-125 Candida parapsilosis Lodderomyces elongisporus NRRL YB-4239 Candida tropicalis MYA-3404 Candida maltosa Xu316 Candida albicans WO-1 Candida dubliniensis CD36 Scheffersomyces stipitis CBS 6054 Millerozyma farinosa CBS7064 Debaryomyces hansenii var. hansenii MTCC 234 Debaryomyces hansenii CBS767 Metschnikowia fructicola 277 Clavispora lusitaniae ATCC 42720 Meyerozyma guilliermondii ATCC 6260 #CUGs at identical alignment positions 316 100 31 10 3 1 Schizosaccharomyces pombe 972hSchizosaccharomyces cryophilus NRRL Y-48691 Schizosaccharomyces octosporus yFS286 Schizosaccharomyces japonicus yFS275 Yarrowia lipolytica CLIB122 Wickerhamomyces anomalus NRRL Y-366 Wickerhamomyces ciferrii NRRL Y-1031 Cyberlindnera jadinii NBRC 0988 Kluyveromyces aestuarii ATCC 18862 Kluyveromyces wickerhamii UCD 54-210 Kluyveromyces marxianus var. marxianus KCTC 17555 Kluyveromyces lactis NRRL Y-1140 Lachancea thermotolerans CBS 6340 Lachancea waltii NCYC 2644 Eremothecium gossypii ATCC 10895 Eremothecium cymbalariae DBVPG#7215 Lachancea kluyveri NRRL Y-12651 Tetrapisispora phaffii CBS 4417 Vanderwaltozyma polyspora DSM 70294 Tetrapisispora blattae CBS 6284 Torulaspora delbrueckii CBS 1146 Zygosaccharomyces bailii ISA1307 Zygosaccharomyces rouxii CBS732 Candida castellii CBS 4332 Nakaseomyces bacillisporus CBS 7720 Candida glabrata CBS138 Nakaseomyces delphensis CBS 2170 Candida nivariensis CBS 9983 Candida bracarensis CBS 10154 Saccharomyces bayanus 623-6C Saccharomyces kudriavzevii IFO 1802 Saccharomyces mikatae IFO 1815 Saccharomyces cerevisiae S288c Saccharomyces paradoxus NRRL Y-17217 Saccharomyces arboricola H-6 Kazachstania africana CBS 2517 Kazachstania naganishii CBS 8797 Naumovozyma dairenensis CBS 421 Naumovozyma castellii NRRL Y-12630 Hanseniaspora vineae T02/19AF Komagataella pastoris GS115 Pichia kudriavzevii M12 Dekkera bruxellensis CBS 2499 Ogataea parapolymorpha DL-1 Kuraishia capsulata CBS 1993 Pachysolen tannophilus NRRL Y-2460 Candida tenuis NRRL Y-1498 Spathaspora passalidarum NRRL Y-27907 Spathaspora arborariae UFMG-19.1A Candida metapsilosis CANME Candida orthopsilosis Co 90-125 Candida parapsilosis Lodderomyces elongisporus NRRL YB-4239 Candida tropicalis MYA-3404 Candida maltosa Xu316 Candida albicans WO-1 Candida dubliniensis CD36 Scheffersomyces stipitis CBS 6054 Millerozyma farinosa CBS7064 Debaryomyces hansenii var. hansenii MTCC 234 Debaryomyces hansenii CBS767 Metschnikowia fructicola 277 Clavispora lusitaniae ATCC 42720 Meyerozyma guilliermondii ATCC 6260 100/100/100 100/100/100 587 Ma 100/100/100 100/100/88 100/100/100 100/100/100 235 Ma 100/100/100 100/100/52 100/100/100 Downloaded from http://gbe.oxfordjournals.org/ by guest on October 6, 2014 87/100/100 100/100/100 56/100/- 235 220 100 175 135 Ma amino acid substitutions per site 0.2 0 190 200 300 400 500 600 700 800 # CTG codons
© Copyright 2024