An expanded Notch-Delta model exhibiting long-range patterning and incorporating MicroRNA regulation

Expression Profiling of Adipogenic and Anti-Adipogenic MicroRNA Sequences following Methylmercury Exposure in Caenorhabditis elegans

MicroRNA (miRNA) are important regulators of gene expression that respond not only to developmental and pathological cues, but also to environmental stimuli. Dyslipidemia is a hallmark of metabolic conditions and has been shown to significantly affect the expression of circulating miRNA sequences. Recently, our lab has shown that the environmental toxicant methylmercury (MeHg) causes dyslipidemia in the Caenorhabditis elegans model organism. While 10 and 20 μM MeHg increases the expression of adipogenic transcription factors and lipid-binding proteins in worms, there is limited information on how the toxicant affects the miRNA regulators of these genes. We hypothesized that MeHg would increase the expression of adipogenic miRNA sequences and/or decrease the expression of anti-adipogenic miRNA sequences. We further hypothesized that the target mRNA sequences for the miRNAs affected by MeHg would be consequently altered. We selected three potentially adipogenic (mir-34, mir-124, and mir-355) and three potentially anti-adipogenic (mir-240, mir-786, and let-7) miRNA sequences homologous to known human miRNA sequences altered in obesity, and quantified their levels 24 h and 48 h post MeHg treatment. At 24 h post exposure, MeHg significantly increased expression of both the adipogenic and anti-adipogenic miRNA sequences 1.5-3x above untreated control. By 48 h post exposure, only the adipogenic miRNA sequences were elevated, while the anti-adipogenic miRNA sequences were decreased by 50% compared to untreated control. These data suggest that there are developmental changes in miRNA expression over time following MeHg exposure. We next selected one target mRNA sequence for each miRNA sequence based on miRNA-mRNA relationships observed in humans. MeHg altered the gene expression of all the target genes assayed. Except for mir-34, all the tested miRNA-mRNA sequences showed a conserved relationship between nematode and humans. To determine whether the selected miRNA sequences were involved in lipid accumulation in response to MeHg, lipid storage was investigated in transgenic worm strains that lacked the specific miRNA strains. Of the six strains investigated, only the mir-124 and let-7 mutant worms had lipid storage levels that were statistically different from wild type, suggesting that these two sequences can be potential mediators of MeHg-induced lipid dysregulation.

Adjusting the range of cell-cell communication enables fine-tuning of cell fate patterns from checkerboard to engulfing

During development, spatio-temporal patterns ranging from checkerboard to engulfing occur with precise proportions of the respective cell fates. Key developmental regulators are intracellular transcriptional interactions and intercellular signaling. We present an analytically tractable mathematical model based on signaling that reliably generates different cell type patterns with specified proportions. Employing statistical mechanics, We derived a cell fate decision model for two cell types. A detailed steady state analysis on the resulting dynamical system yielded necessary conditions to generate spatially heterogeneous patterns. This allows the cell type proportions to be controlled by a single model parameter. Cell-cell communication is realized by local and global signaling mechanisms. These result in different cell type patterns. A nearest neighbor signal yields checkerboard patterns. Increasing the signal dispersion, cell fate clusters and an engulfing pattern can be generated. Altogether, the presented model allows us to reliably generate heterogeneous cell type patterns of different kinds as well as desired proportions.

Coupling dynamics of 2D Notch-Delta signalling

Understanding pattern formation driven by cell-cell interactions has been a significant theme in cellular biology for many years. In particular, due to their implications within many biological contexts, lateral-inhibition mechanisms present in the Notch-Delta signalling pathway led to an extensive discussion between biologists and mathematicians. Deterministic and stochastic models have been developed as a consequence of this discussion, some of which address long-range signalling by considering cell protrusions reaching non-neighbouring cells. The dynamics of such signalling systems reveal intricate properties of the coupling terms involved in these models. In this work, we investigate the advantages and drawbacks of a single-parameter long-range signalling model across diverse scenarios. By employing linear and multi-scale analyses, we discover that pattern selection is not only partially explained but also depends on nonlinear effects that extend beyond the scope of these analytical techniques.

Delineating cooperative effects of Notch and biomechanical signals on patterned liver differentiation

Controlled in vitro multicellular culture systems with defined biophysical microenvironment have been used to elucidate the role of Notch signaling in the spatiotemporal regulation of stem and progenitor cell differentiation. In addition, computational models incorporating features of Notch ligand-receptor interactions have provided important insights into Notch pathway signaling dynamics. However, the mechanistic relationship between Notch-mediated intercellular signaling and cooperative microenvironmental cues is less clear. Here, liver progenitor cell differentiation patterning was used as a model to systematically evaluate the complex interplay of cellular mechanics and Notch signaling along with identifying combinatorial mechanisms guiding progenitor fate. We present an integrated approach that pairs a computational intercellular signaling model with defined microscale culture configurations provided within a cell microarray platform. Specifically, the cell microarray-based experiments were used to validate and optimize parameters of the intercellular Notch signaling model. This model incorporated the experimentally established multicellular dimensions of the cellular microarray domains, mechanical stress-related activation parameters, and distinct Notch receptor-ligand interactions based on the roles of the Notch ligands Jagged-1 and Delta-like-1. Overall, these studies demonstrate the spatial control of mechanotransduction-associated components, key growth factor and Notch signaling interactions, and point towards a possible role of E-Cadherin in translating intercellular mechanical gradients to downstream Notch signaling.

Stochastic fluctuations promote ordered pattern formation of cells in the Notch-Delta signaling pathway

The Notch-Delta signaling pathway mediates cell differentiation implicated in many regulatory processes including spatiotemporal patterning in tissues by promoting alternate cell fates between neighboring cells. At the multicellular level, this "lateral inhibition” principle leads to checkerboard patterns with alternation of Sender and Receiver cells. While it is well known that stochasticity modulates cell fate specification, little is known about how stochastic fluctuations at the cellular level propagate during multicell pattern formation. Here, we model stochastic fluctuations in the Notch-Delta pathway in the presence of two different noise types–shot and white–for a multicell system. Our results show that intermediate fluctuations reduce disorder and guide the multicell lattice toward checkerboard-like patterns. By further analyzing cell fate transition events, we demonstrate that intermediate noise amplitudes provide enough perturbation to facilitate “proofreading” of disordered patterns and cause cells to switch to the correct ordered state (Sender surrounded by Receivers, and vice versa). Conversely, high noise can override environmental signals coming from neighboring cells and lead to switching between ordered and disordered patterns. Therefore, in analogy with spin glass systems, intermediate noise levels allow the multicell Notch system to escape frustrated patterns and relax towards the lower energy checkerboard pattern while at large noise levels the system is unable to find this ordered base of attraction.

The Notch pathway is involved in many biological processes and is known to form precise spatial patterns alternating Sender and Receiver cell states. Quantifying the implications of stochastic fluctuations provided insight that patterns formed in Notch-mediated pathways must follow a predetermined path towards checkerboard or exist in a noisy environment which promotes order through error correction.

We model Notch pattern formation stochastically and analyze the spatiotemporal dynamics. Our results show multicellular systems equilibrate towards ordered systems, but mistakes in the initial lattice propagate causing the systems to relax into frustrated systems. Only through existing in a noisy environment are the systems able to relax into the checkerboard pattern. Analyzing the temporal dynamics confirms, in environments with intermediate noise, the “incorrect” cells (Sender in a Sender environment, and vice versa) can be flipped to the correct state (Sender in a Receiver environment, and vice versa). Comparing with the spin glass energy landscape, we suggest the multicellular model follows a rugged landscape to form patterns with stochastic fluctuations required to enforce order throughout the system.

Mechano-regulated cell-cell signaling in the context of cardiovascular tissue engineering

Cardiovascular tissue engineering (CVTE) aims to create living tissues, with the ability to grow and remodel, as replacements for diseased blood vessels and heart valves. Despite promising results, the (long-term) functionality of these engineered tissues still needs improvement to reach broad clinical application. The functionality of native tissues is ensured by their specific mechanical properties directly arising from tissue organization. We therefore hypothesize that establishing a native-like tissue organization is vital to overcome the limitations of current CVTE approaches. To achieve this aim, a better understanding of the growth and remodeling (G&R) mechanisms of cardiovascular tissues is necessary. Cells are the main mediators of tissue G&R, and their behavior is strongly influenced by both mechanical stimuli and cell-cell signaling. An increasing number of signaling pathways has also been identified as mechanosensitive. As such, they may have a key underlying role in regulating the G&R of tissues in response to mechanical stimuli. A more detailed understanding of mechano-regulated cell-cell signaling may thus be crucial to advance CVTE, as it could inspire new methods to control tissue G&R and improve the organization and functionality of engineered tissues, thereby accelerating clinical translation. In this review, we discuss the organization and biomechanics of native cardiovascular tissues; recent CVTE studies emphasizing the obtained engineered tissue organization; and the interplay between mechanical stimuli, cell behavior, and cell-cell signaling. In addition, we review past contributions of computational models in understanding and predicting mechano-regulated tissue G&R and cell-cell signaling to highlight their potential role in future CVTE strategies.

STEADY STATE AND SENSITIVITY ANALYSIS OF A NOTCH–DELTA SIGNALING SYSTEM OF ONE SINGLE CELL INTERACTING WITH FIXED ENVIRONMENT

The Notch–Delta signaling pathway is a highly conserved signaling system that partakes in a diverse process of growth, patterns and differentiation. Experiments have shown that Delta from different cells activates this pathway (trans-activation) while Delta from the same cell inhibits this pathway (cis-inhibition). The Notch–Delta interactions could switch a cell to one of the two opposite fates: either Sender (high Delta/low Notch) or Receiver (low Delta/high Notch). We studied a Notch–Delta signaling model from Sprinzak et al., (2010), to investigate the cell fate through steady state analysis. The focus was placed on a fundamental case of one single cell with fixed external Delta and Notch supplies. First, we proved there exists a unique steady state which is asymptotically stable. Second, we derived the increasing/decreasing and asymptotic properties of the steady state with respect to all the parameters. Third, we studied the sensitivity and discovered the cell fate is only sensitive to the production rates of Notch and Delta under strong cis-inhibition. Finally, we applied this model to multi-cellular cases and found that the lateral inhibition pattern could be created with the spatially varied Delta production rate. The Hopf bifurcation is not observed in the current model.

Quantitative Characteristic of ncRNA Regulation in Gene Regulatory Networks

RNA is mostly known for its role in protein synthesis, where it encodes information for protein sequence in its messenger RNA (mRNA) form (translation). Yet, RNA molecules regulate several cellular processes other than translation. Here, we present an overview of several mathematical models that help understanding and characterizing the role of noncoding RNA molecules (ncRNAs) in regulating gene expression and protein synthesis. First, we discuss relatively simple models where ncRNAs can modulate protein synthesis via targeting a mRNA. Then, we consider the case of feedback interactions between ncRNAs and their target proteins, and discuss several biological applications where these feedback architectures modulate a cellular phenotype and control the levels of intrinsic and extrinsic noise. Building from these simple circuit motifs, we examine feed-forward circuit motifs involving ncRNAs that generate precise spatial and temporal patterns of protein expression. Further, we investigate the competition between ncRNAs and other endogenous RNA molecules and show that the cross talk between coding and noncoding RNAs can form large genetic circuits that involve up to hundreds of chemical species. Finally, we discuss the role of ncRNAs in modulating cell-cell signaling pathways and therefore the dynamics of spatiotemporal pattern formation in a tissue.

Enhanced Delta-Notch Lateral Inhibition Model Incorporating Intracellular Notch Heterogeneity and Tension-Dependent Rate of Delta-Notch Binding that Reproduces Sprouting Angiogenesis Patterns

Endothelial cells adopt unique cell fates during sprouting angiogenesis, differentiating into tip or stalk cells. The fate selection process is directed by Delta-Notch lateral inhibition pathway. Classical Delta-Notch models produce a spatial pattern of tip cells separated by a single stalk cell, or the salt-and-pepper pattern. However, classical models cannot explain alternative tip-stalk patterning, such as tip cells that are separated by two or more stalk cells. We show that lateral inhibition models involving only Delta and Notch proteins can also recapitulate experimental tip-stalk patterns by invoking two mechanisms, specifically, intracellular Notch heterogeneity and tension-dependent rate of Delta-Notch binding. We introduce our computational model and analysis where we establish that our enhanced Delta-Notch lateral inhibition model can recapitulate a greater variety of tip-stalk patterning which is previously not possible using classical lateral inhibition models. In our enhanced Delta-Notch lateral inhibition model, we observe the existence of a hybrid cell type displaying intermediate tip and stalk cells’ characteristics. We validate the existence of such hybrid cells by immuno-staining of endothelial cells with tip cell markers, Delta and CD34, which substantiates our enhanced model.

Nonprotein-Coding RNAs as Regulators of Development in Tunicates

Tunicates, or urochordates, are a group of small marine organisms that are found widely throughout the seas of the world. As most plausible sister group of the vertebrates, they are of utmost importance for a comprehensive understanding of chordate evolution; hence, they have served as model organisms for many aspects of the developmental biology. Current genomic analysis of tunicates indicates that their genomes evolved with a fast rate not only at the level of nucleotide substitutions but also in terms of genomic organization. The latter involves genome reduction, rearrangements, as well as the loss of some important coding and noncoding RNA (ncRNAs) elements and even entire genomic regions that are otherwise well conserved. These observations are largely based on evidence from comparative genomics resulting from the analysis of well-studied gene families such as the Hox genes and their noncoding elements. In this chapter, the focus lies on the ncRNA complement of tunicates, with a particular emphasis on microRNAs, which have already been studied extensively for other animal clades. MicroRNAs are known as important regulators of key genes in animal development, and they are intimately related to the increase morphological complexity in higher metazoans. Here we review the discovery, evolution, and genome organization of the miRNA repertoire, which has been drastically reduced and restructured in tunicates compared to the chordate ancestor. Known functions of microRNAs as regulators of development in tunicates are a central topic. For instance, we consider the role of miRNAs as regulators of the muscle development and their importance in the regulation of the differential expression during the oral siphon regeneration. Beyond microRNAs, we touch upon the functions of some other ncRNAs such as yellow crescent RNA, moRNAs, RMST lncRNAs, or spliced-leader (SL) RNAs, which have diverse functions associated with the embryonic development, neurogenesis, and mediation of mRNA stability in general.

Role of domain in pattern formation

Pattern formation during development is one of the elegant self-organized phenomena that allow cells to regulate their functions. At all levels, from DNA to a tissue or organ, many developmental processes include the determination of cellular functions through pattern formation. To elucidate the mechanism underlying pattern formation, numerous mathematical models have been developed and applied. However, model simplification has resulted in the role of domains not being seriously considered in pattern formation. Here, we introduce a novel application of the phase-field method for analysis of chromatin dynamics, and a mathematical approach that includes domain information into a biochemical model of pattern formation. Using this new modeling method, here, we consider the role of nuclear and cellular cell shapes on pattern formation.

Lateral inhibition-induced pattern formation controlled by the size and geometry of the cell

Pattern formation in development biology is one of the fundamental processes by which cells change their functions. It is based on the communication of cells via intra- and intercellular dynamics of biochemicals. Thus, the cell is directly involved in biochemical interactions. However, many theoretical approaches describing biochemical pattern formation have usually neglected the cell's role or have simplified the subcellular process without considering cellular aspects despite the cell being the environment where biochemicals interact. On the other hand, recent experimental observations suggest that a change in the physical conditions of cell-to-cell contact can result in a change in cell fate and tissue patterning in a lateral inhibition system. Here we develop a mathematical model by which biochemical dynamics can be directly observed with explicitly expressed cell structure and geometry in higher dimensions, and reconsider pattern formation by lateral inhibition of the Notch-Delta signaling pathway. We explore how the physical characteristic of cell, such as cell geometry or size, influences the biochemical pattern formation in a multi-cellular system. Our results suggest that a property based on cell geometry can be a novel mechanism for symmetry breaking inducing cell asymmetry. We show that cell volume can critically influence cell fate determination and pattern formation at the tissue level, and the surface area of the cell-to-cell contact can directly affect the spatial range of patterning.

Pattern formation in discrete cell tissues under long range filopodia-based direct cell to cell contact

Pattern formation via direct cell to cell contact has received considerable attention over the years. In particular the lateral-inhibition mechanism observed in the Notch signalling pathway can generate a regular periodic pattern of differential cell activity, and has been proposed to explain the emergence of patterns in various tissues and organs. The majority of models of this system have focussed on short-range contacts: a cell signals only to its nearest neighbours and the resulting patterns tend to be of fine-scale "salt and pepper" nature. The capacity of certain cells to extend signalling filopodia (cytonemes) over multiple cell lengths, however, inserts a long-range or non-local component into this process. Here we explore how long range signalling can impact on pattern formation. Specifically, we extend a standard model for Notch-like lateral inhibition to include cytoneme-mediated signalling, and investigate how pattern formation depends on the spatial distribution of signal from the signalling cell. We show that a variety of patterns can be obtained, ranging from a sparse pattern of single isolated cells to larger clusters or stripes.

Development of hair cells in inner ear is associated with expression and promoter methylation of Notch-1 in postnatal mice

The present study was designed to investigate the correlation among the number of hair cells in inner ear, Notch-1 gene expression levels and its methylation status of the promoter region in the postnatal mice. The hair cells in inner ear were collected from postnatal mice at day 0, 4, 8 and 16 and counted by immunofluorescence. Notch-1 mRNA expression were measured by real-time quantitative polymerize chain reaction (PCR). Methylation levels of CpG islands in Notch-1 promoters were analyzed by matrix-assisted laser desorption/ionization-time of flight mass spectrometry. The results showed that the number of hair cells in the inner ear increased gradually after birth, which were positively correlated to Notch-1 mRNA expression. However, analysis on methylation of CpG sites in Notch-1 promoter showed that the methylation rates increased gradually after births, which were correlated with the decreased expression of Notch-1. Drug lesion induced the loss of hair cells, and stimulated the expression of Notch-1 mRNA expression, but didn't influence the methylation rates of Notch-1 promoter. We concluded that the Notch-1 mRNA expression level in inner ear tissues is correlated with the development of hair cells. CpG islands in Notch-1 promoter region manifest hypermethylation status when hair cells in inner ear are mature.

Jagged-Delta asymmetry in Notch signaling can give rise to a Sender/Receiver hybrid phenotype

Notch signaling pathway mediates cell-fate determination during embryonic development, wound healing, and tumorigenesis. This pathway is activated when the ligand Delta or the ligand Jagged of one cell interacts with the Notch receptor of its neighboring cell, releasing the Notch Intracellular Domain (NICD) that activates many downstream target genes. NICD affects ligand production asymmetrically--it represses Delta, but activates Jagged. Although the dynamical role of Notch-Jagged signaling remains elusive, it is widely recognized that Notch-Delta signaling behaves as an intercellular toggle switch, giving rise to two distinct fates that neighboring cells adopt--Sender (high ligand, low receptor) and Receiver (low ligand, high receptor). Here, we devise a specific theoretical framework that incorporates both Delta and Jagged in Notch signaling circuit to explore the functional role of Jagged in cell-fate determination. We find that the asymmetric effect of NICD renders the circuit to behave as a three-way switch, giving rise to an additional state--a hybrid Sender/Receiver (medium ligand, medium receptor). This phenotype allows neighboring cells to both send and receive signals, thereby attaining similar fates. We also show that due to the asymmetric effect of the glycosyltransferase Fringe, different outcomes are generated depending on which ligand is dominant: Delta-mediated signaling drives neighboring cells to have an opposite fate; Jagged-mediated signaling drives the cell to maintain a similar fate to that of its neighbor. We elucidate the role of Jagged in cell-fate determination and discuss its possible implications in understanding tumor-stroma cross-talk, which frequently entails Notch-Jagged communication.