Hydrodynamics of diatom chains and semiflexible fibres

Diatoms are non-motile, unicellular phytoplankton that have the ability to form colonies in the form of chains. Depending upon the species of diatoms and the linking structures that hold the cells together, these chains can be quite stiff or very flexible. Recently, the bending rigidities of some species of diatom chains have been quantified. In an effort to understand the role of flexibility in nutrient uptake and aggregate formation, we begin by developing a three-dimensional model of the coupled elastic–hydrodynamic system of a diatom chain moving in an incompressible fluid. We find that simple beam theory does a good job of describing diatom chain deformation in a parabolic flow when its ends are tethered, but does not tell the whole story of chain deformations when they are subjected to compressive stresses in shear. While motivated by the fluid dynamics of diatom chains, our computational model of semiflexible fibres illustrates features that apply widely to other systems. The use of an adaptive immersed boundary framework allows us to capture complicated buckling and recovery dynamics of long, semiflexible fibres in shear.     

Disparity and convergence in bipedal archosaur locomotion

This study aims to investigate functional disparity in the locomotor apparatus of bipedal archosaurs. We use reconstructions of hindlimb myology of extant and extinct archosaurs to generate musculoskeletal biomechanical models to test hypothesized convergence between bipedal crocodile-line archosaurs and dinosaurs. Quantitative comparison of muscle leverage supports the inference that bipedal crocodile-line archosaurs and non-avian theropods had highly convergent hindlimb myology, suggesting similar muscular mechanics and neuromuscular control of locomotion. While these groups independently evolved similar musculoskeletal solutions to the challenges of parasagittally erect bipedalism, differences also clearly exist, particularly the distinct hip and crurotarsal ankle morphology characteristic of many pseudosuchian archosaurs. Furthermore, comparative analyses of muscle design in extant archosaurs reveal that muscular parameters such as size and architecture are more highly adapted or optimized for habitual locomotion than moment arms. The importance of these aspects of muscle design, which are not directly retrievable from fossils, warns against over-extrapolating the functional significance of anatomical convergences. Nevertheless, links identified between posture, muscle moments and neural control in archosaur locomotion suggest that functional interpretations of osteological changes in limb anatomy traditionally linked to postural evolution in Late Triassic archosaurs could be constrained through musculoskeletal modelling.     

A consistent muscle activation strategy underlies crawling and swimming in Caenorhabditis elegans

Although undulatory swimming is observed in many organisms, the neuromuscular basis for undulatory movement patterns is not well understood. To better understand the basis for the generation of these movement patterns, we studied muscle activity in the nematode Caenorhabditis elegans. Caenorhabditis elegans exhibits a range of locomotion patterns: in low viscosity fluids the undulation has a wavelength longer than the body and propagates rapidly, while in high viscosity fluids or on agar media the undulatory waves are shorter and slower. Theoretical treatment of observed behaviour has suggested a large change in force–posture relationships at different viscosities, but analysis of bend propagation suggests that short-range proprioceptive feedback is used to control and generate body bends. How muscles could be activated in a way consistent with both these results is unclear. We therefore combined automated worm tracking with calcium imaging to determine muscle activation strategy in a variety of external substrates. Remarkably, we observed that across locomotion patterns spanning a threefold change in wavelength, peak muscle activation occurs approximately 45° (1/8th of a cycle) ahead of peak midline curvature. Although the location of peak force is predicted to vary widely, the activation pattern is consistent with required force in a model incorporating putative length- and velocity-dependence of muscle strength. Furthermore, a linear combination of local curvature and velocity can match the pattern of activation. This suggests that proprioception can enable the worm to swim effectively while working within the limitations of muscle biomechanics and neural control.     

Logic programming to predict cell fate patterns and retrodict genotypes in organogenesis

Caenorhabditis elegans vulval development is a paradigm system for understanding cell differentiation in the process of organogenesis. Through temporal and spatial controls, the fate pattern of six cells is determined by the competition of the LET-23 and the Notch signalling pathways. Modelling cell fate determination in vulval development using state-based models, coupled with formal analysis techniques, has been established as a powerful approach in predicting the outcome of combinations of mutations. However, computing the outcomes of complex and highly concurrent models can become prohibitive. Here, we show how logic programs derived from state machines describing the differentiation of C. elegans vulval precursor cells can increase the speed of prediction by four orders of magnitude relative to previous approaches. Moreover, this increase in speed allows us to infer, or ‘retrodict’, compatible genomes from cell fate patterns. We exploit this technique to predict highly variable cell fate patterns resulting from dig-1 reduced-function mutations and let-23 mosaics. In addition to the new insights offered, we propose our technique as a platform for aiding the design and analysis of experimental data.     

Analytical model for instantaneous lift and shape deformation of an insect-scale flapping wing in hover

In the analysis of flexible flapping wings of insects, the aerodynamic outcome depends on the combined structural dynamics and unsteady fluid physics. Because the wing shape and hence the resulting effective angle of attack are a priori unknown, predicting aerodynamic performance is challenging. Here, we show that a coupled aerodynamics/structural dynamics model can be established for hovering, based on a linear beam equation with the Morison equation to account for both added mass and aerodynamic damping effects. Lift strongly depends on the instantaneous angle of attack, resulting from passive pitch associated with wing deformation. We show that both instantaneous wing deformation and lift can be predicted in a much simplified framework. Moreover, our analysis suggests that resulting wing kinematics can be explained by the interplay between acceleration-related and aerodynamic damping forces. Interestingly, while both forces combine to create a high angle of attack resulting in high lift around the midstroke, they offset each other for phase control at the end of the stroke.     

Scaling law and enhancement of lift generation of an insect-size hovering flexible wing

We report a comprehensive scaling law and novel lift generation mechanisms relevant to the aerodynamic functions of structural flexibility in insect flight. Using a Navier–Stokes equation solver, fully coupled to a structural dynamics solver, we consider the hovering motion of a wing of insect size, in which the dynamics of fluid–structure interaction leads to passive wing rotation. Lift generated on the flexible wing scales with the relative shape deformation parameter, whereas the optimal lift is obtained when the wing deformation synchronizes with the imposed translation, consistent with previously reported observations for fruit flies and honeybees. Systematic comparisons with rigid wings illustrate that the nonlinear response in wing motion results in a greater peak angle compared with a simple harmonic motion, yielding higher lift. Moreover, the compliant wing streamlines its shape via camber deformation to mitigate the nonlinear lift-degrading wing–wake interaction to further enhance lift. These bioinspired aeroelastic mechanisms can be used in the development of flapping wing micro-robots.     

The management of fluid and wave resistances by whirligig beetles

Whirligig beetles (Coleoptera: Gyrinidae) are semi-aquatic insects with a morphology and propulsion system highly adapted to their life at the air–water interface. When swimming on the water surface, beetles are subject to both fluid resistance and wave resistance.The purpose of this study was to analyse swimming speed, leg kinematics and the capillarity waves produced by whirligig beetles on the water surface in a simple environment. Whirligig beetles of the species Gyrinus substriatus were filmed in a large container, with a high-speed camera. Resistance forces were also estimated.These beetles used three types of leg kinematics, differing in the sequence of leg strokes: two for swimming at low speed and one for swimming at high speed. Four main speed patterns were produced by different combinations of these types of leg kinematics, and the minimum speed for the production of surface waves (23 cm s⁻¹) corresponded to an upper limit when beetles used low-speed leg kinematics. Each type of leg kinematics produced characteristic capillarity waves, even if the beetles moved at a speed below 23 cm s⁻¹. Our results indicate that whirligig beetles use low- and high-speed leg kinematics to avoid maximum drag and swim at speed corresponding to low resistances.     

The interplay of cell–cell and cell–substrate adhesion in collective cell migration

Collective cell migration often involves notable cell–cell and cell–substrate adhesions and highly coordinated motion of touching cells. We focus on the interplay between cell–substrate adhesion and cell–cell adhesion. We show that the loss of cell-surface contact does not significantly alter the dynamic pattern of protrusions and retractions of fast migrating amoeboid cells (Dictyostelium discoideum), but significantly changes their ability to adhere to other cells. Analysis of the dynamics of cell shapes reveals that cells that are adherent to a surface may coordinate their motion with neighbouring cells through protrusion waves that travel across cell–cell contacts. However, while shape waves exist if cells are detached from surfaces, they do not couple cell to cell. In addition, our investigation of actin polymerization indicates that loss of cell-surface adhesion changes actin polymerization at cell–cell contacts. To further investigate cell–cell/cell–substrate interactions, we used optical micromanipulation to form cell–substrate contact at controlled locations. We find that both cell-shape dynamics and cytoskeletal activity respond rapidly to the formation of cell–substrate contact.     

Remote automated multi-generational growth and observation of an animal in low Earth orbit

The ultimate survival of humanity is dependent upon colonization of other planetary bodies. Key challenges to such habitation are (patho)physiologic changes induced by known, and unknown, factors associated with long-duration and distance space exploration. However, we currently lack biological models for detecting and studying these changes. Here, we use a remote automated culture system to successfully grow an animal in low Earth orbit for six months. Our observations, over 12 generations, demonstrate that the multi-cellular soil worm Caenorhabditis elegans develops from egg to adulthood and produces progeny with identical timings in space as on the Earth. Additionally, these animals display normal rates of movement when fully fed, comparable declines in movement when starved, and appropriate growth arrest upon starvation and recovery upon re-feeding. These observations establish C. elegans as a biological model that can be used to detect changes in animal growth, development, reproduction and behaviour in response to environmental conditions during long-duration spaceflight. This experimental system is ready to be incorporated on future, unmanned interplanetary missions and could be used to study cost-effectively the effects of such missions on these biological processes and the efficacy of new life support systems and radiation shielding technologies.     

Network-based vaccination improves prospects for disease control in wild chimpanzees

Many endangered wildlife populations are vulnerable to infectious diseases for which vaccines exist; yet, pragmatic considerations often preclude large-scale vaccination efforts. These barriers could be reduced by focusing on individuals with the highest contact rates. However, the question then becomes whether targeted vaccination is sufficient to prevent large outbreaks. To evaluate the efficacy of targeted wildlife vaccinations, we simulate pathogen transmission and control on monthly association networks informed by behavioural data from a wild chimpanzee community (Kanyawara N = 37, Kibale National Park, Uganda). Despite considerable variation across monthly networks, our simulations indicate that targeting the most connected individuals can prevent large outbreaks with up to 35% fewer vaccines than random vaccination. Transmission heterogeneities might be attributed to biological differences among individuals (e.g. sex, age, dominance and family size). Thus, we also evaluate the effectiveness of a trait-based vaccination strategy, as trait data are often easier to collect than interaction data. Our simulations indicate that a trait-based strategy can prevent large outbreaks with up to 18% fewer vaccines than random vaccination, demonstrating that individual traits can serve as effective estimates of connectivity. Overall, these results suggest that fine-scale behavioural data can help optimize pathogen control efforts for endangered wildlife.     

Planar covariation of limb elevation angles during bipedal walking in the Japanese macaque

We investigated the planar covariation of lower limb segment elevation angles during bipedal walking in macaques to elucidate the mechanisms underlying the origin and evolution of the planar law in human walking. Two Japanese macaques and four adult humans walking on a treadmill were recorded, and the time course of the elevation angles at the thigh, shank and foot segments relative to the vertical axis were calculated. Our analyses indicated that the planar law also applies to macaque bipedal walking. However, planarity was much lower in macaques, and orientations of the plane differed between the two species because of differences in the foot elevation angle. The human foot is rigidly structured to form a longitudinal arch, whereas the macaque''s foot is more flexible and bends at the midtarsal region in the stance phase. This difference in midfoot flexibility between the two species studied was the main source of the difference in the planar law. Thus, the evolution of a stable midfoot in early hominins may have preceded the acquisition of the strong planar intersegmental coordination and possibly facilitated the subsequent emergence of habitual bipedal walking in the human lineage.     

Searching for motifs in the behaviour of larval Drosophila melanogaster and Caenorhabditis elegans reveals continuity between behavioural states

We present a novel method for the unsupervised discovery of behavioural motifs in larval Drosophila melanogaster and Caenorhabditis elegans. A motif is defined as a particular sequence of postures that recurs frequently. The animal''s changing posture is represented by an eigenshape time series, and we look for motifs in this time series. To find motifs, the eigenshape time series is segmented, and the segments clustered using spline regression. Unlike previous approaches, our method can classify sequences of unequal duration as the same motif. The behavioural motifs are used as the basis of a probabilistic behavioural annotator, the eigenshape annotator (ESA). Probabilistic annotation avoids rigid threshold values and allows classification uncertainty to be quantified. We apply eigenshape annotation to both larval Drosophila and C. elegans and produce a good match to hand annotation of behavioural states. However, we find many behavioural events cannot be unambiguously classified. By comparing the results with ESA of an artificial agent''s behaviour, we argue that the ambiguity is due to greater continuity between behavioural states than is generally assumed for these organisms.     

The Use of Structural Templates in Protein Backbone Modeling

The procedures used to model a protein structure are well established when the novel protein has high sequence similarity to a protein of known structure. Many proteins of interest have low (i.e. <50%) sequence similarity to any known structure. In these cases new approaches to prediction of structure are required.The use of sequence profiles which relate sequence to known structure has been proposed as one method to assign local regions of structure. As a first stage, templates or “icons” of the many relevant substructural motifs found in proteins must be defined. The sequences which gave rise to these structures are then aligned and a weighted profile obtained.Average structures of the 8 and 12 residue helix-turn and turn-helix motifs have been prepared. These coordinate templates were then used to scan through the Brookhaven protein structural database for similar, superimposable fragments. A composite template of 100 similar fragments for each element was found to be internally consistent to a rmsd=0.92 Å for HT8, 1.54 Å for HT12, 0.41 Å for TH8 and 1.40 Å for TH12. All of the sequences, from these structures, were then used to create an overall sequence profile.The four sequence profiles were scanned against the amino acid sequences of the proteins in the Brookhaven database: tertiary structure was correctly identified only about 10% of the time. This value is too low for predictive purposes. However, it could be increased by checking for multiple occurrences of the template in one protein.     

Experimental demonstration of a two-phase population extinction hazard

Population extinction is a fundamental biological process with applications to ecology, epidemiology, immunology, conservation biology and genetics. Although a monotonic relationship between initial population size and mean extinction time is predicted by virtually all theoretical models, attempts at empirical demonstration have been equivocal. We suggest that this anomaly is best explained with reference to the transient properties of ensembles of populations. Specifically, we submit that under experimental conditions, many populations escape their initially vulnerable state to reach quasi-stationarity, where effects of initial conditions are erased. Thus, extinction of populations initialized far from quasi-stationarity may be exposed to a two-phase extinction hazard. An empirical prediction of this theory is that the fit Cox proportional hazards regression model for the observed survival time distribution of a group of populations will be shown to violate the proportional hazards assumption early in the experiment, but not at later times. We report results of two experiments with the cladoceran zooplankton Daphnia magna designed to exhibit this phenomenon. In one experiment, habitat size was also varied. Statistical analysis showed that in one of these experiments a transformation occurred so that very early in the experiment there existed a transient phase during which the extinction hazard was primarily owing to the initial population size, and that this was gradually replaced by a more stable quasi-stationary phase. In the second experiment, only habitat size unambiguously displayed an effect. Analysis of data pooled from both experiments suggests that the overall extinction time distribution in this system results from the mixture of extinctions during the initial rapid phase, during which the effects of initial population size can be considerable, and a longer quasi-stationary phase, during which only habitat size has an effect. These are the first results, to our knowledge, of a two-phase population extinction process.     

Constraint-based network model of pathogen–immune system interactions

Pathogenic bacteria such as Bordetella bronchiseptica modulate host immune responses to enable their establishment and persistence; however, the immune response is generally successful in clearing these bacteria. Here, we model the dynamic outcome of the interplay between host immune components and B. bronchiseptica virulence factors. The model extends our previously published interaction network of B. bronchiseptica and includes the existing experimental information on the relative timing of IL10 and IFNγ activation in the form of qualitative inequalities. The current model improves the previous one in two directions: (i) by augmenting the network with new nodes with specific function in T helper cell differentiation and effector mechanisms and (ii) by using a dynamic approach that allows us to quantify node states and mechanisms revealed to be important from our previous model. The model makes predictions about the time scales of each process, the activity thresholds of each node and novel regulatory interactions. For example, the model predicts that the activity threshold of IL4 is higher than that of IL12 and that pro-inflammatory cytokines regulate the activity of Th2 cells. Some of these predictions are supported by the literature, and many can serve as targets of future experiments.     

Single cell profiling of surface carbohydrates on Bacillus cereus

Cell surface carbohydrates are important to various bacterial activities and functions. It is well known that different types of Bacillus display heterogeneity of surface carbohydrate compositions, but detection of their presence, quantitation and estimation of variation at the single cell level have not been previously solved. Here, using atomic force microscopy (AFM)-based recognition force mapping coupled with lectin probes, the specific carbohydrate distributions of N-acetylglucosamine and mannose/glucose were detected, mapped and quantified on single B. cereus surfaces at the nanoscale across the entire cell. Further, the changes of the surface carbohydrate compositions from the vegetative cell to spore were shown. These results demonstrate AFM-based ‘recognition force mapping’ as a versatile platform to quantitatively detect and spatially map key bacterial surface biomarkers (such as carbohydrate compositions), and monitor in situ changes in surface biochemical properties during intracellular activities at the single cell level.     

Variation in host susceptibility and infectiousness generated by co-infection: the myxoma-Trichostrongylus retortaeformis case in wild rabbits.

One of the conditions that can affect host susceptibility and parasite transmission is the occurrence of concomitant infections. Parasites interact directly or indirectly within an individual host and often these interactions are modulated by the host immune response. We used a free-living rabbit population co-infected with the nematode Trichostrongylus retortaeformis, which appears to stimulate an acquired immune response, and the immunosuppressive poxvirus myxoma. Modelling was used to examine how myxoma infection alters the immune-mediated establishment and death/expulsion of T. retortaeformis, and consequently affects parasite intensity and duration of the infection. Simulations were based on the general TH1-TH2 immunological paradigm that proposes the polarization of the host immune response towards one of the two subsets of T helper cells. Our findings suggest that myxoma infections contribute to alter host susceptibility to the nematode, as co-infected rabbits showed higher worm intensity compared with virus negative hosts. Results also suggest that myxoma disrupts the ability of the host to clear T. retortaeformis as worm intensities were consistently high and remained high in old rabbits. However, the co-infection model has to include some immune-mediated nematode regulation to be consistent with field data, indicating that the TH1-TH2 dichotomy is not complete. We conclude that seasonal myxoma outbreaks enhance host susceptibility to the nematode and generate highly infected hosts that remain infectious for a longer time. Finally, the virus-nematode co-infection increases heterogeneities among individuals and potentially has a large effect on parasite transmission.     

Proposed extension of the SINTAP/FITNET thin wall option based on a simple method for reference load determination

A “correct” limit or yield load is an essential element of flaw assessment procedures of the R6 or SINTAP/FITNET type. In this paper the authors propose a definition of this quantity which is based on the SINTAP option 3 failure assessment function. This “reference load” can be determined for any component geometry by finite element analysis. The method is applied to two kinds of thin walled structures (notched plates and curved stiffened panels) and the results demonstrate that the approach is a suitable extension of the existing thin wall module.     

The attachment strategy of English ivy: a complex mechanism acting on several hierarchical levels

English ivy (Hedera helix L.) is able to grow on vertical substrates such as trees, rocks and house plaster, thereby attaching so firmly to the surface that when removed by force typically whole pieces of the climbing substrate are torn off. The structural details of the attachment process are not yet entirely understood. We studied the attachment process of English ivy in detail and suggest a four-phase process to describe the attachment strategy: (i) initial physical contact, (ii) form closure of the root with the substrate, (iii) chemical adhesion, and (iv) shape changes of the root hairs and form-closure with the substrate. These four phases and their variations play an important role in the attachment to differently structured surfaces. We demonstrate that, in English ivy, different mechanisms work together to allow the plant''s attachment to various climbing substrates and reveal the importance of micro-fibril orientation in the root hairs for the attachment based on structural changes at the subcellular level.