The full Bayesian significance test for mixture models: results in gene expression clustering

Gene clustering is a useful exploratory technique to group together genes with similar expression levels under distinct cell cycle phases or distinct conditions. It helps the biologist to identify potentially meaningful relationships between genes. In this study, we propose a clustering method based on multivariate normal mixture models, where the number of clusters is predicted via sequential hypothesis tests: at each step, the method considers a mixture model of m components (m = 2 in the first step) and tests if in fact it should be m - 1. If the hypothesis is rejected, m is increased and a new test is carried out. The method continues (increasing m) until the hypothesis is accepted. The theoretical core of the method is the full Bayesian significance test, an intuitive Bayesian approach, which needs no model complexity penalization nor positive probabilities for sharp hypotheses. Numerical experiments were based on a cDNA microarray dataset consisting of expression levels of 205 genes belonging to four functional categories, for 10 distinct strains of Saccharomyces cerevisiae. To analyze the method's sensitivity to data dimension, we performed principal components analysis on the original dataset and predicted the number of classes using 2 to 10 principal components. Compared to Mclust (model-based clustering), our method shows more consistent results.

Ventilation behavior in trained and untrained men during incremental test: evidence of one metabolic transition point

This study aimed to describe and compare the ventilation behavior during an incremental test utilizing three mathematical models and to compare the feature of ventilation curve fitted by the best mathematical model between aerobically trained (TR) and untrained ( UT) men. Thirty five subjects underwent a treadmill test with 1 km.h(-1) increases every minute until exhaustion. Ventilation averages of 20 seconds were plotted against time and fitted by: bi-segmental regression model (2SRM); three-segmental regression model (3SRM); and growth exponential model (GEM). Residual sum of squares (RSS) and mean square error (MSE) were calculated for each model. The correlations between peak VO2 (VO2PEAK), peak speed (Speed(PEAK)), ventilatory threshold identified by the best model (VT2SRM) and the first derivative calculated for workloads below (moderate intensity) and above (heavy intensity) VT2SRM were calculated. The RSS and MSE for GEM were significantly higher (p < 0.01) than for 2SRM and 3SRM in pooled data and in UT, but no significant difference was observed among the mathematical models in TR. In the pooled data, the first derivative of moderate intensities showed significant negative correlations with VT2SRM (r = -0.58; p < 0.01) and Speed(PEAK) (r = -0.46; p < 0.05) while the first derivative of heavy intensities showed significant negative correlation with VT2SRM (r = -0.43; p < 0.05). In UT group the first derivative of moderate intensities showed significant negative correlations with VT2SRM (r = -0.65; p < 0.05) and Speed(PEAK) (r = -0.61; p < 0.05), while the first derivative of heavy intensities showed significant negative correlation with VT2SRM (r= -0.73; p < 0.01), Speed(PEAK) (r = -0.73; p < 0.01) and VO2PEAK (r = -0.61; p < 0.05) in TR group. The ventilation behavior during incremental treadmill test tends to show only one threshold. UT subjects showed a slower ventilation increase during moderate intensities while TR subjects showed a slower ventilation increase during heavy intensities.

Analysis of the tip roundness effects on the micro- and macroindentation response of elastic-plastic materials

In this work, the effects of indenter tip roundness oil the load-depth indentation curves were analyzed using finite element modeling. The tip roundness level was Studied based on the ratio between tip radius and maximum penetration depth (R/h(max)), which varied from 0.02 to 1. The proportional Curvature constant (C), the exponent of depth during loading (alpha), the initial unloading slope (S), the correction factor (beta), the level of piling-up or sinking-in (h(c)/h(max)), and the ratio h(max)/h(f) are shown to be strongly influenced by the ratio R/h(max). The hardness (H) was found to be independent of R/h(max) in the range studied. The Oliver and Pharr method was successful in following the variation of h(c)/h(max) with the ratio R/h(max) through the variation of S with the ratio R/h(max). However, this work confirmed the differences between the hardness values calculated using the Oliver-Pharr method and those obtained directly from finite element calculations; differences which derive from the error in area calculation that Occurs when given combinations of indented material properties are present. The ratio of plastic work to total work (W(p)/W(t)) was found to be independent of the ratio R/h(max), which demonstrates that the methods for the Calculation of mechanical properties based on the *indentation energy are potentially not Susceptible to errors caused by tip roundness.

Dynamical models for computer viruses propagation

Nowadays, digital computer systems and networks are the main engineering tools, being used in planning, design, operation, and control of all sizes of building, transportation, machinery, business, and life maintaining devices. Consequently, computer viruses became one of the most important sources of uncertainty, contributing to decrease the reliability of vital activities. A lot of antivirus programs have been developed, but they are limited to detecting and removing infections, based on previous knowledge of the virus code. In spite of having good adaptation capability, these programs work just as vaccines against diseases and are not able to prevent new infections based on the network state. Here, a trial on modeling computer viruses propagation dynamics relates it to other notable events occurring in the network permitting to establish preventive policies in the network management. Data from three different viruses are collected in the Internet and two different identification techniques, autoregressive and Fourier analyses, are applied showing that it is possible to forecast the dynamics of a new virus propagation by using the data collected from other viruses that formerly infected the network. Copyright (c) 2008 J. R. C. Piqueira and F. B. Cesar. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Assessing influence in survival data with a cure fraction and covariates

Diagnostic methods have been an important tool in regression analysis to detect anomalies, such as departures from error assumptions and the presence of outliers and influential observations with the fitted models. Assuming censored data, we considered a classical analysis and Bayesian analysis assuming no informative priors for the parameters of the model with a cure fraction. A Bayesian approach was considered by using Markov Chain Monte Carlo Methods with Metropolis-Hasting algorithms steps to obtain the posterior summaries of interest. Some influence methods, such as the local influence, total local influence of an individual, local influence on predictions and generalized leverage were derived, analyzed and discussed in survival data with a cure fraction and covariates. The relevance of the approach was illustrated with a real data set, where it is shown that, by removing the most influential observations, the decision about which model best fits the data is changed.

NEAR INFRARED SPECTROSCOPY FOR ESTIMATING SUGARCANE BAGASSE CONTENT IN MEDIUM DENSITY FIBERBOARD

Medium density fiberboard (MDF) is an engineered wood product formed by breaking down selected lignin-cellulosic material residuals into fibers, combining it with wax and a resin binder, and then forming panels by applying high temperature and pressure. Because the raw material in the industrial process is ever-changing, the panel industry requires methods for monitoring the composition of their products. The aim of this study was to estimate the ratio of sugarcane (SC) bagasse to Eucalyptus wood in MDF panels using near infrared (NIR) spectroscopy. Principal component analysis (PCA) and partial least square (PLS) regressions were performed. MDF panels having different bagasse contents were easily distinguished from each other by the PCA of their NIR spectra with clearly different patterns of response. The PLS-R models for SC content of these MDF samples presented a strong coefficient of determination (0.96) between the NIR-predicted and Lab-determined values and a low standard error of prediction (similar to 1.5%) in the cross-validations. A key role of resins (adhesives), cellulose, and lignin for such PLS-R calibrations was shown. PLS-DA model correctly classified ninety-four percent of MDF samples by cross-validations and ninety-eight percent of the panels by independent test set. These NIR-based models can be useful to quickly estimate sugarcane bagasse vs. Eucalyptus wood content ratio in unknown MDF samples and to verify the quality of these engineered wood products in an online process.

Evaluation of Antiproliferative, Anti-Type 2 Diabetes, and Antihypertension Potentials of Ellagitannins from Strawberries (Fragaria x ananassa Duch.) Using In Vitro Models

Strawberries represent the main source of ellagic acid derivatives in the Brazilian diet, corresponding to more than 50% of all phenolic compounds found in the fruit. There is a particular interest in the determination of the ellagic acid content in fruits because of possible chemopreventive benefits. In the present study, the potential health benefits of purified ellagitannins from strawberries were evaluated in relation to the antiproliferative activity and in vitro inhibition of alpha-amylase, alpha-glucosidase, and angiotensin I-converting enzyme (ACE) relevant for potential management of hyperglycemia and hypertension. Therefore, a comparison among ellagic acid, purified ellagitannins, and a strawberry extract was done to evaluate the possible synergistic effects of phenolics. In relation to the antiproliferative activity, it was observed that ellagic acid had the highest percentage inhibition of cell proliferation. The strawberry extract had lower efficacy in inhibiting the cell proliferation, indicating that in the case of this fruit there is no synergism. Purified ellagitannins had high alpha-amylase and ACE inhibitory activities. However, these compounds had low alpha-glucosidase inhibitory activity. These results suggested that the ellagitannins and ellagic acid have good potential for the management of hyperglycemia and hypertension linked to type 2 diabetes. However, further studies with animal and human models are needed to advance the in vitro assay-based biochemical rationale from this study.

Evaluation of Antihyperglycemia and Antihypertension Potential of Native Peruvian Fruits Using In Vitro Models

Local food diversity and traditional crops are essential for cost-effective management of the global epidemic of type 2 diabetes and associated complications of hypertension. Water and 12% ethanol extracts of native Peruvian fruits such as Lucuma (Pouteria lucuma), Pacae (Inga feuille), Papayita arequipena (Carica pubescens), Capuli (Prunus capuli), Aguaymanto (Physalis peruviana), and Algarrobo (Prosopis pallida) were evaluated for total phenolics, antioxidant activity based on 2, 2-diphenyl-1-picrylhydrazyl radical scavenging assay, and functionality such as in vitro inhibition of alpha-amylase, alpha-glucosidase, and angiotensin I-converting enzyme (ACE) relevant for potential management of hyperglycemia and hypertension linked to type 2 diabetes. The total phenolic content ranged from 3.2 (Aguaymanto) to 11.4 (Lucuma fruit) mg/g of sample dry weight. A significant positive correlation was found between total phenolic content and antioxidant activity for the ethanolic extracts. No phenolic compound was detected in Lucuma (fruit and powder) and Pacae. Aqueous extracts from Lucuma and Algarrobo had the highest alpha-glucosidase inhibitory activities. Papayita arequipena and Algarrobo had significant ACE inhibitory activities reflecting antihypertensive potential. These in vitro results point to the excellent potential of Peruvian fruits for food-based strategies for complementing effective antidiabetes and antihypertension solutions based on further animal and clinical studies.

Antidiabetes and antihypertension potential of commonly consumed carbohydrate sweeteners using in vitro models

Commonly consumed carbohydrate sweeteners derived from sugar cane, palm, and corn (syrups) were investigated to determine their potential to inhibit key enzymes relevant to Type 2 diabetes and hypertension based on the total phenolic content and antioxidant activity using in vitro models. Among sugar cane derivatives, brown sugars showed higher antidiabetes potential than white sugars; nevertheless, no angiotensin I-converting enzyme (ACE) inhibition was detected in both sugar classes. Brown sugar from Peru and Mauritius (dark muscovado) had the highest total phenolic content and 1,1-diphenyl-2-picrylhydrazyl radical scavenging activity, which correlated with a moderate inhibition of yeast alpha-glucosidase without showing a significant effect on porcine pancreatic alpha-amylase activity. In addition, chlorogenic acid quantified by high-performance liquid chromatography was detected in these sugars (128 +/- 6 and 144 +/- 2 mu g/g of sample weight, respectively). Date sugar exhibited high alpha-glucosidase, alpha-amylase, and ACE inhibitory activities that correlated with high total phenolic content and antioxidant activity. Neither phenolic compounds or antioxidant activity was detected in corn syrups, indicating that nonphenolic factors may be involved in their significant ability to inhibit alpha-glucosidase, alpha-amylase, and ACE. This study provides a strong biochemical rationale for further in vivo studies and useful information to make better dietary sweetener choices for Type 2 diabetes and hypertension management.

Data collapse, scaling functions, and analytical solutions of generalized growth models

We consider a nontrivial one-species population dynamics model with finite and infinite carrying capacities. Time-dependent intrinsic and extrinsic growth rates are considered in these models. Through the model per capita growth rate we obtain a heuristic general procedure to generate scaling functions to collapse data into a simple linear behavior even if an extrinsic growth rate is included. With this data collapse, all the models studied become independent from the parameters and initial condition. Analytical solutions are found when time-dependent coefficients are considered. These solutions allow us to perceive nontrivial transitions between species extinction and survival and to calculate the transition's critical exponents. Considering an extrinsic growth rate as a cancer treatment, we show that the relevant quantity depends not only on the intensity of the treatment, but also on when the cancerous cell growth is maximum.

An empirical evaluation of imputation accuracy for association statistics reveals increased type-I error rates in genome-wide associations

Background: Genome wide association studies (GWAS) are becoming the approach of choice to identify genetic determinants of complex phenotypes and common diseases. The astonishing amount of generated data and the use of distinct genotyping platforms with variable genomic coverage are still analytical challenges. Imputation algorithms combine directly genotyped markers information with haplotypic structure for the population of interest for the inference of a badly genotyped or missing marker and are considered a near zero cost approach to allow the comparison and combination of data generated in different studies. Several reports stated that imputed markers have an overall acceptable accuracy but no published report has performed a pair wise comparison of imputed and empiric association statistics of a complete set of GWAS markers. Results: In this report we identified a total of 73 imputed markers that yielded a nominally statistically significant association at P < 10(-5) for type 2 Diabetes Mellitus and compared them with results obtained based on empirical allelic frequencies. Interestingly, despite their overall high correlation, association statistics based on imputed frequencies were discordant in 35 of the 73 (47%) associated markers, considerably inflating the type I error rate of imputed markers. We comprehensively tested several quality thresholds, the haplotypic structure underlying imputed markers and the use of flanking markers as predictors of inaccurate association statistics derived from imputed markers. Conclusions: Our results suggest that association statistics from imputed markers showing specific MAF (Minor Allele Frequencies) range, located in weak linkage disequilibrium blocks or strongly deviating from local patterns of association are prone to have inflated false positive association signals. The present study highlights the potential of imputation procedures and proposes simple procedures for selecting the best imputed markers for follow-up genotyping studies.

Differential expression of microRNAs in mouse pain models

Background: MicroRNAs (miRNAs) are short non-coding RNAs that inhibit translation of target genes by binding to their mRNAs. The expression of numerous brain-specific miRNAs with a high degree of temporal and spatial specificity suggests that miRNAs play an important role in gene regulation in health and disease. Here we investigate the time course gene expression profile of miR-1, -16, and -206 in mouse dorsal root ganglion (DRG), and spinal cord dorsal horn under inflammatory and neuropathic pain conditions as well as following acute noxious stimulation. Results: Quantitative real-time polymerase chain reaction analyses showed that the mature form of miR-1, -16 and -206, is expressed in DRG and the dorsal horn of the spinal cord. Moreover, CFA-induced inflammation significantly reduced miRs-1 and -16 expression in DRG whereas miR-206 was downregulated in a time dependent manner. Conversely, in the spinal dorsal horn all three miRNAs monitored were upregulated. After sciatic nerve partial ligation, miR-1 and -206 were downregulated in DRG with no change in the spinal dorsal horn. On the other hand, axotomy increases the relative expression of miR-1, -16, and 206 in a time-dependent fashion while in the dorsal horn there was a significant downregulation of miR-1. Acute noxious stimulation with capsaicin also increased the expression of miR-1 and -16 in DRG cells but, on the other hand, in the spinal dorsal horn only a high dose of capsaicin was able to downregulate miR-206 expression. Conclusions: Our results indicate that miRNAs may participate in the regulatory mechanisms of genes associated with the pathophysiology of chronic pain as well as the nociceptive processing following acute noxious stimulation. We found substantial evidence that miRNAs are differentially regulated in DRG and the dorsal horn of the spinal cord under different pain states. Therefore, miRNA expression in the nociceptive system shows not only temporal and spatial specificity but is also stimulus-dependent.

Class II Correction with the Cantilever Bite Jumper

Objective: To identify the skeletal, dentoalveolar, and soft tissue changes that occur during Class II correction with the Cantilever Bite Jumper (CBJ). Materials and Methods: This prospective cephalometric study was conducted on 26 subjects with Class II division 1 malocclusion treated with the CBJ appliance. A comparison was made with 26 untreated subjects with Class II malocclusion. Lateral head films from before and after CBJ therapy were analyzed through conventional cephalometric and Johnston analyses. Results: Class II correction was accomplished by means of 2.9 mm apical base change, 1.5 mm distal movement of the maxillary molars, and 1.1 mm mesial movement of the mandibular molars. The CBJ exhibited good control of the vertical dimension. The main side effect of the CBJ is that the vertical force vectors of the telescope act as lever arms and can produce mesial tipping of the mandibular molars. Conclusions: The Cantilever Bite Jumper corrects Class II malocclusions with similar percentages of skeletal and dentoalveolar effects. (Angle Orthod. 2009:79;)

Additive genetic relationship of longevity with fertility and production traits in Nellore cattle based on bivariate models

Survival or longevity is an economically important trait in beef cattle. The main inconvenience for its inclusion in selection criteria is delayed recording of phenotypic data and the high computational demand for including survival in proportional hazard models. Thus, identification of a longevity-correlated trait that could be recorded early in life would be very useful for selection purposes. We estimated the genetic relationship of survival with productive and reproductive traits in Nellore cattle, including weaning weight (WW), post-weaning growth (PWG), muscularity (MUSC), scrotal circumference at 18 months (SC18), and heifer pregnancy (HP). Survival was measured in discrete time intervals and modeled through a sequential threshold model. Five independent bivariate Bayesian analyses were performed, accounting for cow survival and the five productive and reproductive traits. Posterior mean estimates for heritability (standard deviation in parentheses) were 0.55 (0.01) for WW, 0.25 (0.01) for PWG, 0.23 (0.01) for MUSC, and 0.48 (0.01) for SC18. The posterior mean estimates (95% confidence interval in parentheses) for the genetic correlation with survival were 0.16 (0.13-0.19), 0.30 (0.25-0.34), 0.31 (0.25-0.36), 0.07 (0.02-0.12), and 0.82 (0.78-0.86) for WW, PWG, MUSC, SC18, and HP, respectively. Based on the high genetic correlation and heritability (0.54) posterior mean estimates for HP, the expected progeny difference for HP can be used to select bulls for longevity, as well as for post-weaning gain and muscle score.

Comparison of different models to estimate genetic parameters for carcass traits in a commercial broiler line

The effect of genetic and non-genetic factors for carcass, breast meat and leg weights, and yields of a commercial broiler line were investigated using the restricted maximum likelihood method, considering four different animal models, including or excluding maternal genetic effect with covariance between direct and maternal genetic effects, and maternal permanent environmental effect. The likelihood ratio test was used to determine the most adequate model for each trait. For carcass, breast, and leg weight, and for carcass and breast yield, maternal genetic and permanent environmental effects as well as the covariance between direct and maternal genetic effects were significant. The estimates of direct and maternal heritability were 0.17 and 0.04 for carcass weight, 0.26 and 0.06 for breast weight, 0.22 and 0.02 for leg weight, 0.32 and 0.02 for carcass yield, and 0.52 and 0.04 for breast yield, respectively. For leg yield, maternal permanent environmental effect was important, in addition to direct genetic effects. For that trait, direct heritability and maternal permanent environmental variance as a proportion of the phenotypic variance were 0.43 and 0.02, respectively. The results indicate that ignoring maternal effects in the models, even though they were of small magnitude (0.02 to 0.06), tended to overestimate direct genetic variance and heritability for all traits.