Thomas D. Montenegro-Johnson

Physics of rheologically enhanced propulsion: Different strokes in generalized Stokes

Shear-thinning is an important rheological property of many biological fluids, such as mucus, whereby the apparent viscosity of the fluid decreases with shear. Certain microscopic swimmers have been shown to progress more rapidly through shear-thinning fluids, but is this behavior generic to all microscopic swimmers, and what are the physics through which shear-thinning rheology affects a swimmers propulsion? We examine swimmers employing prescribed stroke kinematics in two-dimensional, inertialess Carreau fluid: shear-thinning “generalized Stokes” flow. Swimmers are modeled, using the method of femlets, by a set of immersed, regularized forces. The equations governing the fluid dynamics are then discretized over a body-fitted mesh and solved with the finite element method. We analyze the locomotion of three distinct classes of microswimmer: (1) conceptual swimmers comprising sliding spheres employing both one- and two-dimensional strokes, (2) slip-velocity envelope models of ciliates commonly referred t...

Left-Right Organizer Flow Dynamics: How Much Cilia Activity Reliably Yields Laterality?

Internal organs are asymmetrically positioned inside the body. Embryonic motile cilia play an essential role in this process by generating a directional fluid flow inside the vertebrate left-right organizer. Detailed characterization of how fluid flow dynamics modulates laterality is lacking. We used zebrafish genetics to experimentally generate a range of flow dynamics. By following the development of each embryo, we show that fluid flow in the left-right organizer is asymmetric and provides a good predictor of organ laterality. This was tested in mosaic organizers composed of motile and immotile cilia generated by dnah7 knockdowns. In parallel, we used simulations of fluid dynamics to analyze our experimental data. These revealed that fluid flow generated by 30 or more cilia predicts 90% situs solitus, similar to experimental observations. We conclude that cilia number, dorsal anterior motile cilia clustering, and left flow are critical to situs solitus via robust asymmetric charon expression.

Modelling the fluid mechanics of cilia and flagella in reproduction and development

AbstractCilia and flagella are actively bending slender organelles, performing functions such as motility, feeding and embryonic symmetry breaking. We review the mechanics of viscous-dominated microscale flow, including time-reversal symmetry, drag anisotropy of slender bodies, and wall effects. We focus on the fundamental force singularity, higher-order multipoles, and the method of images, providing physical insight and forming a basis for computational approaches. Two biological problems are then considered in more detail: 1) left-right symmetry breaking flow in the node, a microscopic structure in developing vertebrate embryos, and 2) motility of microswimmers through non-Newtonian fluids. Our model of the embryonic node reveals how particle transport associated with morphogenesis is modulated by the gradual emergence of cilium posterior tilt. Our model of swimming makes use of force distributions within a body-conforming finite-element framework, allowing the solution of nonlinear inertialess Carreau flow. We find that a three-sphere model swimmer and a model sperm are similarly affected by shear-thinning; in both cases swimming due to a prescribed beat is enhanced by shear-thinning, with optimal Deborah number around 0.8. The sperm exhibits an almost perfect linear relationship between velocity and the logarithm of the ratio of zero to infinite shear viscosity, with shear-thickening hindering cell progress.

Digital Single-Cell Analysis of Plant Organ Development Using 3DCellAtlas

A tool was developed to extract biologically relevant information from quantitative 3D image data, allowing the relationship between diverse regulatory networks and 3D organ growth to be revealed. Diverse molecular networks underlying plant growth and development are rapidly being uncovered. Integrating these data into the spatial and temporal context of dynamic organ growth remains a technical challenge. We developed 3DCellAtlas, an integrative computational pipeline that semiautomatically identifies cell types and quantifies both 3D cellular anisotropy and reporter abundance at single-cell resolution across whole plant organs. Cell identification is no less than 97.8% accurate and does not require transgenic lineage markers or reference atlases. Cell positions within organs are defined using an internal indexing system generating cellular level organ atlases where data from multiple samples can be integrated. Using this approach, we quantified the organ-wide cell-type-specific 3D cellular anisotropy driving Arabidopsis thaliana hypocotyl elongation. The impact ethylene has on hypocotyl 3D cell anisotropy identified the preferential growth of endodermis in response to this hormone. The spatiotemporal dynamics of the endogenous DELLA protein RGA, expansin gene EXPA3, and cell expansion was quantified within distinct cell types of Arabidopsis roots. A significant regulatory relationship between RGA, EXPA3, and growth was present in the epidermis and endodermis. The use of single-cell analyses of plant development enables the dynamics of diverse regulatory networks to be integrated with 3D organ growth.

Organized chaos in Kupffer's vesicle: how a heterogeneous structure achieves consistent left-right patterning.

Successful establishment of left-right asymmetry is crucial to healthy vertebrate development. In many species this process is initiated in a ciliated, enclosed cavity, for example Kupffers vesicle (KV) in zebrafish. The microarchitecture of KV is more complex than that present in the left-right organizer of many other species. While swirling flow in KV is recognized as essential for left-right patterning, its generation, nature and conversion to asymmetric gene expression are only beginning to be fully understood. We recently [Sampaio, P et al. Dev Cell 29:716–728] combined imaging, genetics and fluid dynamics simulation to characterize normal and perturbed ciliary activity, and their correlation to asymmetric charon expression and embryonic organ fate. Randomness in cilia number and length have major implications for robust flow generation; even a modest change in mean cilia length has a major effect on flow speed to due to nonlinear scaling arising from fluid mechanics. Wildtype, and mutant embryos with normal liver laterality, exhibit stronger flow on the left prior to asymmetric inhibition of charon. Our discovery of immotile cilia, taken with data on morphant embryos with very few cilia, further support the role of mechanosensing in initiating and/or enhancing flow conversion into gene expression.

Geometric pumping in autophoretic channels

Many microfluidic devices use macroscopic pressure differentials to overcome viscous friction and generate flows in microchannels. In this work, we investigate how the chemical and geometric properties of the channel walls can drive a net flow by exploiting the autophoretic slip flows induced along active walls by local concentration gradients of a solute species. We show that chemical patterning of the wall is not required to generate and control a net flux within the channel, rather channel geometry alone is sufficient. Using numerical simulations, we determine how geometric characteristics of the wall influence channel flow rate, and confirm our results analytically in the asymptotic limit of lubrication theory.

Spermatozoa scattering by a microchannel feature: an elastohydrodynamic model

Sperm traverse their microenvironment through viscous fluid by propagating flagellar waves; the waveform emerges as a consequence of elastic structure, internal active moments and low Reynolds number fluid dynamics. Engineered microchannels have recently been proposed as a method of sorting and manipulating motile cells; the interaction of cells with these artificial environments therefore warrants investigation. A numerical method is presented for large-amplitude elastohydrodynamic interaction of active swimmers with domain features. This method is employed to examine hydrodynamic scattering by a model microchannel backstep feature. Scattering is shown to depend on backstep height and the relative strength of viscous and elastic forces in the flagellum. In a ‘high viscosity’ parameter regime corresponding to human sperm in cervical mucus analogue, this hydrodynamic contribution to scattering is comparable in magnitude to recent data on contact effects, being of the order of 5°–10°. Scattering can be positive or negative depending on the relative strength of viscous and elastic effects, emphasizing the importance of viscosity on the interaction of sperm with their microenvironment. The modulation of scattering angle by viscosity is associated with variations in flagellar asymmetry induced by the elastohydrodynamic interaction with the boundary feature.

A regularised singularity approach to phoretic problems

Abstract.An efficient, accurate, and flexible numerical method is proposed for the solution of the swimming problem of one or more autophoretic particles in the purely diffusive limit. The method relies on successive boundary element solutions of the Laplacian and the Stokes flow equations using regularised Green’s functions for swift, simple implementations, an extension of the well-known method of “regularised stokeslets” for Stokes flow problems. The boundary element method is particularly suitable for phoretic problems, since no quantities in the domain bulk are required to compute the swimming velocity. For time-dependent problems, the method requires no re-meshing and simple boundaries such as a plane wall may be added at no increase to the size of the linear system through the method of images. The method is validated against two classical examples for which an analytical or semi-analytical solution is known, a two-sphere system and a Janus particle, and provides a rigorous computational pipeline to address further problems with complex geometry and multiple bodies.Graphical abstract

Dynamics of cilia length in left–right development

Reduction in the length of motile cilia in the zebrafish left–right organizer (LRO), also known as Kupffers vesicle, has a large impact on left–right development. Here we demonstrate through genetic overexpression in zebrafish embryos and mathematical modelling that the impact of increased motile cilia length in embryonic LRO fluid flow is milder than that of short cilia. Through Arl13b overexpression, which increases cilia length without impacting cilia beat frequency, we show that the increase in cilium length is associated with a decrease in beat amplitude, resulting in similar flow strengths for Arl13b overexpression and wild-type (WT) embryos, which were not predicted by current theory. Longer cilia exhibit pronounced helical beat patterns and, consequently, lower beat amplitudes relative to WT, a result of an elastohydrodynamic shape transition. For long helical cilia, fluid dynamics modelling predicts a mild (approx. 12%) reduction in the torque exerted on the fluid relative to the WT, resulting in a proportional reduction in flow generation. This mild reduction is corroborated by experiments, providing a mechanism for the mild impact on organ situs.

Flow analysis of the low Reynolds number swimmer C. elegans

Swimming cells and microorganisms are a critical component of many biological processes. In order to better interpret experimental studies of low Reynolds number swimming, we combine experimental and numerical methods to perform an analysis of the flow-field around the swimming nematode Caenorhabditis elegans. We first use image processing and particle tracking velocimetry to extract the body shape, kinematics, and flow-fields around the nematode. We then construct a three-dimensional model using the experimental swimming kinematics and employ a boundary element method to simulate flow-fields, obtaining very good quantitative agreement with experiment. We use this numerical model to show that calculation of flow shear rates using purely planar data results in significant underestimates of the true three-dimensional value. Applying symmetry arguments, validated against numerics, we demonstrate that the out-of-plane contribution can be accounted for via incompressibility and therefore simply calculated from particle tracking velocimetry. Our results show how fundamental fluid mechanics considerations may be used to improve the accuracy of measurements in biofluiddynamics.