Response analysis of the lumbar spine during regular daily activities—A finite element analysis

A non-linear poroelastic finite element model of the lumbar spine was developed to investigate spinal response during daily dynamic physiological activities. Swelling was simulated by imposing a boundary pore pressure of 0.25 MPa at all external surfaces. Partial saturation of the disc was introduced to circumvent the negative pressures otherwise computed upon unloading. The loading conditions represented a pre-conditioning full day followed by another day of loading: 8 h rest under a constant compressive load of 350 N, followed by 16 h loading phase under constant or cyclic compressive load varying in between 1000 and 1600 N. In addition, the effect of one or two short resting periods in the latter loading phase was studied.The model yielded fairly good agreement with in-vivo and in-vitro measurements. Taking the partial saturation of the disc into account, no negative pore pressures were generated during unloading and recovery phase. Recovery phase was faster than the loading period with equilibrium reached in only ～3 h. With time and during the day, the axial displacement, fluid loss, axial stress and disc radial strain increased whereas the pore pressure and disc collagen fiber strains decreased. The fluid pressurization and collagen fiber stiffening were noticeable early in the morning, which gave way to greater compression stresses and radial strains in the annulus bulk as time went by. The rest periods dampened foregoing differences between the early morning and late in the afternoon periods. The forgoing diurnal variations have profound effects on lumbar spine biomechanics and risk of injury.     

Can dog-ear formation be decreased when an S-shaped skin resection is used instead of a spindle skin resection? A three-dimensional analysis of skin surgery techniques using the finite element method

Dog-ear formation is often unavoidable with resection and suturing of the skin, including spindle excision. Regarding dog-ear formation after basic spindle skin resection during removal of a round tumor of the skin, we quantitatively analyzed the frequency of dog-ear formation with respect to the following three techniques: previous spindle skin resection, S-shaped skin resection, which has been experientially considered to induce limited deformity, and mosque-shaped skin resection for control. To date, by using paper models or sponges, various techniques of skin resection have been simulated in the field of plastic surgery. In the present study, we performed three-dimensional simulation and analyzed three different techniques of skin resection by using the finite element method. As a result, image simulation demonstrated that the frequency of dog-ear formation was limited by S-shaped, spindle, and mosque-shaped skin resection, in descending order.     

Investigation of different cage designs and mechano-regulation algorithms in the lumbar interbody fusion process – A finite element analysis

Lumbar interbody fusion cages are commonly used to treat painful spinal degeneration and instability by achieving bony fusion. Many different cage designs exist, however the effect of cage morphology and material properties on the fusion process remains largely unknown. This finite element model study aims to investigate the influence of different cage designs on bone fusion using two mechano-regulation algorithms of tissue formation. It could be observed that different cages play a distinct key role in the mechanical conditions within the fusion region and therefore regulate the time course of the fusion process.     

Does the graft-tunnel friction influence knee joint kinematics and biomechanics after anterior cruciate ligament reconstruction? A finite element study

Graft tissues within bone tunnels remain mobile for a long time after anterior cruciate ligament (ACL) reconstruction. However, whether the graft-tunnel friction affects the finite element (FE) simulation of the ACL reconstruction is still unclear. Four friction coefficients (from 0 to 0.3) were simulated in the ACL-reconstructed joint model as well as two loading levels of anterior tibial drawer. The graft-tunnel friction did not affect joint kinematics and the maximal principal strain of the graft. By contrast, both the relative graft-tunnel motion and equivalent strain for the bone tunnels were altered, which corresponded to different processes of graft-tunnel integration and bone remodeling, respectively. It implies that the graft-tunnel friction should be defined properly for studying the graft-tunnel integration or bone remodeling after ACL reconstruction using numerical simulation.     

Do TRPV1 antagonists increase the risk for skin tumourigenesis? A collaborative in vitro and in vivo assessment

A recent hypothesis suggesting that the pharmacological target TRPV1 (transient receptor potential vanilloid subfamily, member 1) may function as a tumour suppressor, which potentially impacts the development of TRPV1 antagonist therapeutics for a range of conditions. However, little is known about the long-term physiologic effects of TRPV1 blockade in the skin. In vitro and in vivo studies suggested that the potent TRPV1 competitive antagonist AMG-9810 promoted proliferation in N/TERT1 cells (telomerase-immortalised primary human keratinocytes 1) and tumour development in mouse skin that was mediated through EGFR/Akt/mTOR signalling. We attempted to reproduce the reported in vitro and in vivo findings to further explore this hypothesis to understand the underlying mechanism and the risk associated with TRPV1 antagonism in the skin. In vitro proliferation studies using multiple methods and topical application with AMG-9810 and structurally similar TRPV1 antagonists such as SB-705498 and PAC-14028 were performed. Although we confirmed expression of TRPV1 in primary human epidermal keratinocytes (HEKn) and spontaneously immortalised human keratinocytes (HaCaT), we were unable to demonstrate cell proliferation in either cell type or any clear evidence of increased expression of proteins in the EGFR/Akt/mTOR signalling pathway with these molecules. We were also unable to demonstrate skin tumour promotion or underlying molecular mechanisms involved in the EGFR/Akt/mTOR signalling pathway in a single-dose and two-stage carcinogenesis mouse study treated with TRPV1 antagonists. In conclusion, our data suggest that inhibiting the pharmacological function of TRPV1 in skin by specific antagonists has not been considered to be indicative of skin tumour development.     

Do predators influence the behaviour of bats?

Many aspects of animal behaviour are affected by real‐time changes in the risk of predation. This conclusion holds for virtually all taxa and ecological systems studied, but does it hold for bats? Bats are poorly represented in the literature on anti‐predator behaviour, which may reflect a lack of nocturnal predators specialized on bats. If bats actually experience a world with minimal anti‐predator concerns, then they will provide a unique contrast within the realm of vertebrate ecology. Alternatively, such predator‐driven behaviour in bats may not yet be fully understood, given the difficulties in working with these highly mobile and nocturnal animals. We provide a wide‐ranging exploration of these issues in bat behaviour. We first cover the basic predator‐prey information available on bats, both on potential predators and the ways in which bats might perceive predators and respond to attacks. We then cover work relevant to key aspects of bat behaviour, such as choice of daytime roosts, the nature of sleep and torpor, evening roost departures, moonlight avoidance, landscape‐related movement patterns, and habitat selection. Overall, the evidence in favour of a strong influence of predators on bat behaviour is equivocal, with the picture clouded by contradictory results and a lack of information on potential predators and the perception of risk by bats. It seems clear that day‐active bats run a considerable risk of being killed by diurnal raptors, which are able to capture bats with relative ease. Thus, bats taking advantage of a pulse of insects just prior to sunset are likely taking risks to gain much‐needed energy. Further, the choice of daytime roosts by bats is probably strongly influenced by roost safety. Few studies, however, have directly addressed either of these topics. As a group, insectivorous temperate‐zone bats show no clear tendency to avoid apparently risky situations, such as activity on moonlit nights. However, some observations are consistent with the idea that predation risk affects choice of movement paths and feeding areas by temperate‐zone bats, as well as the timing of roost departures. The behaviour of tropical bats, on the other hand, seems more generally influenced by predators; this is especially true for tropical nectarivores and frugivores, but also for insectivorous bats. Presumably there are more serious predators on bats in the tropics (e.g. specialized raptors or carnivorous bats), but the identity of these predators is unclear. More information is needed to assess fully the influence of predators on bat behaviour. There is much need for work on the ways in which bats perceive predators via auditory, visual, and olfactory cues, and whether bats have some knowledge of the risks posed by different predators. Also needed is information on how predators attack bats and how bats react to attacking predators. Difficult to obtain, but of critical value, will be information on the nature of the predation risk experienced by bats while away from roosts and during the full darkness of night.     

A common element in epidermal expression?

Members of an epidermally expressed gene cluster on human chromosome 21 each contain a short sequence similar to an element that can drive ectoderm-specific gene expression.     

Can the modeling for simplification of a dental implant surface affect the accuracy of 3D finite element analysis?

The aim of this study was to assess stress/strain of different implant modeling simplifications by 3D-FEA. Three variation of external hexagon implant (Ø3.75?×?10 mm) supporting one molar crown were simulated: A (no threads); B (slightly threads simplification); C (original design). 200 N (axial) and 100 N (oblique) were applied. Cortical bone was evaluated by maximum principal stress and microstrain qualitatively and quantitatively (ANOVA and Tukey post hoc (p < 0.05)). Higher stress levels (p < 0.05) were observed in model A. Models B and C presented similar stress transmission. It was possible to conclude that slightly simplification should be used for studies evaluating stress transferring for bone tissue.     

Evaporative heat loss in Bos taurus: Do different cattle breeds cope with heat stress in the same way?

The aim of this study was to compare two Portuguese (Alentejana and Mertolenga) and two exotic (Frisian and Limousine) cattle breeds in terms of the relationship between the increase in ambient temperature and the responses of the evaporative heat loss pathways and the effects on homeothermy. In the experiment, six heifers of the Alentejana, Frisian, and Mertolenga breeds and four heifers of the Limousine breed were used. The animals were placed in four temperature levels, the first one under thermoneutral conditions and the other ones with increase levels of thermal stress. When submitted to severe heat stress, the Frisian developed high thermal tachypnea (125 mov/min) and moderate sweating rates (117 g m⁻² h⁻¹), which did not prevent an increase in the rectal temperature (from 38.4 °C to 40.0 °C). Moderate increases in rectal temperature were observed in the Alentejana (from 38.8 °C to 39.4 °C) and Limousine (from 38.6 °C to 39.4 °C), especially in the period of highest heat stress. The Limousine showed moderate levels of tachypnea (101 mov/min) while showing the lowest sweating rates. The Alentejana showed significant increases in sweating rate (156 g m⁻² h⁻¹) that played a major role in homeothermy. The Mertolenga showed a superior stability of body temperature, even in the period of highest heat stress (from 38.5 °C to 39.1 °C). Uncommonly, the maintenance of homeothermy during moderate heat stress was achieved primarily by intense tachypnea (122 mov/min). The sweating rate remained abnormally low under conditions of moderate heat stress, rising significantly (110 g m⁻² h⁻¹) without evidence of stabilization, only when tendency for heat storage occurred. This unusual response of the evaporative heat loss pathways infers a different thermoregulatory strategy, suggesting a different adaptation to semi-arid environment and strong association with water metabolism.     

Does screw–bone interface modelling matter in finite element analyses?

The effect of screw-bone interface modelling strategies was evaluated in the setting of a tibial mid-shaft fracture stabilised using locking plates. Three interface models were examined: fully bonded interface; screw with sliding contact with bone; and screw with sliding contact with bone in an undersized pilot hole. For the simulation of the last interface condition we used a novel thermal expansion approach to generate the pre-stress that the bone would be exposed to during screw insertion. The study finds that the global load-deformation response is not influenced by the interface modelling approach employed; the deformation varied by less than 1% between different interaction models. However, interface modelling is found to have a considerable impact on the local stress-strain environment within the bone in the vicinity of the screws. Frictional and tied representations did not have significantly different peak strain values (<5% difference); the frictional interface had higher peak compressive strains while the tied interface had higher tensile strains. The undersized pilot hole simulation produced the largest strains. The peak minimum principal strains for the frictional interface were 26% of those for the undersized pilot hole simulation at a load of 770 N. It is concluded that the commonly used tie constraint can be used effectively when the only interest is the global load-deformation behaviour. Different contact interface models, however, alter the mechanical response around screw holes leading to different predictions for screw loosening, bone damage and stress shielding.     

Search for critical loading condition of the spine–A meta analysis of a nonlinear viscoelastic finite element model

The relative vulnerability of spinal motion segments to different loading combinations remains unknown. The meta-analysis described here using the results of a validated L2–L3 nonlinear viscoelastic finite element model was designed to investigate the critical loading and its effect on the internal mechanics of the human lumbar spine. A Box-Behnken experimental design was used to design the magnitude of seven independent variables associated with loads, rotations and velocity of motion. Subsequently, an optimization method was used to find the primary and secondary variables that influence spine mechanical output related to facet forces, disc pressure, ligament forces, annulus matrix compressive/shear stresses and anulus fibers strain. The mechanical responses with respect to the two most-relevant variables were then regressed linearly using the response surface quadratic model. Axial force and sagittal rotation were identified as the most-relevant variables for mechanical responses. The procedure developed can be used to find the critical loading for finite element models with multi input variables. The derived meta-models can be used to predict the risk associated with various loading parameters and in setting safer load limits.     

Is the sheep a suitable model to study the mechanical alterations of disc degeneration in humans? A probabilistic finite element model study

Intervertebral disc degeneration is one major source of low back pain, which because of its complex multifactorial nature renders the treatment challenging and thus necessitates extensive research. Experimental animal models have proven valuable in improving our understanding of degenerative processes and potentially promising therapies. Currently, the sheep is the most frequently used large animal in vivo model in intervertebral disc research. However, despite its undoubted value for investigations of the complex biological and cellular aspects, to date, it is unclear whether the sheep is also suited to study the mechanical aspects of disc degeneration in humans.A parametric finite element (FE) model of the L4–5 spinal motion segment was developed. Using this model, the geometry and the material properties of both the human and the ovine spinal segment as well as different appearances of disc degeneration can be depicted. Under pure and combined loads, it was investigated whether degenerative changes to both the human and the ovine model equivalent caused the same mechanical response.Different patterns of degeneration resulted in large variations in the ranges of motion, intradiscal pressure, ligament and facet loads. In the human, but not in the ovine model, all these results differed significantly between different degrees of degeneration.This FE model study highlighted possible differences in the mechanical response to disc degeneration between human and ovine intervertebral discs and indicates the necessity of further, more detailed, investigations.     

The sharpest tools in the box? Quantitative analysis of conodont element functional morphology

Conodonts have been considered the earliest skeletonizing vertebrates and their mineralized feeding apparatus interpreted as having performed a tooth function. However, the absence of jaws in conodonts and the small size of their oropharyngeal musculature limits the force available for fracturing food items, presenting a challenge to this interpretation. We address this issue quantitatively using engineering approaches previously applied to mammalian dentitions. We show that the morphology of conodont food-processing elements was adapted to overcome size limitations through developing dental tools of unparalleled sharpness that maximize applied pressure. Combined with observations of wear, we also show how this morphology was employed, demonstrating how Wurmiella excavata used rotational kinematics similar to other conodonts, suggesting that this occlusal style is typical for the clade. Our work places conodont elements within a broader dental framework, providing a phylogenetically independent system for examining convergence and scaling in dental tools.     

Do rivers influence fine-scale population genetic structure of tigers in the Sundarbans?

Global tiger Panthera tigris populations mostly survive within the geographically fragmented forest patches, thereby limited genetic exchange between isolated populations. Assessing the genetic status of these populations can reveal the effects of dispersal barriers and provide critical insights to guide future conservation actions. Using non-invasively collected biological samples, we investigated fine-scale genetic structure of tigers in the Sundarbans mangrove forests intersected by the complex river systems, and which holds one of the largest global tiger populations. We genotyped 52 tiger samples at 10 polymorphic microsatellite loci, and sequenced 33 of them for a total of 1263 base-pairs at four mitochondrial gene fragments. Microsatellite analyses exhibit a signature of fine-scale genetic structure, which might have been the consequence of limited tiger dispersal due to wide rivers across the Sundarbans. Similarly, mitochondrial data show a historic pattern of population isolation that might be due to wider rivers across the entire Sundarbans shared by Bangladesh and India. Given the intrinsic nature of the mangrove habitat embedded with numerous rivers, increased commercial traffic and human activities may further impede tiger dispersal across wide rivers, escalating further genetic isolation of the Sundarbans tigers.     

Why do humans have chins? Testing the mechanical significance of modern human symphyseal morphology with finite element analysis

The modern human mandibular symphysis differs from those of all other primates in being vertically orientated and possessing a chin, but the functional significance of this unique morphology is not well understood. Some hypotheses propose that it is an adaptation to specific loads occurring during masticatory function. This study uses finite element analysis to examine these symphyseal loads in a model of a modern human mandible. By modifying the symphyseal cross-sectional form, the mechanical significance of the presence of the chin and symphyseal orientation is tested, and modern human and Neanderthal symphyseal cross-sections are compared with regard to their ability to withstand different loads. The results show that changes in symphyseal form have profound effects on the strains. The presence of a chin leads to lower symphyseal strains overall, whereas a vertical orientation of the symphysis results in higher strains under wishboning, but not under vertical bending in the coronal plane and dorsoventral shear. Compared to Neanderthals, the modern human symphysis shows higher strains during dorsoventral shear and wishboning, but is as effective as the Neanderthal symphysis in resisting vertical bending in the coronal plane and the loads resulting from simulated incision and unilateral molar biting. In general, the results of this study corroborate prior hypotheses about the mechanical effects of the human chin and vertical symphyseal orientation and support the idea that the relative importance of wishboning and vertical bending in the coronal plane might have played a role in the evolution of modern human symphyseal morphology.     

Water-stress-induced heat tolerance in geranium leaf tissues: A possible linkage through stress proteins?

Evidence is accumulating in favor of a linkage at the cellular level between various abiotic stresses. We conducted a study to evaluate the effect of water stress on the heat tolerance of zonal geraniums, Pelargonium × hortorum cv. Evening Glow. Water stress was imposed by withholding irrigation until pots reached 30% (by weight) of well‐watered controls, and by maintaining the pot weight by additions of water for another 7 days. Leaf xylem water potential (XWP, MPa), relative water content (RWC. %), and heat‐stress tolerance (HST; LT₅₀, defined as the temperature causing half‐maximal % injury based on electrolyte leakage) were measured in control, stressed, and recovered plants. Proteins were extracted from the leaves following the above treatments, and SDS‐PAGE and immunoblotting were performed by using standard procedures. Immunoblots were probed with antibodies to dehydrin and 70‐kDa heat shock cognate (HSC70) proteins. Data indicate that XWP and RWC, respectively, were −0.378 MPa and 92.3% for control plants and −0.804 MPa and 78.6% for stressed plants. Water‐stressed plants exhibited a significant increase in HST compared to control (LT₅₀ of 55°C vs 51°C). Water‐stress‐induced HST was not due to heat acclimation (leaf warming in stressed plants). Data also indicate that water‐stress treatment did not increase freezing tolerance of geranium leaves. Increased HST was associated with the accumulation of several heat‐stable, dehydrin proteins (25–60 kDa), and both cytosolic and ER luminal (BiP) HSC70 proteins. Leaf XWP, RWC, and HST reversed to control levels concomitant with the disappearance/reduction of dehydrins and HSC70 proteins in water‐stress‐relieved plants. The possibility of a cellular linkage between water stress and heat‐stress tolerance is discussed.     

Sensitivity analysis of the Poisson Nernst–Planck equations: a finite element approximation for the sensitive analysis of an electrodiffusion model

Biological structures exhibiting electric potential fluctuations such as neuron and neural structures with complex geometries are modelled using an electrodiffusion or Poisson Nernst–Planck system of equations. These structures typically depend upon several parameters displaying a large degree of variation or that cannot be precisely inferred experimentally. It is crucial to understand how the mathematical model (and resulting simulations) depend on specific values of these parameters. Here we develop a rigorous approach based on the sensitivity equation for the electrodiffusion model. To illustrate the proposed methodology, we investigate the sensitivity of the electrical response of a node of Ranvier with respect to ionic diffusion coefficients and the membrane dielectric permittivity.

     

Do gymnosperm needles pull water through the xylem produced in the same year as the needle?

This research investigated the longevity of functional connections between leaf traces and stem xylem in 16 species of conifers to better understand the spatial use of sapwood for water transport. The first question was which ring(s) stained when a vacuum was applied to the distal end of the cut surface of a needle attached to a short stem segment. The vacuum was applied to either 1- or 2-yr-old foliage taken from 4-6-yr-old saplings. The 16 species were then categorized based on the growth ring that most consistently stained when the vacuum was applied to the 2-yr-old needles. There were three distinct stain patterns for the evergreen conifers and one pattern for the deciduous conifers. Three evergreen species could not be categorized. The second question used needle fall data to ask whether the leaf trace appeared to break because of age or stem diameter in 5-yr-old Pseudotsuga menziesii saplings. An apparent threshold stem diameter at which needles tended to be shed was more related to diameter than age. These xylem connection patterns could affect leaf cohort physiology and the spatial pattern of water flux in sapwood.     

Three-dimensional finite element analysis of stress distribution in retention screws of different crown–implant ratios

The retaining screw of the implant-supported dental prosthesis is the weakest point of the crown/implant system. Furthermore, crown height is another important factor that may increase the lever arm. Therefore, the aim of this study was to assess the stress distribution in implant prosthetic screws with different heights of the clinical crown of the prosthesis using the method of three-dimensional finite element analysis. Three models were created with implants (3.75 mm × 10 mm) and crowns (heights of 10, 12.5 and 15 mm). The results were visualised by means of von Mises stress maps that increased the crown heights. The screw structure exhibited higher levels of stresses in the oblique load. The oblique loading resulted in higher stress concentration when compared with the axial loading. It is concluded that the increase of the crown was damaging to the stress distribution on the screw, mainly in oblique loading.