Modelling of Yeast Mating Reveals Robustness Strategies for Cell-Cell Interactions

Mating of budding yeast cells is a model system for studying cell-cell interactions. Haploid yeast cells secrete mating pheromones that are sensed by the partner which responds by growing a mating projection toward the source. The two projections meet and fuse to form the diploid. Successful mating relies on precise coordination of dynamic extracellular signals, signaling pathways, and cell shape changes in a noisy background. It remains elusive how cells mate accurately and efficiently in a natural multi-cell environment. Here we present the first stochastic model of multiple mating cells whose morphologies are driven by pheromone gradients and intracellular signals. Our novel computational framework encompassed a moving boundary method for modeling both a-cells and α-cells and their cell shape changes, the extracellular diffusion of mating pheromones dynamically coupled with cell polarization, and both external and internal noise. Quantification of mating efficiency was developed and tested for different model parameters. Computer simulations revealed important robustness strategies for mating in the presence of noise. These strategies included the polarized secretion of pheromone, the presence of the α-factor protease Bar1, and the regulation of sensing sensitivity; all were consistent with data in the literature. In addition, we investigated mating discrimination, the ability of an a-cell to distinguish between α-cells either making or not making α-factor, and mating competition, in which multiple a-cells compete to mate with one α-cell. Our simulations were consistent with previous experimental results. Moreover, we performed a combination of simulations and experiments to estimate the diffusion rate of the pheromone a-factor. In summary, we constructed a framework for simulating yeast mating with multiple cells in a noisy environment, and used this framework to reproduce mating behaviors and to identify strategies for robust cell-cell interactions.     

Role of α-Factor and the Mfα1 α-Factor Precursor in Mating in Yeast

The peptide pheromones secreted by a and α cells (called a-factor and α-factor, respectively) are each encoded by two structural genes. For strains of either mating type, addition of exogenous pheromone does not alleviate the mating defect of mutants with disruptions of both structural genes. In addition, a particular insertion mutation in the major α-factor structural gene (MFα1) that should result in an altered product inhibits α mating. These results suggested that the pheromone precursors (the MFα1 pro region in particular) might play a second role in mating separate from the role of pheromone production. To analyze the role of α-factor and the MFα1 precursor in α mating, we have constructed two classes of mutants. The mating defects of mutants that should produce the MFα1 pro region peptide but no α-factor could not be alleviated by addition of exogenous α-factor in crosses to a wild-type a strain, indicating that the previous results were not due to an inability of the disruption mutants to produce the pro region peptide. Mutants able to produce α-factor, but with a variety of alterations in MFα1 precursor structure, mated at levels proportional to the levels of α-factor produced, suggesting that the only role of the α-factor precursor in mating is to produce α-factor. Both of these results argue against a role for the MFα1 pro region separate from its role in α-factor production. We also describe results that show that in vivo production of α-factor'' (the form of α-factor encoded by one of the two α-factor repeats of MFα2) is equivalent to the major form of α-factor for induction of all responses necessary for mating. We discuss the implications of these results on the role of the pheromones in mating.     

Reverse Engineering a Signaling Network Using Alternative Inputs

One of the goals of systems biology is to reverse engineer in a comprehensive fashion the arrow diagrams of signal transduction systems. An important tool for ordering pathway components is genetic epistasis analysis, and here we present a strategy termed Alternative Inputs (AIs) to perform systematic epistasis analysis. An alternative input is defined as any genetic manipulation that can activate the signaling pathway instead of the natural input. We introduced the concept of an “AIs-Deletions matrix” that summarizes the outputs of all combinations of alternative inputs and deletions. We developed the theory and algorithms to construct a pairwise relationship graph from the AIs-Deletions matrix capturing both functional ordering (upstream, downstream) and logical relationships (AND, OR), and then interpreting these relationships into a standard arrow diagram. As a proof-of-principle, we applied this methodology to a subset of genes involved in yeast mating signaling. This experimental pilot study highlights the robustness of the approach and important technical challenges. In summary, this research formalizes and extends classical epistasis analysis from linear pathways to more complex networks, facilitating computational analysis and reconstruction of signaling arrow diagrams.     

Plasma membrane aminoglycerolipid flippase function is required for signaling competence in the yeast mating pheromone response pathway

The class 4 P-type ATPases (“flippases”) maintain membrane asymmetry by translocating phosphatidylethanolamine and phosphatidylserine from the outer leaflet to the cytosolic leaflet of the plasma membrane. In Saccharomyces cerevisiae, five related gene products (Dnf1, Dnf2, Dnf3, Drs2, and Neo1) are implicated in flipping of phosphatidylethanolamine, phosphatidylserine, and phosphatidylcholine. In MATa cells responding to α-factor, we found that Dnf1, Dnf2, and Dnf3, as well as the flippase-activating protein kinase Fpk1, localize at the projection (“shmoo”) tip where polarized growth is occurring and where Ste5 (the central scaffold protein of the pheromone-initiated MAPK cascade) is recruited. Although viable, a MATa dnf1∆ dnf2∆ dnf3∆ triple mutant exhibited a marked decrease in its ability to respond to α-factor, which we could attribute to pronounced reduction in Ste5 stability resulting from an elevated rate of its Cln2⋅Cdc28-initiated degradation. Similarly, a MATa dnf1∆ dnf3∆ drs2∆ triple mutant also displayed marked reduction in its ability to respond to α-factor, which we could attribute to inefficient recruitment of Ste5 to the plasma membrane due to severe mislocalization of the cellular phosphatidylinositol 4-phosphate and phosphatidylinositol 4,5-bisphosphate pools. Thus proper remodeling of plasma membrane aminoglycerolipids and phosphoinositides is necessary for efficient recruitment, stability, and function of the pheromone signaling apparatus.     

Yeast G-proteins mediate directional sensing and polarization behaviors in response to changes in pheromone gradient direction

Yeast cells polarize by projecting up mating pheromone gradients, a classic cell polarity behavior. However, these chemical gradients may shift direction. We examine how yeast cells sense and respond to a 180^o switch in the direction of microfluidically generated pheromone gradients. We identify two behaviors: at low concentrations of α-factor, the initial projection grows by bending, whereas at high concentrations, cells form a second projection toward the new source. Mutations that increase heterotrimeric G-protein activity expand the bending-growth morphology to high concentrations; mutations that increase Cdc42 activity result in second projections at low concentrations. Gradient-sensing projection bending requires interaction between Gβγ and Cdc24, whereas gradient-nonsensing projection extension is stimulated by Bem1 and hyperactivated Cdc42. Of interest, a mutation in Gα affects both bending and extension. Finally, we find a genetic perturbation that exhibits both behaviors. Overexpression of the formin Bni1, a component of the polarisome, makes both bending-growth projections and second projections at low and high α-factor concentrations, suggesting a role for Bni1 downstream of the heterotrimeric G-protein and Cdc42 during gradient sensing and response. Thus we demonstrate that G-proteins modulate in a ligand-dependent manner two fundamental cell-polarity behaviors in response to gradient directional change.     

Robust Spatial Sensing of Mating Pheromone Gradients by Yeast Cells

Projecting or moving up a chemical gradient is a universal behavior of living organisms. We tested the ability of S. cerevisiae a-cells to sense and respond to spatial gradients of the mating pheromone α-factor produced in a microfluidics chamber; the focus was on bar1Δ strains, which do not degrade the pheromone input. The yeast cells exhibited good accuracy with the mating projection typically pointing in the correct direction up the gradient (∼80% under certain conditions), excellent sensitivity to shallow gradients, and broad dynamic range so that gradient-sensing was relatively robust over a 1000-fold range of average α-factor concentrations. Optimal directional sensing occurred at lower concentrations (5 nM) close to the K_d of the receptor and with steeper gradient slopes. Pheromone supersensitive mutations (sst2Δ and ste2^300Δ) that disrupt the down-regulation of heterotrimeric G-protein signaling caused defects in both sensing and response. Interestingly, yeast cells employed adaptive mechanisms to increase the robustness of the process including filamentous growth (i.e. directional distal budding) up the gradient at low pheromone concentrations, bending of the projection to be more aligned with the gradient, and forming a more accurate second projection when the first projection was in the wrong direction. Finally, the cells were able to amplify a shallow external gradient signal of α-factor to produce a dramatic polarization of signaling proteins at the front of the cell. Mathematical modeling revealed insights into the mechanism of this amplification and how the supersensitive mutants can disrupt accurate polarization. Together, these data help to specify and elucidate the abilities of yeast cells to sense and respond to spatial gradients of pheromone.     

Mutants of SACCHAROMYCES CEREVISIAE Resistant to the α Mating-Type Factor

Mutants that are resistant to α-factor have been isolated from a mating-type haploid strains of yeast by direct selection on agar medium containing partially purified α-factor. All resistant mutants isolated were found to be sterile. They were characterized and compared with mutants previously isolated as nonmating. Among 93 able to mate at low frequency and to sporulate, none showed linkage to the mating-type locus. The results support the hypothesis that the response to α-factor by cells of mating-type a is essential for mating.     

Deletion of the Schizophyllum Commune Aα Locus: The Roles of Aα Y and Z Mating-Type Genes

The Aα locus is one of four master regulatory loci that determine mating type and regulate sexual development in Schizophyllum commune. We have made a plasmid containing a URA1 gene disruption of the Aα Y1 gene. Y1 is the sole Aα gene in Aα1 strains. We used the plasmid construction to produce an Aα null (i.e., AαΔ) strain by replacing the genomic Y1 gene with URA1 in an Aα1 strain. To characterize the role of the Aα genes in the regulation of sexual development, we transformed various Aα Y and Z alleles into AαΔ strains and examined the acquired mating types and mating abilities of the transformants. These experiments demonstrate that the Aα Y gene is not essential for fungal viability and growth, that a solitary Z Aα mating-type gene does not itself activate development, that Aβ proteins are sufficient to activate the A developmental pathway in the absence of Aα proteins and confirm that Y and Z genes are the sole determinants of Aα mating type. The data from these experiments support and refine our model of the regulation of A-pathway events by Y and Z proteins.     

Nonsex Genes in the Mating Type Locus of Candida albicans Play Roles in a/α Biofilm Formation,Including Impermeability and Fluconazole Resistance

The mating type locus (MTL) of Candida albicans contains the mating type genes and has, therefore, been assumed to play an exclusive role in the mating process. In mating-incompetent a/α cells, two of the mating type genes, MTL a1 and MTLα2, encode components of the a1-α2 corepressor that suppresses mating and switching. But the MTL locus of C. albicans also contains three apparently unrelated “nonsex” genes (NSGs), PIK, PAP and OBP, the first two essential for growth. Since it had been previously demonstrated that deleting either the a/α copy of the entire MTL locus, or either MTLa1 or MTLα2, affected virulence, we hypothesized that the NSGs in the MTL locus may also play a role in pathogenesis. Here by mutational analysis, it is demonstrated that both the mating type and nonsex genes in the MTL locus play roles in a/α biofilm formation, and that OBP is essential for impermeability and fluconazole resistance.     

The Stem Cell-Expressed Receptor Lgr5 Possesses Canonical and Functionally Active Molecular Determinants Critical to β-arrestin-2 Recruitment

Lgr5 is a membrane protein related to G protein-coupled receptors (GPCR)s whose expression identifies stem cells in multiple tissues and is strongly correlated with cancer. Despite the recent identification of endogenous ligands for Lgr5, its mode of signaling remains enigmatic. The ability to couple to G proteins and βarrestins are classical molecular behaviors of GPCRs that have yet to be observed for Lgr5. Therefore, the goal of this study was to determine if Lgr5 can engage a classical GPCR behavior and elucidate the molecular determinants of this process. Structural analysis of Lgr5 revealed several motifs consistent with its ability to recruit βarr2. Among them, a “SSS” serine cluster located at amino acid position 873-875 within the C-terminal tail (C-tail), is in a region consistent with other GPCRs that bind βarr2 with high-affinity. To test its functionality, a ligand-independent βarr2 translocation assay was implemented. We show that Lgr5 recruits βarr2 and that the “SSS” amino acids (873-875) are absolutely critical to this process. We also demonstrate that for full efficacy, this cluster requires other Lgr5 C-tail serines that were previously shown to be important for constitutive and βarr2 independent internalization of Lgr5. These data are proof of principle that a classical GPCR behavior can be manifested by Lgr5. The existence of alternative ligands or missing effectors of Lgr5 that scaffold this classical GPCR behavior and the downstream signaling pathways engaged should be considered. Characterizing Lgr5 signaling will be invaluable for assessing its role in tissue maintenance, repair, and disease.     

Localized secretion of acid phosphatase reflects the pattern of cell surface growth in saccharomyces cerevisiae

Secretion of cell wall-bound acid phosphatase by Saccharomyces cerevisiae occurs along a restricted portion of the cell surface. Acid phosphatase activity produced during derepressed synthesis on a phosphate-limited growth medium is detected with an enzyme-specific stain and is localized initially to the bud portion of a dividing cell. After two to three generations of phosphate-limited growth, most of the cells can be stained; if further phosphatase synthesis is repressed by growth in excess phosphate, dividing cells are produced in which the parent but not the bud can be stained. Budding growth is interrupted in α-mating-type cells by a pheromone (α-factor) secreted by the opposite mating type; cell surface growth continues in the presence of α-factor and produces a characteristic cell tip. When acid phosphatase synthesis is initiated during α-factor treatment, only the cell tip can br stained; when phosphate synthesis is repressed during α-factor treatment, the cell body but not the tip can be stained. A mixture of derepressed α cells and phosphatase-negative α cells form zygotes in which mainly one parent cell surface can be stained. The cell cycle mutant, cdc 24 (Hartwell, L.H. 1971. Exp. Cell Res. 69:265-276), fails to bud and, instead, expands symmetrically as a sphere at a nonpermissive temperature (37 degrees C). This mutant does not form a cell tip during α-factor treatment at 37 degrees C, and although acid phosphatade secretion occurs at this temperature, it is not localized. These results suggest that secretion reflects a polar mode of yeast cell- surface growth, and that this organization requires the cdc 24 gene product.     

Pheromone-regulated Genes Required for Yeast Mating Differentiation

Yeast cells mate by an inducible pathway that involves agglutination, mating projection formation, cell fusion, and nuclear fusion. To obtain insight into the mating differentiation of Saccharomyces cerevisiae, we carried out a large-scale transposon tagging screen to identify genes whose expression is regulated by mating pheromone. 91,200 transformants containing random lacZ insertions were screened for β-galactosidase (β-gal) expression in the presence and absence of α factor, and 189 strains containing pheromone-regulated lacZ insertions were identified. Transposon insertion alleles corresponding to 20 genes that are novel or had not previously been known to be pheromone regulated were examined for effects on the mating process. Mutations in four novel genes, FIG1, FIG2, KAR5/ FIG3, and FIG4 were found to cause mating defects. Three of the proteins encoded by these genes, Fig1p, Fig2p, and Fig4p, are dispensible for cell polarization in uniform concentrations of mating pheromone, but are required for normal cell polarization in mating mixtures, conditions that involve cell–cell communication. Fig1p and Fig2p are also important for cell fusion and conjugation bridge shape, respectively. The fourth protein, Kar5p/Fig3p, is required for nuclear fusion. Fig1p and Fig2p are likely to act at the cell surface as Fig1:: β-gal and Fig2::β-gal fusion proteins localize to the periphery of mating cells. Fig4p is a member of a family of eukaryotic proteins that contain a domain homologous to the yeast Sac1p. Our results indicate that a variety of novel genes are expressed specifically during mating differentiation to mediate proper cell morphogenesis, cell fusion, and other steps of the mating process.     

Non-Canonical Interleukin 23 Receptor Complex Assembly: P40 PROTEIN RECRUITS INTERLEUKIN 12 RECEPTOR β1 VIA SITE II AND INDUCES P19/INTERLEUKIN 23 RECEPTOR INTERACTION VIA SITE III*

IL-23, composed of the cytokine subunit p19 and the soluble α receptor subunit p40, binds to a receptor complex consisting of the IL-23 receptor (IL-23R) and the IL-12 receptor β1 (IL-12Rβ1). Complex formation was hypothesized to follow the “site I-II-III” architectural paradigm, with site I of p19 being required for binding to p40, whereas sites II and III of p19 mediate binding to IL-12Rβ1 and IL-23R, respectively. Here we show that the binding mode of p19 to p40 and of p19 to IL-23R follow the canonical site I and III paradigm but that interaction of IL-23 to IL-12Rβ1 is independent of site II in p19. Instead, binding of IL-23 to the cytokine binding module of IL-12Rβ1 is mediated by domains 1 and 2 of p40 via corresponding site II amino acids of IL-12Rβ1. Moreover, domains 2 and 3 of p40 were sufficient for complex formation with p19 and to induce binding of p19 to IL-23R. The Fc-tagged fusion protein of p40_D2D3/p19 did, however, not act as a competitive IL-23 antagonist but, at higher concentrations, induced proliferation via IL-23R but independent of IL-12Rβ1. On the basis of our experimental validation, we propose a non-canonical topology of the IL-23·IL-23R·IL-12Rβ1 complex. Furthermore, our data help to explain why p40 is an antagonist of IL-23 and IL-12 signaling and show that site II of p19 is dispensable for IL-23 signaling.     

Pleiotropic Mutations at the TUP1 Locus That Affect the Expression of Mating-Type-Dependent Functions in SACCHAROMYCES CEREVISIAE

The umr7–1 mutation, previously identified in a set of mutants that had been selected for defective UV-induced mutagenesis at CAN1, affects other cellular functions, including many of those regulated by the mating-type locus (MAT) in heterothallic Saccharomyces cerevisiae. The recessive umr7–1 allele, mapping approximately 20 cM distal to thr4 on chromosome III, causes clumpy growth in both a and α cells and has no apparent effect on a mating functions. However, α umr7 meiotic segregants fail to express several α-specific functions (e.g., high-frequency conjugation with a strains, secretion of the hormone α-factor and response to the hormone a-factor). In addition, α umr7 cells exhibit some a-specific characteristics, such as the barrier phenotype (Bar⁺) that prevents diffusion of α-factor and an increased mating frequency with α strains. The most striking property of α umr7 strains is their altered morphology, in which mitotic cells develop an asymmetric pear shape, like that of normal a cells induced to form "shmoos" by interaction with α-factor. Some a/α-specific diploid functions are also affected by umr7; instead of polar budding patterns, a/α umr7/umr7 diploids have medial budding like a/a, α/α and haploid strains. Moreover, a/α umr7/umr7 diploids have lost the ability to sporulate and are Bar⁺ like a or a/a strains. Revertant studies indicate that umr7–1 is a single point mutation. The umr7 mutant fails to complement mutants of both tup1 (selected for deoxythymidine monophosphate utilization) and cyc9 (selected for high iso-2-cytochrome c levels), and all three isolates have similar genetic and phenotypic properties. It is suggested that the product of this gene plays some common central role in the complex regulation of the expression of both MAT-dependent and MAT-independent functions.     

Ste18p Is a Positive Control Element in the Mating Process of Candida albicans

Heterotrimeric G proteins are an important class of eukaryotic signaling molecules that have been identified as central elements in the pheromone response pathways of many fungi. In the fungal pathogen Candida albicans, the STE18 gene (ORF19.6551.1) encodes a potential γ subunit of a heterotrimeric G protein; this protein contains the C-terminal CAAX box characteristic of γ subunits and has sequence similarity to γ subunits implicated in the mating pathways of a variety of fungi. Disruption of this gene was shown to cause sterility of MTLa mating cells and to block pheromone-induced gene expression and shmoo formation; deletion of just the CAAX box residues is sufficient to inactivate Ste18 function in the mating process. Intriguingly, ectopic expression behind the strong ACT1 promoter of either the Gα or the Gβ subunit of the heterotrimeric G protein is able to suppress the mating defect caused by deletion of the Gγ subunit and restore both pheromone-induced gene expression and morphology changes.     

Pheromone-encoding mRNA is transported to the yeast mating projection by specific RNP granules

Association of messenger RNAs with large complexes such as processing bodies (PBs) plays a pivotal role in regulating their translation and decay. Little is known about other possible functions of these assemblies. Exposure of haploid yeast cells, carrying mating type “a,” to “α pheromone” stimulates polarized growth resulting in a “shmoo” projection; it also induces synthesis of “a pheromone,” encoded by MFA2. In this paper, we show that, in response to α pheromone, MFA2 mRNA is assembled with two types of granules; both contain some canonical PB proteins, yet they differ in size, localization, motility, and sensitivity to cycloheximide. Remarkably, one type is involved in mRNA transport to the tip of the shmoo, whereas the other—in local translation in the shmoo. Normal assembly of these granules is critical for their movement, localization, and for mating. Thus, MFA2 mRNAs are transported to the shmoo tip, in complex with PB-like particles, where they are locally translated.     

HIF1-alpha functions as a tumor promoter in cancer associated fibroblasts,and as a tumor suppressor in breast cancer cells: Autophagy drives compartment-specific oncogenesis

Our recent studies have mechanistically implicated a loss of stromal Cav-1 expression and HIF1α-activation in driving the cancer-associated fibroblast phenotype, through the paracrine production of nutrients via autophagy and aerobic glycolysis. However, it remains unknown if HIF1α-activation is sufficient to confer the cancer-associated fibroblast phenotype. To test this hypothesis directly, we stably-expressed activated HIF1α in fibroblasts and then examined their ability to promote tumor growth using a xenograft model employing human breast cancer cells (MDA-MB-231). Fibroblasts harboring activated HIF1α showed a dramatic reduction in Cav-1 levels and a shift towards aerobic glycolysis, as evidenced by a loss of mitochondrial activity, and an increase in lactate production. Activated HIF1α also induced BNIP3 and BNIP3L expression, markers for the autophagic destruction of mitochondria. Most importantly, fibroblasts expressing activated HIF1α increased tumor mass by ∼2-fold and tumor volume by ∼3-fold, without a significant increase in tumor angiogenesis. In this context, HIF1α also induced an increase in the lymph node metastasis of cancer cells. Similar results were obtained by driving NFκB activation in fibroblasts, another inducer of autophagy. Thus, activated HIF1α is sufficient to functionally confer the cancer-associated fibroblast phenotype. It is also known that HIF1α expression is required for the induction of autophagy in cancer cells. As such, we next directly expressed activated HIF1α in MDA-MB-231 cells and assessed its effect on tumor growth via xenograft analysis. Surprisingly, activated HIF1α in cancer cells dramatically suppressed tumor growth, resulting in a 2-fold reduction in tumor mass and a three-fold reduction in tumor volume. We conclude that HIF1α activation in different cell types can either promote or repress tumorigenesis. Based on these studies, we suggest that autophagy in cancer-associated fibroblasts promotes tumor growth via the paracrine production of recycled nutrients, which can directly “feed” cancer cells. Conversely, autophagy in cancer cells represses tumor growth via their “self-digestion.” Thus, we should consider that the activities of various known oncogenes and tumor-suppressors may be compartment and cell-type specific, and are not necessarily an intrinsic property of the molecule itself. As such, other “classic” oncogenes and tumor suppressors will have to be re-evaluated to determine their compartment specific effects on tumor growth and metastasis. Lastly, our results provide direct experimental support for the recently proposed “autophagic tumor stroma model of cancer.”Key words: caveolin-1, autophagy, mitophagy, the Warburg effect, tumor stroma, hypoxia, HIF1A, NFκB, compartment-specific oncogenesis, cancer-associated fibroblasts