Automatic monitoring of neural activity with single-cell resolution in behaving Hydra.

The ability to record every spike from every neuron in a behaving animal is one of the holy grails of neuroscience. Here, we report coming one step closer towards this goal with the development of an end-to-end pipeline that automatically tracks and extracts calcium signals from individual neurons in the cnidarian Hydra vulgaris. We imaged dually labeled (nuclear tdTomato and cytoplasmic GCaMP7s) transgenic Hydra and developed an open-source Python platform (TraSE-IN) for the Tracking and Spike Estimation of Individual Neurons in the animal during behavior. The TraSE-IN platform comprises a series of modules that segments and tracks each nucleus over time and extracts the corresponding calcium activity in the GCaMP channel. Another series of signal processing modules allows robust prediction of individual spikes from each neuron's calcium signal. This complete pipeline will facilitate the automatic generation and analysis of large-scale datasets of single-cell resolution neural activity in Hydra, and potentially other model organisms, paving the way towards deciphering the neural code of an entire animal.

Automatic monitoring of whole-body neural activity in behaving Hydra.

The ability to record every spike from every neuron in a behaving animal is one of the holy grails of neuroscience. Here, we report coming one step closer towards this goal with the development of an end-to-end pipeline that automatically tracks and extracts calcium signals from individual neurons in the cnidarian Hydra vulgaris. We imaged dually labeled (nuclear tdTomato and cytoplasmic GCaMP7s) transgenic Hydra and developed an open-source Python platform (TraSE-IN) for the Tracking and Spike Estimation of Individual Neurons in the animal during behavior. The TraSE-IN platform comprises a series of modules that segments and tracks each nucleus over time and extracts the corresponding calcium activity in the GCaMP channel. Another series of signal processing modules allows robust prediction of individual spikes from each neuron's calcium signal. This complete pipeline will facilitate the automatic generation and analysis of large-scale datasets of single-cell resolution neural activity in Hydra, and potentially other model organisms, paving the way towards deciphering the neural code of an entire animal.

Tracking receptor motions at the plasma membrane reveals distinct effects of ligands on CCR5 dynamics depending on its dimerization status.

G-protein-coupled receptors (GPCR) are present at the cell surface in different conformational and oligomeric states. However, how these states impact GPCRs biological function and therapeutic targeting remains incompletely known. Here, we investigated this issue in living cells for the CC chemokine receptor 5 (CCR5), a major receptor in inflammation and the principal entry co-receptor for Human Immunodeficiency Viruses type 1 (HIV-1). We used TIRF microscopy and a statistical method to track and classify the motion of different receptor subpopulations. We showed a diversity of ligand-free forms of CCR5 at the cell surface constituted of various oligomeric states and exhibiting transient Brownian and restricted motions. These forms were stabilized differently by distinct ligands. In particular, agonist stimulation restricted the mobility of CCR5 and led to its clustering, a feature depending on ß-arrestin, while inverse agonist stimulation exhibited the opposite effect. These results suggest a link between receptor activation and immobilization. Applied to HIV-1 envelope glycoproteins gp120, our quantitative analysis revealed agonist-like properties of gp120s. Distinct gp120s influenced CCR5 dynamics differently, suggesting that they stabilize different CCR5 conformations. Then, using a dimerization-compromized mutant, we showed that dimerization (i) impacts CCR5 precoupling to G proteins, (ii) is a pre-requisite for the immobilization and clustering of receptors upon activation, and (iii) regulates receptor endocytosis, thereby impacting the fate of activated receptors. This study demonstrates that tracking the dynamic behavior of a GPCR is an efficient way to link GPCR conformations to their functions, therefore improving the development of drugs targeting specific receptor conformations.

Tracking calcium dynamics from individual neurons in behaving animals.

Measuring the activity of neuronal populations with calcium imaging can capture emergent functional properties of neuronal circuits with single cell resolution. However, the motion of freely behaving animals, together with the intermittent detectability of calcium sensors, can hinder automatic monitoring of neuronal activity and their subsequent functional characterization. We report the development and open-source implementation of a multi-step cellular tracking algorithm (Elastic Motion Correction and Concatenation or EMC2) that compensates for the intermittent disappearance of moving neurons by integrating local deformation information from detectable neurons. We demonstrate the accuracy and versatility of our algorithm using calcium imaging data from two-photon volumetric microscopy in visual cortex of awake mice, and from confocal microscopy in behaving Hydra, which experiences major body deformation during its contractions. We quantify the performance of our algorithm using ground truth manual tracking of neurons, along with synthetic time-lapse sequences, covering a wide range of particle motions and detectability parameters. As a demonstration of the utility of the algorithm, we monitor for several days calcium activity of the same neurons in layer 2/3 of mouse visual cortex in vivo, finding significant turnover within the active neurons across days, with only few neurons that remained active across days. Also, combining automatic tracking of single neuron activity with statistical clustering, we characterize and map neuronal ensembles in behaving Hydra, finding three major non-overlapping ensembles of neurons (CB, RP1 and RP2) whose activity correlates with contractions and elongations. Our results show that the EMC2 algorithm can be used as a robust and versatile platform for neuronal tracking in behaving animals.

Cytokine receptor cluster size impacts its endocytosis and signaling.

The interleukin-2 receptor (IL-2R) is a cytokine receptor essential for immunity that transduces proliferative signals regulated by its uptake and degradation. IL-2R is a well-known marker of clathrin-independent endocytosis (CIE), a process devoid of any coat protein, raising the question of how the CIE vesicle is generated. Here, we investigated the impact of IL-2Rγ clustering in its endocytosis. Combining total internal reflection fluorescence (TIRF) live imaging of a CRISPR-edited T cell line endogenously expressing IL-2Rγ tagged with green fluorescent protein (GFP), with multichannel imaging, single-molecule tracking, and quantitative analysis, we were able to decipher IL-2Rγ stoichiometry at the plasma membrane in real time. We identified three distinct IL-2Rγ cluster populations. IL-2Rγ is secreted to the cell surface as a preassembled small cluster of three molecules maximum, rapidly diffusing at the plasma membrane. A medium-sized cluster composed of four to six molecules is key for IL-2R internalization and is promoted by interleukin 2 (IL-2) binding, while larger clusters (more than six molecules) are static and inefficiently internalized. Moreover, we identified membrane cholesterol and the branched actin cytoskeleton as key regulators of IL-2Rγ clustering and IL-2-induced signaling. Both cholesterol depletion and Arp2/3 inhibition lead to the assembly of large IL-2Rγ clusters, arising from the stochastic interaction of receptor molecules in close correlation with their enhanced lateral diffusion at the membrane, thus resulting in a default in IL-2R endocytosis. Despite similar clustering outcomes, while cholesterol depletion leads to a sustained IL-2-dependent signaling, Arp2/3 inhibition prevents signal initiation. Taken together, our results reveal the importance of cytokine receptor clustering for CIE initiation and signal transduction.

Evaluating the Stability of Spatial Keypoints via Cluster Core Correspondence Index.

Detection and analysis of informative keypoints is a fundamental problem in image analysis and computer vision. Keypoint detectors are omnipresent in visual automation tasks, and recent years have witnessed a significant surge in the number of such techniques. Evaluating the quality of keypoint detectors remains a challenging task owing to the inherent ambiguity over what constitutes a good keypoint. In this context, we introduce a reference based keypoint quality index which is based on the theory of spatial pattern analysis. Unlike traditional correspondence-based quality evaluation which counts the number of feature matches within a specified neighborhood, we present a rigorous mathematical framework to compute the statistical correspondence of the detections inside a set of salient zones (cluster cores) defined by the spatial distribution of a reference set of keypoints. We leverage the versatility of the level sets to handle hypersurfaces of arbitrary geometry, and develop a mathematical framework to estimate the model parameters analytically to reflect the robustness of a feature detection algorithm. Extensive experimental studies involving several keypoint detectors tested under different imaging scenarios demonstrate efficacy of our method to evaluate keypoint quality for generic applications in computer vision and image analysis.

Electrodiffusion models of synaptic potentials in dendritic spines.

The biophysical properties of dendritic spines play a critical role in neuronal integration but are still poorly understood, due to experimental difficulties in accessing them. Spine biophysics has been traditionally explored using theoretical models based on cable theory. However, cable theory generally assumes that concentration changes associated with ionic currents are negligible and, therefore, ignores electrodiffusion, i.e. the interaction between electric fields and ionic diffusion. This assumption, while true for large neuronal compartments, could be incorrect when applied to femto-liter size structures such as dendritic spines. To extend cable theory and explore electrodiffusion effects, we use here the Poisson (P) and Nernst-Planck (NP) equations, which relate electric field to charge and Fick's law of diffusion, to model ion concentration dynamics in spines receiving excitatory synaptic potentials (EPSPs). We use experimentally measured voltage transients from spines with nanoelectrodes to explore these dynamics with realistic parameters. We find that (i) passive diffusion and electrodiffusion jointly affect the dynamics of spine EPSPs; (ii) spine geometry plays a key role in shaping EPSPs; and, (iii) the spine-neck resistance dynamically decreases during EPSPs, leading to short-term synaptic facilitation. Our formulation, which complements and extends cable theory, can be easily adapted to model ionic biophysics in other nanoscale bio-compartments.

Real-time analysis of the behaviour of groups of mice via a depth-sensing camera and machine learning.

Preclinical studies of psychiatric disorders use animal models to investigate the impact of environmental factors or genetic mutations on complex traits such as decision-making and social interactions. Here, we introduce a method for the real-time analysis of the behaviour of mice housed in groups of up to four over several days and in enriched environments. The method combines computer vision through a depth-sensing infrared camera, machine learning for animal and posture identification, and radio-frequency identification to monitor the quality of mouse tracking. It tracks multiple mice accurately, extracts a list of behavioural traits of both individuals and the groups of mice, and provides a phenotypic profile for each animal. We used the method to study the impact of Shank2 and Shank3 gene mutations-mutations that are associated with autism-on mouse behaviour. Characterization and integration of data from the behavioural profiles of Shank2 and Shank3 mutant female mice revealed their distinctive activity levels and involvement in complex social interactions.

Histamine releasing factor and elongation factor 1 alpha secreted via malaria parasites extracellular vesicles promote immune evasion by inhibiting specific T cell responses.

Protozoan pathogens secrete nanosized particles called extracellular vesicles (EVs) to facilitate their survival and chronic infection. Here, we show the inhibition by Plasmodium berghei NK65 blood stage-derived EVs of the proliferative response of CD4+ T cells in response to antigen presentation. Importantly, these results were confirmed in vivo by the capacity of EVs to diminish the ovalbumin-specific delayed type hypersensitivity response. We identified two proteins associated with EVs, the histamine releasing factor (HRF) and the elongation factor 1α (EF-1α) that were found to have immunosuppressive activities. Interestingly, in contrast to WT parasites, EVs from genetically HRF- and EF-1α-deficient parasites failed to inhibit T cell responses in vitro and in vivo. At the level of T cells, we demonstrated that EVs from WT parasites dephosphorylate key molecules (PLCγ1, Akt, and ERK) of the T cell receptor signalling cascade. Remarkably, immunisation with EF-1α alone or in combination with HRF conferred a long-lasting antiparasite protection and immune memory. In conclusion, we identified a new mechanism by which P. berghei-derived EVs exert their immunosuppressive functions by altering T cell responses. The identification of two highly conserved immune suppressive factors offers new conceptual strategies to overcome EV-mediated immune suppression in malaria-infected individuals.

Mapping molecular assemblies with fluorescence microscopy and object-based spatial statistics.

Elucidating protein functions and molecular organisation requires to localise precisely single or aggregated molecules and analyse their spatial distributions. We develop a statistical method SODA (Statistical Object Distance Analysis) that uses either micro- or nanoscopy to significantly improve on standard co-localisation techniques. Our method considers cellular geometry and densities of molecules to provide statistical maps of isolated and associated (coupled) molecules. We use SODA with three-colour structured-illumination microscopy (SIM) images of hippocampal neurons, and statistically characterise spatial organisation of thousands of synapses. We show that presynaptic synapsin is arranged in asymmetric triangle with the 2 postsynaptic markers homer and PSD95, indicating a deeper localisation of homer. We then determine stoichiometry and distance between localisations of two synaptic vesicle proteins with 3D-STORM. These findings give insights into the protein organisation at the synapse, and prove the efficiency of SODA to quantitatively assess the geometry of molecular assemblies.

Quantitative and Statistical Study of the Dynamics of Clathrin-Dependent and -Independent Endocytosis Reveal a Differential Role of EndophilinA2.

Eukaryotic cells internalize cargos specifically through clathrin-mediated endocytosis (CME) or clathrin-independent endocytosis (CIE). EndophilinA2 was shown as preferentially implicated in CIE, although initially involved in CME. Here, we investigated the native interplay of endophilinA2 and dynamin2 during CME as compared to CIE. We developed an unbiased integrative approach based on genome engineering, robust tracking methodology, and advanced analytics. We statistically identified CME and CIE subpopulations corresponding to abortive, active, and static endocytic events. Depletion of dynamin2 strongly affected active CME and CIE events, whereas the absence of endophilinA2 impacted only CIE. Accordingly, we demonstrated that endophilinA2 is needed for dynamin2 recruitment during CIE, but not in CME. Despite these differences, endophilinA2 and dynamin2 acted at the latest stage of endocytosis within a similar stoichiometry in both mechanisms. Thus, we propose a conserved function of dynamin2 and endophilinA2 in vesicle scission, but a differential regulation of their recruitment during CME and CIE.

Legionella pneumophila Modulates Mitochondrial Dynamics to Trigger Metabolic Repurposing of Infected Macrophages.

The intracellular bacteria Legionella pneumophila encodes a type IV secretion system (T4SS) that injects effector proteins into macrophages in order to establish and replicate within the Legionella-containing vacuole (LCV). Once generated, the LCV interacts with mitochondria through unclear mechanisms. We show that Legionella uses both T4SS-independent and T4SS-dependent mechanisms to respectively interact with mitochondria and induce mitochondrial fragmentation that ultimately alters mitochondrial metabolism. The T4SS effector MitF, a Ran GTPase activator, is required for fission of the mitochondrial network. These effects of MitF occur through accumulation of mitochondrial DNM1L, a GTPase critical for fission. Furthermore mitochondrial respiration is abruptly halted in a T4SS-dependent manner, while T4SS-independent upregulation of cellular glycolysis remains elevated. Collectively, these alterations in mitochondrial dynamics promote a Warburg-like phenotype in macrophages that favors bacterial replication. Hence the rewiring of cellular bioenergetics to create a replication permissive niche in host cells is a virulence strategy of L. pneumophila.

A role for septin 2 in Drp1-mediated mitochondrial fission.

Mitochondria are essential eukaryotic organelles often forming intricate networks. The overall network morphology is determined by mitochondrial fusion and fission. Among the multiple mechanisms that appear to regulate mitochondrial fission, the ER and actin have recently been shown to play an important role by mediating mitochondrial constriction and promoting the action of a key fission factor, the dynamin-like protein Drp1. Here, we report that the cytoskeletal component septin 2 is involved in Drp1-dependent mitochondrial fission in mammalian cells. Septin 2 localizes to a subset of mitochondrial constrictions and directly binds Drp1, as shown by immunoprecipitation of the endogenous proteins and by pulldown assays with recombinant proteins. Depletion of septin 2 reduces Drp1 recruitment to mitochondria and results in hyperfused mitochondria and delayed FCCP-induced fission. Strikingly, septin depletion also affects mitochondrial morphology in Caenorhabditis elegans, strongly suggesting that the role of septins in mitochondrial dynamics is evolutionarily conserved.

Membrane protrusion powers clathrin-independent endocytosis of interleukin-2 receptor.

Endocytosis controls many functions including nutrient uptake, cell division, migration and signal transduction. A clathrin- and caveolin-independent endocytosis pathway is used by important physiological cargos, including interleukin-2 receptors (IL-2R). However, this process lacks morphological and dynamic data. Our electron microscopy (EM) and tomography studies reveal that IL-2R-pits and vesicles are initiated at the base of protrusions. We identify the WAVE complex as a specific endocytic actor. The WAVE complex interacts with IL-2R, via a WAVE-interacting receptor sequence (WIRS) present in the receptor polypeptide, and allows for receptor clustering close to membrane protrusions. In addition, using total internal reflection fluorescent microscopy (TIRF) and automated analysis we demonstrate that two timely distinct bursts of actin polymerization are required during IL-2R uptake, promoted first by the WAVE complex and then by N-WASP. Finally, our data reveal that dynamin acts as a transition controller for the recruitment of Arp2/3 activators required for IL-2R endocytosis. Altogether, our work identifies the spatio-temporal specific role of factors initiating clathrin-independent endocytosis by a unique mechanism that does not depend on the deformation of a flat membrane, but rather on that of membrane protrusions.

Imaging intraflagellar transport in trypanosomes.

Trypanosoma brucei is a flagellated eukaryotic pathogen responsible for sleeping sickness in central Africa. Because of the presence of a long motile flagellum (>20 µm) and its amenity to genetic manipulation, it is becoming an attractive model to study the assembly and the functions of cilia and flagella. In recent years, several aspects have been investigated, especially intraflagellar transport (IFT) that has been exhaustively characterized at the light microscopy level. In this manuscript, we review various methods to express fluorescent fusion proteins and to record IFT in living trypanosomes in normal or mutant contexts. We present an approach for separating anterograde and retrograde IFT, hence facilitating quantification of train speed, frequency, and size. A statistical analysis to discriminate different subpopulations of IFT trains is also summarized. These methods have proven their efficiency for the study of IFT in trypanosomes and could be applied to any other organism.

KDEL Receptors Assist Dengue Virus Exit from the Endoplasmic Reticulum.

Membrane receptors at the surface of target cells are key host factors for virion entry; however, it is unknown whether trafficking and secretion of progeny virus requires host intracellular receptors. In this study, we demonstrate that dengue virus (DENV) interacts with KDEL receptors (KDELR), which cycle between the ER and Golgi apparatus, for vesicular transport from ER to Golgi. Depletion of KDELR by siRNA reduced egress of both DENV progeny and recombinant subviral particles (RSPs). Coimmunoprecipitation of KDELR with dengue structural protein prM required three positively charged residues at the N terminus, whose mutation disrupted protein interaction and inhibited RSP transport from the ER to the Golgi. Finally, siRNA depletion of class II Arfs, which results in KDELR accumulation in the Golgi, phenocopied results obtained with mutagenized prME and KDELR knockdown. Our results have uncovered a function for KDELR as an internal receptor involved in DENV trafficking.

Statistical analysis of molecule colocalization in bioimaging.

The quantitative analysis of molecule interactions in bioimaging is key for understanding the molecular orchestration of cellular processes and is generally achieved through the study of the spatial colocalization between the different populations of molecules. Colocalization methods are traditionally divided into pixel-based methods that measure global correlation coefficients from the overlap between pixel intensities in different color channels, and object-based methods that first segment molecule spots and then analyze their spatial distributions with second-order statistics. Here, we present a review of such colocalization methods and give a quantitative comparison of their relative merits in different types of biological applications and contexts. We show on synthetic and biological images that object-based methods are more robust statistically than pixel-based methods, and allow moreover to quantify accurately the number of colocalized molecules.

Membrane trafficking. Nucleoside diphosphate kinases fuel dynamin superfamily proteins with GTP for membrane remodeling.

Dynamin superfamily molecular motors use guanosine triphosphate (GTP) as a source of energy for membrane-remodeling events. We found that knockdown of nucleoside diphosphate kinases (NDPKs) NM23-H1/H2, which produce GTP through adenosine triphosphate (ATP)-driven conversion of guanosine diphosphate (GDP), inhibited dynamin-mediated endocytosis. NM23-H1/H2 localized at clathrin-coated pits and interacted with the proline-rich domain of dynamin. In vitro, NM23-H1/H2 were recruited to dynamin-induced tubules, stimulated GTP-loading on dynamin, and triggered fission in the presence of ATP and GDP. NM23-H4, a mitochondria-specific NDPK, colocalized with mitochondrial dynamin-like OPA1 involved in mitochondria inner membrane fusion and increased GTP-loading on OPA1. Like OPA1 loss of function, silencing of NM23-H4 but not NM23-H1/H2 resulted in mitochondrial fragmentation, reflecting fusion defects. Thus, NDPKs interact with and provide GTP to dynamins, allowing these motor proteins to work with high thermodynamic efficiency.

Early sequence of events triggered by the interaction of Neisseria meningitidis with endothelial cells.

Neisseria meningitidis is a bacterium responsible for severe sepsis and meningitis. Following type IV pilus-mediated adhesion to endothelial cells, bacteria proliferating on the cellular surface trigger a potent cellular response that enhances the ability of adhering bacteria to resist the mechanical forces generated by the blood flow. This response is characterized by the formation of numerous 100 nm wide membrane protrusions morphologically related to filopodia. Here, a high-resolution quantitative live-cell fluorescence microscopy procedure was designed and used to study this process. A farnesylated plasma membrane marker was first detected only a few seconds after bacterial contact, rapidly followed by actin cytoskeleton reorganization and bulk cytoplasm accumulation. The bacterial type IV pili-associated minor pilin PilV is necessary for the initiation of this cascade. Plasma membrane composition is a key factor as cholesterol depletion with methyl-ß-cyclodextrin completely blocks the initiation of the cellular response. In contrast membrane deformation does not require the actin cytoskeleton. Strikingly, plasma membrane remodelling undermicrocolonies is also independent of common intracellular signalling pathways as cellular ATP depletion is not inhibitory. This study shows that bacteria-induced plasma membrane reorganization is a rapid event driven by a direct cross-talk between type IV pili and the plasma membrane rather than by the activation of an intracellular signalling pathway that would lead to actin remodelling.

Analysis of the spatial organization of molecules with robust statistics.

One major question in molecular biology is whether the spatial distribution of observed molecules is random or organized in clusters. Indeed, this analysis gives information about molecules' interactions and physical interplay with their environment. The standard tool for analyzing molecules' distribution statistically is the Ripley's K function, which tests spatial randomness through the computation of its critical quantiles. However, quantiles' computation is very cumbersome, hindering its use. Here, we present an analytical expression of these quantiles, leading to a fast and robust statistical test, and we derive the characteristic clusters' size from the maxima of the Ripley's K function. Subsequently, we analyze the spatial organization of endocytic spots at the cell membrane and we report that clathrin spots are randomly distributed while clathrin-independent spots are organized in clusters with a radius of 2 µm, which suggests distinct physical mechanisms and cellular functions for each pathway.