Comparison of eight published static finite element models of the intact lumbar spine : predictive power of models improves when combined together

Finite element (FE) model studies have made important contributions to our understanding of functional biomechanics of the lumbar spine. However, if a model is used to answer clinical and biomechanical questions over a certain population, their inherently large inter-subject variability has to be considered. Current FE model studies, however, generally account only for a single distinct spinal geometry with one set of material properties. This raises questions concerning their predictive power, their range of results and on their agreement with in vitro and in vivo values. Eight well-established FE models of the lumbar spine (L1-5) of different research centres around the globe were subjected to pure and combined loading modes and compared to in vitro and in vivo measurements for intervertebral rotations, disc pressures and facet joint forces. Under pure moment loading, the predicted L1-5 rotations of almost all models fell within the reported in vitro ranges, and their median values differed on average by only 2° for flexion-extension, 1° for lateral bending and 5° for axial rotation. Predicted median facet joint forces and disc pressures were also in good agreement with published median in vitro values. However, the ranges of predictions were larger and exceeded those reported in vitro, especially for the facet joint forces. For all combined loading modes, except for flexion, predicted median segmental intervertebral rotations and disc pressures were in good agreement with measured in vivo values. In light of high inter-subject variability, the generalization of results of a single model to a population remains a concern. This study demonstrated that the pooled median of individual model results, similar to a probabilistic approach, can be used as an improved predictive tool in order to estimate the response of the lumbar spine.

The cytochrome P-450 isoenzyme CYP2E1 in the biological processing of industrial chemicals : consequences for occupational and environmental medicine

The importance of the isoform CYP2E1 of the human cytochrome P-450 superfamily of enzymes for occupational and environmental medicine is derived from its unique substrate spectrum that includes a number of highly important high-production chemicals, such as aliphatic and aromatic hydrocarbons, solvents and industrial monomers (i.a. alkanes, alkenes, aromatic and halogenated hydrocarbons). Many polymorphic genes, such as CYP2E1, show considerable differences in allelic distribution between different human populations. The polymorphic nature of the human CYP2E1 gene is significant for inter-individual differences in toxicity of its substrates. Since the substrate spectrum of CYP2E1 includes many compounds of basic relevance to industrial toxicology, a rationale for metabolic interactions of different CYP2E1 substrates is provided. In-depth research into the inter-individual phenotypic differences of human CYP2E1 enzyme activities was enabled by the recognition that the 6-hydroxylation of the drug chlorzoxazone is mediated by CYP2E1. Studies on CYP2E1 phenotyping have pointed to inter-individual variations in enzyme activities. There are consistent ethnic differences in CYP2E1 enzyme expression, mostly demonstrated between European and Japanese populations, which point to a major impact of genetic factors. The most frequently studied genetic polymorphisms are the restriction fragment length polymorphisms PstI/RsaI (mutant allele: CYP2E1*5B) located in the 5′-flanking region of the gene, as well as the DraI polymorphism (mutant allele: CYP2E1*6) located in intron 6. These polymorphisms are partly related, as they form the common allele designated CYP2E1*5A. Striking inter-ethnic differences between Europeans and Asians appear with respect to the frequencies of the CYP2E1*5A allele (only approximately 5% of Europeans are heterozygous, but 37% of Asians are, whilst 6% of Asians are homozygous). Available studies indicate a wide variation in human CYP2E1 expression, which are very likely based on complex gene-environment interactions. Major inter-ethnic differences are apparent on the genotyping and the phenotyping levels. Selected cases are presented where inter-ethnic variations of CYP2E1 may provide likely explanations for unexplained findings concerning industrial chemicals that are CYP2E1 substrates. Possible consequences of differential inter-individual and inter-ethnic susceptibilities are related to individual expressions of clinical symptoms of chemical toxicity, to results of biological monitoring of exposed workers, and to the interpretation of results of epidemiological or molecular-epidemiological studies.

Detoxification enzyme activities (CYP1A1 and GST) in the skin of humpback whales as a function of organochlorine burdens and migration status

The activities of glutathione-s-transferase (GST) and cytochrome P-450 1A1 (CYP1A1) enzymes were measured in freshly extracted epidermis of live-biopsied, migrating, southern hemisphere humpback whales (Megaptera novaeangliae). The two quantified enzyme activities did not correlate strongly with each other. Similarly, neither correlated strongly with any of the organochlorine compound groups previously measured in the superficial blubber of the sample biopsy core, likely reflecting the anticipated low levels of typical aryl-hydrocarbon receptor ligands. GST activity did not differ significantly between genders or between northward (early migration) or southward (late migration) migrating cohorts. Indeed, the inter-individual variability in GST measurements was relatively low. This observation raises the possibility that measured activities were basal activities and that GST function was inherently impacted by the fasting state of the sampled animals, as seen in other species. These results do not support the implementation of CYP1A1 or GST as effective biomarkers of organochlorine contaminant burdens in southern hemisphere populations of humpback whales as advocated for other cetacean species. Further investigation of GST activity in feeding versus fasting cohorts may, however, provide some insight into the fasting metabolism of these behaviourally adapted populations.

Comparison of eight published lumbar spine finite element models

Due to its ability to represent intricate systems with material nonlinearities as well as irregular loading, boundary, geometrical and material domains, the finite element (FE) method has been recognized as an important computational tool in spinal biomechanics. Current FE models generally account for a single distinct spinal geometry with one set of material properties despite inherently large inter-subject variability. The uncertainty and high variability in tissue material properties, geometry, loading and boundary conditions has cast doubt on the reliability of their predictions and comparability with reported in vitro and in vivo values. A multicenter study was undertaken to compare the results of eight well-established models of the lumbar spine that have been developed, validated and applied for many years. Models were subjected to pure and combined loading modes and their predictions were compared to in vitro and in vivo measurements for intervertebral rotations, disc pressures and facet joint forces. Under pure moment loading, the predicted L1-5 rotations of almost all models fell within the reported in vitro ranges; their median values differed on average by only 2° for flexion-extension, 1° for lateral bending and 5° for axial rotation. Predicted median facet joint forces and disc pressures were also in good agreement with previously published median in vitro values. However, the ranges of predictions were larger and exceeded the in vitro ranges, especially for facet joint forces. For all combined loading modes, except for flexion, predicted median segmental intervertebral rotations and disc pressures were in good agreement with in vivo values. The simulations yielded median facet joint forces of 0 N in flexion, 38 N in extension, 14 N in lateral bending and 60 N in axial rotation that could not be validated due to the paucity of in vivo facet joint forces. In light of high inter-subject variability, one must be cautious when generalizing predictions obtained from one deterministic model. This study demonstrates however that the predictive power increases when FE models are combined together. The median of individual numerical results can hence be used as an improved tool in order to estimate the response of the lumbar spine.

Bi-modal biometric authentication on mobile phones in challenging conditions

This paper examines the issue of face, speaker and bi-modal authentication in mobile environments when there is significant condition mismatch. We introduce this mismatch by enrolling client models on high quality biometric samples obtained on a laptop computer and authenticating them on lower quality biometric samples acquired with a mobile phone. To perform these experiments we develop three novel authentication protocols for the large publicly available MOBIO database. We evaluate state-of-the-art face, speaker and bi-modal authentication techniques and show that inter-session variability modelling using Gaussian mixture models provides a consistently robust system for face, speaker and bi-modal authentication. It is also shown that multi-algorithm fusion provides a consistent performance improvement for face, speaker and bi-modal authentication. Using this bi-modal multi-algorithm system we derive a state-of-the-art authentication system that obtains a half total error rate of 6.3% and 1.9% for Female and Male trials, respectively.

Robust face clustering for real-world data

This thesis has investigated how to cluster a large number of faces within a multi-media corpus in the presence of large session variation. Quality metrics are used to select the best faces to represent a sequence of faces; and session variation modelling improves clustering performance in the presence of wide variations across videos. Findings from this thesis contribute to improving the performance of both face verification systems and the fully automated clustering of faces from a large video corpus.

Persistent organic pollutants in matched breast milk and infant faeces samples

Assessing blood concentration of persistent organic pollutants (POPs) in infants is difficult due to the ethical and practical difficulties in obtaining sufficient quantities of blood. To determine whether measuring POPs in faeces might reflect blood concentration during infancy, we measured the concentrations of a range of POPs (i.e. polychlorinated biphenyls (PCBs), polybrominated diphenyl ethers (PBDEs) and organochlorine pesticides (OCPs)) in a pilot study using matched breast milk and infant faecal samples obtained from ten mother-child pairs. All infants were breast fed, with 8 of them also receiving solid food at the time of faecal sampling. In this small dataset faecal concentrations (range 0.01-41ngg-1 lipid) are strongly associated with milk concentrations (range 0.02-230ngg-1 lipid). Associations with other factors generally could not be detected in this dataset, with the exception of a small effect of age or growth. Different sources (external or internal) of exposure appeared to directly influence faecal concentrations of different chemicals based on different inter-individual variability in the faeces-to-milk concentration ratio Rfm. Overall, the matrix of faeces as an external measure of internal exposure in infants looks promising for some chemicals and is worth assessing further in larger datasets.

A framework for fully decentralised cycle stealing

Ordinary desktop computers continue to obtain ever more resources – in-creased processing power, memory, network speed and bandwidth – yet these resources spend much of their time underutilised. Cycle stealing frameworks harness these resources so they can be used for high-performance computing. Traditionally cycle stealing systems have used client-server based architectures which place significant limits on their ability to scale and the range of applica-tions they can support. By applying a fully decentralised network model to cycle stealing the limits of centralised models can be overcome. Using decentralised networks in this manner presents some difficulties which have not been encountered in their previous uses. Generally decentralised ap-plications do not require any significant fault tolerance guarantees. High-performance computing on the other hand requires very stringent guarantees to ensure correct results are obtained. Unfortunately mechanisms developed for traditional high-performance computing cannot be simply translated because of their reliance on a reliable storage mechanism. In the highly dynamic world of P2P computing this reliable storage is not available. As part of this research a fault tolerance system has been created which provides considerable reliability without the need for a persistent storage. As well as increased scalability, fully decentralised networks offer the ability for volunteers to communicate directly. This ability provides the possibility of supporting applications whose tasks require direct, message passing style communication. Previous cycle stealing systems have only supported embarrassingly parallel applications and applications with limited forms of communication so a new programming model has been developed which can support this style of communication within a cycle stealing context. In this thesis I present a fully decentralised cycle stealing framework. The framework addresses the problems of providing a reliable fault tolerance sys-tem and supporting direct communication between parallel tasks. The thesis includes a programming model for developing cycle stealing applications with direct inter-process communication and methods for optimising object locality on decentralised networks.

Whole of life cycle cost analysis in bridge rehabilitation

The report presents a methodology for whole of life cycle cost analysis of alternative treatment options for bridge structures, which require rehabilitation. The methodology has been developed after a review of current methods and establishing that a life cycle analysis based on a probabilistic risk approach has many advantages including the essential ability to consider variability of input parameters. The input parameters for the analysis are identified as initial cost, maintenance, monitoring and repair cost, user cost and failure cost. The methodology utilizes the advanced simulation technique of Monte Carlo simulation to combine a number of probability distributions to establish the distribution of whole of life cycle cost. In performing the simulation, the need for a powerful software package, which would work with spreadsheet program, has been identified. After exploring several products on the market, @RISK software has been selected for the simulation. In conclusion, the report presents a typical decision making scenario considering two alternative treatment options.

Variability, regularity and coupling measures distinguish PD tremor from voluntary 5Hz tremor

A characteristic of Parkinson's disease (PD) is the development of tremor within the 4–6 Hz range. One method used to better understand pathological tremor is to compare the responses to tremor-type actions generated intentionally in healthy adults. This study was designed to investigate the similarities and differences between voluntarily generated 4–6 Hz tremor and PD tremor in regards to their amplitude, frequency and coupling characteristics. Tremor responses for 8 PD individuals (on- and off-medication) and 12 healthy adults were assessed under postural and resting conditions. Results showed that the voluntary and PD tremor were essentially identical with regards to the amplitude and peak frequency. However, differences between the groups were found for the variability (SD of peak frequency, proportional power) and regularity (Approximate Entropy, ApEn) of the tremor signal. Additionally, coherence analysis revealed strong inter-limb coupling during voluntary conditions while no bilateral coupling was seen for the PD persons. Overall, healthy participants were able to produce a 5 Hz tremulous motion indistinguishable to that of PD patients in terms of peak frequency and amplitude. However, differences in the structure of variability and level of inter-limb coupling were found for the tremor responses of the PD and healthy adults. These differences were preserved irrespective of the medication state of the PD persons. The results illustrate the importance of assessing the pattern of signal structure/variability to discriminate between different tremor forms, especially where no differences emerge in standard measures of mean amplitude as traditionally defined.

Inter-vertebral rotational deformity after endoscopic anterior scoliosis correction may contribute to rib hump recurrence after two years

Introduction. Endoscopic anterior scoliosis correction has been employed recently as a less invasive and level-sparing approach compared with open surgical techniques. We have previously demonstrated that during the two-year post-operative period, there was a mean loss of rib hump correction by 1.4 degrees. The purpose of this study was to determine whether intra- or inter-vertebral rotational deformity during the post-operative period could account for the loss of rib hump correction. Materials and Methods. Ten consecutive patients diagnosed with adolescent idiopathic scoliosis were treated with an endoscopic anterior scoliosis correction. Low-dose computed tomography scans of the instrumented segment were obtained post-operatively at 6 and 24 months following institutional ethical approval and patient consent. Three-dimensional multi-planar reconstruction software (Osirix Imaging Software, Pixmeo, Switzerland) was used to create axial slices of each vertebral level, corrected in both coronal and sagittal planes. Vertebral rotation was measured using Ho’s method for every available superior and inferior endplate at 6 and 24 months. Positive changes in rotation indicate a reduction and improvement in vertebral rotation. Intra-observer variability analysis was performed on a subgroup of images. Results. Mean change in rotation for vertebral endplates between 6 and 24 months post-operatively was -0.26˚ (range -3.5 to 4.9˚) within the fused segment and +1.26˚ (range -7.2 to 15.1˚) for the un-instrumented vertebrae above and below the fusion. Mean change in clinically measured rib hump for the 10 patients was -1.6˚ (range -3 to 0˚). The small change in rotation within the fused segment accounts for only 16.5% of the change in rib hump measured clinically whereas the change in rotation between the un-instrumented vertebrae above and below the construct accounts for 78.8%. There was no clear association between rib hump recurrence and intra- or inter-vertebral rotation in individual patients. Intra-rater variability was ± 3˚. Conclusions. Intra- and inter-vertebral rotation continues post-operatively both within the instrumented and un-instrumented segments of the immature spine. Rotation between the un-instrumented vertebrae above and below the fusion was +1.26˚, suggesting that the un-instrumented vertebrae improved and de-rotated slightly after surgery. This may play a role in rib hump recurrence, however this remains clinically insignificant.

Investigations on road noise level spatial variability within a specially designed acoustic balcony

An investigation into the spatial distribution of road traffic noise levels on a balcony is conducted. A balcony constructed to a special acoustic design due to its elevation above an 8 lane motorway is selected for detailed measurements. The as-constructed balcony design includes solid parapets, side walls, ceiling shields and highly absorptive material placed on the ceiling. Road traffic noise measurements are conducted spatially using a five channel acoustic analyzer, where four microphones are located at various positions within the balcony space and one microphone placed outside the parapet at a reference position. Spatial distributions in both vertical and horizontal planes are measured. A theoretical model and prediction configuration is presented that assesses the acoustic performance of the balcony under existing traffic flow conditions. The prediction model implements a combined direct path, specular reflection path and diffuse reflection path utilizing image source and radiosity techniques. Results obtained from the prediction model are presented and compared to the measurement results. The predictions are found to correlate well with measurements with some minor differences that are explained. It is determined that the prediction methodology is acceptable to assess a wider range of street and balcony configuration scenarios.

Seasonal and spatial variations in rare earth elements to identify inter-aquifer linkages and recharge processes in an Australian catchment

With the aim of elucidating the seasonal behaviour of rare earth elements (REEs), surface and groundwaters were collected under dry and wet conditions in different hydrological units of the Teviot Brook catchment (Southeast Queensland, Australia). Sampled waters showed a large degree of variability in both REE abundance and normalised patterns. Overall REE abundance ranged over nearly three orders of magnitude, and was consistently lower in the sedimentary bedrock aquifer (18ppt<∑REE<477ppt) than in the other hydrological systems studied. Abundance was greater in springs draining rhyolitic rocks (∑REE=300 and 2054ppt) than in springs draining basalt ranges (∑REE=25 and 83ppt), yet was highly variable in the shallow alluvial groundwater (16ppt<∑REE<5294ppt) and, to a lesser extent, in streamwater (85ppt<∑REE<2198ppt). Generally, waters that interacted with different rock types had different REE patterns. In order to obtain an unbiased characterisation of REE patterns, the ratios between light and middle REEs (R(M/L)) and the ratios between middle and heavy REEs (R(H/M)) were calculated for each sample. The sedimentary bedrock aquifer waters had highly evolved patterns depleted in light REEs and enriched in middle and heavy REEs (0.17

Population of computational rabbit-specific ventricular action potential models for investigating sources of variability in cellular repolarisation

Variability is observed at all levels of cardiac electrophysiology. Yet, the underlying causes and importance of this variability are generally unknown, and difficult to investigate with current experimental techniques. The aim of the present study was to generate populations of computational ventricular action potential models that reproduce experimentally observed intercellular variability of repolarisation (represented by action potential duration) and to identify its potential causes. A systematic exploration of the effects of simultaneously varying the magnitude of six transmembrane current conductances (transient outward, rapid and slow delayed rectifier K(+), inward rectifying K(+), L-type Ca(2+), and Na(+)/K(+) pump currents) in two rabbit-specific ventricular action potential models (Shannon et al. and Mahajan et al.) at multiple cycle lengths (400, 600, 1,000 ms) was performed. This was accomplished with distributed computing software specialised for multi-dimensional parameter sweeps and grid execution. An initial population of 15,625 parameter sets was generated for both models at each cycle length. Action potential durations of these populations were compared to experimentally derived ranges for rabbit ventricular myocytes. 1,352 parameter sets for the Shannon model and 779 parameter sets for the Mahajan model yielded action potential duration within the experimental range, demonstrating that a wide array of ionic conductance values can be used to simulate a physiological rabbit ventricular action potential. Furthermore, by using clutter-based dimension reordering, a technique that allows visualisation of multi-dimensional spaces in two dimensions, the interaction of current conductances and their relative importance to the ventricular action potential at different cycle lengths were revealed. Overall, this work represents an important step towards a better understanding of the role that variability in current conductances may play in experimentally observed intercellular variability of rabbit ventricular action potential repolarisation.