Next Article in Journal
Metabolic Reprogramming and Predominance of Solute Carrier Genes during Acquired Enzalutamide Resistance in Prostate Cancer
Next Article in Special Issue
Symmetry Breaking as an Interdisciplinary Concept Unifying Cell and Developmental Biology
Previous Article in Journal
Diazepam Promotes Translocation of Human Constitutive Androstane Receptor (CAR) via Direct Interaction with the Ligand-Binding Domain
Previous Article in Special Issue
Activation of Cdc42 GTPase upon CRY2-Induced Cortical Recruitment Is Antagonized by GAPs in Fission Yeast
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

The Path towards Predicting Evolution as Illustrated in Yeast Cell Polarity

Department of Bionanoscience, TU Delft, 2629 HZ Delft, The Netherlands
*
Author to whom correspondence should be addressed.
Cells 2020, 9(12), 2534; https://doi.org/10.3390/cells9122534
Submission received: 30 September 2020 / Revised: 18 November 2020 / Accepted: 21 November 2020 / Published: 24 November 2020
(This article belongs to the Special Issue Symmetry Breaking in Cells and Tissues)

Abstract

:
A bottom-up route towards predicting evolution relies on a deep understanding of the complex network that proteins form inside cells. In a rapidly expanding panorama of experimental possibilities, the most difficult question is how to conceptually approach the disentangling of such complex networks. These can exhibit varying degrees of hierarchy and modularity, which obfuscate certain protein functions that may prove pivotal for adaptation. Using the well-established polarity network in budding yeast as a case study, we first organize current literature to highlight protein entrenchments inside polarity. Following three examples, we see how alternating between experimental novelties and subsequent emerging design strategies can construct a layered understanding, potent enough to reveal evolutionary targets. We show that if you want to understand a cell’s evolutionary capacity, such as possible future evolutionary paths, seemingly unimportant proteins need to be mapped and studied. Finally, we generalize this research structure to be applicable to other systems of interest.

1. Introduction

How cells work and how they evolve is at the heart of cell biology. In this work we will review how cellular architecture (“how cells work”) and its evolutionary properties (“how they evolve”) are related to each other. Understanding evolution and possible mutational paths of protein networks, and especially the cell polarity network, is not only satisfying our curiosity but may also help us understand and possibly predict cancer progression [1].
Every cell consists of many different interconnected functional protein networks (for definitions, see Table 1), such as transcription, translation, or polarity establishment [2]. The network’s architecture, (for example: which protein binds to/reacts with which other protein), impacts the evolutionary possibilities of a network in multiple ways. For example, hubs, proteins with many binding partners, tend to evolve slower [3]. Less connected proteins, that may be deleted in a cell without a detectable change in cell physiology, can permit duplication of other genes and thus promote evolution [4]: duplicates of a gene enable new options for diversification, which facilitate further evolution of a gene/protein and the surrounding network [5,6,7]. Interestingly, many mutations (from 3% of nonsilent mutations in bacteria to 30% in hominids [8]) in a cell show very weak, or no effect on the cell’s function, a phenomenon called neutrality [9]. Thus, proteins that may seem unimportant for how the cell works now, in this environment, may become important when changes occur in the network architecture due to a mutation or a switch in environment [10,11].
A well-studied model organism to concretize how these proteins, without a detectable phenotype, shape a network is Saccharomyces cerevisiae, or budding yeast. The organism generally exhibits many of the network properties defined in Table 1, such as hierarchy and the presence of hubs [2]. Here, neutrality is also pervasive, as only 40% of homozygous gene deletions for the entire organism initially had obvious phenotypes [12]. Moreover, the environment has been shown to have a notable influence on neutrality, as lethal heterozygous deletions can be compensated by poor medium [13]. As a general rule, both the network architecture and the environment can mask the function of many proteins.
Within this organism, a, to some extent, representative example of a protein network is the polarity network, which governs how the yeast chooses a direction in which to divide and involves directing dozens of proteins in a process of breaking its internal spherical symmetry (see, e.g., [14]). As required, we observe the presence of the common network properties demonstrated in Section 2, such as hierarchy and redundancy. Polarity is also one of many biological functions in yeast for which a subdivision of many proteins into a quasi-modular network proved possible [15]. Within polarity, even more detailed submodules can be distinguished [16]. Neutrality is also exhibited by several polarity proteins discussed in Section 3, in part responsible for the difficulty in determining the role of each protein. Lastly, polarity is a pattern formation process where, by definition, spatiotemporal dependencies are important, and understanding evolution generally relies on understanding this type of dependencies [17].
However, the polarity network is not a prime example of the sort of networks with abundant transcriptional regulation. Other templates are better suited for learning about the evolution of gene regulation, such as interaction networks centered around transcriptional regulators (e.g., Mcm1 [18]), ribosomal regulation (e.g., [19]) the stress response (e.g., [20]) or metabolic response (e.g., GAL pathway [21]). The existence of established regulatory templates thus conveniently complements our focus on symmetry breaking during polarity as a model for the protein interaction network, which is also a topic for evolutionary studies.
Concretely, in [22] a mutant strongly defective in polarity establishment was experimentally evolved and found to recover remarkably reproducibly, e.g., the first rescuing mutation to sweep the population was always the same. Because of this exhibited tractability of the adaptations, network structures within the polarity network that facilitated evolution could be concretely interpreted in terms of redundancies [23]. In another approach to determine the flexibility of the polarity network, historical evolution was studied for 40+ proteins in almost 300 fungal species in [11]. Again, the polarity network exhibited sufficient modularity so that studying its evolution separately from other functions still yielded interpretable results. For example, authors showed that polarity network size was shaped in part by fungal lifestyle (e.g., uni- or multicellular).
Nevertheless, clear justification for observed evolutionary trajectories in this network remain difficult to make. It is sometimes possible to generate more abstract predictions by linking network architecture to the evolvability of the polarity network using classical regulatory motifs. For example, the presence of positive and negative feedback in polarity establishment confers robustness to the network [24,25], which in turn may facilitate evolution. But in order to make more concrete and detailed predictions about evolutionary trajectories, an important insight is that we need mechanistic information of (parts of) the polarity network.
To illustrate this point, a bottom-up model was constructed in [26] where molecular details were coarse-grained following analysis on numerical simulations of multiple polarity mechanisms [23]. This approach was successful in quantitatively describing the fitness along the evolutionary trajectory in [22]. Furthermore, the predictions on epistasis, an important bottleneck for predicting evolution [27], can be extended to other modules, and although use of full mechanistic understanding is superior in quantitative assessments of epistasis, biofunctional information (viz. from GO-terms, in agreement with [28]) as input to the model suffices for epistasis sign predictions.
Instructive is the Nrp1 case where full information is absent, but some phenomenological information is available. Based on the latter, inclusion of Nrp1 into the bottom-up coarse-grained model of [26] is still worthwhile, but leaves room for improvement. This marks the importance of continuing to investigate protein networks until molecular mechanisms have been elucidated. In this review we advocate that obtaining this knowledge is (soon) feasible, motivating the use of the yeast polarity system for studies in network organization (with properties such as hierarchy, modularity/connectivity and redundancy that can “hide” proteins) as well as evolution.
We present three case studies in our review to reflect different stages in this knowledge quest; while literature on the first case, bud scar proteins, has been quite advanced for several years, only recently the GAP mechanism [23] (second case) has been revealed and for Nrp1 (third case) more work remains. These cases illustrate how to move from detecting the strong and more obvious phenotypes to unveiling evolutionary important hidden features (such as for Nrp1). The cases build on literature of the polarity network summarized in the form of a Venn diagram, ideal for depicting a hierarchical and semimodular protein grouping. The feasibility of deciphering a protein’s role in this Venn diagram turns out to depend on how deep the protein is embedded inside the network, which can impose neutrality on mutations to hinder unique attribution of genotypes to phenotypes. In the sections thereafter, our three examples with varying depths of embedding serve to show how improved understanding of the nontrivial parts of the networks can elucidate evolutionary trajectories. Ultimately, we believe that in all the work done in yeast for many decades by many researchers there is a common recipe applicable and useful for many protein networks, as expanded upon in the outlook.

2. Polarity Overview

Within the yeast polarity network, four pathways to polarization exist which cannot easily be considered modular. Their interconnectivity can be conveniently visualized in the form of a Venn diagram (Figure 1, with references found in Appendix Table A1). These pathways are hierarchically set up, in the following order: the mating pathway at the top, then the bud scar pathway, followed by the reaction−diffusion pathway and finally the actin pathway. In short, their function boils down to the act of condensing the GTPase Cdc42 bound to GTP molecules (i.e., active Cdc42) to one point on the plasma membrane, which can signal downstream effectors to proceed the cell cycle [29]. To prevent premature or overdue localization of active Cdc42 and allow some influence on the hierarchy of pathways, a fifth pathway exists to control the previous four, namely the timing pathway. The next section summarizes the most important interactions in and across all pathways, starting with the timing cue, before expanding upon the three examples. As a site note, we have done our best to include all relevant papers, but apologize for important papers we have missed.

2.1. Timing: The Control Knob

During isotropic growth in G1, active Cdc42 localization is suppressed by overactivity of its associated GTPase activating proteins (GAPs) and sequestration of its guanine nucleotide exchange factor (GEF). The consequence of both circumstances is the vast abundance of inactive Cdc42, which is bound to GDP instead of GTP [30], rendering it impossible to signal the polarity cue. The purpose of this pathway (see top dark purple region in Figure 1) is hence to timely reduce GAP activity and release the GEF, which must be in response to important physiological parameters that indicate the readiness of the cell: sufficient protein production, a sufficient size, and sufficient nutrition.
The physiological state of the cell enters the equation through nuclear levels of cyclin Cln3. Upon sufficient nutrition and size, Cln3 levels rise either more directly through higher Cln3 mRNA abundance [31], or more indirectly through Ydj1 disturbing Cln3 localization by Whi3 [32], the latter also being an inhibitor of Cln3 mRNA translation [33]. The arrival of nuclear Cln3 allows binding partner and cyclin-dependent kinase Cdc28 [34] to phosphorylate Whi5, which had inhibited expression of Cln2, another cyclin [35,36]. Cln2 can then reinforce its own expression, consolidating the original Cln3 signal [37].
Now, the Cdc28-Cln2 complex can distribute the physiological signal to the aforementioned targets, the GAPs and the GEF. The kinase Cdc28 phosphorylates all four GAPs Bem2, Bem3, Rga1 and Rga2 [38,39,40,41] and Far1, which was keeping the GEF Cdc24 in the nucleus [42]. Now cytoplasmic levels of active Cdc42 can rise, leading to polarity establishment through subsequent pathways.
Importantly, the completion of the timing pathway causes the hierarchy of the subsequent pathways to change. While the mating pathway is otherwise dominant, the kinase Cdc28 phosphorylates Ste5, a crucial hub in the mating pathway, to stop the mating in its tracks [43,44,45]. In the following discussion of the mating pathway, the situation is considered where the timing pathway did not overwrite its behavior.

2.2. Mating: Heavily Cross-Linked

The mating pathway is the dominant force across the four symmetry-breaking pathways. While polarization in a random orientation is possible after the timing cue (see the section on reaction−diffusion further on), the presence of pheromones of the opposite mating type (a or α) should redirect the Cdc42 localization to the side of the pheromone signal. This process revolves around Ste5, as also depicted in the left, blue-grey circle of Figure 1.
Briefly put, once pheromones bind membrane proteins Ste2 and Ste3 [46,47], Ste4 is released from the membrane [48,49] and binds Ste20 and scaffold Ste5. This scaffold binds Ste7, Ste11 and Fus3, which are activated by sequential phosphorylation [50,51,52]. Fus3 may inhibit the GAPs Bem2 and Bem3 [39], while Ste5 binds the GEF Cdc24 [53], replacing the absence of the timing pathway result to stimulate activity of Cdc42.
While this simplified view would suffice to redirect the Cdc42 localization, the mating pathway is much more intertwined with the other pathways than seemingly necessary, particularly with the actin pathway. The abundant mechanistic redundancies result in a more complex picture, obfuscating the role of the proteins involved. For example, active Cdc42 stimulates Ste11 phosphorylation/activation [54]. Another form of positive feedback, as well as a bridge to the actin pathway, is the Cdc42 recruitment of formin Bni1 [55]. The resulting nucleation of actin cables may transport Ste5-GEF Cdc24 complexes [56], possibly also through Bem1. This scaffold co-immunoprecipitates with Act1 [57] and Far1 [58], which is bound to the GEF Cdc24 [59], but is itself in turn also bound to Bem1 [60]. Another actin cross-link is the phosphorylation and localization of Bni1 through Fus3 [61]. Clearly, care must be taken in assigning roles to different proteins, as many are overloaded.

2.3. Bud Scar: Mostly Modular and Ordered

In the absence of a mating cue, the timing pathway reduces GAP activity and releases the GEF, while the mating pathway is repressed. Under the new hierarchy, the bud scar pathway is normally dominant. The scar refers to leftover proteins from the previous division, named septins [62,63]. This spatial cue can be exploited for polarity establishment; a new bud forms adjacent to the scar (axial budding, haploids) or also at the opposing side (bipolar budding, diploids) [64,65]. The bottom, dark blue circle of Figure 1 represents this path from septins to Cdc42 recruitment graphically. More background information about the core bud scar protein group Bud1 to Bud5 is discussed separately in one of the three case studies, and only a brief overview of the pathway as a whole is discussed in this section.
An important bud scar localization target is Bud5, which activates and recruits Bud1 [66,67], a protein also known as Rsr1 [65,68]. In haploids, where Axl1 is specifically expressed [69], the Bud5 localization occurs by relay of the septin signal to a protein complex of Bud3, Bud4, Axl1, Axl2 and Bud5 [70,71,72,73]. In diploids, functionality of Rax1 and Rax2 is not impaired, presumably by blocked expression of Axl1 [74], so these can localize Bud9 and Bud8 to the bud scar or the opposite end of the scar respectively [75], which in turn recruit Bud5 [76].
After Bud5, localization follows of Bud2 [77], the GAP for Bud1 [78], to complete the control of the GTPase cycle of Bud1. Finally, Bud1 binds GEF Cdc24 and Bem1 [79] (although Cdc24 has the strongest affinity with Bem1 [80]), to redirect the pattern formation made possible after the timing cue. As linkage of Cdc42 GAP Rga1 to septins prevents reuse of the previous location [81], the new bud forms adjacent to the bud scar.
As a whole, the bud scar pathway is not completely modular either. Aside from nudging the reaction−diffusion pathway (see next section), an example of a cross-link is that Bud8 and Bud9 are delivered by actin transport [82]. The highest position in the hierarchy in the absence of the mating cue is also not absolute; multiple ways to promote the subsequent reaction−diffusion pathway exist, such as deletion of Bud1 and Bud8 [83], Axl2 and Rax1 [84], or Bem1 [22]. Therefore, it has been possible to retrieve the information discussed in the following section.

2.4. Reaction−Diffusion: Ample Redundancy

Even in the absence of chemical or spatial cues, the shift in balance towards activation of Cdc42 induced by the timing pathway still provides the conditions for swift symmetry breaking. Theoretical models concerning this pathway have been subject to updates as more molecular details have been revealed. Central is the strong positive feedback generated by the Bem1-GEF Cdc24 complex, as modelled in, e.g., [85,86] with further refinement in [87]. This feedback is sufficient for polarity success, which becomes rather insensitive to GAP abundance. More details on the GAPs were uncovered in [23] and are placed in a broader context in the case study further on.
What makes this pathway special is the limited number of proteins that are unique to this pathway, as seen from the central, emerald circle in Figure 1, namely only Cla4 and Rdi1. The latter is the least cross-linked of the two, providing a possible justification for referring to the WT mechanism as the Rdi1 polarity mechanism, as in [24]. Cla4 is more context-dependent, possibly having two opposing roles, promoting and inhibiting polarity [23,88,89]. Yet both Rdi1 and Cla4 are dispensable for polarity [88,90].
An even stronger addition to the redundancy within this pathway is on the positive feedback side. Without Bem1, generic rescuing feedbacks suffice [23], among which Cla4 could account for 20% of their function [91]. More feedbacks may be found in the GAPs (see GAP case study) through actin transport as described in, e.g., [24]. This brings us to the actin pathway as the final layer to discuss.

2.5. Actin: The Mysterious Auxiliary Layer

The actin pathway (rightmost green circle in Figure 1) has featured several times already in the previously discussed pathways, but its individual role is still quite uncertain. Yeast formin Bni1 which nucleates actin cables, binds active Cdc42 [55], and is known to be involved in exo- and endocytosis [92,93]. This suggests transport of polarity proteins from and to the presumptive bud site. The resulting actin pathway has been confusingly implicated in two opposing roles; promoting Cdc42 polarization, see, e.g., [24,94,95], as well as negatively impacting Cdc42 polarization [96,97].
A way to reconcile these findings is that actin transport contributes to a process promoting Cdc42 polarization but without relying on significant transport of Cdc42 itself. As mentioned in the mating pathway, Bem1 and Act1 co-immunoprecipitate [57], suggesting that Bem1, and concordantly its multiple binding partners, might get transported through the actin pathway. However, in the absence of Bem1, 80% of the positive feedback is still unidentified [91], which may very well be actin-related. Instead, a prime candidate is the GAP group, which is known to bind the epsin-coating of actin cables involved in endocytosis [98]. This is further discussed in the case study on the GAPs. In any case, it is quite difficult to decipher the actin pathway, in large part due to its low positioning in the yeast polarity hierarchy.

3. Case Studies

In the introduction the need for determining the (potential) functions of a protein inside a complex network even when these are normally hidden, was explained as, evolution may exploit these later. In the following sections, the reconstitution of the mechanistic details involving three protein(s) (classes), that are currently at different stages in their discovery (see graphical mechanistic summary in Figure 2), is illustrated; the Bud1-5 bud scar protein subset, the GAP proteins for Cdc42, and finally Nrp1.
In the first case, detailed knowledge has been established the longest, and has also been put in an evolutionary perspective as forming a “weak regulatory link” that facilitates phenotypic variation [99,100]. Therefore, relatively small gains are expected from the remaining open questions. For the GAP proteins, detailed knowledge is only very recently obtained. Together, the purpose of putting the literature of these cases into historical context is that these two examples form convenient templates to delineate the route towards complete understanding for proteins as in the final case, Nrp1, where knowledge is scarce. On the one hand, technical advances that improve detection of phenotypes were the most prominent method of progress for the bud scar proteins. On the other hand, the GAP protein and the Nrp1 case show that progress can rely more on very specific designs of experiments than on new technology. All three cases combined in turn illustrate the foundation under the general strategy to approach networks put forward in the outlook.

3.1. Bud Scar Proteins

One of the strongest phenotypes to observe for budding yeast is the location of the next bud. Already more than half a century ago, it was documented that S. cerevisiae exhibited two possible budding patterns [64]. Normally for haploids, the next bud grows next to the previous division site (axial budding), whereas diploids can also pick the location opposite the previous site (bipolar budding). Changes in this pattern are clear phenotypes, and are therefore a useful detection tool in bulk mutagenesis screens. In 1991, five genes were identified that affected (seemingly as their sole phenotype) the bud site selection, and were therefore named BUD1 (formerly known as RSR1 [68]), BUD2, BUD3, BUD4 and BUD5 [65,66]. These genes and their associated proteins are therefore straightforwardly localized in the bud scar circle in the Venn diagram (Figure 1).
The next step, retrieving information on physical interactions, is more elaborate. For example, authors in [66] reasoned the physical mechanism of Bud5 (amongst others a GEF for Bud1) based on sequence information, which was further substantiated in [78], together with the GAP action of Bud2. Bud2 and Bud5 hence underlie the relevant GTPase cycling of Bud1, as in vitro studies established that Bud1 when GTP-bound, mostly binds Cdc24, while binding and recruiting Bem1 when GDP-bound [79,101]. This provides the bridge to the previously described reaction−diffusion pathway. With the advent of GFP [102], it also became possible to get in vivo spatiotemporal information from yeast proteins, later even in bulk [103]. Now armed with the tool of fluorescence microscopy, substantial progress was made for validating the roles for Bud2 and Bud5 [77].
The final confirmation on the mechanisms comes when zooming into the domains of a protein, by expressing truncated versions or intelligently chosen point mutations. For example, the domains responsible for how Bud5 affects whether budding is axial or bipolar are described in [104]. For other proteins, however, it was necessary to consider in greater detail the protein network that the bud genes compose. Therefore, it has taken more time to reverse-engineer the details for Bud3 and Bud4. Authors in [72] combine their findings with earlier literature to demonstrate that septins from the previous division localize Bud4, which can then bind Bud3. Subsequently, Bud4 binds Axl2, and after this Axl1 as well. Domain analysis of Bud3 further uncovered more details [105] to show Bud3 can serve as a GEF for Cdc42 (but early in G1, before Start), which is important for the Bud4−Axl1 link.
This leaves only a few mechanistic unknowns. For example, the exact recruitment of Bud2 as depicted in Figure 2 is putative, as it is unclear how Bud2 is recruited even in the absence of Bud1 or Bud5 [74]. Bud1 can recruit itself by dimerization and Cdc42 at the bud site of the membrane [67]. In the same paper, it is shown that while Bud5 also recruits Bud1, the former may only indirectly be involved in the Bud1 dimerization, considering the full GTP-GDP cycle is critical for appropriate dimerization. This suggests guidance of Bud2 for the self-recruitment, but the exact order for Bud1-binding and hydrolysis of its GTP is open for interpretation, given the alternative model in [74].
As an alternative representation of the current state of research of the bud scar proteins, Figure 3 depicts an ordering in the form of a mind map. For clarity, only Bud1 is considered. As can be seen, the bud scar protein case is relatively well advanced, with ample coverage of all categories.

3.2. GAPs

Central to the function of GTPase Cdc42 are proteins that promote GTP hydrolysis, i.e., its GTPase activating proteins. Bem3, was the first to be identified [106], followed by Rga1 and Rga2 [107,108,109]. A summary of Bem3 studies/information is given in Figure 4. The specific molecular function of GAPs allowed validation in vitro, where Cdc42-GTP was incubated with a GAP to determine whether the amount of GTP indeed decayed faster than without a GAP. In [106] Bem2 was not found to exhibit detectable GAP activity, and not until [110] could Bem2 be convincingly considered a GAP as well.
After their establishment as a GAP, the localization of Bem3, Rga1 and Rga2 was determined using either antibody staining or GFP-tagging [111]. These GAPs colocalized with Cdc42 at the bud site (although Rga1 is slightly more dispersed). Bem3 was also found in Spitzenkörper-like structures, but seemingly after polarity establishment during polarized growth [112]. The GAP localizations affirmed their role in polarity, although they also seemed somewhat redundant, as a triple mutant was still viable [109].
More information on GAPs was gathered through their interactions. Bem2 and Bem3 were found to be associated with the mating pathway [39], Bem3, Rga1 and Rga2 related to the actin pathway [98,113], while [81] suggested a link between Rga1 and the bud scar pathway. However, several open questions remain. For example, why is there more than one GAP? Why and how are these distributed across multiple pathways?
More clues were to follow from the tedious deciphering of the mechanistic action of the GAPs. The difficulty behind this problem was clearly elucidated by the model of [87]. There, authors showed the dominance of the Bem1-mediated positive feedback minimized the phenotypical influence of varying GAP concentration. Consequently, this led to the realization that GAP details only emerged in a Δbem1 background.
In [23] a proposed mechanistic model for the GAPs was validated against this background. The idea was that GAPs were temporarily retained on the membrane by Cdc42 during the GTP hydrolysis process. Only in the location in the cell with high membrane concentrations of Cdc42-GTP would this lead to a local depletion of available GAPs, which cannot be compensated by the cytosolic diffusive flux. Elsewhere on the membrane, this was not a problem and Cdc42 was promptly inactivated and recycled, leading to only one spot on the plasma membrane where active Cdc42 accumulated. A generic positive feedback mechanism for Cdc42 (due to, e.g., Cla4) completed the symmetry breaking. This provides a good example of the added value of revealing nonobvious mechanistic details. This GAP model, which is obscured in the presence of Bem1, allowed the authors in [23] to provide a detailed explanation of the evolutionary trajectory observed in [22].
Yet, this may not be the complete role of the GAPs. Their mechanism partly supplements the need for the Bem1-mediated positive feedback, but a further deletion of CLA4 and BEM3 reveals that 80% of the positive feedback originates elsewhere [91] (p. 101). Here the interconnectivity with another pathway can surface, which initially seemed elusive and redundant.
As discussed in the actin pathway, there could be a critical link with actin and the GAPs. The aforementioned local depletion of unbound GAPs may be reinforced by actin transport [91] (p. 34). Epsin coatings of endocytic vesicles colocalize with polarized growth and bind GAPs [98,113]. Moreover, active Cdc42 releases the auto-inhibition of kinases Ste20 [54] and Cla4 [114], both of which phosphorylate myosins 3 and 5 to ultimately lead to activation of the Arp2/3 complex [115], critical for endocytosis. In this way, sites of active Cdc42 can promote the endocytosis which may reinforce the stability of the site, providing the feedback needed to establish polarity (model F in [116]). The combination of interactions of recyclable GAPs and active Cdc42 would also fulfil the requirement of actin-mediated recruitment formulated in [117], provided the GAPs diffuse slowly on the membrane [96].
Figure 4 graphically summarizes the available information for one of the GAPs, Bem3. However, there is still room for improvement with experiments whose design can just now be established. As with the bud scar pathway, subtle information may be retrieved through experiments at the domain level, to test, e.g., the GAP trafficking hypothesis. One could remove the link between GAPs and actin through deletion of epsins ENT1 and ENT2 (in the Δbem1 Δcla4 background) and replacing this by only a weakly expressed ENTH-domain of Ent1 (truncation). Modulating this expression should show how strong this effect is. More information on GAP interactors in this role may also be retrieved by using this mutant as the crippled starting point in an evolution experiment, akin to [22].
Furthermore, the resulting scatter of the GAPs across the other polarity pathways currently leaves room for interpretation and speculation. Given Figure 1, the components that are most shared also seem the most critical. This is most obvious when noting that actin pathway components are not just essential for polarity establishment. For example, Rho1 is needed for cell wall integrity synthesis later on during polarized growth [118,119,120]. Extrapolating, the location (wedged between pathways) of the GAPs in the Venn diagram may suggest an important (but not essential) function for each of them in establishing polarity.
Hypothetically, the GAPs might serve as an evolutionary control knob to mediate the relative hierarchy between pathways. This could be favorable in situations where different hierarchies are optimal, such as when mating is infrequent (e.g., the diploid state becomes the default), or when the bud scar is not often used (frequent sporulation). Strategic dispersal of multiple GAPs may therefore provide more handles for the cell to optimize the pathways then simply having one GAP in larger copy numbers.

3.3. Nrp1

After discussing two well-studied protein classes, we address an underexposed protein, namely Nrp1. With this we would like to show the difficulties and possibilities that are still open in a case when the most straightforward experiments do not provide obvious, interpretable phenotypes. Although its deletion does not have a detectable phenotype in standard lab conditions, Nrp1 is an important evolutionary [22,121]. Usually, Nrp1 is mentioned merely peripherally in articles as a bycatch in studies with an alternative focus. Therefore, a chronologically ordered literature overview does not make sense here, as very little of the research findings actually builds on previous work. An overview of the Nrp1 knowledge is given in Figure 5, where it is apparent that there are some gaps in our understanding.
Nrp1 was first described by [122], who have given the protein its name. NRP1 stands for ‘Asparagine rich protein’, the name refers to the region of the protein sequence that has many asparagines (short name: “N”). Genes are often named for their defining characteristics or functions. Reynaud and coworkers [122] did not find a phenotype or function for Nrp1, thus the seemingly nondescript name.
Nrp1 has been linked to stress response and stress granules. For example, it has been implicated in the response to glucose and oxygen [123,124,125], although the precise mechanism or function of Nrp1 in this response is not known. A hypothesis is that Nrp1 forms an aggregate or prion (like a stress granule) by its low complexity domains, because the repeated asparagine sequence in Nrp1 is often found to form prions for other proteins [126]. Nrp1 itself also seems to form a prion. In addition, Nrp1 can potentially bind and regulate mRNA. The mRNA regulation occurs through another documented domain, namely an RNA binding motif [126]. This implies that Nrp1 can bind specific mRNA sequences, however no specific mRNAs have been identified until now [127].
The supposed link between Nrp1 and the polarity network was found by authors in [22]. They showed that null mutations in NRP1 could rescue a bem1Δ. This prompted interest in a search for the function of Nrp1. Another connection to the polarity network and a possible explanation for what was found in [22] is the synthetic lethality of NRP1 with CLA4 [128].
Diepeveen et al. [11] have found that Nrp1 is highly conserved within the Ascomycota which hints to a function of some importance. One typically expects that essential genes are the most conserved parts as they cannot easily be mutated [3]. However this is not always the case, for example CDC42 is conserved in most fungal species (and also outside of fungi [129,130]) but it is not present in others [11]. Interestingly, NRP1 is more conserved than several essential genes [11]. What the reason for this conservation is, is unclear. A conserved sequence does not mean conserved function or interactions. So, although a highly similar Nrp1 protein is present in a species, this does not mean it has the same function in that species.
An interesting look into the kind of experiments and research done for a protein that is similar in sequence to Nrp1 and also probably functionally related to Nrp1 is given by Whi3 [131]. Whi3 like Nrp1 has an RNA recognition motif and a low complexity/repeated region. One difference is that Whi3 has repeated glutamines and Nrp1 has repeated asparagines, though these amino acids are chemically similar. In the paper of [131] the authors removed these domains and checked functionality and localization of Whi3 in Ashbya gossypii. Whi3 in A.gossypii and budding yeast share important functionality [33,131]. Whi3 localizes to stress granules, indirectly regulates the G1/S phase transition via Cln3 and affects many mRNAs in yeast [33]. A similar role may be hypothesized for Nrp1. It also localizes to stress granules [124], affects the G1 exit [22], and may bind some mRNAs [127]. Different from Whi3, Nrp1 does not have a clear RNA target (like Cln3) that explains its functioning.
Here we will provide some research paths from different areas of research to test our hypothesis from the previous paragraph. First, from a more chemistry perspective one strategy is to look at the structure of Nrp1. The order of the domains is known, but it is unclear whether the low-complexity domain will fold, as they have been shown to be disordered [132]. This unstructured part of the protein may move about freely. It would be interesting to see if the C-terminal part of the protein does fold specifically again after the unstructured part. One would be able to visualize the structure of the protein by means of NMR [133], however it might be difficult in this case to determine the structure if it is indeed unstructured/moving.
Second, looking at Nrp1 from a cell biological perspective, deletion studies are often used. Unfortunately, this does not work as easily for NRP1 as the single deletion does not yield any phenotype. However, in a different background a phenotype can be found, the bem1Δ [22]. Analyzing what happens in these cells will help understand Nrp1. This can be done on a population level by doing a fitness assay. On a single cell level one can use microscopy. A previous high-throughput study shows Nrp1 present in the cytoplasm [103]. A more detailed study focusing on polarity establishment can give more insight into the functioning of Nrp1, especially when combined with the previous approach (different genetic backgrounds).
By the same token, a more in-depth analysis of the RNA binding ability of Nrp1 is also relevant. Again, using Nrp1 in different environments may give different results and find different specific RNAs that are bound. RNA chip-seq has become easier to execute over the last few years [134] and thus now it might be possible to do the proposed experiments.
The ultimate goal is to find a molecular mechanism for Nrp1, but this goal cannot be reached without knowledge of other aspects of the network. Nrp1 is a good example of a protein that is buried deep in the network, which makes the investigation challenging. However, it is worthwhile to dig deep and find the hidden functionality of proteins like Nrp1 that seem neutral at first glance, but have a significant evolutionary role.

4. Outlook

Protein−protein interaction networks have diverse properties that influence its evolution, such as redundancy, hierarchy and neutrality. We have advocated that studying yeast polarity provides a suitable starting point for general studies on the evolution of these complex networks, as it exhibits many of the aforementioned traits and includes the evolutionary relevant spatiotemporal dynamics [17]. From work in [26], it had become clear that predicting epistasis, which precedes predicting evolution, required deep understanding of the interactions and reactions constituting a protein network like yeast polarity.
However, our case studies demonstrate that obtaining this knowledge in a complex network architecture is far from trivial due to the aforementioned redundancy, hierarchy and neutrality. While some genes generate obvious phenotypes from which the roles of the corresponding gene products are easily deduced, for others the required information is only revealed after multiple rounds of precisely designed experiments. The latter situation can be due to various reasons related to the network architecture, for example the associated protein is found deep down in the hierarchy, forms part of a redundancy in the network or is currently otherwise peripheral to the core function. Regardless of the origin, encountered neutrality is a complication, which is a particularly unresolved feature in the third case study (Nrp1).
By exploring the past and present of the polarity network, we have aimed to determine how to advance for future research, to answer the open questions in the field and underline research opportunities. Although the yeast polarity network has been studied for many years, there still remain many unsolved mysteries. For example, the study by [22] gave much insight into the evolution and possible back-up mechanisms of the polarity network, but gave rise to the question of how Nrp1 is involved.
Much research has been done under perfect lab conditions, which results in specific results. It would be interesting to also design experiments that explore the genotype−phenotype map in different conditions. This would make it possible to find previously hidden components, that do not show up under standard lab conditions. Another possibility is changing the environment together with the expression level of specific genes, which can affect fitness [135], as has been shown for Cdc42 [23].
In vitro work is also an important next step to isolate parts of the network and see if these parts can independently perform a function [136]. As an example in budding yeast, in vitro reconstitution of a She-protein mediated mechanism of asymmetric mRNA transport revealed subtle details that were hard to demonstrate in vivo [137].
From an evolution standpoint it is interesting to see how the network as it is in current strains came to be and what the variation is that can be found in the wild. The variation within a wild population shows the spread that is available in the genotype map. It provides insight into what genotypes are preferable in certain environments. Apart from the genotype showing the history of a population, also, for example, the expression levels of proteins can be inherited [138]. Such epigenetic inheritance can be important for the reaction of an organism to stress and other environmental factors [139].
The yeast polarity examples discussed also delineate a general route forward in dismantling complexly connected protein networks. A graphical overview is depicted in Figure 6. The arrow heads indicate the type of information obtained from proteins inside the network, and higher degrees of information are successively more difficult to obtain, and usually rely on first reaching the previous level. In this way, we work our way deeper into the protein network and slowly but steadily elucidate beyond the obvious phenotypes.
In the top category, much work has been done determining the effects of simple deletions. For budding yeast, a large knock-out database has been present for more than two decades [12]. The ease of experimentally parallelizing the deletion construction even allows for ample double deletion data constituting genetic interactions, which for multiple model systems is bundled in the BioGRID database [140]. Yet, there is a myriad of combinatorial possibilities for gene deletions, and it has become clear from the GAP and Nrp1 examples that these need to be explored intelligently, rather than by brute force. Well-designed starting points for evolution experiments, for example to find that GAPs and Nrp1 genetically interact with Bem1 in [22], or SATAY assays [141] can elucidate interactions of genes of interest with important domains that only surface with strong phenotypes in the right genetic background.
However, genetic interactions can be very indirect. Epistasis is known to act globally even on unrelated networks during adaptation [142], so gathering physical information is a welcome next step. The efficient two-hybrid screens date back more than three decades [143], where two proteins fused to a DNA binding domain and transcriptional activator, respectively, promote transcription of a reporter gene when interacting. To follow interactions across the cell (see also next paragraph), FRET imaging, where two fluorescent fusion proteins cause the emission spectrum to shift when in close proximity, has also been extensively used [144], but this method also relies on the proteins of interest to be tolerant to protein fusions. More subtle modifications for tagging are used for co-immunoprecipitation [145], and are still heavily used in budding yeast [146]. Once physical interactors have been established, these can be further confirmed with their relevant binding sites by point mutation to influence the binding, as in, e.g., [84] for Bem1. This level of precision also means a move towards low-throughput data gathering, but can be useful to conjecture how cells after the deletion of BEM1 evolved [23].
While the physical interactions constitute a rudimentary form of the protein network, the next step in understanding originates from adding spatiotemporal information. For example, Bud3 is a GEF for Cdc42, but only during early G1 phase [105], before the timing pathway gives the cue for symmetry breaking. If the function of interest is symmetry breaking, this can be excluded from the network overview as in Figure 1. Generally, a high-throughput manner for establishing protein localization in vivo has been with fluorescent protein fusions, particularly with GFP variants [103]. As aforementioned, immunostaining provides a similar option for tagging and hence localization, but requires fixation of the cells, making the temporally transient contributions of components more difficult to trace. If temporal rather than spatial information on the importance of a protein is of the essence, ingenious solutions exist that conditionally disable the protein of interest, after which the results can be swiftly observed. Examples include degron systems [147,148] and optogenetic tools [149]. Finally, as with Nrp1 localization of mRNAs may be the important function, options to trace these are also present with in situ hybridization, as described in, e.g., [150].
Ultimately, when all previous steps have been performed, it becomes possible to make the next step towards complete understanding, which would be mechanistic understanding. As shown in the case of actin, if insufficient information is available, modelling can lead to uncertain and contradictory results, such as is the case with the role of actin in polarity establishment [24,95,96]. If possible, the most unambiguous results would come from bottom-up approaches, such as modelling Michaelis−Menten kinetics for metabolism (e.g., [151] (pp. 165–180)) or solving reaction−diffusion systems, as for polarity done in [23,87] in case proteins cannot be assumed to be uniformly distributed. The more complete understanding of the protein network is then put to the test by the prediction of observed adaptive trajectories in historical or experimental evolution, as done, for example, in the latter case for the path of [22] in [23,26]. Considering the many options and solutions evolution can explore, accurate prediction of evolution is the best guarantee that the full extent of the knowledge of a network has been explored.
In summary, the general path to full network understanding as outlined in Figure 6 brings us from genetic and physical interactions to visualizing precise protein dynamics, modeling and full reconstitution. While initially, even the genetic interaction map was a tedious chore, high-throughput studies and bioinformatics tools continue to facilitate the gathering of information. Considering the speed with which the technological advances occur, the necessary data for network understanding becomes feasible for many more functions and organisms. This marks the relevance of establishing a generalizable and efficient workflow to obtain the right network data and use it for understanding and predicting its evolution.

Author Contributions

Conceptualization, W.K.-G.D., E.S. and L.L.; investigation, W.K.-G.D. and E.S.; writing—original draft preparation, W.K.-G.D. and E.S.; writing—review and editing, W.K.-G.D., E.S. and L.L.; visualization, W.K.-G.D. and E.S.; supervision, L.L. All authors have read and agreed to the published version of the manuscript.

Funding

L.L. and W.K-G.D. gratefully acknowledge support from the Netherlands Organization for Scientific Research (NWO/OCW), as part of the Gravitation Program: Frontiers of Nanoscience. L.L. and E.S. gratefully acknowledge funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (Grant agreement No. [758132]).

Acknowledgments

We would like to thank Enzo Kingma and Leila Iñigo de la Cruz for careful reading of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

Table A1 shows the corresponding literature references for the information in the Venn diagram in Figure 1. In addition, gene essentiality was obtained from [152].
Table A1. Literature supporting information (e.g., bonds and reactions) underlying Figure 1.
Table A1. Literature supporting information (e.g., bonds and reactions) underlying Figure 1.
DescriptionReference
Nrp1 binds Pub1[153]
Nrp1-Pub1 promotes Whi3-Cln3 bond (putative)
Pub1 binds Whi3[154]
Nrp1-Pub1 inhibits Whi3-Cla4 (mRNA) binding (putative)
Whi3 binds and inhibits Cln3[33,155]
Whi3 binds Cdc28[34]
Cln3 binds Cdc28[34]
Cdc28-Cln3 inhibits activation of Whi5[36]
Whi5 binds and inhibits activation of Swi4/Swi6/Mbp1[36]
Cdc28-Cln2 inhibits activation of Whi5 as well[36,44]
Swi4/Swi6/Mbp1 promotes activity of Cln2[35]
Whi3 binds Cla4 (mRNA)[156]
Whi3 can aggregate in stress granules[33]
Cdc28 binds Cln2[157]
Cdc28-Cln2 inhibits Ste5[43]
Cdc28-Cln2 inhibits Far1-Cdc24 binding[42]
Cdc28-Cln2 inhibits Bem2[39]
Cdc28-Cln2 inhibits Rga1[38]
Cdc28-Cln2 inhibits Bem3[39]
Cdc28-Cln2 inhibits Rga2[40,41]
Far1 inhibits Cln2[44,45]
Swi4/Swi6/Mbp1 promotes activity of Bud9[82]
Ydj1 inhibits Whi3[32]
Protein synthesis inhibits Ydj1[32]
Size promotes Ydj1[32]
Nutrition promotes Cln3[31]
Pheromone activates Ste2/Ste3[46,47]
Ste2/Ste3 binds Ste4-Ste18[48]
Ste4-Ste18 binds Ste20[158]
Ste4-Ste18 binds Cdc24-Far1[59]
Ste20 activates Ste11[51]
Ste5 binds Ste11[50]
Ste11 activates Ste7[51]
Ste5 binds Ste7[50]
Ste7 activates Fus3[52]
Ste5 binds Fus3[50]
Ste5 binds Bem1[57]
Fus3 inhibits Bem3 activity (putative)[39]
Fus3 inhibits Bem2 activity (putative)[39]
Far1 binds Cdc24[59]
Far1 binds and promotes Cdc24 activity[159]
Cdc42 binds and promotes Ste20 activity[54]
Bem1 binds Far1[58]
Bem1 binds Ste20[57]
Fus3 binds and promotes Bni1[61]
Ste5 binds the exocytosis network[56]
Ste4 binds Ste5[51]
Ste5 promotes Ste20 activating Ste11[51]
Ste5 promotes Ste11 activating Ste7[51]
Ste5 promotes Ste7 activating Fus3[51,52]
Septins bind Bud4[73]
Bud3 binds Bud4[72]
Bud3-Bud4 binds Axl1[72]
Bud3-Bud4 binds Axl2[72]
Axl2 binds Bud5[72]
Septins bind Rga1 (putative)[81]
Septins bind Bud9[82]
Axl1 inhibits Rax2 (putative)[74]
Axl1 inhibits Rax1 (putative)[74]
Rax1 binds Bud8[75,76]
Rax1 binds Bud9[75,76]
Bud8 binds exocytosis network[82]
Bud2 binds and inhibits Bud1[78]
Bud5 binds and activates Bud1[78,160]
Bud1 binds Cdc24[80]
Bud1 binds Bem1[79,101]
Bud9 binds exocytosis network[82]
Bud1 binds Cdc42[67,161]
Bud3 binds Cdc42[105]
Bud5 bind Bud8[76]
Bud5 binds Bud9[76]
Rdi1 binds Cdc42[162]
Rga2 binds and deactivates Cdc42[109]
Rga2 binds Bem1[40]
Rga1 binds and deactivates Cdc42[109]
Cdc42 binds Bem1[163]
Cla4 binds Bem1[164]
Bem1 binds Cdc24, activates it[60,84,164]
Cla4 inhibits binding Bem1-Cdc24[89,163,165]
Cdc24 binds and activates Cdc42[106,164]
Bem2 binds and deactivates Cdc42[110]
Cla4 binds and activates Cdc42[88]
Bem3 binds and deactivates Cdc42[106]
Cla4 inhibits Cdc42-Rdi1 bond[88]
Cdc24 binds Bem1-Cla4[164]
Bni1 binds Act1[166]
Act1 self-organizes[167]
Bnr1 binds Act1[166]
Exocytosis network binds Cdc42[94]
Act1 binds Arp2/Arp3[115]
Arp2/Arp3 interacts with intermediate[168,169]
Intermediates interact with Ent1/Ent2[168,169]
Ent2 inhibits Bem3 deactivating Cdc42[113]
Ent1/Ent2 inhibits Rga2 deactivating Cdc42[98]
Cdc42-Cla4 activates Arp2/3[115]
Cdc42-Ste20 activates Arp2/3[115]
Ent1/Ent2 inhibits Rga1 deactivating Cdc42[98]
Bni1 binds Cdc42[55]
Act1 binds Bem1[57]
Bud6 binds Bni1[93]
Rho1 binds Bni1[93]
Spa2 binds Bni1[93]

References

  1. Greaves, M.; Maley, C.C. Clonal evolution in cancer. Nature 2012, 481, 306–313. [Google Scholar] [CrossRef] [PubMed]
  2. Costanzo, M.; VanderSluis, B.; Koch, E.N.; Baryshnikova, A.; Pons, C.; Tan, G.; Wang, W.; Usaj, M.; Hanchard, J.; Lee, S.D. A global genetic interaction network maps a wiring diagram of cellular function. Science 2016, 353. [Google Scholar] [CrossRef] [PubMed]
  3. Fraser, H.B.; Wall, D.P.; Hirsh, A.E. A simple dependence between protein evolution rate and the number of protein-protein interactions. BMC Evol. Biol. 2003, 3, 1–6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Ratmann, O.; Wiuf, C.; Pinney, J.W. From evidence to inference: Probing the evolution of protein interaction networks. HFSP J. 2009, 3, 290–306. [Google Scholar] [CrossRef] [PubMed]
  5. Zhang, J. Evolution by gene duplication: An update. Trends Ecol. Evol. 2003, 18, 292–298. [Google Scholar] [CrossRef] [Green Version]
  6. Näsvall, J.; Sun, L.; Roth, J.R.; Andersson, D.I. Real-time evolution of new genes by innovation, amplification, and divergence. Science 2012, 338, 384–387. [Google Scholar] [CrossRef] [Green Version]
  7. Magadum, S.; Banerjee, U.; Murugan, P.; Gangapur, D.; Ravikesavan, R. Gene duplication as a major force in evolution. J. Genet. 2013, 92, 155–161. [Google Scholar] [CrossRef]
  8. Eyre-Walker, A.; Keightley, P.D. The distribution of fitness effects of new mutations. Nat. Rev. Genet. 2007, 8, 610–618. [Google Scholar] [CrossRef]
  9. Thoday, J.M. Non-Darwinian “evolution” and biological progress. Nature 1975, 255, 675–677. [Google Scholar] [CrossRef]
  10. Khan, A.I.; Dinh, D.M.; Schneider, D.; Lenski, R.E.; Cooper, T.F. Negative epistasis between beneficial mutations in an evolving bacterial population. Science 2011, 332, 1193–1196. [Google Scholar] [CrossRef]
  11. Diepeveen, E.T.; Gehrmann, T.; Pourquié, V.; Abeel, T.; Laan, L. Patterns of Conservation and Diversification in the Fungal Polarization Network. Genome Biol. Evol. 2018, 10, 1765–1782. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Winzeler, E.A. Functional Characterization of the S. cerevisiae Genome by Gene Deletion and Parallel Analysis. Science 1999, 285, 901–906. [Google Scholar] [CrossRef] [Green Version]
  13. Deutschbauer, A.M.; Jaramillo, D.F.; Proctor, M.; Kumm, J.; Hillenmeyer, M.E.; Davis, R.W.; Nislow, C.; Giaever, G. Mechanisms of Haploinsufficiency Revealed by Genome-Wide Profiling in Yeast. Genetics 2005, 169, 1915–1925. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Martin, S.G.; Arkowitz, R.A. Cell polarization in budding and fission yeasts. FEMS Microbiol. Rev. 2014, 38, 228–253. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Costanzo, M.; Baryshnikova, A.; Bellay, J.; Kim, Y.; Spear, E.D.; Sevier, C.S.; Ding, H.; Koh, J.L.; Toufighi, K.; Mostafavi, S.; et al. The genetic landscape of a cell. Science 2010, 327, 425–431. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Gao, J.T.; Guimera, R.; Li, H.; Pinto, I.M.; Sales-Pardo, M.; Wai, S.C.; Rubinstein, B.; Li, R. Modular coherence of protein dynamics in yeast cell polarity system. Proc. Natl. Acad. Sci. USA 2011, 108, 7647–7652. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Yamada, T.; Bork, P. Evolution of biomolecular networks—Lessons from metabolic and protein interactions. Nat. Rev. Mol. Cell Biol. 2009, 10, 791–803. [Google Scholar] [CrossRef] [PubMed]
  18. Tuch, B.B.; Galgoczy, D.J.; Hernday, A.D.; Li, H.; Johnson, A.D. The evolution of combinatorial gene regulation in fungi. PLoS Biol. 2008, 6, e38. [Google Scholar] [CrossRef]
  19. Tanay, A.; Regev, A.; Shamir, R. Conservation and evolvability in regulatory networks: The evolution of ribosomal regulation in yeast. Proc. Natl. Acad. Sci. USA 2005, 102, 7203–7208. [Google Scholar] [CrossRef] [Green Version]
  20. Dhar, R.; Sagesser, R.; Weikert, C.; Wagner, A. Yeast Adapts to a Changing Stressful Environment by Evolving Cross-Protection and Anticipatory Gene Regulation. Mol. Biol. Evol. 2013, 30, 573–588. [Google Scholar] [CrossRef] [Green Version]
  21. Slot, J.C.; Rokas, A. Multiple GAL pathway gene clusters evolved independently and by different mechanisms in fungi. Proc. Natl. Acad. Sci. USA 2010, 107, 10136–10141. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Laan, L.; Koschwanez, J.H.; Murray, A.W. Evolutionary adaptation after crippling cell polarization follows reproducible trajectories. eLife 2015, 4, e09638. [Google Scholar] [CrossRef] [PubMed]
  23. Brauns, F.; De la Cruz, L.M.I.; Daalman, W.K.-G.; De Bruin, I.; Halatek, J.; Laan, L.; Frey, E. Adaptability and evolution of the cell polarization machinery in budding yeast. bioRxiv 2020. [Google Scholar] [CrossRef]
  24. Freisinger, T.; Klünder, B.; Johnson, J.; Müller, N.; Pichler, G.; Beck, G.; Costanzo, M.; Boone, C.; Cerione, R.A.; Frey, E.; et al. Establishment of a robust single axis of cell polarity by coupling multiple positive feedback loops. Nat. Commun. 2013, 4, 1807. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Howell, A.S.; Jin, M.; Wu, C.-F.; Zyla, T.R.; Elston, T.C.; Lew, D.J. Negative Feedback Enhances Robustness in the Yeast Polarity Establishment Circuit. Cell 2012, 149, 322–333. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Daalman, W.K.-G.; Laan, L. Predicting an epistasis-rich genotype-phenotype map with a coarse-grained bottom-up model of budding yeast polarity. bioRxiv 2020. [Google Scholar] [CrossRef]
  27. Breen, M.S.; Kemena, C.; Vlasov, P.K.; Notredame, C.; Kondrashov, F.A. Epistasis as the primary factor in molecular evolution. Nature 2012, 490, 535–538. [Google Scholar] [CrossRef]
  28. Yu, M.K.; Kramer, M.; Dutkowski, J.; Srivas, R.; Licon, K.; Kreisberg, J.F.; Ng, C.T.; Krogan, N.; Sharan, R.; Ideker, T. Translation of Genotype to Phenotype by a Hierarchy of Cell Subsystems. Cell Syst. 2016, 2, 77–88. [Google Scholar] [CrossRef] [Green Version]
  29. Park, H.-O.; Bi, E. Central Roles of Small GTPases in the Development of Cell Polarity in Yeast and Beyond. Microbiol. Mol. Biol. Rev. 2007, 71, 48–96. [Google Scholar] [CrossRef] [Green Version]
  30. Hall, A. The cellular functions of small GTP-binding proteins. Science 1990, 249, 635–640. [Google Scholar] [CrossRef]
  31. Gallego, C.; Garí, E.; Colomina, N.; Herrero, E.; Aldea, M. The Cln3 cyclin is down-regulated by translational repression and degradation during the G1 arrest caused by nitrogen deprivation in budding yeast. EMBO J. 1997, 16, 7196–7206. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Vergés, E.; Colomina, N.; Garí, E.; Gallego, C.; Aldea, M. Cyclin Cln3 Is Retained at the ER and Released by the J Chaperone Ydj1 in Late G1 to Trigger Cell Cycle Entry. Mol. Cell 2007, 26, 649–662. [Google Scholar] [CrossRef]
  33. Cai, Y.; Futcher, B. Effects of the Yeast RNA-Binding Protein Whi3 on the Half-Life and Abundance of CLN3 mRNA and Other Targets. PLoS ONE 2013, 8, e84630. [Google Scholar] [CrossRef] [Green Version]
  34. Wang, H.; Garí, E.; Vergés, E.; Gallego, C.; Aldea, M. Recruitment of Cdc28 by Whi3 restricts nuclear accumulation of the G1 cyclin–Cdk complex to late G1. EMBO J. 2004, 23, 180–190. [Google Scholar] [CrossRef] [Green Version]
  35. Koch, C.; Schleiffer, A.; Ammerer, G.; Nasmyth, K. Switching transcription on and off during the yeast cell cycle: Cln/Cdc28 kinases activate bound transcription factor SBF (Swi4/Swi6) at start, whereas Clb/Cdc28 kinases displace it from the promoter in G2. Genes Dev. 1996, 10, 129–141. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Costanzo, M.; Nishikawa, J.L.; Tang, X.; Millman, J.S.; Schub, O.; Breitkreuz, K.; Dewar, D.; Rupes, I.; Andrews, B.; Tyers, M. CDK activity antagonizes Whi5, an inhibitor of G1/S transcription in yeast. Cell 2004, 117, 899–913. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Stuart, D.; Wittenberg, C. CLN3, not positive feedback, determines the timing of CLN2 transcription in cycling cells. Genes Dev. 1995, 9, 2780–2794. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Archambault, V.; Chang, E.J.; Drapkin, B.J.; Cross, F.R.; Chait, B.T.; Rout, M.P. Targeted proteomic study of the cyclin-Cdk module. Mol. Cell 2004, 14, 699–711. [Google Scholar] [CrossRef]
  39. Knaus, M.; Pelli-Gulli, M.-P.; Van Drogen, F.; Springer, S.; Jaquenoud, M.; Peter, M. Phosphorylation of Bem2p and Bem3p may contribute to local activation of Cdc42p at bud emergence. EMBO J. 2007, 26, 4501–4513. [Google Scholar] [CrossRef] [Green Version]
  40. McCusker, D.; Denison, C.; Anderson, S.; Egelhofer, T.A.; Yates, J.R.; Gygi, S.P.; Kellogg, D.R. Cdk1 coordinates cell-surface growth with the cell cycle. Nat. Cell Biol. 2007, 9, 506–515. [Google Scholar] [CrossRef]
  41. Sopko, R.; Huang, D.; Smith, J.C.; Figeys, D.; Andrews, B.J. Activation of the Cdc42p GTPase by cyclin-dependent protein kinases in budding yeast. EMBO J. 2007, 26, 4487–4500. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Shimada, Y.; Gulli, M.-P.; Peter, M. Nuclear sequestration of the exchange factor Cdc24 by Far1 regulates cell polarity during yeast mating. Nat. Cell Biol. 2000, 2, 117–124. [Google Scholar] [CrossRef] [PubMed]
  43. Garrenton, L.S.; Braunwarth, A.; Irniger, S.; Hurt, E.; Kunzler, M.; Thorner, J. Nucleus-Specific and Cell Cycle-Regulated Degradation of Mitogen-Activated Protein Kinase Scaffold Protein Ste5 Contributes to the Control of Signaling Competence. Mol. Cell. Biol. 2009, 29, 582–601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Doncic, A.; Falleur-Fettig, M.; Skotheim, J.M. Distinct Interactions Select and Maintain a Specific Cell Fate. Mol. Cell 2011, 43, 528–539. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Doncic, A.; Skotheim, J.M. Feedforward Regulation Ensures Stability and Rapid Reversibility of a Cellular State. Mol. Cell 2013, 50, 856–868. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Blumer, K.J.; Reneke, J.E.; Thorner, J. The STE2 gene product is the ligand-binding component of the alpha-factor receptor of Saccharomyces cerevisiae. J. Biol. Chem. 1988, 263, 10836–10842. [Google Scholar] [PubMed]
  47. Hagen, D.C.; McCaffrey, G.; Sprague, G.F. Evidence the yeast STE3 gene encodes a receptor for the peptide pheromone a factor: Gene sequence and implications for the structure of the presumed receptor. Proc. Natl. Acad. Sci. USA 1986, 83, 1418–1422. [Google Scholar] [CrossRef] [Green Version]
  48. Whiteway, M.; Hougan, L.; Dignard, D.; Thomas, D.Y.; Bell, L.; Saari, G.C.; Grant, F.J.; O’Hara, P.; MacKay, V.L. The STE4 and STE18 genes of yeast encode potential β and γ subunits of the mating factor receptor-coupled G protein. Cell 1989, 56, 467–477. [Google Scholar] [CrossRef]
  49. Hirschman, J.E.; De Zutter, G.S.; Simonds, W.F.; Jenness, D.D. The Gβγ Complex of the Yeast Pheromone Response Pathway SUBCELLULAR FRACTIONATION AND PROTEIN-PROTEIN INTERACTIONS. J. Biol. Chem. 1997, 272, 240–248. [Google Scholar] [CrossRef] [Green Version]
  50. Chol, K.-Y.; Satterberg, B.; Lyons, D.M.; Elion, E.A. Ste5 tethers multiple protein kinases in the MAP kinase cascade required for mating in S. cerevisiae. Cell 1994, 78, 499–512. [Google Scholar] [CrossRef]
  51. Feng, Y.; Song, L.Y.; Kincaid, E.; Mahanty, S.K.; Elion, E.A. Functional binding between Gβ and the LIM domain of Ste5 is required to activate the MEKK Ste11. Curr. Biol. 1998, 8, 267–282. [Google Scholar] [CrossRef] [Green Version]
  52. Good, M.; Tang, G.; Singleton, J.; Reményi, A.; Lim, W.A. The Ste5 Scaffold Directs Mating Signaling by Catalytically Unlocking the Fus3 MAP Kinase for Activation. Cell 2009, 136, 1085–1097. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Wang, Y.; Chen, W.; Simpson, D.M.; Elion, E.A. Cdc24 Regulates Nuclear Shuttling and Recruitment of the Ste5 Scaffold to a Heterotrimeric G Protein in Saccharomyces cerevisiae. J. Biol. Chem. 2005, 280, 13084–13096. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Lamson, R.E.; Winters, M.J.; Pryciak, P.M. Cdc42 Regulation of Kinase Activity and Signaling by the Yeast p21-Activated Kinase Ste20. Mol. Cell. Biol. 2002, 22, 2939–2951. [Google Scholar] [CrossRef] [Green Version]
  55. Evangelista, M.; Blundell, K.; Longtine, M.S.; Chow, C.J.; Adames, N.; Pringle, J.R.; Peter, M.; Boone, C. Bni1p, a yeast formin linking cdc42p and the actin cytoskeleton during polarized morphogenesis. Science 1997, 276, 118–122. [Google Scholar] [CrossRef]
  56. Qi, M.; Elion, E.A. Formin-induced actin cables are required for polarized recruitment of the Ste5 scaffold and high level activation of MAPK Fus3. J. Cell Sci. 2005, 118, 2837–2848. [Google Scholar] [CrossRef] [Green Version]
  57. Leeuw, T.; Fourest-Lieuvin, A.; Wu, C.; Chenevert, J.; Clark, K.; Whiteway, M.; Thomas, D.Y.; Leberer, E. Pheromone Response in Yeast: Association of Bem1p with Proteins of the MAP Kinase Cascade and Actin. Science 1995, 1210–1213. [Google Scholar] [CrossRef]
  58. Lyons, D.M.; Mahanty, S.K.; Choi, K.-Y.; Manandhar, M.; Elion, E.A. The SH3-domain protein Bem1 coordinates mitogen-activated protein kinase cascade activation with cell cycle control in Saccharomyces cerevisiae. Mol. Cell. Biol. 1996, 16, 4095–4106. [Google Scholar] [CrossRef] [Green Version]
  59. Butty, A.-C.; Pryciak, P.M.; Huang, L.S.; Herskowitz, I.; Peter, M. The role of Far1p in linking the heterotrimeric G protein to polarity establishment proteins during yeast mating. Science 1998, 282, 1511–1516. [Google Scholar] [CrossRef]
  60. Peterson, J. Interactions between the bud emergence proteins Bem1p and Bem2p and Rho- type GTPases in yeast. J. Cell Biol. 1994, 127, 1395–1406. [Google Scholar] [CrossRef] [Green Version]
  61. Matheos, D.; Metodiev, M.; Muller, E.; Stone, D.; Rose, M.D. Pheromone-induced polarization is dependent on the Fus3p MAPK acting through the formin Bni1p. J. Cell Biol. 2004, 165, 99–109. [Google Scholar] [CrossRef] [Green Version]
  62. Kim, H.B. Cellular morphogenesis in the Saccharomyces cerevisiae cell cycle: Localization of the CDC3 gene product and the timing of events at the budding site. J. Cell Biol. 1991, 112, 535–544. [Google Scholar] [CrossRef]
  63. Flescher, E.G. Components required for cytokinesis are important for bud site selection in yeast. J. Cell Biol. 1993, 122, 373–386. [Google Scholar] [CrossRef] [Green Version]
  64. Freifelder, D. Bud position in Saccharomyces cerevisiae. J. Bacteriol. 1960, 80, 567. [Google Scholar] [CrossRef] [Green Version]
  65. Chant, J.; Herskowitz, I. Genetic control of bud site selection in yeast by a set of gene products that constitute a morphogenetic pathway. Cell 1991, 65, 1203–1212. [Google Scholar] [CrossRef]
  66. Chant, J.; Corrado, K.; Pringle, J.R.; Herskowitz, I. Yeast BUD5, encoding a putative GDP-GTP exchange factor, is necessary for bud site selection and interacts with bud formation gene BEM1. Cell 1991, 65, 1213–1224. [Google Scholar] [CrossRef] [Green Version]
  67. Kang, P.J.; Béven, L.; Hariharan, S.; Park, H.-O. The Rsr1/Bud1 GTPase interacts with itself and the Cdc42 GTPase during bud-site selection and polarity establishment in budding yeast. Mol. Biol. Cell 2010, 21, 3007–3016. [Google Scholar] [CrossRef] [Green Version]
  68. Bender, A.; Pringle, J.R. Multicopy suppression of the cdc24 budding defect in yeast by CDC42 and three newly identified genes including the ras-related gene RSR1. Proc. Natl. Acad. Sci. USA 1989, 86, 9976–9980. [Google Scholar] [CrossRef] [Green Version]
  69. Fujita, A.; Oka, C.; Arikawa, Y.; Katagai, T.; Tonouchi, A.; Kuhara, S.; Misumi, Y. A yeast gene necessary for bud-site selection encodes a protein similar to insulin-degrading enzymes. Nature 1994, 372, 567. [Google Scholar] [CrossRef]
  70. Chant, J.; Mischke, M.; Mitchell, E.; Herskowitz, I.; Pringle, J.R. Role of Bud3p in producing the axial budding pattern of yeast. J. Cell Biol. 1995, 129, 767–778. [Google Scholar] [CrossRef] [Green Version]
  71. Sanders, S.L.; Herskowitz, I. The BUD4 protein of yeast, required for axial budding, is localized to the mother/BUD neck in a cell cycle-dependent manner. J. Cell Biol. 1996, 134, 413–427. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Kang, P.J.; Angerman, E.; Jung, C.-H.; Park, H.-O. Bud4 mediates the cell-type-specific assembly of the axial landmark in budding yeast. J. Cell Sci. 2012, 125, 3840–3849. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Kang, P.J.; Hood-DeGrenier, J.K.; Park, H.-O. Coupling of septins to the axial landmark by Bud4 in budding yeast. J. Cell Sci. 2013, 126, 1218–1226. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Bi, E.; Park, H.-O. Cell Polarization and Cytokinesis in Budding Yeast. Genetics 2012, 191, 347–387. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Kang, P.J.; Angerman, E.; Nakashima, K.; Pringle, J.R.; Park, H.-O. Interactions among Rax1p, Rax2p, Bud8p, and Bud9p in marking cortical sites for bipolar bud-site selection in yeast. Mol. Biol. Cell 2004, 15, 5145–5157. [Google Scholar] [CrossRef]
  76. Krappmann, A.-B.; Taheri, N.; Heinrich, M.; Mösch, H.-U. Distinct domains of yeast cortical tag proteins Bud8p and Bud9p confer polar localization and functionality. Mol. Biol. Cell 2007, 18, 3323–3339. [Google Scholar] [CrossRef] [Green Version]
  77. Marston, A.L.; Chen, T.; Yang, M.C.; Belhumeur, P.; Chant, J. A localized GTPase exchange factor, Bud5, determines the orientation of division axes in yeast. Curr. Biol. 2001, 11, 803–807. [Google Scholar] [CrossRef] [Green Version]
  78. Bender, A. Genetic evidence for the roles of the bud-site-selection genes BUD5 and BUD2 in control of the Rsr1p (Bud1p) GTPase in yeast. Proc. Natl. Acad. Sci. USA 1993, 90, 9926–9929. [Google Scholar] [CrossRef] [Green Version]
  79. Park, H.-O.; Bi, E.; Pringle, J.R.; Herskowitz, I. Two active states of the Ras-related Bud1/Rsr1 protein bind to different effectors to determine yeast cell polarity. Proc. Natl. Acad. Sci. USA 1997, 94, 4463–4468. [Google Scholar] [CrossRef] [Green Version]
  80. Zheng, Y.; Bender, A.; Cerione, R.A. Interactions among proteins involved in bud-site selection and bud-site assembly in Saccharomyces cerevisiae. J. Biol. Chem. 1995, 270, 626–630. [Google Scholar] [CrossRef] [Green Version]
  81. Tong, Z.; Gao, X.-D.; Howell, A.S.; Bose, I.; Lew, D.J.; Bi, E. Adjacent positioning of cellular structures enabled by a Cdc42 GTPase-activating protein mediated zone of inhibition. J. Cell Biol. 2007, 179, 1375–1384. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Schenkman, L.R.; Caruso, C.; Pagé, N.; Pringle, J.R. The role of cell cycle–regulated expression in the localization of spatial landmark proteins in yeast. J. Cell Biol. 2002, 156, 829–841. [Google Scholar] [CrossRef] [PubMed]
  83. Irazoqui, J.E.; Gladfelter, A.S.; Lew, D.J. Scaffold-mediated symmetry breaking by Cdc42p. Nat. Cell Biol. 2003, 5, 1062–1070. [Google Scholar] [CrossRef] [PubMed]
  84. Smith, S.E.; Rubinstein, B.; Mendes Pinto, I.; Slaughter, B.D.; Unruh, J.R.; Li, R. Independence of symmetry breaking on Bem1-mediated autocatalytic activation of Cdc42. J. Cell Biol. 2013, 202, 1091–1106. [Google Scholar] [CrossRef] [Green Version]
  85. Ozbudak, E.M.; Becskei, A.; Van Oudenaarden, A. A System of Counteracting Feedback Loops Regulates Cdc42p Activity during Spontaneous Cell Polarization. Dev. Cell 2005, 9, 565–571. [Google Scholar] [CrossRef] [Green Version]
  86. Goryachev, A.B.; Pokhilko, A.V. Dynamics of Cdc42 network embodies a Turing-type mechanism of yeast cell polarity. FEBS Lett. 2008, 582, 1437–1443. [Google Scholar] [CrossRef] [Green Version]
  87. Klünder, B.; Freisinger, T.; Wedlich-Söldner, R.; Frey, E. GDI-Mediated Cell Polarization in Yeast Provides Precise Spatial and Temporal Control of Cdc42 Signaling. PLoS Comput. Biol. 2013, 9, e1003396. [Google Scholar] [CrossRef] [Green Version]
  88. Tiedje, C.; Sakwa, I.; Just, U.; Höfken, T. The rho gdi rdi1 regulates rho gtpases by distinct mechanisms. Mol. Biol. Cell 2008, 19, 2885–2896. [Google Scholar] [CrossRef] [Green Version]
  89. Kuo, C.-C.; Savage, N.S.; Chen, H.; Wu, C.-F.; Zyla, T.R.; Lew, D.J. Inhibitory GEF Phosphorylation Provides Negative Feedback in the Yeast Polarity Circuit. Curr. Biol. 2014, 24, 753–759. [Google Scholar] [CrossRef] [Green Version]
  90. Cvrckova, F.; De Virgilio, C.; Manser, E.; Pringle, J.R.; Nasmyth, K. Ste20-like protein kinases are required for normal localization of cell growth and for cytokinesis in budding yeast. Genes Dev. 1995, 9, 1817–1830. [Google Scholar] [CrossRef] [Green Version]
  91. Daalman, W.K.-G. Coupling of Genotype-Phenotype Maps to Noise-Driven Adaptation, Showcased in Yeast Polarity; TU Delft: Delft, The Netherlands, 2020. [Google Scholar] [CrossRef]
  92. Pruyne, D.; Gao, L.; Bi, E.; Bretscher, A. Stable and dynamic axes of polarity use distinct formin isoforms in budding yeast. Mol. Biol. Cell 2004, 15, 4971–4989. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  93. Prosser, D.C.; Drivas, T.G.; Maldonado-Báez, L.; Wendland, B. Existence of a novel clathrin-independent endocytic pathway in yeast that depends on Rho1 and formin. J. Cell Biol. 2011, 195, 657–671. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  94. Wedlich-Soldner, R.; Altschuler, S.; Wu, L.; Li, R. Spontaneous cell polarization through actomyosin-based delivery of the Cdc42 GTPase. Science 2003, 299, 1231–1235. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  95. Slaughter, B.D.; Unruh, J.R.; Das, A.; Smith, S.E.; Rubinstein, B.; Li, R. Non-uniform membrane diffusion enables steady-state cell polarization via vesicular trafficking. Nat. Commun. 2013, 4, 1380. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Layton, A.T.; Savage, N.S.; Howell, A.S.; Carroll, S.Y.; Drubin, D.G.; Lew, D.J. Modeling Vesicle Traffic Reveals Unexpected Consequences for Cdc42p-Mediated Polarity Establishment. Curr. Biol. 2011, 21, 184–194. [Google Scholar] [CrossRef] [Green Version]
  97. Watson, L.J.; Rossi, G.; Brennwald, P. Quantitative Analysis of Membrane Trafficking in Regulation of Cdc42 Polarity. Traffic 2014, 15, 1330–1343. [Google Scholar] [CrossRef] [Green Version]
  98. Aguilar, R.C.; Longhi, S.A.; Shaw, J.D.; Yeh, L.-Y.; Kim, S.; Schön, A.; Freire, E.; Hsu, A.; McCormick, W.K.; Watson, H.A.; et al. Epsin N-terminal homology domains perform an essential function regulating Cdc42 through binding Cdc42 GTPase-activating proteins. Proc. Natl. Acad. Sci. USA 2006, 103, 4116–4121. [Google Scholar] [CrossRef] [Green Version]
  99. Gerhart, J.; Kirschner, M. The theory of facilitated variation. Proc. Natl. Acad. Sci. USA 2007, 104, 8582–8589. [Google Scholar] [CrossRef] [Green Version]
  100. Wedlich-Soldner, R.; Li, R. Yeast and fungal morphogenesis from an evolutionary perspective. Semin. Cell Dev. Biol. 2008, 19, 224–233. [Google Scholar] [CrossRef] [Green Version]
  101. Miller, K.E.; Lo, W.-C.; Chou, C.-S.; Park, H.-O. Temporal regulation of cell polarity via the interaction of the Ras GTPase Rsr1 and the scaffold protein Bem1. Mol. Biol. Cell 2019, 30, 2543–2557. [Google Scholar] [CrossRef]
  102. Niedenthal, R.K.; Riles, L.; Johnston, M.; Hegemann, J.H. Green fluorescent protein as a marker for gene expression and subcellular localization in budding yeast. Yeast 1996, 12, 773–786. [Google Scholar] [CrossRef]
  103. Huh, W.-K.; Falvo, J.V.; Gerke, L.C.; Carroll, A.S.; Howson, R.W.; Weissman, J.S.; O’Shea, E.K. Global analysis of protein localization in budding yeast. Nature 2003, 425, 686–691. [Google Scholar] [CrossRef] [PubMed]
  104. Kang, P.J.; Lee, B.; Park, H.-O. Specific Residues of the GDP/GTP Exchange Factor Bud5p Are Involved in Establishment of the Cell Type-specific Budding Pattern in Yeast. J. Biol. Chem. 2004, 279, 27980–27985. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  105. Kang, P.J.; Lee, M.E.; Park, H.-O. Bud3 activates Cdc42 to establish a proper growth site in budding yeast. J. Cell Biol. 2014, 206, 19–28. [Google Scholar] [CrossRef] [Green Version]
  106. Zheng, Y.; Cerione, R.; Bender, A. Control of the yeast bud-site assembly GTPase Cdc42. Catalysis of guanine nucleotide exchange by Cdc24 and stimulation of GTPase activity by Bem3. J. Biol. Chem. 1994, 269, 2369–2372. [Google Scholar]
  107. Stevenson, B.J.; Ferguson, B.; De Virgilio, C.; Bi, E.; Pringle, J.R.; Ammerer, G.; Sprague, G.F. Mutation of RGA1, which encodes a putative GTPase-activating protein for the polarity-establishment protein Cdc42p, activates the pheromone-response pathway in the yeast Saccharomyces cerevisiae. Genes Dev. 1995, 9, 2949–2963. [Google Scholar] [CrossRef] [Green Version]
  108. Gladfelter, A.S.; Bose, I.; Zyla, T.R.; Bardes, E.S.G.; Lew, D.J. Septin ring assembly involves cycles of GTP loading and hydrolysis by Cdc42p. J. Cell Biol. 2002, 156, 315–326. [Google Scholar] [CrossRef]
  109. Smith, G.R.; Givan, S.A.; Cullen, P.; Sprague, G.F. GTPase-Activating Proteins for Cdc42. Eukaryot. Cell 2002, 1, 469–480. [Google Scholar] [CrossRef] [Green Version]
  110. Marquitz, A.R.; Harrison, J.C.; Bose, I.; Zyla, T.R.; McMillan, J.N.; Lew, D.J. The Rho-GAP Bem2p plays a GAP-independent role in the morphogenesis checkpoint. EMBO J. 2002, 21, 4012–4025. [Google Scholar] [CrossRef] [Green Version]
  111. Caviston, J.P.; Longtine, M.; Pringle, J.R.; Bi, E. The role of Cdc42p GTPase-activating proteins in assembly of the septin ring in yeast. Mol. Biol. Cell 2003, 14, 4051–4066. [Google Scholar] [CrossRef] [Green Version]
  112. Mukherjee, D.; Sen, A.; Boettner, D.R.; Fairn, G.D.; Schlam, D.; Bonilla Valentin, F.J.; McCaffery, J.M.; Hazbun, T.; Staiger, C.J.; Grinstein, S.; et al. Bem3, a Cdc42 GTPase-activating protein, traffics to an intracellular compartment and recruits the secretory Rab GTPase Sec4 to endomembranes. J. Cell Sci. 2013, 126, 4560–4571. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  113. Mukherjee, D.; Coon, B.G.; Edwards, D.F.; Hanna, C.B.; Longhi, S.A.; McCaffery, J.M.; Wendland, B.; Retegui, L.A.; Bi, E.; Aguilar, R.C. The yeast endocytic protein Epsin 2 functions in a cell-division signaling pathway. J. Cell Sci. 2009, 122, 2453–2463. [Google Scholar] [CrossRef] [Green Version]
  114. Wild, A.C.; Yu, J.W.; Lemmon, M.A.; Blumer, K.J. The p21-activated Protein Kinase-related Kinase Cla4 Is a Coincidence Detector of Signaling by Cdc42 and Phosphatidylinositol 4-Phosphate. J. Biol. Chem. 2004, 279, 17101–17110. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  115. Lechler, T.; Shevchenko, A.; Shevchenko, A.; Li, R. Direct Involvement of Yeast Type I Myosins in Cdc42-Dependent Actin Polymerization. J. Cell Biol. 2000, 148, 363–374. [Google Scholar] [CrossRef] [Green Version]
  116. Goryachev, A.B.; Leda, M. Many roads to symmetry breaking: Molecular mechanisms and theoretical models of yeast cell polarity. Mol. Biol. Cell 2017, 28, 370–380. [Google Scholar] [CrossRef]
  117. Howell, A.S.; Savage, N.S.; Johnson, S.A.; Bose, I.; Wagner, A.W.; Zyla, T.R.; Nijhout, H.F.; Reed, M.C.; Goryachev, A.B.; Lew, D.J. Singularity in Polarization: Rewiring Yeast Cells to Make Two Buds. Cell 2009, 139, 731–743. [Google Scholar] [CrossRef] [Green Version]
  118. Kamada, Y.; Qadota, H.; Python, C.P.; Anraku, Y.; Ohya, Y.; Levin, D.E. Activation of yeast protein kinase C by Rho1 GTPase. J. Biol. Chem. 1996, 271, 9193–9196. [Google Scholar] [CrossRef] [Green Version]
  119. Qadota, H.; Python, C.P.; Inoue, S.B.; Arisawa, M.; Anraku, Y.; Zheng, Y.; Watanabe, T.; Levin, D.E.; Ohya, Y. Identification of yeast Rho1p GTPase as a regulatory subunit of 1, 3-β-glucan synthase. Science 1996, 272, 279–281. [Google Scholar] [CrossRef]
  120. Saka, A.; Abe, M.; Okano, H.; Minemura, M.; Qadota, H.; Utsugi, T.; Mino, A.; Tanaka, K.; Takai, Y.; Ohya, Y. Complementing Yeast rho1 Mutation Groups with Distinct Functional Defects. J. Biol. Chem. 2001, 276, 46165–46171. [Google Scholar] [CrossRef] [Green Version]
  121. Diepeveen, E.T.; De la Cruz, L.I.; Laan, L. Evolutionary dynamics in the fungal polarization network, a mechanistic perspective. Biophys. Rev. 2017, 9, 375–387. [Google Scholar] [CrossRef] [Green Version]
  122. Reynaud, A.; Facca, C.; Sor, F.; Faye, G. Disruption and functional analysis of six ORFs of chromosome IV: YDL103c (QRI1), YDL105w (QRI2), YDL112w (TRM3), YDL113c, YDL116w (NUP84) and YDL167c (NRP1). Yeast 2001, 18, 273–282. [Google Scholar] [CrossRef]
  123. Outten, C.E.; Falk, R.L.; Culotta, V.C. Cellular factors required for protection from hyperoxia toxicity in Saccharomyces cerevisiae. Biochem. J. 2005, 388, 93–101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  124. Buchan, J.R.; Muhlrad, D.; Parker, R. P bodies promote stress granule assembly in Saccharomyces cerevisiae. J. Cell Biol. 2008, 183, 441–455. [Google Scholar] [CrossRef] [PubMed]
  125. Hoepfner, D.; Helliwell, S.B.; Sadlish, H.; Schuierer, S.; Filipuzzi, I.; Brachat, S.; Bhullar, B.; Plikat, U.; Abraham, Y.; Altorfer, M.; et al. High-resolution chemical dissection of a model eukaryote reveals targets, pathways and gene functions. Microbiol. Res. 2014, 169, 107–120. [Google Scholar] [CrossRef]
  126. 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, 146–158. [Google Scholar] [CrossRef] [Green Version]
  127. Hogan, D.J.; Riordan, D.P.; Gerber, A.P.; Herschlag, D.; Brown, P.O. Diverse RNA-Binding Proteins Interact with Functionally Related Sets of RNAs, Suggesting an Extensive Regulatory System. PLoS Biol. 2008, 6, e255. [Google Scholar] [CrossRef] [Green Version]
  128. Sharifpoor, S.; Van Dyk, D.; Costanzo, M.; Baryshnikova, A.; Friesen, H.; Douglas, A.C.; Youn, J.-Y.; VanderSluis, B.; Myers, C.L.; Papp, B.; et al. Functional wiring of the yeast kinome revealed by global analysis of genetic network motifs. Genome Res. 2012, 22, 791–801. [Google Scholar] [CrossRef] [Green Version]
  129. Etienne-Manneville, S. Cdc42—The centre of polarity. J. Cell Sci. 2004, 117, 1291–1300. [Google Scholar] [CrossRef] [Green Version]
  130. Johnson, D.I. Cdc42: An Essential Rho-Type GTPase Controlling Eukaryotic Cell Polarity. Microbiol. Mol. Biol. Rev. 1999, 63, 54–105. [Google Scholar] [CrossRef] [Green Version]
  131. Lee, C.; Occhipinti, P.; Gladfelter, A.S. PolyQ-dependent RNA–protein assemblies control symmetry breaking. J. Cell Biol. 2015, 208, 533–544. [Google Scholar] [CrossRef]
  132. Kato, M.; McKnight, S.L. Cross-β Polymerization of Low Complexity Sequence Domains. Cold Spring Harb. Perspect. Biol. 2017, 9, a023598. [Google Scholar] [CrossRef] [Green Version]
  133. Wuthrich, K. Protein structure determination in solution by nuclear magnetic resonance spectroscopy. Science 1989, 243, 45–50. [Google Scholar] [CrossRef] [Green Version]
  134. Hrdlickova, R.; Toloue, M.; Tian, B. RNA-Seq methods for transcriptome analysis: RNA-Seq. Wiley Interdiscip. Rev. RNA 2017, 8, e1364. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  135. Keren, L.; Hausser, J.; Lotan-Pompan, M.; Vainberg Slutskin, I.; Alisar, H.; Kaminski, S.; Weinberger, A.; Alon, U.; Milo, R.; Segal, E. Massively Parallel Interrogation of the Effects of Gene Expression Levels on Fitness. Cell 2016, 166, 1282–1294.e18. [Google Scholar] [CrossRef]
  136. Vendel, K.J.A.; Tschirpke, S.; Shamsi, F.; Dogterom, M.; Laan, L. Minimal in vitro systems shed light on cell polarity. J. Cell Sci. 2019, 132, jcs217554. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  137. Heym, R.G.; Zimmermann, D.; Edelmann, F.T.; Israel, L.; Ökten, Z.; Kovar, D.R.; Niessing, D. In vitro reconstitution of an mRNA-transport complex reveals mechanisms of assembly and motor activation. J. Cell Biol. 2013, 203, 971–984. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  138. Cerulus, B.; New, A.M.; Pougach, K.; Verstrepen, K.J. Noise and Epigenetic Inheritance of Single-Cell Division Times Influence Population Fitness. Curr. Biol. 2016, 26, 1138–1147. [Google Scholar] [CrossRef] [Green Version]
  139. Xue, Y.; Acar, M. Mechanisms for the epigenetic inheritance of stress response in single cells. Curr. Genet. 2018, 64, 1221–1228. [Google Scholar] [CrossRef]
  140. Stark, C. BioGRID: A general repository for interaction datasets. Nucleic Acids Res. 2006, 34, D535–D539. [Google Scholar] [CrossRef] [Green Version]
  141. Michel, A.H.; Hatakeyama, R.; Kimmig, P.; Arter, M.; Peter, M.; Matos, J.; De Virgilio, C.; Kornmann, B. Functional mapping of yeast genomes by saturated transposition. Elife 2017, 6, e23570. [Google Scholar] [CrossRef]
  142. Kryazhimskiy, S.; Rice, D.P.; Jerison, E.R.; Desai, M.M. Global epistasis makes adaptation predictable despite sequence-level stochasticity. Science 2014, 344, 1519–1522. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  143. Fields, S.; Song, O. A novel genetic system to detect protein–protein interactions. Nature 1989, 340, 245–246. [Google Scholar] [CrossRef] [PubMed]
  144. Jares-Erijman, E.A.; Jovin, T.M. FRET imaging. Nat. Biotechnol. 2003, 21, 1387–1395. [Google Scholar] [CrossRef] [PubMed]
  145. Phizicky, E.M.; Fields, S. Protein-protein interactions: Methods for detection and analysis. Microbiol. Rev. 1995, 59, 94–123. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  146. Foltman, M.; Sanchez-Diaz, A. Studying Protein–Protein Interactions in Budding Yeast Using Co-immunoprecipitation. In Yeast Cytokinesis: Methods and Protocols; Sanchez-Diaz, A., Perez, P., Eds.; Springer: New York, NY, USA, 2016; pp. 239–256. ISBN 978-1-4939-3145-3. [Google Scholar]
  147. Nishimura, K.; Fukagawa, T.; Takisawa, H.; Kakimoto, T.; Kanemaki, M. An auxin-based degron system for the rapid depletion of proteins in nonplant cells. Nat. Methods 2009, 6, 917–922. [Google Scholar] [CrossRef] [PubMed]
  148. Papagiannakis, A.; De Jonge, J.J.; Zhang, Z.; Heinemann, M. Quantitative characterization of the auxin-inducible degron: A guide for dynamic protein depletion in single yeast cells. Sci. Rep. 2017, 7, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  149. Jost, A.P.-T.; Weiner, O.D. Probing Yeast Polarity with Acute, Reversible, Optogenetic Inhibition of Protein Function. ACS Synth. Biol. 2015, 4, 1077–1085. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  150. Jensen, E. Technical review: In situ hybridization. Anat. Rec. 2014, 297, 1349–1353. [Google Scholar] [CrossRef]
  151. Pelillo, M.; Poli, I.; Roli, A.; Serra, R.; Slanzi, D.; Villani, M. (Eds.) Artificial Life and Evolutionary Computation; Communications in Computer and Information Science; Springer International Publishing: Cham, Switzerland, 2018; Volume 830, ISBN 978-3-319-78657-5. [Google Scholar]
  152. Giaever, G.; Chu, A.M.; Ni, L.; Connelly, C.; Riles, L.; Véronneau, S.; Dow, S.; Lucau-Danila, A.; Anderson, K.; Andre, B. Functional profiling of the Saccharomyces cerevisiae genome. Nature 2002, 418, 387–391. [Google Scholar] [CrossRef]
  153. Gavin, A.-C.; Bösche, M.; Krause, R.; Grandi, P.; Marzioch, M.; Bauer, A.; Schultz, J.; Rick, J.M.; Michon, A.-M.; Cruciat, C.-M. Functional organization of the yeast proteome by systematic analysis of protein complexes. Nature 2002, 415, 141. [Google Scholar] [CrossRef]
  154. Tarassov, K.; Messier, V.; Landry, C.R.; Radinovic, S.; Molina, M.M.S.; Shames, I.; Malitskaya, Y.; Vogel, J.; Bussey, H.; Michnick, S.W. An in vivo map of the yeast protein interactome. Science 2008, 320, 1465–1470. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  155. Garí Marsol, E.; Volpe, T.; Wang, H.; Gallego, C.; Futcher, B.; Aldea, M. Whi3 binds the mRNA of the G (1) cyclin CLN3 to modulate cell fate in budding yeast. Genes Dev. 2001, 15, 2803–2808. [Google Scholar]
  156. Colomina, N.; Ferrezuelo, F.; Wang, H.; Aldea, M.; Garí, E. Whi3, a Developmental Regulator of Budding Yeast, Binds a Large Set of mRNAs Functionally Related to the Endoplasmic Reticulum. J. Biol. Chem. 2008, 283, 28670–28679. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  157. Wittenberg, C.; Sugimoto, K.; Reed, S.I. G1-specific cyclins of S. cerevisiae: Cell cycle periodicity, regulation by mating pheromone, and association with the p34CDC28 protein kinase. Cell 1990, 62, 225–237. [Google Scholar] [CrossRef]
  158. Leeuw, T.; Wu, C.; Schrag, J.D.; Whiteway, M.; Thomas, D.Y.; Leberer, E. Interaction of a G-protein β-subunit with a conserved sequence in Ste20/PAK family protein kinases. Nature 1998, 391, 191. [Google Scholar] [CrossRef]
  159. Wiget, P.; Shimada, Y.; Butty, A.-C.; Bi, E.; Peter, M. Site-specific regulation of the GEF Cdc24p by the scaffold protein Far1p during yeast mating. EMBO J. 2004, 23, 1063–1074. [Google Scholar] [CrossRef] [Green Version]
  160. Camus, C.; Boy-Marcotte, E.; Jacquet, M. Two subclasses of guanine exchange factor (GEF) domains revealed by comparison of activities of chimeric genes constructed from CDC25, SDC25 and BUD5 in Saccharomyces cerevisiae. Mol. Gen. Genet. MGG 1994, 245, 167–176. [Google Scholar] [CrossRef]
  161. Kozminski, K.G.; Beven, L.; Angerman, E.; Tong, A.H.Y.; Boone, C.; Park, H.-O. Interaction between a Ras and a Rho GTPase couples selection of a growth site to the development of cell polarity in yeast. Mol. Biol. Cell 2003, 14, 4958–4970. [Google Scholar] [CrossRef] [Green Version]
  162. Koch, G.; Tanaka, K.; Masuda, T.; Yamochi, W.; Nonaka, H.; Takai, Y. Association of the Rho family small GTP-binding proteins with Rho GDP dissociation inhibitor (Rho GDI) in Saccharomyces cerevisiae. Oncogene 1997, 15, 417. [Google Scholar] [CrossRef]
  163. Bose, I.; Irazoqui, J.E.; Moskow, J.J.; Bardes, E.S.G.; Zyla, T.R.; Lew, D.J. Assembly of Scaffold-mediated Complexes Containing Cdc42p, the Exchange Factor Cdc24p, and the Effector Cla4p Required for Cell Cycle-regulated Phosphorylation of Cdc24p. J. Biol. Chem. 2001, 276, 7176–7186. [Google Scholar] [CrossRef] [Green Version]
  164. Kozubowski, L.; Saito, K.; Johnson, J.M.; Howell, A.S.; Zyla, T.R.; Lew, D.J. Symmetry-Breaking Polarization Driven by a Cdc42p GEF-PAK Complex. Curr. Biol. 2008, 18, 1719–1726. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  165. Gulli, M.-P.; Jaquenoud, M.; Shimada, Y.; Niederhäuser, G.; Wiget, P.; Peter, M. Phosphorylation of the Cdc42 exchange factor Cdc24 by the PAK-like kinase Cla4 may regulate polarized growth in yeast. Mol. Cell 2000, 6, 1155–1167. [Google Scholar] [CrossRef]
  166. Sagot, I.; Klee, S.K.; Pellman, D. Yeast formins regulate cell polarity by controlling the assembly of actin cables. Nat. Cell Biol. 2002, 4, 42. [Google Scholar] [CrossRef] [PubMed]
  167. Adams, A.E. Relationship of actin and tubulin distribution to bud growth in wild- type and morphogenetic-mutant Saccharomyces cerevisiae. J. Cell Biol. 1984, 98, 934–945. [Google Scholar] [CrossRef]
  168. Wendland, B. Epsins: Adaptors in endocytosis? Nat. Rev. Mol. Cell Biol. 2002, 3, 971. [Google Scholar] [CrossRef]
  169. Goode, B.L.; Eskin, J.A.; Wendland, B. Actin and Endocytosis in Budding Yeast. Genetics 2015, 199, 315–358. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Venn diagram representation of the hierarchical, nonmodular pathway structure for protein−protein interactions in budding yeast polarity establishment (Table A1 for references).
Figure 1. Venn diagram representation of the hierarchical, nonmodular pathway structure for protein−protein interactions in budding yeast polarity establishment (Table A1 for references).
Cells 09 02534 g001
Figure 2. Mechanistic model for the haploid bud scar pathway, the reaction−diffusion pathway including the GAPs, and the hypothesized role of Nrp1 in the timing pathway.
Figure 2. Mechanistic model for the haploid bud scar pathway, the reaction−diffusion pathway including the GAPs, and the hypothesized role of Nrp1 in the timing pathway.
Cells 09 02534 g002
Figure 3. Mind map representation of the bud scar protein Bud1 summarizing its most important properties from literature, such as localization, domains and interactions.
Figure 3. Mind map representation of the bud scar protein Bud1 summarizing its most important properties from literature, such as localization, domains and interactions.
Cells 09 02534 g003
Figure 4. Mind map representation of the GTPase activating protein (GAP) Bem3 summarizing its most important properties from literature, such as localization, domains and interactions.
Figure 4. Mind map representation of the GTPase activating protein (GAP) Bem3 summarizing its most important properties from literature, such as localization, domains and interactions.
Cells 09 02534 g004
Figure 5. Mind map representation of Nrp1 summarizing its most important properties from literature, such as localization, domains and interactions.
Figure 5. Mind map representation of Nrp1 summarizing its most important properties from literature, such as localization, domains and interactions.
Cells 09 02534 g005
Figure 6. Flow chart for understanding the roles of proteins inside a complex network. Level of detail increases from top to bottom, while the enumerations of accompanying experiments corresponding to a level are ordered in decreasing feasibility of high-throughput data generation.
Figure 6. Flow chart for understanding the roles of proteins inside a complex network. Level of detail increases from top to bottom, while the enumerations of accompanying experiments corresponding to a level are ordered in decreasing feasibility of high-throughput data generation.
Cells 09 02534 g006
Table 1. Definitions of important terms used throughout the introduction.
Table 1. Definitions of important terms used throughout the introduction.
TermDefinition
Protein networkGroup of proteins with physical interactions together performing a function
ConnectivityDegree to which parts of the network are embedded with other parts in the network. In this sense, it can be received as the reciprocal of modularity.
ModularityPotential to group parts of a protein network given a certain representation of the protein network (e.g., in terms of mechanisms, genetic or physical interactions)
Hub proteinHighly connected protein in a network (often essential)
NeutralityNo consequence of a mutation to phenotype (in current environment)
HierarchyClear layering of pathways inside a protein network
RedundancyMultiple mechanisms that can to some extent interchangeably contribute to the same function
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Daalman, W.K.-G.; Sweep, E.; Laan, L. The Path towards Predicting Evolution as Illustrated in Yeast Cell Polarity. Cells 2020, 9, 2534. https://doi.org/10.3390/cells9122534

AMA Style

Daalman WK-G, Sweep E, Laan L. The Path towards Predicting Evolution as Illustrated in Yeast Cell Polarity. Cells. 2020; 9(12):2534. https://doi.org/10.3390/cells9122534

Chicago/Turabian Style

Daalman, Werner Karl-Gustav, Els Sweep, and Liedewij Laan. 2020. "The Path towards Predicting Evolution as Illustrated in Yeast Cell Polarity" Cells 9, no. 12: 2534. https://doi.org/10.3390/cells9122534

APA Style

Daalman, W. K. -G., Sweep, E., & Laan, L. (2020). The Path towards Predicting Evolution as Illustrated in Yeast Cell Polarity. Cells, 9(12), 2534. https://doi.org/10.3390/cells9122534

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop