Simulations demonstrate a simple network to be sufficient to control branch point selection, smooth muscle and vasculature formation during lung branching morphogenesis

Simple Models of Lung Development

Models are essential to further our understanding of lung development and regeneration and to facilitate identification and testing of potential treatments for lung diseases. A wide variety of rodent and human models are available that recapitulate one or more stages of lung development. This chapter describes the existing 'simple' in vitro, in silico and ex vivo models of lung development. We define which stage(s) of development each model recapitulates and highlight their pros and cons.

Embryology and anatomy of congenital diaphragmatic hernia

Prenatal and postnatal treatment modalities for congenital diaphragmatic hernia (CDH) continue to improve, however patients still face high rates of morbidity and mortality caused by severe underlying persistent pulmonary hypertension and pulmonary hypoplasia. Though the majority of CDH cases are idiopathic, it is believed that CDH is a polygenic developmental defect caused by interactions between candidate genes, as well as environmental and epigenetic factors. However, the origin and pathogenesis of these developmental insults are poorly understood. Further, connections between disrupted lung development and the failure of diaphragmatic closure during embryogenesis have not been fully elucidated. Though several animal models have been useful in identifying candidate genes and disrupted signalling pathways, more studies are required to understand the pathogenesis and to develop effective preventative care. In this article, we summarize the most recent litterature on disrupted embryological lung and diaphragmatic development associated with CDH.

Instability mechanisms of repelling peak solutions in a multi-variable activator-inhibitor system

We study the linear stability properties of spatially localized single- and multi-peak states generated in a subcritical Turing bifurcation in the Meinhardt model of branching. In one spatial dimension, these states are organized in a foliated snaking structure owing to peak-peak repulsion but are shown to be all linearly unstable, with the number of unstable modes increasing with the number of peaks present. Despite this, in two spatial dimensions, direct numerical simulations reveal the presence of stable single- and multi-spot states whose properties depend on the repulsion from nearby spots as well as the shape of the domain and the boundary conditions imposed thereon. Front propagation is shown to trigger the growth of new spots while destabilizing others. The results indicate that multi-variable models may support new types of behavior that are absent from typical two-variable models.

Tracheal Ring Formation

The trachea is a long tube that enables air passage between the larynx and the bronchi. C-shaped cartilage rings on the ventral side stabilise the structure. On its esophagus-facing dorsal side, deformable smooth muscle facilitates the passage of food in the esophagus. While the symmetry break along the dorsal-ventral axis is well understood, the molecular mechanism that results in the periodic Sox9 expression pattern that translates into the cartilage rings has remained elusive. Here, we review the molecular regulatory interactions that have been elucidated, and discuss possible patterning mechanisms. Understanding the principles of self-organisation is important, both to define biomedical interventions and to enable tissue engineering.

Negative Transpulmonary Pressure Disrupts Airway Morphogenesis by Suppressing Fgf10

Mechanical forces are increasingly recognized as important determinants of cell and tissue phenotype and also appear to play a critical role in organ development. During the fetal stages of lung morphogenesis, the pressure of the fluid within the lumen of the airways is higher than that within the chest cavity, resulting in a positive transpulmonary pressure. Several congenital defects decrease or reverse transpulmonary pressure across the developing airways and are associated with a reduced number of branches and a correspondingly underdeveloped lung that is insufficient for gas exchange after birth. The small size of the early pseudoglandular stage lung and its relative inaccessibility in utero have precluded experimental investigation of the effects of transpulmonary pressure on early branching morphogenesis. Here, we present a simple culture model to explore the effects of negative transpulmonary pressure on development of the embryonic airways. We found that negative transpulmonary pressure decreases branching, and that it does so in part by altering the expression of fibroblast growth factor 10 (Fgf10). The morphogenesis of lungs maintained under negative transpulmonary pressure can be rescued by supplementing the culture medium with exogenous FGF10. These data suggest that Fgf10 expression is regulated by mechanical stress in the developing airways. Understanding the mechanical signaling pathways that connect transpulmonary pressure to FGF10 can lead to the establishment of novel non-surgical approaches for ameliorating congenital lung defects.

A review of regulated self-organizing approaches for tissue regeneration

Tissue and organ regeneration is the dynamic process by which a population of cells rearranges into a specific form with specific functions. Traditional tissue regeneration utilizes tissue grafting, cell implantation, and structured scaffolds to achieve clinical efficacy. However, tissue grafting methods face a shortage of donor tissue, while cell implantation may involve leakage of the implanted cells without a supportive 3D matrix. Cell migration, proliferation, and differentiation in structured scaffolds may disorganize and frustrate the artificially pre-designed structures, and sometimes involve immunogenic reactions. To overcome this limitation, the self-organizing properties and innate regenerative capability of tissue/organism formation in the absence of guidance by structured scaffolds has been investigated. This review emphasizes the growing subfield of the regulated self-organizing approach for neotissue formation and describes advances in the subfield using diverse, cutting-edge, inter-disciplinarity technologies. We cohesively summarize the directed self-organization of cells in the micro-engineered cell-ECM system and 3D/4D cell printing. Mathematical modeling of cellular self-organization is also discussed for providing rational guidance to intractable problems in tissue regeneration. It is envisioned that future self-organization approaches integrating biomathematics, micro-nano engineering, and gene circuits developed from synthetic biology will continue to work in concert with self-organizing morphogenesis to enhance rational control during self-organizing in tissue and organ regeneration.

Organ-Specific Branching Morphogenesis

A common developmental process, called branching morphogenesis, generates the epithelial trees in a variety of organs, including the lungs, kidneys, and glands. How branching morphogenesis can create epithelial architectures of very different shapes and functions remains elusive. In this review, we compare branching morphogenesis and its regulation in lungs and kidneys and discuss the role of signaling pathways, the mesenchyme, the extracellular matrix, and the cytoskeleton as potential organ-specific determinants of branch position, orientation, and shape. Identifying the determinants of branch and organ shape and their adaptation in different organs may reveal how a highly conserved developmental process can be adapted to different structural and functional frameworks and should provide important insights into epithelial morphogenesis and developmental disorders.

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.

The control of lung branching morphogenesis

Branching morphogenesis generates epithelial trees which facilitate gas exchange, filtering, as well as secretion processes with their large surface to volume ratio. In this review, we focus on the developmental mechanisms that control the early stages of lung branching morphogenesis. Lung branching morphogenesis involves the stereotypic, recurrent definition of new branch points, subsequent epithelial budding, and lung tube elongation. We discuss current models and experimental evidence for each of these steps. Finally, we discuss the role of the mesenchyme in determining the organ-specific shape.

Perturbation analysis of a multi-morphogen Turing reaction-diffusion stripe patterning system reveals key regulatory interactions

Periodic patterning is widespread in development and can be modelled by reaction-diffusion (RD) processes. However, minimal two-component RD descriptions are vastly simpler than the multi-molecular events that actually occur and are often hard to relate to real interactions measured experimentally. Addressing these issues, we investigated the periodic striped patterning of the rugae (transverse ridges) in the mammalian oral palate, focusing on multiple previously implicated pathways: FGF, Hh, Wnt and BMP. For each, we experimentally identified spatial patterns of activity and distinct responses of the system to inhibition. Through numerical and analytical approaches, we were able to constrain substantially the number of network structures consistent with the data. Determination of the dynamics of pattern appearance further revealed its initiation by 'activators' FGF and Wnt, and 'inhibitor' Hh, whereas BMP and mesenchyme-specific-FGF signalling were incorporated once stripes were formed. This further limited the number of possible networks. Experimental constraint thus limited the number of possible minimal networks to 154, just 0.004% of the number of possible diffusion-driven instability networks. Together, these studies articulate the principles of multi-morphogen RD patterning and demonstrate the utility of perturbation analysis for constraining RD systems.This article has an associated 'The people behind the papers' interview.

Image-based modeling of kidney branching morphogenesis reveals GDNF-RET based Turing-type mechanism and pattern-modulating WNT11 feedback

Branching patterns and regulatory networks differ between branched organs. It has remained unclear whether a common regulatory mechanism exists and how organ-specific patterns can emerge. Of all previously proposed signalling-based mechanisms, only a ligand-receptor-based Turing mechanism based on FGF10 and SHH quantitatively recapitulates the lung branching patterns. We now show that a GDNF-dependent ligand-receptor-based Turing mechanism quantitatively recapitulates branching of cultured wildtype and mutant ureteric buds, and achieves similar branching patterns when directing domain outgrowth in silico. We further predict and confirm experimentally that the kidney-specific positive feedback between WNT11 and GDNF permits the dense packing of ureteric tips. We conclude that the ligand-receptor based Turing mechanism presents a common regulatory mechanism for lungs and kidneys, despite the differences in the molecular implementation. Given its flexibility and robustness, we expect that the ligand-receptor-based Turing mechanism constitutes a likely general mechanism to guide branching morphogenesis and other symmetry breaks during organogenesis.

Many organs develop through branching morphogenesis, but whether the underlying mechanisms are shared is unknown. Here, the authors show that a ligand-receptor based Turing mechanisms, similar to that observed in lung development, likely underlies branching morphogenesis of the kidney.

Mathematical Approaches of Branching Morphogenesis

Many organs require a high surface to volume ratio to properly function. Lungs and kidneys, for example, achieve this by creating highly branched tubular structures during a developmental process called branching morphogenesis. The genes that control lung and kidney branching share a similar network structure that is based on ligand-receptor reciprocal signalling interactions between the epithelium and the surrounding mesenchyme. Nevertheless, the temporal and spatial development of the branched epithelial trees differs, resulting in organs of distinct shape and size. In the embryonic lung, branching morphogenesis highly depends on FGF10 signalling, whereas GDNF is the driving morphogen in the kidney. Knockout of Fgf10 and Gdnf leads to lung and kidney agenesis, respectively. However, FGF10 plays a significant role during kidney branching and both the FGF10 and GDNF pathway converge on the transcription factors ETV4/5. Although the involved signalling proteins have been defined, the underlying mechanism that controls lung and kidney branching morphogenesis is still elusive. A wide range of modelling approaches exists that differ not only in the mathematical framework (e.g., stochastic or deterministic) but also in the spatial scale (e.g., cell or tissue level). Due to advancing imaging techniques, image-based modelling approaches have proven to be a valuable method for investigating the control of branching events with respect to organ-specific properties. Here, we review several mathematical models on lung and kidney branching morphogenesis and suggest that a ligand-receptor-based Turing model represents a potential candidate for a general but also adaptive mechanism to control branching morphogenesis during development.

FGF10 and Human Lung Disease Across the Life Spectrum

Lung diseases impact patients across the lifespan, from infants in the first minutes of life through the aged population. Congenital abnormalities of lung structure can cause lung disease at birth or make adults more susceptible to chronic disease. Continuous inhalation of atmospheric components also requires the lung to be resilient to cellular injury. Fibroblast growth factor 10 (FGF10) regulates multiple stages of structural lung morphogenesis, cellular differentiation, and the response to injury. As a driver of lung airway branching morphogenesis, FGF10 signaling defects during development lead to neonatal lung disease. Alternatively, congenital airway abnormalities attributed to FGF10 mutations increase the risk of chronic airway disease in adulthood. FGF10 also maintains progenitor cell populations in the airway and promotes alveolar type 2 cell expansion and differentiation following injury. Here we review the cellular and molecular mechanisms linking FGF10 to multiple lung diseases, from bronchopulmonary dysplasia in extremely preterm neonates, cystic fibrosis in children, and chronic adult lung disorders. Understanding the connections between FGF10 and lung diseases may lead to exciting new therapeutic strategies.

Turing Instability-Driven Biofabrication of Branching Tissue Structures: A Dynamic Simulation and Analysis Based on the Reaction⁻Diffusion Mechanism †

Four-dimensional (4D) biofabrication techniques aim to dynamically produce and control three-dimensional (3D) biological structures that would transform their shapes or functionalities with time, when a stimulus is imposed or cell post-printing self-assembly occurs. The evolution of 3D branching patterns via self-assembly of cells is critical for the 4D biofabrication of artificial organs or tissues with branched geometry. However, it is still unclear how the formation and evolution of these branching patterns are biologically encoded. Here, we study the biofabrication of lung branching structures utilizing a simulation model based on Turing instability that raises a dynamic reaction⁻diffusion (RD) process of the biomolecules and cells. The simulation model incorporates partial differential equations of four variables, describing the tempo-spatial distribution of the variables in 3D over time. The simulation results present the formation and evolution process of 3D branching patterns over time and also interpret both the behaviors of side-branching and tip-splitting as the stalk grows and the fabrication style under an external concentration gradient of morphogen, through 3D visualization. This provides a theoretical framework for rationally guiding the 4D biofabrication of lung airway grafts via cellular self-organization, which would potentially reduce the complexity of future experimental research and number of trials.

Computational models of airway branching morphogenesis

The bronchial network of the mammalian lung consists of millions of dichotomous branches arranged in a highly complex, space-filling tree. Recent computational models of branching morphogenesis in the lung have helped uncover the biological mechanisms that construct this ramified architecture. In this review, we focus on three different theoretical approaches - geometric modeling, reaction-diffusion modeling, and continuum mechanical modeling - and discuss how, taken together, these models have identified the geometric principles necessary to build an efficient bronchial network, as well as the patterning mechanisms that specify airway geometry in the developing embryo. We emphasize models that are integrated with biological experiments and suggest how recent progress in computational modeling has advanced our understanding of airway branching morphogenesis.

What We Talk About When We Talk About Evolution

Currently, the biologic sciences are a Tower of Babel, having become so highly specialized that one discipline cannot effectively communicate with another. A mechanism for evolution that integrates development and physiologic homeostasis phylogenetically has been identified—cell-cell interactions. By reducing this process to ligand-receptor interactions and their intermediate down-stream signaling partners, it is possible, for example, to envision the functional homologies between such seemingly disparate structures and functions as the lung alveolus and kidney glomerulus, the skin and brain, or the skin and lung. For example, by showing the continuum of the lung phenotype for gas exchange at the cell-molecular level, being selected for increased surface area by augmenting lung surfactant production and function in lowering surface tension, we have determined an unprecedented structural-functional continuum from proximate to ultimate causation in evolution. It is maintained that tracing the changes in structure and function that have occurred over both the short-term history of the organism (as ontogeny), and the long-term history of the organism (as phylogeny), and how the mechanisms shared in common can account for both biologic stability and novelty, will provide the key to understanding the mechanisms of evolution. We need to better understand evolution from its unicellular origins as the Big Bang of biology.

LBIBCell: a cell-based simulation environment for morphogenetic problems

MOTIVATION

The simulation of morphogenetic problems requires the simultaneous and coupled simulation of signalling and tissue dynamics. A cellular resolution of the tissue domain is important to adequately describe the impact of cell-based events, such as cell division, cell-cell interactions and spatially restricted signalling events. A tightly coupled cell-based mechano-regulatory simulation tool is therefore required.

RESULTS

We developed an open-source software framework for morphogenetic problems. The environment offers core functionalities for the tissue and signalling models. In addition, the software offers great flexibility to add custom extensions and biologically motivated processes. Cells are represented as highly resolved, massless elastic polygons; the viscous properties of the tissue are modelled by a Newtonian fluid. The Immersed Boundary method is used to model the interaction between the viscous and elastic properties of the cells, thus extending on the IBCell model. The fluid and signalling processes are solved using the Lattice Boltzmann method. As application examples we simulate signalling-dependent tissue dynamics.

AVAILABILITY AND IMPLEMENTATION

The documentation and source code are available on http://tanakas.bitbucket.org/lbibcell/index.html

Simulating tissue morphogenesis and signaling

During embryonic development tissue morphogenesis and signaling are tightly coupled. It is therefore important to simulate both tissue morphogenesis and signaling simultaneously in in silico models of developmental processes. The resolution of the processes depends on the questions of interest. As part of this chapter we introduce different descriptions of tissue morphogenesi s. In the simplest approximation tissue is a continuous domain and tissue expansion is described according to a predefined function of time (and possibly space). In a slightly more advanced version the expansion speed and direction of the tissue may depend on a signaling variable that evolves on the domain. Both versions will be referred to as "prescribed growth." Alternatively tissue can be regarded as incompressible fluid and can be described with Navier-Stokes equations. Local cell expansion, proliferation, and death are then incorporated by a source term. In other applications the cell boundaries may be important and cell-based models must be introduced. Finally, cells may move within the tissue, a process best described by agent-based models.

An interplay of geometry and signaling enables robust lung branching morphogenesis

Early branching events during lung development are stereotyped. Although key regulatory components have been defined, the branching mechanism remains elusive. We have now used a developmental series of 3D geometric datasets of mouse embryonic lungs as well as time-lapse movies of cultured lungs to obtain physiological geometries and displacement fields. We find that only a ligand-receptor-based Turing model in combination with a particular geometry effect that arises from the distinct expression domains of ligands and receptors successfully predicts the embryonic areas of outgrowth and supports robust branch outgrowth. The geometry effect alone does not support bifurcating outgrowth, while the Turing mechanism alone is not robust to noisy initial conditions. The negative feedback between the individual Turing modules formed by fibroblast growth factor 10 (FGF10) and sonic hedgehog (SHH) enlarges the parameter space for which the embryonic growth field is reproduced. We therefore propose that a signaling mechanism based on FGF10 and SHH directs outgrowth of the lung bud via a ligand-receptor-based Turing mechanism and a geometry effect.

Feedback, receptor clustering, and receptor restriction to single cells yield large Turing spaces for ligand-receptor-based Turing models

Turing mechanisms can yield a large variety of patterns from noisy, homogenous initial conditions and have been proposed as patterning mechanism for many developmental processes. However, the molecular components that give rise to Turing patterns have remained elusive, and the small size of the parameter space that permits Turing patterns to emerge makes it difficult to explain how Turing patterns could evolve. We have recently shown that Turing patterns can be obtained with a single ligand if the ligand-receptor interaction is taken into account. Here we show that the general properties of ligand-receptor systems result in very large Turing spaces. Thus, the restriction of receptors to single cells, negative feedbacks, regulatory interactions among different ligand-receptor systems, and the clustering of receptors on the cell surface all greatly enlarge the Turing space. We further show that the feedbacks that occur in the FGF10-SHH network that controls lung branching morphogenesis are sufficient to result in large Turing spaces. We conclude that the cellular restriction of receptors provides a mechanism to sufficiently increase the size of the Turing space to make the evolution of Turing patterns likely. Additional feedbacks may then have further enlarged the Turing space. Given their robustness and flexibility, we propose that receptor-ligand-based Turing mechanisms present a general mechanism for patterning in biology.

Lung development: orchestrating the generation and regeneration of a complex organ

The respiratory system, which consists of the lungs, trachea and associated vasculature, is essential for terrestrial life. In recent years, extensive progress has been made in defining the temporal progression of lung development, and this has led to exciting discoveries, including the derivation of lung epithelium from pluripotent stem cells and the discovery of developmental pathways that are targets for new therapeutics. These discoveries have also provided new insights into the regenerative capacity of the respiratory system. This Review highlights recent advances in our understanding of lung development and regeneration, which will hopefully lead to better insights into both congenital and acquired lung diseases.

Branching patterns emerge in a mathematical model of the dynamics of lung development

Recent experimental work has described an elegant pattern of branching in the development of the lung. Multiple forms of branching have been identified, including side branching and tip bifurcation. A particularly interesting feature is the phenomenon of 'orthogonal rotation of the branching plane'. The lung must fill 3D space with the essentially 2D phenomenon of branching. It accomplishes this by rotating the branching plane by 90° with each generation. The mechanisms underlying this rotation are not understood. In general, the programmes that underlie branching have been hypothetically attributed to genetic 'subroutines' under the control of a 'global master routine' to invoke particular subroutines at the proper time and location, but the mechanisms of these routines are not known. Here, we demonstrate that fundamental mechanisms, the reaction and diffusion of biochemical morphogens, can create these patterns. We used a partial differential equation model that postulates three morphogens, which we identify with specific molecules in lung development. We found that cascades of branching events, including side branching, tip splitting and orthogonal rotation of the branching plane, all emerge immediately from the model, without further assumptions. In addition, we found that one branching mode can be easily switched to another, by increasing or decreasing the values of key parameters. This shows how a 'global master routine' could work by the alteration of a single parameter. Being able to simulate cascades of branching events is necessary to understand the critical features of branching, such as orthogonal rotation of the branching plane between successive generations, and branching mode switch during lung development. Thus, our model provides a paradigm for how genes could possibly act to produce these spatial structures. Our low-dimensional model gives a qualitative understanding of how generic physiological mechanisms can produce branching phenomena, and how the system can switch from one branching pattern to another using low-dimensional 'control knobs'. The model provides a number of testable predictions, some of which have already been observed (though not explained) in experimental work.

The control of branching morphogenesis

Many organs of higher organisms are heavily branched structures and arise by an apparently similar process of branching morphogenesis. Yet the regulatory components and local interactions that have been identified differ greatly in these organs. It is an open question whether the regulatory processes work according to a common principle and how far physical and geometrical constraints determine the branching process. Here, we review the known regulatory factors and physical constraints in lung, kidney, pancreas, prostate, mammary gland and salivary gland branching morphogenesis, and describe the models that have been formulated to analyse their impacts.

Computational modelling of bovine ovarian follicle development

Background

The development of ovarian follicles hinges on the timely exposure to the appropriate combination of hormones. Follicle stimulating hormone (FSH) and luteinizing hormone (LH) are both produced in the pituitary gland and are transported via the blood circulation to the thecal layer surrounding the follicle. From there both hormones are transported into the follicle by diffusion. FSH-receptors are expressed mainly in the granulosa while LH-receptors are expressed in a gradient with highest expression in the theca. How this spatial organization is achieved is not known. Equally it is not understood whether LH and FSH trigger distinct signalling programs or whether the distinct spatial localization of their G-protein coupled receptors is sufficient to convey their distinct biological function.

Results

We have developed a data-based computational model of the spatio-temporal signalling processes within the follicle and (i) predict that FSH and LH form a gradient inside the follicle, (ii) show that the spatial distribution of FSH- and LH-receptors can arise from the well known regulatory interactions, and (iii) find that the differential activity of FSH and LH may well result from the distinct spatial localisation of their receptors, even when both receptors respond with the same intracellular signalling cascade to their ligand.

Conclusion

The model integrates the large amount of published data into a consistent framework that can now be used to better understand how observed defects translate into failed follicle maturation.

Kidney branching morphogenesis under the control of a ligand-receptor-based Turing mechanism

The main signalling proteins that control early kidney branching have been defined. Yet the underlying mechanism is still elusive. We have previously shown that a Schnakenberg-type Turing mechanism can recapitulate the branching and protein expression patterns in wild-type and mutant lungs, but it is unclear whether this mechanism would extend to other branched organs that are regulated by other proteins. Here, we show that the glial cell line-derived neurotrophic factor-RET regulatory interaction gives rise to a Schnakenberg-type Turing model that reproduces the observed budding of the ureteric bud from the Wolffian duct, its invasion into the mesenchyme and the observed branching pattern. The model also recapitulates all relevant protein expression patterns in wild-type and mutant mice. The lung and kidney models are both based on a particular receptor-ligand interaction and require (1) cooperative binding of ligand and receptor, (2) a lower diffusion coefficient for the receptor than for the ligand and (3) an increase in the receptor concentration in response to receptor-ligand binding (by enhanced transcription, more recycling or similar). These conditions are met also by other receptor-ligand systems. We propose that ligand-receptor-based Turing patterns represent a general mechanism to control branching morphogenesis and other developmental processes.