Fluid mechanics of nodal flow due to embryonic primary cilia

Breaking of left–right symmetry is crucial in vertebrate development. The role of cilia-driven flow has been the subject of many recent publications, but the underlying mechanisms remain controversial. At approximately 8 days post-fertilization, after the establishment of the dorsal–ventral and anterior–posterior axes, a depressed structure is found on the ventral side of mouse embryos, termed the ventral node. Within the node, ‘whirling’ primary cilia, tilted towards the posterior, drive a flow implicated in the initial left–right signalling asymmetry. However, the underlying fluid mechanics have not been fully and correctly explained until recently and accurate characterization is required in determining how the flow triggers the downstream signalling cascades. Using the approximation of resistive force theory, we show how the flow is produced and calculate the optimal configuration to cause maximum flow, showing excellent agreement with in vitro measurements and numerical simulation, and paralleling recent analogue experiments. By calculating numerical solutions of the slender body theory equations, we present time-dependent physically based fluid dynamics simulations of particle pathlines in flows generated by large arrays of beating cilia, showing the far-field radial streamlines predicted by the theory

Mathematical modelling of cilia driven transport of biological fluids

Cilia-driven flow occurs in the airway surface liquid, in the female and male reproductive tracts and enables symmetry-breaking in the embryonic node. Viscoelastic rheology is found in healthy states in some systems, whereas in others may characterise disease, motivating the development of mathematical models that take this effect into account. We derive the fundamental solution for linear viscoelastic flow, which is subsequently used as a basis for slender-body theory. Our numerical algorithm allows efficient computation of three-dimensional time-dependent flow, bending moments, power and particle transport. We apply the model to the large-amplitude motion of a single cilium in a linear Maxwell liquid. A relatively short relaxation time of just 0.032 times the beat period significantly reduces forces, bending moments, power and particle transport, the last variable exhibiting exponential decay with relaxation time. A test particle is propelled approximately one-fifth as quickly along the direction of cilia beating for scaled relaxation time 0.032 as in the Newtonian case, and mean volume flow is abolished, emphasizing the sensitivity of cilia function to fluid rheology. These results may have implications for flow in the airways, where the transition from Newtonian to viscoelastic rheology in the peri-ciliary fluid may reduce clearance

A viscoelastic traction layer model of mucociliary flow

A new mathematical model of the transport of mucus and periciliary liquid (PCL) in the airways by cilia is presented. Mucus is represented by a linearly viscoelastic fluid, the mat of cilia is modelled as an ‘active porous medium.’ The propulsive effect of the cilia is modelled by a time-dependent force acting in a shear-thinned ‘traction layer’ between the mucus and the PCL. The effects of surface and interface tension are modelled by constraining the mucus free surface and mucus–PCL interface to be flat. It is assumed that the epithelium is impermeable to fluid. Using Fourier series, the system is converted into ODEs and solved numerically. We calculate values for mean mucus speed close to those observed by Matsui et~al. [{J. Clin. Invest.}, 102(6):1125’1131, 1998], (~40 μms−1). We obtain more detail regarding the dynamics of the flow and the nonlinear relationships between physical parameters in healthy and diseased states than in previously published models. Pressure gradients in the PCL caused by interface and surface tension are vital to ensuring efficient transport of mucus, and the role of the mucus–PCL interface appears to be to support such pressure gradients, ensuring efficient transport. Mean transport of PCL is found to be very small, consistent with previous analyses, providing insight into theories regarding the normal tonicity of PCL

The influence of receptor-mediated interactions on reaction-diffusion mechanisms of cellular self-organisation

Understanding the mechanisms governing and regulating self-organisation in the developing embryo is a key challenge that has puzzled and fascinated scientists for decades. Since its conception in 1952 the Turing model has been a paradigm for pattern formation, motivating numerous theoretical and experimental studies, though its verification at the molecular level in biological systems has remained elusive. In this work, we consider the influence of receptor-mediated dynamics within the framework of Turing models, showing how non-diffusing species impact the conditions for the emergence of self-organisation. We illustrate our results within the framework of hair follicle pre-patterning, showing how receptor interaction structures can be constrained by the requirement for patterning, without the need for detailed knowledge of the network dynamics. Finally, in the light of our results, we discuss the ability of such systems to pattern outside the classical limits of the Turing model, and the inherent dangers involved in model reduction

Mammalian Sperm Motility: Observation and Theory

Mammalian spermatozoa motility is a subject of growing importance because of rising human infertility and the possibility of improving animal breeding. We highlight opportunities for fluid and continuum dynamics to provide novel insights concerning the mechanics of these specialized cells, especially during their remarkable journey to the egg. The biological structure of the motile sperm appendage, the flagellum, is described and placed in the context of the mechanics underlying the migration of mammalian sperm through the numerous environments of the female reproductive tract. This process demands certain specific changes to flagellar movement and motility for which further mechanical insight would be valuable, although this requires improved modeling capabilities, particularly to increase our understanding of sperm progression in vivo. We summarize current theoretical studies, highlighting the synergistic combination of imaging and theory in exploring sperm motility, and discuss the challenges for future observational and theoretical studies in understanding the underlying mechanics.\ud Acronyms and Definitions\ud Acrosome: the cap of the sperm head containing enzymes allowing penetration of the zona pellucida via the acrosome reaction\ud Adenosine triphosphate (ATP): the currency unit of chemical energy transfer in living cells\ud Axoneme: a phylogenetically conserved structure within the eukaryotic flagellum consisting of a ring of nine microtubule doublets and a central pair, frequently referred to as 9 + 2\ud Bending moment density: the moment per unit length associated with flagellar bending; it can be divided into a hydrodynamic moment, an elastic moment (from the flagellar bending stiffness), an active moment (generated by dyneins exerting forces between adjacent microtubule doublets), and a passive moment resisting shear\ud Capacitation: the physiological state of a sperm required for fertilization, which is accompanied by the motility patterns associated with hyperactivation, characterized in saline by high-amplitude asymmetric beating\ud Central pair: a pair of microtubules along the length of the axoneme, symmetrically and slightly offset from the axoneme centerline\ud Cumulus oophorus: the outer vestment of the mammalian egg consisting of hundreds of cells radiating out from the egg embedded within a non-Newtonian hyaluronic acid gel\ud Dynein: a molecular motor within the axoneme, attached between adjacent microtubule doublets, that exerts a shearing force to induce axonemal bending\ud Flagellum: a motile cellular appendage that drives the swimming of sperm and other cells; this article focuses on the eukaryotic flagellum\ud Microtubule doublet: a pair of proteinaceous filament structures running the length of the axoneme; dyneins drive their bending, which induces flagellar motion\ud Mid-piece: the region of a sperm flagellum with a mitochondrial sheath, where ATP is generated\ud Oocyte: the egg\ud Outer dense fibers and fibrous sheath: accessory structures reinforcing the mammalian sperm flagellum; the combined axoneme and accessory structures are referred to as 9+9+2\ud Resistive-force theory: an approximation for the local drag of a slender filament element in Stokes flow (or a viscoelastic generalization thereof)\ud Rheotaxis: directed motility in response to the influence of fluid flow\ud Shear: in the context of the flagellum, the relative movement of adjacent microtubule doublets\ud Slender-body theory: an improved approximation for the local drag on a slender filament element in Stokes flow (or a viscoelastic generalization thereof)\ud Zona pellucida: a tough glycoprotein coat between the human egg and the cumulus oophorus, which a sperm must penetrate for successful fertilizatio

Incorporating spatial correlations into multispecies mean-field models

In biology, we frequently observe different species existing within the same environment. For example, there are many cell types in a tumour, or different animal species may occupy a given habitat. In modeling interactions between such species, we often make use of the mean-field approximation, whereby spatial correlations between the locations of individuals are neglected. Whilst this approximation holds in certain situations, this is not always the case, and care must be taken to ensure the mean-field approximation is only used in appropriate settings. In circumstances where the mean-field approximation is unsuitable, we need to include information on the spatial distributions of individuals, which is not a simple task. In this paper, we provide a method that overcomes many of the failures of the mean-field approximation for an on-lattice volume-excluding birth-death-movement process with multiple species. We explicitly take into account spatial information on the distribution of individuals by including partial differential equation descriptions of lattice site occupancy correlations. We demonstrate how to derive these equations for the multispecies case and show results specific to a two-species problem. We compare averaged discrete results to both the mean-field approximation and our improved method, which incorporates spatial correlations. We note that the mean-field approximation fails dramatically in some cases, predicting very different behavior from that seen upon averaging multiple realizations of the discrete system. In contrast, our improved method provides excellent agreement with the averaged discrete behavior in all cases, thus providing a more reliable modeling framework. Furthermore, our method is tractable as the resulting partial differential equations can be solved efficiently using standard numerical techniques

Functional Imaging of the Outer Retinal Complex using High Fidelity Imaging Retinal Densitometry

We describe a new technique, high fidelity Imaging Retinal Densitometry (IRD), which probes the functional integrity of the outer retinal complex. We demonstrate the ability of the technique to map visual pigment optical density and synthesis rates in eyes with and without macular disease. A multispectral retinal imaging device obtained precise measurements of retinal reflectance over space and time. Data obtained from healthy controls and 5 patients with intermediate AMD, before and after photopigment bleaching, were used to quantify visual pigment metrics. Heat maps were plotted to summarise the topography of rod and cone pigment kinetics and descriptive statistics conducted to highlight differences between those with and without AMD. Rod and cone visual pigment synthesis rates in those with AMD (v = 0.043 SD 0.019 min-1 and v = 0.119 SD 0.046 min-1, respectively) were approximately half those observed in healthy controls (v = 0.079 SD 0.024 min-1 for rods and v = 0.206 SD 0.069 min-1 for cones). By mapping visual pigment kinetics across the central retina, high fidelity IRD provides a unique insight into outer retinal complex function. This new technique will improve the phenotypic characterisation, diagnosis and treatment monitoring of various ocular pathologies, including AMD

Atrazine and Alachlor Adsorption Characteristics to Benchmark Soil Series in Eastern South Dakota

Corn, grain sorghum, and soybean are grown on about six million acres in eastern South Dakota each year. Two herbicides used routinely for weed control are atrazine(6-chloro-N-ethyl-N’-(1-methylethyl)-1,3,5-triazine-2,4-diamine) in corn and grain sorghum and alachlor (2-chloro-N-(2,6-diethylphenyl)-N-(methoxymethyl)acetamide) in all three crops. Six benchmark soil series that include a majority of the cropped acres treated with these herbicides are the Egan, Moody, Nora, and Brandt silty clay loams, and Clarno and Enet loams. Batch adsorption studies determined atrazine and alachlor binding characteristics to these soils and aids in assessing the amount of herbicide available for movement. These data also provide a basis for future use and management decisions for these and other related herbicides in similar soils. Soils from three horizons (A, B, and C) for each soil type were treated with atrazine or alachlor at four herbicide concentrations. Atrazine and alachlor sorption partition coefficients differed most in A horizon soils and ranged from 2.16 to 5.35 µmol1-1/n L1/nKg-1 for atrazine and 1.95 to 5.78 µmol1-1/n L1/nKg-1 for alachlor. Atrazine binding to A horizon soils ranked as Brandt \u3eEgan = Moody \u3e Enet = Clarno \u3e Nora. Alachlor binding to A horizon soils ranked as Brandt \u3eMoody \u3e Nora \u3e Enet \u3e Clarno. B and C horizon soils had lower binding for both herbicides; the sorption partition coefficient for atrazine ranged from 0.12 to 1.9 µmol1-1/n L1/nKg-1 while alachlor ranged from 0.43 to 1.64 µmol1-1/n L1/nKg-1. These data indicate that some soil types would be more susceptible to herbicide leaching than others. Once the herbicide moves through the A horizon, it may move rapidly through the lower soil profile (because of the decrease in binding capacity), and therefore, increase the vulnerability of the aquifer to contamination. Best management practices for these herbicides are being investigated to limit their movement through soil