Time course decomposition of cell heterogeneity in TFEB signaling states reveals homeostatic mechanisms restricting the magnitude and duration of TFEB responses to mTOR activity modulation.

BACKGROUND: TFEB (transcription factor EB) regulates metabolic homeostasis through its activation of lysosomal biogenesis following its nuclear translocation. TFEB activity is inhibited by mTOR phosphorylation, which signals its cytoplasmic retention. To date, the temporal relationship between alterations to mTOR activity states and changes in TFEB subcellular localization and concentration has not been sufficiently addressed. METHODS: mTOR was activated by renewed addition of fully-supplemented medium, or inhibited by Torin1 or nutrient deprivation. Single-cell TFEB protein levels and subcellular localization in HeLa and MCF7 cells were measured over a time course of 15 hours by multispectral imaging cytometry. To extract single-cell level information on heterogeneous TFEB activity phenotypes, we developed a framework for identification of TFEB activity subpopulations. Through unsupervised clustering, cells were classified according to their TFEB nuclear concentration, which corresponded with downstream lysosomal responses. RESULTS: Bulk population results revealed that mTOR negatively regulates TFEB protein levels, concomitantly to the regulation of TFEB localization. Subpopulation analysis revealed maximal sensitivity of HeLa cells to mTOR activity stimulation, leading to inactivation of 100 % of the cell population within 0.5 hours, which contrasted with a lower sensitivity in MCF7 cells. Conversely, mTOR inhibition increased the fully active subpopulation only fractionally, and full activation of 100 % of the population required co-inhibition of mTOR and the proteasome. Importantly, mTOR inhibition activated TFEB for a limited duration of 1.5 hours, and thereafter the cell population was progressively re-inactivated, with distinct kinetics for Torin1 and nutrient deprivation treatments. CONCLUSION: TFEB protein levels and subcellular localization are under control of a short-term rheostat, which is highly responsive to negative regulation by mTOR, but under conditions of mTOR inhibition, restricts TFEB activation in a manner dependent on the proteasome. We further identify a long-term, mTOR-independent homeostatic control negatively regulating TFEB upon prolonged mTOR inhibition. These findings are of relevance for developing strategies to target TFEB activity in disease treatment. Moreover, our quantitative approach to decipher phenotype heterogeneity in imaging datasets is of general interest, as shifts between subpopulations provide a quantitative description of single cell behaviour, indicating novel regulatory behaviors and revealing differences between cell types.

A dynamic and expandable digital 3D-atlas maker for monitoring the temporal changes in tissue growth during hindbrain morphogenesis.

Reconstruction of prototypic three-dimensional (3D) atlases at the scale of whole tissues or organs requires specific methods to be developed. We have established a digital 3D-atlas maker (DAMAKER) and built a digital 3D-atlas to monitor the changes in the growth of the neuronal differentiation domain in the zebrafish hindbrain upon time. DAMAKER integrates spatial and temporal data of cell populations, neuronal differentiation and brain morphogenesis, through in vivo imaging techniques paired with image analyses and segmentation tools. First, we generated a 3D-reference from several imaged hindbrains and segmented them using a trainable tool; these were aligned using rigid registration, revealing distribution of neuronal differentiation growth patterns along the axes. Second, we quantified the dynamic growth of the neuronal differentiation domain by in vivo neuronal birthdating experiments. We generated digital neuronal birthdating 3D-maps and revealed that the temporal order of neuronal differentiation prefigured the spatial distribution of neurons in the tissue, with an inner-outer differentiation gradient. Last, we applied it to specific differentiated neuronal populations such as glutamatergic and GABAergic neurons, as proof-of-concept that the digital birthdating 3D-maps could be used as a proxy to infer neuronal birthdate. As this protocol uses open-access tools and algorithms, it can be shared for standardized, accessible, tissue-wide cell population atlas construction.

The brain, like most other organs, is formed by the coordinated growth of a few unspecialized cells in the embryo, which give rise to billions of neurons. For the brain to work properly, it is crucial that, during embryonic development, each neuron ends up in the correct location. This migration to the right spot has to happen while the brain grows and changes shape, which affects how and how far neurons and their precursor cells need to move to reach their final position. If these movements and changes in shape are not coordinated correctly, neurons can end up in the wrong place, form the wrong connections, and ultimately impact how the brain works. Previous work done in fruit flies and zebrafish resulted in three-dimensional maps of these animals' healthy brains, which allowed scientists to have a holistic view of how brains are organized. Although these maps are a valuable resource to study the structure of the brain, they do not provide information on how the brain transforms over time, especially during embryonic development. To get a clearer picture of how a few precursor cells give rise to the incredibly complex tissue that is the brain, a three-dimensional map spanning the entire developmental process is needed. To fill this gap in knowledge, Blanc et al. developed a digital atlas-maker pipeline (DAMAKER) that allows scientists to generate three-dimensional models of the embryonic brain from microscopy images of several individuals. They then used this pipeline to construct a three-dimensional digital atlas of how a part of the brain called the hindbrain develops in the zebrafish embryo. First, they collected images of the hindbrain showing neurons born at different times and matched these images to the existing static maps. Next, DAMAKER was used to follow neurons from the time of their birth to their final location, allowing Blanc et al. to create a map showing where neurons born at different stages during development end up. This type of map allows users to accurately determine when different populations of mature neurons are born, which allows scientists to estimate when different defects in brain development might originate. Based on these data, Blanc et al. concluded that in zebrafish most of the cells that will end up forming the hindbrain acquire their specialized neuronal identities very early in development, between 24 and 48 hours post fertilization. These temporal maps of healthy hindbrains were then compared to maps of brains in which the birth of neurons was disrupted, thus changing the final number of neurons in the brain. This experiment showed that changing the number of neurons that are born early in development alters the final positions of neurons and the overall shape of the brain. Therefore, for the brain to grow to its correct size, there must be a balance between the number of unspecialized cells in the developing brain, and the rate at which these cells become neurons. The DAMAKER pipeline not only provides scientists with a tool to study neurodevelopmental disorders, but also serves as a method that can be adjusted to map growth and shaping of other organs.

4D reconstruction of murine developmental trajectories using spherical harmonics.

Normal organogenesis cannot be recapitulated in vitro for mammalian organs, unlike in species including Drosophila and zebrafish. Available 3D data in the form of ex vivo images only provide discrete snapshots of the development of an organ morphology. Here, we propose a computer-based approach to recreate its continuous evolution in time and space from a set of 3D volumetric images. Our method is based on the remapping of shape data into the space of the coefficients of a spherical harmonics expansion where a smooth interpolation over time is simpler. We tested our approach on mouse limb buds and embryonic hearts. A key advantage of this method is that the resulting 4D trajectory can take advantage of all the available data while also being able to interpolate well through time intervals for which there are little or no data. This allows for a quantitative, data-driven 4D description of mouse limb morphogenesis.

Robustness of frequency vs. amplitude coding of calcium oscillations during changing temperatures.

Intracellular calcium oscillations have been widely studied. It is assumed that information is conveyed in the frequency, amplitude and shape of these oscillations. In particular, calcium signalling in mammalian liver cells has repeatedly been reported to display frequency coding so that an increasing amount of stimulus is translated into an increasing frequency of the oscillations. However, recently, we have shown that calcium oscillations in fish liver cells rather exhibit amplitude coding with increasing stimuli being translated into increasing amplitudes. Practical consequences of this difference are unknown so far. Here we investigated advantages and disadvantages of frequency vs. amplitude coding, in particular in environments with substantially changing temperatures (e.g. 10-20 degrees). For this purpose, we use computational modelling and a new approach to generate a calcium model exactly displaying a specific frequency and/or amplitude. We conclude that despite the advantages in flexibility that frequencies might offer for the transmission of information in the cell, amplitude coding is obviously more robust with respect to changes in environmental temperatures. This potentially explains the observed differences between two classes of organisms, one operating at constant temperatures whereas the other is not.

Agent-Based Modeling of Mitochondria Links Sub-Cellular Dynamics to Cellular Homeostasis and Heterogeneity.

Mitochondria are semi-autonomous organelles that supply energy for cellular biochemistry through oxidative phosphorylation. Within a cell, hundreds of mobile mitochondria undergo fusion and fission events to form a dynamic network. These morphological and mobility dynamics are essential for maintaining mitochondrial functional homeostasis, and alterations both impact and reflect cellular stress states. Mitochondrial homeostasis is further dependent on production (biogenesis) and the removal of damaged mitochondria by selective autophagy (mitophagy). While mitochondrial function, dynamics, biogenesis and mitophagy are highly-integrated processes, it is not fully understood how systemic control in the cell is established to maintain homeostasis, or respond to bioenergetic demands. Here we used agent-based modeling (ABM) to integrate molecular and imaging knowledge sets, and simulate population dynamics of mitochondria and their response to environmental energy demand. Using high-dimensional parameter searches we integrated experimentally-measured rates of mitochondrial biogenesis and mitophagy, and using sensitivity analysis we identified parameter influences on population homeostasis. By studying the dynamics of cellular subpopulations with distinct mitochondrial masses, our approach uncovered system properties of mitochondrial populations: (1) mitochondrial fusion and fission activities rapidly establish mitochondrial sub-population homeostasis, and total cellular levels of mitochondria alter fusion and fission activities and subpopulation distributions; (2) restricting the directionality of mitochondrial mobility does not alter morphology subpopulation distributions, but increases network transmission dynamics; and (3) maintaining mitochondrial mass homeostasis and responding to bioenergetic stress requires the integration of mitochondrial dynamics with the cellular bioenergetic state. Finally, (4) our model suggests sources of, and stress conditions amplifying, cell-to-cell variability of mitochondrial morphology and energetic stress states. Overall, our modeling approach integrates biochemical and imaging knowledge, and presents a novel open-modeling approach to investigate how spatial and temporal mitochondrial dynamics contribute to functional homeostasis, and how subcellular organelle heterogeneity contributes to the emergence of cell heterogeneity.