Nfinder: automatic inference of cell neighborhood in 2D and 3D using nuclear markers.

BACKGROUND: In tissues and organisms, the coordination of neighboring cells is essential to maintain their properties and functions. Therefore, knowing which cells are adjacent is crucial to understand biological processes that involve physical interactions among them, e.g. cell migration and proliferation. In addition, some signaling pathways, such as Notch or extrinsic apoptosis, are highly dependent on cell-cell communication. While this is straightforward to obtain from membrane images, nuclei labelling is much more ubiquitous for technical reasons. However, there are no automatic and robust methods to find neighboring cells based only on nuclear markers. RESULTS: In this work, we describe Nfinder, a method to assess the cell's local neighborhood from images with nuclei labeling. To achieve this goal, we approximate the cell-cell interaction graph by the Delaunay triangulation of nuclei centroids. Then, links are filtered by automatic thresholding in cell-cell distance (pairwise interaction) and the maximum angle that a pair of cells subtends with shared neighbors (non-pairwise interaction). We systematically characterized the detection performance by applying Nfinder to publicly available datasets from Drosophila melanogaster, Tribolium castaneum, Arabidopsis thaliana and C. elegans. In each case, the result of the algorithm was compared to a cell neighbor graph generated by manually annotating the original dataset. On average, our method detected 95% of true neighbors, with only 6% of false discoveries. Remarkably, our findings indicate that taking into account non-pairwise interactions might increase the Positive Predictive Value up to + 11.5%. CONCLUSION: Nfinder is the first robust and automatic method for estimating neighboring cells in 2D and 3D based only on nuclear markers and without any free parameters. Using this tool, we found that taking non-pairwise interactions into account improves the detection performance significantly. We believe that using our method might improve the effectiveness of other workflows to study cell-cell interactions from microscopy images. Finally, we also provide a reference implementation in Python and an easy-to-use napari plugin.

Robust and unbiased estimation of the background distribution for automated quantitative imaging.

Background estimation is the first step in quantitative analysis of images. It has an impact on all subsequent analyses, in particular for segmentation and calculation of ratiometric quantities. Most methods recover only a single value such as the median or yield a biased estimation in non-trivial cases. We introduce, to our knowledge, the first method to recover an unbiased estimation of background distribution. It leverages the lack of local spatial correlation in background pixels to robustly select a subset that accurately represents the background. The resulting background distribution can be used to test for foreground membership of individual pixels or estimate confidence intervals in derived quantities.

Smaug1 membrane-less organelles respond to AMPK and mTOR and affect mitochondrial function.

Smaug is a conserved translational regulator that binds numerous mRNAs, including nuclear transcripts that encode mitochondrial enzymes. Smaug orthologs form cytosolic membrane-less organelles (MLOs) in several organisms and cell types. We have performed single-molecule fluorescence in situ hybridization (FISH) assays that revealed that SDHB and UQCRC1 mRNAs associate with Smaug1 bodies in U2OS cells. Loss of function of Smaug1 and Smaug2 (also known as SAMD4A and SAMD4B, respectively) affected both mitochondrial respiration and morphology of the mitochondrial network. Phenotype rescue by Smaug1 transfection depends on the presence of its RNA-binding domain. Moreover, we identified specific Smaug1 domains involved in MLO formation, and found that impaired Smaug1 MLO condensation correlates with mitochondrial defects. Mitochondrial complex I inhibition upon exposure to rotenone, but not strong mitochondrial uncoupling upon exposure to CCCP, rapidly induced the dissolution of Smaug1 MLOs. Metformin and rapamycin elicited similar effects, which were blocked by pharmacological inhibition of AMP-activated protein kinase (AMPK). Finally, we found that Smaug1 MLO dissolution weakens the interaction with target mRNAs, thus enabling their release. We propose that mitochondrial respiration and the AMPK-mTOR balance controls the condensation and dissolution of Smaug1 MLOs, thus regulating nuclear mRNAs that encode key mitochondrial proteins. This article has an associated First Person interview with the first authors of the paper.

A first approach to the use of upconversion nanoparticles to measure fluorescent tracers in water: a proof of concept.

In this work we use lanthanide based NaYF4:Er3+, Yb3+upconversion nanoparticles (UCNP) to detect ppb-level sensitibity of a xanthene dye, Rhodamine B (RB) dye, under NIR excitation. A static energy transfer was observed between the luminescent UCNP energy donors and RB acceptor in aqueous solution for three different sizes of UCNP. No specific covalent functionalization of the UCNPs was performed providing a direct method of detection, particularly promising in natural systems where the interfering fluorescence background is a detrimental limitation to the performance of the detection method. This procedure is a first approach to be applied in estuarine and coastal zone where the high content of suspended particulate matter prevents the detection of tracers.

CASPAM: A Triple-Modality Biosensor for Multiplexed Imaging of Caspase Network Activity.

Understanding signal propagation across biological networks requires to simultaneously monitor the dynamics of several nodes to uncover correlations masked by inherent intercellular variability. To monitor the enzymatic activity of more than two components over short time scales has proven challenging. Exploiting the narrow spectral width of homo-FRET-based biosensors, up to three activities can be imaged through fluorescence polarization anisotropy microscopy. We introduce Caspase Activity Sensor by Polarization Anisotropy Multiplexing (CASPAM) a single-plasmid triple-modality reporter of key nodes of the apoptotic network. Apoptosis provides an ideal molecular framework to study interactions between its three composing pathways (intrinsic, extrinsic, and effector). We characterized the biosensor performance and demonstrated the advantages that equimolar expression has in both simplifying experimental procedure and reducing observable variation, thus enabling robust data-driven modeling. Tools like CASPAM become essential to analyze molecular pathways where multiple nodes need to be simultaneously monitored.

Furry is required for cell movements during gastrulation and functionally interacts with NDR1.

Gastrulation is a key event in animal embryogenesis during which germ layer precursors are rearranged and the embryonic axes are established. Cell polarization is essential during gastrulation, driving asymmetric cell division, cell movements, and cell shape changes. The furry (fry) gene encodes an evolutionarily conserved protein with a wide variety of cellular functions, including cell polarization and morphogenesis in invertebrates. However, little is known about its function in vertebrate development. Here, we show that in Xenopus, Fry plays a role in morphogenetic processes during gastrulation, in addition to its previously described function in the regulation of dorsal mesoderm gene expression. Using morpholino knock-down, we demonstrate a distinct role for Fry in blastopore closure and dorsal axis elongation. Loss of Fry function drastically affects the movement and morphological polarization of cells during gastrulation and disrupts dorsal mesoderm convergent extension, responsible for head-to-tail elongation. Finally, we evaluate a functional interaction between Fry and NDR1 kinase, providing evidence of an evolutionarily conserved complex required for morphogenesis.

Interaction Between the Angiotensin-(1-7) Mas Receptor and the Dopamine D2 Receptor: Implications in Inflammation.

[Figure: see text].

Effective description of bistability and irreversibility in apoptosis.

Apoptosis is a mechanism of programmed cell death in which cells engage in a controlled demolition and prepare to be digested without damaging their environment. In normal conditions, apoptosis is repressed until it is irreversibly induced by an appropriate signal. In adult organisms, apoptosis is a natural way to dispose of damaged cells and its disruption or excess is associated with cancer and autoimmune diseases. Apoptosis is regulated by a complex signaling network controlled by caspases, specialized enzymes that digest essential cellular components and promote the degradation of genomic DNA. In this work, we propose an effective description of the signaling network focused on caspase-3 as a readout of cell fate. We integrate intermediate network interactions into a nonlinear feedback function acting on caspase-3 and introduce the effect of pro-apoptotic stimuli and regulatory elements as a saturating activation function. We show that activation dynamics in the theory is similar to previously reported experimental results. We compute bifurcation diagrams and obtain cell fate maps describing how stimulus intensity and feedback strength affect cell survival and death fates. These fates overlap within a bistable region that depends on total caspase concentration, regulatory elements, and feedback nonlinearity. We study a strongly nonlinear regime to obtain analytical expressions for bifurcation curves and fate map boundaries. For a broad range of parameters, strong stimuli can induce an irreversible switch to the death fate. We use the theory to explore dynamical stimulation conditions and determine how cell fate depends on stimulation temporal patterns. This analysis predicts a critical relation between transient stimuli intensity and duration to trigger irreversible apoptosis. We derive an analytical expression for this critical relation, valid for short stimuli. Our description provides distinct predictions and offers a framework to study how this signaling network processes different stimuli to make a cell fate decision.

Compact and reflective light-sheet microscopy for long-term imaging of living embryos.

The development of light-sheet fluorescence microscopy has been a revolution for developmental biology as it allows long-term imaging during embryonic development. An important reason behind the quick adoption has been the availability of open hardware alternatives. In this work, we present a robust and compact version of a light-sheet fluorescence microscope that is easy to assemble and requires little to no maintenance. An important aspect of the design is that the illumination unit consists of reflective elements, thereby reducing chromatic aberrations an order of magnitude as compared to refractive counterparts.

Robustness in spatially driven bistability in signaling systems.

Biological systems are spatially organized. This microscopic heterogeneity has been shown to produce emergent complex behaviors such as bistability. Even though the connection between spatiality and dynamic response is essential to understand biological output, its robustness and extent has not been sufficiently explored. This work focuses on a previously described system which is composed of two monostable modules acting on different cellular compartments and sharing species through linear shuttling reactions. One of the two main purposes of this paper is to quantify the frequency of occurrence of bistability throughout the parameter space and to identify which parameters and in which value ranges control the emergence and the properties of bistability. We found that a very small fraction of the sampled parameter space produced a bistable response. Most importantly, shuttling parameters were among the most influential ones to control this property. The other goal of this paper is to simplify the same system as much as possible without losing compartment-induced bistability. This procedure provided a simplified model that still connects two monostable systems by a reduced set of linear shuttling reactions that circulates all the species around the two compartments. Bistable systems are one of the main building blocks of more complex behaviors such as oscillations, memory, and digitalization. Therefore, we expect that the proposed minimal system provides insight into how these behaviors can arise from compartmentalization.

An Integrative and Modular Framework to Recapitulate Emergent Behavior in Cell Migration.

Cell migration has been a subject of study in a broad variety of biological systems, from morphogenetic events during development to cancer progression. In this work, we describe single-cell movement in a modular framework from which we simulate the collective behavior of glioblastoma cells, the most prevalent and malignant primary brain tumor. We used the U87 cell line, which can be grown as a monolayer or spatially closely packed and organized in 3D structures called spheroids. Our integrative model considers the most relevant mechanisms involved in cell migration: chemotaxis of attractant factor, mechanical interactions and random movement. The effect of each mechanism is integrated into the overall probability of the cells to move in a particular direction, in an automaton-like approach. Our simulations fit and reproduced the emergent behavior of the spheroids in a set of migration assays where single-cell trajectories were tracked. We also predicted the effect of migration inhibition on the colonies from simple experimental characterization of single treated cell tracks. The development of tools that allow complementing molecular knowledge in migratory cell behavior is relevant for understanding essential cellular processes, both physiological (such as organ formation, tissue regeneration among others) and pathological perspectives. Overall, this is a versatile tool that has been proven to predict individual and collective behavior in U87 cells, but that can be applied to a broad variety of scenarios.

Co-imaging extrinsic, intrinsic and effector caspase activity by fluorescence anisotropy microscopy.

In order to overcome intercellular variability and thereby effectively assess signal propagation in biological networks it is imperative to simultaneously quantify multiple biological observables in single living cells. While fluorescent biosensors have been the tool of choice to monitor the dynamics of protein interaction and enzymatic activity, co-measuring more than two of them has proven challenging. In this work, we designed three spectrally separated anisotropy-based Förster Resonant Energy Transfer (FRET) biosensors to overcome this difficulty. We demonstrate this principle by monitoring the activation of extrinsic, intrinsic and effector caspases upon apoptotic stimulus. Together with modelling and simulations we show that time of maximum activity for each caspase can be derived from the anisotropy of the corresponding biosensor. Such measurements correlate relative activation times and refine existing models of biological signalling networks, providing valuable insight into signal propagation.

Interdependence between EGFR and Phosphatases Spatially Established by Vesicular Dynamics Generates a Growth Factor Sensing and Responding Network.

The proto-oncogenic epidermal growth factor receptor (EGFR) is a tyrosine kinase whose sensitivity to growth factors and signal duration determines cellular behavior. We resolve how EGFR's response to epidermal growth factor (EGF) originates from dynamically established recursive interactions with spatially organized protein tyrosine phosphatases (PTPs). Reciprocal genetic PTP perturbations enabled identification of receptor-like PTPRG/J at the plasma membrane and ER-associated PTPN2 as the major EGFR dephosphorylating activities. Imaging spatial-temporal PTP reactivity revealed that vesicular trafficking establishes a spatially distributed negative feedback with PTPN2 that determines signal duration. On the other hand, single-cell dose-response analysis uncovered a reactive oxygen species-mediated toggle switch between autocatalytically activated monomeric EGFR and the tumor suppressor PTPRG that governs EGFR's sensitivity to EGF. Vesicular recycling of monomeric EGFR unifies the interactions with these PTPs on distinct membrane systems, dynamically generating a network architecture that can sense and respond to time-varying growth factor signals.

Oxygen-Insensitive Aggregates of Pt(II) Complexes as Phosphorescent Labels of Proteins with Luminescence Lifetime-Based Readouts.

The synthesis and photophysical properties of a tailored Pt(II) complex are presented. The phosphorescence of its monomeric species in homogeneous solutions is quenched by interaction with the solvent and therefore absent even upon deoxygenation. However, aggregation-induced shielding from the environment and suppression of rotovibrational degrees of freedom trigger a phosphorescence turn-on that is not suppressed by molecular oxygen, despite possessing an excited-state lifetime ranging in the microsecond scale. Thus, the photoinduced production of reactive oxygen species is avoided by the suppression of diffusion-controlled Dexter-type energy transfer to triplet molecular oxygen. These aggregates emit with the characteristic green luminescence profile of monomeric complexes, indicating that Pt-Pt or excimeric interactions are negligible. Herein, we show that these aggregates can be used to label a model biomolecule (bovine serum albumin) with a microsecond-range luminescence. The protein stabilizes the aggregates, acting as a carrier in aqueous environments. Despite spectral overlaps, the green phosphorescence can be separated by time-gated detection from the dominant autofluorescence of the protein arising from a covalently bound green fluorophore that emits in the nanosecond range. Interestingly, the aggregates also acted as energy donors able to sensitize the emission of a fraction of the fluorophores bound to the protein. This resulted in a microsecond-range luminescence of the fluorescent acceptors and a shortening of the excited-state lifetime of the phosphorescent aggregates. The process that can be traced by a 1000-fold increase in the acceptor's lifetime mirrors the donor's triplet character. The implications for phosphorescence lifetime imaging are discussed.

Diverse patterns of molecular changes in the mechano-responsiveness of focal adhesions.

Focal adhesions anchor contractile actin fibers with the extracellular matrix, sense the generated tension and respond to it by changing their morphology and composition. Here we ask how this mechanosensing is enabled at the protein-network level, given the modular assembly and multitasking of focal adhesions. To address this, we applied a sensitive 4-color live cell imaging approach, enabling monitoring patterns of molecular changes in single focal adhesions. Co-imaging zyxin, FAK, vinculin and paxillin revealed heterogeneities in their responses to Rho-associated kinase (ROCK)-mediated perturbations of actomyosin contractility. These responses were rather weakly correlated between the proteins, reflecting diverse compositional changes in different focal adhesions. This diversity is partially attributable to the location of focal adhesions, their area, molecular content and previous contractility perturbations, suggesting that integration of multiple local cues shapes differentially focal adhesion mechano-responsiveness. Importantly, the compositional changes upon ROCK perturbations exhibited distinct paths in different focal adhesions. Moreover, the protein exhibiting the strongest response to ROCK perturbations varied among different focal adhesions. The diversity in response patterns is plausibly enabled by the modular mode of focal adhesions assembly and can provide them the needed flexibility to perform multiple tasks by combining optimally a common set of multifunctional components.

MAS1 Receptor Trafficking Involves ERK1/2 Activation Through a ß-Arrestin2-Dependent Pathway.

The MAS1 receptor (R) exerts protective effects in the brain, heart, vessels, and kidney. R trafficking plays a critical function in signal termination and propagation and in R resensitization. We examined MAS1R internalization and trafficking on agonist stimulation and the role of ß-arrestin2 in the activation of ERK1/2 (extracellular signal-regulated kinase 1/2) and Akt after MAS1R stimulation. Human embryonic kidney 293T cells were transfected with the coding sequence for MAS1R-YFP (MAS1R fused to yellow fluorescent protein). MAS1R internalization was evaluated by measuring the MAS1R present in the plasma membrane after agonist stimulation using a ligand-binding assay. MAS1R trafficking was evaluated by its colocalization with trafficking markers. MAS1R internalization was blocked in the presence of shRNAcaveolin-1 and with dominant negatives for Eps15 (a protein involved in endocytosed Rs by clathrin-coated pits) and for dynamin. After stimulation, MAS1R colocalized with Rab11-a slow recycling vesicle marker-and not with Rab4-a fast recycling vesicle marker-or LysoTracker-a lysosome marker. Cells transfected with MAS1R showed an increase in Akt and ERK1/2 activation on angiotensin-(1-7) stimulation, which was blocked when the clathrin-coated pits pathway was blocked. Suppression of ß-arrestin2 by shRNA reduced the angiotensin-(1-7)-induced ERK1/2 activation, whereas Akt activation was not modified. We conclude that on agonist stimulation, MAS1R is internalized through clathrin-coated pits and caveolae in a dynamin-dependent manner and is then slowly recycled back to the plasma membrane. MAS1R induced Akt and ERK1/2 activation from early endosomes, and the activation of ERK1/2 was mediated by ß-arrestin2. Thus, MAS1R activity and density may be tightly controlled by the cell.

pawFLIM: reducing bias and uncertainty to enable lower photon count in FLIM experiments.

Förster resonant energy transfer measured by fluorescence lifetime imaging microscopy (FRET-FLIM) is the method of choice for monitoring the spatio-temporal dynamics of protein interactions in living cells. To obtain an accurate estimate of the molecular fraction of interacting proteins requires a large number of photons, which usually precludes the observation of a fast process, particularly with time correlated single photon counting (TCSPC) based FLIM. In this work, we propose a novel method named pawFLIM (phasor analysis via wavelets) that allows the denoising of FLIM datasets by adaptively and selectively adjusting the desired compromise between spatial and molecular resolution. The method operates by applying a weighted translational-invariant Haar-wavelet transform denoising algorithm to phasor images. This results in significantly less bias and mean square error than other existing methods. We also present a new lifetime estimator (named normal lifetime) with a smaller mean squared error and overall bias as compared to frequency domain phase and modulation lifetimes. Overall, we present an approach that will enable the observation of the dynamics of biological processes at the molecular level with better temporal and spatial resolution.

Heteromerization Between the Bradykinin B2 Receptor and the Angiotensin-(1-7) Mas Receptor: Functional Consequences.

Bradykinin B2 receptor (B2R) and angiotensin-(1-7) Mas receptor (MasR)-mediated effects are physiologically interconnected. The molecular basis for such cross talk is unknown. It is hypothesized that the cross talk occurs at the receptor level. We investigated B2R-MasR heteromerization and the functional consequences of such interaction. B2R fused to the cyan fluorescent protein and MasR fused to the yellow fluorescent protein were transiently coexpressed in human embryonic kidney293T cells. Fluorescence resonance energy transfer analysis showed that B2R and MasR formed a constitutive heteromer, which was not modified by their agonists. B2R or MasR antagonists decreased fluorescence resonance energy transfer efficiency, suggesting that the antagonist promoted heteromer dissociation. B2R-MasR heteromerization induced an 8-fold increase in the MasR ligand-binding affinity. On agonist stimulation, the heteromer was internalized into early endosomes with a slower sequestration rate from the plasma membrane, compared with single receptors. B2R-MasR heteromerization induced a greater increase in arachidonic acid release and extracellular signal-regulated kinase phosphorylation after angiotensin-(1-7) stimulation, and this effect was blocked by the B2R antagonist. Concerning serine/threonine kinase Akt activity, a significant bradykinin-promoted activation was detected in B2R-MasR but not in B2R-expressing cells. Angiotensin-(1-7) and bradykinin elicited antiproliferative effects only in cells expressing B2R-MasR heteromers, but not in cells expressing each receptor alone. Proximity ligation assay confirmed B2R-MasR interaction in human glomerular endothelial cells supporting the interaction between both receptors in vivo. Our findings provide an explanation for the cross talk between bradykinin B2R and angiotensin-(1-7) MasR-mediated effects. B2R-MasR heteromerization induces functional changes in the receptor that may lead to long-lasting protective properties.

Highly Multiplexed Imaging Uncovers Changes in Compositional Noise within Assembling Focal Adhesions.

Integrin adhesome proteins bind each other in alternative manners, forming within the cell diverse cell-matrix adhesion sites with distinct properties. An intriguing question is how such modular assembly of adhesion sites is achieved correctly solely by self-organization of their components. Here we address this question using high-throughput multiplexed imaging of eight proteins and two phosphorylation sites in a large number of single focal adhesions. We found that during the assembly of focal adhesions the variances of protein densities decrease while the correlations between them increase, suggesting reduction in the noise levels within these structures. These changes correlate independently with the area and internal density of focal adhesions, but not with their age or shape. Artificial neural network analysis indicates that a joint consideration of multiple components improves the predictability of paxillin and zyxin levels in internally dense focal adhesions. This suggests that paxillin and zyxin densities in focal adhesions are fine-tuned by integrating the levels of multiple other components, thus averaging-out stochastic fluctuations. Based on these results we propose that increase in internal protein densities facilitates noise suppression in focal adhesions, while noise suppression enables their stable growth and further density increase-hence forming a feedback loop giving rise to a quality-controlled assembly.

Multiplexed imaging of intracellular protein networks.

Cellular functions emerge from the collective action of a large number of different proteins. Understanding how these protein networks operate requires monitoring their components in intact cells. Due to intercellular and intracellular molecular variability, it is important to monitor simultaneously multiple components at high spatiotemporal resolution. However, inherent trade-offs narrow the boundaries of achievable multiplexed imaging. Pushing these boundaries is essential for a better understanding of cellular processes. Here the motivations, challenges and approaches for multiplexed imaging of intracellular protein networks are discussed. © 2016 International Society for Advancement of Cytometry.