Cell proliferation and Notch signaling coordinate the formation of epithelial folds in the Drosophila leg.

The formation of complex three-dimensional organs during development requires precise coordination between patterning networks and mechanical forces. In particular, tissue folding is a crucial process that relies on a combination of local and tissue-wide mechanical forces. Here, we investigate the contribution of cell proliferation to epithelial morphogenesis using the Drosophila leg tarsal folds as a model. We reveal that tissue-wide compression forces generated by cell proliferation, in coordination with the Notch signaling pathway, are essential for the formation of epithelial folds in precise locations along the proximo-distal axis of the leg. As cell numbers increase, compressive stresses arise, promoting the folding of the epithelium and reinforcing the apical constriction of invaginating cells. Additionally, the Notch target dysfusion plays a key function specifying the location of the folds, through the apical accumulation of F-actin and the apico-basal shortening of invaginating cells. These findings provide new insights into the intricate mechanisms involved in epithelial morphogenesis, highlighting the crucial role of tissue-wide forces in shaping a three-dimensional organ in a reproducible manner.

An ETFDH-driven metabolon supports OXPHOS efficiency in skeletal muscle by regulating coenzyme Q homeostasis.

Coenzyme Q (Q) is a key lipid electron transporter, but several aspects of its biosynthesis and redox homeostasis remain undefined. Various flavoproteins reduce ubiquinone (oxidized form of Q) to ubiquinol (QH2); however, in eukaryotes, only oxidative phosphorylation (OXPHOS) complex III (CIII) oxidizes QH2 to Q. The mechanism of action of CIII is still debated. Herein, we show that the Q reductase electron-transfer flavoprotein dehydrogenase (ETFDH) is essential for CIII activity in skeletal muscle. We identify a complex (comprising ETFDH, CIII and the Q-biosynthesis regulator COQ2) that directs electrons from lipid substrates to the respiratory chain, thereby reducing electron leaks and reactive oxygen species production. This metabolon maintains total Q levels, minimizes QH2-reductive stress and improves OXPHOS efficiency. Muscle-specific Etfdh-/- mice develop myopathy due to CIII dysfunction, indicating that ETFDH is a required OXPHOS component and a potential therapeutic target for mitochondrial redox medicine.

The effects of Hh morphogen source movement on signaling dynamics.

Morphogens of the Hh family trigger gene expression changes in receiving cells in a concentration-dependent manner to regulate their identity, proliferation, death or metabolism, depending on the tissue or organ. This variety of responses relies on a conserved signaling pathway. Its logic includes a negative-feedback loop involving the Hh receptor Ptc. Here, using experiments and computational models we study and compare the different spatial signaling profiles downstream of Hh in several developing Drosophila organs. We show that the spatial distributions of Ptc and the activator transcription factor CiA in wing, antenna and ocellus show similar features, but are markedly different from that in the compound eye. We propose that these two profile types represent two time points along the signaling dynamics, and that the interplay between the spatial displacement of the Hh source in the compound eye and the negative-feedback loop maintains the receiving cells effectively in an earlier stage of signaling. These results show how the interaction between spatial and temporal dynamics of signaling and differentiation processes may contribute to the informational versatility of the conserved Hh signaling pathway.

A coarse-grained approach to model the dynamics of the actomyosin cortex.

BACKGROUND: The dynamics of the actomyosin machinery is at the core of many important biological processes. Several relevant cellular responses such as the rhythmic compression of the cell cortex are governed, at a mesoscopic level, by the nonlinear interaction between actin monomers, actin crosslinkers, and myosin motors. Coarse-grained models are an optimal tool to study actomyosin systems, since they can include processes that occur at long time and space scales, while maintaining the most relevant features of the molecular interactions. RESULTS: Here, we present a coarse-grained model of a two-dimensional actomyosin cortex, adjacent to a three-dimensional cytoplasm. Our simplified model incorporates only well-characterized interactions between actin monomers, actin crosslinkers and myosin, and it is able to reproduce many of the most important aspects of actin filament and actomyosin network formation, such as dynamics of polymerization and depolymerization, treadmilling, network formation, and the autonomous oscillatory dynamics of actomyosin. CONCLUSIONS: We believe that the present model can be used to study the in vivo response of actomyosin networks to changes in key parameters of the system, such as alterations in the attachment of actin filaments to the cell cortex.

Quantitative Approaches to Study Retinal Neurogenesis.

The study of the development of the vertebrate retina can be addressed from several perspectives: from a purely qualitative to a more quantitative approach that takes into account its spatio-temporal features, its three-dimensional structure and also the regulation and properties at the systems level. Here, we review the ongoing transition toward a full four-dimensional characterization of the developing vertebrate retina, focusing on the challenges at the experimental, image acquisition, image processing and quantification. Using the developing zebrafish retina, we illustrate how quantitative data extracted from these type of highly dense, three-dimensional tissues depend strongly on the image quality, image processing and algorithms used to segment and quantify. Therefore, we propose that the scientific community that focuses on developmental systems could strongly benefit from a more detailed disclosure of the tools and pipelines used to process and analyze images from biological samples.

FGF2 modulates simultaneously the mode, the rate of division and the growth fraction in cultures of radial glia.

Radial glial progenitors in the mammalian developing neocortex have been shown to follow a deterministic differentiation program restricted to an asymmetric-only mode of division. This feature seems incompatible with their well-known ability to increase in number when cultured in vitro, driven by fibroblast growth factor 2 and other mitogenic signals. The changes in their differentiation dynamics that allow this transition from in vivo asymmetric-only division mode to an in vitro self-renewing culture have not been fully characterized. Here, we combine experiments of radial glia cultures with numerical models and a branching process theoretical formalism to show that fibroblast growth factor 2 has a triple effect by simultaneously increasing the growth fraction, promoting symmetric divisions and shortening the length of the cell cycle. These combined effects partner to establish and sustain a pool of rapidly proliferating radial glial progenitors in vitro We also show that, in conditions of variable proliferation dynamics, the branching process tool outperforms other commonly used methods based on thymidine analogs, such as BrdU and EdU, in terms of accuracy and reliability.

Midbody Remnant Inheritance Is Regulated by the ESCRT Subunit CHMP4C.

The inheritance of the midbody remnant (MBR) breaks the symmetry of the two daughter cells, with functional consequences for lumen and primary cilium formation by polarized epithelial cells, and also for development and differentiation. However, despite its importance, neither the relationship between the plasma membrane and the inherited MBR nor the mechanism of MBR inheritance is well known. Here, the analysis by correlative light and ultra-high-resolution scanning electron microscopy reveals a membranous stalk that physically connects the MBR to the apical membrane of epithelial cells. The stalk, which derives from the uncleaved side of the midbody, concentrates the ESCRT machinery. The ESCRT CHMP4C subunit enables MBR inheritance, and its depletion dramatically reduces the percentage of ciliated cells. We demonstrate (1) that MBRs are physically connected to the plasma membrane, (2) how CHMP4C helps maintain the integrity of the connection, and (3) the functional importance of the connection.

Human-iPSC-Derived Cardiac Stromal Cells Enhance Maturation in 3D Cardiac Microtissues and Reveal Non-cardiomyocyte Contributions to Heart Disease.

Cardiomyocytes (CMs) from human induced pluripotent stem cells (hiPSCs) are functionally immature, but this is improved by incorporation into engineered tissues or forced contraction. Here, we showed that tri-cellular combinations of hiPSC-derived CMs, cardiac fibroblasts (CFs), and cardiac endothelial cells also enhance maturation in easily constructed, scaffold-free, three-dimensional microtissues (MTs). hiPSC-CMs in MTs with CFs showed improved sarcomeric structures with T-tubules, enhanced contractility, and mitochondrial respiration and were electrophysiologically more mature than MTs without CFs. Interactions mediating maturation included coupling between hiPSC-CMs and CFs through connexin 43 (CX43) gap junctions and increased intracellular cyclic AMP (cAMP). Scaled production of thousands of hiPSC-MTs was highly reproducible across lines and differentiated cell batches. MTs containing healthy-control hiPSC-CMs but hiPSC-CFs from patients with arrhythmogenic cardiomyopathy strikingly recapitulated features of the disease. Our MT model is thus a simple and versatile platform for modeling multicellular cardiac diseases that will facilitate industry and academic engagement in high-throughput molecular screening.

Control of size, fate and time by the Hh morphogen in the eyes of flies.

Molecules of the hedgehog (hh) family are involved in the specification and patterning of eyes in vertebrates and invertebrates. These organs, though, are of very different sizes, raising the question of how Hh molecules operate at such different scales. In this paper we discuss the strategies used by Hh to control the development of the two eye types in Drosophila: the large compound eye and the small ocellus. We first describe the distinct ways in which these two eyes develop and the evidence for the key role played by Hh in both; then we consider the potential for variation in the range of action of a "typical" morphogen and measure this range ("characteristic length") for Hh in different organs, including the compound eye and the ocellus. Finally, we describe how different feedback mechanisms are used to extend the Hh range of action to pattern the large and even the small eye. In the ocellus, the basic Hh signaling pathway adds to its dynamics the attenuation of its receptor as cell differentiate. This sole regulatory change can result in the decoding of the Hh gradient by receiving cells as a wave of constant speed. Therefore, in the fly ocellus, the Hh morphogen adds to its spatial patterning role a novel one: patterning along a time axis.

Dynamic Hh signalling can generate temporal information during tissue patterning.

The differentiation of tissues and organs requires that cells exchange information in space and time. Spatial information is often conveyed by morphogens: molecules that disperse across receiving cells to generate signalling gradients. Cells translate such concentration gradients into space-dependent patterns of gene expression and cellular behaviour. But could morphogen gradients also convey developmental time? Here, by investigating the developmental role of Hh on a component of the Drosophila visual system, the ocellar retina, we have discovered that ocellar cells use the non-linear gradient of Hh as a temporal cue, collectively performing the biological equivalent of a mathematical logarithmic transformation. In this way, a morphogen diffusing from a non-moving source is decoded as a wave of differentiating photoreceptors that travels at constant speed throughout the retinal epithelium.

Drug treatment efficiency depends on the initial state of activation in nonlinear pathways.

An accurate prediction of the outcome of a given drug treatment requires quantitative values for all parameters and concentrations involved as well as a detailed characterization of the network of interactions where the target molecule is embedded. Here, we present a high-throughput in silico screening of all potential networks of three interacting nodes to study the effect of the initial conditions of the network in the efficiency of drug inhibition. Our study shows that most network topologies can induce multiple dose-response curves, where the treatment has an enhanced, reduced or even no effect depending on the initial conditions. The type of dual response observed depends on how the potential bistable regimes interplay with the inhibition of one of the nodes inside a nonlinear pathway architecture. We propose that this dependence of the strength of the drug on the initial state of activation of the pathway may be affecting the outcome and the reproducibility of drug studies and clinical trials.

Novel role for the midbody in primary ciliogenesis by polarized epithelial cells.

The primary cilium is a membrane protrusion that is crucial for vertebrate tissue homeostasis and development. Here, we investigated the uncharacterized process of primary ciliogenesis in polarized epithelial cells. We show that after cytokinesis, the midbody is inherited by one of the daughter cells as a remnant that initially locates peripherally at the apical surface of one of the daughter cells. The remnant then moves along the apical surface and, once proximal to the centrosome at the center of the apical surface, enables cilium formation. The physical removal of the remnant greatly impairs ciliogenesis. We developed a probabilistic cell population-based model that reproduces the experimental data. In addition, our model explains, solely in terms of cell area constraints, the various observed transitions of the midbody, the beginning of ciliogenesis, and the accumulation of ciliated cells. Our findings reveal a biological mechanism that links the three microtubule-based organelles-the midbody, the centrosome, and the cilium-in the same cellular process.

Myosin light-chain phosphatase regulates basal actomyosin oscillations during morphogenesis.

Contractile actomyosin networks generate forces that drive tissue morphogenesis. Actomyosin contractility is controlled primarily by reversible phosphorylation of the myosin-II regulatory light chain through the action of myosin kinases and phosphatases. While the role of myosin light-chain kinase in regulating contractility during morphogenesis has been largely characterized, there is surprisingly little information on myosin light-chain phosphatase (MLCP) function in this context. Here, we use live imaging of Drosophila follicle cells combined with mathematical modelling to demonstrate that the MLCP subunit flapwing (flw) is a key regulator of basal myosin oscillations and cell contractions underlying egg chamber elongation. Flw expression decreases specifically on the basal side of follicle cells at the onset of contraction and flw controls the initiation and periodicity of basal actomyosin oscillations. Contrary to previous reports, basal F-actin pulsates similarly to myosin. Finally, we propose a quantitative model in which periodic basal actomyosin oscillations arise in a cell-autonomous fashion from intrinsic properties of motor assemblies.

A Branching Process to Characterize the Dynamics of Stem Cell Differentiation.

The understanding of the regulatory processes that orchestrate stem cell maintenance is a cornerstone in developmental biology. Here, we present a mathematical model based on a branching process formalism that predicts average rates of proliferative and differentiative divisions in a given stem cell population. In the context of vertebrate neurogenesis, the model predicts complex non-monotonic variations in the rates of pp, pd and dd modes of division as well as in cell cycle length, in agreement with experimental results. Moreover, the model shows that the differentiation probability follows a binomial distribution, allowing us to develop equations to predict the rates of each mode of division. A phenomenological simulation of the developing spinal cord informed with the average cell cycle length and division rates predicted by the mathematical model reproduces the correct dynamics of proliferation and differentiation in terms of average numbers of progenitors and differentiated cells. Overall, the present mathematical framework represents a powerful tool to unveil the changes in the rate and mode of division of a given stem cell pool by simply quantifying numbers of cells at different times.

Synergistic interaction between selective drugs in cell populations models.

The design of selective drugs and combinatorial drug treatments are two of the main focuses in modern pharmacology. In this study we use a mathematical model of chimeric ligand-receptor interaction to show that the combination of selective drugs is synergistic in nature, providing a way to gain optimal selective potential at reduced doses compared to the same drugs when applied individually. We use a cell population model of proliferating cells expressing two different amounts of a target protein to show that both selectivity and synergism are robust against variability and heritability in the cell population. The reduction in the total drug administered due to the synergistic performance of the selective drugs can potentially result in reduced toxicity and off-target interactions, providing a mechanism to improve the treatment of cell-based diseases caused by aberrant gene overexpression, such as cancer and diabetes.

A tunable coarse-grained model for ligand-receptor interaction.

Cell-surface receptors are the most common target for therapeutic drugs. The design and optimization of next generation synthetic drugs require a detailed understanding of the interaction with their corresponding receptors. Mathematical approximations to study ligand-receptor systems based on reaction kinetics strongly simplify the spatial constraints of the interaction, while full atomistic ligand-receptor models do not allow for a statistical many-particle analysis, due to their high computational requirements. Here we present a generic coarse-grained model for ligand-receptor systems that accounts for the essential spatial characteristics of the interaction, while allowing statistical analysis. The model captures the main features of ligand-receptor kinetics, such as diffusion dependence of affinity and dissociation rates. Our model is used to characterize chimeric compounds, designed to take advantage of the receptor over-expression phenotype of certain diseases to selectively target unhealthy cells. Molecular dynamics simulations of chimeric ligands are used to study how selectivity can be optimized based on receptor abundance, ligand-receptor affinity and length of the linker between both ligand subunits. Overall, this coarse-grained model is a useful approximation in the study of systems with complex ligand-receptor interactions or spatial constraints.

Smad2 and Smad3 cooperate and antagonize simultaneously in vertebrate neurogenesis.

The transforming growth factor beta (TGF-ß) pathway plays key roles in development and cancer. TGF-ß signaling converges on the Smad2 and Smad3 effectors, which can either cooperate or antagonize to regulate their transcriptional targets. Here we performed in vivo and in silico experiments to study how such cooperativity and antagonism might function during neurogenesis. In vivo electroporation experiments in the chick embryo neural tube show that Smad2 and Smad3 cooperate to promote neurogenesis, as well as the transcription of Smad3-specific targets. Knockdown of Smad2 enhances neurogenesis and the transcription of Smad3-specific targets. A mathematical model of the TGF-ß pathway fits the experimental results and predicts that the proportions of the three different trimeric complexes formed dictates the transcriptional responses of the R-Smad proteins. As such, Smad2 targets are activated solely by the Smad2-Smad2-Smad4 complex, whereas Smad3 targets are activated both by Smad2-Smad3-Smad4 and Smad3-Smad3-Smad4 trimers. We have modeled the Smad responses onto arbitrary genes and propose that this mechanism might be extended to additional activities of TGF-ß in development and disease.

Synthetic approaches to study transcriptional networks and noise in mammalian systems.

Synthetic biology aims to build new functional organisms and to rationally re-design existing ones by applying the engineering principle of modularity. Apart from building new life forms to perform technical applications, the approach of synthetic biology is useful to dissect complex biological phenomena into simple and easy to understand synthetic modules. Synthetic gene networks have been successfully implemented in prokaryotes and lower eukaryotes, with recent approaches moving ahead towards the mammalian environment. However, synthetic circuits in higher eukaryotes present a more challenging scenario, since its reliability is compromised because of the strong stochastic nature of transcription. Here, the authors review recent approaches that take advantage of the noisy response of synthetic regulatory circuits to learn key features of the complex machinery that orchestrates transcription in higher eukaryotes. Understanding the causes and consequences of biological noise will allow us to design more reliable mammalian synthetic circuits with revolutionary medical applications.

Sonic hedgehog signaling switches the mode of division in the developing nervous system.

The different modes of stem cell division are tightly regulated to balance growth and differentiation during organ development and homeostasis, and these regulatory processes are subverted in tumor formation. Here, we developed markers that provided the single-cell resolution necessary to quantify the three modes of division taking place in the developing nervous system in vivo: self-expanding, PP; self-replacing, PN; and self-consuming, NN. Using these markers and a mathematical model that predicts the dynamics of motor neuron progenitor division, we identify a role for the morphogen Sonic hedgehog in the maintenance of stem cell identity in the developing spinal cord. Moreover, our study provides insight into the process linking lineage commitment to neurogenesis with changes in cell-cycle parameters. As a result, we propose a challenging model in which the external Sonic hedgehog signal dictates stem cell identity, reflected in the consequent readjustment of cell-cycle parameters.

Network nonlinearities in drug treatment.

Despite major achievements in the understanding of human disease, there is a general perception that the drug development industry has failed to meet the expectations that recent advances in biotechnology should drive. One of the potential sources of failure of many next generation drugs is that their targets are embedded in highly nonlinear signaling pathways and gene networks with multiple negative and positive feedback loops of regulation. There is increasing evidence that this complex network shapes the response to external perturbations in the form of drug treatment, originating bistability, hypersensitivity, robustness, complex dose-response curves or schedule dependent activity. This review focuses on the effect of nonlinearities on signaling and gene networks involved in human disease, using tools from Nonlinear Dynamics to discuss the implications and to overcome the effects of the nonlinearities on regulatory networks.