- Research article
- Open Access
Emergence and evolution of yeast prion and prion-like proteins
BMC Evolutionary Biology volume 16, Article number: 24 (2016)
Prions are transmissible, propagating alternative states of proteins, and are usually made from the fibrillar, beta-sheet-rich assemblies termed amyloid. Prions in the budding yeast Saccharomyces cerevisiae propagate heritable phenotypes, uncover hidden genetic variation, function in large-scale gene regulation, and can act like diseases. Almost all these amyloid prions have asparagine/glutamine-rich (N/Q–rich) domains. Other proteins, that we term here ‘prionogenic amyloid formers’ (PAFs), have been shown to form amyloid in vivo, and to have N/Q-rich domains that can propagate heritable states in yeast cells. Also, there are >200 other S.cerevisiae proteins with prion-like N/Q-rich sequence composition. Furthermore, human proteins with such N/Q-rich composition have been linked to the pathomechanisms of neurodegenerative amyloid diseases.
Here, we exploit the increasing abundance of complete fungal genomes to examine the ancestry of prions/PAFs and other N/Q-rich proteins across the fungal kingdom. We find distinct evolutionary behavior for Q-rich and N-rich prions/PAFs; those of ancient ancestry (outside the budding yeasts, Saccharomycetes) are Q-rich, whereas N-rich cases arose early in Saccharomycetes evolution. This emergence of N-rich prion/PAFs is linked to a large-scale emergence of N-rich proteins during Saccharomycetes evolution, with Saccharomycetes showing a distinctive trend for population sizes of prion-like proteins that sets them apart from all the other fungi. Conversely, some clades, e.g. Eurotiales, have much fewer N/Q-rich proteins, and in some cases likely lose them en masse, perhaps due to greater amyloid intolerance, although they contain relatively more non-N/Q-rich predicted prions. We find that recent mutational tendencies arising during Saccharomycetes evolution (i.e., increased numbers of N residues and a tendency to form more poly-N tracts), contributed to the expansion/development of the prion phenomenon. Variation in these mutational tendencies in Saccharomycetes is correlated with the population sizes of prion-like proteins, thus implying that selection pressures on N/Q-rich protein sequences against amyloidogenesis are not generally maintained in budding yeasts.
These results help to delineate further the limits and origins of N/Q-rich prions, and provide insight as a case study of the evolution of compositionally-defined protein domains.
Yeast prions are propagating alternative states of proteins. These states can be transmitted sustainably in the yeast Saccharomyces cerevisiae during budding, mating or laboratory infection protocols. The first well-characterized yeast prions, that underlie the [PSI+] and [URE3] prion states, are propagating amyloids (i.e., fibrillar beta-sheet aggregates) of the proteins Sup35p and Ure2p. The protein Sup35p is part of the translation termination complex. Formation of [PSI+] prions reduces the efficiency of translation termination and increases levels of nonsense-codon read-through [1, 2]. Such read-through has been shown to be a potential mechanism to uncover cryptic genetic variation [3, 4]. [URE3] causes upregulation of poor nitrogen source usage, even when rich sources are available [5–7]. Prion variants may be considered as diseases of S. cerevisiae in some contexts [8, 9]. A more recently discovered example, the [MOT3+] prion, has been shown to govern acquisition of multicellularity in S. cerevisiae . There are now at least 10 known prions of S. cerevisiae that are propagated by amyloids [11, 12].
A common compositional feature of almost all amyloid-based yeast prions is bias for asparagine (N) and/or glutamine (Q) residues [11, 12]. A majority of them are N-rich (6/10 at the time of this analysis), rather than Q-rich. Bioinformatic surveys have revealed the existence of hundreds of proteins with such N/Q-richness in S. cerevisiae and diverse other fungi [13–15]. Evolutionary analysis showed that the [PSI+] prion N/Q bias is conserved across fungal clades that diverged >1 billion years ago, with only eight other S. cerevisiae proteins showing similar, phylogenetically deep patterns of N/Q bias conservation . The [URE3] prion domain is unique to Saccharomycetes, with different parts of the domain demonstrating purifying selection (i.e., significant conservation of amino-acid identity from examination of codon position mutation rates), and variation in N/Q bias between clades [14, 16].
The peculiar composition of known prions has been exploited to computationally detect likely prions that were then tested experimentally for prion-forming ability . Tests for in vitro and in vivo amyloid formation were combined with a Sup35 prion assay, wherein predicted prion-forming domains were fused to the C-terminal part of the Sup35p protein, and these constructs were then tested for the ability to produce [PSI+]-like states in yeast cells . About twenty novel ‘prionogenic’ proteins were identified. The results from this survey have been used to train other algorithms to predict prion domains bioinformatically  (PrionW, PLAAC [19, 20]). On a related note, ‘scrambled’ forms of the Ure2p and Sup35p prion domains, that maintain the same amino-acid composition, can form prions in S. cerevisiae, indicating that prion formation is primarily determined by composition but not by specific sequence features [21, 22]. Building on these analyses, an amino-acid propensity scale for prion formation was developed and incorporated into the PAPA method for prion prediction [23, 24].
Putative prion domains from other Saccharomycetes (but not from fungal clades outside of this one) can make prions in S. cerevisiae or in their own cells, although this ability is sporadic [25–30], and can rely on small changes in the protein sequence . Conversely, the full-length non-yeast protein CPEB from the sea hare Aplysia californica can form prions in S. cerevisiae, albeit with much less efficiency than native prions [31, 32]. Mutational experiments indicate that many N/Q-rich domains in S. cerevisiae may only be a small number of sequence mutations away from prion-forming ability, implying that natural selection may only act to keep aggregation propensities sufficiently low ; this may be an under-appreciated effect in the analysis of mammalian prion disease mutations [34, 35].
Several human proteins have prion-like N/Q-rich domains that have been directly linked to neurodegenerative diseases. Cytoplasmic aggregates of the RNA-binding protein FUS, which contains a Q-rich domain, are implicated in amyotrophic lateral sclerosis, and its aggregation has been re-capitulated in an induced S. cerevisiae proteinopathy . Mutations in two yeast-prion-like proteins hnRNPA2B1 and hnRNPA1 initiate neurodegenerative disease in humans through amyloid formation . HNRPDL has a yeast-prion-like domain, and has been linked to development of limb-girdle muscular dystrophy 1G . Also, pathogenic proteins in at least nine other neurodegenerative disorders have disease-linked poly-Q expansions. Thus, the degree of conservation and variation of yeast prion domains has implications not just in fungi, but for human diseases as well.
Here, we probe how prion and prion-like proteins have evolved across the fungal kingdom. We discover that the ancestors of N-rich prion formers emerged during Saccharomycetes speciation, in tandem with a general dramatic increase in the number of N-rich proteins. Conversely, more ancient prion biases are Q-rich, at least back to the last common ancestor of fungi. Some fungal clades have very few N/Q-rich proteins, and in some cases likely lose them en masse. We find evidence that recent emergence of large populations of N/Q-rich proteins in Saccharomycetes may be partly due to mutational tendencies leading to more frequent initiation and elongation of poly-N runs. Variation in these mutational tendencies in Saccharomycetes is correlated with the population sizes of prion-like proteins, thus implying that selection pressure on N/Q-rich protein sequences to prevent their formation of amyloids is not generally maintained across Saccharomycetes.
We downloaded complete proteomes of 169 fungi from various sources as listed in Additional file 1: Table S1.
A phylogeny of fungi was obtained from NCBI Taxonomy (http://www.ncbi.nlm.nih.gov). Organismal phylogenetic trees were drawn using phyloT (http://phylot.biobyte.de) to generate a Newick format file, which was then input into Phylodendron (http://iubio.bio.indiana.edu/treeapp/treeprint-form.html). Orthologs for all of the S. cerevisiae proteins in all of the other fungi were calculated using the bi-directional best hits method. Families of paralogous proteins were determined using the CDHIT program (http://cd-hit.org). Removing the small numbers of putative paralogous duplications (identified using the CDHIT program) for N/Q-rich proteins has no effect on these observed trends for reported in this paper.
Prion and prionogenic proteins
Prion and prionogenic sequence sets for S. cerevisiae were taken from the PrionHome database [11, 12]. Here, we analyze as groups: (i) the set of known prions, and (ii) a larger set made from these known prions plus other ‘prionogenic’ amyloid-forming proteins (PAFs). We initially included the non-N/Q-rich prion protein Mod5p in our analysis, which underlies the [MOD+] prion state , to check whether it acquires N/Q-rich domains in other clades (Table 1); but discovered that it does not gain any. The PAFs set includes the prionogenic proteins from the analysis of Alberti, et al.  that have been shown to form prions, through a SUP35C prion assay in conjunction with evidence for in vivo amyloid formation by the full-length proteins from the other assays. The list of PAFs is as follows: (UniProt IDs and standard names, well-characterized ‘known’ prion proteins have an asterisk): P05453, Sup35p*; P07884, Mod5p*; P09547, Swi1p*; P14922, Cyc8p*; P23202, Ure2p*; P25367, Rnq1p*; P32432, Sfp1p*; Q08972, New1p*; P54785, Mot3p*; Q02629, NUP100p*; P32588, Pub1p*; P40070, Lsm4p; P14907, Nsp1p; P18494, Gln3p; P32770, Nrp1p; P38180, YBL081W; P38216, YBR016W; P38429, Sap30p; P38691, Ksp1p; P40356, Pgd1p; P53894, Cbk1p; Q05166, Asm4p; Q08925, Mrn1p; Q12139, YPR022C; Q12221, Puf2p; Q12224, Rlm1p; Q12361, Gpr1p. Further evidence for the importance/relevance of many of the PAF proteins to the prion phenomenon and other aggregation-dependent phenomena in yeast is continuing to accumulate. For example, Lsm4 amyloids can both act as a [PSI+] prion inducer and prion clearer (the latter when overexpressed) [40–43], and can underlie the aggregation of P-bodies . Also, fragments of Sap30p and Gpr1p have been shown to act as prion inducers . The N/Q-rich regions of Nsp1p are important in mediating nucleoporin hydrogel formation, and interact in trans with the Sup35p prion domain .
Prion-like proteins are defined in two ways: (a) through identification of compositional bias for N and Q residues, and (b) through application of algorithms designed to predict prion-forming domains.
N/Q-rich Proteins (NQPs) have N/Q-rich domains in them. For S. cerevisiae and all other fungal proteomes, N/Q-rich domains were determined using the LPS compositional-bias binomial probability minimization algorithm (with a maximum binomial P-value threshold of 1×10−10, derived from analyzing known prion-determinant domains [12, 14, 15, 47–49]), but testing three different criteria for expected amino-acid composition: (i) using the amino-acid composition for structured protein domains for each proteome (determined using blastp search against the ASTRAL database  with e-value threshold <1e–04 ); (ii) using equal expected amino-acid composition (i.e. =0.05), and (iii) using the average amino-acid composition of all the proteomes examined. Results are reported for criterion (ii), but the same dominant trend for large-scale emergence of N-rich proteins in Saccharomycetes is observed regardless of the criterion used.
Also, we applied the PAPA and PLAAC prion prediction programs [19, 23, 24] to all the complete proteomes. The PAPA algorithm using an experimentally derived prion propensity score combined with explicit consideration of the intrinsic disorder. For PAPA, the default threshold for prion prediction was used. PLAAC uses a Hidden Markov Model trained on the composition of known prion-forming domains, which all have a pronounced bias for N and/or Q residues, and were all known or predicted intrinsically-disordered domains. For PLAAC, we used as a threshold the lowest COREscore value for the known prion-forming proteins (20.58 for Sfp1p). Also, the PrionW webserver  was applied to the PAF data set and their orthologs. Intrinsic disorder was annotated using IUPRED .
Results and discussion
Firstly, we examine the evolutionary ancestry of the prions and other prionogenic amyloid-forming (PAF) protein sequences across the fungal kingdom. Then, we describe a dramatic large-scale emergence of N-rich prion-like proteins in the budding yeasts (Saccharomycetes), and how this contrasts with the notable lacks/losses of prion-like proteins in particular fungal clades/species. We show how the N-rich protein emergence in Saccharomycetes species is a striking trend that sets them apart from other fungi. We analyze how this trend is linked to recent mutational tendencies in this clade, and discuss the implications of our observations for prion formation.
In general, we examine trends at three evolutionary depths: (i) within the class of budding yeasts Saccharomycetes (also known as Hemiascomycetes), (ii) within the phylum of the sac fungi Ascomycota but beyond the Saccharomycetes, and (iii) outside of Ascomycota in other fungi (Fig. 1). We examine the conservation of N/Q-richness and predicted prion status (annotated as described in Methods).
Evolutionary origins of the ancestral sequences of prions and other prionogenic amyloid formers
We find conspicuously distinct evolutionary ancestry for prions/PAFs if we separate them into N-rich and Q-rich cases. They are designated N-rich if they have smaller P-value from the LPS algorithm for N bias than for Q bias, and vice versa for Q-rich cases (see Methods for details). N-rich sequences dominate the set of prions/PAFs (18/26 cases), with almost all of these arising evolutionarily within the Saccharomycetes (Fig. 2). If these N-rich prions/PAFs arose earlier in fungal evolution, they had a Q-rich sequence which subsequently became N-rich within Saccharomycetes (Fig. 2). All but one of the N-rich prion/PAF domains arose within an evolutionarily short time frame, after the last common ancestor of Saccharomycetes, and before the whole genome duplication (WGD) event that occurred within Saccharomycetes. (A small minority (5/26) of the PAFs have ohnologs (i.e. WGD gene duplications), of which four maintain N/Q-rich predicted prion status.) The ancestor of the prion protein Rnq1p appeared in the same evolutionary epoch. Rnq1p is one of three mostly Q-rich prions that arose within Saccharomycetes, whose [PIN+]/[RNQ+] prion is required for the induction in wild strains of the [PSI+] prion made from Q-rich Sup35p [40, 53].
All originally Q-rich prion/PAF sequences arose before the last common ancestor of the Saccharomycetes, either within Ascomycota, or further back in fungal evolution (Fig. 2). We can see that ancestral Q-rich composition, back beyond the last common ancestor of the Saccharomycetes, occurs for 6/10 of the N/Q-rich prions (bold names in Fig. 2) and 12/26 of the PAFs overall (names in bold or italics, Fig. 2). For such Q-rich sequences of ancient origin, switching to N bias arises only rarely outside of Saccharomycetes yeasts (in 5-6 % of all orthologs) (Fig. 2). A particularly notable case is the Q-rich prion protein CYC8, which is a part of a global transcription repressor complex that controls the expression of ~7 % of yeast genes . It is deeply conserved as a prion-like domain across diverse Ascomycota and in a few other fungal clades. Almost all of the putative prion-like domains for orthologs of CYC8 are Q-rich (66/67, i.e. with a single N-rich case). Such deep conservation of Q-rich prion-like domains (to before the last common ancestor of Ascomycota) for the prion proteins Sup35p, Cyc8p, Swi1p and others may be linked to a function other than prion formation. The prion domains of Sup35p and Ure2p have been shown to also have non-prion-forming functions [55, 56]. However, Pub1/Tia1 functions in stress granule assembly in mammals through aggregation mediated by its prion-like domain, a phenomenon that also arises in single-celled yeast [57, 58]. It is particularly intriguing that the Q-richness of such prion sequences has been deeply conserved across diverse fungal clades, since it has been shown through mutation experiments that Q-richness tends to lead to formation of toxic non-amyloid conformations in S. cerevisiae, whereas N-richness tends to produce benign propagating amyloids . The chromatin remodelling factor Swi1p has distinct N-rich and Q-rich domains . The N-rich domain, which has arisen recently in the evolution of S. cerevisiae and closely related yeasts, is required for the formation of the [SWI+] prion, which causes a partial loss of function phenotype [61, 62]; the deeply conserved Q-rich domain can modify aggregation patterns . The N-rich and Q-rich regions function in causing a gain in sensitivity to Na+/Li + ions .
The conservation of the prion-like character of prions/PAFs for the three studied taxonomic groups is summarized in Table 1. For most prions/PAFs with orthologs ‘beyond Ascomycota’, or ‘within Ascomycota beyond Saccharomycetes’, prion-like domains are observed, either as N/Q-rich annotations or algorithmic prion predictions (Table 1). There is substantial agreement between prion/PAF ortholog annotations for N/Q-rich domains and algorithmic prion predictions (discussed in Additional file 2: Text S1).
Emergence of N-rich proteins in Saccharomycetes yeasts
Is there any link between patterns of emergence of prion ancestors and prion-like proteins across the fungal kingdom? To answer this, we annotated all N/Q-rich proteins across the whole proteomes of the >160 fungal species under study. In doing this, we discovered a dramatic expansion in the number of N/Q-rich proteins in the Saccharomycetes clade, with all other clades having on average substantially fewer (Additional file 3: Figure S1 in detail, with a schematic summary in Fig. 3a). This is due to emergence of large numbers of N-rich prion-like domains, during Saccharomyetes evolution (Additional file 4: Figure S2). The trend for the evolutionary emergence of N-rich domains in prions and other PAFs is thus linked to a more general large-scale trend during Saccharomycetes evolution. This evolutionary trend for N/Q-rich domain genesis is replicated for numbers of prion predictions by the PAPA and PLAAC programs (Fig. 3b). N/Q-rich domains which arose within Saccharomycetes have significant functional linkage to transcription regulation, as determined by analysis of Gene Ontology process category enrichments (Additional file 5: Table S2, corrected P-values <0.05). N-rich prion-like domains (for example the one in Swi1) thus may have a specific functional influence on the recent evolutionary dynamics of transcriptional regulation pathways in the Saccharomyces genus. The Gene Ontology category enrichments (amongst many others) are also observed for the N/Q-rich proteins that occur beyond Saccharomycetes (with corrected P-values ≤1e–26). At least seven of the known prions/PAFs likewise function in regulation of transcription (the prions MOT3 , SWI1 , CYC8 , SFP1  and the PAFs GLN3 , PGD1  and RLM1 ). N-bias is not the only bias to become prominent in Saccharomycetes; there is also an emergence of more D-, E- and K-rich proteins (Additional file 6: Table S3, discussed in more detail below).
Many fungal clades and species have few prion-like proteins and in some cases have likely lost them in their recent evolution
Some clades have very few N/Q-rich proteins or predicted prions and thus fewer possible N/Q-rich prions (Additional file 3: Figure S1 and Fig. 3). These include the Eurotiales (containing the filamentous fungi genus Aspergillus), the fission yeasts Schizosaccharomycetes and the Agaricomycetes (the class including the mushrooms). Also, species in these clades contain few orthologs of the known prions/PAFs (Table 2). The dearth of likely prion-forming and prion-like proteins may be perhaps due to some mechanistic intolerance to their aggregation/propagation. Indeed, they may be too easily propagated to daughter cells in some fungal species, and thus subject to greater selection pressure on their sequences against formation of prion-forming domains. For Eurotiales, the most parsimonious explanation is that N/Q-rich domains and possible prions have been lost since their last common ancestor with other Ascomycota relatives.
Notably, two of the clades with fewest overall N/Q-rich proteins or prion predictions (Eurotiales and Schizosaccharomycetes) have some of the largest numbers of non-N/Q-rich prion predictions (Additional file 7: Figure S3 and Fig. 3c; a strict N/Q bias threshold of 1×10−5 is used). Although these are speculative predictions, this may indicate that these clades harbour differently composed cohorts of functional amyloid-forming proteins.
Substantial losses of prion-like proteins also occur in two individual Saccharomycetes species (Additional file 3: Figure S1). These are the non-WGD species Ogataea parapolymorpha and Ashbya gossypii. The species O. parapolymorpha has the lowest level of conservation of known prions and PAFs in Saccharomycetes (Table 2). O. parapolymorpha is an atypical thermotolerant yeast with a relatively high GC% (percentage guanidine + cytidine) genome that can grow on methanol, acquiring large numbers of cellular peroxisomes in the process; it also has an unusual thermotolerance mechanism linked to production of trehalose, a sugar normally found in insects. The filamentous yeast A. gossypii has undergone substantial genome evolution since divergence from its close relative E. cymbalariae, gaining higher GC% and losing transposons and 10 % of its genome size . The high GC% of these two genomes (51 % for A. gossypii and 48 % for O. parapolymorpha; the highest and third-highest GC% values of the Saccharomycetes species examined) may be a contributing factor to the loss of N-rich prion-like domains in particular: N codons are one sixth guanidine/cytidine (whereas Q codons are half). Correlation with variation in GC% is discussed in more detail below. O. parapolymorpha and A. gossypii share three conserved prion proteins (Cyc8p, Sup35p and Sfp1p). Sup35p and Sfp1p are functionally linked prions that exert control over translation accuracy . This is thus evidence for selection to maintain a small core of prions, despite many others being lost.
A distinctive trend for formation of prion-like proteins in Saccharomycetes yeasts
Since prion and prion-like proteins are intrinsically disordered, we surmised that maybe the trends in variation of N/Q-rich proteins or predicted prions are due to a more general trend for variation in the number of intrinsically-disordered proteins (IDPs) across fungal evolution. Thus, we compared the numbers of IDPs with the numbers of N/Q-rich proteins and prion predictions, for each proteome (Fig. 4). In general, we find some degree of correlation between numbers of IDPs and prion-like proteins. This may be because many intrinsically disordered regions (including those that contain N/Q-rich regions) are evolving neutrally or nearly neutrally, with little negative selection pressure. Also, such intrinsically disordered regions (including those that contain N/Q-rich regions) may have some organizational function in the cell that makes their precise amino-acid composition unimportant. In Fig. 4, we find a distinct trend for Saccharomycetes yeasts that sets them apart from non-Saccharomycetes. Saccharomycetes occupy the lowest portion of the scatter plots with a shallow correlation, where they segregate from all the other fungi. The highest numbers of N/Q-rich proteins (400+, or 150+ for those also predicted as prions by the PAPA and PLAAC programs), are in the genera Candida and Tetrapisispora, and in the species N. dairensis and L. elongisporus. This is the case whether we consider compositionally-defined N/Q-rich proteins, or PAPA/PLAAC-predicted prions that are N/Q-rich (Fig. 4a-b). The shallow distinct correlation for Saccharomycetes implies that additional N/Q-rich domains tend to arise by mutation without formation of many additional IDPs, and this is at a rate that sets them apart from other fungi. Notably, considering only non-N/Q-rich algorithmic prion predictions makes the trend for Saccharomycetes less distinct, implying that the trend observed is primarily due to N/Q-richness, and that numbers of non-N/Q-rich prion predictions are more correlated with numbers of IDPs generally (Fig. 4c). Algorithmically predicted non-N/Q-rich yeast prions are a largely untested cohort, and their exact design principles have yet to be determined.
Motivated by this distinct trend for Saccharomycetes yeasts, and given that N/Q-rich proteins (NQPs) are compositionally defined, we checked whether their numbers in fungal proteomes are correlated with other compositional characteristics (Table 3). We defined two special types of N and Q composition, % ‘lone’ N and Q and % ‘run’ N and Q. ‘Lone’ N and Q do not occur in homopeptide runs and are surrounded on either side by ≥2 non-N/Q residues. ‘Run’ N and Q occur in runs of 3–5 residues. Both ‘lone’ and ‘run’ N and Q residues are counted only from proteins that are not N/Q-biased (using a strict P-value threshold of <1×10−5), and that are not predicted to be prions by PAPA or PLAAC. Overall in fungi, and in other fungal clades examined for comparison, we find significant correlations for NQP numbers arising out of Q percentages (Table 3). However, when we specifically examine the Saccharomycetes clade, we see a different situation. There is a prominent correlation for % of run N, with other compositional traits having less significant or non-significant correlations (Table 3). These results indicate that the sizes of populations of N/Q-rich proteins in Saccharomycetes yeasts is directly linked to a mutational tendencies for more N residues, particularly in poly-N runs. Lower %GC (% guanidine + cytidine in the DNA) may lead to a higher proportion of Ns for initiation of runs (since N codons are AT-rich/GC-poor). For Saccharomycetes, we see also an increase in K-, D- and E-rich proteins, and a depletion of A-, G-, P- and R-rich proteins compared to other fungal clades (Additional file 6: Table S3). Notably, K residues are encoded by the most AT-rich set of codons in the genetic code, while A, G, P and R comprise the amino-acid residues encoded by GC-rich codons. Also, in line with this observation, we see that within Saccharomycetes, GC% has high positive correlation with the occurrence of alanine-rich proteins, and high negative correlation with the occurrence of isoleucine- and lysine-rich proteins (Additional file 8: Table S4). These latter two amino acids are encoded by AT-rich codons. Thus, GC% is an important contributor to the occurrence of compositionally biased regions in Saccharomycetes proteins, including N/Q-rich regions.
However, %GC is not correlated with numbers of NQPs across the fungal kingdom (Table 3). Indeed, there are several notable clades that have similar %GC but drastically different % NQPs. For example, Schizosaccharomyces and Saccharomyces species have similar %GC (~37 % versus ~39 %), but Schizosaccharomyces have much fewer NQPs and prion predictions (Fig. 3 and Additional file 3: Figure S1). A similar situation arises in the Basidiomycota, where Ustilaginomycotina have much more NQPs and prion predictions than Agaricomycotina, but have similar %GC (~54 % versus ~52 %). Thus, the precise nature of the selection pressures that contribute to the populations of NQPs remains to be elucidated fully. It has been shown that part of the Ure2p prion-forming domain is under purifying selection in Saccharomycetes, whereas another part of the domain varies widely in its N and Q composition . Whether in some sequences this variation is partly caused by diversifying or positive selection (i.e., significantly increased amino-acid mutation rate from examination of codon position mutation rates) will require further developments in molecular evolution models for biased sequences.
The evolutionary vista for the ancestors of prion and prion-like proteins changed substantially in Saccharomycetes budding yeasts. During Saccharomycetes evolution, large-scale formation of N-rich regions occurred. This thus may have provided a trigger for the expansion and development of the prion phenomenon and so consequently were born the ancestral sequences of the prions Ure2p, Mot3p and New1p, and other N-rich PAF proteins of S. cerevisiae. Thus, new prion domains could have initially arisen from the formation of sufficiently long poly-N/Q tracts (particularly poly-N tracts) . Certain individual newly-formed N-rich domains subsequently have been maintained to perform a function that may or may not be related to prion formation . Other factors being conducive, these evolutionarily novel N-rich domains could evolve to produce benign propagating amyloids in S. cerevisiae . Also within the same epoch (before the whole genome duplication in budding yeasts), the Rnq1p protein required for [PSI+] induction has arisen as a novel protein. Variation in recent mutational tendencies for more N residues, particularly in the form of poly-N tracts, is correlated with population sizes of N/Q-rich proteins in individual Saccharomycetes yeast species. Given the correlation that we see between numbers of N/Q-rich proteins and numbers of short poly-N tracts in other proteins, these results suggest that there is no clade-wide maintenance of selection pressure on N/Q-rich protein sequences to prevent N/Q-rich protein aggregation. This may be either because in many species there is cellular machinery to prevent/handle them effectively, or because they do not often enough tend to aggregate. In the amoeba Dictyostelium, there are large numbers of N/Q-rich proteins, and experiments on Sup35p aggregation indicate that there are cellular mechanisms preventing their aggregation generally . Such mechanisms may also allow larger populations of N/Q-rich proteins in the Tetrapisispora and Candida clades, and are of interest for the analysis of diseases in humans that are linked to prion-like proteins or poly-Q repeat expansions, such as Huntington’s disease. Also, the tendency to form poly-N homopeptide runs per se may be under selection variably in different lineages of budding yeasts, to control the evolution of functional N/Q-rich domains. Indeed, the relative selective burden on the protein sequences per se against harmful aggregation may vary as the potency of anti-aggregation cellular mechanisms varies. Assessment of these latter hypotheses would require experimental evolution investigations in tandem with novel theoretical developments. The evolution of mutation rates and the heterogeneity of rates for different types of mutation is a current area of interest in experimental evolution analysis [70, 71].
Availability of supporting data
Proteomes analyzed can be downloaded from the links listed in Additional file 1: Table S1. All other data sets supporting the results of this article are included within the article (and its additional files).
adenine + thymidine
guanidine + cytidine
intrinsically disordered protein
prionogenic amyloid former
whole genome duplication
Cox B. [PSI], a cytoplasmic suppressor of super-suppression in yeast. Heredity. 1965;20:505–21.
Shorter J, Lindquist S. Prions as adaptive conduits of memory and inheritance. Nat Rev Genets. 2005;6:435–50.
True H, Berlin I, Lindquist S. Epigenetic regulation of translation reveals hidden genetic variation to produce comlex traits. Nature. 2004;431:184–7.
True H, Lindquist S. A yeast prion provides a mechanism for genetic variation and phenotypic diversity. Nature. 2000;407:477–83.
Lacroute F. Non-Mendelian mutation allowing ureidosuccinic acid uptake in yeast. J Bacteriol. 1971;106:519–22.
Wickner R. [URE3] as an altered URE2 protein: evidence for a prion analog in Saccharomyces cerevisiae. Science. 1994;264:528–30.
Wickner R, Edskes H, Roberts B, Baxa U, Pierce M, Ross E, et al. Prions: proteins as genes and infectious entities. Genes Dev. 2004;18:470–85.
McGlinchey RP, Kryndushkin D, Wickner RB. Suicidal [PSI+] is a lethal yeast prion. Proc Natl Acad Sci U S A. 2011;108(13):5337–41.
Nakayashiki T, Kurtzman C, Edskes H, Wickner R. Yeast prions [URE3] and [PSI+] are diseases. PNAS. 2005;102:10575–80.
Holmes DL, Lancaster AK, Lindquist S, Halfmann R. Heritable remodeling of yeast multicellularity by an environmentally responsive prion. Cell. 2013;153(1):153–65.
Harbi D, Harrison PM: Classifying prion and prion-like phenomena. Prion 2014, 8(2):161-165.
Harbi D, Parthiban M, Gendoo DM, Ehsani S, Kumar M, Schmitt-Ulms G, et al. PrionHome: a database of prions and other sequences relevant to prion phenomena. PLoS One. 2012;7(2):e31785.
Michelitsch MD, Weissman JS. A census of glutamine/asparagine-rich regions: implications for their conserved function and the prediction of novel prions. Proc Natl Acad Sci U S A. 2000;97(22):11910–5.
Harrison LB, Yu Z, Stajich JE, Dietrich FS, Harrison PM. Evolution of budding yeast prion-determinant sequences across diverse fungi. J Mol Biol. 2007;368(1):273–82.
Harrison PM, Gerstein M. A method to assess compositional bias in biological sequences and its application to prion-like glutamine/asparagine-rich domains in eukaryotic proteomes. Genome Biol. 2003;4(6):R40.
Medina EM, Jones GW, Fitzpatrick DA. Reconstructing the fungal tree of life using phylogenomics and a preliminary investigation of the distribution of yeast prion-like proteins in the fungal kingdom. J Mol Evol. 2011;73(3-4):116–33.
Alberti S, Halfmann R, King O, Kapila A, Lindquist S. A systematic survey identifies prions and illuminates sequence features of prionogenic proteins. Cell. 2009;137(1):146–58.
Espinosa Angarica V, Ventura S, Sancho J. Discovering putative prion sequences in complete proteomes using probabilistic representations of Q/N-rich domains. BMC Genomics. 2013;14:316.
Lancaster AK, Nutter-Upham A, Lindquist S, King OD. PLAAC: a web and command-line application to identify proteins with prion-like amino acid composition. Bioinformatics. 2014;30(17):2501–2.
Zambrano R, Conchillo-Sole O, Iglesias V, Illa R, Rousseau F, Schymkowitz J, et al. PrionW: a server to identify proteins containing glutamine/asparagine rich prion-like domains and their amyloid cores. Nucleic Acids Res. 2015;43(W1):W331–7.
Ross E, Edskes H, Terry M, Wickner R. Primary sequence independence for prion formation. PNAS. 2005;102:12825–30.
Ross ED, Baxa U, Wickner RB. Scrambled prion domains form prions and amyloid. Mol Cell Biol. 2004;24(16):7206–13.
Toombs JA, Petri M, Paul KR, Kan GY, Ben-Hur A, Ross ED. De novo design of synthetic prion domains. Proc Natl Acad Sci U S A. 2012;109(17):6519–24.
Ross ED, Maclea KS, Anderson C, Ben-Hur A. A bioinformatics method for identifying Q/N-rich prion-like domains in proteins. Methods Mol Biol. 2013;1017:219–28.
Edskes HK, Wickner RB. The [URE3] prion in Candida. Eukaryot Cell. 2013;12(4):551–8.
Edskes HK, Khamar HJ, Winchester CL, Greenler AJ, Zhou A, McGlinchey RP, et al. Sporadic distribution of prion-forming ability of Sup35p from yeasts and fungi. Genetics. 2014;198(2):605–16.
Tanaka M, Chien P, Yonekura K, Weissman JS. Mechanism of cross-species prion transmission: an infectious conformation compatible with two highly divergent yeast prion proteins. Cell. 2005;121(1):49–62.
Kushnirov VV, Ter-Avanesyan MD, Didichenko SA, Smirnov VN, Chernoff YO, Derkach IL, et al. Divergence and conservation of SUP2 (SUP35) gene of yeast Pichia pinus and Saccharomyces cerevisiae. Yeast. 1990;6(6):461–72.
Talarek N, Maillet L, Cullin C, Aigle M. The [URE3] prion is not conserved among Saccharomyces species. Genetics. 2005;171(1):23–34.
Nakayashiki T, Ebihara K, Bannai H, Nakamura Y. Yeast [PSI+] “prions” that are crosstransmissible and susceptible beyond a species barrier through a quasi-prion state. Mol Cell. 2001;7(6):1121–30.
Si K, Lindquist S, Kandel ER. A neuronal isoform of the aplysia CPEB has prion-like properties. Cell. 2003;115(7):879–91.
Si K, Choi YB, White-Grindley E, Majumdar A, Kandel ER. Aplysia CPEB can form prion-like multimers in sensory neurons that contribute to long-term facilitation. Cell. 2010;140(3):421–35.
Paul KR, Hendrich CG, Waechter A, Harman MR, Ross ED. Generating new prions by targeted mutation or segment duplication. Proc Natl Acad Sci U S A. 2015;112(28):8584–9.
Gendoo DM, Harrison PM. The landscape of the prion protein’s structural response to mutation revealed by principal component analysis of multiple NMR ensembles. PLoS Comput Biol. 2012;8(8):e1002646.
Gendoo DM, Harrison PM. Discordant and chameleon sequences: their distribution and implications for amyloidogenicity. Protein Sci. 2011;20(3):567–79.
Sun Z, Diaz Z, Fang X, Hart MP, Chesi A, Shorter J, et al. Molecular determinants and genetic modifiers of aggregation and toxicity for the ALS disease protein FUS/TLS. PLoS Biol. 2011;9(4):e1000614.
Kim HJ, Kim NC, Wang YD, Scarborough EA, Moore J, Diaz Z, et al. Mutations in prion-like domains in hnRNPA2B1 and hnRNPA1 cause multisystem proteinopathy and ALS. Nature. 2013;495(7442):467–73.
Navarro S, Marinelli P, Diaz-Caballero M, Ventura S. The prion-like RNA-processing protein HNRPDL forms inherently toxic amyloid-like inclusion bodies in bacteria. Microb Cell Fact. 2015;14:102.
Suzuki G, Shimazu N, Tanaka M. A yeast prion, Mod5, promotes acquired drug resistance and cell survival under environmental stress. Science. 2012;336(6079):355–9.
Derkatch I, Bradley M, Zhou P, Chernoff Y, Liebman S. Prions affect the appearance of other prions: The story of [PIN+]. Cell. 2001;106:171–82.
Klucevsek KM, Braun MA, Arndt KM. The Paf1 complex subunit Rtf1 buffers cells against the toxic effects of [PSI+] and defects in Rkr1-dependent protein quality control in Saccharomyces cerevisiae. Genetics. 2012;191(4):1107–18.
Kurahashi H, Oishi K, Nakamura Y. A bipolar personality of yeast prion proteins. Prion. 2011;5(4):305–10.
Oishi K, Kurahashi H, Pack CG, Sako Y, Nakamura Y. A bipolar functionality of Q/N-rich proteins: Lsm4 amyloid causes clearance of yeast prions. MicrobiologyOpen. 2013;2(3):415–30.
Decker CJ, Teixeira D, Parker R. Edc3p and a glutamine/asparagine-rich domain of Lsm4p function in processing body assembly in Saccharomyces cerevisiae. J Cell Biol. 2007;179(3):437–49.
Ross CD, McCarty BR, Hamilton M, Ben-Hur A, Ross ED. A promiscuous prion: efficient induction of [URE3] prion formation by heterologous prion domains. Genetics. 2009;183(3):929–40.
Ader C, Frey S, Maas W, Schmidt HB, Gorlich D, Baldus M. Amyloid-like interactions within nucleoporin FG hydrogels. Proc Natl Acad Sci U S A. 2010;107(14):6281–5.
Harrison PM. Exhaustive assignment of compositional bias reveals universally prevalent biased regions: analysis of functional associations in human and Drosophila. BMC Bioinformatics. 2006;7:441.
Harbi D, Kumar M, Harrison PM. LPS-annotate: complete annotation of compositionally biased regions in the protein knowledgebase. Database (Oxford). 2011;2011:baq031.
Harbi D, Harrison PM. Interaction networks of prion, prionogenic and prion-like proteins in budding yeast, and their role in gene regulation. PLoS One. 2014;9(6):e100615.
Fox NK, Brenner SE, Chandonia JM. SCOPe: Structural Classification of Proteins: extended integrating SCOP and ASTRAL data and classification of new structures. Nucleic Acids Res. 2014;42:D304–9.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402.
Dosztanyi Z, Csizmok V, Tompa P, Simon I. IUPred: web server for the prediction of intrinsically unstructured regions of proteins based on estimated energy content. Bioinformatics. 2005;21(16):3433–4.
Halfmann R, Jarosz DF, Jones SK, Chang A, Lancaster AK, Lindquist S. Prions are a common mechanism for phenotypic inheritance in wild yeasts. Nature. 2012;482(7385):363–8.
Patel BK, Gavin-Smyth J, Liebman SW. The yeast global transcriptional co-repressor protein Cyc8 can propagate as a prion. Nat Cell Biol. 2009;11(3):344–9.
Shewmaker F, Mull L, Nakayashiki T, Masison DC, Wickner RB. Ure2p function is enhanced by its prion domain in Saccharomyces cerevisiae. Genetics. 2007;176(3):1557–65.
Wickner RB, Edskes HK, Shewmaker FP, Kryndushkin D, Nemecek J, McGlinchey R, et al. The relationship of prions and translation. Wiley Interdiscip Rev RNA. 2010;1(1):81–9.
Gilks N, Kedersha N, Ayodele M, Shen L, Stoecklin G, Dember LM, et al. Stress granule assembly is mediated by prion-like aggregation of TIA-1. Mol Biol Cell. 2004;15(12):5383–98.
Li X, Rayman JB, Kandel ER, Derkatch IL. Functional role of Tia1/Pub1 and Sup35 prion domains: directing protein synthesis machinery to the tubulin cytoskeleton. Mol Cell. 2014;55(2):305–18.
Halfmann R, Alberti S, Krishnan R, Lyle N, O’Donnell CW, King OD, et al. Opposing effects of glutamine and asparagine govern prion formation by intrinsically disordered proteins. Mol Cell. 2011;43(1):72–84.
Du Z, Crow ET, Kang HS, Li L. Distinct subregions of Swi1 manifest striking differences in prion transmission and SWI/SNF function. Mol Cell Biol. 2010;30(19):4644–55.
Du Z, Park KW, Yu H, Fan Q, Li L. Newly identified prion linked to the chromatin-remodeling factor Swi1 in Saccharomyces cerevisiae. Nat Genet. 2008;40(4):460–5.
Crow ET, Du Z, Li L. A small, glutamine-free domain propagates the [SWI(+)] prion in budding yeast. Mol Cell Biol. 2011;31(16):3436–44.
Rogoza T, Goginashvili A, Rodionova S, Ivanov M, Viktorovskaya O, Rubel A, et al. Non-Mendelian determinant [ISP+] in yeast is a nuclear-residing prion form of the global transcriptional regulator Sfp1. Proc Natl Acad Sci U S A. 2010;107(23):10573–7.
Mitchell AP, Magasanik B. Regulation of glutamine-repressible gene products by the GLN3 function in Saccharomyces cerevisiae. Mol Cell Biol. 1984;4(12):2758–66.
Myers LC, Gustafsson CM, Hayashibara KC, Brown PO, Kornberg RD. Mediator protein mutations that selectively abolish activated transcription. Proc Natl Acad Sci U S A. 1999;96(1):67–72.
Watanabe Y, Takaesu G, Hagiwara M, Irie K, Matsumoto K. Characterization of a serum response factor-like protein in Saccharomyces cerevisiae, Rlm1, which has transcriptional activity regulated by the Mpk1 (Slt2) mitogen-activated protein kinase pathway. Mol Cell Biol. 1997;17(5):2615–23.
Wendland J, Walther A. Genome evolution in the eremothecium clade of the Saccharomyces complex revealed by comparative genomics. G3. 2011;1(7):539–48.
Alexandrov AI, Ter-Avanesyan MD. Could yeast prion domains originate from polyQ/N tracts? Prion. 2013;7(3):209–14.
Malinovska L, Palm S, Gibson K, Verbavatz JM, Alberti S. Dictyostelium discoideum has a highly Q/N-rich proteome and shows an unusual resilience to protein aggregation. Proc Natl Acad Sci U S A. 2015;112(20):E2620–9.
Farlow A, Long H, Arnoux S, Sung W, Doak TG, Nordborg M, et al. The spontaneous mutation rate in the fission yeast schizosaccharomyces pombe. Genetics. 2015;201(2):737–44.
Ness RW, Morgan AD, Vasanthakrishnan RB, Colegrave N, Keightley PD. Extensive de novo mutation rate variation between individuals and across the genome of Chlamydomonas reinhardtii. Genome Res. 2015;25(11):1739–49.
This work was funded by the Natural Sciences and Engineering Research Council of Canada.
The authors declare they have no competing interests.
LA performed data analysis, made figures, and wrote some article text. DF performed data analysis. PH conceived the project, performed data analysis, made figures and wrote some article text. All authors read and approved the final manuscript.
List of fungal proteomes and their sources (weblinks). (XLS 37 kb)
Comparison of annotations of N/Q-rich proteins and prion predictions by the PAPA, PLAAC and PrionW programs. (DOCX 163 kb)
Large phylogenetic tree showing the trend in numbers of prion-like N/Q-rich proteins (NQPs). Colour-coding is according to a heatmap with green for low N/Q-rich numbers and red for high. The heatmap scale is indicated in the figure. The numbers of N-, Q-, N/Q- and Q/N-rich regions are listed for each species. Q/N-rich are regions that have a mingled bias of Qs and Ns, but mostly Q; similarly, for N/Q-rich. Clades are labelled where they branch off in the tree. (PNG 737 kb)
Large phylogenetic tree for N-rich/Q-rich ratio. Colour-coding is according to a heatmap with green for low N-rich/Q-rich ratio and red for high. The heatmap scale is indicated in the figure. Listed for each species are the total number of N-, Q-, N/Q- and Q/N-rich proteins and the N-rich/Q-rich ratio, which is the number of N-rich divided by the number of Q-rich proteins. Q/N-rich are regions that have a mingled bias of Qs and Ns, but mostly Q; similarly, for N/Q-rich. (PNG 737 kb)
Gene Ontology process category enrichments. (TXT 1 kb)
Other biases, that become prominent or depleted in Saccharomycetes. (DOCX 17 kb)
Large phylogenetic tree showing numbers of PLAAC, PAPA and combined PLAAC and PAPA (union of the two sets), and non-N/Q-rich PLAAC/PAPA predictions. Colour-coding is according to a heatmap with green for low percentage of non-N/Q-rich prion predictions and red for high. The heatmap scale is indicated in the figure. For this tree, for counting non-N/Q-rich PLAAC/PAPA predictions we use a strict threshold for N/Q bias (P = 1×10−5). (PNG 790 kb)
Further correlations referenced in the manuscript. (DOC 43 kb)