Single-cell profiling screen identifies microtubule-dependent reduction of variability in signaling.

Populations of isogenic cells often respond coherently to signals, despite differences in protein abundance and cell state. Previously, we uncovered processes in the Saccharomyces cerevisiae pheromone response system (PRS) that reduced cell-to-cell variability in signal strength and cellular response. Here, we screened 1,141 non-essential genes to identify 50 "variability genes". Most had distinct, separable effects on strength and variability of the PRS, defining these quantities as genetically distinct "axes" of system behavior. Three genes affected cytoplasmic microtubule function: BIM1, GIM2, and GIM4 We used genetic and chemical perturbations to show that, without microtubules, PRS output is reduced but variability is unaffected, while, when microtubules are present but their function is perturbed, output is sometimes lowered, but its variability is always high. The increased variability caused by microtubule perturbations required the PRS MAP kinase Fus3 and a process at or upstream of Ste5, the membrane-localized scaffold to which Fus3 must bind to be activated. Visualization of Ste5 localization dynamics demonstrated that perturbing microtubules destabilized Ste5 at the membrane signaling site. The fact that such microtubule perturbations cause aberrant fate and polarity decisions in mammals suggests that microtubule-dependent signal stabilization might also operate throughout metazoans.

Utilization of extracellular information before ligand-receptor binding reaches equilibrium expands and shifts the input dynamic range.

Cell signaling systems sense and respond to ligands that bind cell surface receptors. These systems often respond to changes in the concentration of extracellular ligand more rapidly than the ligand equilibrates with its receptor. We demonstrate, by modeling and experiment, a general "systems level" mechanism cells use to take advantage of the information present in the early signal, before receptor binding reaches a new steady state. This mechanism, pre-equilibrium sensing and signaling (PRESS), operates in signaling systems in which the kinetics of ligand-receptor binding are slower than the downstream signaling steps, and it typically involves transient activation of a downstream step. In the systems where it operates, PRESS expands and shifts the input dynamic range, allowing cells to make different responses to ligand concentrations so high as to be otherwise indistinguishable. Specifically, we show that PRESS applies to the yeast directional polarization in response to pheromone gradients. Consideration of preexisting kinetic data for ligand-receptor interactions suggests that PRESS operates in many cell signaling systems throughout biology. The same mechanism may also operate at other levels in signaling systems in which a slow activation step couples to a faster downstream step.

Phasor plotting with frequency-domain flow cytometry.

Interest in time resolved flow cytometry is growing. In this paper, we collect time-resolved flow cytometry data and use it to create polar plots showing distributions that are a function of measured fluorescence decay rates from individual fluorescently-labeled cells and fluorescent microspheres. Phasor, or polar, graphics are commonly used in fluorescence lifetime imaging microscopy (FLIM). In FLIM measurements, the plotted points on a phasor graph represent the phase-shift and demodulation of the frequency-domain fluorescence signal collected by the imaging system for each image pixel. Here, we take a flow cytometry cell counting system, introduce into it frequency-domain optoelectronics, and process the data so that each point on a phasor plot represents the phase shift and demodulation of an individual cell or particle. In order to demonstrate the value of this technique, we show that phasor graphs can be used to discriminate among populations of (i) fluorescent microspheres, which are labeled with one fluorophore type; (ii) Chinese hamster ovary (CHO) cells labeled with one and two different fluorophore types; and (iii) Saccharomyces cerevisiae cells that express combinations of fluorescent proteins with different fluorescence lifetimes. The resulting phasor plots reveal differences in the fluorescence lifetimes within each sample and provide a distribution from which we can infer the number of cells expressing unique single or dual fluorescence lifetimes. These methods should facilitate analysis time resolved flow cytometry data to reveal complex fluorescence decay kinetics.

Negative feedback that improves information transmission in yeast signalling.

Haploid Saccharomyces cerevisiae yeast cells use a prototypic cell signalling system to transmit information about the extracellular concentration of mating pheromone secreted by potential mating partners. The ability of cells to respond distinguishably to different pheromone concentrations depends on how much information about pheromone concentration the system can transmit. Here we show that the mitogen-activated protein kinase Fus3 mediates fast-acting negative feedback that adjusts the dose response of the downstream system response to match the dose response of receptor-ligand binding. This 'dose-response alignment', defined by a linear relationship between receptor occupancy and downstream response, can improve the fidelity of information transmission by making downstream responses corresponding to different receptor occupancies more distinguishable and reducing amplification of stochastic noise during signal transmission. We also show that one target of the feedback is a previously uncharacterized signal-promoting function of the regulator of G-protein signalling protein Sst2. Our work suggests that negative feedback is a general mechanism used in signalling systems to align dose responses and thereby increase the fidelity of information transmission.

Scaffold number in yeast signaling system sets tradeoff between system output and dynamic range.

Although the proteins comprising many signaling systems are known, less is known about their numbers per cell. Existing measurements often vary by more than 10-fold. Here, we devised improved quantification methods to measure protein abundances in the Saccharomyces cerevisiae pheromone response pathway, an archetypical signaling system. These methods limited variation between independent measurements of protein abundance to a factor of two. We used these measurements together with quantitative models to identify and investigate behaviors of the pheromone response system sensitive to precise abundances. The difference between the maximum and basal signaling output (dynamic range) of the pheromone response MAPK cascade was strongly sensitive to the abundance of Ste5, the MAPK scaffold protein, and absolute system output depended on the amount of Fus3, the MAPK. Additional analysis and experiment suggest that scaffold abundance sets a tradeoff between maximum system output and system dynamic range, a prediction supported by recent experiments.

Autorepression-Based Conditional Gene Expression System in Yeast for Variation-Suppressed Control of Protein Dosage.

Conditional control of gene expression allows an experimenter to investigate many aspects of a gene's function. In the model organism Saccharomyces cerevisiae, a number of methods to control gene expression are widely practiced, including induction by metabolites, small molecules, and even light. However, all current methods suffer from at least one of a set of drawbacks, including need for specialized growth conditions, leaky expression, or requirement of specialized equipment. Here we describe protocols using two transformations to construct strains that carry a new controller in which all these drawbacks are overcome. In these strains, the expression of a controlled gene of interest is repressed by the bacterial repressor TetR and induced by anhydrotetracycline. TetR also regulates its own expression, creating an autorepression loop. This autorepression allows tight control of gene expression and protein dosage with low cell-to-cell variation in expression. A second repressor, TetR-Tup1, prevents any leaky expression. We also present a protocol showing a particular workhorse application of such strains to generate synchronized cell populations. We turn off expression of the cell cycle regulator CDC20 completely, arresting the cell population, and then we turn it back on so that the synchronized cells resume cell cycle progression. This control system can be applied to any endogenous or exogenous gene for precise expression. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Generating a parent WTC846 strain Basic Protocol 2: Generating a WTC846 strain with controlled expression of the targeted gene Alternate Protocol: CRISPR-mediated promoter replacement Basic Protocol 3: Cell cycle synchronization/arrest and release using the WTC846- K3 ::CDC20 strain.

Detailed simulations of cell biology with Smoldyn 2.1.

Most cellular processes depend on intracellular locations and random collisions of individual protein molecules. To model these processes, we developed algorithms to simulate the diffusion, membrane interactions, and reactions of individual molecules, and implemented these in the Smoldyn program. Compared to the popular MCell and ChemCell simulators, we found that Smoldyn was in many cases more accurate, more computationally efficient, and easier to use. Using Smoldyn, we modeled pheromone response system signaling among yeast cells of opposite mating type. This model showed that secreted Bar1 protease might help a cell identify the fittest mating partner by sharpening the pheromone concentration gradient. This model involved about 200,000 protein molecules, about 7000 cubic microns of volume, and about 75 minutes of simulated time; it took about 10 hours to run. Over the next several years, as faster computers become available, Smoldyn will allow researchers to model and explore systems the size of entire bacterial and smaller eukaryotic cells.

Regulated cell-to-cell variation in a cell-fate decision system.

Here we studied the quantitative behaviour and cell-to-cell variability of a prototypical eukaryotic cell-fate decision system, the mating pheromone response pathway in yeast. We dissected and measured sources of variation in system output, analysing thousands of individual, genetically identical cells. Only a small proportion of total cell-to-cell variation is caused by random fluctuations in gene transcription and translation during the response ('expression noise'). Instead, variation is dominated by differences in the capacity of individual cells to transmit signals through the pathway ('pathway capacity') and to express proteins from genes ('expression capacity'). Cells with high expression capacity express proteins at a higher rate and increase in volume more rapidly. Our results identify two mechanisms that regulate cell-to-cell variation in pathway capacity. First, the MAP kinase Fus3 suppresses variation at high pheromone levels, while the MAP kinase Kss1 enhances variation at low pheromone levels. Second, pathway capacity and expression capacity are negatively correlated, suggesting a compensatory mechanism that allows cells to respond more precisely to pheromone in the presence of a large variation in expression capacity.

A precisely adjustable, variation-suppressed eukaryotic transcriptional controller to enable genetic discovery.

Conditional expression of genes and observation of phenotype remain central to biological discovery. Current methods enable either on/off or imprecisely controlled graded gene expression. We developed a 'well-tempered' controller, WTC846, for precisely adjustable, graded, growth condition independent expression of genes in Saccharomyces cerevisiae. Controlled genes are expressed from a strong semisynthetic promoter repressed by the prokaryotic TetR, which also represses its own synthesis; with basal expression abolished by a second, 'zeroing' repressor. The autorepression loop lowers cell-to-cell variation while enabling precise adjustment of protein expression by a chemical inducer. WTC846 allelic strains in which the controller replaced the native promoters recapitulated known null phenotypes (CDC42, TPI1), exhibited novel overexpression phenotypes (IPL1), showed protein dosage-dependent growth rates and morphological phenotypes (CDC28, TOR2, PMA1 and the hitherto uncharacterized PBR1), and enabled cell cycle synchronization (CDC20). WTC846 defines an 'expression clamp' allowing protein dosage to be adjusted by the experimenter across the range of cellular protein abundances, with limited variation around the setpoint.

Quantifying small numbers of antibodies with a 'near-universal' protein-DNA chimera.

We present general means to greatly increase the sensitivity of antibody-based assays. Augmentation relies on a 'tadpole' protein-DNA chimera whose protein moiety binds most classes of mammalian antibodies but not avian immunoglobulin Y (IgY). We used this tadpole in affinity capture assays followed by real-time PCR to quantify numerous molecules, including prostate-specific antigen (PSA) in human serum, with great sensitivity and accuracy.

Recipes and Tools for Culture of Escherichia coli.

In this article, we provide information about culture media, including minimal liquid media, rich liquid media, solid media, top agar, and stab agar. We also provide descriptions and useful information about tools used with growth media such as inoculating loops, sterile toothpicks, and spreaders. © 2018 by John Wiley & Sons, Inc.

Growth of E. coli in Liquid Medium.

We describe the procedure for inoculating overnight (starter) cultures of E. coli from a single colony, along with considerations for growing larger cultures. We also include two methods for monitoring the number of cells per unit volume (density) of liquid cultures using a spectrophotometer and a hemacytometer or "count slide." © 2018 by John Wiley & Sons, Inc.

Growth of E. coli on Solid Media.

We provide protocols for titering and isolating bacterial colonies from single cells by serial dilutions, for streaking agar plates, and for spreading suspensions of cells on plates. Support protocols describe replica plating and methods for storing strains as agar stabs and frozen stocks. © 2018 by John Wiley & Sons, Inc.

In vivo measurements reveal a single 5'-intron is sufficient to increase protein expression level in Caenorhabditis elegans.

Introns can increase gene expression levels using a variety of mechanisms collectively referred to as Intron Mediated Enhancement (IME). IME has been measured in cell culture and plant models by quantifying expression of intronless and intron-bearing reporter genes in vitro. We developed hardware and software to implement microfluidic chip-based gene expression quantification in vivo. We altered position, number and sequence of introns in reporter genes controlled by the hsp-90 promoter. Consistent with plant and mammalian studies, we determined a single, natural or synthetic, 5'-intron is sufficient for the full IME effect conferred by three synthetic introns, while a 3'-intron is not. We found coding sequence can affect IME; the same three synthetic introns that increase mcherry protein concentration by approximately 50%, increase mEGFP by 80%. We determined IME effect size is not greatly affected by the stronger vit-2 promoter. Our microfluidic imaging approach should facilitate screens for factors affecting IME and other intron-dependent processes.

Chaperone biomarkers of lifespan and penetrance track the dosages of many other proteins.

Many traits vary among isogenic individuals in homogeneous environments. In microbes, plants and animals, variation in the protein chaperone system affects many such traits. In the animal model C. elegans, the expression level of hsp-16.2 chaperone biomarkers correlates with or predicts the penetrance of mutations and lifespan after heat shock. But the physiological mechanisms causing cells to express different amounts of the biomarker were unknown. Here, we used an in vivo microscopy approach to dissect different contributions to cell-to-cell variation in hsp-16.2 expression in the intestines of young adult animals, which generate the most lifespan predicting signal. While we detected both cell autonomous intrinsic noise and signaling noise, we found both contributions were relatively unimportant. The major contributor to cell-to-cell variation in biomarker expression was general differences in protein dosage. The hsp-16.2 biomarker reveals states of high or low effective dosage for many genes.

Automatic generation of cellular reaction networks with Moleculizer 1.0.

Accurate simulation of intracellular biochemical networks is essential to furthering our understanding of biological system behavior. The number of protein complexes and of chemical interactions among them has traditionally posed significant problems for simulation algorithms. Here we describe an approach to the exact stochastic simulation of biochemical networks that emphasizes the contribution of protein complexes to these systems. This simulation approach starts from a description of monomeric proteins and specifications for binding, unbinding and other reactions. This manageable specification is reasonably intuitive for biologists. Rather than requiring the inclusion of all possible complexes and reactions from the outset, our approach incorporates new complexes and reactions only when needed as the simulation proceeds. As a result, the simulation generates much smaller reaction networks, which can be exported to other simulators for further analysis. We apply this approach to the automatic generation of reaction systems for the study of signal transduction networks.

Transferring information without distortion.

Despite employing diverse molecular mechanisms, many different cell signaling systems avoid losing information by transmitting it in a linear manner.

Caenorhabditis elegans Genes Affecting Interindividual Variation in Life-span Biomarker Gene Expression.

Genetically identical organisms grown in homogenous environments differ in quantitative phenotypes. Differences in one such trait, expression of a single biomarker gene, can identify isogenic cells or organisms that later manifest different fates. For example, in isogenic populations of young adult Caenorhabditis elegans, differences in Green Fluorescent Protein (GFP) expressed from the hsp-16.2 promoter predict differences in life span. Thus, it is of interest to determine how interindividual differences in biomarker gene expression arise. Prior reports showed that the thermosensory neurons and insulin signaling systems controlled the magnitude of the heat shock response, including absolute expression of hsp-16.2. Here, we tested whether these regulatory signals might also influence variation in hsp-16.2 reporter expression. Genetic experiments showed that the action of AFD thermosensory neurons increases interindividual variation in biomarker expression. Further genetic experimentation showed the insulin signaling system acts to decrease interindividual variation in life-span biomarker expression; in other words, insulin signaling canalizes expression of the hsp-16.2-driven life-span biomarker. Our results show that specific signaling systems regulate not only expression level, but also the amount of interindividual expression variation for a life-span biomarker gene. They raise the possibility that manipulation of these systems might offer means to reduce heterogeneity in the aging process.