Single-cell Time-lapse Imaging Research Articles

Tools and methods for high-throughput single-cell imaging with the mother machine

Despite much progress, image processing remains a significant bottleneck for high-throughput analysis of microscopy data. One popular platform for single-cell time-lapse imaging is the mother machine, which enables long-term tracking of microbial cells under precisely controlled growth conditions. While several mother machine image analysis pipelines have been developed in the past several years, adoption by a non-expert audience remains a challenge. To fill this gap, we implemented our own software, MM3, as a plugin for the multidimensional image viewer napari. napari-MM3 is a complete and modular image analysis pipeline for mother machine data, which takes advantage of the high-level interactivity of napari. Here, we give an overview of napari-MM3 and test it against several well-designed and widely used image analysis pipelines, including BACMMAN and DeLTA. Researchers often analyze mother machine data with custom scripts using varied image analysis methods, but a quantitative comparison of the output of different pipelines has been lacking. To this end, we show that key single-cell physiological parameter correlations and distributions are robust to the choice of analysis method. However, we also find that small changes in thresholding parameters can systematically alter parameters extracted from single-cell imaging experiments. Moreover, we explicitly show that in deep learning-based segmentation, ‘what you put is what you get’ (WYPIWYG) – that is, pixel-level variation in training data for cell segmentation can propagate to the model output and bias spatial and temporal measurements. Finally, while the primary purpose of this work is to introduce the image analysis software that we have developed over the last decade in our lab, we also provide information for those who want to implement mother machine-based high-throughput imaging and analysis methods in their research.

Tools and methods for high-throughput single-cell imaging with the mother machine.

Despite much progress, image processing remains a significant bottleneck for high-throughput analysis of microscopy data. One popular platform for single-cell time-lapse imaging is the mother machine, which enables long-term tracking of microbial cells under precisely controlled growth conditions. While several mother machine image analysis pipelines have been developed in the past several years, adoption by a non-expert audience remains a challenge. To fill this gap, we implemented our own software, MM3, as a plugin for the multidimensional image viewer napari. napari-MM3 is a complete and modular image analysis pipeline for mother machine data, which takes advantage of the high-level interactivity of napari. Here, we give an overview of napari-MM3 and test it against several well-designed and widely used image analysis pipelines, including BACMMAN and DeLTA. Researchers often analyze mother machine data with custom scripts using varied image analysis methods, but a quantitative comparison of the output of different pipelines has been lacking. To this end, we show that key single-cell physiological parameter correlations and distributions are robust to the choice of analysis method. However, we also find that small changes in thresholding parameters can systematically alter parameters extracted from single-cell imaging experiments. Moreover, we explicitly show that in deep learning-based segmentation, 'what you put is what you get' (WYPIWYG) - that is, pixel-level variation in training data for cell segmentation can propagate to the model output and bias spatial and temporal measurements. Finally, while the primary purpose of this work is to introduce the image analysis software that we have developed over the last decade in our lab, we also provide information for those who want to implement mother machine-based high-throughput imaging and analysis methods in their research.

The oscillation of mitotic kinase governs cell cycle latches in mammalian cells.

The mammalian cell cycle alternates between two phases - S-G2-M with high levels of A- and B-type cyclins (CycA and CycB, respectively) bound to cyclin-dependent kinases (CDKs), and G1 with persistent degradation of CycA and CycB by an activated anaphase promoting complex/cyclosome (APC/C) bound to Cdh1 (also known as FZR1 in mammals; denoted APC/C:Cdh1). Because CDKs phosphorylate and inactivate Cdh1, these two phases are mutually exclusive. This 'toggle switch' is flipped from G1 to S by cyclin-E bound to a CDK (CycE:CDK), which is not degraded by APC/C:Cdh1, and from M to G1 by Cdc20-bound APC/C (APC/C:Cdc20), which is not inactivated by CycA:CDK or CycB:CDK. After flipping the switch, cyclin E is degraded and APC/C:Cdc20 is inactivated. Combining mathematical modelling with single-cell timelapse imaging, we show that dysregulation of CycB:CDK disrupts strict alternation of the G1-S and M-G1 switches. Inhibition of CycB:CDK results in Cdc20-independent Cdh1 'endocycles', and sustained activity of CycB:CDK drives Cdh1-independent Cdc20 endocycles. Our model provides a mechanistic explanation for how whole-genome doubling can arise, a common event in tumorigenesis that can drive tumour evolution.

The intensities of canonical senescence biomarkers integrate the duration of cell-cycle withdrawal

Senescence, a state of irreversible cell-cycle withdrawal, is difficult to distinguish from quiescence, a state of reversible cell-cycle withdrawal. This difficulty arises because quiescent and senescent cells are defined by overlapping biomarkers, raising the question of whether these states are truly distinct. To address this, we use single-cell time-lapse imaging to distinguish slow-cycling cells that spend long periods in quiescence from cells that never cycle after recovery from senescence-inducing treatments, followed by staining for various senescence biomarkers. We find that the staining intensity of multiple senescence biomarkers is graded rather than binary and reflects the duration of cell-cycle withdrawal, rather than senescence per se. Together, our data show that quiescent and apparent senescent cells are nearly molecularly indistinguishable from each other at a snapshot in time. This suggests that cell-cycle withdrawal itself is graded rather than binary, where the intensities of senescence biomarkers integrate the duration of past cell-cycle withdrawal.

Cyclin E/CDK2 and feedback from soluble histone protein regulate the S phase burst of histone biosynthesis

Faithful DNA replication requires that cells fine-tune their histone pool in coordination with cell-cycle progression. Replication-dependent histone biosynthesis is initiated at a low level upon cell-cycle commitment, followed by a burst at the G1/S transition, but it remains unclear how exactly the cell regulates this burst in histone biosynthesis as DNA replication begins. Here, we use single-cell time-lapse imaging to elucidate the mechanisms by which cells modulate histone production during different phases of the cell cycle. We find that CDK2-mediated phosphorylation of NPAT at the restriction point triggers histone transcription, which results in a burst of histone mRNA precisely at the G1/S phase boundary. Excess soluble histone protein further modulates histone abundance by promoting the degradation of histone mRNA for the duration of S phase. Thus, cells regulate their histone production in strict coordination with cell-cycle progression by two distinct mechanisms acting in concert.

Multiparameter analysis of timelapse imaging reveals kinetics of megakaryocytic erythroid progenitor clonal expansion and differentiation

Single-cell assays have enriched our understanding of hematopoiesis and, more generally, stem and progenitor cell biology. However, these single-end-point approaches provide only a static snapshot of the state of a cell. To observe and measure dynamic changes that may instruct cell fate, we developed an approach for examining hematopoietic progenitor fate specification using long-term (> 7-day) single-cell time-lapse imaging for up to 13 generations with in situ fluorescence staining of primary human hematopoietic progenitors followed by algorithm-assisted lineage tracing. We analyzed progenitor cell dynamics, including the division rate, velocity, viability, and probability of lineage commitment at the single-cell level over time. We applied a Markov probabilistic model to predict progenitor division outcome over each generation in culture. We demonstrated the utility of this methodological pipeline by evaluating the effects of the cytokines thrombopoietin and erythropoietin on the dynamics of self-renewal and lineage specification in primary human bipotent megakaryocytic-erythroid progenitors (MEPs). Our data support the hypothesis that thrombopoietin and erythropoietin support the viability and self-renewal of MEPs, but do not affect fate specification. Thus, single-cell tracking of time-lapse imaged colony-forming unit assays provides a robust method for assessing the dynamics of progenitor self-renewal and lineage commitment.

YB-1 Knockdown Inhibits the Proliferation of Mesothelioma Cells through Multiple Mechanisms.

Y-box binding protein-1 (YB-1) is a multifunctional oncoprotein that has been shown to regulate proliferation, invasion and metastasis in a variety of cancer types. We previously demonstrated that YB-1 is overexpressed in mesothelioma cells and its knockdown significantly reduces tumour cell proliferation, migration, and invasion. However, the mechanisms driving these effects are unclear. Here, we utilised an unbiased RNA-seq approach to characterise the changes to gene expression caused by loss of YB-1 knockdown in three mesothelioma cell lines (MSTO-211H, VMC23 and REN cells). Bioinformatic analysis showed that YB-1 knockdown regulated 150 common genes that were enriched for regulators of mitosis, integrins and extracellular matrix organisation. However, each cell line also displayed unique gene expression signatures, that were differentially enriched for cell death or cell cycle control. Interestingly, deregulation of STAT3 and p53-pathways were a key differential between each cell line. Using flow cytometry, apoptosis assays and single-cell time-lapse imaging, we confirmed that MSTO-211H, VMC23 and REN cells underwent either increased cell death, G1 arrest or aberrant mitotic division, respectively. In conclusion, this data indicates that YB-1 knockdown affects a core set of genes in mesothelioma cells. Loss of YB-1 causes a cascade of events that leads to reduced mesothelioma proliferation, dependent on the underlying functionality of the STAT3/p53-pathways and the genetic landscape of the cell.

Dynamic subcellular translocation of V-type H+ -ATPase is essential for biomineralization of the diatom silica cell wall.

Diatom cell walls, called frustules, are main sources of biogenic silica in the ocean and their intricate morphology is an inspiration for nanoengineering. Here we show dynamic aspects of frustule biosynthesis involving acidification of the silica deposition vesicle (SDV) by V-type H+ ATPase (VHA). Transgenic Thalassiosira pseudonana expressing the VHA B subunit tagged with enhanced green fluorescent protein (VHAB -eGFP) enabled subcellular protein localization in live cells. In exponentially growing cultures, VHAB -eGFP was present in various subcellular localizations including the cytoplasm, SDVs and vacuoles. We studied the role of VHA during frustule biosynthesis in synchronized cell cultures of T.pseudonana. During the making of new biosilica components, VHAB -eGFP first localized in the girdle band SDVs, and subsequently in valve SDVs. In single cell time-lapse imaging experiments, VHAB -eGFP localization in SDVs precluded accumulation of the acidotropic silica biomineralization marker PDMPO. Furthermore, pharmacological VHA inhibition prevented PDMPO accumulation in the SDV, frustule biosynthesis and cell division, as well as insertion of the silicalemma-associated protein SAP1 into the SDVs. Finally, partial inhibition of VHA activity affected the nanoscale morphology of the valve. Altogether, these results indicate that VHA is essential for frustule biosynthesis by acidifying the SDVs and regulating the insertion of other structural proteins into the SDV.

Single-cell kinetics of siRNA-mediated mRNA degradation

RNA interference (RNAi) enables the therapeutic use of small interfering RNAs (siRNAs) to silence disease-related genes. The efficiency of silencing is commonly assessed by measuring expression levels of the target protein at a given time point post-transfection. Here, we determine the siRNA-induced fold change in mRNA degradation kinetics from single-cell fluorescence time-courses obtained using live-cell imaging on single-cell arrays (LISCA). After simultaneous transfection of mRNAs encoding eGFP (target) and CayRFP (reference), the eGFP expression is silenced by siRNA. The single-cell time-courses are fitted using a mathematical model of gene expression. Analysis yields best estimates of related kinetic rate constants, including mRNA degradation constants. We determine the siRNA-induced changes in kinetic rates and their correlations between target and reference protein expression. Assessment of mRNA degradation constants using single-cell time-lapse imaging is fast (<30 h) and returns an accurate, time-independent measure of siRNA-induced silencing, thus allowing the exact evaluation of siRNA therapeutics.

Investigating Age-Related Changes in Proliferation and the Cell Division Repertoire of Parenchymal Reactive Astrocytes.

Reactive gliosis is a complicated process involving all types of glial cells and is the therapeutic target of efforts to treat several types of neuropathologies. Parenchymal astrocytes continuously survey their microenvironment to identify even tiny abnormalities in the central nervous system (CNS) homeostasis and react rapidly to brain damage, such as following ischemia, trauma, or neurodegenerative diseases, to prevent propagation of tissue damage. Aging can play causal roles in certain astroglial dysfunctions, however, still little is known to what extent the heterogeneous reaction of astrocytes at the injury site might be impaired over the course of aging. Based on our experience with both in vitro and in vivo experimental paradigms, we describe here in detail the analysis of age-related changes in (1) proliferative response of parenchymal astrocytes within the posttraumatic cerebral cortex grey matter (GM), and (2) repertoire of their cell divisions in adherent cell culture prepared from the injured GM of young and old double transgenic GFAP-mRFP1/(FUCCI)-S/G2/M-mAG-hGeminin mice by single cell time-lapse imaging.

Random migration of induced pluripotent stem cell-derived human gastrulation-stage mesendoderm.

Gastrulation is the initial systematic deformation of the embryo to form germ layers, which is characterized by the placement of appropriate cells in their destined locations. Thus, gastrulation, which occurs at the beginning of the second month of pregnancy, is a critical stage in human body formation. Although histological analyses indicate that human gastrulation is similar to that of other amniotes (birds and mammals), much of human gastrulation dynamics remain unresolved due to ethical and technical limitations. We used human induced pluripotent stem cells (hiPSCs) to study the migration of mesendodermal cells through the primitive streak to form discoidal germ layers during gastrulation. Immunostaining results showed that hiPSCs differentiated into mesendodermal cells and that epithelial–mesenchymal transition occurred through the activation of the Activin/Nodal and Wnt/beta-catenin pathways. Single-cell time-lapse imaging of cells adhered to cover glass showed that mesendodermal differentiation resulted in the dissociation of cells and an increase in their migration speed, thus confirming the occurrence of epithelial–mesenchymal transition. These results suggest that mesendodermal cells derived from hiPSCs may be used as a model system for studying migration during human gastrulation in vitro. Using random walk analysis, we found that random migration occurred for both undifferentiated hiPSCs and differentiated mesendodermal cells. Two-dimensional random walk simulation showed that homogeneous dissociation of particles may form a discoidal layer, suggesting that random migration might be suitable to effectively disperse cells homogeneously from the primitive streak to form discoidal germ layers during human gastrulation.

Feedback-mediated signal conversion promotes viral fitness

A fundamental signal-processing problem is how biological systems maintain phenotypic states (i.e., canalization) long after degradation of initial catalyst signals. For example, to efficiently replicate, herpesviruses (e.g., human cytomegalovirus, HCMV) rapidly counteract cell-mediated silencing using transactivators packaged in the tegument of the infecting virion particle. However, the activity of these tegument transactivators is inherently transient-they undergo immediate proteolysis but delayed synthesis-and how transient activation sustains lytic viral gene expression despite cell-mediated silencing is unclear. By constructing a two-color, conditional-feedback HCMV mutant, we find that positive feedback in HCMV's immediate-early 1 (IE1) protein is of sufficient strength to sustain HCMV lytic expression. Single-cell time-lapse imaging and mathematical modeling show that IE1 positive feedback converts transient transactivation signals from tegument pp71 proteins into sustained lytic expression, which is obligate for efficient viral replication, whereas attenuating feedback decreases fitness by promoting a reversible silenced state. Together, these results identify a regulatory mechanism enabling herpesviruses to sustain expression despite transient activation signals-akin to early electronic transistors-and expose a potential target for therapeutic intervention.

Quantification and modeling of signalling dynamics in single hematopoietic stem and progenitor cells

Signalling pathways influence fate decisions of hematopoietic cells. However, quantifying signalling dynamics over time and linking these to future fate decisions of the same individual cells have been technically infeasible. Here, we establish continuous single-cell time-lapse imaging to quantify the activity of various signalling pathways over time in living hematopoietic stem and progenitor cells (HSPCs). Live signalling biosensors reporting the activity of the ERK, PI3K/AKT, STAT3, WNT and NOTCH signalling pathways are combined to generate a transgenic mouse line simultaneously displaying their activity. Isolating known single cells after live observation allows to link signalling dynamics to live functional assays or to destructive high-dimensional RNA or protein quantification. We expect to identify novel HSPC sub-populations with specific signalling responses, and a better understanding of the relevance of signalling pathways in controlling transcriptional programs and HSPC fates, of molecular HSC control in general.

Single-cell time-lapse imaging of intracellular O2 in response to metabolic inhibition and mitochondrial cytochrome-c release

The detection of intracellular molecular oxygen (O2) levels is important for understanding cell physiology, cell death, and drug effects, and has recently been improved with the development of oxygen-sensitive probes that are compatible with live cell time-lapse microscopy. We here provide a protocol for the use of the nanoparticle probe MitoImage-MM2 to monitor intracellular oxygen levels by confocal microscopy under baseline conditions, in response to mitochondrial toxins, and following mitochondrial cytochrome-c release. We demonstrate that the MitoImage-MM2 probe, which embeds Pt(II)-5,10,15,20-tetrakis-(2,3,4,5,6–pentafluorophenyl)-porphyrin as oxygen sensor and poly(9,9-dioctylfluorene) as an O2-independent component, enables quantitative, ratiometric time-lapse imaging of intracellular O2. Multiplexing with tetra-methyl-rhodamine-methyl ester in HeLa cervical cancer cells showed significant increases in intracellular O2 accompanied by strong mitochondrial depolarization when respiratory chain complexes III or IV were inhibited by Antimycin A or sodium azide, respectively, and when cells were maintained at ‘physiological’ tissue O2 levels (5% O2). Multiplexing also allowed us to monitor intracellular O2 during the apoptotic signaling process of mitochondrial outer membrane permeabilization in HeLa expressing cytochrome-c-eGFP, and demonstrated that mitochondria post cytochrome-c release are able to retain their capacity to respire at physiological O2 despite a decrease in mitochondrial membrane potential.

Endogenous Replication Stress in Mother Cells Leads to Quiescence of Daughter Cells

Mammalian cells have two fundamentally different states, proliferative and quiescent, but our understanding of how and why cells switch between these states is limited. We previously showed that actively proliferating populations contain a subpopulation that enters quiescence (G0) in an apparently stochastic manner. Using single-cell time-lapse imaging of CDK2 activity and DNA damage, we now show that unresolved endogenous replication stress in the previous (mother) cell cycle prompts p21-dependent entry of daughter cells into quiescence immediately after mitosis. Furthermore, the amount of time daughter cells spend in quiescence is correlated with the extent of inherited damage. Our study thus links replication errors in one cell cycle to the fate of daughter cells in the subsequent cell cycle. More broadly, this work reveals that entry into quiescence is not purely stochastic but has a strong deterministic component arising from a memory of events that occurred in the previous generation(s).

Abstract 4224: Multiplexing single cell time lapse imaging of intra-cellular oxygen in cancer cells

Abstract Oxygen plays an important role in aerobic energy metabolism and signal transduction. The detection of intracellular or subcellular molecular oxygen (O2) levels is important for the understanding of cell physiology and mitochondrial energy metabolism, and the alterations correlated with cancer. Detection of O2 has recently been improved significantly with the development of live cell time lapse microscopy compatible oxygen sensitive probes. We here utilized nanoparticle probe MitoImageTM-MM2 with embedded poly(9,9-dioctylfluorene) (PFO) as an O2-independent component and Pt(II)-5,10,15,20-tetrakis-(2,3,4,5,6-pentafluorophenyl)-porphyrin (PtTFPP) to perform time lapse ratiometric measurements of intra-cellular oxygen levels. Loading of HeLa cervical cancer cells with MM2 enables us to image kinetics of intracellular O2 concentration in response to the mitochondrial respiratory chain inhibitors Antimycin A (complex III) or sodium azide (complex IV). These changes were detected in respiration-supporting medium equilibrated with 5% (53 μM) and 20% (225 μM)O2. Also exposing cells to varying extracellular oxygen concentrations (21 - 225 μM) showed strong dependence of the intracellular oxygen levels. The mitochondrial membrane potential was monitored by simultaneous imaging of Tetra-Methyl-Rhodamine-Methyl ester (TMRM) in the same sample. Multiplexing with GFP finally allowed us to monitor intracellular O2 during the apoptotic signaling process of mitochondrial outer membrane permeabilisation (MOMP) in HeLa expressing cytochrome-c-eGFP and stained with TMRM. In conclusion, the MitoImageTM-MM2 nanoparticle probe enables single-cell, time-lapse imaging of alterations in intracellular oxygen levels in cancer cells, and can be employed in multiplex assays to investigate alterations in intracellular oxygen levels relative to protein dynamics and mitochondrial function. Keywords: Intracellular oxygen; phosphorescent oxygen-sensitive nanoparticles; HeLa cells; MitoImage ratiometric probes; confocal microscopy. Citation Format: Heiko Düssmann, Sergio Perez, Ujval Anil Kumar, Dmitri B. Papkovsky, Jochen HM Prehn. Multiplexing single cell time lapse imaging of intra-cellular oxygen in cancer cells. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 4224.

Combinatorial gene regulation by modulation of relative pulse timing

Studies of individual living cells have revealed that many transcription factors activate in dynamic, and often stochastic, pulses within the same cell. However, it has remained unclear whether cells might modulate the relative timing of these pulses to control gene expression. Here, using quantitative single-cell time-lapse imaging of Saccharomyces cerevisiae, we show that the pulsatile transcription factors Msn2 and Mig1 combinatorially regulate their target genes through modulation of their relative pulse timing. The activator Msn2 and repressor Mig1 pulsed in either a temporally overlapping or non-overlapping manner during their transient response to different inputs, with only the non-overlapping dynamics efficiently activating target gene expression. Similarly, under constant environmental conditions, where Msn2 and Mig1 exhibit sporadic pulsing, glucose concentration modulated the temporal overlap between pulses of the two factors. Together, these results reveal a time-based mode of combinatorial gene regulation. Regulation through relative signal timing is common in engineering and neurobiology, and these results suggest that it could also function broadly within the signaling and regulatory systems of the cell.

Imaging Pluripotency: Time-Lapse Analysis of Mouse Embryonic Stem Cells.

The current view of the pluripotent state is that of a transient, dynamic state, maintained by the balance between opposing cues. Understanding how this dynamic state is established in pluripotent cells and how it relates to gene expression is essential to obtain a more detailed description of the pluripotent state.In this chapter, we describe how to study the dynamic expression of a core pluripotency gene regulator-Nanog-by exploiting single-cell time-lapse imaging of a reporter mESC line grown in different cell culture media. We further describe an automated image analysis method and discuss how to extract information from the generated quantitative time-course data.

Proteasome Inhibition Can Impair Caspase-8 Activation upon Submaximal Stimulation of Apoptotic Tumor Necrosis Factor-related Apoptosis Inducing Ligand (TRAIL) Signaling

Tumor necrosis factor-related apoptosis inducing ligand (TRAIL) can induce extrinsic apoptosis, resulting in caspase-8 activation, but may also initiate transcription-dependent prosurvival signaling. Proteasome inhibitors were suggested to promote TRAIL signal transduction through the death-inducing signaling complex (DISC) by modulating the relative abundance of core DISC components, thereby enhancing caspase-8 activation and apoptosis. To test this hypothesis, we quantified the changes in DISC protein levels as an early consequence of proteasome inhibition in HeLa cervical cancer cells and, based on these data, mathematically modeled the proapoptotic TRAIL signaling toward caspase-8 activation. Modeling results surprisingly suggested that caspase-8 activation might be delayed in presence of proteasome inhibitors, in particular at submaximal TRAIL doses. Subsequent FRET-based single cell time-lapse imaging at conditions where transcription dependent prosurvival signaling was blocked confirmed this hypothesis: caspase-8 activity was delayed by hours in the presence of proteasome inhibitors epoxomicin or bortezomib. Corresponding delays were detected for effector caspase processing and cell death. Contrary to current models, we therefore provide evidence that synergies between TRAIL and proteasome inhibitors do not result from changes in the levels of core DISC signaling proteins.

Understanding the Mechanotransductional Basis of Intravascular Bubble Injury

Intravascular bubbles, commonly introduced during surgery or precipitated from dissolved gases in decompression sickness, can occlude vessels, diminish perfusion, and initiate thrombotic and inflammatory pathways. The specific interactions between bubbles and the endothelium that lead to cell injury and death are poorly understood. We report live single cell time-lapse imaging of human umbilical vein endothelial cells (HUVECs) perturbed by microbubbles produced by injecting air through a pulled glass capillary ground so that its bevel opposes the buoyant force of the bubble. The bubble is moved with a three-stage micromanipulator into contact with a cell containing Fluo-4, a calcium sensitive dye, and imaged using phase-contrast and fluorescence microscopy. A significant transient elevation in intracellular calcium is observed in response to bubble impact. Membrane depolarization (external K+, 130 mM) does not block the calcium response, indicating that voltage-dependent calcium channels are not involved. Extracellular calcium and an intact actin cytoskeleton are required for the elevation of intracellular calcium upon bubble impact, which is reduced with the addition of gadolinium, a general inhibitor of mechanosensitive channels. These results suggest that a combination of the mechanical deformation and the air-water interface of the bubble activate a plasmalemma ion channel triggering further release of internal calcium stores. The calcium response is ameliorated with addition of surfactant, likely competing for occupancy of the air-liquid bubble interface. Supported by ONR N00014-08-1-0436 and NIH R01 HL67986.