Cooperative motion of intrinsic and actuated semiflexible swimmers

We examine the phenomenon of hydrodynamic-induced cooperativity for pairs of flagellated micro-organism swimmers, of which spermatozoa cells are an example. We consider semiflexible swimmers, where inextensible filaments are driven by an internal intrinsic force and torque-free mechanism (intrinsic swimmers). The velocity gain for swimming cooperatively, which depends on both the geometry and the driving, develops as a result of the near-field coupling of bending and hydrodynamic stresses. We identify the regimes where hydrodynamic cooperativity is advantageous and quantify the change in efficiency. When the filaments' axes are parallel, hydrodynamic interaction induces a directional instability that causes semiflexible swimmers that profit from swimming together to move apart from each other. Biologically, this implies that flagella need to select different synchronized collective states and to compensate for directional instabilities (e.g., by binding) in order to profit from swimming together. By analyzing the cooperative motion of pairs of externally actuated filaments, we assess the impact that stress distribution along the filaments has on their collective displacements.

Soft tissue artifact distribution on lower limbs during treadmill gait: Influence of skin markers' location on cluster design.

Segment poses and joint kinematics estimated from skin markers are highly affected by soft tissue artifact (STA) and its rigid motion component (STARM). While four marker-clusters could decrease the STA non-rigid motion during gait activity, other data, such as marker location or STARM patterns, would be crucial to compensate for STA in clinical gait analysis. The present study proposed 1) to devise a comprehensive average map illustrating the spatial distribution of STA for the lower limb during treadmill gait and 2) to analyze STARM from four marker-clusters assigned to areas extracted from spatial distribution. All experiments were realized using a stereophotogrammetric system to track the skin markers and a bi-plane fluoroscopic system to track the knee prosthesis. Computation of the spatial distribution of STA was realized on 19 subjects using 80 markers apposed on the lower limb. Three different areas were extracted from the distribution map of the thigh. The marker displacement reached a maximum of 24.9mm and 15.3mm in the proximal areas of thigh and shank, respectively. STARM was larger on thigh than the shank with RMS error in cluster orientations between 1.2° and 8.1°. The translation RMS errors were also large (3.0mm to 16.2mm). No marker-cluster correctly compensated for STARM. However, the coefficient of multiple correlations exhibited excellent scores between skin and bone kinematics, as well as for STARM between subjects. These correlations highlight dependencies between STARM and the kinematic components. This study provides new insights for modeling STARM for gait activity.

Towards high-quality simultaneous EEG-fMRI at 7T: Detection and reduction of EEG artifacts due to head motion.

The enhanced functional sensitivity offered by ultra-high field imaging may significantly benefit simultaneous EEG-fMRI studies, but the concurrent increases in artifact contamination can strongly compromise EEG data quality. In the present study, we focus on EEG artifacts created by head motion in the static B0 field. A novel approach for motion artifact detection is proposed, based on a simple modification of a commercial EEG cap, in which four electrodes are non-permanently adapted to record only magnetic induction effects. Simultaneous EEG-fMRI data were acquired with this setup, at 7T, from healthy volunteers undergoing a reversing-checkerboard visual stimulation paradigm. Data analysis assisted by the motion sensors revealed that, after gradient artifact correction, EEG signal variance was largely dominated by pulse artifacts (81-93%), but contributions from spontaneous motion (4-13%) were still comparable to or even larger than those of actual neuronal activity (3-9%). Multiple approaches were tested to determine the most effective procedure for denoising EEG data incorporating motion sensor information. Optimal results were obtained by applying an initial pulse artifact correction step (AAS-based), followed by motion artifact correction (based on the motion sensors) and ICA denoising. On average, motion artifact correction (after AAS) yielded a 61% reduction in signal power and a 62% increase in VEP trial-by-trial consistency. Combined with ICA, these improvements rose to a 74% power reduction and an 86% increase in trial consistency. Overall, the improvements achieved were well appreciable at single-subject and single-trial levels, and set an encouraging quality mark for simultaneous EEG-fMRI at ultra-high field.

Intracranial Aneurysms: Wall Motion Analysis for Prediction of Rupture.

Intracranial aneurysms are a common pathologic condition with a potential severe complication: rupture. Effective treatment options exist, neurosurgical clipping and endovascular techniques, but guidelines for treatment are unclear and focus mainly on patient age, aneurysm size, and localization. New criteria to define the risk of rupture are needed to refine these guidelines. One potential candidate is aneurysm wall motion, known to be associated with rupture but difficult to detect and quantify. We review what is known about the association between aneurysm wall motion and rupture, which structural changes may explain wall motion patterns, and available imaging techniques able to analyze wall motion.

Effect of humeral stem design on humeral position and range of motion in reverse shoulder arthroplasty.

PURPOSE: The impacts of humeral offset and stem design after reverse shoulder arthroplasty (RSA) have not been well-studied, particularly with regard to newer stems which have a lower humeral inclination. The purpose of this study was to analyze the effect of different humeral stem designs on range of motion and humeral position following RSA. METHODS: Using a three-dimensional computer model of RSA, a traditional inlay Grammont stem was compared to a short curved onlay stem with different inclinations (155°, 145°, 135°) and offset (lateralised vs medialised). Humeral offset, the acromiohumeral distance (AHD), and range of motion were evaluated for each configuration. RESULTS: Altering stem design led to a nearly 7-mm change in humeral offset and 4 mm in the AHD. Different inclinations of the onlay stems had little influence on humeral offset and larger influence on decreasing the AHD. There was a 10° decrease in abduction and a 5° increase in adduction between an inlay Grammont design and an onlay design with the same inclination. Compared to the 155° model, the 135° model improved adduction by 28°, extension by 24° and external rotation of the elbow at the side by 15°, but led to a decrease in abduction of 9°. When the tray was placed medially, on the 145° model, a 9° loss of abduction was observed. CONCLUSIONS: With varus inclination prostheses (135° and 145°), elevation remains unchanged, abduction slightly decreases, but a dramatic improvement in adduction, extension and external rotation with the elbow at the side are observed.

On the computation of reducible invariant tori on a parallel computer

We present an algorithm for the computation of reducible invariant tori of discrete dynamical systems that is suitable for tori of dimensions larger than 1. It is based on a quadratically convergent scheme that approximates, at the same time, the Fourier series of the torus, its Floquet transformation, and its Floquet matrix. The Floquet matrix describes the linearization of the dynamics around the torus and, hence, its linear stability. The algorithm presents a high degree of parallelism, and the computational effort grows linearly with the number of Fourier modes needed to represent the solution. For these reasons it is a very good option to compute quasi-periodic solutions with several basic frequencies. The paper includes some examples (flows) to show the efficiency of the method in a parallel computer. In these flows we compute invariant tori of dimensions up to 5, by taking suitable sections.

Flexible Multibody Simulation Approach in the Dynamic Analysis of Bone Strains Activity

The objective of this study is to show that bone strains due to dynamic mechanical loading during physical activity can be analysed using the flexible multibody simulation approach. Strains within the bone tissue play a major role in bone (re)modeling. Based on previous studies, it has been shown that dynamic loading seems to be more important for bone (re)modeling than static loading. The finite element method has been used previously to assess bone strains. However, the finite element method may be limited to static analysis of bone strains due to the expensive computation required for dynamic analysis, especially for a biomechanical system consisting of several bodies. Further, in vivo implementation of strain gauges on the surfaces of bone has been used previously in order to quantify the mechanical loading environment of the skeleton. However, in vivo strain measurement requires invasive methodology, which is challenging and limited to certain regions of superficial bones only, such as the anterior surface of the tibia. In this study, an alternative numerical approach to analyzing in vivo strains, based on the flexible multibody simulation approach, is proposed. In order to investigate the reliability of the proposed approach, three 3-dimensional musculoskeletal models where the right tibia is assumed to be flexible, are used as demonstration examples. The models are employed in a forward dynamics simulation in order to predict the tibial strains during walking on a level exercise. The flexible tibial model is developed using the actual geometry of the subject’s tibia, which is obtained from 3 dimensional reconstruction of Magnetic Resonance Images. Inverse dynamics simulation based on motion capture data obtained from walking at a constant velocity is used to calculate the desired contraction trajectory for each muscle. In the forward dynamics simulation, a proportional derivative servo controller is used to calculate each muscle force required to reproduce the motion, based on the desired muscle contraction trajectory obtained from the inverse dynamics simulation. Experimental measurements are used to verify the models and check the accuracy of the models in replicating the realistic mechanical loading environment measured from the walking test. The predicted strain results by the models show consistency with literature-based in vivo strain measurements. In conclusion, the non-invasive flexible multibody simulation approach may be used as a surrogate for experimental bone strain measurement, and thus be of use in detailed strain estimation of bones in different applications. Consequently, the information obtained from the present approach might be useful in clinical applications, including optimizing implant design and devising exercises to prevent bone fragility, accelerate fracture healing and reduce osteoporotic bone loss.

Application of intravoxel incoherent motion perfusion imaging to shoulder muscles after a lift-off test of varying duration.

Intravoxel incoherent motion (IVIM) MRI is a method to extract microvascular blood flow information out of diffusion-weighted images acquired at multiple b-values. We hypothesized that IVIM can identify the muscles selectively involved in a specific task, by measuring changes in activity-induced local muscular perfusion after exercise. We tested this hypothesis using a widely used clinical maneuver, the lift-off test, which is known to assess specifically the subscapularis muscle functional integrity. Twelve shoulders from six healthy male volunteers were imaged at 3 T, at rest, as well as after a lift-off test hold against resistance for 30 s, 1 and 2 min respectively, in three independent sessions. IVIM parameters, consisting of perfusion fraction (f), diffusion coefficient (D), pseudo-diffusion coefficient D* and blood flow-related fD*, were estimated within outlined muscles of the rotator cuff and the deltoid bundles. The mean values at rest and after the lift-off tests were compared in each muscle using a one-way ANOVA. A statistically significant increase in fD* was measured in the subscapularis, after a lift-off test of any duration, as well as in D. A fD* increase was the most marked (30 s, +103%; 1 min, +130%; 2 min, +156%) and was gradual with the duration of the test (in 10(-3) mm(2) /s: rest, 1.41 ± 0.50; 30 s, 2.86 ± 1.17; 1 min, 3.23 ± 1.22; 2 min, 3.60 ± 1.21). A significant increase in fD* and D was also visible in the posterior bundle of the deltoid. No significant change was consistently visible in the other investigated muscles of the rotator cuff and the other bundles of the deltoid. In conclusion, IVIM fD* allows the demonstration of a task-related microvascular perfusion increase after a specific task and suggests a direct relationship between microvascular perfusion and the duration of the effort. It is a promising method to investigate non-invasively skeletal muscle physiology and clinical perfusion-related muscular disorders.

Fast and Rigorous Computation of Gene and Pathway Scores from SNP-Based Summary Statistics.

Integrating single nucleotide polymorphism (SNP) p-values from genome-wide association studies (GWAS) across genes and pathways is a strategy to improve statistical power and gain biological insight. Here, we present Pascal (Pathway scoring algorithm), a powerful tool for computing gene and pathway scores from SNP-phenotype association summary statistics. For gene score computation, we implemented analytic and efficient numerical solutions to calculate test statistics. We examined in particular the sum and the maximum of chi-squared statistics, which measure the strongest and the average association signals per gene, respectively. For pathway scoring, we use a modified Fisher method, which offers not only significant power improvement over more traditional enrichment strategies, but also eliminates the problem of arbitrary threshold selection inherent in any binary membership based pathway enrichment approach. We demonstrate the marked increase in power by analyzing summary statistics from dozens of large meta-studies for various traits. Our extensive testing indicates that our method not only excels in rigorous type I error control, but also results in more biologically meaningful discoveries.

A minimal set of coordinates for describing humanoid shoulder motion

The kinematics of the anatomical shoulder are analysed and modelled as a parallel mechanism similar to a Stewart platform. A new method is proposed to describe the shoulder kinematics with minimal coordinates and solve the indeterminacy. The minimal coordinates are defined from bony landmarks and the scapulothoracic kinematic constraints. Independent from one another, they uniquely characterise the shoulder motion. A humanoid mechanism is then proposed with identical kinematic properties. It is then shown how minimal coordinates can be obtained for this mechanism and how the coordinates simplify both the motion-planning task and trajectory-tracking control. Lastly, the coordinates are also shown to have an application in the field of biomechanics where they can be used to model the scapulohumeral rhythm.

Prospective motion correction with FID-triggered image navigators

In this work, we propose a method for prospective motion correction in MRI using a novel image navigator module, which is triggered by a free induction decay (FID) navigator. Only when motion occurs, the image navigator is run and new positional information is obtained through image registration. The image navigator was specifically designed to match the impact on the magnetization and the acoustic noise of the host sequence. This detection-correction scheme was implemented for an MP-RAGE sequence and 5 healthy volunteers were scanned at 3T while performing various head movements. The correction performance was demonstrated through automated brain segmentation and an image quality index whose results are sensitive to motion artifacts.

Apnea-like suppression of respiratory motion: First evaluation in radiotherapy.

BACKGROUND AND PURPOSE: Compensation for respiratory motion is needed while administering radiotherapy (RT) to tumors that are moving with respiration to reduce the amount of irradiated normal tissues and potentially decrease radiation-induced collateral damages. The purpose of this study was to test a new ventilation system designed to induce apnea-like suppression of respiratory motion and allow long enough breath hold durations to deliver complex RT. MATERIAL AND METHODS: The High Frequency Percussive Ventilation system was initially tested in a series of 10 volunteers and found to be well tolerated, allowing a median breath hold duration of 11.6min (range 3.9-16.5min). An evaluation of this system was subsequently performed in 4 patients eligible for adjuvant breast 3D conformal RT, for lung stereotactic body RT (SBRT), lung volumetric modulated arc therapy (VMAT), and VMAT for palliative pleural metastases. RESULTS: When compared to free breathing (FB) and maximal inspiration (MI) gating, this Percussion Assisted RT (PART) offered favorable dose distribution profiles in 3 out of the 4 patients tested. PART was applied in these 3 patients with good tolerance, without breaks during the "beam on time period" throughout the overall courses of RT. The mean duration of the apnea-like breath hold that was necessary for delivering all the RT fractions was 7.61min (SD=2.3). CONCLUSIONS: This first clinical implementation of PART was found to be feasible, tolerable and offers new opportunities in the field of RT for suppressing respiratory motion.

Monte Carlo derivation of filtered tungsten anode X-ray spectra for dose computation in digital mammography

Abstract Objective: Derive filtered tungsten X-ray spectra used in digital mammography systems by means of Monte Carlo simulations. Materials and Methods: Filtered spectra for rhodium filter were obtained for tube potentials between 26 and 32 kV. The half-value layer (HVL) of simulated filtered spectra were compared with those obtained experimentally with a solid state detector Unfors model 8202031-H Xi R/F & MAM Detector Platinum and 8201023-C Xi Base unit Platinum Plus w mAs in a Hologic Selenia Dimensions system using a direct radiography mode. Results: Calculated HVL values showed good agreement as compared with those obtained experimentally. The greatest relative difference between the Monte Carlo calculated HVL values and experimental HVL values was 4%. Conclusion: The results show that the filtered tungsten anode X-ray spectra and the EGSnrc Monte Carlo code can be used for mean glandular dose determination in mammography.