Multiscale dynamics of branching morphogenesis

A tridimensional atlas of the developing human head

The evolution and development of the head have long captivated researchers due to the crucial role of the head as the gateway for sensory stimuli and the intricate structural complexity of the head. Although significant progress has been made in understanding head development in various vertebrate species, our knowledge of early human head ontogeny remains limited. Here, we used advanced whole-mount immunostaining and 3D imaging techniques to generate a comprehensive 3D cellular atlas of human head embryogenesis. We present detailed developmental series of diverse head tissues and cell types, including muscles, vasculature, cartilage, peripheral nerves, and exocrine glands. These datasets, accessible through a dedicated web interface, provide insights into human embryogenesis. We offer perspectives on the branching morphogenesis of human exocrine glands and unknown features of the development of neurovascular and skeletomuscular structures. These insights into human embryology have important implications for understanding craniofacial defects and neurological disorders and advancing diagnostic and therapeutic strategies.

Topology recapitulates morphogenesis of neuronal dendrites

Branching allows neurons to make synaptic contacts with large numbers of other neurons, facilitating the high connectivity of nervous systems. Neuronal arbors have geometric properties such as branch lengths and diameters that are optimal in that they maximize signaling speeds while minimizing construction costs. In this work, we asked whether neuronal arbors have topological properties that may also optimize their growth or function. We discovered that for a wide range of invertebrate and vertebrate neurons the distributions of their subtree sizes follow power laws, implying that they are scale invariant. The power-law exponent distinguishes different neuronal cell types. Postsynaptic spines and branchlets perturb scale invariance. Through simulations, we show that the subtree-size distribution depends on the symmetry of the branching rules governing arbor growth and that optimal morphologies are scale invariant. Thus, the subtree-size distribution is a topological property that recapitulates the functional morphology of dendrites.

Soliton approximation in continuum models of leader-follower behavior

Complex biological processes involve collective behavior of entities (bacteria, cells, animals) over many length and time scales and can be described by discrete models that track individuals or by continuum models involving densities and fields. We consider hybrid stochastic agent-based models of branching morphogenesis and angiogenesis (new blood vessel creation from preexisting vasculature), which treat cells as individuals that are guided by underlying continuous chemical and/or mechanical fields. In these descriptions, leader (tip) cells emerge from existing branches and follower (stalk) cells build the new sprout in their wake. Vessel branching and fusion (anastomosis) occur as a result of tip and stalk cell dynamics. Coarse graining these hybrid models in appropriate limits produces continuum partial differential equations (PDEs) for endothelial cell densities that are more analytically tractable. While these models differ in nonlinearity, they produce similar equations at leading order when chemotaxis is dominant. We analyze this leading order system in a simple quasi-one-dimensional geometry and show that the numerical solution of the leading order PDE is well described by a soliton wave that evolves from vessel to source. This wave is an attractor for intermediate times until it arrives at the hypoxic region releasing the growth factor. The mathematical techniques used here thus identify common features of discrete and continuum approaches and provide insight into general biological mechanisms governing their collective dynamics.

Fourteenth Annual ENBDC Workshop: Methods in Mammary Gland Biology and Breast Cancer

The fourteenth annual workshop of the European Network for Breast Development and Cancer (ENBDC) on Methods in Mammary Gland Biology and Breast Cancer was held on April 26th - 29th in Weggis, Switzerland. For the first time, early career researchers organised and took part in an additional ECR workshop on the 26th of April, which was received with great enthusiasm. The topics of the main workshop included mammary branching and morphogenesis, novel experimental systems (model organisms), systemic influences on tumour progression and the tumour microenvironment. Novel and recent findings were shared across excellent oral and poster presentations.

Self-organized and directed branching results in optimal coverage in developing dermal lymphatic networks

Branching morphogenesis is a ubiquitous process that gives rise to high exchange surfaces in the vasculature and epithelial organs. Lymphatic capillaries form branched networks, which play a key role in the circulation of tissue fluid and immune cells. Although mouse models and correlative patient data indicate that the lymphatic capillary density directly correlates with functional output, i.e., tissue fluid drainage and trafficking efficiency of dendritic cells, the mechanisms ensuring efficient tissue coverage remain poorly understood. Here, we use the mouse ear pinna lymphatic vessel network as a model system and combine lineage-tracing, genetic perturbations, whole-organ reconstructions and theoretical modeling to show that the dermal lymphatic capillaries tile space in an optimal, space-filling manner. This coverage is achieved by two complementary mechanisms: initial tissue invasion provides a non-optimal global scaffold via self-organized branching morphogenesis, while VEGF-C dependent side-branching from existing capillaries rapidly optimizes local coverage by directionally targeting low-density regions. With these two ingredients, we show that a minimal biophysical model can reproduce quantitatively whole-network reconstructions, across development and perturbations. Our results show that lymphatic capillary networks can exploit local self-organizing mechanisms to achieve tissue-scale optimization.

Inflationary theory of branching morphogenesis in the mouse salivary gland

The mechanisms that regulate the patterning of branched epithelia remain a subject of long-standing debate. Recently, it has been proposed that the statistical organization of multiple ductal tissues can be explained through a local self-organizing principle based on the branching-annihilating random walk (BARW) in which proliferating tips drive a process of ductal elongation and stochastic bifurcation that terminates when tips encounter maturing ducts. Here, applied to mouse salivary gland, we show the BARW model struggles to explain the large-scale organization of tissue. Instead, we propose that the gland develops as a tip-driven branching-delayed random walk (BDRW). In this framework, a generalization of the BARW, tips inhibited through steric interaction with proximate ducts may continue their branching program as constraints become alleviated through the persistent expansion of the surrounding tissue. This inflationary BDRW model presents a general paradigm for branching morphogenesis when the ductal epithelium grows cooperatively with the domain into which it expands.

Deep into the niche: Deciphering local endoderm-microenvironment interactions in development, homeostasis, and disease of pancreas and intestine

Unraveling molecular and functional heterogeneity of niche cells within the developing endoderm could resolve mechanisms of tissue formation and maturation. Here, we discuss current unknowns in molecular mechanisms underlying key developmental events in pancreatic islet and intestinal epithelial formation. Recent breakthroughs in single-cell and spatial transcriptomics, paralleled with functional studies in vitro, reveal that specialized mesenchymal subtypes drive the formation and maturation of pancreatic endocrine cells and islets via local interactions with epithelium, neurons, and microvessels. Analogous to this, distinct intestinal niche cells regulate both epithelial development and homeostasis throughout life. We propose how this knowledge can be used to progress research in the human context using pluripotent stem cell-derived multilineage organoids. Overall, understanding the interactions between the multitude of microenvironmental cells and how they drive tissue development and function could help us make more therapeutically relevant in vitro models.

A Guide Toward Multi-scale and Quantitative Branching Analysis in the Mammary Gland

The mammary gland consists of a bilayered epithelial structure with an extensively branched morphology. The majority of this epithelial tree is laid down during puberty, during which actively proliferating terminal end buds repeatedly elongate and bifurcate to form the basic structure of the ductal tree. Mammary ducts consist of a basal and luminal cell layer with a multitude of identified sub-lineages within both layers. The understanding of how these different cell lineages are cooperatively driving branching morphogenesis is a problem of crossing multiple scales, as this requires information on the macroscopic branched structure of the gland, as well as data on single-cell dynamics driving the morphogenic program. Here we describe a method to combine genetic lineage tracing with whole-gland branching analysis. Quantitative data on the global organ structure can be used to derive a model for mammary gland branching morphogenesis and provide a backbone on which the dynamics of individual cell lineages can be simulated and compared to lineage-tracing approaches. Eventually, these quantitative models and experiments allow to understand the couplings between the macroscopic shape of the mammary gland and the underlying single-cell dynamics driving branching morphogenesis.

Matrix viscoelasticity controls spatiotemporal tissue organization

Biomolecular and physical cues of the extracellular matrix environment regulate collective cell dynamics and tissue patterning. Nonetheless, how the viscoelastic properties of the matrix regulate collective cell spatial and temporal organization is not fully understood. Here we show that the passive viscoelastic properties of the matrix encapsulating a spheroidal tissue of breast epithelial cells guide tissue proliferation in space and in time. Matrix viscoelasticity prompts symmetry breaking of the spheroid, leading to the formation of invading finger-like protrusions, YAP nuclear translocation and epithelial-to-mesenchymal transition both in vitro and in vivo in a Arp2/3-complex-dependent manner. Computational modelling of these observations allows us to establish a phase diagram relating morphological stability with matrix viscoelasticity, tissue viscosity, cell motility and cell division rate, which is experimentally validated by biochemical assays and in vitro experiments with an intestinal organoid. Altogether, this work highlights the role of stress relaxation mechanisms in tissue growth dynamics, a fundamental process in morphogenesis and oncogenesis.

Human branching cholangiocyte organoids recapitulate functional bile duct formation

Human cholangiocyte organoids show great promise for regenerative therapies and in vitro modeling of bile duct development and diseases. However, the cystic organoids lack the branching morphology of intrahepatic bile ducts (IHBDs). Here, we report establishing human branching cholangiocyte organoid (BRCO) cultures. BRCOs self-organize into complex tubular structures resembling the IHBD architecture. Single-cell transcriptomics and functional analysis showed high similarity to primary cholangiocytes, and importantly, the branching growth mimics aspects of tubular development and is dependent on JAG1/NOTCH2 signaling. When applied to cholangiocarcinoma tumor organoids, the morphology changes to an in vitro morphology like primary tumors. Moreover, these branching cholangiocarcinoma organoids (BRCCAOs) better match the transcriptomic profile of primary tumors and showed increased chemoresistance to gemcitabine and cisplatin. In conclusion, BRCOs recapitulate a complex process of branching morphogenesis in vitro. This provides an improved model to study tubular formation, bile duct functionality, and associated biliary diseases.

Theory of branching morphogenesis by local interactions and global guidance

Branching morphogenesis governs the formation of many organs such as lung, kidney, and the neurovascular system. Many studies have explored system-specific molecular and cellular regulatory mechanisms, as well as self-organizing rules underlying branching morphogenesis. However, in addition to local cues, branched tissue growth can also be influenced by global guidance. Here, we develop a theoretical framework for a stochastic self-organized branching process in the presence of external cues. Combining analytical theory with numerical simulations, we predict differential signatures of global vs. local regulatory mechanisms on the branching pattern, such as angle distributions, domain size, and space-filling efficiency. We find that branch alignment follows a generic scaling law determined by the strength of global guidance, while local interactions influence the tissue density but not its overall territory. Finally, using zebrafish innervation as a model system, we test these key features of the model experimentally. Our work thus provides quantitative predictions to disentangle the role of different types of cues in shaping branched structures across scales.

Mechanical Pressure Driving Proteoglycan Expression in Mammographic Density: a Self-perpetuating Cycle?

Regions of high mammographic density (MD) in the breast are characterised by a proteoglycan (PG)-rich fibrous stroma, where PGs mediate aligned collagen fibrils to control tissue stiffness and hence the response to mechanical forces. Literature is accumulating to support the notion that mechanical stiffness may drive PG synthesis in the breast contributing to MD. We review emerging patterns in MD and other biological settings, of a positive feedback cycle of force promoting PG synthesis, such as in articular cartilage, due to increased pressure on weight bearing joints. Furthermore, we present evidence to suggest a pro-tumorigenic effect of increased mechanical force on epithelial cells in contexts where PG-mediated, aligned collagen fibrous tissue abounds, with implications for breast cancer development attributable to high MD. Finally, we summarise means through which this positive feedback mechanism of PG synthesis may be intercepted to reduce mechanical force within tissues and thus reduce disease burden.

Emerging roles of MT-MMPs in embryonic development

Membrane-type matrix metalloproteinases (MT-MMPs) are cell membrane-tethered proteinases that belong to the family of the MMPs. Apart from their roles in degradation of the extracellular milieu, MT-MMPs are able to activate through proteolytic processing at the cell surface distinct molecules such as receptors, growth factors, cytokines, adhesion molecules, and other pericellular proteins. Although most of the information regarding these enzymes comes from cancer studies, our current knowledge about their contribution in distinct developmental processes occurring in the embryo is limited. In this review, we want to summarize the involvement of MT-MMPs in distinct processes during embryonic morphogenesis, including cell migration and proliferation, epithelial-mesenchymal transition, cell polarity and branching, axon growth and navigation, synapse formation, and angiogenesis. We also considered information about MT-MMP functions from studies assessed in pathological conditions and compared these data with those relevant for embryonic development.

The nonlinear initiation of side-branching by activator-inhibitor-substrate (Turing) morphogenesis

An understanding of the underlying mechanism of side-branching is paramount in controlling and/or therapeutically treating mammalian organs, such as lungs, kidneys, and glands. Motivated by an activator-inhibitor-substrate approach that is conjectured to dominate the initiation of side-branching in a pulmonary vascular pattern, I demonstrate a distinct transverse front instability in which new fingers grow out of an oscillatory breakup dynamics at the front line without any typical length scale. These two features are attributed to unstable peak solutions in 1D that subcritically emanate from Turing bifurcation and that exhibit repulsive interactions. The results are based on a bifurcation analysis and numerical simulations and provide a potential strategy toward also developing a framework of side-branching for other biological systems, such as plant roots and cellular protrusions.

An Intravital Microscopy Toolbox to Study Mammary Gland Dynamics from Cellular Level to Organ Scale

The architecture of the mouse mammary gland is highly dynamic and constantly remodeled during pubertal development and estrous cycle-driven sprouting and regression of alveolar side branches. During each of these developmental stages, turnover is driven by distinct subsets of mammary epithelial cells. Extensive previous research has shed light on the unique morphological and cell biological characteristics of each stage. However, technological shortcomings failed to capture the dynamics and single-cell contributions to mammary remodeling. Here, we developed in vivo imaging strategies to follow the same mammary ducts over time and quantify the dynamics of mammary gland growth and remodeling from single-cell level to organ scale. Using a combination of intravital microscopy and genetic reporter systems we show how proliferative heterogeneity drives ductal morphogenesis during different developmental stages. To visualize pubertal growth at the cellular level, we performed long-term time-lapse imaging of extending terminal end buds through a mammary imaging window. We show that single-cells within the terminal end buds are extremely motile and continuously exchange position whilst the duct is elongating. To visualize short-term remodeling in the adult mammary gland at the single cell level, we performed multi-day intravital imaging in photoconvertible Kikume Green-Red mice and fluorescent ubiquitination-based cell cycle indicator mice. We demonstrate that the contribution of single-cells to estrous-driven remodeling is highly variable between cells in the same micro-environment. To assess the effects of this dynamic proliferative contribution on the long-term stability of tissue architecture, we developed a repeated skin flap method to assess mammary gland morphology by intravital microscopy over extended time spans for up to six months. Interestingly, in contrast to the short-term dynamic remodeling, the long-term morphology of the mammary gland remains remarkably stable. Together, our tool box of imaging strategies allows to identify and map transient and continuing dynamics of single cells to the architecture of the mammary gland.

Multi-scale simulation of early kidney branching morphogenesis

An important feature of the branch morphogenesis during kidney development is the termination of the tips on the outer surface of a kidney. This feature requires the avoidance of the intersection between the tips and existing ducts inside the kidney. Here, we started from a continuous model and implemented the coarse grained rules into a fast and discrete simulations. The ligand-receptor-based Turing mechanism suggests a repulsion that decreases exponentially with distance between interacting branches, preventing the intersection between neighboring branches. We considered this repulsive effect in numerical simulations and successfully reproduce the key features of the experimentally observed branch morphology for an E15.5 kidney. We examine the similarity of several geometrical parameters between the simulation results and experimental observations. The good agreement between the simulations and experiments suggests that the concentration decay caused by the absorption of glial cell line derived neurotrophic factor might be the key factor to affect the geometry in early kidney development.

Evaluating snail-trail frameworks for leader-follower behavior with agent-based modeling

Branched networks constitute a ubiquitous structure in biology, arising in plants, lungs, and the circulatory system; however, the mechanisms behind their creation are not well understood. A commonly used model for network morphogenesis proposes that sprouts develop through interactions between leader (tip) cells and follower (stalk) cells. In this description, tip cells emerge from existing structures, travel up chemoattractant gradients, and form new networks by guiding the movement of stalk cells. Such dynamics have been mathematically represented by continuum "snail-trail" models in which the tip cell flux contributes to the stalk cell proliferation rate. Although snail-trail models constitute a classical depiction of leader-follower behavior, their accuracy has yet to be evaluated in a rigorous quantitative setting. Here, we extend the snail-trail modeling framework to two spatial dimensions by introducing a novel multiplicative factor to the stalk cell rate equation, which corrects for neglected network creation in directions other than that of the migrating front. Our derivation of this factor demonstrates that snail-trail models are valid descriptions of cell dynamics when chemotaxis dominates cell movement. We confirm that our snail-trail model accurately predicts the dynamics of tip and stalk cells in an existing agent-based model (ABM) for network formation [Pillay et al., Phys. Rev. E 95, 012410 (2017)10.1103/PhysRevE.95.012410]. We also derive conditions for which it is appropriate to use a reduced, one-dimensional snail-trail model to analyze ABM results. Our analysis identifies key metrics for cell migration that may be used to anticipate when simple snail-trail models will accurately describe experimentally observed cell dynamics in network formation.

Exploring the interdependence between self-organization and functional morphology in cellular systems

All living matter is subject to continuous adaptation and functional optimization via natural selection. Consequentially, structures with close morphological resemblance repeatedly appear across the phylogenetic tree. How these designs emerge at the cellular level is not fully understood. Here, we explore core concepts of functional morphology and discuss its cause and consequences, with a specific focus on emerging properties of self-organizing systems as the potential driving force. We conclude with open questions and limitations that are present when studying shape-function interdependence in single cells and cellular ensembles.