Hydraulic integration and shrub growth form linked across continental aridity gradients.

Both engineered hydraulic systems and plant hydraulic systems are protected against failure by resistance, reparability, and redundancy. A basic rule of reliability engineering is that the level of independent redundancy should increase with increasing risk of fatal system failure. Here we show that hydraulic systems of plants function as predicted by this engineering rule. Hydraulic systems of shrubs sampled along two transcontinental aridity gradients changed with increasing aridity from highly integrated to independently redundant modular designs. Shrubs in humid environments tend to be hydraulically integrated, with single, round basal stems, whereas dryland shrubs typically have modular hydraulic systems and multiple, segmented basal stems. Modularity is achieved anatomically at the vessel-network scale or developmentally at the whole-plant scale through asymmetric secondary growth, which results in a semiclonal or clonal shrub growth form that appears to be ubiquitous in global deserts.

Primary motor cortex tuning to intended movement kinematics in humans with tetraplegia.

The relationship between spiking activities in motor cortex and movement kinematics has been well studied in neurologically intact nonhuman primates. We examined the relationship between spiking activities in primary motor cortex (M1) and intended movement kinematics (position and velocity) using 96-microelectrode arrays chronically implanted in two humans with tetraplegia. Study participants were asked to perform two different tasks: imagined pursuit tracking of a cursor moving on a computer screen and a "neural cursor center-out" task in which cursor position was controlled by the participant's neural activity. In the pursuit tracking task, the majority of neurons were significantly tuned: 90% were tuned to velocity and 86% were tuned to position in one participant; 95% and 84%, respectively, in the other. Additionally, velocity and position of the tracked cursor could be decoded from the ensemble of neurons. In the neural cursor center-out task, tuning to direction of the intended target was well captured by a log-linear cosine function. Neural spiking soon after target appearance could be used to classify the intended target with an accuracy of 95% in one participant, and 80% in the other. It was also possible to extract information about the direction of the difference vector between the target position and the instantaneous neural cursor position. Our results indicate that correlations between spiking activity and intended movement velocity and position are present in human M1 after the loss of descending motor pathways, and that M1 spiking activities share many kinematic tuning features whether movement is imagined by humans with tetraplegia, or is performed as shown previously in able-bodied nonhuman primates.

Biomechanics of the vibrissa motor plant in rat: rhythmic whisking consists of triphasic neuromuscular activity.

The biomechanics of a motor plant constrain the behavioral strategies that an animal has available to extract information from its environment. We used the rat vibrissa system as a model for active sensing and determined the pattern of muscle activity that drives rhythmic exploratory whisking. Our approach made use of electromyography to measure the activation of all relevant muscles in both head-fixed and unrestrained rats and two-dimensional imaging to monitor the position of the vibrissae in head-fixed rats. Our essential finding is that the periodic motion of the vibrissae and mystacial pad during whisking results from three phases of muscle activity. First, the vibrissae are thrust forward as the rostral extrinsic muscle, musculus (m.) nasalis, contracts to pull the pad and initiate protraction. Second, late in protraction, the intrinsic muscles pivot the vibrissae farther forward. Third, retraction involves the cessation of m. nasalis and intrinsic muscle activity and the contraction of the caudal extrinsic muscles m. nasolabialis and m. maxillolabialis to pull the pad and the vibrissae backward. We developed a biomechanical model of the whisking motor plant that incorporates the measured muscular mechanics along with movement vectors observed from direct muscle stimulation in anesthetized rats. The results of simulations of the model quantify how the combination of extrinsic and intrinsic muscle activity leads to an enhanced range of vibrissa motion than would be available from the intrinsic muscles alone.

Virtual lesions of the anterior intraparietal area disrupt goal-dependent on-line adjustments of grasp.

Adaptive motor behavior requires efficient error detection and correction. The posterior parietal cortex is critical for on-line control of reach-to-grasp movements. Here we show a causal relationship between disruption of cortical activity within the anterior intraparietal sulcus (aIPS) by transcranial magnetic stimulation (TMS) and disruption of goal-directed prehensile actions (either grip size or forearm rotation, depending on the task goal, with reaching preserved in either case). Deficits were elicited by applying TMS within 65 ms after object perturbation, which attributes a rapid control process on the basis of visual feedback to aIPS. No aperture deficits were produced when TMS was applied to a more caudal region within the intraparietal sulcus, to the parieto-occipital complex (putative V6, V6A) or to the hand area of primary motor cortex. We contend that aIPS is critical for dynamic error detection during goal-dependent reach-to-grasp action that is visually guided.

Effects of hydration on mechanical properties of a highly sclerotized tissue.

The jaws of the bloodworm Glycera dibranchiata consist principally of protein and melanin scaffolds with small amounts of unmineralized copper (Cu) and mineralized atacamite (Cu(2)Cl(OH)(3)) fibers in distinct regions. Remarkably, when tested in air, the regions containing unmineralized Cu are the hardest, stiffest, and most abrasion resistant. To establish the functions of jaw constituents in physiologically relevant environments, this study examines the effects of hydration on their response to indentation, scratching, and wear. Although all jaw regions are degraded by the presence of water, the ones containing unmineralized Cu are affected least. Notably, scratch depths in the bulk and the atacamite-containing regions double when wet, whereas the corresponding increase in the regions with unmineralized Cu is approximately 20%. The results support the view that Cu ions are involved in the formation of intermolecular coordination complexes, creating a cross-linked molecular network that is both mechanically robust and resistant to water ingress. Hydration effects are greatest during wear testing, rates of material removal in water being about three times those in air. The mechanism underlying accelerated wear is suspected to involve coupled effects of near-surface damage and enhanced water ingress, resulting in increased plasticization and susceptibility to plastic plowing.

Cytoskeletal bundle mechanics.

The mechanical properties of cytoskeletal actin bundles play an essential role in numerous physiological processes, including hearing, fertilization, cell migration, and growth. Cells employ a multitude of actin-binding proteins to actively regulate bundle dimensions and cross-linking properties to suit biological function. The mechanical properties of actin bundles vary by orders of magnitude depending on diameter and length, cross-linking protein type and concentration, and constituent filament properties. Despite their importance to cell function, the molecular design principles responsible for this mechanical behavior remain unknown. Here, we examine the mechanics of cytoskeletal bundles using a molecular-based model that accounts for the discrete nature of constituent actin filaments and their distinct cross-linking proteins. A generic competition between filament stretching and cross-link shearing determines three markedly different regimes of mechanical response that are delineated by the relative values of two simple design parameters, revealing the universal nature of bundle-bending mechanics. In each regime, bundle-bending stiffness displays distinct scaling behavior with respect to bundle dimensions and molecular composition, as observed in reconstituted actin bundles in vitro. This mechanical behavior has direct implications on the physiological bending, buckling, and entropic stretching behavior of cytoskeletal processes, as well as reconstituted actin systems. Results are used to predict the bending regimes of various in vivo cytoskeletal bundles that are not easily accessible to experiment and to generate hypotheses regarding implications of the isolated behavior on in vivo bundle function.

Quantitative mechanical evaluation and analysis of Drosophila embryos through the stages of embryogenesis.

The fruit fly Drosophila embryo is one of the most important model organisms in genetics and developmental biology research. To better understand the biomechanical properties involved in Drosophila embryo research, this work presents a mechanical characterization of living Drosophila embryos through the stages of embryogenesis. Measurements of the mechanical forces of Drosophila embryos are implemented using a novel, in situ, and minimally invasive force sensing tool with a resolution in the range of microN. The measurements offer an essential understanding of penetration force profiles during the microinjection of Drosophila embryos. Sequentially quantitative evaluation and analysis of the mechanical properties, such as Young's modulus, stiffness, and mechanical impedance of living Drosophila embryos are performed by extracting the force measurements throughout the stages of embryogenesis. Experimental results illustrate the changing mechanical properties of Drosophila embryos during development, and thus mathematical models are proposed. The evaluation provides a critical step toward better understanding of the biomechanical properties of Drosophila embryos during embryogenesis, and could contribute to more efficient and significant genetic and embryonic development research on Drosophila.

Choice of contact points during multidigit grasping: effect of predictability of object center of mass location.

It has been shown that when subjects can predict object properties [e.g., weight or center of mass (CM)], fingertip forces are appropriately scaled before the object is lifted, i.e., before somatosensory feedback can be processed. However, it is not known whether subjects, in addition to these anticipatory force mechanisms, exploit the ability to choose where digits can be placed to facilitate object manipulation. We addressed this question by asking subjects to reach and grasp an object whose CM was changed to the left, center, or right of the object in either a predictable or unpredictable manner. The only task requirement was to minimize object roll during lift. We hypothesized that subjects would modulate contact points but only when object CM location could be predicted. As expected, object roll was significantly smaller in the predictable condition. This experimental condition was also associated with statistically distinct spatial distributions of contact points as a function of object CM location but primarily when large torques had to be counteracted, i.e., for right and left CM locations. In contrast, when subjects could not anticipate CM location, a "default" distribution of contact points was used, this being statistically indistinguishable from that adopted for the center CM location in the predictable condition. We conclude that choice of contact points is integrated with anticipatory force control mechanisms to facilitate object manipulation. These results demonstrate that planning of digit placement is an important component of grasp control.

Modular control of limb movements during human locomotion.

The idea that the CNS may control complex interactions by modular decomposition has received considerable attention. We explored this idea for human locomotion by examining limb kinematics. The coordination of limb segments during human locomotion has been shown to follow a planar law for walking at different speeds, directions, and levels of body unloading. We compared the coordination for different gaits. Eight subjects were asked to walk and run on a treadmill at different speeds or to walk, run, and hop over ground at a preferred speed. To explore various constraints on limb movements, we also recorded stepping over an obstacle, walking with the knees flexed, and air-stepping with body weight support. We found little difference among covariance planes that depended on speed, but there were differences that depended on gait. In each case, we could fit the planar trajectories with a weighted sum of the limb length and orientation trajectories. This suggested that limb length and orientation might provide independent predictors of limb coordination. We tested this further by having the subjects step, run, and hop in place, thereby varying only limb length and maintaining limb orientation fixed, and also by marching with knees locked to maintain limb length constant while varying orientation. The results were consistent with a modular control of limb kinematics where limb movements result from a superposition of separate length- and orientation-related angular covariance. The hypothesis finds support in the animal findings that limb proprioception may also be encoded in terms of these global limb parameters.

Cancellous bone lamellae strongly affect microcrack propagation and apparent mechanical properties: separation of patients with osteoporotic fracture from normal controls using a 2D nonlinear finite element method (biomechanical stereology).

Biomechanical stereology is proposed as a two-dimensional (2D) finite element (FE) method to estimate the ability of bone tissue to sustain damage and to separate patients with osteoporotic fracture from normal controls. Briefly, 2D nonlinear compact tension FE models were created from quantitative back scattered electron images taken of iliac crest bone specimens collected from the individuals with or without osteoporotic fracture history. The effects of bone mineral microstructure on predicted bone fracture toughness and microcrack propagation were examined. The 2D FE models were used as surrogates for the real bone tissues. The calculated microcrack propagation results and bone mechanical properties were examined as surrogates for measurements from mechanical testing of actual specimens. The results for the 2D FE simulation separated patients with osteoporotic fracture from normal controls even though only the variability in tissue mineral microstructure was used to build the models. The models were deliberately created to ignore all differences in mean mineralization. Hence, the current results support the following hypotheses: (1) that material heterogeneity is important to the separation of patients with osteoporotic fracture from normal controls; and (2) that 2D nonlinear finite element modeling can produce surrogate mechanical parameters that separate patients with fracture from normal controls.

Relating osteon diameter to strain.

Osteon diameter is generally smaller in bone regions that experience larger strains. A mechanism relating osteon diameter to strain is as yet unknown. We propose that strain-induced osteocyte signals inhibit osteoclastic bone resorption. This mechanism was previously shown to produce load-aligned osteons in computer simulations. Now we find that it also predicts smaller osteon diameter for higher loads. Additionally, we find that our model predicts osteon development with two cutting cones, one moving up and one moving down the loading axis. Such 'double-ended osteons' were reported in literature as a common type of osteon development. Further, we find that a steep gradient in strain magnitude can result in an osteonal tunnel with continuous resorption along the less strained side, which corresponds to 'drifting osteons' reported in literature.

Three-dimensional microstructure of the bone in a hamster model of senile osteoporosis.

Age-related bone loss, which is poorly characterized, is a major underlying cause of osteoporotic fractures in the elderly. In order to identify the morphological feature of age-related bone loss, we investigated sex and site (tibia, femur and vertebra) dependence of bone microstructure in aging hamsters from 3 to 24 months of age using micro-CT. In the proximal tibia and distal femur, trabecular bone volume (BV/TV), trabecular number (Tb.N), trabecular thickness (Tb.Th) and bone mineral density (BMD) increased to a maximum at 6 or 12 months and then declined progressively from 12 to 24 months of age. Trabecular separation (Tb.Sp), trabecular bone pattern factor (TBPf) and structure model index (SMI) increased with age. As compared with male hamsters, BV/TV and Tb.N were significantly lower in females at 18 and 24 months of age. Age-related decrease of trabecular BV/TV in the vertebral body was less than that of the femoral and tibial metaphyses. In the mid-femoral diaphysis, cortical bone area remained constant from 3 to 24 months of age. Cortical thickness decreased from 12 to 24 months and cortical BMD declined significantly from 18 to 24 months of age. These findings indicate that skeletal site and sex differences exist in hamster bone structure. Age-related bone changes in hamsters resemble those in humans. We conclude that hamster may be a useful model to study at least some aspects of bone loss during human aging.

Differences in matrix composition between calvaria and long bone in mice suggest differences in biomechanical properties and resorption: Special emphasis on collagen.

The mammalian skeleton consists of bones that are formed in two different ways: long bones via endochondral ossification and flat bones via intramembranous ossification. These different formation modes may result in differences in the composition of the two bone types. Using the 2D-difference in gel electrophoresis technique and mass spectrometry, we analyzed the composition of murine mineral-associated proteins of calvaria and long bone. Considerable differences in protein composition were observed. Flat bones (calvariae) contained more soluble collagen (8x), pigment epithelium derived factor (3x) and osteoglycin (4x); whereas long bones expressed more chondrocalcin (3x), thrombospondin- 1 (4x), fetuin (4x), secreted phosphoprotein 24 (3x), and thrombin (7x). Although cystatin motifs containing proteins, such as secreted phosphoprotein 24 and fetuin are highly expressed in long bone, they did not inhibit the activity of the cysteine proteinases cathepsin B and K. The solubility of collagen differed which coincided with differences in collagen crosslinking, long bone containing 3x more (hydroxylysine)-pyridinoline. The degradation of long bone collagen by MMP2 (but not by cathepsin K) was impaired. These differences in collagen crosslinking may explain the differences in the proteolytic pathways osteoclasts use to degrade bone. Our data demonstrate considerable differences in protein composition of flat and long bones and strongly suggest functional differences in formation, resorption, and mechanical properties of these bone types.

Do the left atrial substrate properties correlate with the left atrial mechanical function? A novel insight from the electromechanical study in patients with atrial fibrillation.

BACKGROUND: The atrial substrate is the determinant of occurrence and maintenance of atrial fibrillation (AF), which can induce remodeling of atrial function and structure. This study investigated the relationship between the left atrial (LA) substrate properties and LA mechanical function. METHODS: Forty-four consecutive patients (50.3 +/- 10.7 years old, 33 men) who presented with sinus rhythm during echocardiographic study before receiving catheter ablation for AF were enrolled. The LA diameter, LA volume, ratio of early and late transmitral filling flow velocities (E/A), LA appendage flow velocity, and transmitral velocity-time integral (VTI) were measured by the echocardiography. The LA empty fraction (LAEF), which was obtained via dividing the difference between maximal and minimal LA volume by maximal LA volume, was calculated as a parameter of the global LA contractile function. The LA global contact voltage mapping (NavX system) was performed before pulmonary vein isolation. RESULTS: Mean LA voltage and LA low voltage zone index (LVZ index, area with voltage < 0.5 mV, divided by total LA surface area) showed significant correlation with LA diameter and volume, but only the LA LVZ index showed significant correlation with A-wave velocity, transmitral A-wave VTI, and LAEF (r =-0.340, -0.411, -0.426; P = 0.024, 0.006, 0.005, respectively). We divided the LA LVZ index into three groups (< 10%, 10-20%, > 20%). The LAEF got worse and the transmitral A-wave VTI percentage (divided by transmitral VTI) decreased as LA LVZ index increased. CONCLUSIONS: The LA substrate properties showed close correlation with LA size, but only the LA LVZ index correlated with the LA mechanical function.

Hydrodynamic simulation of multicellular embryo invagination.

The mechanical aspects of embryonic morphogenesis have been widely analysed by numerical simulations of invagination in sea urchins and Drosophila gastrulation. Finite element models, which describe the tissue as a continuous medium, lead to the global invagination morphogenesis observed in vivo. Here we develop a simulation of multicellular embryo invagination that allows access to both cellular and multicellular mechanical behaviours of the embryo. In this model, the tissue is composed of adhesive individual cells, in which shape change dynamics is governed by internal acto-myosin forces and the hydrodynamic flow associated with membrane movements. We investigated the minimal structural and force elements sufficient to phenocopy mesoderm invagination. The minimal structures are cell membranes characterized by an acto-myosin cortical tension and connected by apical and basal junctions and an acto-myosin contractile ring connected to the apical junctions. An increase in the apical-cortical surface tension is the only control parameter change required to phenocopy most known multicellular and cellular shape changes of Drosophila gastrulation. Specifically, behaviours observed in vivo, including apical junction movements at the onset of gastrulation, cell elongation and subsequent shortening during invagination, and the development of a dorso-ventral gradient of thickness of the embryo, are predicted by this model as passive mechanical consequences of the genetically controlled increase in the apical surface tension in invaginating mesoderm cells, thus demonstrating the accurate description of structures at both global and single cell scales.

The forces that shape embryos: physical aspects of convergent extension by cell intercalation.

We discuss the physical aspects of the morphogenic process of convergence (narrowing) and extension (lengthening) of tissues by cell intercalation. These movements, often referred to as 'convergent extension', occur in both epithelial and mesenchymal tissues during embryogenesis and organogenesis of invertebrates and vertebrates, and they play large roles in shaping the body plan during development. Our focus is on the presumptive mesodermal and neural tissues of the Xenopus (frog) embryo, tissues for which some physical measurements have been made. We discuss the physical aspects of how polarized cell motility, oriented along future tissue axes, generate the forces that drive oriented cell intercalation and how this intercalation results in convergence and extension or convergence and thickening of the tissue. Our goal is to identify aspects of these morphogenic movements for further biophysical, molecular and cell biological, and modeling studies.

Rayleigh instability of the inverted one-cell amphibian embryo.

The one-cell amphibian embryo is modeled as a rigid spherical shell containing equal volumes of two immiscible fluids with different densities and viscosities and a surface tension between them. The fluids represent denser yolk in the bottom hemisphere and clearer cytoplasm and the germinal vesicle in the top hemisphere. The unstable equilibrium configuration of the inverted system (the heavier fluid on top) depends on the value of the contact angle. The theoretically calculated normal modes of perturbation and the instability of each mode are in agreement with the results from ComFlo computational fluid dynamic simulations of the same system. The two dominant types of modes of perturbation give rise to axisymmetric and asymmetric sloshing of the cytoplasm of the inverted embryos, respectively. This work quantifies our hypothesis that the axisymmetric mode corresponds to failure of development, and the asymmetric sloshing mode corresponds to development proceeding normally, but with reversed pigmentation, for inverted embryos.

Comparative anatomy and muscle architecture of selected hind limb muscles in the Quarter Horse and Arab.

The Quarter Horse (bred for acceleration) and the Arab (bred for endurance) are situated at either end of the equine athletic spectrum. Studies into the form and function of the leg muscles in human sprint and endurance runners have demonstrated that differences exist in their muscle architecture. It is not known whether similar differences exist in the horse. Six Quarter Horse and six Arab fresh hind limb cadavers were dissected to gain information on the muscle mass and architecture of the following muscles: gluteus medius; biceps femoris; semitendinosus; vastus lateralis; gastrocnemius; tibialis cranialis and extensor digitorum longus. Specifically, muscle mass, fascicle length and pennation angle were quantified and physiological cross-sectional area (PCSA) and maximum isometric force were estimated. The hind limb muscles of the Quarter Horse were of a significantly greater mass, but had similar fascicle lengths and pennation angles when compared with those of the Arab; this resulted in the Quarter Horse hind limb muscles having greater PCSAs and hence greater isometric force potential. This study suggests that Quarter Horses as a breed inherently possess large strong hind limb muscles, with the potential to accelerate their body mass more rapidly than those of the Arab.

Constellation of phase singularities in a speckle-like pattern for optical vortex metrology applied to biological kinematic analysis.

A novel technique for biological kinematic analysis is proposed that makes use of the pseudophase singularities in a complex signal generated from a speckle-like pattern. In addition to the information about the locations and the anisotropic core structures of the pseudophase singularities, we also detect the spatial structures of a cluster of phase singularities, which serves as a unique constellation characterizing the mutual position relation between the individual pseudophase singularities. Experimental results of in vivo measurements for a swimming fish along with its kinematic analysis are presented, which demonstrate the validity of the proposed technique.

Predictions of hepatic disposition properties using a mechanistically realistic, physiologically based model.

Quantitative mappings were established between drug physicochemical properties (PCPs) and parameter values of a physiologically based, mechanistically realistic, in silico liver (ISL). The ISL plugs together autonomous software objects that represent hepatic components at different scales and levels of detail. Microarchitectural features are represented separately from the mechanisms that influence drug metabolism. The same ISL has been validated against liver perfusion data for sucrose and four cationic drugs: antipyrine, atenolol, labetalol, and diltiazem. Parameters sensitive to drug-specific PCPs were tuned so that ISL outflow profiles from a single ISL matched in situ perfused rat liver outflow profiles of all five compounds. Quantitative relationships were then established between the four sets of drug PCPs and the corresponding four sets of PCP-sensitive, ISL parameter values; those relationships were used to predict PCP-sensitive, ISL parameter values for prazosin and propranolol given only their PCPs. Relationships were established using three different methods: 1) a simple linear correlation method, 2) the fuzzy c-means algorithm, and 3) a simple artificial neural network. Each relationship was used separately to predict ISL parameter values for prazosin and propranolol, given their PCPs. Those values were applied in the ISL used earlier to predict the hepatic disposition details for each drug. Although we had only sparse data available, all predicted disposition profiles were judged reasonable (within a factor of 2 of referent profile data). The order of precision, based on a similarity measure, was 3 > 2 > 1.