rcell2: Microscopy-Based Cytometry in R.

This article describes a method for quantifying various cellular features (e.g., volume, curvature, total and sub-cellular fluorescence localization) of individual cells from sets of microscope images, and for tracking them over time-course microscopy experiments. One purposely defocused transmission image (sometimes referred to as bright-field or BF) is used to segment the image and locate each cell. Fluorescence images (one for each of the color channels or z-stacks to be analyzed) may be acquired by conventional wide-field epifluorescence or confocal microscopy. This method uses a set of R packages called rcell2. Relative to the original release of Rcell (Bush et al., 2012), the updated version bundles, into a single software suite, the image-processing capabilities of Cell-ID, offers new data analysis tools for cytometry, and relies on the widely used data analysis and visualization tools of the statistical programming framework R. © 2023 Wiley Periodicals LLC. Basic Protocol: Extracting quantitative information from single cells Support Protocol 1: Obtaining and installing Cell-ID and R Support Protocol 2: Preparing cells for imaging.

Modulation of transcription factor dynamics allows versatile information transmission.

Cells detect changes in their environment and generate responses, often involving changes in gene expression. In this paper we use information theory and a simple transcription model to analyze whether the resulting gene expression serves to identify extracellular stimuli and assess their intensity when they are encoded in the amplitude, duration or frequency of pulses of a transcription factor's nuclear concentration (or activation state). We find, for all cases, that about three ranges of input strengths can be distinguished and that maximum information transmission occurs for fast and high activation threshold promoters. The three input modulation modes differ in the sensitivity to changes in the promoters parameters. Frequency modulation is the most sensitive and duration modulation, the least. This is key for signal identification: there are promoter parameters that yield a relatively high information transmission for duration or amplitude modulation and a much smaller value for frequency modulation. The reverse situation cannot be found with a single promoter transcription model. Thus, pulses of transcription factors can selectively activate the "frequency-tuned" promoter while prolonged nuclear accumulation would activate promoters of all three modes simultaneously. Frequency modulation is therefore highly selective and better suited than the other encoding modes for signal identification without requiring other mediators of the transduction process.

Characterization of cell-to-cell variation in nuclear transport rates and identification of its sources.

Nuclear transport is an essential part of eukaryotic cell function. Here, we present scFRAP, a model-assisted fluorescent recovery after photobleaching (FRAP)- based method to determine nuclear import and export rates independently in individual live cells. To overcome the inherent noise of single-cell measurements, we performed sequential FRAPs on the same cell. We found large cell-to-cell variation in transport rates within isogenic yeast populations. For passive transport, the variability in NPC number might explain most of the variability. Using this approach, we studied mother-daughter cell asymmetry in the active nuclear shuttling of the transcription factor Ace2, which is specifically concentrated in daughter cell nuclei in early G1. Rather than reduced export in the daughter cell, as previously hypothesized, we found that this asymmetry is mainly due to an increased import in daughters. These results shed light on cell-to-cell variation in cellular dynamics and its sources.

Akt Is S-Palmitoylated: A New Layer of Regulation for Akt.

The protein kinase Akt/PKB participates in a great variety of processes, including translation, cell proliferation and survival, as well as malignant transformation and viral infection. In the last few years, novel Akt posttranslational modifications have been found. However, how these modification patterns affect Akt subcellular localization, target specificity and, in general, function is not thoroughly understood. Here, we postulate and experimentally demonstrate by acyl-biotin exchange (ABE) assay and 3H-palmitate metabolic labeling that Akt is S-palmitoylated, a modification related to protein sorting throughout subcellular membranes. Mutating cysteine 344 into serine blocked Akt S-palmitoylation and diminished its phosphorylation at two key sites, T308 and T450. Particularly, we show that palmitoylation-deficient Akt increases its recruitment to cytoplasmic structures that colocalize with lysosomes, a process stimulated during autophagy. Finally, we found that cysteine 344 in Akt1 is important for proper its function, since Akt1-C344S was unable to support adipocyte cell differentiation in vitro. These results add an unexpected new layer to the already complex Akt molecular code, improving our understanding of cell decision-making mechanisms such as cell survival, differentiation and death.

A new tool to sense pH changes at the neuromuscular junction synaptic cleft.

Synaptic transmission triggers transient acidification of the synaptic cleft. Recently, it has been shown that pH affects the opening of postsynaptic channels and therefore the production of tools that allow to study these behaviors should result of paramount value. We fused α-bungarotoxin, a neurotoxin derived from the snake Bungarus multicinctus that binds irreversibly to the acetylcholine receptor extracellular domain, to the pH sensitive GFP Super Ecliptic pHluorin, and efficiently expressed it in Pichia pastoris. This sensor allows synaptic changes in pH to be measured without the need of incorporating transgenes into animal cells. Here, we show that incubation of the mouse levator auris muscle with a solution containing this recombinant protein is enough to fluorescently label the endplate post synaptic membrane. Furthermore, we could physiologically alter and measure post synaptic pH by evaluating changes in the fluorescent signal of pHluorin molecules bound to acetylcholine receptors. In fact, using this tool we were able to detect a drop in 0.01 to 0.05 pH units in the vicinity of the acetylcholine receptors following vesicle exocytosis triggered by nerve electrical stimulation. Further experiments will allow to learn the precise changes in pH during and after synaptic activation.

GPCR receptor phosphorylation and endocytosis are not necessary to switch polarized growth between internal cues during pheromone response in S. cerevisiae.

Chemotactic/chemotropic cells follow accurately the direction of gradients of regulatory molecules. Many G-protein-coupled receptors (GPCR) function as chemoattractant receptors to guide polarized responses. In "a" mating type yeast, the GPCR Ste2 senses the α-cell's pheromone. Previously, phosphorylation and trafficking of this receptor have been implicated in the process of gradient sensing, where cells dynamically correct growth. Correction is often necessary since yeast have intrinsic polarity sites that interfere with a correct initial gradient decoding. We have recently showed that when actively dividing (not in G1) yeast are exposed to a uniform pheromone concentration, they initiate a pheromone-induced polarization next to the mother-daughter cytokinesis site. Then, they reorient their growth to the intrinsic polarity site. Here, to study if Ste2 phosphorylation and internalization are involved in this process, we generated receptor variants combining three types of mutated signals for the first time: phosphorylation, ubiquitylation and the NPFX1,2D Sla1-binding motif. We first characterized their effect on endocytosis and found that these processes regulate internalization in a more complex manner than previously shown. Interestingly, we showed that receptor phosphorylation can drive internalization independently of ubiquitylation and the NPFX1,2D motif. When tested in our assays, cells expressing either phosphorylation or endocytosis-deficient receptors were able to switch away from the cytokinesis site to find the intrinsic polarity site as efficiently as their WT counterparts. Thus, we conclude that these processes are not necessary for the reorientation of polarization.

Mitotic and pheromone-specific intrinsic polarization cues interfere with gradient sensing in Saccharomyces cerevisiae.

Polarity decisions are central to many processes, including mitosis and chemotropism. In Saccharomyces cerevisiae, budding and mating projection (MP) formation use an overlapping system of cortical landmarks that converges on the small G protein Cdc42. However, pheromone-gradient sensing must override the Rsr1-dependent internal polarity cues used for budding. Using this model system, we asked what happens when intrinsic and extrinsic spatial cues are not aligned. Is there competition, or collaboration? By live-cell microscopy and microfluidics techniques, we uncovered three previously overlooked features of this signaling system. First, the cytokinesis-associated polarization patch serves as a polarity landmark independently of all known cues. Second, the Rax1-Rax2 complex functions as a pheromone-promoted polarity cue in the distal pole of the cells. Third, internal cues remain active during pheromone-gradient tracking and can interfere with this process, biasing the location of MPs. Yeast defective in internal-cue utilization align significantly better than wild type with artificially generated pheromone gradients.

Properties of cell signaling pathways and gene expression systems operating far from steady-state.

Ligand-receptor systems, covalent modification cycles, and transcriptional networks are basic units of signaling systems and their steady-state properties are well understood. However, the behavior of such systems before steady-state is poorly characterized. Here, we analyzed the properties of input-output curves for each of these systems as they approach steady-state. In ligand-receptor systems, the EC50 (concentration of the ligand that occupies 50% of the receptors) is higher before the system reaches steady-state. Based on this behavior, we have previously defined PRESS (for pre-equilibrium sensing and signaling), a general "systems level" mechanism cells may use to overcome input saturation. Originally, we showed that, given a step stimulation, PRESS operates when the kinetics of ligand-receptor binding are slower than the downstream signaling steps. Now, we show that, provided the input increases slowly, it is not essential for the ligand binding reaction itself to be slow. In addition, we demonstrate that covalent modification cycles and gene expression systems may also operate in PRESS mode. Thus, nearly all biochemical processes may operate in PRESS mode, suggesting that this mechanism may be ubiquitous in cell signaling systems.

Heat-stress triggers MAPK crosstalk to turn on the hyperosmotic response pathway.

Cells make decisions based on a combination of external and internal signals. In yeast, the high osmolarity response (HOG) is a mitogen-activated protein kinase (MAPK) pathway that responds to a variety of stimuli, and it is central to the general stress response. Here we studied the effect of heat-stress (HS) on HOG. Using live-cell reporters and genetics, we show that HS promotes Hog1 phosphorylation and Hog1-dependent gene expression, exclusively via the Sln1 phosphorelay branch, and that the strength of the activation is larger in yeast adapted to high external osmolarity. HS stimulation of HOG is indirect. First, we show that HS causes glycerol loss, necessary for HOG activation. Preventing glycerol efflux by deleting the glyceroporin FPS1 or its regulators RGC1 and ASK10/RGC2, or by increasing external glycerol, greatly reduced HOG activation. Second, we found that HOG stimulation by HS depended on the operation of a second MAPK pathway, the cell-wall integrity (CWI), a well-known mediator of HS, since inactivating Pkc1 or deleting the MAPK SLT2 greatly reduced HOG activation. Our data suggest that the main role of the CWI in this process is to stimulate glycerol loss. We found that in yeast expressing the constitutively open channel mutant (Fps1-Δ11), HOG activity was independent of Slt2. In summary, we suggest that HS causes a reduction in turgor due to the loss of glycerol and the accompanying water, and that this is what actually stimulates HOG. Thus, taken together, our findings highlight a central role for Fps1, and the metabolism of glycerol, in the communication between the yeast MAPK pathways, essential for survival and reproduction in changing environments.

CDK and MAPK Synergistically Regulate Signaling Dynamics via a Shared Multi-site Phosphorylation Region on the Scaffold Protein Ste5.

We report an unanticipated system of joint regulation by cyclin-dependent kinase (CDK) and mitogen-activated protein kinase (MAPK), involving collaborative multi-site phosphorylation of a single substrate. In budding yeast, the protein Ste5 controls signaling through a G1 arrest pathway. Upon cell-cycle entry, CDK inhibits Ste5 via multiple phosphorylation sites, disrupting its membrane association. Using quantitative time-lapse microscopy, we examined Ste5 membrane recruitment dynamics at different cell-cycle stages. Surprisingly, in S phase, where Ste5 recruitment should be blocked, we observed an initial recruitment followed by a steep drop-off. This delayed inhibition revealed a requirement for both CDK activity and negative feedback from the pathway MAPK Fus3. Mutagenesis, mass spectrometry, and electrophoretic analyses suggest that the CDK and MAPK modify shared sites, which are most extensively phosphorylated when both kinases are active and able to bind their docking sites on Ste5. Such collaborative phosphorylation can broaden regulatory inputs and diversify output dynamics of signaling pathways.

Ultrasensitivity in signaling cascades revisited: Linking local and global ultrasensitivity estimations.

Ultrasensitive response motifs, capable of converting graded stimuli into binary responses, are well-conserved in signal transduction networks. Although it has been shown that a cascade arrangement of multiple ultrasensitive modules can enhance the system's ultrasensitivity, how a given combination of layers affects a cascade's ultrasensitivity remains an open question for the general case. Here, we introduce a methodology that allows us to determine the presence of sequestration effects and to quantify the relative contribution of each module to the overall cascade's ultrasensitivity. The proposed analysis framework provides a natural link between global and local ultrasensitivity descriptors and it is particularly well-suited to characterize and understand mathematical models used to study real biological systems. As a case study, we have considered three mathematical models introduced by O'Shaughnessy et al. to study a tunable synthetic MAPK cascade, and we show how our methodology can help modelers better understand alternative models.

Yeast GPCR signaling reflects the fraction of occupied receptors, not the number.

According to receptor theory, the effect of a ligand depends on the amount of agonist-receptor complex. Therefore, changes in receptor abundance should have quantitative effects. However, the response to pheromone in Saccharomyces cerevisiae is robust (unaltered) to increases or reductions in the abundance of the G-protein-coupled receptor (GPCR), Ste2, responding instead to the fraction of occupied receptor. We found experimentally that this robustness originates during G-protein activation. We developed a complete mathematical model of this step, which suggested the ability to compute fractional occupancy depends on the physical interaction between the inhibitory regulator of G-protein signaling (RGS), Sst2, and the receptor. Accordingly, replacing Sst2 by the heterologous hsRGS4, incapable of interacting with the receptor, abolished robustness. Conversely, forcing hsRGS4:Ste2 interaction restored robustness. Taken together with other results of our work, we conclude that this GPCR pathway computes fractional occupancy because ligand-bound GPCR-RGS complexes stimulate signaling while unoccupied complexes actively inhibit it. In eukaryotes, many RGSs bind to specific GPCRs, suggesting these complexes with opposing activities also detect fraction occupancy by a ratiometric measurement. Such complexes operate as push-pull devices, which we have recently described.

Messages do diffuse faster than messengers: reconciling disparate estimates of the morphogen bicoid diffusion coefficient.

The gradient of Bicoid (Bcd) is key for the establishment of the anterior-posterior axis in Drosophila embryos. The gradient properties are compatible with the SDD model in which Bcd is synthesized at the anterior pole and then diffuses into the embryo and is degraded with a characteristic time. Within this model, the Bcd diffusion coefficient is critical to set the timescale of gradient formation. This coefficient has been measured using two optical techniques, Fluorescence Recovery After Photobleaching (FRAP) and Fluorescence Correlation Spectroscopy (FCS), obtaining estimates in which the FCS value is an order of magnitude larger than the FRAP one. This discrepancy raises the following questions: which estimate is "correct''; what is the reason for the disparity; and can the SDD model explain Bcd gradient formation within the experimentally observed times? In this paper, we use a simple biophysical model in which Bcd diffuses and interacts with binding sites to show that both the FRAP and the FCS estimates may be correct and compatible with the observed timescale of gradient formation. The discrepancy arises from the fact that FCS and FRAP report on different effective (concentration dependent) diffusion coefficients, one of which describes the spreading rate of the individual Bcd molecules (the messengers) and the other one that of their concentration (the message). The latter is the one that is more relevant for the gradient establishment and is compatible with its formation within the experimentally observed times.

Modification of Akt by SUMO conjugation regulates alternative splicing and cell cycle.

Akt/PKB is a key signaling molecule in higher eukaryotes and a crucial protein kinase in human health and disease. Phosphorylation, acetylation, and ubiquitylation have been reported as important regulatory post-translational modifications of this kinase. We describe here that Akt is modified by SUMO conjugation, and show that lysine residues 276 and 301 are the major SUMO attachment sites within this protein. We found that phosphorylation and SUMOylation of Akt appear as independent events. However, decreasing Akt SUMOylation levels severely affects the role of this kinase as a regulator of fibronectin and Bcl-x alternative splicing. Moreover, we observed that the Akt mutant (Akt E17K) found in several human tumors displays increased levels of SUMOylation and also an enhanced capacity to regulate fibronectin splicing patterns. This splicing regulatory activity is completely abolished by decreasing Akt E17K SUMO conjugation levels. Additionally, we found that SUMOylation controls Akt regulatory function at G1/S transition during cell cycle progression. These findings reveal SUMO conjugation as a novel level of regulation for Akt activity, opening new areas of exploration related to the molecular mechanisms involved in the diverse cellular functions of this kinase.

Modulation of the Akt pathway reveals a novel link with PERK/eIF2α, which is relevant during hypoxia.

The unfolded protein response (UPR) and the Akt signaling pathway share several regulatory functions and have the capacity to determine cell outcome under specific conditions. However, both pathways have largely been studied independently. Here, we asked whether the Akt pathway regulates the UPR. To this end, we used a series of chemical compounds that modulate PI3K/Akt pathway and monitored the activity of the three UPR branches: PERK, IRE1 and ATF6. The antiproliferative and antiviral drug Akt-IV strongly and persistently activated all three branches of the UPR. We present evidence that activation of PERK/eIF2α requires Akt and that PERK is a direct Akt target. Chemical activation of this novel Akt/PERK pathway by Akt-IV leads to cell death, which was largely dependent on the presence of PERK and IRE1. Finally, we show that hypoxia-induced activation of eIF2α requires Akt, providing a physiologically relevant condition for the interaction between Akt and the PERK branch of the UPR. These data suggest the UPR and the Akt pathway signal to one another as a means of controlling cell fate.

Pheromone-induced morphogenesis improves osmoadaptation capacity by activating the HOG MAPK pathway.

Environmental and internal conditions expose cells to a multiplicity of stimuli whose consequences are difficult to predict. We investigate the response to mating pheromone of yeast cells adapted to high osmolarity. Events downstream of pheromone binding involve two mitogen-activated protein kinase (MAPK) cascades: the pheromone response (PR) and the cell wall integrity (CWI) response. Although the PR MAPK pathway shares components with a third MAPK pathway, the high osmolarity (HOG) response, each one is normally only activated by its cognate stimulus, a phenomenon called insulation. We found that in cells adapted to high osmolarity, PR activated the HOG pathway in a pheromone- and osmolarity-dependent manner. Activation of HOG by the PR was not due to loss of insulation, but rather a response to a reduction in internal osmolarity, which resulted from an increase in glycerol release caused by the PR. By analyzing single-cell time courses, we found that stimulation of HOG occurred in discrete bursts that coincided with the "shmooing" morphogenetic process. Activation required the polarisome, the CWI MAPK Slt2, and the aquaglyceroporin Fps1. HOG activation resulted in high glycerol turnover, which improved adaptability to rapid changes in osmolarity. Our work shows how a differentiation signal can recruit a second, unrelated sensory pathway to fine-tune yeast response in a complex environment.

Quantitative measurement of protein relocalization in live cells.

Microscope cytometry provides a powerful means to study signaling in live cells. Here we present a quantitative method to measure protein relocalization over time, which reports the absolute fraction of a tagged protein in each compartment. Using this method, we studied an essential step in the early propagation of the pheromone signal in Saccharomyces cerevisiae: recruitment to the membrane of the scaffold Ste5 by activated Gßγ dimers. We found that the dose response of Ste5 recruitment is graded (EC50 = 0.44 ± 0.08 nM, Hill coefficient = 0.8 ± 0.1). Then, we determined the effective dissociation constant (K(de)) between Ste5 and membrane sites during the first few minutes when the negative feedback from the MAPK Fus3 is first activated. K(de) changed during the first minutes from a high affinity of < 0.65 nM to a steady-state value of 17 ± 9 nM. During the same period, the total number of binding sites decreased slightly, from 1940 ± 150 to 1400 ± 200. This work shows how careful quantification of a protein relocalization dynamic can give insight into the regulation mechanisms of a biological system.

Using Cell-ID 1.4 with R for microscope-based cytometry.

This unit describes a method for quantifying various cellular features (e.g., volume, total and subcellular fluorescence localization) from sets of microscope images of individual cells. It includes procedures for tracking cells over time. One purposely defocused transmission image (sometimes referred to as bright-field or BF) is acquired to segment the image and locate each cell. Fluorescence images (one for each of the color channels to be analyzed) are then acquired by conventional wide-field epifluorescence or confocal microscopy. This method uses the image-processing capabilities of Cell-ID and data analysis by the statistical programming framework R, which is supplemented with a package of routines for analyzing Cell-ID output. Both Cell-ID and the analysis package are open-source.

PI3K/AKT pathway regulates phosphorylation of steroid receptors, hormone independence and tumor differentiation in breast cancer.

Using a model of medroxyprogesterone acetate (MPA)-induced mouse mammary tumors that transit through different stages of hormone dependence, we previously reported that the activation of the phosphatidylinositol 3-kinase (PI3K)/AKT (protein kinase B) pathway is critical for the growth of hormone-independent (HI) mammary carcinomas but not for the growth of hormone-dependent (HD) mammary carcinomas. The objective of this work was to explore whether the activation of the PI3K/AKT pathway is responsible for the changes in tumor phenotype and for the transition to autonomous growth. We found that the inhibition of the PI3K/AKT/mTOR (mammalian target of rapamycin) pathway suppresses HI tumor growth. In addition, we were able to induce mammary tumors in mice in the absence of MPA by inoculating HD tumor cells expressing a constitutively active form of AKT1, myristoylated AKT1 (myrAKT1). These tumors were highly differentiated and displayed a ductal phenotype with laminin-1 and cytokeratin 8 expression patterns typical of HI tumors. Furthermore, myrAKT1 increased the tumor growth of IBH-6 and IBH-7 human breast cancer cell lines. In the estrogen-dependent IBH-7 cell line, myrAKT1 induced estrogen-independent growth accompanied by the expression of the adhesion markers focal adhesion kinase and E-cadherin. Finally, we found that cells expressing myrAKT1 exhibited increased phosphorylation of the progesterone receptor at Ser190 and Ser294 and of the estrogen receptor α at Ser118 and Ser167, independently of exogenous MPA or estrogen supply. Our results indicate that the activation of the PI3K/AKT/mTOR pathway promotes tissue architecture remodeling and the activation of steroid receptors.