LATENT RANK CHANGE DETECTION FOR ANALYSIS OF SPLICE-JUNCTION MICROARRAYS WITH NONLINEAR EFFECTS

Alternative splicing of gene transcripts greatly expands the functional capacity of the genome, and certain splice isoforms may indicate specific disease states such as cancer. Splice junction microarrays interrogate thousands of splice junctions, but data analysis is difficult and error prone because of the increased complexity compared to differential gene expression analysis. We present Rank Change Detection (RCD) as a method to identify differential splicing events based upon a straightforward probabilistic model comparing the over-or underrepresentation of two or more competing isoforms. RCD has advantages over commonly used methods because it is robust to false positive errors due to nonlinear trends in microarray measurements. Further, RCD does not depend on prior knowledge of splice isoforms, yet it takes advantage of the inherent structure of mutually exclusive junctions, and it is conceptually generalizable to other types of splicing arrays or RNA-Seq. RCD specifically identifies the biologically important cases when a splice junction becomes more or less prevalent compared to other mutually exclusive junctions. The example data is from different cell lines of glioblastoma tumors assayed with Agilent microarrays.

Diagnostic errors and repetitive sequential classifications in on-line process control by attributes

The procedure of on-line process control by attributes, known as Taguchi`s on-line process control, consists of inspecting the mth item (a single item) at every m produced items and deciding, at each inspection, whether the fraction of conforming items was reduced or not. If the inspected item is nonconforming, the production is stopped for adjustment. As the inspection system can be subject to diagnosis errors, one develops a probabilistic model that classifies repeatedly the examined item until a conforming or b non-conforming classification is observed. The first event that occurs (a conforming classifications or b non-conforming classifications) determines the final classification of the examined item. Proprieties of an ergodic Markov chain were used to get the expression of average cost of the system of control, which can be optimized by three parameters: the sampling interval of the inspections (m); the number of repeated conforming classifications (a); and the number of repeated non-conforming classifications (b). The optimum design is compared with two alternative approaches: the first one consists of a simple preventive policy. The production system is adjusted at every n produced items (no inspection is performed). The second classifies the examined item repeatedly r (fixed) times and considers it conforming if most classification results are conforming. Results indicate that the current proposal performs better than the procedure that fixes the number of repeated classifications and classifies the examined item as conforming if most classifications were conforming. On the other hand, the preventive policy can be averagely the most economical alternative rather than those ones that require inspection depending on the degree of errors and costs. A numerical example illustrates the proposed procedure. (C) 2009 Elsevier B. V. All rights reserved.

Economic design of an X chart for short-run production

The aim of this paper is to present an economical design of an X chart for a short-run production. The process mean starts equal to mu(0) (in-control, State I) and in a random time it shifts to mu(1) > mu(0) (out-of-control, State II). The monitoring procedure consists of inspecting a single item at every m produced ones. If the measurement of the quality characteristic does not meet the control limits, the process is stopped, adjusted, and additional (r - 1) items are inspected retrospectively. The probabilistic model was developed considering only shifts in the process mean. A direct search technique is applied to find the optimum parameters which minimizes the expected cost function. Numerical examples illustrate the proposed procedure. (C) 2009 Elsevier B.V. All rights reserved.

A practical approach to assess depression risk and to guide risk reduction strategies in later life

Background: Many factors have been associated with the onset and maintenance of depressive symptoms in later life, although this knowledge is yet to be translated into significant health gains for the population. This study gathered information about common modifiable and non-modifiable risk factors for depression with the aim of developing a practical probabilistic model of depression that can be used to guide risk reduction strategies. \Methods: A cross-sectional study was undertaken of 20,677 community-dwelling Australians aged 60 years or over in contact with their general practitioner during the preceding 12 months. Prevalent depression (minor or major) according to the Patient Health Questionnaire (PHQ-9) assessment was the main outcome of interest. Other measured exposures included self-reported age, gender, education, loss of mother or father before age 15 years, physical or sexual abuse before age 15 years, marital status, financial stress, social support, smoking and alcohol use, physical activity, obesity, diabetes, hypertension, and prevalent cardiovascular diseases, chronic respiratory diseases and cancer. Results: The mean age of participants was 71.7 +/- 7.6 years and 57.9% were women. Depression was present in 1665 (8.0%) of our subjects. Multivariate logistic regression showed depression was independently associated with age older than 75 years, childhood adverse experiences, adverse lifestyle practices (smoking, risk alcohol use, physical inactivity), intermediate health hazards (obesity, diabetes and hypertension), comorbid medical conditions (clinical history of coronary heart disease, stroke, asthma, chronic obstructive pulmonary disease, emphysema or cancers), and social or financial strain. We stratified the exposures to build a matrix that showed that the probability of depression increased progressively with the accumulation of risk factors, from less than 3% for those with no adverse factors to more than 80% for people reporting the maximum number of risk factors. Conclusions: Our probabilistic matrix can be used to estimate depression risk and to guide the introduction of risk reduction strategies. Future studies should now aim to clarify whether interventions designed to mitigate the impact of risk factors can change the prevalence and incidence of depression in later life.

Coupled reliability and boundary element model for probabilistic fatigue life assessment in mixed mode crack propagation

Fatigue and crack propagation are phenomena affected by high uncertainties, where deterministic methods fail to predict accurately the structural life. The present work aims at coupling reliability analysis with boundary element method. The latter has been recognized as an accurate and efficient numerical technique to deal with mixed mode propagation, which is very interesting for reliability analysis. The coupled procedure allows us to consider uncertainties during the crack growth process. In addition, it computes the probability of fatigue failure for complex structural geometry and loading. Two coupling procedures are considered: direct coupling of reliability and mechanical solvers and indirect coupling by the response surface method. Numerical applications show the performance of the proposed models in lifetime assessment under uncertainties, where the direct method has shown faster convergence than response surface method. (C) 2010 Elsevier Ltd. All rights reserved.

Monte Carlo investigation of the critical behavior of Stavskaya's probabilistic cellular automaton

Stavskaya's model is a one-dimensional probabilistic cellular automaton (PCA) introduced in the end of the 1960s as an example of a model displaying a nonequilibrium phase transition. Although its absorbing state phase transition is well understood nowadays, the model never received a full numerical treatment to investigate its critical behavior. In this Brief Report we characterize the critical behavior of Stavskaya's PCA by means of Monte Carlo simulations and finite-size scaling analysis. The critical exponents of the model are calculated and indicate that its phase transition belongs to the directed percolation universality class of critical behavior, as would be expected on the basis of the directed percolation conjecture. We also explicitly establish the relationship of the model with the Domany-Kinzel PCA on its directed site percolation line, a connection that seems to have gone unnoticed in the literature so far.

An extinction-survival-type phase transition in the probabilistic cellular automaton p182-q200

We investigate the critical behaviour of a probabilistic mixture of cellular automata (CA) rules 182 and 200 (in Wolfram`s enumeration scheme) by mean-field analysis and Monte Carlo simulations. We found that as we switch off one CA and switch on the other by the variation of the single parameter of the model, the probabilistic CA (PCA) goes through an extinction-survival-type phase transition, and the numerical data indicate that it belongs to the directed percolation universality class of critical behaviour. The PCA displays a characteristic stationary density profile and a slow, diffusive dynamics close to the pure CA 200 point that we discuss briefly. Remarks on an interesting related stochastic lattice gas are addressed in the conclusions.

Validation of an experimental polyurethane model for biomechanical studies on implant supported prosthesis - tension tests

OBJECTIVES: The complexity and heterogeneity of human bone, as well as ethical issues, frequently hinder the development of clinical trials. The purpose of this in vitro study was to determine the modulus of elasticity of a polyurethane isotropic experimental model via tension tests, comparing the results to those reported in the literature for mandibular bone, in order to validate the use of such a model in lieu of mandibular bone in biomechanical studies. MATERIAL AND METHODS: Forty-five polyurethane test specimens were divided into 3 groups of 15 specimens each, according to the ratio (A/B) of polyurethane reagents (PU-1: 1/0.5, PU-2: 1/1, PU-3: 1/1.5). RESULTS: Tension tests were performed in each experimental group and the modulus of elasticity values found were 192.98 MPa (SD=57.20) for PU-1, 347.90 MPa (SD=109.54) for PU-2 and 304.64 MPa (SD=25.48) for PU-3. CONCLUSION: The concentration of choice for building the experimental model was 1/1.

Validation of an experimental polyurethane model for biomechanical studies on implant-supported prosthesis: compression tests

OBJECTIVES: The complexity and heterogeneity of human bone, as well as ethical issues, most always hinder the performance of clinical trials. Thus, in vitro studies become an important source of information for the understanding of biomechanical events on implant-supported prostheses, although study results cannot be considered reliable unless validation studies are conducted. The purpose of this work was to validate an artificial experimental model based on its modulus of elasticity, to simulate the performance of human bone in vivo in biomechanical studies of implant-supported prostheses. MATERIAL AND METHODS: In this study, fast-curing polyurethane (F16 polyurethane, Axson) was used to build 40 specimens that were divided into five groups. The following reagent ratios (part A/part B) were used: Group A (0.5/1.0), Group B (0.8/1.0), Group C (1.0/1.0), Group D (1.2/1.0), and Group E (1.5/1.0). A universal testing machine (Kratos model K - 2000 MP) was used to measure modulus of elasticity values by compression. RESULTS: Mean modulus of elasticity values were: Group A - 389.72 MPa, Group B - 529.19 MPa, Group C - 571.11 MPa, Group D - 470.35 MPa, Group E - 437.36 MPa. CONCLUSION: The best mechanical characteristics and modulus of elasticity value comparable to that of human trabecular bone were obtained when A/B ratio was 1:1.

Lesion of the subthalamic nucleus reverses motor deficits but not death of nigrostriatal dopaminergic neurons in a rat 6-hydroxydopamine-lesion model of Parkinson's disease

The objective of the present study was to determine whether lesion of the subthalamic nucleus (STN) promoted by N-methyl-D-aspartate (NMDA) would rescue nigrostriatal dopaminergic neurons after unilateral 6-hydroxydopamine (6-OHDA) injection into the medial forebrain bundle (MFB). Initially, 16 mg 6-OHDA (6-OHDA group) or vehicle (artificial cerebrospinal fluid - aCSF; Sham group) was infused into the right MFB of adult male Wistar rats. Fifteen days after surgery, the 6-OHDA and SHAM groups were randomly subdivided and received ipsilateral injection of either 60 mM NMDA or aCSF in the right STN. Additionally, a control group was not submitted to stereotaxic surgery. Five groups of rats were studied: 6-OHDA/NMDA, 6-OHDA/Sham, Sham/NMDA, Sham/Sham, and Control. Fourteen days after injection of 6-OHDA, rats were submitted to the rotational test induced by apomorphine (0.1 mg/kg, ip) and to the open-field test. The same tests were performed again 14 days after NMDA-induced lesion of the STN. The STN lesion reduced the contralateral turns induced by apomorphine and blocked the progression of motor impairment in the open-field test in 6-OHDA-treated rats. However, lesion of the STN did not prevent the reduction of striatal concentrations of dopamine and metabolites or the number of nigrostriatal dopaminergic neurons after 6-OHDA lesion. Therefore, STN lesion is able to reverse motor deficits after severe 6-OHDA-induced lesion of the nigrostriatal pathway, but does not protect or rescue dopaminergic neurons in the substantia nigra pars compacta.

Viability of using enamel and dentin from bovine origin as a substitute for human counterparts in an intraoral erosion model

This study ascertained whether under dental erosion models that closely mimics the real-life situation enamel and root dentin from bovine origin would be reliable substitutes for human counterparts. Through a 2x2 crossover design, in a first trial, 14 volunteers wore a palatal device containing slabs of bovine and human enamel. Half of the participants ingested (4x daily, for 10 days) orange juice first, crossing over to mineral water, while the remainder received the reverse sequence. In a second trial, volunteers wore devices with slabs of bovine and human root dentin. Except for the duration of each intraoral phase, which lasted 2 rather 10 days, the experiment with root dentin run exactly as for enamel. Dental substrates were analyzed for surface microhardness. Two-way ANOVAs (α=0.05) indicated no difference between the microhardness values recorded for human and bovine enamel (p=0.1350), but bovine root dentin had lower microhardness compared to its human counterpart (p=0.0432). While bovine enamel can reliably substitute its human counterpart in in situ dental erosion models, bovine root dentin does not seem to be a viable alternative to the corresponding human tissue.

An experimental model for the study of craniofacial deformities

PURPOSE: To develop an experimental surgical model in rats for the study of craniofacial abnormalities. METHODS: Full thickness calvarial defects with 10x10-mm and 5x8-mm dimensions were created in 40 male NIS Wistar rats, body weight ranging from 320 to 420 g. The animals were equally divided into two groups. The periosteum was removed and dura mater was left intact. Animals were killed at 8 and 16 weeks postoperatively and cranial tissue samples were taken from the defects for histological analysis. RESULTS: Cranial defects remained open even after 16 weeks postoperatively. CONCLUSION: The experimental model with 5x8-mm defects in the parietal region with the removal of the periosteum and maintenance of the integrity of the dura mater are critical and might be used for the study of cranial bone defects in craniofacial abnormalities.

Validation of the Deardorff model for estimating energy balance components for a sugarcane crop

The quantification of the available energy in the environment is important because it determines photosynthesis, evapotranspiration and, therefore, the final yield of crops. Instruments for measuring the energy balance are costly and indirect estimation alternatives are desirable. This study assessed the Deardorff's model performance during a cycle of a sugarcane crop in Piracicaba, State of São Paulo, Brazil, in comparison to the aerodynamic method. This mechanistic model simulates the energy fluxes (sensible, latent heat and net radiation) at three levels (atmosphere, canopy and soil) using only air temperature, relative humidity and wind speed measured at a reference level above the canopy, crop leaf area index, and some pre-calibrated parameters (canopy albedo, soil emissivity, atmospheric transmissivity and hydrological characteristics of the soil). The analysis was made for different time scales, insolation conditions and seasons (spring, summer and autumn). Analyzing all data of 15 minute intervals, the model presented good performance for net radiation simulation in different insolations and seasons. The latent heat flux in the atmosphere and the sensible heat flux in the atmosphere did not present differences in comparison to data from the aerodynamic method during the autumn. The sensible heat flux in the soil was poorly simulated by the model due to the poor performance of the soil water balance method. The Deardorff's model improved in general the flux simulations in comparison to the aerodynamic method when more insolation was available in the environment.

Study of the equatorial Atlantic Ocean mixing layer using a one-dimensional turbulence model

The General Ocean Turbulence Model (GOTM) is applied to the diagnostic turbulence field of the mixing layer (ML) over the equatorial region of the Atlantic Ocean. Two situations were investigated: rainy and dry seasons, defined, respectively, by the presence of the intertropical convergence zone and by its northward displacement. Simulations were carried out using data from a PIRATA buoy located on the equator at 23º W to compute surface turbulent fluxes and from the NASA/GEWEX Surface Radiation Budget Project to close the surface radiation balance. A data assimilation scheme was used as a surrogate for the physical effects not present in the one-dimensional model. In the rainy season, results show that the ML is shallower due to the weaker surface stress and stronger stable stratification; the maximum ML depth reached during this season is around 15 m, with an averaged diurnal variation of 7 m depth. In the dry season, the stronger surface stress and the enhanced surface heat balance components enable higher mechanical production of turbulent kinetic energy and, at night, the buoyancy acts also enhancing turbulence in the first meters of depth, characterizing a deeper ML, reaching around 60 m and presenting an average diurnal variation of 30 m.

Immune cells and oxidative stress in the endotoxin tolerance mouse model

Sepsis is a systemic inflammatory response that can lead to tissue damage and death. In order to increase our understanding of sepsis, experimental models are needed that produce relevant immune and inflammatory responses during a septic event. We describe a lipopolysaccharide tolerance mouse model to characterize the cellular and molecular alterations of immune cells during sepsis. The model presents a typical lipopolysaccharide tolerance pattern in which tolerance is related to decreased production and secretion of cytokines after a subsequent exposure to a lethal dose of lipopolysaccharide. The initial lipopolysaccharide exposure also altered the expression patterns of cytokines and was followed by an 8- and a 1.5-fold increase in the T helper 1 and 2 cell subpopulations. Behavioral data indicate a decrease in spontaneous activity and an increase in body temperature following exposure to lipopolysaccharide. In contrast, tolerant animals maintained production of reactive oxygen species and nitric oxide when terminally challenged by cecal ligation and puncture (CLP). Survival study after CLP showed protection in tolerant compared to naive animals. Spleen mass increased in tolerant animals followed by increases of B lymphocytes and subpopulation Th1 cells. An increase in the number of stem cells was found in spleen and bone marrow. We also showed that administration of spleen or bone marrow cells from tolerant to naive animals transfers the acquired resistance status. In conclusion, lipopolysaccharide tolerance is a natural reprogramming of the immune system that increases the number of immune cells, particularly T helper 1 cells, and does not reduce oxidative stress.