Parametric correlation functions to model the structure of permanent environmental (co)variances in milk yield random regression models

The objective of the present study was to estimate milk yield genetic parameters applying random regression models and parametric correlation functions combined with a variance function to model animal permanent environmental effects. A total of 152,145 test-day milk yields from 7,317 first lactations of Holstein cows belonging to herds located in the southeastern region of Brazil were analyzed. Test-day milk yields were divided into 44 weekly classes of days in milk. Contemporary groups were defined by herd-test-day comprising a total of 2,539 classes. The model included direct additive genetic, permanent environmental, and residual random effects. The following fixed effects were considered: contemporary group, age of cow at calving (linear and quadratic regressions), and the population average lactation curve modeled by fourth-order orthogonal Legendre polynomial. Additive genetic effects were modeled by random regression on orthogonal Legendre polynomials of days in milk, whereas permanent environmental effects were estimated using a stationary or nonstationary parametric correlation function combined with a variance function of different orders. The structure of residual variances was modeled using a step function containing 6 variance classes. The genetic parameter estimates obtained with the model using a stationary correlation function associated with a variance function to model permanent environmental effects were similar to those obtained with models employing orthogonal Legendre polynomials for the same effect. A model using a sixth-order polynomial for additive effects and a stationary parametric correlation function associated with a seventh-order variance function to model permanent environmental effects would be sufficient for data fitting.

Antenatal steroid and tracheal occlusion restore vascular endothelial growth factor receptors in congenital diaphragmatic hernia rat model

OBJECTIVE: Investigate the effects of antenatal steroids and tracheal occlusion on pulmonary expression of vascular endothelial growth factor receptors in rats with nitrofen-induced congenital diaphragmatic hernia. STUDY DESIGN: Fetuses were exposed to nitrofen at embryonic day 9.5. Subgroups received dexamethasone or were operated on for tracheal occlusion, or received combined treatment. Morphologic variables were recorded. To analyze vascular endothelial growth factor receptor 1 and vascular endothelial growth factor receptor 2 expression, we performed Western blotting and immunohistochemistry. Morphologic variables were analyzed by analysis of variance and immunohistochemistry by Kruskal-Wallis test. RESULTS: Congenital diaphragmatic hernia decreased body weight, total lung weight, and lung-to-body weight ratio. Tracheal occlusion increased total lung weight and lung-to-body weight ratio (P < .05). Fetuses with congenital diaphragmatic hernia had reduced vascular endothelial growth factor receptor 1 and vascular endothelial growth factor receptor 2 expression, whereas steroids and tracheal occlusion increased their expression. Combined treatment increased expression of receptors, but had no additive effect. CONCLUSION: Vascular endothelial growth factor signaling disruption may be associated with pulmonary hypertension in congenital diaphragmatic hernia. Tracheal occlusion and steroids provide a pathway for restoring expression of vascular endothelial growth factor receptors.

The IWA Anaerobic Digestion Model No 1 (ADM1)

The IWA Anaerobic Digestion Modelling Task Group was established in 1997 at the 8th World Congress on,Anaerobic Digestion (Sendai, Japan) with the goal of developing a generalised anaerobic digestion model. The structured model includes multiple steps describing biochemical as well as physicochemical processes. The biochemical steps include disintegration from homogeneous particulates to carbohydrates, proteins and lipids; extracellular hydrolysis of these particulate substrates to sugars, amino acids, and long chain fatty acids (LCFA), respectively; acidogenesis from sugars and amino acids to volatile fatty acids (VFAs) and hydrogen; acetogenesis of LCFA and VFAs to acetate; and separate methanogenesis steps from acetate and hydrogen/CO2. The physico-chemical equations describe ion association and dissociation, and gas-liquid transfer. Implemented as a differential and algebraic equation (DAE) set, there are 26 dynamic state concentration variables, and 8 implicit algebraic variables per reactor vessel or element. Implemented as differential equations (DE) only, there are 32 dynamic concentration state variables.

Analysis of the interaction of polycyclic aromatic compounds in a model organism: integration of genotoxic and histopathological effects

Due to their toxicity, especially their carcinogenic potential, polycyclic aromatic hydrocarbons (PAHs) became priority pollutants in biomonitoring programmes and environmental policy, such as the European Water Framework Directive. The model substances tested in this study, namely benzo[b]fluoranthene (B[b]F), considered potentially carcinogenic to humans and an effector carcinogenic PAH to wildlife, and phenanthrene (Phe), deemed a non-carcinogenic PAH, are common PAHs in coastal waters, owning distinct properties reflected in different, albeit overlapping, mechanisms of toxicity. Still, as for similar PAHs, their interaction effects remain largely unknown. In order to study the genotoxic effects of caused by the interaction of carcinogenic and non-carcinogenic PAHs, and their relation to histopathological alterations, juvenile sea basses, Dicentrarchus labrax, a highly ecologically- and economically-relevant marine fish, were injected with different doses (5 and 10 μg.g-1 fish ww) of the two PAHs, isolated or in mixture, and incubated for 48 h. Individuals injected with B[b]F and the PAH mixture exhibited higher clastogenic/aneugenic effects and DNA strand breakage in blood cells, determined through the erythrocytic nuclear abnormalities (ENA) and Comet assays, respectively. Also, hepatic histopathological alterations were found in all animals, especially those injected with B[b]F and the PAH mixture, relating especially to inflammation. Still, Phe also exhibited genotoxic effects in sea bass, especially in higher doses, revealing a very significant acute effect that was accordant with the Microtox test performed undergone in parallel. Overall, sea bass was sensitive to B[b]F (a higher molecular weight PAH), likely due to efficient bioactivation of the pollutant (yielding genotoxic metabolites and reactive oxygen species), when compared to Phe, the latter revealing a more significant acute effect. The results indicate no significant additive effect between the substances, under the current experimental conditions. The present study highlights the importance of understanding PAH interactions in aquatic organisms, since they are usually present in the aquatic environment in complex mixtures.

A Multisectorial model of prices: the SAM approach

Social Accounting Matrices (SAM) are normally used to analyse the income generation process. They are also useful, however, for analysing the cost transmission and price formation mechanisms. For price contributions, Roland-Holst and Sancho (1995) used the SAM structure to analyse the price and cost linkages through a representation of the interdependence between activities, households and factors. This paper is a further analysis of the cost transmission mechanisms, in which I add the capital account to the endogenous components of the Roland-Holst and Sancho approach. By doing this I reflect the responses of prices to the exogenous shocks in savings and investment. I also present an additive decomposition of the global price effects into categories of interdependence that isolates the impact on price levels of shocks in the capital account. I use a 1994 Social Accounting Matrix to make an empirical application of the Catalan economy. Keywords: social accounting matrix, cost linkages, price transmission, capital account. JEL Classification: C63, C69, D59.

The additive genetic variance after bottlenecks is affected by the number of loci involved in epistatic interactions.

We investigated the role of the number of loci coding for a neutral trait on the release of additive variance for this trait after population bottlenecks. Different bottleneck sizes and durations were tested for various matrices of genotypic values, with initial conditions covering the allele frequency space. We used three different types of matrices. First, we extended Cheverud and Routman's model by defining matrices of "pure" epistasis for three and four independent loci; second, we used genotypic values drawn randomly from uniform, normal, and exponential distributions; and third we used two models of simple metabolic pathways leading to physiological epistasis. For all these matrices of genotypic values except the dominant metabolic pathway, we find that, as the number of loci increases from two to three and four, an increase in the release of additive variance is occurring. The amount of additive variance released for a given set of genotypic values is a function of the inbreeding coefficient, independently of the size and duration of the bottleneck. The level of inbreeding necessary to achieve maximum release in additive variance increases with the number of loci. We find that additive-by-additive epistasis is the type of epistasis most easily converted into additive variance. For a wide range of models, our results show that epistasis, rather than dominance, plays a significant role in the increase of additive variance following bottlenecks.

Additive functions in boolean models of gene regulatory network modules.

Gene-on-gene regulations are key components of every living organism. Dynamical abstract models of genetic regulatory networks help explain the genome's evolvability and robustness. These properties can be attributed to the structural topology of the graph formed by genes, as vertices, and regulatory interactions, as edges. Moreover, the actual gene interaction of each gene is believed to play a key role in the stability of the structure. With advances in biology, some effort was deployed to develop update functions in Boolean models that include recent knowledge. We combine real-life gene interaction networks with novel update functions in a Boolean model. We use two sub-networks of biological organisms, the yeast cell-cycle and the mouse embryonic stem cell, as topological support for our system. On these structures, we substitute the original random update functions by a novel threshold-based dynamic function in which the promoting and repressing effect of each interaction is considered. We use a third real-life regulatory network, along with its inferred Boolean update functions to validate the proposed update function. Results of this validation hint to increased biological plausibility of the threshold-based function. To investigate the dynamical behavior of this new model, we visualized the phase transition between order and chaos into the critical regime using Derrida plots. We complement the qualitative nature of Derrida plots with an alternative measure, the criticality distance, that also allows to discriminate between regimes in a quantitative way. Simulation on both real-life genetic regulatory networks show that there exists a set of parameters that allows the systems to operate in the critical region. This new model includes experimentally derived biological information and recent discoveries, which makes it potentially useful to guide experimental research. The update function confers additional realism to the model, while reducing the complexity and solution space, thus making it easier to investigate.

Use of Two-Part Regression Calibration Model to Correct for Measurement Error in Episodically Consumed Foods in a Single-Replicate Study Design: EPIC Case Study.

In epidemiologic studies, measurement error in dietary variables often attenuates association between dietary intake and disease occurrence. To adjust for the attenuation caused by error in dietary intake, regression calibration is commonly used. To apply regression calibration, unbiased reference measurements are required. Short-term reference measurements for foods that are not consumed daily contain excess zeroes that pose challenges in the calibration model. We adapted two-part regression calibration model, initially developed for multiple replicates of reference measurements per individual to a single-replicate setting. We showed how to handle excess zero reference measurements by two-step modeling approach, how to explore heteroscedasticity in the consumed amount with variance-mean graph, how to explore nonlinearity with the generalized additive modeling (GAM) and the empirical logit approaches, and how to select covariates in the calibration model. The performance of two-part calibration model was compared with the one-part counterpart. We used vegetable intake and mortality data from European Prospective Investigation on Cancer and Nutrition (EPIC) study. In the EPIC, reference measurements were taken with 24-hour recalls. For each of the three vegetable subgroups assessed separately, correcting for error with an appropriately specified two-part calibration model resulted in about three fold increase in the strength of association with all-cause mortality, as measured by the log hazard ratio. Further found is that the standard way of including covariates in the calibration model can lead to over fitting the two-part calibration model. Moreover, the extent of adjusting for error is influenced by the number and forms of covariates in the calibration model. For episodically consumed foods, we advise researchers to pay special attention to response distribution, nonlinearity, and covariate inclusion in specifying the calibration model.

Assessing an intuitive condition for stability under a range of traffic conditions via a generalised Lu-Kumar network

We argue the importance both of developing simple sufficientconditions for the stability of general multiclass queueing networks and also of assessing such conditions under a range of assumptions on the weight of the traffic flowing between service stations. To achieve the former, we review a peak-rate stability condition and extend its range of application and for the latter, we introduce a generalisation of the Lu-Kumar network on which the stability condition may be tested for a range of traffic configurations. The peak-rate condition is close to exact when the between-station traffic is light, but degrades as this traffic increases.

Mometasone and desloratadine additive effect on eosinophil survival and cytokine secretion from epithelial cells

Although antihistamines and topical corticosteroids are used in combination to treat allergic rhinitis, their additive effect has not been yet demonstrated. The aim was investigate the antiinflammatory additive effect of mometasone and desloratadine on cytokine and sICAM-1 secretion by epithelial cells, and on eosinophil survival stimulated by human epithelial cells secretions from nasal mucosa and polyps. Methods Epithelial cells obtained from nasal mucosa or polyps were stimulated with 10% fetal bovine serum in presence of mometasone (10-11M-10-5M) with/without desloratadine (10-5M). Cytokine and sICAM-1 concentrations in supernatants were measured by ELISA. Peripheral blood eosinophils were incubated during 4 days with epithelial cell secretions with (10-11M-10-5M) and/or desloratadine (10-5M) and survival assessed by Trypan blue. Results are expressed as percentage (mean ± SEM) compared to control. Results Fetal bovine serum stimulated IL-6, IL-8, GM-CSF and sICAM-1 secretion. In mucosa and polyp epithelial cells, mometasone inhibited this induced secretion while desloratadine inhibited IL-6 and IL-8. The combination of 10-5M desloratadine and 10-9M mometasone reduced IL-6 secretion (48 ± 11%, p < 0.05) greater extent than mometasone alone (68 ± 10%) compared to control (100%). Epithelial cell secretions induced eosinophil survival from day 1 to 4, this effect being inhibited by mometasone. At day 4, the combination of mometasone (10-11M) and desloratadine (10-5M) provoked an increased inhibition of eosinophil survival induced by cell secretions (27 ± 5%, p < 0.01) than mometasone (44 ± 7%) or desloratadine (46 ± 7%) alone. Conclusions These results suggest that the combination of desloratadine and mometasone furoate have a greater antinflammatory effect in an in vitro model of eosinophil inflammation than those drugs administered alone.

Finite mixture analysis of beauty-contest data using generalised beta distributions

This paper introduces a mixture model based on the beta distribution, without preestablishedmeans and variances, to analyze a large set of Beauty-Contest data obtainedfrom diverse groups of experiments (Bosch-Domenech et al. 2002). This model gives a bettert of the experimental data, and more precision to the hypothesis that a large proportionof individuals follow a common pattern of reasoning, described as iterated best reply (degenerate),than mixture models based on the normal distribution. The analysis shows thatthe means of the distributions across the groups of experiments are pretty stable, while theproportions of choices at dierent levels of reasoning vary across groups.

Genotype by environment interaction for seed yield per plant in rapeseed using AMMI model

The objective of this study was to assess genotype by environment interaction for seed yield per plant in rapeseed cultivars grown in Northern Serbia by the AMMI (additive main effects and multiplicative interaction) model. The study comprised 19 rapeseed genotypes, analyzed in seven years through field trials arranged in a randomized complete block design, with three replicates. Seed yield per plant of the tested cultivars varied from 1.82 to 19.47 g throughout the seven seasons, with an average of 7.41 g. In the variance analysis, 72.49% of the total yield variation was explained by environment, 7.71% by differences between genotypes, and 19.09% by genotype by environment interaction. On the biplot, cultivars with high yield genetic potential had positive correlation with the seasons with optimal growing conditions, while the cultivars with lower yield potential were correlated to the years with unfavorable conditions. Seed yield per plant is highly influenced by environmental factors, which indicates the adaptability of specific genotypes to specific seasons.

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.

Derivation of an observer model adapted to irregular signals based on convolution channels.

Anthropomorphic model observers are mathe- matical algorithms which are applied to images with the ultimate goal of predicting human signal detection and classification accuracy across varieties of backgrounds, image acquisitions and display conditions. A limitation of current channelized model observers is their inability to handle irregularly-shaped signals, which are common in clinical images, without a high number of directional channels. Here, we derive a new linear model observer based on convolution channels which we refer to as the "Filtered Channel observer" (FCO), as an extension of the channelized Hotelling observer (CHO) and the nonprewhitening with an eye filter (NPWE) observer. In analogy to the CHO, this linear model observer can take the form of a single template with an external noise term. To compare with human observers, we tested signals with irregular and asymmetrical shapes spanning the size of lesions down to those of microcalfications in 4-AFC breast tomosynthesis detection tasks, with three different contrasts for each case. Whereas humans uniformly outperformed conventional CHOs, the FCO observer outperformed humans for every signal with only one exception. Additive internal noise in the models allowed us to degrade model performance and match human performance. We could not match all the human performances with a model with a single internal noise component for all signal shape, size and contrast conditions. This suggests that either the internal noise might vary across signals or that the model cannot entirely capture the human detection strategy. However, the FCO model offers an efficient way to apprehend human observer performance for a non-symmetric signal.

Image quality in CT: From physical measurements to model observers.

Evaluation of image quality (IQ) in Computed Tomography (CT) is important to ensure that diagnostic questions are correctly answered, whilst keeping radiation dose to the patient as low as is reasonably possible. The assessment of individual aspects of IQ is already a key component of routine quality control of medical x-ray devices. These values together with standard dose indicators can be used to give rise to 'figures of merit' (FOM) to characterise the dose efficiency of the CT scanners operating in certain modes. The demand for clinically relevant IQ characterisation has naturally increased with the development of CT technology (detectors efficiency, image reconstruction and processing), resulting in the adaptation and evolution of assessment methods. The purpose of this review is to present the spectrum of various methods that have been used to characterise image quality in CT: from objective measurements of physical parameters to clinically task-based approaches (i.e. model observer (MO) approach) including pure human observer approach. When combined together with a dose indicator, a generalised dose efficiency index can be explored in a framework of system and patient dose optimisation. We will focus on the IQ methodologies that are required for dealing with standard reconstruction, but also for iterative reconstruction algorithms. With this concept the previously used FOM will be presented with a proposal to update them in order to make them relevant and up to date with technological progress. The MO that objectively assesses IQ for clinically relevant tasks represents the most promising method in terms of radiologist sensitivity performance and therefore of most relevance in the clinical environment.