The Influence of Osteoblast Differentiation Stage on Bone Formation in Autogenously Implanted Cell-Based Poly(Lactide-Co-Glycolide) and Calcium Phosphate Constructs

We tested the hypothesis that the osteoblast differentiation status of bone marrow stem cells (BMSCs) combined with a three-dimensional (3D) structure modulates bone formation when autogenously implanted. Rat BMSCs were aspirated, expanded, and seeded into a 3D composite of poly(lactide-co-glycolide) and calcium phosphate (PLGA/CaP) to produce a hybrid biomaterial. Calvarial defects were implanted with (1) scaffold without cells (SC/NC), (2) scaffold and BMSCs (SC + BMSC), (3) scaffold and osteoblasts differentiated for 7 days (SC + OB7), and (4) for 14 days (SC + OB14). After 4 weeks, there was more bone formation in groups combining scaffold and cells, SC + BMSC and SC + OB7. A nonsignificant higher amount of bone formation was observed on SC + OB14 compared with SC/NC. Additionally, more blood vessels were counted within all hybrid biomaterials, without differences among them, than into SC/NC. These findings provide evidences that the cell differentiation status affects in vivo bone formation in autogenously implanted cell-based constructs. Undifferentiated BMSCs or osteoblasts in early stage of differentiation combined with PLGA/CaP scaffold favored bone formation compared with plain scaffold and that one associated with more mature osteoblasts.

Potential of Bioactive Glass Particles of Different Size Ranges to Affect Bone Formation in Interproximal Periodontal Defects in Dogs

Background: The aim of this study was to compare the potential of bioactive glass particles of different size ranges to affect bone formation in periodontal defects, using the guided tissue regeneration model in dogs. Methods: In six dogs, 2-wall intrabony periodontal defects were surgically created and chronified on the mesial surfaces of mandibular third premolars and first molars bilaterally. After 1 month, each defect was randomly assigned to treatment with bioabsorbable membrane in association with bioactive glass with particle sizes between 300 and 355 mu m (group 1) or between 90 and 710 mu m (group 2), membrane alone (group 3), or negative control (group 4). The dogs were sacrificed 12 weeks after surgeries, and histomorphometric measurements were made of the areas of newly formed bone, new mineralized bone, and bioactive glass particle remnants. Results: With regard to the area of bioactive glass particle remnants, there was a statistically significant difference between groups 1 and 2, favoring group 1. There were greater areas of mineralized bone in groups 1 and 2 compared to groups 3 and 4 (P<0.05). Conclusion: The bioactive glass particles of small size range underwent faster resorption and substitution by new bone than the larger particles, and the use of bioactive glass particles favored the formation of mineralized bone. J Periodontol 2009;80:808-815.

Parameter recovery for a skew-normal IRT model under a Bayesian approach: hierarchical framework, prior and kernel sensitivity and sample size

Item response theory (IRT) comprises a set of statistical models which are useful in many fields, especially when there is an interest in studying latent variables (or latent traits). Usually such latent traits are assumed to be random variables and a convenient distribution is assigned to them. A very common choice for such a distribution has been the standard normal. Recently, Azevedo et al. [Bayesian inference for a skew-normal IRT model under the centred parameterization, Comput. Stat. Data Anal. 55 (2011), pp. 353-365] proposed a skew-normal distribution under the centred parameterization (SNCP) as had been studied in [R. B. Arellano-Valle and A. Azzalini, The centred parametrization for the multivariate skew-normal distribution, J. Multivariate Anal. 99(7) (2008), pp. 1362-1382], to model the latent trait distribution. This approach allows one to represent any asymmetric behaviour concerning the latent trait distribution. Also, they developed a Metropolis-Hastings within the Gibbs sampling (MHWGS) algorithm based on the density of the SNCP. They showed that the algorithm recovers all parameters properly. Their results indicated that, in the presence of asymmetry, the proposed model and the estimation algorithm perform better than the usual model and estimation methods. Our main goal in this paper is to propose another type of MHWGS algorithm based on a stochastic representation (hierarchical structure) of the SNCP studied in [N. Henze, A probabilistic representation of the skew-normal distribution, Scand. J. Statist. 13 (1986), pp. 271-275]. Our algorithm has only one Metropolis-Hastings step, in opposition to the algorithm developed by Azevedo et al., which has two such steps. This not only makes the implementation easier but also reduces the number of proposal densities to be used, which can be a problem in the implementation of MHWGS algorithms, as can be seen in [R.J. Patz and B.W. Junker, A straightforward approach to Markov Chain Monte Carlo methods for item response models, J. Educ. Behav. Stat. 24(2) (1999), pp. 146-178; R. J. Patz and B. W. Junker, The applications and extensions of MCMC in IRT: Multiple item types, missing data, and rated responses, J. Educ. Behav. Stat. 24(4) (1999), pp. 342-366; A. Gelman, G.O. Roberts, and W.R. Gilks, Efficient Metropolis jumping rules, Bayesian Stat. 5 (1996), pp. 599-607]. Moreover, we consider a modified beta prior (which generalizes the one considered in [3]) and a Jeffreys prior for the asymmetry parameter. Furthermore, we study the sensitivity of such priors as well as the use of different kernel densities for this parameter. Finally, we assess the impact of the number of examinees, number of items and the asymmetry level on the parameter recovery. Results of the simulation study indicated that our approach performed equally as well as that in [3], in terms of parameter recovery, mainly using the Jeffreys prior. Also, they indicated that the asymmetry level has the highest impact on parameter recovery, even though it is relatively small. A real data analysis is considered jointly with the development of model fitting assessment tools. The results are compared with the ones obtained by Azevedo et al. The results indicate that using the hierarchical approach allows us to implement MCMC algorithms more easily, it facilitates diagnosis of the convergence and also it can be very useful to fit more complex skew IRT models.

Exact field-driven interface dynamics in the two-dimensional stochastic Ising model with helicoidal boundary conditions

We investigate the interface dynamics of the two-dimensional stochastic Ising model in an external field under helicoidal boundary conditions. At sufficiently low temperatures and fields, the dynamics of the interface is described by an exactly solvable high-spin asymmetric quantum Hamiltonian that is the infinitesimal generator of the zero range process. Generally, the critical dynamics of the interface fluctuations is in the Kardar-Parisi-Zhang universality class of critical behavior. We remark that a whole family of RSOS interface models similar to the Ising interface model investigated here can be described by exactly solvable restricted high-spin quantum XXZ-type Hamiltonians. (C) 2012 Elsevier B.V. All rights reserved.

Connections between the Sznajd model with general confidence rules and graph theory

The Sznajd model is a sociophysics model that is used to model opinion propagation and consensus formation in societies. Its main feature is that its rules favor bigger groups of agreeing people. In a previous work, we generalized the bounded confidence rule in order to model biases and prejudices in discrete opinion models. In that work, we applied this modification to the Sznajd model and presented some preliminary results. The present work extends what we did in that paper. We present results linking many of the properties of the mean-field fixed points, with only a few qualitative aspects of the confidence rule (the biases and prejudices modeled), finding an interesting connection with graph theory problems. More precisely, we link the existence of fixed points with the notion of strongly connected graphs and the stability of fixed points with the problem of finding the maximal independent sets of a graph. We state these results and present comparisons between the mean field and simulations in Barabasi-Albert networks, followed by the main mathematical ideas and appendices with the rigorous proofs of our claims and some graph theory concepts, together with examples. We also show that there is no qualitative difference in the mean-field results if we require that a group of size q > 2, instead of a pair, of agreeing agents be formed before they attempt to convince other sites (for the mean field, this would coincide with the q-voter model).

A Simple Artificial Life Model Explains Irrational Behavior in Human Decision-Making

Although praised for their rationality, humans often make poor decisions, even in simple situations. In the repeated binary choice experiment, an individual has to choose repeatedly between the same two alternatives, where a reward is assigned to one of them with fixed probability. The optimal strategy is to perseverate with choosing the alternative with the best expected return. Whereas many species perseverate, humans tend to match the frequencies of their choices to the frequencies of the alternatives, a sub-optimal strategy known as probability matching. Our goal was to find the primary cognitive constraints under which a set of simple evolutionary rules can lead to such contrasting behaviors. We simulated the evolution of artificial populations, wherein the fitness of each animat (artificial animal) depended on its ability to predict the next element of a sequence made up of a repeating binary string of varying size. When the string was short relative to the animats' neural capacity, they could learn it and correctly predict the next element of the sequence. When it was long, they could not learn it, turning to the next best option: to perseverate. Animats from the last generation then performed the task of predicting the next element of a non-periodical binary sequence. We found that, whereas animats with smaller neural capacity kept perseverating with the best alternative as before, animats with larger neural capacity, which had previously been able to learn the pattern of repeating strings, adopted probability matching, being outperformed by the perseverating animats. Our results demonstrate how the ability to make predictions in an environment endowed with regular patterns may lead to probability matching under less structured conditions. They point to probability matching as a likely by-product of adaptive cognitive strategies that were crucial in human evolution, but may lead to sub-optimal performances in other environments.

Characterization of depolymerized residues from extremely low acid hydrolysis (ELA) of sugarcane bagasse cellulose: Effects of degree of polymerization, crystallinity and crystallite size on thermal decomposition

Sugarcane bagasse cellulose was subjected to the extremely low acid (ELA) hydrolysis in 0.07% H2SO4 at 190, 210 and 225 degrees C for various times. The cellulose residues from this process were characterized by TGA, XRD, GPC, FIR and SEM. A kinetic study of thermal decomposition of the residues was also carried out, using the ASTM and Kissinger methods. The thermal studies revealed that residues of cellulose hydrolyzed at 190, 210 and 225 degrees C for 80,40 and 8 min have initial decomposition temperature and activation energy for the main decomposition step similar to those of Avicel PH-101. XRD studies confirmed this finding by showing that these cellulose residues are similar to Avicel in crystallinity index and crystallite size in relation to the 110 and 200 planes. FTIR spectra revealed no significant changes in the cellulose chemical structure and analysis of SEM micrographs demonstrated that the particle size of the cellulose residues hydrolyzed at 190 and 210 degrees C were similar to that of Avicel. (C) 2011 Elsevier B.V. All rights reserved.

Biological evaluation of the bone healing process after application of two potentially osteogenic proteins: an animal experimental model

doi: 10.1111/j.1741-2358.2011.00526.x Biological evaluation of the bone healing process after application of two potentially osteogenic proteins: an animal experimental model Objective: The aim of this work was to analyse qualitatively and quantitatively the newly formed bone after insertion of rhBMP-2 and protein extracted from Hevea brasiliensis (P-1), associated or not with a carrier in critical bone defects created in Wistar rat calvarial bone, using histological and histomorphometrical analyses. Materials and methods: Eighty-four male Wistar rats were used, divided into two groups, according to the period of time until the sacrifice (2 and 6 weeks). Each one of these groups was subdivided into six groups with seven animals each, according to the treatments: (1) 5 mu g of pure rhBMP-2, (2) 5 mu g of rhBMP-2/monoolein gel, (3) pure monoolein gel, (4) 5 mu g of pure P-1, (5) 5 mu g of P-1/monoolein gel and (6) critical bone defect controls. The animals were euthanised and the calvarial bone tissue removed for histological and histomorphometrical analyses. Result and conclusion: The results showed an improvement in the bone healing process using the rhBMP-2 protein, associated or not with a material carrier in relation to the other groups, and this process demonstrated to be time dependent.

DEVELOPMENT OF AN EXPERIMENTAL MODEL OF INFECTED BONE VOID IN THE ULNA OF RABBITS

Objective: Develop a model that allowed the study of bone regeneration in infection conditions. Method: A 15 mm defect was surgically created in the rabbit ulna and inoculated with 5x10(8) colony-forming units (CFU) of S. aureus. Surgical debridement was performed two weeks after and systemic gentamicin was administered for four weeks. Animals were followed up to 12 weeks to evaluate infection control and bone regeneration. Result: Bone regeneration was inferior to 25% of the defect in radiological and histological analysis. Conclusion: Infected bone defect of 15 mm in the rabbit ulna was unable to achieve full regeneration without further treatment. Level of Evidence V, Experimental Study.

Four-gene expression model predictive of lymph node metastases in oral squamous cell carcinoma

Background. Previous knowledge of cervical lymph node compromise may be crucial to choose the best treatment strategy in oral squamous cell carcinoma (OSCC). Here we propose a set four genes, whose mRNA expression in the primary tumor predicts nodal status in OSCC, excluding tongue. Material and methods. We identified differentially expressed genes in OSCC with and without compromised lymph nodes using Differential Display RT-PCR. Known genes were chosen to be validated by means of Northern blotting or real time RT-PCR (qRT-PCR). Thereafter we constructed a Nodal Index (NI) using discriminant analysis in a learning set of 35 patients, which was further validated in a second independent group of 20 patients. Results. Of the 63 differentially expressed known genes identified comparing three lymph node positive (pN+) and three negative (pN0) primary tumors, 23 were analyzed by Northern analysis or RT-PCR in 49 primary tumors. Six genes confirmed as differentially expressed were used to construct a NI, as the best set predictive of lymph nodal status, with the final result including four genes. The NI was able to correctly classify 32 of 35 patients comprising the learning group (88.6%; p = 0.009). Casein kinase 1alpha1 and scavenger receptor class B, member 2 were found to be up regulated in pN + group in contrast to small proline-rich protein 2B and Ras-GTPase activating protein SH3 domain-binding protein 2 which were upregulated in the pN0 group. We validated further our NI in an independent set of 20 primary tumors, 11 of them pN0 and nine pN+ with an accuracy of 80.0% (p = 0.012). Conclusions. The NI was an independent predictor of compromised lymph nodes, taking into the consideration tumor size and histological grade. The genes identified here that integrate our "Nodal Index" model are predictive of lymph node metastasis in OSCC.

Forced Desynchronization of Activity Rhythms in a Model of Chronic Jet Lag in Mice

We studied locomotor activity rhythms of C57/Bl6 mice under a chronic jet lag (CJL) protocol (ChrA(6/2)), which consisted of 6-hour phase advances of the light-dark schedule (LD) every 2 days. Through periodogram analysis, we found 2 components of the activity rhythm: a short-period component (21.01 +/- 0.04 h) that was entrained by the LD schedule and a long-period component (24.68 +/- 0.26 h). We developed a mathematical model comprising 2 coupled circadian oscillators that was tested experimentally with different CJL schedules. Our simulations suggested that under CJL, the system behaves as if it were under a zeitgeber with a period determined by (24 -[phase shift size/days between shifts]). Desynchronization within the system arises according to whether this effective zeitgeber is inside or outside the range of entrainment of the oscillators. In this sense, ChrA(6/2) is interpreted as a (24 - 6/2 = 21 h) zeitgeber, and simulations predicted the behavior of mice under other CJL schedules with an effective 21-hour zeitgeber. Animals studied under an asymmetric T = 21 h zeitgeber (carried out by a 3-hour shortening of every dark phase) showed 2 activity components as observed under ChrA(6/2): an entrained short-period (21.01 +/- 0.03 h) and a long-period component (23.93 +/- 0.31 h). Internal desynchronization was lost when mice were subjected to 9-hour advances every 3 days, a possibility also contemplated by the simulations. Simulations also predicted that desynchronization should be less prevalent under delaying than under advancing CJL. Indeed, most mice subjected to 6-hour delay shifts every 2 days (an effective 27-hour zeitgeber) displayed a single entrained activity component (26.92 +/- 0.11 h). Our results demonstrate that the disruption provoked by CJL schedules is not dependent on the phase-shift magnitude or the frequency of the shifts separately but on the combination of both, through its ratio and additionally on their absolute values. In this study, we present a novel model of forced desynchronization in mice under a specific CJL schedule; in addition, our model provides theoretical tools for the evaluation of circadian disruption under CJL conditions that are currently used in circadian research.

A model of road effect using line integrals and a test of the performance of two new road indices using the distribution of small mammals in an Atlantic Forest landscape

Effects of roads on wildlife and its habitat have been measured using metrics, such as the nearest road distance, road density, and effective mesh size. In this work we introduce two new indices: (1) Integral Road Effect (IRE), which measured the sum effects of points in a road at a fixed point in the forest; and (2) Average Value of the Infinitesimal Road Effect (AVIRE), which measured the average of the effects of roads at this point. IRE is formally defined as the line integral of a special function (the infinitesimal road effect) along the curves that model the roads, whereas AVIRE is the quotient of IRE by the length of the roads. Combining tools of ArcGIS software with a numerical algorithm, we calculated these and other road and habitat cover indices in a sample of points in a human-modified landscape in the Brazilian Atlantic Forest, where data on the abundance of two groups of small mammals (forest specialists and habitat generalists) were collected in the field. We then compared through the Akaike Information Criterion (AIC) a set of candidate regression models to explain the variation in small mammal abundance, including models with our two new road indices (AVIRE and IRE) or models with other road effect indices (nearest road distance, mesh size, and road density), and reference models (containing only habitat indices, or only the intercept without the effect of any variable). Compared to other road effect indices, AVIRE showed the best performance to explain abundance of forest specialist species, whereas the nearest road distance obtained the best performance to generalist species. AVIRE and habitat together were included in the best model for both small mammal groups, that is, higher abundance of specialist and generalist small mammals occurred where there is lower average road effect (less AVIRE) and more habitat. Moreover, AVIRE was not significantly correlated with habitat cover of specialists and generalists differing from the other road effect indices, except mesh size, which allows for separating the effect of roads from the effect of habitat on small mammal communities. We suggest that the proposed indices and GIS procedures could also be useful to describe other spatial ecological phenomena, such as edge effect in habitat fragments. (C) 2012 Elsevier B.V. All rights reserved.

Stripe-tetragonal phase transition in the two-dimensional Ising model with dipole interactions: Partition function zeros approach

We have performed multicanonical simulations to study the critical behavior of the two-dimensional Ising model with dipole interactions. This study concerns the thermodynamic phase transitions in the range of the interaction delta where the phase characterized by striped configurations of width h = 1 is observed. Controversial results obtained from local update algorithms have been reported for this region, including the claimed existence of a second-order phase transition line that becomes first order above a tricritical point located somewhere between delta = 0.85 and 1. Our analysis relies on the complex partition function zeros obtained with high statistics from multicanonical simulations. Finite size scaling relations for the leading partition function zeros yield critical exponents. that are clearly consistent with a single second-order phase transition line, thus excluding such a tricritical point in that region of the phase diagram. This conclusion is further supported by analysis of the specific heat and susceptibility of the orientational order parameter.

Why do we not find polyps in the lungs? Bronchial mucosa as a model in the treatment of polyposis

The link between lower and upper airways has been reported since the beginning of 1800s. They share the same pseudostratified ciliated columnar epithelium lining and the concept of one airway, one disease is quite well widespread. Nasal polyposis and asthma share basically the same inflammatory process: predominant infiltration of eosinophils, mucus cell hyperplasia, edema, thickened basal membrane, polarization for Th2 cell immune response, similar pro-inflammatory mediators are increased, for example cysteinyl leukotrienes. If the lower and upper airways share a lot of common epithelial structural features so why is the edema in the nasal mucosa able to increase so much the size of the mucosa to the point of developing polyps? The article tries to underline some differences between the nasal and the bronchial mucosa that could be implicated in this aberrant change from normal mucosa to polyps. This paper creates the concept that there are no polyps with the features of nasal polyposis disease in the lower airway and through it is developed the hypothesis of the nasal polyps origin could partially lie on the difference between the upper and lower airway histology. (C) 2012 Elsevier Ltd. All rights reserved.

Finite-size corrections of the entanglement entropy of critical quantum chains

Using the density matrix renormalization group, we calculated the finite-size corrections of the entanglement alpha-Renyi entropy of a single interval for several critical quantum chains. We considered models with U(1) symmetry such as the spin-1/2 XXZ and spin-1 Fateev-Zamolodchikov models, as well as models with discrete symmetries such as the Ising, the Blume-Capel, and the three-state Potts models. These corrections contain physically relevant information. Their amplitudes, which depend on the value of a, are related to the dimensions of operators in the conformal field theory governing the long-distance correlations of the critical quantum chains. The obtained results together with earlier exact and numerical ones allow us to formulate some general conjectures about the operator responsible for the leading finite-size correction of the alpha-Renyi entropies. We conjecture that the exponent of the leading finite-size correction of the alpha-Renyi entropies is p(alpha) = 2X(epsilon)/alpha for alpha > 1 and p(1) = nu, where X-epsilon denotes the dimensions of the energy operator of the model and nu = 2 for all the models.