Statistical model-based segmentation of the proximal femur in digital antero-posterior (AP) pelvic radiographs.

PURPOSE    Segmentation of the proximal femur in digital antero-posterior (AP) pelvic radiographs is required to create a three-dimensional model of the hip joint for use in planning and treatment. However, manually extracting the femoral contour is tedious and prone to subjective bias, while automatic segmentation must accommodate poor image quality, anatomical structure overlap, and femur deformity. A new method was developed for femur segmentation in AP pelvic radiographs. METHODS    Using manual annotations on 100 AP pelvic radiographs, a statistical shape model (SSM) and a statistical appearance model (SAM) of the femur contour were constructed. The SSM and SAM were used to segment new AP pelvic radiographs with a three-stage approach. At initialization, the mean SSM model is coarsely registered to the femur in the AP radiograph through a scaled rigid registration. Mahalanobis distance defined on the SAM is employed as the search criteria for each annotated suggested landmark location. Dynamic programming was used to eliminate ambiguities. After all landmarks are assigned, a regularized non-rigid registration method deforms the current mean shape of SSM to produce a new segmentation of proximal femur. The second and third stages are iteratively executed to convergence. RESULTS    A set of 100 clinical AP pelvic radiographs (not used for training) were evaluated. The mean segmentation error was [Formula: see text], requiring [Formula: see text] s per case when implemented with Matlab. The influence of the initialization on segmentation results was tested by six clinicians, demonstrating no significance difference. CONCLUSIONS    A fast, robust and accurate method for femur segmentation in digital AP pelvic radiographs was developed by combining SSM and SAM with dynamic programming. This method can be extended to segmentation of other bony structures such as the pelvis.

Characterization of the 9L gliosarcoma implanted in the Fischer rat: an orthotopic model for a grade IV brain tumor.

Among rodent models for brain tumors, the 9L gliosarcoma is one of the most widely used. Our 9L-European Synchrotron Radiation Facility (ESRF) model was developed from cells acquired at the Brookhaven National Laboratory (NY, USA) in 1997 and implanted in the right caudate nucleus of syngeneic Fisher rats. It has been largely used by the user community of the ESRF during the last decade, for imaging, radiotherapy, and chemotherapy, including innovative treatments based on particular irradiation techniques and/or use of new drugs. This work presents a detailed study of its characteristics, assessed by magnetic resonance imaging (MRI), histology, immunohistochemistry, and cytogenetic analysis. The data used for this work were from rats sampled in six experiments carried out over a 3-year period in our lab (total number of rats = 142). The 9L-ESRF tumors were induced by a stereotactic inoculation of 10(4) 9L cells in the right caudate nucleus of the brain. The assessment of vascular parameters was performed by MRI (blood volume fraction and vascular size index) and by immunostaining of vessels (rat endothelial cell antigen-1 and type IV collagen). Immunohistochemistry and regular histology were used to describe features such as tumor cell infiltration, necrosis area, nuclear pleomorphism, cellularity, mitotic characteristics, leukocytic infiltration, proliferation, and inflammation. Moreover, for each of the six experiments, the survival of the animals was assessed and related to the tumor growth observed by MRI or histology. Additionally, the cytogenetic status of the 9L cells used at ESRF lab was investigated by comparative genomics hybridization analysis. Finally, the response of the 9L-ESRF tumor to radiotherapy was estimated by plotting the survival curves after irradiation. The median survival time of 9L-ESRF tumor-bearing rats was highly reproducible (19-20 days). The 9L-ESRF tumors presented a quasi-exponential growth, were highly vascularized with a high cellular density and a high proliferative index, accompanied by signs of inflammatory responses. We also report an infiltrative pattern which is poorly observed on conventional 9 L tumor. The 9L-ESRF cells presented some cytogenetic specificities such as altered regions including CDK4, CDKN2A, CDKN2B, and MDM2 genes. Finally, the lifespan of 9L-ESRF tumor-bearing rats was enhanced up to 28, 35, and 45 days for single doses of 10, 20, and 2 × 20 Gy, respectively. First, this report describes an animal model that is used worldwide. Second, we describe few features typical of our model if compared to other 9L models worldwide. Altogether, the 9L-ESRF tumor model presents characteristics close to the human high-grade gliomas such as high proliferative capability, high vascularization and a high infiltrative pattern. Its response to radiotherapy demonstrates its potential as a tool for innovative radiotherapy protocols.

Non-exponential return time distributions for vorticity extremes explained by fractional Poisson processes

Serial correlation of extreme midlatitude cyclones observed at the storm track exits is explained by deviations from a Poisson process. To model these deviations, we apply fractional Poisson processes (FPPs) to extreme midlatitude cyclones, which are defined by the 850 hPa relative vorticity of the ERA interim reanalysis during boreal winter (DJF) and summer (JJA) seasons. Extremes are defined by a 99% quantile threshold in the grid-point time series. In general, FPPs are based on long-term memory and lead to non-exponential return time distributions. The return times are described by a Weibull distribution to approximate the Mittag–Leffler function in the FPPs. The Weibull shape parameter yields a dispersion parameter that agrees with results found for midlatitude cyclones. The memory of the FPP, which is determined by detrended fluctuation analysis, provides an independent estimate for the shape parameter. Thus, the analysis exhibits a concise framework of the deviation from Poisson statistics (by a dispersion parameter), non-exponential return times and memory (correlation) on the basis of a single parameter. The results have potential implications for the predictability of extreme cyclones.

Lepton Flavor Violation in the Standard Model with general Dimension-Six Operators

We study lepton flavor observables in the Standard Model (SM) extended with all dimension-6 operators which are invariant under the SM gauge group. We calculate the complete one-loop predictions to the radiative lepton decays μ → eγ, τ → μγ and τ → eγ as well as to the closely related anomalous magnetic moments and electric dipole moments of charged leptons, taking into account all dimension-6 operators which can generate lepton flavor violation. Also the 3-body flavor violating charged lepton decays τ ± → μ ± μ + μ −, τ ± → e ± e + e −, τ ± → e ± μ + μ −, τ ± → μ ± e + e −, τ ± → e ∓ μ ± μ ±, τ ± → μ ∓ e ± e ± and μ ± → e ± e + e − and the Z 0 decays Z 0 → ℓ+iℓ−j are considered, taking into account all tree-level contributions.

Lesion location predicts transient and extended risk of aspiration after supratentorial ischemic stroke

BACKGROUND AND PURPOSE To assess the association of lesion location and risk of aspiration and to establish predictors of transient versus extended risk of aspiration after supratentorial ischemic stroke. METHODS Atlas-based localization analysis was performed in consecutive patients with MRI-proven first-time acute supratentorial ischemic stroke. Standardized swallowing assessment was carried out within 8±18 hours and 7.8±1.2 days after admission. RESULTS In a prospective, longitudinal analysis, 34 of 94 patients (36%) were classified as having acute risk of aspiration, which was extended (≥7 days) or transient (<7 days) in 17 cases. There were no between-group differences in age, sex, cause of stroke, risk factors, prestroke disability, lesion side, or the degree of age-related white-matter changes. Correcting for stroke volume and National Institutes of Health Stroke Scale with a multiple logistic regression model, significant adjusted odds ratios in favor of acute risk of aspiration were demonstrated for the internal capsule (adjusted odds ratio, 6.2; P<0.002) and the insular cortex (adjusted odds ratio, 4.8; P<0.003). In a multivariate model of extended versus transient risk of aspiration, combined lesions of the frontal operculum and insular cortex was the only significant independent predictor of poor recovery (adjusted odds ratio, 33.8; P<0.008). CONCLUSIONS Lesions of the insular cortex and the internal capsule are significantly associated with acute risk of aspiration after stroke. Combined ischemic infarctions of the frontal operculum and the insular cortex are likely to cause extended risk of aspiration in stroke patients, whereas risk of aspiration tends to be transient in subcortical stroke.

A new model construction by making a detour via intuitionistic theories I: Operational set theory without choice is Π1-equivalent to KP

We introduce a version of operational set theory, OST−, without a choice operation, which has a machinery for Δ0Δ0 separation based on truth functions and the separation operator, and a new kind of applicative set theory, so-called weak explicit set theory WEST, based on Gödel operations. We show that both the theories and Kripke–Platek set theory KPKP with infinity are pairwise Π1Π1 equivalent. We also show analogous assertions for subtheories with ∈-induction restricted in various ways and for supertheories extended by powerset, beta, limit and Mahlo operations. Whereas the upper bound is given by a refinement of inductive definition in KPKP, the lower bound is by a combination, in a specific way, of realisability, (intuitionistic) forcing and negative interpretations. Thus, despite interpretability between classical theories, we make “a detour via intuitionistic theories”. The combined interpretation, seen as a model construction in the sense of Visser's miniature model theory, is a new way of construction for classical theories and could be said the third kind of model construction ever used which is non-trivial on the logical connective level, after generic extension à la Cohen and Krivine's classical realisability model.

The Utilization of Classical Spin Monte Carlo Methods to Simulate the Magnetic Behavior of Extended Three-Dimensional Cubic Networks Incorporating M(II) Ions with an S = 5/2 Ground State Spin

The numerical simulations of the magnetic properties of extended three-dimensional networks containing M(II) ions with an S = 5/2 ground-state spin have been carried out within the framework of the isotropic Heisenberg model. Analytical expressions fitting the numerical simulations for the primitive cubic, diamond, together with (10−3) cubic networks have all been derived. With these empirical formulas in hands, we can now extract the interaction between the magnetic ions from the experimental data for these networks. In the case of the primitive cubic network, these expressions are directly compared with those from the high-temperature expansions of the partition function. A fit of the experimental data for three complexes, namely [(N(CH3)4][Mn(N3)] 1, [Mn(CN4)]n 2, and [FeII(bipy)3][MnII2(ox)3] 3, has been carried out. The best fits were those obtained using the following parameters, J = −3.5 cm-1, g = 2.01 (1); J = −8.3 cm-1, g = 1.95 (2); and J = −2.0 cm-1, g = 1.95 (3).

CODE’s new solar radiation pressure model for GNSS orbit determination

The Empirical CODE Orbit Model (ECOM) of the Center for Orbit Determination in Europe (CODE), which was developed in the early 1990s, is widely used in the International GNSS Service (IGS) community. For a rather long time, spurious spectral lines are known to exist in geophysical parameters, in particular in the Earth Rotation Parameters (ERPs) and in the estimated geocenter coordinates, which could recently be attributed to the ECOM. These effects grew creepingly with the increasing influence of the GLONASS system in recent years in the CODE analysis, which is based on a rigorous combination of GPS and GLONASS since May 2003. In a first step we show that the problems associated with the ECOM are to the largest extent caused by the GLONASS, which was reaching full deployment by the end of 2011. GPS-only, GLONASS-only, and combined GPS/GLONASS solutions using the observations in the years 2009–2011 of a global network of 92 combined GPS/GLONASS receivers were analyzed for this purpose. In a second step we review direct solar radiation pressure (SRP) models for GNSS satellites. We demonstrate that only even-order short-period harmonic perturbations acting along the direction Sun-satellite occur for GPS and GLONASS satellites, and only odd-order perturbations acting along the direction perpendicular to both, the vector Sun-satellite and the spacecraft’s solar panel axis. Based on this insight we assess in the third step the performance of four candidate orbit models for the future ECOM. The geocenter coordinates, the ERP differences w. r. t. the IERS 08 C04 series of ERPs, the misclosures for the midnight epochs of the daily orbital arcs, and scale parameters of Helmert transformations for station coordinates serve as quality criteria. The old and updated ECOM are validated in addition with satellite laser ranging (SLR) observations and by comparing the orbits to those of the IGS and other analysis centers. Based on all tests, we present a new extended ECOM which substantially reduces the spurious signals in the geocenter coordinate z (by about a factor of 2–6), reduces the orbit misclosures at the day boundaries by about 10 %, slightly improves the consistency of the estimated ERPs with those of the IERS 08 C04 Earth rotation series, and substantially reduces the systematics in the SLR validation of the GNSS orbits.

Random Sampling with Interspike-Intervals of the Exponential Integrate and Fire Neuron: A Computational Interpretation of UP-States

Oscillations between high and low values of the membrane potential (UP and DOWN states respectively) are an ubiquitous feature of cortical neurons during slow wave sleep and anesthesia. Nevertheless, a surprisingly small number of quantitative studies have been conducted only that deal with this phenomenon’s implications for computation. Here we present a novel theory that explains on a detailed mathematical level the computational benefits of UP states. The theory is based on random sampling by means of interspike intervals (ISIs) of the exponential integrate and fire (EIF) model neuron, such that each spike is considered a sample, whose analog value corresponds to the spike’s preceding ISI. As we show, the EIF’s exponential sodium current, that kicks in when balancing a noisy membrane potential around values close to the firing threshold, leads to a particularly simple, approximative relationship between the neuron’s ISI distribution and input current. Approximation quality depends on the frequency spectrum of the current and is improved upon increasing the voltage baseline towards threshold. Thus, the conceptually simpler leaky integrate and fire neuron that is missing such an additional current boost performs consistently worse than the EIF and does not improve when voltage baseline is increased. For the EIF in contrast, the presented mechanism is particularly effective in the high-conductance regime, which is a hallmark feature of UP-states. Our theoretical results are confirmed by accompanying simulations, which were conducted for input currents of varying spectral composition. Moreover, we provide analytical estimations of the range of ISI distributions the EIF neuron can sample from at a given approximation level. Such samples may be considered by any algorithmic procedure that is based on random sampling, such as Markov Chain Monte Carlo or message-passing methods. Finally, we explain how spike-based random sampling relates to existing computational theories about UP states during slow wave sleep and present possible extensions of the model in the context of spike-frequency adaptation.

Transient Earth system responses to cumulative carbon dioxide emissions: linearities, uncertainties, and probabilities in an observation-constrained model ensemble

Information on the relationship between cumulative fossil CO2 emissions and multiple climate targets is essential to design emission mitigation and climate adaptation strategies. In this study, the transient response of a climate or environmental variable per trillion tonnes of CO2 emissions, termed TRE, is quantified for a set of impact-relevant climate variables and from a large set of multi-forcing scenarios extended to year 2300 towards stabilization. An  ∼ 1000-member ensemble of the Bern3D-LPJ carbon–climate model is applied and model outcomes are constrained by 26 physical and biogeochemical observational data sets in a Bayesian, Monte Carlo-type framework. Uncertainties in TRE estimates include both scenario uncertainty and model response uncertainty. Cumulative fossil emissions of 1000 Gt C result in a global mean surface air temperature change of 1.9 °C (68 % confidence interval (c.i.): 1.3 to 2.7 °C), a decrease in surface ocean pH of 0.19 (0.18 to 0.22), and a steric sea level rise of 20 cm (13 to 27 cm until 2300). Linearity between cumulative emissions and transient response is high for pH and reasonably high for surface air and sea surface temperatures, but less pronounced for changes in Atlantic meridional overturning, Southern Ocean and tropical surface water saturation with respect to biogenic structures of calcium carbonate, and carbon stocks in soils. The constrained model ensemble is also applied to determine the response to a pulse-like emission and in idealized CO2-only simulations. The transient climate response is constrained, primarily by long-term ocean heat observations, to 1.7 °C (68 % c.i.: 1.3 to 2.2 °C) and the equilibrium climate sensitivity to 2.9 °C (2.0 to 4.2 °C). This is consistent with results by CMIP5 models but inconsistent with recent studies that relied on short-term air temperature data affected by natural climate variability.

Exposure assessment of extended-spectrum beta-lactamases/AmpC beta-lactamases-producing Escherichia coli in meat in Denmark

INTRODUCTION Extended-spectrum beta-lactamases (ESBL) and AmpC beta-lactamases (AmpC) are of concern for veterinary and public health because of their ability to cause treatment failure due to antimicrobial resistance in Enterobacteriaceae. The main objective was to assess the relative contribution (RC) of different types of meat to the exposure of consumers to ESBL/AmpC and their potential importance for human infections in Denmark. MATERIAL AND METHODS The prevalence of each genotype of ESBL/AmpC-producing E. coli in imported and nationally produced broiler meat, pork and beef was weighted by the meat consumption patterns. Data originated from the Danish surveillance program for antibiotic use and antibiotic resistance (DANMAP) from 2009 to 2011. DANMAP also provided data about human ESBL/AmpC cases in 2011, which were used to assess a possible genotype overlap. Uncertainty about the occurrence of ESBL/AmpC-producing E. coli in meat was assessed by inspecting beta distributions given the available data of the genotypes in each type of meat. RESULTS AND DISCUSSION Broiler meat represented the largest part (83.8%) of the estimated ESBL/AmpC-contaminated pool of meat compared to pork (12.5%) and beef (3.7%). CMY-2 was the genotype with the highest RC to human exposure (58.3%). However, this genotype is rarely found in human infections in Denmark. CONCLUSION The overlap between ESBL/AmpC genotypes in meat and human E. coli infections was limited. This suggests that meat might constitute a less important source of ESBL/AmpC exposure to humans in Denmark than previously thought - maybe because the use of cephalosporins is restricted in cattle and banned in poultry and pigs. Nonetheless, more detailed surveillance data are required to determine the contribution of meat compared to other sources, such as travelling, pets, water resources, community and hospitals in the pursuit of a full source attribution model.

Establishment of a confluent monolayer model with human primary trophoblast cells: novel insights into placental glucose transport

STUDY HYPOTHESIS Using optimized conditions, primary trophoblast cells isolated from human term placenta can develop a confluent monolayer in vitro, which morphologically and functionally resembles the microvilli structure found in vivo. STUDY FINDING We report the successful establishment of a confluent human primary trophoblast monolayer using pre-coated polycarbonate inserts, where the integrity and functionality was validated by cell morphology, biophysical features, cellular marker expression and secretion, and asymmetric glucose transport. WHAT IS KNOWN ALREADY Human trophoblast cells form the initial barrier between maternal and fetal blood to regulate materno-fetal exchange processes. Although the method for isolating pure human cytotrophoblast cells was developed almost 30 years ago, a functional in vitro model with primary trophoblasts forming a confluent monolayer is still lacking. STUDY DESIGN, SAMPLES/MATERIALS, METHODS Human term cytotrophoblasts were isolated by enzymatic digestion and density gradient separation. The purity of the primary cells was evaluated by flow cytometry using the trophoblast-specific marker cytokeratin 7, and vimentin as an indicator for potentially contaminating cells. We screened different coating matrices for high cell viability to optimize the growth conditions for primary trophoblasts on polycarbonate inserts. During culture, cell confluency and polarity were monitored daily by determining transepithelial electrical resistance (TEER) and permeability properties of florescent dyes. The time course of syncytia-related gene expression and hCG secretion during syncytialization were assessed by quantitative RT-PCR and enzyme-linked immunosorbent assay, respectively. The morphology of cultured trophoblasts after 5 days was determined by light microscopy, scanning electron microscopy (SEM) and transmission electron microscopy (TEM). Membrane makers were visualized using confocal microscopy. Additionally, glucose transport studies were performed on the polarized trophoblasts in the same system. MAIN RESULTS AND THE ROLE OF CHANCE During 5-day culture, the highly pure trophoblasts were cultured on inserts coated with reconstituted basement membrane matrix . They exhibited a confluent polarized monolayer, with a modest TEER and a size-dependent apparent permeability coefficient (Papp) to fluorescently labeled compounds (MW ∼400-70 000 Da). The syncytialization progress was characterized by gradually increasing mRNA levels of fusogen genes and elevating hCG secretion. SEM analyses confirmed a confluent trophoblast layer with numerous microvilli, and TEM revealed a monolayer with tight junctions. Immunocytochemistry on the confluent trophoblasts showed positivity for the cell-cell adhesion molecule E-cadherin, the tight junction protein 1 (ZO-1) and the membrane proteins ATP-binding cassette transporter A1 (ABCA1) and glucose transporter 1 (GLUT1). Applying this model to study the bidirectional transport of a non-metabolizable glucose derivative indicated a carrier-mediated placental glucose transport mechanism with asymmetric kinetics. LIMITATIONS, REASONS FOR CAUTION The current study is only focused on primary trophoblast cells isolated from healthy placentas delivered at term. It remains to be evaluated whether this system can be extended to pathological trophoblasts isolated from diverse gestational diseases. WIDER IMPLICATIONS OF THE FINDINGS These findings confirmed the physiological properties of the newly developed human trophoblast barrier, which can be applied to study the exchange of endobiotics and xenobiotics between the maternal and fetal compartment, as well as intracellular metabolism, paracellular contributions and regulatory mechanisms influencing the vectorial transport of molecules. LARGE-SCALE DATA Not applicable. STUDY FUNDING AND COMPETING INTERESTS This study was supported by the Swiss National Center of Competence in Research, NCCR TransCure, University of Bern, Switzerland, and the Swiss National Science Foundation (grant no. 310030_149958, C.A.). All authors declare that their participation in the study did not involve factual or potential conflicts of interests.

Use of Extended Characteristics of Locomotion and Feeding Behavior for Automated Identification of Lame Dairy Cows.

This study was carried out to detect differences in locomotion and feeding behavior in lame (group L; n = 41; gait score ≥ 2.5) and non-lame (group C; n = 12; gait score ≤ 2) multiparous Holstein cows in a cross-sectional study design. A model for automatic lameness detection was created, using data from accelerometers attached to the hind limbs and noseband sensors attached to the head. Each cow's gait was videotaped and scored on a 5-point scale before and after a period of 3 consecutive days of behavioral data recording. The mean value of 3 independent experienced observers was taken as a definite gait score and considered to be the gold standard. For statistical analysis, data from the noseband sensor and one of two accelerometers per cow (randomly selected) of 2 out of 3 randomly selected days was used. For comparison between group L and group C, the T-test, the Aspin-Welch Test and the Wilcoxon Test were used. The sensitivity and specificity for lameness detection was determined with logistic regression and ROC-analysis. Group L compared to group C had significantly lower eating and ruminating time, fewer eating chews, ruminating chews and ruminating boluses, longer lying time and lying bout duration, lower standing time, fewer standing and walking bouts, fewer, slower and shorter strides and a lower walking speed. The model considering the number of standing bouts and walking speed was the best predictor of cows being lame with a sensitivity of 90.2% and specificity of 91.7%. Sensitivity and specificity of the lameness detection model were considered to be very high, even without the use of halter data. It was concluded that under the conditions of the study farm, accelerometer data were suitable for accurately distinguishing between lame and non-lame dairy cows, even in cases of slight lameness with a gait score of 2.5.