The role of the pericardium in the valveless, tubular heart of the tunicate Ciona savignyi

Mixed uncertainty analysis on pumping by peristaltic hearts using Dempster-Shafer theory

In this paper, we introduce the numerical strategy for mixed uncertainty propagation based on probability and Dempster-Shafer theories, and apply it to the computational model of peristalsis in a heart-pumping system. Specifically, the stochastic uncertainty in the system is represented with random variables while epistemic uncertainty is represented using non-probabilistic uncertain variables with belief functions. The mixed uncertainty is propagated through the system, resulting in the uncertainty in the chosen quantities of interest (QoI, such as flow volume, cost of transport and work). With the introduced numerical method, the uncertainty in the statistics of QoIs will be represented using belief functions. With three representative probability distributions consistent with the belief structure, global sensitivity analysis has also been implemented to identify important uncertain factors and the results have been compared between different peristalsis models. To reduce the computational cost, physics constrained generalized polynomial chaos method is adopted to construct cheaper surrogates as approximations for the full simulation.

Pulsation waves along the Ciona heart tube reverse by bimodal rhythms expressed by a remote pair of pacemakers

The heart of ascidians (marine invertebrate chordates) has a tubular structure, and heartbeats propagate from one end to the other. The direction of pulsation waves intermittently reverses in the heart of ascidians and their relatives; however, the underlying mechanisms remain unclear. We herein performed a series of experiments to characterize the pacemaker systems in isolated hearts and their fragments, and applied a mathematical model to examine the conditions leading to heart reversals. The isolated heart of Ciona robusta autonomously generated pulsation waves at ∼20 to 25 beats min-1 with reversals at ∼1 to 10 min intervals. Experimental bisections of isolated hearts revealed that independent pacemakers resided on each side and also that their beating frequencies periodically changed as they expressed bimodal rhythms, which comprised an ∼1.25 to 5.5 min acceleration/deceleration cycle of a beating rate of between 0 and 25 beats min-1. Only fragments including 5% or shorter terminal regions of the heart tube maintained autonomous pulsation rhythms, whereas other regions did not. Our mathematical model, based on FitzHugh-Nagumo equations applied to a one-dimensional alignment of cells, demonstrated that the difference between frequencies expressed by the two independent terminal pacemakers determined the direction of propagated waves. Changes in the statuses of terminal pacemakers between the excitatory and oscillatory modes as well as in their endogenous oscillation frequencies were sufficient to lead to heart reversals. These results suggest that the directions of pulsation waves in the Ciona heart reverse according to the changing rhythms independently expressed by remotely coupled terminal pacemakers.

Central nervous system regeneration in ascidians: cell migration and differentiation

Adult ascidians have the capacity to regenerate the central nervous system (CNS) and are therefore excellent models for studies on neuroregeneration. The possibility that undifferentiated blood cells are involved in adult neuroregeneration merits investigation. We analyzed the migration, circulation, and role of hemocytes of the ascidian Styela plicata in neuroregeneration. Hemocytes were removed and incubated with superparamagnetic iron oxide nanoparticles (SPION), and these SPION-labeled hemocytes were injected back into the animals (autologous transplant), followed by neurodegeneration with the neurotoxin 3-acetylpyridine (3AP). Magnetic resonance imaging showed that 1, 5, and 10 days after injury, hemocytes migrated to the intestinal region, siphons, and CNS. Immunohistochemistry revealed that the hemocytes that migrated to the CNS were putative stem cells (P-element-induced wimpy testis + or PIWI + cells). In the cortex of the neural ganglion, migrated hemocytes started to lose their PIWI labeling 5 days after injury, and 10 days later started to show β-III tubulin labeling. In the neural gland, however, the hemocytes remained undifferentiated during the entire experimental period. Transmission electron microscopy revealed regions in the neural gland with characteristics of neurogenic niches, not previously reported in ascidians. These results showed that migration of hemocytes to the hematopoietic tissue and to the 3AP-neurodegenerated region is central to the complex mechanism of neuroregeneration.

Modeling action potential reversals in tunicate hearts

Tunicates are small invertebrates which possess a unique ability to reverse flow in their hearts. Scientists have debated various theories regarding how and why flow reversals occur. Here we explore the electrophysiological basis for reversals by simulating action potential propagation in an idealized model of the tubelike tunicate heart. Using asymptotic formulas for action potential duration and conduction velocity, we propose tunicate-specific parameters for a two-current ionic model of the action potential. Then, using a kinematic model, we derive analytical criteria for reversals to occur. These criteria inform subsequent numerical simulations of action potential propagation in a fiber paced at both ends. In particular, we explore the role that variability of pacemaker firing rates plays in generating reversals, and we identify various favorable conditions for triggering retrograde propagation. Our analytical framework extends to other species; for instance, it can be used to model competition between the sinoatrial node and abnormal ectopic foci in human heart tissue.

Transcription Factors of the bHLH Family Delineate Vertebrate Landmarks in the Nervous System of a Simple Chordate

Tunicates are marine invertebrates whose tadpole-like larvae feature a highly simplified version of the chordate body plan. Similar to their distant vertebrate relatives, tunicate larvae develop a regionalized central nervous system and form distinct neural structures, which include a rostral sensory vesicle, a motor ganglion, and a caudal nerve cord. The sensory vesicle contains a photoreceptive complex and a statocyst, and based on the comparable expression patterns of evolutionarily conserved marker genes, it is believed to include proto-hypothalamic and proto-retinal territories. The evolutionarily conserved molecular fingerprints of these landmarks of the vertebrate brain consist of genes encoding for different transcription factors, and of the gene batteries that they control, and include several members of the bHLH family. Here we review the complement of bHLH genes present in the streamlined genome of the tunicate Ciona robusta and their current classification, and summarize recent studies on proneural bHLH transcription factors and their expression territories. We discuss the possible roles of bHLH genes in establishing the molecular compartmentalization of the enticing nervous system of this unassuming chordate.

Uncertainty quantification reveals the physical constraints on pumping by peristaltic hearts

Most biological functional systems are complex, and this complexity is a fundamental driver of diversity. Because input parameters interact in complex ways, a holistic understanding of functional systems is key to understanding how natural selection produces diversity. We present uncertainty quantification (UQ) as a quantitative analysis tool on computational models to study the interplay of complex systems and diversity. We investigate peristaltic pumping in a racetrack circulatory system using a computational model and analyse the impact of three input parameters (Womersley number, compression frequency, compression ratio) on flow and the energetic costs of circulation. We employed two models of peristalsis (one that allows elastic interactions between the heart tube and fluid and one that does not), to investigate the role of elastic interactions on model output. A computationally cheaper surrogate of the input parameter space was created with generalized polynomial chaos expansion to save computational resources. Sobol indices were then calculated based on the generalized polynomial chaos expansion and model output. We found that all flow metrics were highly sensitive to changes in compression ratio and insensitive to Womersley number and compression frequency, consistent across models of peristalsis. Elastic interactions changed the patterns of parameter sensitivity for energetic costs between the two models, revealing that elastic interactions are probably a key physical metric of peristalsis. The UQ analysis created two hypotheses regarding diversity: favouring high flow rates (where compression ratio is large and highly conserved) and minimizing energetic costs (which avoids combinations of high compression ratios, high frequencies and low Womersley numbers).

Surprising simplicities and syntheses in limbless self-propulsion in sand

Animals moving on and in fluids and solids move their bodies in diverse ways to generate propulsion and lift forces. In fluids, animals can wiggle, stroke, paddle or slap, whereas on hard frictional terrain, animals largely engage their appendages with the substrate to avoid slip. Granular substrates, such as desert sand, can display complex responses to animal interactions. This complexity has led to locomotor strategies that make use of fluid-like or solid-like features of this substrate, or combinations of the two. Here, we use examples from our work to demonstrate the diverse array of methods used and insights gained in the study of both surface and subsurface limbless locomotion in these habitats. Counterintuitively, these seemingly complex granular environments offer certain experimental, theoretical, robotic and computational advantages for studying terrestrial movement, with the potential for providing broad insights into morphology and locomotor control in fluids and solids, including neuromechanical control templates and morphological and behavioral evolution. In particular, granular media provide an excellent testbed for a locomotion framework called geometric mechanics, which was introduced by particle physicists and control engineers in the last century, and which allows quantitative analysis of alternative locomotor patterns and morphology to test for control templates, optimality and evolutionary alternatives. Thus, we posit that insights gained from movement in granular environments can be translated into principles that have broader applications across taxa, habitats and movement patterns, including those at microscopic scales.

The ventral peptidergic system of the adult ascidian Ciona robusta (Ciona intestinalis Type A) insights from a transgenic animal model

Ascidians are the sister group of vertebrates and occupy a critical position in explorations of the evolution of the endocrine and nervous systems of chordates. Here, we describe the complete ventral peptidergic system in adult transgenic Ciona robusta (Ciona intestinalis Type A) which expresses the Kaede reporter gene driven by the prohormone convertase 2 (PC2) gene promoter. Numerous PC2 promoter-driven fluorescent (Kaede-positive) non-neural cells were distributed in the blood sinus located at the anterior end of the pharynx, suggesting the acquisition of a peptidergic circulatory system in Ciona. Kaede-positive ciliated columnar cells, rounded cells, and tall ciliated cells were observed in the alimentary organs, including the endostyle, pharynx, esophagus, stomach, and intestine, suggesting that digestive functions are regulated by multiple peptidergic systems. In the heart, Kaede-positive neurons were located in the ring-shaped plexus at both ends of the myocardium. Nerve fiber-like tracts ran along the raphe and appeared to be connected with the plexuses. Such unique structures suggest a role for the peptidergic system in cardiac function. Collectively, the present anatomic analysis revealed the major framework of the ventral peptidergic system of adult Ciona, which could facilitate investigations of peptidergic regulation of the pharynx, endostyle, alimentary tissues, and heart.

Pharmacological analysis of the transmembrane action potential configuration in myoepithelial cells of the spontaneously beating heart of the ascidian Styela rustica in vitro

The mechanisms of action potential (AP) generation in the myoepithelial cells of the tunicate heart are not yet well understood. Here, an attempt was made to elucidate these mechanisms by analyzing the effects of specific blockers of K⁺, Na⁺ and Ca²⁺ currents on the configuration of transmembrane APs and their frequency in the spontaneously beating ascidian heart. In addition, an immunocytochemical analysis of heart myoepithelial cells was performed. Staining with anti-FMRF-amide and anti-tubulin antibodies did not reveal any nerve elements within the heart tube. Treatment with 1 mmol l^-1 TEA (I_K blocker) resulted in depolarization of heart cell sarcolemma by 10 mV, and inhibition of APs generation was recorded after 3 min of exposure. Prior to this moment, the frequency of AP generation in a burst decreased from 16-18 to 2 beats min^-1 owing to prolongation of the diastole. After application of ivabradine (3 or 10 µmol l^-1), the spontaneous APs generation frequency decreased by 24%. Based on these results and published data, it is concluded that the key role in the automaticity of the ascidian heart is played by the outward K⁺ currents, Na⁺ currents, activated hyperpolarization current I_f and a current of unknown nature I_X.

IB2d: a Python and MATLAB implementation of the immersed boundary method

The development of fluid-structure interaction (FSI) software involves trade-offs between ease of use, generality, performance, and cost. Typically there are large learning curves when using low-level software to model the interaction of an elastic structure immersed in a uniform density fluid. Many existing codes are not publicly available, and the commercial software that exists usually requires expensive licenses and may not be as robust or allow the necessary flexibility that in house codes can provide. We present an open source immersed boundary software package, IB2d, with full implementations in both MATLAB and Python, that is capable of running a vast range of biomechanics models and is accessible to scientists who have experience in high-level programming environments. IB2d contains multiple options for constructing material properties of the fiber structure, as well as the advection-diffusion of a chemical gradient, muscle mechanics models, and artificial forcing to drive boundaries with a preferred motion.