Genetic Analysis of Default Mating Behavior in Saccharomyces Cerevisiae

Haploid Saccharomyces cerevisiae cells find each other during conjugation by orienting their growth toward each other along pheromone gradients (chemotropism). However, when their receptors are saturated for pheromone binding, yeast cells must select a mate by executing a default pathway in which they choose a mating partner at random. We previously demonstrated that this default pathway requires the SPA2 gene. In this report we show that the default mating pathway also requires the AXL1, FUS1, FUS2, FUS3, PEA2, RVS161, and BNI1 genes. These genes, including SPA2, are also important for efficient cell fusion during chemotropic mating. Cells containing null mutations in these genes display defects in cell fusion that subtly affect mating efficiency. In addition, we found that the defect in default mating caused by mutations in SPA2 is partially suppressed by multiple copies of two genes, FUS2 and MFA2. These findings uncover a molecular relationship between default mating and cell fusion. Moreover, because axl1 mutants secrete reduced levels of a-factor and are defective at both cell fusion and default mating, these results reveal an important role for a-factor in cell fusion and default mating. We suggest that default mating places a more stringent requirement on some aspects of cell fusion than does chemotropic mating.     

Epistatic interaction maps relative to multiple metabolic phenotypes

An epistatic interaction between two genes occurs when the phenotypic impact of one gene depends on another gene, often exposing a functional association between them. Due to experimental scalability and to evolutionary significance, abundant work has been focused on studying how epistasis affects cellular growth rate, most notably in yeast. However, epistasis likely influences many different phenotypes, affecting our capacity to understand cellular functions, biochemical networks adaptation, and genetic diseases. Despite its broad significance, the extent and nature of epistasis relative to different phenotypes remain fundamentally unexplored. Here we use genome-scale metabolic network modeling to investigate the extent and properties of epistatic interactions relative to multiple phenotypes. Specifically, using an experimentally refined stoichiometric model for Saccharomyces cerevisiae, we computed a three-dimensional matrix of epistatic interactions between any two enzyme gene deletions, with respect to all metabolic flux phenotypes. We found that the total number of epistatic interactions between enzymes increases rapidly as phenotypes are added, plateauing at approximately 80 phenotypes, to an overall connectivity that is roughly 8-fold larger than the one observed relative to growth alone. Looking at interactions across all phenotypes, we found that gene pairs interact incoherently relative to different phenotypes, i.e. antagonistically relative to some phenotypes and synergistically relative to others. Specific deletion-deletion-phenotype triplets can be explained metabolically, suggesting a highly informative role of multi-phenotype epistasis in mapping cellular functions. Finally, we found that genes involved in many interactions across multiple phenotypes are more highly expressed, evolve slower, and tend to be associated with diseases, indicating that the importance of genes is hidden in their total phenotypic impact. Our predictions indicate a pervasiveness of nonlinear effects in how genetic perturbations affect multiple metabolic phenotypes. The approaches and results reported could influence future efforts in understanding metabolic diseases and the role of biochemical regulation in the cell.     

Derivation of genetic interaction networks from quantitative phenotype data

We have generalized the derivation of genetic-interaction networks from quantitative phenotype data. Familiar and unfamiliar modes of genetic interaction were identified and defined. A network was derived from agar-invasion phenotypes of mutant yeast. Mutations showed specific modes of genetic interaction with specific biological processes. Mutations formed cliques of significant mutual information in their large-scale patterns of genetic interaction. These local and global interaction patterns reflect the effects of gene perturbations on biological processes and pathways.     

The landscape of microbial phenotypic traits and associated genes

Bacteria and Archaea display a variety of phenotypic traits and can adapt to diverse ecological niches. However, systematic annotation of prokaryotic phenotypes is lacking. We have therefore developed ProTraits, a resource containing ∼545 000 novel phenotype inferences, spanning 424 traits assigned to 3046 bacterial and archaeal species. These annotations were assigned by a computational pipeline that associates microbes with phenotypes by text-mining the scientific literature and the broader World Wide Web, while also being able to define novel concepts from unstructured text. Moreover, the ProTraits pipeline assigns phenotypes by drawing extensively on comparative genomics, capturing patterns in gene repertoires, codon usage biases, proteome composition and co-occurrence in metagenomes. Notably, we find that gene synteny is highly predictive of many phenotypes, and highlight examples of gene neighborhoods associated with spore-forming ability. A global analysis of trait interrelatedness outlined clusters in the microbial phenotype network, suggesting common genetic underpinnings. Our extended set of phenotype annotations allows detection of 57 088 high confidence gene-trait links, which recover many known associations involving sporulation, flagella, catalase activity, aerobicity, photosynthesis and other traits. Over 99% of the commonly occurring gene families are involved in genetic interactions conditional on at least one phenotype, suggesting that epistasis has a major role in shaping microbial gene content.     

Codependency and mutual exclusivity for gene community detection from sparse single-cell transcriptome data

Single-cell RNA-seq (scRNA-seq) can be used to characterize cellular heterogeneity in thousands of cells. The reconstruction of a gene network based on coexpression patterns is a fundamental task in scRNA-seq analyses, and the mutual exclusivity of gene expression can be critical for understanding such heterogeneity. Here, we propose an approach for detecting communities from a genetic network constructed on the basis of coexpression properties. The community-based comparison of multiple coexpression networks enables the identification of functionally related gene clusters that cannot be fully captured through differential gene expression-based analysis. We also developed a novel metric referred to as the exclusively expressed index (EEI) that identifies mutually exclusive gene pairs from sparse scRNA-seq data. EEI quantifies and ranks the exclusive expression levels of all gene pairs from binary expression patterns while maintaining robustness against a low sequencing depth. We applied our methods to glioblastoma scRNA-seq data and found that gene communities were partially conserved after serum stimulation despite a considerable number of differentially expressed genes. We also demonstrate that the identification of mutually exclusive gene sets with EEI can improve the sensitivity of capturing cellular heterogeneity. Our methods complement existing approaches and provide new biological insights, even for a large, sparse dataset, in the single-cell analysis field.     

Pheromone-induced polarization is dependent on the Fus3p MAPK acting through the formin Bni1p

During mating, budding yeast cells reorient growth toward the highest concentration of pheromone. Bni1p, a formin homologue, is required for this polarized growth by facilitating cortical actin cable assembly. Fus3p, a pheromone-activated MAP kinase, is required for pheromone signaling and cell fusion. We show that Fus3p phosphorylates Bni1p in vitro, and phosphorylation of Bni1p in vivo during the pheromone response is dependent on Fus3p. fus3 mutants exhibited multiple phenotypes similar to bni1 mutants, including defects in actin and cell polarization, as well as Kar9p and cytoplasmic microtubule localization. Disruption of the interaction between Fus3p and the receptor-associated Galpha subunit caused similar mutant phenotypes. After pheromone treatment, Bni1p-GFP and Spa2p failed to localize to the cortex of fus3 mutants, and cell wall growth became completely unpolarized. Bni1p overexpression suppressed the actin assembly, cell polarization, and cell fusion defects. These data suggest a model wherein activated Fus3p is recruited back to the cortex, where it activates Bni1p to promote polarization and cell fusion.     

Genome shuffling: Progress and applications for phenotype improvement

Although rational method and global technique have been successfully applied in strain improvement respectively, the demand for engineering complex phenotypes required combinatorial approach. The technology of genome shuffling has been presented as a novel whole genome engineering approach for the rapid improvement of cellular phenotypes. This approach using recursive protoplast fusion with multi-parental strains offers the advantage of recombination throughout the entire genome without the necessity for genome sequence data or network information. Genome shuffling has been demonstrated as an effective method, which is not only for producing improved strain but also for providing information on complex phenotype. In this review we attempt to present the advantage of genome shuffling, introduce the procedure of this technology, summarize the applications of this approach for phenotype improvement and then give perspective on the development of this method in the future.     

A Network of Genes Antagonistic to the LIN-35 Retinoblastoma Protein of Caenorhabditis elegans

The Caenorhabditis elegans pRb ortholog, LIN-35, functions in a wide range of cellular and developmental processes. This includes a role of LIN-35 in nutrient utilization by the intestine, which it carries out redundantly with SLR-2, a zinc-finger protein. This and other redundant functions of LIN-35 were identified in genetic screens for mutations that display synthetic phenotypes in conjunction with loss of lin-35. To explore the intestinal role of LIN-35, we conducted a genome-wide RNA-interference-feeding screen for suppressors of lin-35; slr-2 early larval arrest. Of the 26 suppressors identified, 17 fall into three functional classes: (1) ribosome biogenesis genes, (2) mitochondrial prohibitins, and (3) chromatin regulators. Further characterization indicates that different categories of suppressors act through distinct molecular mechanisms. We also tested lin-35; slr-2 suppressors, as well as suppressors of the synthetic multivulval phenotype, to determine the spectrum of lin-35-synthetic phenotypes that could be suppressed following inhibition of these genes. We identified 19 genes, most of which are evolutionarily conserved, that can suppress multiple unrelated lin-35-synthetic phenotypes. Our study reveals a network of genes broadly antagonistic to LIN-35 as well as genes specific to the role of LIN-35 in intestinal and vulval development. Suppressors of multiple lin-35 phenotypes may be candidate targets for anticancer therapies. Moreover, screening for suppressors of phenotypically distinct synthetic interactions, which share a common altered gene, may prove to be a novel and effective approach for identifying genes whose activities are most directly relevant to the core functions of the shared gene.     

Expression profiling reveals differential gene induction underlying specific and non-specific memory for pheromones in mice

Memory for the mating male’s pheromones in female mice is thought to require synaptic changes in the accessory olfactory bulb (AOB). Induction of this memory depends on release of glutamate in response to pheromonal exposure coincident with release of norepinephrine (NE) in the AOB following mating. A similar memory for pheromones can also be induced artificially by local infusion of the GABA_A receptor antagonist bicuculline into the AOB. The natural memory formed by exposure to pheromones during mating is specific to the pheromones sensed by the female during mating. In contrast, the artificial memory induced by bicuculline is non-specific and results in the female mice recognizing all pheromones as if they were from the mating male. Although protein synthesis has been shown to be essential for development of pheromone memory, the gene expression cascades critical for memory formation are not known. We investigated changes in gene expression in the AOB using oligonucleotide microarrays during mating-induced pheromone memory (MIPM) as well as bicuculline-induced pheromone memory (BIPM). We found the set of genes induced during MIPM and BIPM are largely non-overlapping and Ingenuity Pathway Analysis revealed that the signaling pathways in MIPM and BIPM also differ. The products of genes induced during MIPM are associated with synaptic function, indicating the possibility of modification at specific synapses, while those induced during BIPM appear to possess neuron-wide functions, which would be consistent with global cellular changes. Thus, these results begin to provide a mechanistic explanation for specific and non-specific memories induced by pheromones and bicuculline infusion respectively.     

A filamentous growth response mediated by the yeast mating pathway.

Haploid cells of the budding yeast Saccharomyces cerevisiae respond to mating pheromones by arresting their cell-division cycle in G1 and differentiating into a cell type capable of locating and fusing with mating partners. Yeast cells undergo chemotactic cell surface growth when pheromones are present above a threshold level for morphogenesis; however, the morphogenetic responses of cells to levels of pheromone below this threshold have not been systematically explored. Here we show that MATa haploid cells exposed to low levels of the alpha-factor mating pheromone undergo a novel cellular response: cells modulate their division patterns and cell shape, forming colonies composed of filamentous chains of cells. Time-lapse analysis of filament formation shows that its dynamics are distinct from that of pseudohyphal growth; during pheromone-induced filament formation, daughter cells are delayed relative to mother cells with respect to the timing of bud emergence. Filament formation requires the RSR1(BUD1), BUD8, SLK1/BCK1, and SPA2 genes and many elements of the STE11/STE7 MAP kinase pathway; this response is also independent of FAR1, a gene involved in orienting cell polarization during the mating response. We suggest that mating yeast cells undergo a complex response to low levels of pheromone that may enhance the ability of cells to search for mating partners through the modification of cell shape and alteration of cell-division patterns.     

The yeast G-protein homolog is involved in the mating pheromone signal transduction system.

I have isolated a new type of sterile mutant of Saccharomyces cerevisiae, carrying a single mutant allele, designated dac1, which was mapped near the centromere on chromosome VIII. The dac1 mutation caused specific defects in the pheromone responsiveness of both a and alpha cells and did not seem to be associated with any pleiotropic phenotypes. Thus, in contrast to the ste4, ste5, ste7, ste11, and ste12 mutations, the dac1 mutation had no significant effect on such constitutive functions of haploid cells as pheromone production and alpha-factor destruction. The characteristics of this phenotype suggest that the DAC1 gene encodes a component of the pheromone response pathway common to both a and alpha cells. Introduction of the GPA1 gene encoding an S. cerevisiae homolog of the alpha subunit of mammalian guanine nucleotide-binding regulatory proteins (G proteins) into sterile dac1 mutants resulted in restoration of pheromone responsiveness and mating competence to both a and alpha cells. These results suggest that the dac1 mutation is an allele of the GPA1 gene and thus provide genetic evidence that the yeast G protein homolog is directly involved in the mating pheromone signal transduction pathway.     

Genetic interactions indicate a role for Mdg1p and the SH3 domain protein Bem1p in linking the G-protein mediated yeast pheromone signalling pathway to regulators of cell polarity

The pheromone signal in the yeastSaccharomyces cerevisiae is transmitted by the and subunits of the mating response G-protein. TheSTE20 gene, encoding a protein kinase required for pheromone signal transduction, has recently been identified in a genetic screen for high-gene-dosage suppressors of a partly defective G mutation. The same genetic screen identifiedBEM1, which encodes an SH3 domain protein required for polarized morphogenesis in response to pheromone, and a novel gene, designatedMDG1 (multicopy suppressor ofdefectiveG-protein). TheMDG1 gene was independently isolated in a search for multicopy suppressors of abem1 mutation. TheMDG1 gene encodes a predicted hydrophilic protein of 364 amino acids with a molecular weight of 41 kDa that has no homology with known proteins. A fusion of Mdg1p with the green fluorescent protein fromAequorea victoria localizes to the plasma membrane, suggesting that Mdg1p is an extrinsically bound membrane protein. Deletion ofMDG1 causes sterility in cells in which the wild-type G has been replaced by partly defective G derivatives but does not cause any other obvious phenotypes. The mating defect of cells deleted forSTE20 is partially suppressed by multiple copies ofBEM1 andCDC42, which encodes a small GTP-binding protein that binds to Ste20p and is necessary for the development of cell polarity. Elevated levels ofSTE20 andBEM1 are capable of suppressing a temperature-sensitive mutation inCDC42. This complex network of genetic interactions points to a role for Bem1p and Mdg1p in G-protein mediated signal transduction and indicates a functional linkage between components of the pheromone signalling pathway and regulators of cell polarity during yeast mating.     

Instability of the mitofusin Fzo1 regulates mitochondrial morphology during the mating response of the yeast Saccharomyces cerevisiae

Mitochondria form a highly dynamic network that is shaped by continuous fission and fusion of these organelles. In the yeast Saccharomyces cerevisiae two machineries are involved in this process, one of which includes the mitochondrial fusion promoting GTPase Fzo1. Although a role for the F-box protein Mdm30 in regulating the stability of Fzo1 has been proposed, the molecular basis for the regulation of the fission to fusion ratio of mitochondria remains unknown. To discern the mechanism of the regulation of mitochondrial morphology, we arrested cells at different stages of the cell cycle and examined mitochondrial morphology as well as the stability of mitochondrial fission and fusion proteins. In response to a G1 arrest evoked by the mating pheromone alpha factor the mitochondrial network fragmented into small pieces, which was accompanied by dramatic down-regulation of Fzo1. Mating pheromone also triggered the degradation of Fzo1 produced under the control of a constitutive promoter, and Fzo1 was stabilized upon proteasome inhibition, indicating a role for the proteasome system in the degradation of Fzo1. However, deletion of MDM30 did not stabilize Fzo1 after mating pheromone treatment, showing a different mechanism from the previously reported process of steady state Fzo1 regulation. We show an example for a regulated change of the mitochondrial fission to fusion ratio during the life cycle of budding yeast. Proteasomal degradation of Fzo1 in response to the mating pheromone is proposed to mediate the remodeling of the mitochondrial network during the process of mating.     

Promiscuity resolves constraints on social mate choice imposed by population viscosity

Population viscosity can have major consequences for adaptive evolution, in particular for phenotypes involved in social interactions. For example, population viscosity increases the probability of mating with close kin, resulting in selection for mechanisms that circumvent the potential negative consequences of inbreeding. Female promiscuity is often suggested to be one such mechanism. However, whether avoidance of genetically similar partners is a major selective force shaping patterns of promiscuity remains poorly supported by empirical data. Here, we show (i) that fine‐scale genetic structure constrains social mate choice in a pair‐bonding lizard, resulting in individuals pairing with genetically similar individuals, (ii) that these constraints are circumvented by multiple mating with less related individuals and (iii) that this results in increased heterozygosity of offspring. Despite this, we did not detect any significant effects of heterozygosity on offspring or adult fitness or a strong relationship between pair relatedness and female multiple mating. We discuss these results within the context of incorporating the genetic context dependence of mating strategies into a holistic understanding of mating system evolution.