Using patient specific finite element models to understand the biomechanics of paediatric spinal deformity surgery

Adolescent idiopathic scoliosis (AIS) is a spinal deformity, which may require surgical correction by attaching rods to the patient’s spine using screws inserted into the vertebrae. Complication rates for deformity correction surgery are unacceptably high. Determining an achievable correction without overloading the adjacent spinal tissues or implants requires an understanding of the mechanical interaction between these components. We have developed novel patient specific modelling software to create individualized finite element models (FEM) representing the thoracolumbar spine and ribcage of scoliosis patients. We are using these models to better understand the biomechanics of spinal deformity correction.

The topography of multiple drug use among adolescent Australians : findings from the National Drug Strategy Household Survey

Introduction and aims: Despite evidence that many Australian adolescents have considerable experience with various drug types, little is known about the extent to which adolescents use multiple substances. The aim of this study was to examine the degree of clustering of drug types within individuals, and the extent to which demographic and psychosocial predictors are related to cluster membership. Design and method: A sample of 1402 adolescents aged 12-17. years were extracted from the Australian 2007 National Drug Strategy Household Survey. Extracted data included lifetime use of 10 substances, gender, psychological distress, physical health, perceived peer substance use, socioeconomic disadvantage, and regionality. Latent class analysis was used to determine clusters, and multinomial logistic regression employed to examine predictors of cluster membership. Result: There were 3 latent classes. The great majority (79.6%) of adolescents used alcohol only, 18.3% were limited range multidrug users (encompassing alcohol, tobacco, and marijuana), and 2% were extended range multidrug users. Perceived peer drug use and psychological distress predicted limited and extended multiple drug use. Psychological distress was a more significant predictor of extended multidrug use compared to limited multidrug use. Discussion and conclusion: In the Australian school-based prevention setting, a very strong focus on alcohol use and the linkages between alcohol, tobacco and marijuana are warranted. Psychological distress may be an important target for screening and early intervention for adolescents who use multiple drugs.

Phage library screening for the rapid identification and in vivo testing of candidate genes for a DNA vaccine against Mycoplasma mycoides subsp. mycoides small colony biotype

A new strategy for rapidly selecting and testing genetic vaccines has been developed, in which a whole genome library is cloned into a bacteriophage λ ZAP Express vector which contains both prokaryotic (Plac) and eukaryotic (PCMV) promoters upstream of the insertion site. The phage library is plated on Escherichia coli cells, immunoblotted, and probed with hyperimmune and/or convalescent-phase antiserum to rapidly identify vaccine candidates. These are then plaque purified and grown as liquid lysates, and whole bacteriophage particles are then used directly to immunize the host, following which PCMV-driven expression of the candidate vaccine gene occurs. In the example given here, a semirandom genome library of the bovine pathogen Mycoplasma mycoides subsp. mycoides small colony (SC) biotype was cloned into λ ZAP Express, and two strongly immunodominant clones, λ-A8 and λ-B1, were identified and subsequently tested for vaccine potential against M. mycoides subsp. mycoides SC biotype-induced mycoplasmemia. Sequencing and immunoblotting indicated that clone λ-A8 expressed an isopropyl-β-d-thiogalactopyranoside (IPTG)-inducible M. mycoides subsp. mycoides SC biotype protein with a 28-kDa apparent molecular mass, identified as a previously uncharacterized putative lipoprotein (MSC_0397). Clone λ-B1 contained several full-length genes from the M. mycoides subsp. mycoides SC biotype pyruvate dehydrogenase region, and two IPTG-independent polypeptides, of 29 kDa and 57 kDa, were identified on immunoblots. Following vaccination, significant anti-M. mycoides subsp. mycoides SC biotype responses were observed in mice vaccinated with clones λ-A8 and λ-B1. A significant stimulation index was observed following incubation of splenocytes from mice vaccinated with clone λ-A8 with whole live M. mycoides subsp. mycoides SC biotype cells, indicating cellular proliferation. After challenge, mice vaccinated with clone λ-A8 also exhibited a reduced level of mycoplasmemia compared to controls, suggesting that the MSC_0397 lipoprotein has a protective effect in the mouse model when delivered as a bacteriophage DNA vaccine. Bacteriophage-mediated immunoscreening using an appropriate vector system offers a rapid and simple technique for the identification and immediate testing of putative candidate vaccines from a variety of pathogens.

Conjugation of carcinogens by 6 class glutathione S-transferases : mechanisms and relevance to variations in human risk

Conjugation of chemicals with glutathione (GSH) can lead to decreased or increased toxicity. A genetic deficiency in the GSH S-transferase μ class gene M1 has been hypothesized to lead to greater risk of lung cancer in smokers. Recently a gene deletion polymorphism involving the human θ enzyme T1 has been described; the enzyme is present in erythrocytes and can be readily assayed. A rat θ class enzyme, 5-5, has structural and catalytic similarity and the protein was expressed in the Salmonella typhimurium tester strain TA1535. Expression of the cDNA vector increased the mutagenicity of ethylene dibromide and several methylene dihalides. Mutations resulting from the known GSH S-transferase substrate 1,2-epoxy-3-(4′nitrophenoxy)propane were decreased in the presence of the transferase. Expression of transferase 5-5 increased mutations when 1,2,3,4-diepoxybutane (butadiene diepoxide), 4-bromo-1,2-epoxybutane, or 1,3-dichloracetone were added. The latter compound is a model for the putative 1,2-dibromo-3-chloropropane oxidation product 1-bromo-3-chloroacetone. These genotoxicity and genotyping assays may be of use in further studies of the roles of GSH S-transferase θ enzymes in bioactivation and detoxication and any changes in risk due to polymorphism.

Enhancement of bacterial mutagenicity of bifunctional alkylating agents by expression of mammalian glutathione s-transferase

Recently, we inserted the plasmid vector pKK233-2 containing rat GSH S-transferase (GST) 5-5 cDNA into Salmonella typhimurium TA1535 and found that these bacteria [GST 5-5(+)] expressed the protein and produced mutations when ethylene or methylene dihalides were added [Thier, R., Taylor, J. B., Pemble, S. E., Ketterer, B., Persmark, M., Humphreys, W. G., and Guengerich, F. P. (1993) Proc. Natl. Acad. Sci. U.S.A. 90, 8576-8580]. After exposure to the known GST 5-5 substrate 1,2-epoxy-3-(4′-nitrophenoxy)propane, the GST 5-5(+) strain showed fewer mutants than the bacteria transfected with the cDNA clone in a reverse orientation [GST 5-5(-)], suggesting a protective role of GST 5-5. However, mutations were considerably enhanced in the GST 5-5(+) strain [as compared to GST 5-5(-)] when 1,2,3,4-diepoxybutane (butadiene diepoxide) or 1,2-epoxy-4-bromobutane was added. The GST 5-5(+) and GST 5-5(-) bacterial stains showed similar responses to 1,2-epoxypropane, 3,4-epoxy-1-butene, and 1,4-dibromobutane. The results suggest that some bifunctional activated butanes are transformed to mutagenic products through GSH conjugation. We also found that the GST 5-5(+) strain showed enhanced mutagenicity with 1,4-dibromo-2,3-epoxybutane, 1,2-epoxy-3-bromopropane (epibromohydrin), and (±)-1,4-dibromo-2,3-dihydroxybutane. The possibility was considered that a 5-membered thialonium ion may be involved in the mutagenicity. Model thialonium compounds were rather stable to hydrolysis in aqueous solution at pH 7.4 and slowly alkylated 4-(4-nitrobenzyl)pyridine. The presence of a hydroxyl group β to the sulfur did not enhance reactivity. Mechanisms involving episulfonium ions are considered more likely. Potential oxidation products of the toxic pesticide 1,2-dibromo-3-chloropropane (DBCP) were also considered in this system. DBCP itself gave rather similar results in the two strains. Others have reported that oxidation of DBCP is required for mutagenicity, along with GST-catalyzed GSH conjugation [Simula, T. P., Glancey, M. J., Söderlund, E. J., Dybing, E., and Wolf, C. R. (1993) Carcinogenesis 14, 2303-2307]. The putative oxidation product 1,2-dibromopropional did not show a difference between the two strains. However, 1,3-dichloroacetone, a model for the putative oxidation product 1-bromo-3-chloroacetone, was considerably more mutagenic in the GST 5-5(+) strain.

Adaptive estimation of HMM transition probabilities

This paper presents new schemes for recursive estimation of the state transition probabilities for hidden Markov models (HMM's) via extended least squares (ELS) and recursive state prediction error (RSPE) methods. Local convergence analysis for the proposed RSPE algorithm is shown using the ordinary differential equation (ODE) approach developed for the more familiar recursive output prediction error (RPE) methods. The presented scheme converges and is relatively well conditioned compared with the ...

On adaptive HMM state estimation

In this paper new online adaptive hidden Markov model (HMM) state estimation schemes are developed, based on extended least squares (ELS) concepts and recursive prediction error (RPE) methods. The best of the new schemes exploit the idempotent nature of Markov chains and work with a least squares prediction error index, using a posterior estimates, more suited to Markov models then traditionally used in identification of linear systems.

Adaptive estimation of hidden semi-Markov chains with parameterised transition probabilities and exponential decaying states

This paper develops maximum likelihood (ML) estimation schemes for finite-state semi-Markov chains in white Gaussian noise. We assume that the semi-Markov chain is characterised by transition probabilities of known parametric from with unknown parameters. We reformulate this hidden semi-Markov model (HSM) problem in the scalar case as a two-vector homogeneous hidden Markov model (HMM) problem in which the state consist of the signal augmented by the time to last transition. With this reformulation we apply the expectation Maximumisation (EM ) algorithm to obtain ML estimates of the transition probabilities parameters, Markov state levels and noise variance. To demonstrate our proposed schemes, motivated by neuro-biological applications, we use a damped sinusoidal parameterised function for the transition probabilities.

Optimal HMM processing of differentially phase modulated communication systems

This paper investigates demodulation of differentially phase modulated signals DPMS using optimal HMM filters. The optimal HMM filter presented in the paper is computationally of order N3 per time instant, where N is the number of message symbols. Previously, optimal HMM filters have been of computational order N4 per time instant. Also, suboptimal HMM filters have be proposed of computation order N2 per time instant. The approach presented in this paper uses two coupled HMM filters and exploits knowledge of ...

Reduced complexity on-line estimation of hidden Markov model parameters

In this paper we propose and study low complexity algorithms for on-line estimation of hidden Markov model (HMM) parameters. The estimates approach the true model parameters as the measurement noise approaches zero, but otherwise give improved estimates, albeit with bias. On a nite data set in the high noise case, the bias may not be signi cantly more severe than for a higher complexity asymptotically optimal scheme. Our algorithms require O(N3) calculations per time instant, where N is the number of states. Previous algorithms based on earlier hidden Markov model signal processing methods, including the expectation-maximumisation (EM) algorithm require O(N4) calculations per time instant.

Optimal HMM filtering and decision feedback equalisation for differential encoded transmission systems

In this paper conditional hidden Markov model (HMM) filters and conditional Kalman filters (KF) are coupled together to improve demodulation of differential encoded signals in noisy fading channels. We present an indicator matrix representation for differential encoded signals and the optimal HMM filter for demodulation. The filter requires O(N3) calculations per time iteration, where N is the number of message symbols. Decision feedback equalisation is investigated via coupling the optimal HMM filter for estimating the message, conditioned on estimates of the channel parameters, and a KF for estimating the channel states, conditioned on soft information message estimates. The particular differential encoding scheme examined in this paper is differential phase shift keying. However, the techniques developed can be extended to other forms of differential modulation. The channel model we use allows for multiplicative channel distortions and additive white Gaussian noise. Simulation studies are also presented.

Primary and substance-induced psychotic disorders in methamphetamine users

This study investigates the rates of primary psychotic disorders (PPD) and substance induced psychotic disorders (SIPDs) in methamphetamine (MA) users accessing needle and syringe programs (NSPs). The aim was to determine if there are systematic differences in the characteristics of MA users with PPDs and SIPDs compared to those with no psychotic disorder. Participants were 198 MA users reporting use in the previous month. Diagnosis was determined using the Psychiatric Research Interview for DSM-IV Substance and Mental Disorders (PRISM-IV). Current psychiatric symptoms and substance use were also measured. Just over half (N=101) of participants met DSM-IV criteria for a lifetime psychotic disorder, including 81 (80%) with a SIPD and 20 (20%) with a PPD. Those with a younger age of onset of weekly MA use were at increased risk of a lifetime SIPD. A current psychotic disorder was found in 62 (39%), comprising 49 SIPDs (79%) and 13 PPDs (21%). MA users with a current PPD were more likely to have received psychiatric treatment in the past month than those with a current SIPD, despite a similar level of psychotic symptom severity. A high proportion of MA users accessing NSPs have psychotic disorders, the majority of which are substance-induced.

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.

Translating biomechanics research into clinical practice

There is a large amount of research conducted each year examining every aspect of the mechanics of the human body and its interaction with medical devices and the environment; from the cellular level through to the whole body. While, as researchers, we obtain great pleasure from conducting studies and creating new knowledge we need to keep in mind that while this is a good thing it is even better if this new knowledge can lead to improvement in the quality of life for individuals suffering from biomechanical disorders. Such that while commercialisation is a good aim, not all research leads to marketable outcomes. However, it can lead to improvements in surgical techniques and clinical practice. It is important for us to identify and promote how the outcomes of research lead to improvements in quality of care, as this is perhaps the most important outcome for individual patients.

Novel loci affecting iron homeostasis and their effects in individuals at risk for hemochromatosis

Variation in body iron is associated with or causes diseases, including anaemia and iron overload. Here, we analyse genetic association data on biochemical markers of iron status from 11 European-population studies, with replication in eight additional cohorts (total up to 48,972 subjects). We find 11 genome-wide-significant (P<5 × 10−8) loci, some including known iron-related genes (​HFE, ​SLC40A1, ​TF, ​TFR2, ​TFRC, ​TMPRSS6) and others novel (​ABO, ​ARNTL, ​FADS2, ​NAT2, ​TEX14). SNPs at ​ARNTL, ​TF, and ​TFR2 affect iron markers in ​HFE C282Y homozygotes at risk for hemochromatosis. There is substantial overlap between our iron loci and loci affecting erythrocyte and lipid phenotypes. These results will facilitate investigation of the roles of iron in disease.