Partial intervertebral fusion secures successful outcomes after thoracoscopic anterior scoliosis correction: A low-dose computed tomography study

Study Design Retrospective review of prospectively collected data. Objectives To analyze intervertebral (IV) fusion after thoracoscopic anterior spinal fusion (TASF) and explore the relationship between fusion scores and key clinical variables. Summary of Background Information TASF provides comparable correction with some advantages over posterior approaches but reported mechanical complications, and their relationship to non-union and graft material is unclear. Similarly, the optimal combination of graft type and implant stiffness for effecting successful radiologic union remains undetermined. Methods A subset of patients from a large single-center series who had TASF for progressive scoliosis underwent low-dose computed tomographic scans 2 years after surgery. The IV fusion mass in the disc space was assessed using the 4-point Sucato scale, where 1 indicates <50% and 4 indicates 100% bony fusion of the disc space. The effects of rod diameter, rod material, graft type, fusion level, and mechanical complications on fusion scores were assessed. Results Forty-three patients with right thoracic major curves (mean age 14.9 years) participated in the study. Mean fusion scores for patient subgroups ranged from 1.0 (IV levels with rod fractures) to 2.2 (4.5-mm rod with allograft), with scores tending to decrease with increasing rod size and stiffness. Graft type (autograft vs. allograft) did not affect fusion scores. Fusion scores were highest in the middle levels of the rod construct (mean 2.52), dropping off by 20% to 30% toward the upper and lower extremities of the rod. IV levels where a rod fractured had lower overall mean fusion scores compared to levels without a fracture. Mean total Scoliosis Research Society (SRS) questionnaire scores were 98.9 from a possible total of 120, indicating a good level of patient satisfaction. Conclusions Results suggest that 100% radiologic fusion of the entire disc space is not necessary for successful clinical outcomes following thoracoscopic anterior selective thoracic fusion.

Quantifying progressive anterior overgrowth in the thoracic vertebrae of adolescent idiopathic scoliosis patients

Study design Anterior and posterior vertebral body heights were measured from sequential MRI scans of adolescent idiopathic scoliosis (AIS) patients and healthy controls. Objective To measure changes in vertebral body height over time during scoliosis progression to assess how vertebral body height discrepancies change during growth. Summary of background data Relative anterior overgrowth has been proposed as a potential driver for AIS initiation and progression. This theory proposes that the anterior column grows faster, and the posterior column slower, in AIS patients when compared to healthy controls. There is disagreement in the literature as to whether the anterior vertebral body heights are proportionally greater than posterior vertebral body heights in AIS patients when compared to healthy controls. To some extent, these discrepancies may be attributed to methodological differences. Methods MRI scans of the major curve of 21 AIS patients (mean age 12.5 ± 1.4 years, mean Cobb 32.2 ± 12.8º) and between T4 and T12 of 21 healthy adolescents (mean age 12.1 ± 0.5 years) were captured for this study. Of the 21 AIS patients, 14 had a second scan on average 10.8 ± 4.7 months after the first. Anterior and posterior vertebral body heights were measured from the true sagittal plane of each vertebra such that anterior overgrowth could be quantified. Results The difference between anterior and posterior vertebral body height in healthy, non-scoliotic children was significantly greater than in AIS patients with mild to moderate scoliosis. However there was no significant relationship between the overall anterior-posterior vertebral body height difference in AIS and either severity of the curve or its progression over time. Conclusions Whilst AIS patients have a proportionally longer anterior column than non-scoliotic controls, the degree of anterior overgrowth was not related to the rate of progression or the severity of the scoliotic curve.

Overfitting Bayesian mixture models with an unknown number of components

This paper proposes solutions to three issues pertaining to the estimation of finite mixture models with an unknown number of components: the non-identifiability induced by overfitting the number of components, the mixing limitations of standard Markov Chain Monte Carlo (MCMC) sampling techniques, and the related label switching problem. An overfitting approach is used to estimate the number of components in a finite mixture model via a Zmix algorithm. Zmix provides a bridge between multidimensional samplers and test based estimation methods, whereby priors are chosen to encourage extra groups to have weights approaching zero. MCMC sampling is made possible by the implementation of prior parallel tempering, an extension of parallel tempering. Zmix can accurately estimate the number of components, posterior parameter estimates and allocation probabilities given a sufficiently large sample size. The results will reflect uncertainty in the final model and will report the range of possible candidate models and their respective estimated probabilities from a single run. Label switching is resolved with a computationally light-weight method, Zswitch, developed for overfitted mixtures by exploiting the intuitiveness of allocation-based relabelling algorithms and the precision of label-invariant loss functions. Four simulation studies are included to illustrate Zmix and Zswitch, as well as three case studies from the literature. All methods are available as part of the R package Zmix, which can currently be applied to univariate Gaussian mixture models.

Accelerating Pseudo-Marginal MCMC using Gaussian Processes

Pseudo-marginal methods such as the grouped independence Metropolis-Hastings (GIMH) and Markov chain within Metropolis (MCWM) algorithms have been introduced in the literature as an approach to perform Bayesian inference in latent variable models. These methods replace intractable likelihood calculations with unbiased estimates within Markov chain Monte Carlo algorithms. The GIMH method has the posterior of interest as its limiting distribution, but suffers from poor mixing if it is too computationally intensive to obtain high-precision likelihood estimates. The MCWM algorithm has better mixing properties, but less theoretical support. In this paper we propose to use Gaussian processes (GP) to accelerate the GIMH method, whilst using a short pilot run of MCWM to train the GP. Our new method, GP-GIMH, is illustrated on simulated data from a stochastic volatility and a gene network model.

The impact of inter-fraction set-up errors on the probability of pulmonary and cardiac complication in left-sided breast cancer patients

Purpose This study evaluated the impact of patient set-up errors on the probability of pulmonary and cardiac complications in the irradiation of left-sided breast cancer. Methods and Materials Using the CMS XiO Version 4.6 (CMS Inc., St Louis, MO) radiotherapy planning system's NTCP algorithm and the Lyman -Kutcher-Burman (LKB) model, we calculated the DVH indices for the ipsilateral lung and heart and the resultant normal tissue complication probabilities (NTCP) for radiation-induced pneumonitis and excess cardiac mortality in 12 left-sided breast cancer patients. Results Isocenter shifts in the posterior direction had the greatest effect on the lung V20, heart V25, mean and maximum doses to the lung and the heart. Dose volume histograms (DVH) results show that the ipsilateral lung V20 tolerance was exceeded in 58% of the patients after 1cm posterior shifts. Similarly, the heart V25 tolerance was exceeded after 1cm antero-posterior and left-right isocentric shifts in 70% of the patients. The baseline NTCPs for radiation-induced pneumonitis ranged from 0.73% - 3.4% with a mean value of 1.7%. The maximum reported NTCP for radiation-induced pneumonitis was 5.8% (mean 2.6%) after 1cm posterior isocentric shift. The NTCP for excess cardiac mortality were 0 % in 100% of the patients (n=12) before and after setup error simulations. Conclusions Set-up errors in left sided breast cancer patients have a statistically significant impact on the Lung NTCPs and DVH indices. However, with a central lung distance of 3cm or less (CLD <3cm), and a maximum heart distance of 1.5cm or less (MHD<1.5cm), the treatment plans could tolerate set-up errors of up to 1cm without any change in the NTCP to the heart.

CBCT image guidance in head and neck irradiation: The impact of daily and weekly imaging protocols

Purpose This study evaluated the impact of a daily and weekly image-guided radiotherapy protocols in reducing setup errors and setting of appropriate margins in head and neck cancer patients. Materials and methods Interfraction and systematic shifts for the hypothetical day 1–3 plus weekly imaging were extrapolated from daily imaging data from 31 patients (964 cone beam computed tomography (CBCT) scans). In addition, residual setup errors were calculated by taking the average shifts in each direction for each patient based on the first three shifts and were presumed to represent systematic setup error. The clinical target volume (CTV) to planning target volume (PTV) margins were calculated using van Herk formula and analysed for each protocol. Results The mean interfraction shifts for daily imaging were 0·8, 0·3 and 0·5 mm in the S-I (superior-inferior), L-R (left-right) and A-P (anterior-posterior) direction, respectively. On the other hand the mean shifts for day 1–3 plus weekly imaging were 0·9, 1·8 and 0·5 mm in the S-I, L-R and A-P direction, respectively. The mean day 1–3 residual shifts were 1·5, 2·1 and 0·7 mm in the S-I, L-R and A-P direction, respectively. No significant difference was found in the mean setup error for the daily and hypothetical day 1–3 plus weekly protocol. However, the calculated CTV to PTV margins for the daily interfraction imaging data were 1·6, 3·8 and 1·4 mm in the S-I, L-R and A-P directions, respectively. Hypothetical day 1–3 plus weekly resulted in CTV–PTV margins of 5, 4·2 and 5 mm in the S-I, L-R and A-P direction. Conclusions The results of this study show that a daily CBCT protocol reduces setup errors and allows setup margin reduction in head and neck radiotherapy compared to a weekly imaging protocol.

Corneal confocal microscopy predicts 4-year incident peripheral neuropathy in Type 1 diabetes

OBJECTIVE This study determined if deficits in corneal nerve fiber length (CNFL) assessed using corneal confocal microscopy (CCM) can predict future onset of diabetic peripheral neuropathy (DPN). RESEARCH DESIGN AND METHODS CNFL and a range of other baseline measures were compared between 90 nonneuropathic patients with type 1 diabetes who did or did not develop DPN after 4 years. The receiver operator characteristic (ROC) curve was used to determine the capability of single and combined measures of neuropathy to predict DPN. RESULTS DPN developed in 16 participants (18%) after 4 years. Factors predictive of 4-year incident DPN were lower CNFL (P = 0.041); longer duration of diabetes (P = 0.002); higher triglycerides (P = 0.023); retinopathy (higher on the Early Treatment of Diabetic Retinopathy Study scale) (P = 0.008); nephropathy (higher albumin-to-creatinine ratio) (P = 0.001); higher neuropathy disability score (P = 0.037); lower cold sensation (P = 0.001) and cold pain (P = 0.027) thresholds; higher warm sensation (P = 0.008), warm pain (P = 0.024), and vibration (P = 0.003) thresholds; impaired monofilament response (P = 0.003); and slower peroneal (P = 0.013) and sural (P = 0.002) nerve conduction velocity. CCM could predict the 4-year incident DPN with 63% sensitivity and 74% specificity for a CNFL threshold cutoff of 14.1 mm/mm2 (area under ROC curve = 0.66, P = 0.041). Combining neuropathy measures did not improve predictive capability. CONCLUSIONS DPN can be predicted by various demographic, metabolic, and conventional neuropathy measures. The ability of CCM to predict DPN broadens the already impressive diagnostic capabilities of this novel ophthalmic marker.

Atmospheric visibility and PM10 as indicators of new particle formation in an urban environment

It is well-known that new particle formation (NPF) in the atmosphere is inhibited by pre-existing particles in the air that act as condensation sinks to decrease the concentration and, thus, the supersaturation of precursor gases. In this study, we investigate the effects of two parameters - atmospheric visibility, expressed as the particle back-scatter coefficient (BSP), and PM10 particulate mass concentration, on the occurrences of NPF events in an urban environment where the majority of precursor gases originate from motor vehicle and industrial sources. This is the first attempt to derive direct relationships between each of these two parameters and the occurrence of NPF. NPF events were identified from data obtained with a neutral cluster and air ion spectrometer over 245 days within a calendar year. Bayesian logistic regression was used to determine the probability of observing NPF as functions of BSP and PM10. We show that the BSP at 08 h on a given day is a reliable indicator of an NPF event later that day. The posterior median probability of observing an NPF event was greater than 0.5 (95%) when the BSP at 08 h was less than 6.8 Mm-1.

Hypoxic corneal changes following 8 hours of scleral contact lens wear

- Purpose To examine the change in corneal thickness and posterior curvature following 8 hours of miniscleral contact lens wear. - Methods Scheimpflug imaging (Pentacam HR, Oculus) was captured before, and immediately following, 8 hours of miniscleral contact lens wear for 15 young (mean age 22 ± 3 years), healthy participants with normal corneae. Natural diurnal variations were considered by measuring baseline corneal changes obtained on a separate control day without contact lens wear. - Results Over the central 6 mm of the cornea, a small, but highly statistically significant amount of edema was observed following 8 hours of miniscleral lens wear, after accounting for normal diurnal fluctuations (mean ± standard deviation percentage swelling 1.70 ± 0.98%, p < 0.0001). Posterior corneal topography remained stable following lens wear (-0.01 ± 0.07 mm steepening over the central 6 mm, p = 0.60). The magnitude of posterior corneal topographical changes following lens wear did not correlate with the extent of lens-related corneal edema (r = -0.16, p = 0.57). Similarly, the initial central corneal vault (maximum post-lens tear layer depth) was not associated with corneal swelling following lens removal (r = 0.27, p = 0.33). - Conclusions While a small amount of corneal swelling was induced following 8 hours of miniscleral lens wear (on average <2%), modern high Dk miniscleral contact lenses that vault the cornea do not induce clinically significant corneal edema or hypoxic related posterior corneal curvature changes during short-term wear. Longer-term studies of compromised eyes (e.g. corneal ectasia) are still required to inform the optimum lens and fitting characteristics for safe scleral lens wear to minimize corneal hypoxia.

Dynamic Bayesian network inferencing for non-homogeneous complex systems

Dynamic Bayesian Networks (DBNs) provide a versatile platform for predicting and analysing the behaviour of complex systems. As such, they are well suited to the prediction of complex ecosystem population trajectories under anthropogenic disturbances such as the dredging of marine seagrass ecosystems. However, DBNs assume a homogeneous Markov chain whereas a key characteristics of complex ecosystems is the presence of feedback loops, path dependencies and regime changes whereby the behaviour of the system can vary based on past states. This paper develops a method based on the small world structure of complex systems networks to modularise a non-homogeneous DBN and enable the computation of posterior marginal probabilities given evidence in forwards inference. It also provides an approach for an approximate solution for backwards inference as convergence is not guaranteed for a path dependent system. When applied to the seagrass dredging problem, the incorporation of path dependency can implement conditional absorption and allows release from the zero state in line with environmental and ecological observations. As dredging has a marked global impact on seagrass and other marine ecosystems of high environmental and economic value, using such a complex systems model to develop practical ways to meet the needs of conservation and industry through enhancing resistance and/or recovery is of paramount importance.

Biometry of eyes in type 1 diabetes

This is a comprehensive study of a large range of biometric and optical parameters in people with type 1 diabetes. The parameters of 74 people with type 1 diabetes and an age matched control group were assessed. Most of the people with diabetes had low levels of neuropathy, retinopathy and nephropathy. Marginal or no significant differences were found between groups for corneal shape, corneal thickness, pupil size, and pupil decentrations. Relative to the control group, the diabetes group demonstrated smaller anterior chamber depths, more curved lenses, greater lens thickness and lower lens equivalent refractive index. While the optics of diabetic eyes make them appear as older eyes than those of people of the same age without diabetes, the differences did not increase significantly with age. Age-related changes in the optics of the eyes of people with diabetes need not be accelerated if the diabetes is well controlled.

Bayesian estimation of the dynamics of pandemic (H1N1) 2009 influenza transmission in Queensland: A space–time SIR-based model

Background A pandemic strain of influenza A spread rapidly around the world in 2009, now referred to as pandemic (H1N1) 2009. This study aimed to examine the spatiotemporal variation in the transmission rate of pandemic (H1N1) 2009 associated with changes in local socio-environmental conditions from May 7–December 31, 2009, at a postal area level in Queensland, Australia. Method We used the data on laboratory-confirmed H1N1 cases to examine the spatiotemporal dynamics of transmission using a flexible Bayesian, space–time, Susceptible-Infected-Recovered (SIR) modelling approach. The model incorporated parameters describing spatiotemporal variation in H1N1 infection and local socio-environmental factors. Results The weekly transmission rate of pandemic (H1N1) 2009 was negatively associated with the weekly area-mean maximum temperature at a lag of 1 week (LMXT) (posterior mean: −0.341; 95% credible interval (CI): −0.370–−0.311) and the socio-economic index for area (SEIFA) (posterior mean: −0.003; 95% CI: −0.004–−0.001), and was positively associated with the product of LMXT and the weekly area-mean vapour pressure at a lag of 1 week (LVAP) (posterior mean: 0.008; 95% CI: 0.007–0.009). There was substantial spatiotemporal variation in transmission rate of pandemic (H1N1) 2009 across Queensland over the epidemic period. High random effects of estimated transmission rates were apparent in remote areas and some postal areas with higher proportion of indigenous populations and smaller overall populations. Conclusions Local SEIFA and local atmospheric conditions were associated with the transmission rate of pandemic (H1N1) 2009. The more populated regions displayed consistent and synchronized epidemics with low average transmission rates. The less populated regions had high average transmission rates with more variations during the H1N1 epidemic period.

Assessing the combined effect of asbestos exposure and smoking on lung cancer: a Bayesian approach

We review the literature on the combined association between lung cancer and two environmental exposures, asbestos exposure and smoking, and explore a Bayesian approach to assess evidence of interaction between the exposures. The meta-analysis combines separate indices of additive and multiplicative relationships and multivariate relative risk estimates. By making inferences on posterior probabilities we can explore both the form and strength of interaction. This analysis may be more informative than providing evidence to support one relation over another on the basis of statistical significance. Overall, we find evidence for a more than additive and less than multiplicative relation.

Morphological changes in the conjunctiva, episclera and sclera following short-term miniscleral contact lens wear in rigid lens neophytes

PURPOSE To quantify the influence of short-term wear of miniscleral contact lenses on the morphology of the corneo-scleral limbus, the conjunctiva, episclera and sclera. METHODS OCT images of the anterior eye were captured before, immediately following 3h of wear and then 3h after removal of a miniscleral contact lens for 10 young (27±5 years) healthy participants (neophyte rigid lens wearers). The region of analysis encompassed 1mm anterior, to 3.5mm posterior to the scleral spur. Natural diurnal variations in thickness were measured on a separate day and compensated for in subsequent analyses. RESULTS Following 3h of lens wear, statistically significant tissue thinning was observed across all quadrants, with a mean decrease in thickness of -24.1±3.6μm (p<0.001), which diminished, but did not return to baseline 3h after lens removal (-16.9±1.9μm, p<0.001). The largest tissue compression was observed in the superior quadrant (-49.9±8.5μm, p<0.01) and in the annular zone 1.5mm from the scleral spur (-48.2±5.7μm), corresponding to the approximate edge of the lens landing zone. Compression of the conjunctiva/episclera accounted for about 70% of the changes. CONCLUSIONS Optimal fitting miniscleral contact lenses worn for three hours resulted in significant tissue compression in young healthy eyes, with the greatest thinning observed superiorly, potentially due to the additional force of the eyelid, with a partial recovery of compression 3h after lens removal. Most of the morphological changes occur in the conjunctiva/episclera layers.

Evidence from molecules and morphology expands Podonomopsis Brundin (Diptera: Chironomidae: Podonominae) to include ‘genus Chile'

The informal taxon ‘genus Chile’ of Brundin, based solely on pupal exuviae of a podonomine Chironomidae, has remained inadequately known for half a century. New collections reveal life associations, and provide molecular data to hypothesise a precise phylogenetic placement in the austral Podonominae. A densely sampled molecular phylogeny based on two nuclear and one mitochondrial DNA markers shows ‘genus Chile’ to be the sister group to Podonomopsis Brundin, 1966. Within Podonomopsis a clade of South American species is sister to all Australian species. We discuss how to rank such a sister group taxon and treat ‘genus Chile’ as a new subgenus Araucanopsis, subg. nov. with the new species, Podonomopsis (Araucanopsis) avelasse, sp. nov. from Chile and Argentina as genotype of the monotypic subgenus. We describe P. (A.) avelasse in all stages and provide an expanded diagnosis and description of Podonomopsis to include Araucanopsis. A dated biogeographic hypothesis (chronogram) infers the most recent common ancestor (tmcra) of expanded Podonomopsis at 95 million years ago (Mya) (68–122 Mya 95% highest posterior density), ‘core’ Podonomopsis at 83 Mya (58–108) and Australian Podonomopsis at 65 Mya (44–87). All dates are before the South America–Australia geological separation through Antarctica, supporting previous conclusions that the taxon distribution is ‘Gondwanan’ in origin. Podonomopsis, even as expanded here, remains unknown from New Zealand or elsewhere on extant Zealandia.