BadiRate: Estimating Family Turnover Rates by Likelihood-Based Methods

Motivation: The comparative analysis of gene gain and loss rates is critical for understanding the role of natural selection and adaptation in shaping gene family sizes. Studying complete genome data from closely related species allows accurate estimation of gene family turnover rates. Current methods and software tools, however, are not well designed for dealing with certain kinds of functional elements, such as microRNAs or transcription factor binding sites. Results: Here, we describe BadiRate, a new software tool to estimate family turnover rates, as well as the number of elements in internal phylogenetic nodes, by likelihood-based methods and parsimony. It implements two stochastic population models, which provide the appropriate statistical framework for testing hypothesis, such as lineage-specific gene family expansions or contractions. We have assessed the accuracy of BadiRate by computer simulations, and have also illustrated its functionality by analyzing a representative empirical dataset.

Molecular phylogenetics and temporal diversification in the genus Aeromonas based on the sequences of five housekeeping genes

Several approaches have been developed to estimate both the relative and absolute rates of speciation and extinction within clades based on molecular phylogenetic reconstructions of evolutionary relationships, according to an underlying model of diversification. However, the macroevolutionary models established for eukaryotes have scarcely been used with prokaryotes. We have investigated the rate and pattern of cladogenesis in the genus Aeromonas (γ-Proteobacteria, Proteobacteria, Bacteria) using the sequences of five housekeeping genes and an uncorrelated relaxed-clock approach. To our knowledge, until now this analysis has never been applied to all the species described in a bacterial genus and thus opens up the possibility of establishing models of speciation from sequence data commonly used in phylogenetic studies of prokaryotes. Our results suggest that the genus Aeromonas began to diverge between 248 and 266 million years ago, exhibiting a constant divergence rate through the Phanerozoic, which could be described as a pure birth process.

Effects of sample size on the performance of species distribution models

A wide range of modelling algorithms is used by ecologists, conservation practitioners, and others to predict species ranges from point locality data. Unfortunately, the amount of data available is limited for many taxa and regions, making it essential to quantify the sensitivity of these algorithms to sample size. This is the first study to address this need by rigorously evaluating a broad suite of algorithms with independent presence-absence data from multiple species and regions. We evaluated predictions from 12 algorithms for 46 species (from six different regions of the world) at three sample sizes (100, 30, and 10 records). We used data from natural history collections to run the models, and evaluated the quality of model predictions with area under the receiver operating characteristic curve (AUC). With decreasing sample size, model accuracy decreased and variability increased across species and between models. Novel modelling methods that incorporate both interactions between predictor variables and complex response shapes (i.e. GBM, MARS-INT, BRUTO) performed better than most methods at large sample sizes but not at the smallest sample sizes. Other algorithms were much less sensitive to sample size, including an algorithm based on maximum entropy (MAXENT) that had among the best predictive power across all sample sizes. Relative to other algorithms, a distance metric algorithm (DOMAIN) and a genetic algorithm (OM-GARP) had intermediate performance at the largest sample size and among the best performance at the lowest sample size. No algorithm predicted consistently well with small sample size (n < 30) and this should encourage highly conservative use of predictions based on small sample size and restrict their use to exploratory modelling.

Association of calcemia and serum vitamin D with 24H-urinary calcium excretion in a swiss population-based study

Background: Elevated urinary calcium excretion is associated with reduced bone mineral density. Population-based data on urinary calcium excretion are scarce. We explored the association of serum calcium and circulating levels of vitamin D (including 25(OH)D2 and 25(OH)D3) with urinary calcium excretion in men and women in a population-based study. Methods: We used data from the "Swiss Survey on Salt" conducted between 2010 and 2012 and including people aged 15 years and over. Twenty-four hour urine collection, blood analysis, clinical examination and anthropometric measures were collected in 11 centres from the 3 linguistic regions of Switzerland. Vitamin D was measured centrally using liquid chromatography - tandem mass spectrometry. Hypercalciuria was defined as urinary calcium excretion >0.1 mmol/kg/24h. Multivariable linear regression was used to explore factors associated with 24-hour urinary calcium excretion (mmol/24h) squared root transformed, taken as the dependant variable. Vitamin D was divided into monthspecific tertiles with the first tertile having the lowest value and the third tertile having the highest value. Results: The 669 men and 624 women had mean (SD) age of 49.2 (18.1) and 47 (17.9) years and a prevalence of hypercalciuria of 8.9% and 8.0%, respectively. In adjusted models, the association of urinary calcium excretion with protein-corrected serum calcium was (β coefficient } standard error, according to urinary calcium squared root transformed) 1.125 } 0.184 mmol/L per square-root (mmol/24h) (P<0.001) in women and 0.374 } 0.224 (P=0.096) in men. Men in the third month-specific vitamin D tertile had higher urinary calcium excretion than men in the first tertile (0.170 } 0.05 nmol/L per mmol/24h, P=0.001) and the corresponding association was 0.048 } 0.043, P= 0.272 in women. Conclusion: About one in eleven person has hypercalciuria in the Swiss population. The positive association of serum calcium with urinary calcium excretion was steeper in women than in men, independently of menopausal status. Circulating vitamin D was associated positively with urinary calcium excretion only in men. The reasons underlying the observed sex differences in the hormonal control of urinary calcium excretion need to be explored in further studies.

Estimating parameters for generalized mass action models using constraint propagation.

As modern molecular biology moves towards the analysis of biological systems as opposed to their individual components, the need for appropriate mathematical and computational techniques for understanding the dynamics and structure of such systems is becoming more pressing. For example, the modeling of biochemical systems using ordinary differential equations (ODEs) based on high-throughput, time-dense profiles is becoming more common-place, which is necessitating the development of improved techniques to estimate model parameters from such data. Due to the high dimensionality of this estimation problem, straight-forward optimization strategies rarely produce correct parameter values, and hence current methods tend to utilize genetic/evolutionary algorithms to perform non-linear parameter fitting. Here, we describe a completely deterministic approach, which is based on interval analysis. This allows us to examine entire sets of parameters, and thus to exhaust the global search within a finite number of steps. In particular, we show how our method may be applied to a generic class of ODEs used for modeling biochemical systems called Generalized Mass Action Models (GMAs). In addition, we show that for GMAs our method is amenable to the technique in interval arithmetic called constraint propagation, which allows great improvement of its efficiency. To illustrate the applicability of our method we apply it to some networks of biochemical reactions appearing in the literature, showing in particular that, in addition to estimating system parameters in the absence of noise, our method may also be used to recover the topology of these networks.

Associations between mood, anxiety or substance use disorders and inflammatory markers after adjustment for multiple covariates in a population-based study.

Inflammation is one possible mechanism underlying the associations between mental disorders and cardiovascular diseases (CVD). However, studies on mental disorders and inflammation have yielded inconsistent results and the majority did not adjust for potential confounding factors. We examined the associations of several pro-inflammatory cytokines (IL-1β, IL-6 and TNF-α) and high sensitive C-reactive protein (hsCRP) with lifetime and current mood, anxiety and substance use disorders (SUD), while adjusting for multiple covariates. The sample included 3719 subjects, randomly selected from the general population, who underwent thorough somatic and psychiatric evaluations. Psychiatric diagnoses were made with a semi-structured interview. Major depressive disorder was subtyped into "atypical", "melancholic", "combined atypical-melancholic" and "unspecified". Associations between inflammatory markers and psychiatric diagnoses were assessed using multiple linear and logistic regression models. Lifetime bipolar disorders and atypical depression were associated with increased levels of hsCRP, but not after multivariate adjustment. After multivariate adjustment, SUD remained associated with increased hsCRP levels in men (β = 0.13 (95% CI: 0.03,0.23)) but not in women. After multivariate adjustment, lifetime combined and unspecified depression were associated with decreased levels of IL-6 (β = -0.27 (-0.51,-0.02); β = -0.19 (-0.34,-0.05), respectively) and TNF-α (β = -0.16 (-0.30,-0.01); β = -0.10 (-0.19,-0.02), respectively), whereas current combined and unspecified depression were associated with decreased levels of hsCRP (β = -0.20 (-0.39,-0.02); β = -0.12 (-0.24,-0.01), respectively). Our data suggest that the significant associations between increased hsCRP levels and mood disorders are mainly attributable to the effects of comorbid disorders, medication as well as behavioral and physical CVRFs.

BadiRate: estimating family turnover rates by likelihood-based methods

Motivation: The comparative analysis of gene gain and loss rates is critical for understanding the role of natural selection and adaptation in shaping gene family sizes. Studying complete genome data from closely related species allows accurate estimation of gene family turnover rates. Current methods and software tools, however, are not well designed for dealing with certain kinds of functional elements, such as microRNAs or transcription factor binding sites. Results: Here, we describe BadiRate, a new software tool to estimate family turnover rates, as well as the number of elements in internal phylogenetic nodes, by likelihood-based methods and parsimony. It implements two stochastic population models, which provide the appropriate statistical framework for testing hypothesis, such as lineage-specific gene family expansions or contractions. We have assessed the accuracy of BadiRate by computer simulations, and have also illustrated its functionality by analyzing a representative empirical dataset.

Nutrient-based dietary patterns and the risk of head and neck cancer: a pooled analysis in the International Head and Neck Cancer Epidemiology consortium.

Background The association between dietary patterns and head and neck cancer has rarely been addressed. Patients and methods We used individual-level pooled data from five case-control studies (2452 cases and 5013 controls) participating in the International Head and Neck Cancer Epidemiology consortium. A posteriori dietary patterns were identified through a principal component factor analysis carried out on 24 nutrients derived from study-specific food-frequency questionnaires. Odds ratios (ORs) and corresponding 95% confidence intervals (CIs) were estimated using unconditional logistic regression models on quintiles of factor scores. Results We identified three major dietary patterns named 'animal products and cereals', 'antioxidant vitamins and fiber', and 'fats'. The 'antioxidant vitamins and fiber' pattern was inversely related to oral and pharyngeal cancer (OR = 0.57, 95% CI 0.43-0.76 for the highest versus the lowest score quintile). The 'animal products and cereals' pattern was positively associated with laryngeal cancer (OR = 1.54, 95% CI 1.12-2.11), whereas the 'fats' pattern was inversely associated with oral and pharyngeal cancer (OR = 0.78, 95% CI 0.63-0.97) and positively associated with laryngeal cancer (OR = 1.69, 95% CI 1.22-2.34). Conclusions These findings suggest that diets rich in animal products, cereals, and fats are positively related to laryngeal cancer, and those rich in fruit and vegetables inversely related to oral and pharyngeal cancer.

Comparison of an EMG-based and a stress-based method to predict shoulder muscle forces.

The estimation of muscle forces in musculoskeletal shoulder models is still controversial. Two different methods are widely used to solve the indeterminacy of the system: electromyography (EMG)-based methods and stress-based methods. The goal of this work was to evaluate the influence of these two methods on the prediction of muscle forces, glenohumeral load and joint stability after total shoulder arthroplasty. An EMG-based and a stress-based method were implemented into the same musculoskeletal shoulder model. The model replicated the glenohumeral joint after total shoulder arthroplasty. It contained the scapula, the humerus, the joint prosthesis, the rotator cuff muscles supraspinatus, subscapularis and infraspinatus and the middle, anterior and posterior deltoid muscles. A movement of abduction was simulated in the plane of the scapula. The EMG-based method replicated muscular activity of experimentally measured EMG. The stress-based method minimised a cost function based on muscle stresses. We compared muscle forces, joint reaction force, articular contact pressure and translation of the humeral head. The stress-based method predicted a lower force of the rotator cuff muscles. This was partly counter-balanced by a higher force of the middle part of the deltoid muscle. As a consequence, the stress-based method predicted a lower joint load (16% reduced) and a higher superior-inferior translation of the humeral head (increased by 1.2 mm). The EMG-based method has the advantage of replicating the observed cocontraction of stabilising muscles of the rotator cuff. This method is, however, limited to available EMG measurements. The stress-based method has thus an advantage of flexibility, but may overestimate glenohumeral subluxation.

Degree of multicollinearity and variables involved in linear dependence in additive-dominant models

The objective of this work was to assess the degree of multicollinearity and to identify the variables involved in linear dependence relations in additive-dominant models. Data of birth weight (n=141,567), yearling weight (n=58,124), and scrotal circumference (n=20,371) of Montana Tropical composite cattle were used. Diagnosis of multicollinearity was based on the variance inflation factor (VIF) and on the evaluation of the condition indexes and eigenvalues from the correlation matrix among explanatory variables. The first model studied (RM) included the fixed effect of dam age class at calving and the covariates associated to the direct and maternal additive and non-additive effects. The second model (R) included all the effects of the RM model except the maternal additive effects. Multicollinearity was detected in both models for all traits considered, with VIF values of 1.03 - 70.20 for RM and 1.03 - 60.70 for R. Collinearity increased with the increase of variables in the model and the decrease in the number of observations, and it was classified as weak, with condition index values between 10.00 and 26.77. In general, the variables associated with additive and non-additive effects were involved in multicollinearity, partially due to the natural connection between these covariables as fractions of the biological types in breed composition.

Rice yield estimation based on weather conditions and on technological level of production systems in Brazil

The objective of this work was to evaluate an estimation system for rice yield in Brazil, based on simple agrometeorological models and on the technological level of production systems. This estimation system incorporates the conceptual basis proposed by Doorenbos & Kassam for potential and attainable yields with empirical adjusts for maximum yield and crop sensitivity to water deficit, considering five categories of rice yield. Rice yield was estimated from 2000/2001 to 2007/2008, and compared to IBGE yield data. Regression analyses between model estimates and data from IBGE surveys resulted in significant coefficients of determination, with less dispersion in the South than in the North and Northeast regions of the country. Index of model efficiency (E1') ranged from 0.01 in the lower yield classes to 0.45 in higher ones, and mean absolute error ranged from 58 to 250 kg ha‑1, respectively.

Determinants of science-based cooperation: evidence in a sample of small and micro firms

We study the determining factors of cience-based cooperation in the case of small and micro firms. In this research, we propose an analytical framework based on the resource-based view of the firm and we identify a set of organisational characteristics, which we classify as internal, external and structural factors. Each factor can be linked to at least one reason, from the firm¿s point of view, to cooperate with universities and public research centres. Each reason can, in turn, be used as an indicator of a firm¿s organisational needs or organisational capacities. In order to validate the theoretical model, we estimate a logistic regression that models the propensity to participate in science-based cooperation activities within a sample of 285 small and micro firms located in Barcelona. The results show the key role played by the absorptive capacity of new and small companies.

Ambiguous nucleotide calls from population-based sequencing of HIV-1 are a marker for viral diversity and the age of infection.

Background. The time passed since the infection of a human immunodeficiency virus (HIV)-infected individual (the age of infection) is an important but often only poorly known quantity. We assessed whether the fraction of ambiguous nucleotides obtained from bulk sequencing as done for genotypic resistance testing can serve as a proxy of this parameter. Methods. We correlated the age of infection and the fraction of ambiguous nucleotides in partial pol sequences of HIV-1 sampled before initiation of antiretroviral therapy (ART). Three groups of Swiss HIV Cohort Study participants were analyzed, for whom the age of infection was estimated on the basis of Bayesian back calculation (n = 3,307), seroconversion (n = 366), or diagnoses of primary HIV infection (n = 130). In addition, we studied 124 patients for whom longitudinal genotypic resistance testing was performed while they were still ART-naive. Results. We found that the fraction of ambiguous nucleotides increased with the age of infection with a rate of .2% per year within the first 8 years but thereafter with a decreasing rate. We show that this pattern is consistent with population-genetic models for realistic parameters. Finally, we show that, in this highly representative population, a fraction of ambiguous nucleotides of >.5% provides strong evidence against a recent infection event < 1 year prior to sampling (negative predictive value, 98.7%). Conclusions. These findings show that the fraction of ambiguous nucleotides is a useful marker for the age of infection.

Environmental data mining and modeling based on machine learning algorithms and geostatistics

The paper presents some contemporary approaches to spatial environmental data analysis. The main topics are concentrated on the decision-oriented problems of environmental spatial data mining and modeling: valorization and representativity of data with the help of exploratory data analysis, spatial predictions, probabilistic and risk mapping, development and application of conditional stochastic simulation models. The innovative part of the paper presents integrated/hybrid model-machine learning (ML) residuals sequential simulations-MLRSS. The models are based on multilayer perceptron and support vector regression ML algorithms used for modeling long-range spatial trends and sequential simulations of the residuals. NIL algorithms deliver non-linear solution for the spatial non-stationary problems, which are difficult for geostatistical approach. Geostatistical tools (variography) are used to characterize performance of ML algorithms, by analyzing quality and quantity of the spatially structured information extracted from data with ML algorithms. Sequential simulations provide efficient assessment of uncertainty and spatial variability. Case study from the Chernobyl fallouts illustrates the performance of the proposed model. It is shown that probability mapping, provided by the combination of ML data driven and geostatistical model based approaches, can be efficiently used in decision-making process. (C) 2003 Elsevier Ltd. All rights reserved.

Development of Craniofacial Bone Defect Reconstruction Implant Based on Fibre-Reinforced Composite with Photopolymerisable Resin Systems. Experimental studies in vitro and in vivo.

Reconstruction of defects in the craniomaxillofacial (CMF) area has mainly been based on bone grafts or metallic fixing plates and screws. Particularly in the case of large calvarial and/or craniofacial defects caused by trauma, tumours or congenital malformations, there is a need for reliable reconstruction biomaterials, because bone grafts or metallic fixing systems do not completely fulfill the criteria for the best possible reconstruction methods in these complicated cases. In this series of studies, the usability of fibre-reinforced composite (FRC) was studied as a biostable, nonmetallic alternative material for reconstructing artificially created bone defects in frontal and calvarial areas of rabbits. The experimental part of this work describes the different stages of the product development process from the first in vitro tests with resin-impregnated fibrereinforced composites to the in vivo animal studies, in which this FRC was tested as an implant material for reconstructing different size bone defects in rabbit frontal and calvarial areas. In the first in vitro study, the FRC was polymerised in contact with bone or blood in the laboratory. The polymerised FRC samples were then incubated in water, which was analysed for residual monomer content by using high performance liquid chromatography (HPLC). It was found that this in vitro polymerisation in contact with bone and blood did not markedly increase the residual monomer leaching from the FRC. In the second in vitro study, different adhesive systems were tested in fixing the implant to bone surface. This was done to find an alternative implant fixing system to screws and pins. On the basis of this study, it was found that the surface of the calvarial bone needed both mechanical and chemical treatments before the resinimpregnated FRC could be properly fixed onto it. In three animal studies performed with rabbit frontal bone defects and critical size calvarial bone defect models, biological responses to the FRC implants were evaluated. On the basis of theseevaluations, it can be concluded that the FRC, based on E-glass (electrical glass) fibres forming a porous fibre veil enables the ingrowth of connective tissues to the inner structures of the material, as well as the bone formation and mineralization inside the fibre veil. Bone formation could be enhanced by using bioactive glass granules fixed to the FRC implants. FRC-implanted bone defects healed partly; no total healing of defects was achieved. Biological responses during the follow-up time, at a maximum of 12 weeks, to resin-impregnated composite implant seemed to depend on the polymerization time of the resin matrix of the FRC. Both of the studied resin systems used in the FRC were photopolymerised and the heat-induced postpolymerisation was used additionally.