Proper heart morphogenesis requires a delicate balance between hemodynamic forces, myocardial activity, morphogen gradients, and epigenetic signaling, all of which are coupled with genetic regulatory networks. Recently both in vivo and in silico studies have tried to better understand hemodynamics at varying stages of vertebrate cardiogenesis. In particular, the intracardial hemodynamics during the onset of trabeculation is notably complex - the inertial and viscous fluid forces are approximately equal at this stage and small perturbations in morphology, scale, and steadiness of the flow can lead to significant changes in bulk flow structures, shear stress distributions, and chemical morphogen gradients. The immersed boundary method was used to solve the computational fluid dynamics problem involving fluid flow moving through the trabeculated ventricles of 72, 80, and 120 hours post fertilization wild type zebrafish embryos and ErbB2-inhibited embryos at 7 days post fertilization. An idealized trabeculated ventricular model was explored to map the bifurcations in flow structure that occur as a result of the unsteadiness of flow, trabeculae height, and fluid scale ($Re$). Vortex formation occurred in intertrabecular regions for biologically relevant parameter spaces, wherein flow velocities increased. This indicates that trabecular morphology may alter intracardial flow patterns and hence ventricular shear stresses and morphogen gradients. A potential implication of this work is that the onset of vortical (disturbed) flows can upregulate Notch1 expression in endothelial cells in vitro and hence impacts chamber morphogensis, valvulogenesis, and the formation of the trabeculae themselves.

Multiple myeloma (MM), a plasma cell cancer, is associated with many health challenges, including damage to the kidney by tubulointerstitial fibrosis. We develop a mathematical model which captures the qualitative behavior of the cell and protein populations involved. Specifically, we model the interaction between cells in the proximal tubule of the kidney, free light chains, renal fibroblasts, and myeloma cells. We analyze the model for steady-state solutions to find a mathematically and biologically relevant stable steady-state solution. This foundational model provides a representation of dynamics between key populations in tubulointerstitial fibrosis that demonstrates how these populations interact to affect patient prognosis in patients with MM and renal impairment.

In animals, gas exchange between blood and tissues occurs in narrow vessels, whose diameter is comparable to that of a red blood cell. Red blood cells must deform to squeeze through these narrow vessels, transiently blocking or occluding the vessels they pass through. Although the dynamics of vessel occlusion have been studied extensively, it remains an open question why microvessels need to be so narrow. We study occlusive dynamics within a model microvascular network: the embryonic zebrafish trunk. We show that pressure feedbacks created when red blood cells enter the finest vessels of the trunk act together to uniformly partition red blood cells through the microvasculature. Using mathematical models as well as direct observation, we show that these occlusive feedbacks are tuned throughout the trunk network to prevent the vessels closest to the heart from short-circuiting the network. Thus occlusion is linked with another open question of microvascular function: how are red blood cells delivered at the same rate to each micro-vessel? Our analysis shows that tuning of occlusive feedbacks increase the total dissipation within the network by a factor of 11, showing that uniformity of flows rather than minimization of transport costs may be prioritized by the microvascular network.

Life and functioning of higher organisms depends on the continuous supply of metabolites to tissues and organs. What are the requirements on the transport network pervading a tissue to provide a uniform supply of nutrients, minerals, or hormones? To theoretically answer this question, we present an analytical scaling argument and numerical simulations on how flow dynamics and network architecture control active spread and uniform supply of metabolites by studying the example of xylem vessels in plants. We identify the fluid inflow rate as the key factor for uniform supply. While at low inflow rates metabolites are already exhausted close to flow inlets, too high inflow flushes metabolites through the network and deprives tissue close to inlets of supply. In between these two regimes, there exists an optimal inflow rate that yields a uniform supply of metabolites. We determine this optimal inflow analytically in quantitative agreement with numerical results. Optimizing network architecture by reducing the supply variance over all network tubes, we identify patterns of tube dilation or contraction that compensate sub-optimal supply for the case of too low or too high inflow rate.

A key process in the life of any multicellular organism is its development from a single egg into a full grown adult. The first step in this process often consists of forming a tissue layer out of randomly placed cells on the surface of the egg. We present a model for generating such a tissue, and find that the resulting cellular pattern corresponds to the Voronoi tessellation of the nuclei of the cells. Experimentally, we obtain the same result in both fruit flies and flour beetles, with a distribution of cell shapes that matches that of the model, without any adjustable parameters. Finally, we show that this pattern is broken when the cells do not all grow at the same rate.

We study dynamics emergent from a two-dimensional reaction--diffusion process modelled via a finite lattice dynamical system, as well as an analogous PDE system, involving spatially nonlocal interactions. These models govern the evolution of cells in a bioactive porous medium, with evolution of the local cell density depending on a coupled quasi--static fluid flow problem. We demonstrate differences emergent from the choice of a discrete lattice or a continuum for the spatial domain of such a process. We find long--time oscillations and steady states in cell density in both lattice and continuum models, but that the continuum model only exhibits solutions with vertical symmetry, independent of initial data, whereas the finite lattice admits asymmetric oscillations and steady states arising from symmetry-breaking bifurcations. We conjecture that it is the structure of the finite lattice which allows for more complicated asymmetric dynamics. Our analysis suggests that the origin of both types of oscillations is a nonlocal reaction-diffusion mechanism mediated by quasi-static fluid flow.

We developed simulation methodology to assess eventual therapeutic efficiency of exogenous multiparametric changes in a four-component cellular system described by the system of ordinary differential equations. The method is numerically implemented to simulate the temporal behavior of a cellular system of multiple myeloma cells. The problem is conceived as an inverse optimization task where the alternative temporal changes of selected parameters of the ordinary differential equations represent candidate solutions and the objective function quantifies the goals of the therapy. The system under study consists of two main cellular components, tumor cells and their cellular environment, respectively. The subset of model parameters closely related to the environment is substituted by exogenous time dependencies - therapeutic pulses combining continuous functions and discrete parameters subordinated thereafter to the optimization. Synergistic interaction of temporal parametric changes has been observed and quantified whereby two or more dynamic parameters show effects that absent if either parameter is stimulated alone. We expect that the theoretical insight into unstable tumor growth provided by the sensitivity and optimization studies could, eventually, help in designing combination therapies.

In order to understand the phenomenon of longevity in biological world, the relationship between the potential of longevity and the structural complexity of an organism is analyzed. I. The potential of longevity is the maximum lifespan of an organism if it lives in an ideal environment. The longevity of an organism includes two parts: the time for development (development time) and the time for structure-maintenance (maintenance time). II. The development time of an organism depends on its structural complexity. The maintenance time is related to two factors: the degree of damage-exposure and the potential of functionality for structure-maintenance. The potential of functionality is built in structural complexity. Since both of development time and maintenance time are determined by structural complexity, the potential of longevity of an organism is hidden in structural complexity. III. The individuals of different species have different potentials of longevity because they have different structural complexities. An animal has limited longevity because it has limited structural complexity. Limited structural complexity and limited longevity are essential for the survival of a species. IV. Despite having the same potential of longevity, the individuals of a species can have different lifespans. The lifespan of an individual is more related to the degree of damage-exposure, which is determined by the living environment and the living habit of the individual. In conclusion, the potential of longevity of an organism is hidden in structural complexity, but the real lifespan of an organism is more related to the living environment.

The production of mechanical stresses in living organisms largely relies on localized, force-generating active units embedded in filamentous matrices. Numerical simulations of discrete fiber networks with fixed boundaries have shown that buckling in the matrix dramatically amplifies the resulting active stresses. Here we extend this result to a bucklable continuum elastic medium subjected to an arbitrary external stress, and derive analytical expressions for the active, nonlinear constitutive relations characterizing the full active medium. Inserting these relations into popular "active gel" descriptions of living tissues and the cytoskeleton will enable investigations into nonlinear regimes previously inaccessible due to the phenomenological nature of these theories.

In this paper we introduce a new mathematical model for the active contraction of the cardiac muscle under different thermo-electric and nonlinear conductivity properties. The passive hyperelastic response of the tissue is described by an orthotropic exponential model, whereas the ionic activity dictates active contraction incorporated through the concept of orthotropic active strain. We use a fully incompressible formulation, and the generated strain modifies directly the conductivity mechanisms in the medium through the pull-back transformation. We also investigate the influence of thermo-electric effects in the onset of multiphysics emergent spatiotemporal dynamics, using nonlinear diffusion. It turns out that these ingredients have a key role in reproducing pathological rhythms such as ventricular fibrillation, during inflammatory events, for instance. The specific structure of the governing equations (written in terms of Kirchhoff stress, displacements, solid pressure, electric potential, activation generation, and ionic variables) suggests to cast the problem in mixed-primal form, and to use a mixed-primal finite element method for its numerical approximation. We finally explore the properties of the model, together with the importance of coupling variables, by means of a few computational experiments.

We numerically investigate the role of mechanical stress in modifying the conductivity properties of the cardiac tissue and its impact in computational models for cardiac electromechanics. We follow a theoretical framework recently proposed in [Cherubini, Filippi, Gizzi, Ruiz-Baier, JTB 2017], in the context of general reaction-diffusion-mechanics systems using multiphysics continuum mechanics and finite elasticity. In the present study, the adapted models are compared against preliminary experimental data of pig right ventricle fluorescence optical mapping. These data contribute to the characterization of the observed inhomogeneity and anisotropy properties that result from mechanical deformation. Our novel approach simultaneously incorporates two mechanisms for mechano-electric feedback (MEF): stretch-activated currents (SAC) and stress-assisted diffusion (SAD); and we also identify their influence into the nonlinear spatiotemporal dynamics. It is found that i) only specific combinations of the two MEF effects allow proper conduction velocity measurement; ii) expected heterogeneities and anisotropies are obtained via the novel stress-assisted diffusion mechanisms; iii) spiral wave meandering and drifting is highly mediated by the applied mechanical loading. We provide an analysis of the intrinsic structure of the nonlinear coupling using computational tests, conducted using a finite element method. In particular, we compare static and dynamic deformation regimes in the onset of cardiac arrhythmias and address other potential biomedical applications.

The behavior of hyperelastic energies depending on an internal parameter, which is a function of the deformation gradient, is discussed. As an example, the analysis of two models where the parameter describes the activation of a tetanized skeletal muscle tissue is presented. In those models, the activation parameter depends on the strain and it is shown the importance of considering the derivative of the parameter with respect to the strain in order to capture the proper stress-strain relations.

Background: The spread of tumors through neural routes is common in several types of cancer in which patients suffer from a moderate-to-severe neuropathy, neural damage and a distorted quality of life. Here we aim to examine the expression profiles of transient receptor potential vanilloid 1 (TRPV1) and of transient receptor potential vanilloid 4 (TRPV4), toll-like receptor 4 (TLR4) and extracellular signal-regulated kinase (ERK1/2), and to assess the possible therapeutic strategies through blockade of transient receptor potential (TRP) channels. Methods: Cancer was induced within the sciatic nerves of male Copenhagen rats, and tissues from dorsal root ganglia (DRG) were collected and used for measurements of immunofluorescence and Western blotting. The TRPV1 antagonist capsazepine, the selective TRPV4 antagonist HC-067047 and the calcium ions inhibitor ruthenium red were used to treat thermal and/or mechanical hyperalgesia. Results: Transient receptor potential vanilloid 1 showed a lower expression in DRGs on days 7 and 14. The expression of TRPV4, TLR4 and ERK1/2 showed an increase on day 3 then a decrease on days 7 and 14. TRPV1 and TLR4 as well as TRPV4 and ERK1/2 co-existed on the same neuronal cells. The neuropathic pain was reversed in dose-dependent manners by using the TRP antagonists and the calcium ions inhibitor. Conclusion: The decreased expression of TRPV1 and TRPV4 is associated with high activation. The increased expression of TLR4 and ERK1/2 reveals earlier immune response and tumor progression, respectively, and their ultimate decrease is an indicator of nerve damage. We studied the possible role of TRPV1 and TRPV4 in transducing cancer-induced hyperalgesia. The possible treatment strategies of cancer-induced thermal and/or mechanical hyperalgesia using capsazepine, HC-067047 and ruthenium red are examined.

We apply three optical coherence tomography (OCT) image analysis techniques to extract morphometric information from OCT images obtained on peripheral nerves of rat. The accuracy of each technique is evaluated against histological measurements accurate to +/-1 um. The three OCT techniques are: 1) average depth resolved profile (ADRP); 2) autoregressive spectral estimation (AR-SE); and, 3) correlation of the derivative spectral estimation (CoD-SE). We introduce a scanning window to the ADRP technique which provides transverse resolution, and improves epineurium thickness estimates - with the number of analysed images showing agreement with histology increasing from 2/10 to 5/10 (Kruskal-Wallis test, {\alpha} = 0.05). A new method of estimating epineurium thickness, using the AR-SE technique, showed agreement with histology in 6/10 analysed images (Kruskal-Wallis test, {\alpha} = 0.05). Using a tissue sample in which histology identified two fascicles with an estimated difference in mean fibre diameter of 4 um, the AR-SE and CoD-SE techniques both correctly identified the fascicle with larger fibre diameter distribution but incorrectly estimated the magnitude of this difference as 0.5um. The ability of OCT signal analysis techniques to extract accurate morphometric details from peripheral nerve is promising but restricted in depth by scattering in adipose and neural tissues.

An example of phase transition in natural complex systems is the qualitative and sudden change in the heart rhythm between sinus rhythm and atrial fibrillation (AF), the most common irregular heart rhythm in humans. While the system behavior is centrally controlled by the behavior of the sinoatrial node in sinus rhythm, the macro-scale collective behavior of the heart causes the micro-scale behavior in AF. To quantitatively analyze this causation shift associated with phase transition in human heart, we evaluated the causal architecture of the human cardiac system using the time series of multi-lead intracardiac unipolar electrograms in a series of spatiotemporal scales by generating a stochastic renormalization group. We found that the phase transition between sinus rhythm and AF is associated with a significant shift of the peak causation from macroscopic to microscopic scales. Causal architecture analysis may improve our understanding of causality in phase transitions in other natural and social complex systems.

Acute lymphoblastic leukemia (ALL) and chronic lymphocytic leukemia (CLL) are two major forms of leukemia that arise from lymphoid cells (LCs). ALL occurs mostly in children and CLL occurs mainly in old people. However, the Philadelphia-chromosome-positive ALL (Ph+-ALL) and the Ph-like ALL occur in both children and adults. To understand childhood leukemia/lymphoma, we have recently proposed two hypotheses on the causes and the mechanism of cell transformation of a LC. Hypothesis A is: repeated bone-remodeling during bone-growth and bone-repair may be a source of cell injuries of marrow cells including hematopoietic stem cells (HSCs), myeloid cells, and LCs. Hypothesis B is: a LC may have three pathways on transformation: a slow, a rapid, and an accelerated. We discuss in the present paper the developing mechanisms of ALL and CLL by these hypotheses. Having a peak incidence in young children, ALL may develop mainly as a result of rapid cell transformation of a lymphoblast (or pro-lymphocyte). Differently, Ph+-ALL and Ph-like ALL may develop as results of transformation of a lymphoblast via accelerated pathway. Occurring mainly in adults, CLL may be a result of transformation of a memory B-cell via slow pathway. By causing cell injuries of HSCs and LCs, repeated bone-remodeling during bone-growth and bone-repair may be related to the cell transformation of a LC. In conclusion, ALL may develop as a result of cell transformation of a lymphoblast via rapid or accelerated pathway; and repeated bone-remodeling during bone-growth may be a trigger for the cell transformation of a lymphoblast in a child.

Lymphomas are a large group of neoplasms developed from lymphoid cells (LCs) in lymph nodes (LNs) or lymphoid tissues (LTs). Some forms of lymphomas, including Burkitt lymphoma (BL), ALK+ anaplastic large cell lymphoma (ALK+-ALCL), and T-cell lymphoblastic lymphoma/leukemia (T-LBL), occur mainly in children and teenagers. Hodgkin's lymphoma (HL) has a peak incidence at age 20s. To understand pediatric lymphoma, we have recently proposed two hypotheses on the causes and the mechanism of cell transformation of a LC. Hypothesis A is: repeated bone-remodeling during bone-growth and bone-repair may be a source of cell injuries of marrow cells including hematopoietic stem cells (HSCs), myeloid cells, and LCs, and thymic involution may be a source of damage to the developing T-cells in thymus. Hypothesis B is: a LC may have three pathways on transformation: a slow, a rapid, and an accelerated. In this paper, we discuss pediatric lymphomas by this hypothesis. Having a peak incidence at young age, BL, T-LBL, ALK+-ALCL, and HL develop more likely as a result of rapid transformation of a LC. In BL, ALK+-ALCL, and HL, the cell transformations may be triggered by severe viral infections. In T-LBL, the cell transformation may be related to thymic involution. Occurring in both adults and children, diffuse large B-cell lymphoma (DLBCL) may develop via slow or accelerated pathway. In conclusion, pediatric lymphoma may develop as a result of "one-step" cell transformation of a LC, and severe viral infections may be the main trigger for the rapid transformation of a LC in a LN/LT.

Acute myeloid leukemia (AML) and chronic myeloid leukemia (CML) are two major formsof leukemia developed from myeloid cells (MCs). To understand why AML and CML occurin children, we analyzed the causes and the mechanism of cell transformation of a MC. I. Forthe MCs in marrow cavity, repeated bone-remodeling during bone-growth may be a source ofcell injuries. II. As a type of blood cell, a MC may have higher survivability from DNAchanges and require obtaining fewer cancerous properties for cell transformation than a tissuecell. III. Point DNA mutations (PDMs) and chromosome changes (CCs) are the two majortypes of DNA changes. CCs have three subtypes by effects on a cell: great effect CCs(GECCs), mild-effect CCs (MECCs), and intermediate-effect CCs (IECCs). A GECC affectsone or more genes and can alone trigger cell transformation. PDMs/MECCs are mostly mildand can accumulate in cells. Some of the PDMs/MECCs contribute to cell transformation. AnIECC affects one or more genes and participates in cell transformation. IV. Based on II andIII, we hypothesize that a MC may have two pathways on transformation: a slow and anaccelerated. Slow pathway is driven by accumulation of PDMs/MECCs. Accelerated pathwayis driven by accumulation of PDMs/MECCs/IECC(s). A transformation via slow pathwayoccurs at old age; whereas that via accelerated pathway occurs at any age. Thus, CML andpediatric AML may develop via accelerated pathway, and adult AML may develop via bothpathways. In conclusion, pediatric AML and CML may develop as a result of transformationof a MC via accelerated pathway; and repeated bone-remodeling for bone-growth may be atrigger for the transformation of a MC in a child.

A contemporary procedure to grow artificial tissue is to seed cells onto a porous biomaterial scaffold and culture it within a perfusion bioreactor to facilitate the transport of nutrients to growing cells. Typical models of cell growth for tissue engineering applications make use of spatially homogeneous or spatially continuous equations to model cell growth, flow of culture medium, nutrient transport, and their interactions. The network structure of the physical porous scaffold is often incorporated through parameters in these models, either phenomenologically or through techniques like mathematical homogenization. We derive a model on a square grid lattice to demonstrate the importance of explicitly modelling the network structure of the porous scaffold, and compare results from this model with those from a modified continuum model from the literature. We capture two-way coupling between cell growth and fluid flow by allowing cells to block pores, and by allowing the shear stress of the fluid to affect cell growth and death. We explore a range of parameters for both models, and demonstrate quantitative and qualitative differences between predictions from each of these approaches, including spatial pattern formation and local oscillations in cell density present only in the lattice model. These differences suggest that for some parameter regimes, corresponding to specific cell types and scaffold geometries, the lattice model gives qualitatively different model predictions than typical continuum models. Our results inform model selection for bioactive porous tissue scaffolds, aiding in the development of successful tissue engineering experiments and eventually clinically successful technologies.

Chronic wounds fail to proceed through an orderly and timely self healing process, resulting in cutaneous damage with full thickness in depth and leading to a major healthcare and economic burden worldwide. In the UK alone, 200,000 patients suffer from a chronic wound, whilst the global advanced wound care market is expected to reach nearly $11 million in 2022. Despite extensive research efforts so far, clinically-approved chronic wound therapies are still time-consuming, economically unaffordable and present restricted customisation. In this chapter, the role of collagen in the extracellular matrix of biological tissues and wound healing will be discussed, together with its use as building block for the manufacture of advanced wound dressings. Commercially-available collagen dressings and respective clinical performance will be presented, followed by an overview on the latest research advances in the context of multifunctional collagen systems for advanced wound care.

Pancreatic cancer has the poorest prognosis among all cancer types. Intraductal Papillary Mucinous Neoplasms (IPMNs) are radiographically identifiable precursors to pancreatic cancer; hence, early detection and precise risk assessment of IPMN are vital. In this work, we propose a Convolutional Neural Network (CNN) based computer aided diagnosis (CAD) system to perform IPMN diagnosis and risk assessment by utilizing multi-modal MRI. In our proposed approach, we use minimum and maximum intensity projections to ease the annotation variations among different slices and type of MRIs. Then, we present a CNN to obtain deep feature representation corresponding to each MRI modality (T1-weighted and T2-weighted). At the final step, we employ canonical correlation analysis (CCA) to perform a fusion operation at the feature level, leading to discriminative canonical correlation features. Extracted features are used for classification. Our results indicate significant improvements over other potential approaches to solve this important problem. The proposed approach doesn't require explicit sample balancing in cases of imbalance between positive and negative examples. To the best of our knowledge, our study is the first to automatically diagnose IPMN using multi-modal MRI.

Lymphoid leukemia (LL) and lymphoma are blood cancers developed from lymphoid cells (LCs). To understand the cause and the mechanism of cell transformation of a LC, we studied the potential sources of cell injuries of LCs and analyzed how DNA changes are generated and accumulate in LCs. I. The DNA changes that contribute to cell transformation of a LC can be generated in the LCs in marrow, thymus, lymph nodes (LNs), and/or lymphoid tissues (LTs). In LNs/LTs, pathogen-infections may be the main cause for cell injuries of LCs. In marrow cavity, repeated bone-remodeling during bone-growth and bone-repair, by producing toxic substances, may be a source of damage to hematopoietic cells, including hematopoietic stem cells (HSCs) and developing LCs. In thymus, thymic involution and death of stromal cells may be a damaging factor for the developing T-cells. II. Point DNA mutation (PDM) and chromosome change (CC) are the two major types of DNA changes. CCs include numerical CCs (NCCs) and structural CCs (SCCs). Generation of a PDM/SCC may be a result of Misrepair of DNA on DNA breaks. Generation of NCC is rather a consequence of dysfunction of cell division. III. Repeated cell injuries and cell proliferation drive the accumulation of DNA changes in LCs and HSCs. However, long-term accumulation of DNA changes occurs mainly in long-living stem cells including HSCs and memory cells. In conclusion, the DNA changes in LCs are generated and accumulate as a consequence of repeated cell injuries and repeated cell proliferation; and three potential sources of cell injuries of LCs may be: repeated bone-remodeling, long-term thymic involution, and repeated pathogen-infections.

We present a theory of pattern formation in growing domains inspired by biological examples of tissue development. Gradients of signaling molecules regulate growth, while growth changes these graded chemical patterns by dilution and advection. We identify a critical point of this feedback dynamics, which is characterized by spatially homogeneous growth and proportional scaling of patterns with tissue length. We apply this theory to the biological model system of the developing wing of the fruit fly \textit{Drosophila melanogaster} and quantitatively identify signatures of the critical point.

Biological growth is often driven by mechanical cues, such as changes in external pressure or tensile loading. Moreover, it is well known that many living tissues actively maintain a preferred level of mechanical internal stress, called the mechanical homeostasis. The tissue-level feedback mechanism by which changes of the local mechanical stresses affect growth is called a growth law within the theory of morphoelasticity, a theory for understanding the coupling between mechanics and geometry in growing and evolving biological materials. The goal of this article is to develop mathematical techniques to analyze growth laws and to explore issues of heterogeneity and growth stability. We discuss the growth dynamics of tubular structures, which are very common in biology (e.g. arteries, plant stems, airways) and model the homeostasis-driven growth dynamics of tubes which produces spatially inhomogeneous residual stress. We show that the stability of the homeostatic state depends nontrivially on the anisotropy of the growth response. The key role of anisotropy may provide a foundation for experimental testing of homeostasis-driven growth laws.

Immunotherapy plays a major role in tumour treatment, in comparison with other methods of dealing with cancer. The Kirschner-Panetta (KP) model of cancer immunotherapy describes the interaction between tumour cells, effector cells and interleukin-2 which are clinically utilized as medical treatment. The model selects a rich concept of immune-tumour dynamics. In this paper, approximate analytical solutions to KP model are represented by using the differential transform and Adomian decomposition. The complicated nonlinearity of the KP system causes the application of these two methods to require more involved calculations. The approximate analytical solutions to the model are compared with the results obtained by numerical fourth order Runge-Kutta method.