CT and radiographic analysis of sagittal profile changes after thoracoscopic anterior scoliosis surgery

Background. Previous studies report an increase in thoracic kyphosis after anterior approaches and a flattening of sagittal contours following posterior approaches. Difficulties with measuring sagittal parameters on radiographs are avoided with reformatted sagittal CT reconstructions due to the superior endplate clarity afforded by this imaging modality. Methods. A prospective study of 30 Lenke 1 adolescent idiopathic scoliosis (AIS) patients receiving selective thoracoscopic anterior spinal fusion (TASF) was performed. Participants had ethically approved low dose CT scans at minimum 24 months after surgery in addition to their standard care following surgery. The change in sagittal contours on supine CT was compared to standing radiographic measurements of the same patients and with previous studies. Inter-observer variability was assessed as well as whether hypokyphotic and normokyphotic patient groups responded differently to the thoracoscopic anterior approach. Results. Mean T5-12 kyphosis Cobb angle increased by 11.8 degrees and lumbar lordosis increased by 5.9 degrees on standing radiographs two years after surgery. By comparison, CT measurements of kyphosis and lordosis increased by 12.3 degrees and 7.0 degrees respectively. 95% confidence intervals for inter-observer variability of sagittal contour measurements on supine CT ranged between 5-8 degrees. TASF had a slightly greater corrective effect on patients who were hypokyphotic before surgery compared with those who were normokyphotic. Conclusions. Restoration of sagittal profile is an important goal of scoliosis surgery, but reliable measurement with radiographs suffers from poor endplate clarity. TASF significantly improves thoracic kyphosis and lumbar lordosis while preserving proximal and distal junctional alignment in thoracic AIS patients. Supine CT allows greater endplate clarity for sagittal Cobb measurements and linear relationships were found between supine CT and standing radiographic measurements. In this study, improvements in sagittal kyphosis and lordosis following surgery were in agreement with prior anterior surgery studies, and add to the current evidence suggesting that anterior correction is more capable than posterior approaches of addressing the sagittal component of both the instrumented and adjacent non instrumented segments following surgical correction of progressive Lenke 1 idiopathic scoliosis.

Bayesian experimental design for models with intractable likelihoods

In this paper we present a methodology for designing experiments for efficiently estimating the parameters of models with computationally intractable likelihoods. The approach combines a commonly used methodology for robust experimental design, based on Markov chain Monte Carlo sampling, with approximate Bayesian computation (ABC) to ensure that no likelihood evaluations are required. The utility function considered for precise parameter estimation is based upon the precision of the ABC posterior distribution, which we form efficiently via the ABC rejection algorithm based on pre-computed model simulations. Our focus is on stochastic models and, in particular, we investigate the methodology for Markov process models of epidemics and macroparasite population evolution. The macroparasite example involves a multivariate process and we assess the loss of information from not observing all variables.

GABAergic modification of myopic and hyperopic ocular tissue effects on scleral fibroblast growth

Purpose: Myopia is a common eye disorder affecting up to 90% of children in South East Asia and 30% of the population worldwide. Myopia of high severity is a leading cause of blindness around the world (4th to 5th most common). Changes and remodelling of the sclera i.e. increase cellular proliferation & increase protein synthesis within scleral cells (↑ scleral DNA) and thinning and lose of extracellular matrix of sclera (↓ scleral GAG synthesis) have been linked to myopic eye growth in animal models. Signals acting on the sclera are thought to originate in the retina, and are modulated by the retinal pigment epithelium (RPE) with limited evidence suggesting that the RPE can modify scleral cell growth in culture. However, the mechanism of retinal signal transmission and the role of posterior eye cup tissue, including the RPE, in mediating changes in scleral fibroblast growth during myopia development are unclear. Retinal transmitter systems are critically involved in pathways regulating eye growth, which ultimately lead to alterations in the sclera if eye size is to change. A dopaminergic agonist and muscarinic antagonists decrease the proliferation of scleral chondrocytes when co-cultured with chick’s retinal pigment epithelium (RPE). GABA receptors have recently been localised to chick sclera. We therefore hypothesised that posterior eye cup tissue from myopic eyes would stimulate and from hyperopic eyes would inhibit growth of scleral fibroblasts in vitro and that GABAergic agents could directly interact with scleral cells or indirectly modify the effects of myopic and hyperopic posterior eye cup tissue on scleral fibroblast growth. Method: Fibroblastic cells obtained from 8-day-old chick sclera were used to establish cell banks. Two major experiments were performed. Experiment 1: To determine if posterior eye cup tissues from myopic eye stimulates and hyperopic eye inhibits scleral cell proliferation, when co-cultured with scleral cells in vitro. This study comprised two linked experiments, i) monocular visual treatments of FDM (form-deprivation myopia), LIM (lens-induced myopia) and LIH (lens-induced hyperopia) with assessment of the effect of full punch eye cup tissue on DNA and GAG synthesis by cultured chick scleral fibroblasts, and ii) binocular visual treatments comprising LIM and LIH with assessment of the effect of individual layers of eye cup tissues (neural retina, RPE and choroid) on cultured chick scleral fibroblasts. Visual treatment was applied for 3 days. Experiment 2: To determine the direct interaction of GABA agents on scleral cell growth and to establish whether GABA agents modify the stimulatory/inhibitory effect of myopic and hyperopic posterior eye cup tissues on cultured scleral cell growth in vitro. Two linked experiments were performed. i) GABA agonists (muscimol and baclofen) and GABA antagonists (bicuculine (-), CGP46381 and TPMPA) were added to scleral cell culture medium to determine their direct effect on scleral cells. ii) GABAergic agents (agonists and antagonists) were administered to scleral fibroblasts co-cultured with posterior eye cup tissue (retina, RPE, retina/RPE, RPE/choroid). Ocular tissues were obtained from chick eyes wearing +15D (LIH) or -15D lenses (LIM) for 3 days. In both experiments, tissues were added to hanging cell culture insert (pore size 1.0ìm) placed over each well of 24 well plates while scleral cells were cultured in DMEM/F12, Glutamax (Gibco) plus 10% FBS and penicillin/streptomycin (50U/ml)) and fungizone (1.25ug/ml) (Gibco), at seeding density of 30,000 cells/well at the bottom of the well and allowed to grow for 3 days. Scleral cells proliferation rate throughout the study was evaluated by determining GAG and DNA content of scleral cells using Dimethylmethylene blue (DMMB) dye and Quant-iTTm Pico Green® dsDNA reagent respectively. Results and analysis: Based on DNA and GAG content, there was no significant difference in tissue effect of LIM and LIH eyes on scleral fibroblast growth (DNA: 8.4 ± 1.1μg versus 9.3 ± 2.3 μg, p=0.23; GAG: 10.13 ± 1.4 μg versus 12.67 ± 1.2 μg, F2,23=6.16, p=0.0005) when tissues were obtained from monocularly treated chick eyes (FDM or +15D lens or -15D lens over right eyes with left eyes untreated) and co-cultured as full punch. When chick eyes were treated binocularly with -15D lens (LIM) right eye and +15D lens (LIH) left eyes and tissue layers were separated, the retina from LIM eyes did not stimulate scleral cell proliferation compared to LIH eyes (DNA: 27.2 ± 6.7 μg versus 23.2 ± 1.5 μg, p=0.23; GAG: 28.1 ±3.7 μg versus 28.7 ± 4.2 μg, p=0.21). Similarly, the LIH and LIM choroid did not produce a differential effect based on DNA (LIM 46.9 ± 6.4 μg versus LIH 53.5 ± 4.7 μg, p=0.18), however the choroid from LIH eyes induced higher scleral GAG content than from LIM eyes (32.5 ± 6.7 μg versus 18.9 ± 1.2 μg, p=0.023). In contrast, the RPE from LIM eyes caused a significant increase in fibroblast proliferation whereas the RPE from LIH eyes was relatively inhibitory (72.4 ± 6.3 μg versus 27.9 ± 2.3 μg, F1, 6=69.99, p=0.0005). GAG data were opposite to DNA data e.g. the RPE from LIH eyes increased (33.7 ± 7.9 μg) while the RPE from LIM eyes decreased (28.2 ± 3.0 μg) scleral cell growth (F1, 6=13.99, p=0.010). Based on DNA content, GABA agents had a small direct effect on scleral cell growth; GABA agonists increased (21.4 ± 1.0% and 18.3 ± 1.0% with muscimol and baclofen, p=0.0021), whereas GABA antagonists decreased fibroblast proliferation (-23.7 ± 0.9% with bicuculine & CGP46381 and -28.1 ± 0.5% with TPMPA, p=0.0004). GABA agents also modified the effect of LIM and LIH tissues (p=0.0005).The increase in proliferation rate of scleral fibroblasts co-cultured with tissues (RPE, retina, RPE/retina and RPE/choroid) from LIM treated eyes was enhanced by GABA agonists (muscimol: 27.4 ± 1.2%, 35.8 ± 1.6%, 8.4 ± 0.3% and 11.9 ± 0.6%; baclofen: 27.0 ± 1.0%, 15.8 ± 1.5%, 16.8 ± 1.2% and 15.4 ± 0.4%, p=0.014) whereas GABA antagonists further reduced scleral fibroblasts growth (bicuculine: -52.5 ± 2.5%, -36.9 ± 1.4%, -37.5 ± 0.6% and -53.7 ± 0.9%; TPMPA: 57.3 ± 1.3%, -15.7 ± 1.2%, -33.5 ± 0.4% and -45.9 ± 1.5%; CGP46381: -51.9 ± 1.6%, -28.5 ± 1.5%, -25.4 ± 2.0% and -45.5 ± 1.9% respectively, p=0.0034). GAG data were opposite to DNA data throughout the experiment e.g. GABA agonists further inhibited while antagonists relatively enhanced scleral fibroblasts growth for both LIM and LIH tissue co-culture. The effect of GABA agents was relatively lower (p=0.0004) for tissue from LIH versus LIM eyes but was in a similar direction. There was a significant drug effect on all four tissue types e.g. RPE, retina, RPE/retina and RPE/choroid for both LIM and LIH tissue co-culture (F20,92=3.928, p=0.0005). However, the effect of GABA agents was greatest in co-culture with RPE tissue (F18,36=4.865, p=0.0005). Summary and Conclusion: 1) Retinal defocus signals are transferred to RPE and choroid which then exert their modifying effect on scleral GAG and DNA synthesis either through growth stimulating factors or directly interacting with scleral cells in process of scleral remodeling during LIM and LIH visual conditions. 2) GABAergic agents affect the proliferation of scleral fibroblasts both directly and when co-cultured with ocular tissues in vitro.

Marginal reversible jump Markov chain Monte Carlo with application to motor unit number estimation

Motor unit number estimation (MUNE) is a method which aims to provide a quantitative indicator of progression of diseases that lead to loss of motor units, such as motor neurone disease. However the development of a reliable, repeatable and fast real-time MUNE method has proved elusive hitherto. Ridall et al. (2007) implement a reversible jump Markov chain Monte Carlo (RJMCMC) algorithm to produce a posterior distribution for the number of motor units using a Bayesian hierarchical model that takes into account biological information about motor unit activation. However we find that the approach can be unreliable for some datasets since it can suffer from poor cross-dimensional mixing. Here we focus on improved inference by marginalising over latent variables to create the likelihood. In particular we explore how this can improve the RJMCMC mixing and investigate alternative approaches that utilise the likelihood (e.g. DIC (Spiegelhalter et al., 2002)). For this model the marginalisation is over latent variables which, for a larger number of motor units, is an intractable summation over all combinations of a set of latent binary variables whose joint sample space increases exponentially with the number of motor units. We provide a tractable and accurate approximation for this quantity and also investigate simulation approaches incorporated into RJMCMC using results of Andrieu and Roberts (2009).

Corneal changes following short-term rigid contact lens wear

Purpose: The aim of this cross-over study was to investigate the changes in corneal thickness, anterior and posterior corneal topography, corneal refractive power and ocular wavefront aberrations, following the short term use of rigid contact lenses. Method: Fourteen participants wore 4 different types of contact lenses (RGP lenses of 9.5 mm and 10.5 mm diameter, and for comparison a PMMA lens of 9.5 mm diameter and a soft silicone hydrogel lens) on 4 different days for a period of 8 h on each day. Measures were collected before and after contact lens wear and additionally on a baseline day. Results: Anterior corneal curvature generally showed a flattening with both of the RGP lenses and a steepening with the PMMA lens. A significant negative correlation was found between the change in corneal swelling and central and peripheral posterior corneal curvature (all p ≤ 0.001). RGP contact lenses caused a significant decrease in corneal refractive power (hyperopic shift) of approximately 0.5 D. The PMMA contact lenses caused the greatest corneal swelling in both the central (27.92 ± 15.49 μm, p < 0.001) and peripheral (17.78 ± 12.11 μm, p = 0.001) corneal regions, a significant flattening of the posterior cornea and an increase in ocular aberrations (all p ≤ 0.05). Conclusion: The corneal swelling associated with RGP lenses was relatively minor, but there was slight central corneal flattening and a clinically significant hyperopic change in corneal refractive power after the first day of lens wear. The PMMA contact lenses resulted in significant corneal swelling and reduced optical performance of the cornea.

Qualitative study of adversity activated development and resilience in adolescents and young adults with congenital heart disease and their parents

What is it like to have a medical condition that few people have ever heard about? How does it feel to have to question whether daily physical activities are dangerous for you, whilst you watch your friends enjoy those activities without a care? Can you imagine that you need to have a complicated heart surgery, with risks such as paralysis or death? Or even imagine facing the painful recovery period and scars after such a surgery? Then imagine that you are a child or teenager dealing with this medical condition when all your friends are simply occupied with school and normal life. Now consider that surgery has been undertaken to extend your lifespan, but the operation is so new that the long-term outcomes are just not known? All you really know is that you might have ‘surgical repairs’ to your heart and symptoms may be relieved or managed by medications or cardiac devices, but you are never going to be cured. What if you had already experienced painful, frightening, lonely and tedious hospitalisations and you were forced to put your life on hold to re-enter that situation, time and time again. This may be your life, as a Congenital Heart Disease or CHD patient. How do such patients cope and in many cases even thrive? This chapter will review current international literature regarding the medical and personal impact of CHD. Our qualitative study of the perspectives of young CHD patients and their parents contributes to the Australian story of CHD, as well as highlighting the potential for CHD related adversity to promote personal development.

Bayesian algorithms with applications

Advances in algorithms for approximate sampling from a multivariable target function have led to solutions to challenging statistical inference problems that would otherwise not be considered by the applied scientist. Such sampling algorithms are particularly relevant to Bayesian statistics, since the target function is the posterior distribution of the unobservables given the observables. In this thesis we develop, adapt and apply Bayesian algorithms, whilst addressing substantive applied problems in biology and medicine as well as other applications. For an increasing number of high-impact research problems, the primary models of interest are often sufficiently complex that the likelihood function is computationally intractable. Rather than discard these models in favour of inferior alternatives, a class of Bayesian "likelihoodfree" techniques (often termed approximate Bayesian computation (ABC)) has emerged in the last few years, which avoids direct likelihood computation through repeated sampling of data from the model and comparing observed and simulated summary statistics. In Part I of this thesis we utilise sequential Monte Carlo (SMC) methodology to develop new algorithms for ABC that are more efficient in terms of the number of model simulations required and are almost black-box since very little algorithmic tuning is required. In addition, we address the issue of deriving appropriate summary statistics to use within ABC via a goodness-of-fit statistic and indirect inference. Another important problem in statistics is the design of experiments. That is, how one should select the values of the controllable variables in order to achieve some design goal. The presences of parameter and/or model uncertainty are computational obstacles when designing experiments but can lead to inefficient designs if not accounted for correctly. The Bayesian framework accommodates such uncertainties in a coherent way. If the amount of uncertainty is substantial, it can be of interest to perform adaptive designs in order to accrue information to make better decisions about future design points. This is of particular interest if the data can be collected sequentially. In a sense, the current posterior distribution becomes the new prior distribution for the next design decision. Part II of this thesis creates new algorithms for Bayesian sequential design to accommodate parameter and model uncertainty using SMC. The algorithms are substantially faster than previous approaches allowing the simulation properties of various design utilities to be investigated in a more timely manner. Furthermore the approach offers convenient estimation of Bayesian utilities and other quantities that are particularly relevant in the presence of model uncertainty. Finally, Part III of this thesis tackles a substantive medical problem. A neurological disorder known as motor neuron disease (MND) progressively causes motor neurons to no longer have the ability to innervate the muscle fibres, causing the muscles to eventually waste away. When this occurs the motor unit effectively ‘dies’. There is no cure for MND, and fatality often results from a lack of muscle strength to breathe. The prognosis for many forms of MND (particularly amyotrophic lateral sclerosis (ALS)) is particularly poor, with patients usually only surviving a small number of years after the initial onset of disease. Measuring the progress of diseases of the motor units, such as ALS, is a challenge for clinical neurologists. Motor unit number estimation (MUNE) is an attempt to directly assess underlying motor unit loss rather than indirect techniques such as muscle strength assessment, which generally is unable to detect progressions due to the body’s natural attempts at compensation. Part III of this thesis builds upon a previous Bayesian technique, which develops a sophisticated statistical model that takes into account physiological information about motor unit activation and various sources of uncertainties. More specifically, we develop a more reliable MUNE method by applying marginalisation over latent variables in order to improve the performance of a previously developed reversible jump Markov chain Monte Carlo sampler. We make other subtle changes to the model and algorithm to improve the robustness of the approach.

Polycaprolactone-based scaffold plus recombinant human bone morphogenic protein rhBMP-2) in a sheep thoracic spine fusion model

Adolescent idiopathic scoliosis is a complex three dimensional deformity affecting 2-3% of the general population. The resulting spinal deformity consists of coronal curvature, hypokyphosis of the thoracic spine and vertebral rotation in the axial plane with posterior elements turned into the curve concavity. The potential for curve progression is heightened during the adolescent growth spurt. Success of scoliosis deformity correction depends on solid bony fusion between adjacent vertebrae after the intervertebral (IV) discs have been surgically cleared and the disc spaces filled with graft material. Recently a bioactive and resorbable scaffold fabricated from medical grade polycaprolactone has been developed for bone regeneration at load bearing sites. Combined with rhBMP-2, this has been shown to be successful in acting as a bone graft substitute in a porcine lumbar interbody fusion model when compared to autologous bone graft alone. The study aimed to establish a large animal thoracic spine interbody fusion model, develop spine biodegradable scaffolds (PCL) in combination with biologics (rhBMP-2) and to establish a platform for research into spine tissue engineering constructs. Preliminary results demonstrate higher grades of radiologically evident bony fusion across all levels when comparing fusion scores between the 3 and 6 month postop groups at the PCL CaP coated scaffold level, which is observed to be a similar grade to autograft, while no fusion is seen at the scaffold only level. Results to date suggest that the combination of rhBMP-2 and scaffold engineering actively promotes bone formation, laying the basis of a viable tissue engineered constructs.

Biomechanical performance of polycaprolactone (PCL)-based scaffold with rhBMP-2 in a sheep thoracic spine fusion model

Adolescent idiopathic scoliosis is a complex three dimensional deformity affecting 2-3% of the general population. Resulting spine deformities include progressive coronal curvature, hypokyphosis, or frank lordosis in the thoracic spine and vertebral rotation in the axial plane with posterior elements turned into the curve concavity. The potential for curve progression is heightened during the adolescent growth spurt. Success of scoliosis deformity correction depends on solid bony fusion between adjacent vertebrae after the intervertebral discs have been surgically cleared and the disc spaces filled with graft material. Problems with bone graft harvest site morbidity as well as limited bone availability have led to the search for bone graft substitutes. Recently, a bioactive and resorbable scaffold fabricated from medical grade polycaprolactone (PCL) has been developed for bone regeneration at load bearing sites. Combined with recombinant human bone morphogenic protein–2 (rhBMP-2), this has been shown to be successful in acting as a bone graft substitute in acting as a bone graft substitute in a porcine lumbar interbody fusion model when compared to autologous bone graft. This in vivo sheep study intends to evaluate the suitability of a custom designed medical grade PCL scaffold in combination with rhBMP-2 as a bone graft substitute in the setting of mini–thoracotomy surgery as a platform for ongoing research to benefit patients with adolescent idiopathic scoliosis.

Methods to estimate the size and shape of the unaccommodated crystalline lens in vivo

Purpose. The purpose of this article was to present methods capable of estimating the size and shape of the human eye lens without resorting to phakometry or magnetic resonance imaging (MRI). Methods. Previously published biometry and phakometry data of 66 emmetropic eyes of 66 subjects (age range [18, 63] years, spherical equivalent range [−0.75, +0.75] D) were used to define multiple linear regressions for the radii of curvature and thickness of the lens, from which the lens refractive index could be derived. MRI biometry was also available for a subset of 30 subjects, from which regressions could be determined for the vertex radii of curvature, conic constants, equatorial diameter, volume, and surface area. All regressions were compared with the phakometry and MRI data; the radii of curvature regressions were also compared with a method proposed by Bennett and Royston et al. Results. The regressions were in good agreement with the original measurements. This was especially the case for the regressions of lens thickness, volume, and surface area, which each had an R2 > 0.6. The regression for the posterior radius of curvature had an R2 < 0.2, making this regression unreliable. For all other regressions we found 0.25 < R2 < 0.6. The Bennett-Royston method also produced a good estimation of the radii of curvature, provided its parameters were adjusted appropriately. Conclusions. The regressions presented in this article offer a valuable alternative in case no measured lens biometry values are available; however care must be taken for possible outliers.

Exact and approximate Bayesian inference for low count time series models with intractable likelihoods

In this paper we present a new simulation methodology in order to obtain exact or approximate Bayesian inference for models for low-valued count time series data that have computationally demanding likelihood functions. The algorithm fits within the framework of particle Markov chain Monte Carlo (PMCMC) methods. The particle filter requires only model simulations and, in this regard, our approach has connections with approximate Bayesian computation (ABC). However, an advantage of using the PMCMC approach in this setting is that simulated data can be matched with data observed one-at-a-time, rather than attempting to match on the full dataset simultaneously or on a low-dimensional non-sufficient summary statistic, which is common practice in ABC. For low-valued count time series data we find that it is often computationally feasible to match simulated data with observed data exactly. Our particle filter maintains $N$ particles by repeating the simulation until $N+1$ exact matches are obtained. Our algorithm creates an unbiased estimate of the likelihood, resulting in exact posterior inferences when included in an MCMC algorithm. In cases where exact matching is computationally prohibitive, a tolerance is introduced as per ABC. A novel aspect of our approach is that we introduce auxiliary variables into our particle filter so that partially observed and/or non-Markovian models can be accommodated. We demonstrate that Bayesian model choice problems can be easily handled in this framework.

The relationship between hip abductor muscle strength and magnitude of pelvic drop following a 3 week strengthening protocol in non-specific low back pain patients.

Purpose: To examine the relationship between hip abductor muscle (HABD) strength and the magnitude of pelvic drop (MPD) for patients with non-specific low back pain (NSLBP) and controls (CON) prior to and following a 3-week HABD strengthening protocol. At baseline, we hypothesized that NSLBP patients would exhibit reduced HABD strength and greater MPD compared to CON. Following the protocol, we hypothesized that strength would increase and MPD would decrease. Relevance: The Trendelenburg test (TT) is a common clinical test used to examine the ability of the HABD to maintain horizontal pelvic position during single limb stance. However, no study has specifically tested this theory. Moreover, no study has investigated the relationship between HABD strength and pelvic motion during walking or tested whether increased HABD strength would reduce the MPD. Methods: Quasi-experimental with 3-week exercise intervention. Fifteen NSLBP patients (32.5yrs,range 21-51yrs; VAS baseline: 5.3cm) and 10 CON (29.5yrs,range 22-47yrs) were recruited. Isometric HABD strength was measured using a force dynamometer and the average of three maximal voluntary contractions were normalized to body mass (N/kg). Two-dimensional MPD (degrees) was measured using a 60 Hz camera and was derived from two retroreflective-markers placed on the posterior superior iliac spines. MPD was measured while performing the static TT and while walking and averaged over 10 consecutive footfalls. NSLBP patients completed a 3-week HABD strengthening protocol consisting of 2 open-kinetic-chain exercises then all measures were repeated. Non-parametric analysis was used for group comparisons and correlation analysis. Results: At baseline, the NSLBP patients demonstrated 31% reduced HABD strength (mean=6.6 N/kg) compared to CON (mean=9.5 N/kg: p=0.03) and no significant differences in maximal pelvic frontal plane excursion while walking (NSLBP:mean=8.1°, CON:mean=7.1°: p=0.72). No significant correlations were measured between left HABD strength and right MPD (r=-0.37, p=0.11), or between right HABD strength and left MPD (r=-0.04, p=0.84) while performing the static TT. Following the 3-week strengthening protocol, NSLBP patients demonstrated a 12% improvement in strength (Post:mean=7.4 N/kg: p=0.02), a reduction in pain (VAS followup: 2.8cm), but no significant decreases in MPD while walking (p=0.92). Conclusions: NSLBP patients demonstrated reduced HABD strength at baseline and were able to increase strength and reduce pain in a 3-week period. However, despite increases in HABD strength, the NSLBP group exhibited similar MPD motion during the static TT and while walking compared to baseline and controls. Implications: The results suggest that the HABD alone may not be primarily responsible for controlling a horizontal pelvic position during static and dynamic conditions. Increasing the strength of the hip abductors resulted in a reduction of pain in NSLBP patients providing evidence for further research to identify specific musculature responsible for controlling pelvic motion.

The use of electronic portal imaging to verify patient position during intensity-modulated radiotherapy delivered by the dynamic MLC technique

Purpose: The precise shape of the three-dimensional dose distributions created by intensity-modulated radiotherapy means that the verification of patient position and setup is crucial to the outcome of the treatment. In this paper, we investigate and compare the use of two different image calibration procedures that allow extraction of patient anatomy from measured electronic portal images of intensity-modulated treatment beams. Methods and Materials: Electronic portal images of the intensity-modulated treatment beam delivered using the dynamic multileaf collimator technique were acquired. The images were formed by measuring a series of frames or segments throughout the delivery of the beams. The frames were then summed to produce an integrated portal image of the delivered beam. Two different methods for calibrating the integrated image were investigated with the aim of removing the intensity modulations of the beam. The first involved a simple point-by-point division of the integrated image by a single calibration image of the intensity-modulated beam delivered to a homogeneous polymethyl methacrylate (PMMA) phantom. The second calibration method is known as the quadratic calibration method and required a series of calibration images of the intensity-modulated beam delivered to different thicknesses of homogeneous PMMA blocks. Measurements were made using two different detector systems: a Varian amorphous silicon flat-panel imager and a Theraview camera-based system. The methods were tested first using a contrast phantom before images were acquired of intensity-modulated radiotherapy treatment delivered to the prostate and pelvic nodes of cancer patients at the Royal Marsden Hospital. Results: The results indicate that the calibration methods can be used to remove the intensity modulations of the beam, making it possible to see the outlines of bony anatomy that could be used for patient position verification. This was shown for both posterior and lateral delivered fields. Conclusions: Very little difference between the two calibration methods was observed, so the simpler division method, requiring only the single extra calibration measurement and much simpler computation, was the favored method. This new method could provide a complementary tool to existing position verification methods, and it has the advantage that it is completely passive, requiring no further dose to the patient and using only the treatment fields.

Pseudo-marginal algorithms with multiple CPUs

The emergence of pseudo-marginal algorithms has led to improved computational efficiency for dealing with complex Bayesian models with latent variables. Here an unbiased estimator of the likelihood replaces the true likelihood in order to produce a Bayesian algorithm that remains on the marginal space of the model parameter (with latent variables integrated out), with a target distribution that is still the correct posterior distribution. Very efficient proposal distributions can be developed on the marginal space relative to the joint space of model parameter and latent variables. Thus psuedo-marginal algorithms tend to have substantially better mixing properties. However, for pseudo-marginal approaches to perform well, the likelihood has to be estimated rather precisely. This can be difficult to achieve in complex applications. In this paper we propose to take advantage of multiple central processing units (CPUs), that are readily available on most standard desktop computers. Here the likelihood is estimated independently on the multiple CPUs, with the ultimate estimate of the likelihood being the average of the estimates obtained from the multiple CPUs. The estimate remains unbiased, but the variability is reduced. We compare and contrast two different technologies that allow the implementation of this idea, both of which require a negligible amount of extra programming effort. The superior performance of this idea over the standard approach is demonstrated on simulated data from a stochastic volatility model.

Development of an ultra-thin fibroin membrane for RPE cell transplantation

Purpose: One of the challenges associated with cell-based therapies for repairing the retina is the development of suitable materials on which to grow and transplant retinal cells. Using the ARPE-19 cell line, we have previously demonstrated the feasibility of growing RPE-derived cells on membranes prepared from the silk protein fibroin. The present study was aimed at developing a porous, ultra-thin fibroin membrane that might better support development of apical-basal polarity in culture, and to extend this work to primary cultures of human RPE cells. Methods: Ultra-thin fibroin membranes were prepared using a highly polished casting table coated with Topas® (a cyclic olefin copolymer) and a 1:0.03 aqueous solution of fibroin and PEO (Mv 900 000 g/mol). Following drying, the membranes were water annealed to make them water-stable, washed in water to remove PEO, sterilised by treatment with 95% ethanol, and washed extensively in saline. Primary cultures containing human RPE cells were established from donor posterior eye cups and maintained in DMEM/F12 medium supplemented with 10% fetal bovine serum and antibiotics. First passage cultures were seeded onto fibroin membranes pre-coated with vitronectin and grown for 6 weeks in medium supplemented with 1% serum. Comparative cultures were established on porous 1.0 µm pore PET membrane (Millipore) and using ARPE-19 cells. Results: The fibroin membranes displayed an average thickness of 3 µm and contained numerous dimples/pore-like structures of up to 3-5 µm in diameter. The primary cultures predominantly contained pigmented epithelial cells, but mesenchymal cells (presumed fibroblasts) were also often present. Passaged cultures appeared to attach equally well to either fibroin or PET membranes. Over time cells on either material adopted a more cobblestoned morphology. Conclusions: Progress has been made towards developing a porous ultra-thin fibroin membrane that supports cultivation of RPE cells. Further studies are required to determine the degree of membrane permeability and RPE polarity.