Spatial statistical models that use flow and stream distance

We develop spatial statistical models for stream networks that can estimate relationships between a response variable and other covariates, make predictions at unsampled locations, and predict an average or total for a stream or a stream segment. There have been very few attempts to develop valid spatial covariance models that incorporate flow, stream distance, or both. The application of typical spatial autocovariance functions based on Euclidean distance, such as the spherical covariance model, are not valid when using stream distance. In this paper we develop a large class of valid models that incorporate flow and stream distance by using spatial moving averages. These methods integrate a moving average function, or kernel, against a white noise process. By running the moving average function upstream from a location, we develop models that use flow, and by construction they are valid models based on stream distance. We show that with proper weighting, many of the usual spatial models based on Euclidean distance have a counterpart for stream networks. Using sulfate concentrations from an example data set, the Maryland Biological Stream Survey (MBSS), we show that models using flow may be more appropriate than models that only use stream distance. For the MBSS data set, we use restricted maximum likelihood to fit a valid covariance matrix that uses flow and stream distance, and then we use this covariance matrix to estimate fixed effects and make kriging and block kriging predictions.

PHOTOELASTIC ANALISYS IN THE LOWER REGION OF VERTEBRAL BODY L4

Objective: To analyze the shear forces on the vertebral body L4 when submitted to a compression force by means of transmission photoelasticity. Methods: Twelve photoelastic models were divided into three groups, with four models per group, according to the positioning of the sagittal section vertebrae L4-L5 (sections A, B and C). The simulation was performed using a 15N compression force, and the fringe orders were evaluated in the vertebral body L4 by the Tardy compensation method. Results: Photoelastic analysis showed, in general, a homogeneous distribution in the vertebral bodies. The shear forces were higher in section C than B, and higher in B than A. Conclusion: The posterior area of L4, mainly in section C, showed higher shear concentrations, corresponding to a more susceptible area for bone fracture and spondylolisthesis. Economic and Decision Analyses Development of an Economic or Decision Model. Level I

A Systematic Mapping on the Role-Permission Relationship in Role Based Access Control Models

Access control is a key component of security in any computer system. In the last two decades, the research on Role Basead Access Control Models was intense. One of the most important components of a Role Based Model is the Role-Permission Relationship. In this paper, the technique of systematic mapping is used to identify, extract and analyze many approaches applied to establish the Role-Permission Relationship. The main goal of this mapping is pointing directions of significant research in the area of Role Based Access Control Models.

Molecular evolution of a malaria resistance gene (DARC) in primates

Genes involved in host-pathogen interactions are often strongly affected by positive natural selection. The Duffy antigen, coded by the Duffy antigen receptor for chemokines (DARC) gene, serves as a receptor for Plasmodium vivax in humans and for Plasmodium knowlesi in some nonhuman primates. In the majority of sub-Saharan Africans, a nucleic acid variant in GATA-1 of the gene promoter is responsible for the nonexpression of the Duffy antigen on red blood cells and consequently resistance to invasion by P. vivax. The Duffy antigen also acts as a receptor for chemokines and is expressed in red blood cells and many other tissues of the body. Because of this dual role, we sequenced a 3,000-bp region encompassing the entire DARC gene as well as part of its 5' and 3' flanking regions in a phylogenetic sample of primates and used statistical methods to evaluate the nature of selection pressures acting on the gene during its evolution. We analyzed both coding and regulatory regions of the DARC gene. The regulatory analysis showed accelerated rates of substitution at several sites near known motifs. Our tests of positive selection in the coding region using maximum likelihood by branch sites and maximum likelihood by codon sites did not yield statistically significant evidence for the action of positive selection. However, the maximum likelihood test in which the gene was subdivided into different structural regions showed that the known binding region for P. vivax/P. knowlesi is under very different selective pressures than the remainder of the gene. In fact, most of the gene appears to be under strong purifying selection, but this is not evident in the binding region. We suggest that the binding region is under the influence of two opposing selective pressures, positive selection possibly exerted by the parasite and purifying selection exerted by chemokines.

Bayesian Analysis for Errors in Variables with Changepoint Models

Changepoint regression models have originally been developed in connection with applications in quality control, where a change from the in-control to the out-of-control state has to be detected based on the avaliable random observations. Up to now various changepoint models have been suggested for differents applications like reliability, econometrics or medicine. In many practical situations the covariate cannot be measured precisely and an alternative model are the errors in variable regression models. In this paper we study the regression model with errors in variables with changepoint from a Bayesian approach. From the simulation study we found that the proposed procedure produces estimates suitable for the changepoint and all other model parameters.

Influence diagnostics in Gaussian spatial linear models

Spatial linear models have been applied in numerous fields such as agriculture, geoscience and environmental sciences, among many others. Spatial dependence structure modelling, using a geostatistical approach, is an indispensable tool to estimate the parameters that define this structure. However, this estimation may be greatly affected by the presence of atypical observations in the sampled data. The purpose of this paper is to use diagnostic techniques to assess the sensitivity of the maximum-likelihood estimators, covariance functions and linear predictor to small perturbations in the data and/or the spatial linear model assumptions. The methodology is illustrated with two real data sets. The results allowed us to conclude that the presence of atypical values in the sample data have a strong influence on thematic maps, changing the spatial dependence structure.

Evolution of form and function: morphophysiological relationships and locomotor performance in tropidurine lizards

Locomotor capacity is often considered an excellent measure of whole animal performance because it requires the integrated functioning of many morphological, physiological (and biochemical) traits. However, because studies tend to focus on either structural or functional suits of traits, we know little on whether and how morphological and physiological traits coevolve to produce adequate locomotor capacities. Hence, we investigate the evolutionary relationships between morphological and physiological parameters related to exercise physiology, using tropidurine lizards as a model. We employ a phylogenetic principal component analysis (PCA) to identify variable clusters (factors) related to morphology, energetic metabolism and muscle metabolism, and then analyze the relationships between these clusters and measures of locomotor performance, using two models (star and hierarchical phylogenies). Our data indicate that sprint performance is enhanced by simultaneous evolutionary tendencies affecting relative limb and tail size and physiological traits. Specifically, the high absolute sprint speeds exhibited by tropidurines from the sand dunes are explained by longer limbs, feet and tails and an increased proportion of glycolytic fibers in the leg muscle, contrasting with their lower capacity for overall oxidative metabolism [principal component (PC1)]. However, when sprint speeds are corrected for body size, performance correlates with a cluster (PC3) composed by moderate loads for activity metabolic rate and body size. The simultaneous measurement of morphological and physiological parameters is a powerful tool for exploring patterns of coadaptation and proposing morphophysiological associations that are not directly predictable from theory. This approach may trigger novel directions for investigating the evolution of form and function, particularly in the context of organismal performance.

Segmental dataset and whole body expression data do not support the hypothesis that non-random movement is an intrinsic property of Drosophila retrogenes

Background: Several studies in Drosophila have shown excessive movement of retrogenes from the X chromosome to autosomes, and that these genes are frequently expressed in the testis. This phenomenon has led to several hypotheses invoking natural selection as the process driving male-biased genes to the autosomes. Metta and Schlotterer (BMC Evol Biol 2010, 10:114) analyzed a set of retrogenes where the parental gene has been subsequently lost. They assumed that this class of retrogenes replaced the ancestral functions of the parental gene, and reported that these retrogenes, although mostly originating from movement out of the X chromosome, showed female-biased or unbiased expression. These observations led the authors to suggest that selective forces (such as meiotic sex chromosome inactivation and sexual antagonism) were not responsible for the observed pattern of retrogene movement out of the X chromosome. Results: We reanalyzed the dataset published by Metta and Schlotterer and found several issues that led us to a different conclusion. In particular, Metta and Schlotterer used a dataset combined with expression data in which significant sex-biased expression is not detectable. First, the authors used a segmental dataset where the genes selected for analysis were less testis-biased in expression than those that were excluded from the study. Second, sex-biased expression was defined by comparing male and female whole-body data and not the expression of these genes in gonadal tissues. This approach significantly reduces the probability of detecting sex-biased expressed genes, which explains why the vast majority of the genes analyzed (parental and retrogenes) were equally expressed in both males and females. Third, the female-biased expression observed by Metta and Schltterer is mostly found for parental genes located on the X chromosome, which is known to be enriched with genes with female-biased expression. Fourth, using additional gonad expression data, we found that autosomal genes analyzed by Metta and Schlotterer are less up regulated in ovaries and have higher chance to be expressed in meiotic cells of spermatogenesis when compared to X-linked genes. Conclusions: The criteria used to select retrogenes and the sex-biased expression data based on whole adult flies generated a segmental dataset of female-biased and unbiased expressed genes that was unable to detect the higher propensity of autosomal retrogenes to be expressed in males. Thus, there is no support for the authors' view that the movement of new retrogenes, which originated from X-linked parental genes, was not driven by selection. Therefore, selection-based genetic models remain the most parsimonious explanations for the observed chromosomal distribution of retrogenes.

Importance of xenarthrans in the eco-epidemiology of Paracoccidioides brasiliensis

Abstract Background Several pathogens that cause important zoonotic diseases have been frequently associated with armadillos and other xenarthrans. This mammal group typically has evolved on the South American continent and many of its extant species are seriously threatened with extinction. Natural infection of armadillos with Paracoccidioides brasiliensis in hyperendemic areas has provided a valuable opportunity for understanding the role of this mammal in the eco-epidemiology of Paracoccidioidomycosis (PCM), one of the most important systemic mycoses in Latin America. Findings This study aimed to detect P. brasiliensis in different xenarthran species (Dasypus novemcinctus, Cabassous spp., Euphractus sexcinctus, Tamandua tetradactyla and Myrmecophaga tridactyla), by molecular and mycological approaches, in samples obtained by one of the following strategies: i) from road-killed animals (n = 6); ii) from naturally dead animals (n = 8); iii) from animals that died in captivity (n = 9); and iv) from living animals captured from the wild (n = 2). Specific P. brasiliensis DNA was detected in several organs among 7/20 nine-banded armadillos (D. novemcinctus) and in 2/2 anteaters (M. tridactyla). The fungus was also cultured in tissue samples from one of two armadillos captured from the wild. Conclusion Members of the Xenarthra Order, especially armadillos, have some characteristics, including a weak cellular immune response and low body temperature, which make them suitable models for studying host-pathogen interaction. P. brasiliensis infection in wild animals, from PCM endemic areas, may be more common than initially postulated and reinforces the use of these animals as sentinels for the pathogen in the environment.

Body weight, metabolism and clock genes

Biological rhythms are present in the lives of almost all organisms ranging from plants to more evolved creatures. These oscillations allow the anticipation of many physiological and behavioral mechanisms thus enabling coordination of rhythms in a timely manner, adaption to environmental changes and more efficient organization of the cellular processes responsible for survival of both the individual and the species. Many components of energy homeostasis exhibit circadian rhythms, which are regulated by central (suprachiasmatic nucleus) and peripheral (located in other tissues) circadian clocks. Adipocyte plays an important role in the regulation of energy homeostasis, the signaling of satiety and cellular differentiation and proliferation. Also, the adipocyte circadian clock is probably involved in the control of many of these functions. Thus, circadian clocks are implicated in the control of energy balance, feeding behavior and consequently in the regulation of body weight. In this regard, alterations in clock genes and rhythms can interfere with the complex mechanism of metabolic and hormonal anticipation, contributing to multifactorial diseases such as obesity and diabetes. The aim of this review was to define circadian clocks by describing their functioning and role in the whole body and in adipocyte metabolism, as well as their influence on body weight control and the development of obesity.

Validation of models for predicting milk urea nitrogen concentrations, estimating dry matter intake by the NRC (2001)

The objective of this study was to validate three different models for predicting milk urea nitrogen using field conditions, attempting to evaluate the nutritional adequacy diets for dairy cows and prediction of nitrogen excreted to the environment. Observations (4,749) from 855 cows were used. Milk yield, body weight (BW), days in milk and parity were recorded on the milk sampling days. Milk was sampled monthly, for analysis of milk urea nitrogen (MUN), fat, protein, lactose and total solids concentration and somatic cells count. Individual dry matter intake was estimated using the NRC (2001). The three models studied were derived from a first one to predict urinary nitrogen (UN). Model 1 was MUN = UN/12.54, model 2 was MUN = UN/17.6 and model 3 was MUN = UN/(0.0259 × BW), adjusted by body weight effect. To evaluate models, they were tested for accuracy, precision and robustness. Despite being more accurate (mean bias = 0.94 mg/dL), model 2 was less precise (residual error = 4.50 mg/dL) than model 3 (mean bias = 1.41 and residual error = 4.11 mg/dL), while model 1 was the least accurate (mean bias = 6.94 mg/dL) and the least precise (residual error = 5.40 mg/dL). They were not robust, because they were influenced by almost all the variables studied. The three models for predicting milk urea nitrogen were different with respect to accuracy, precision and robustness.

Relationship between body image disturbance and incidence of depression: the SUN prospective cohort

[EN] Background: Body image disturbance is an increasing problem in Western societies and is associated with a number of mental health outcomes including anorexia, bulimia, body dysmorphia, and depression. The aim of this study was to assess the association between body image disturbance and the incidence of depression. Methods: This study included 10,286 participants from a dynamic prospective cohort of Spanish university graduates, who were followed-up for a median period of 4.2 years (Seguimiento Universidad de Navarra – the SUN study). The key characteristic of the study is the permanently open recruitment that started in 1999. The baseline questionnaire included information about body mass index (BMI) and the nine figure schemes that were used to assess body size perception. These variables were grouped according to recommended classifications and the difference between BMI and body size perception was considered as a proxy of body image disturbance. A subject was classified as an incident case of depression if he/she was initially free of depression and reported a physician-made diagnosis of depression and/or the use of antidepressant medication in at least one of the follow-up questionnaires. The association between body image disturbance and the incidence of depression was estimated by calculating the multivariable adjusted Odds Ratio (OR) and its 95% Confidence Interval (95% CI), using logistic regression models. Results: The cumulative incidence of depression during follow-up in the cohort was 4.8%. Men who underestimated their body size had a high percentage of overweight and obesity (50.1% and 12.6%, respectively), whereas women who overestimated their body size had a high percentage of underweight (87.6%). The underestimation exhibited a negative association with the incidence of depression among women (OR: 0.72, 95% CI: 0.54 – 0.95), but this effect disappeared after adjusting for possible confounding variables. The proportion of participants who correctly perceived their body size was high (53.3%) and gross misperception was seldom found, with most cases selecting only one silhouette below (42.7%) or above (2.6%) their actual BMI. Conclusion: We found no association between body image disturbance and subsequent depression in a cohort of university graduates in Spain.

The use of image analysis systems for the assesment of zooplankton physiological rates and carbon fluxes

Trabajo realizado por: Garijo, J. C., Hernández León, S.

Development of musculoskeletal models for the design and the pre-clinical validation of hip resurfacing prosthesis

Background. The surgical treatment of dysfunctional hips is a severe condition for the patient and a costly therapy for the public health. Hip resurfacing techniques seem to hold the promise of various advantages over traditional THR, with particular attention to young and active patients. Although the lesson provided in the past by many branches of engineering is that success in designing competitive products can be achieved only by predicting the possible scenario of failure, to date the understanding of the implant quality is poorly pre-clinically addressed. Thus revision is the only delayed and reliable end point for assessment. The aim of the present work was to model the musculoskeletal system so as to develop a protocol for predicting failure of hip resurfacing prosthesis. Methods. Preliminary studies validated the technique for the generation of subject specific finite element (FE) models of long bones from Computed Thomography data. The proposed protocol consisted in the numerical analysis of the prosthesis biomechanics by deterministic and statistic studies so as to assess the risk of biomechanical failure on the different operative conditions the implant might face in a population of interest during various activities of daily living. Physiological conditions were defined including the variability of the anatomy, bone densitometry, surgery uncertainties and published boundary conditions at the hip. The protocol was tested by analysing a successful design on the market and a new prototype of a resurfacing prosthesis. Results. The intrinsic accuracy of models on bone stress predictions (RMSE < 10%) was aligned to the current state of the art in this field. The accuracy of prediction on the bone-prosthesis contact mechanics was also excellent (< 0.001 mm). The sensitivity of models prediction to uncertainties on modelling parameter was found below 8.4%. The analysis of the successful design resulted in a very good agreement with published retrospective studies. The geometry optimisation of the new prototype lead to a final design with a low risk of failure. The statistical analysis confirmed the minimal risk of the optimised design over the entire population of interest. The performances of the optimised design showed a significant improvement with respect to the first prototype (+35%). Limitations. On the authors opinion the major limitation of this study is on boundary conditions. The muscular forces and the hip joint reaction were derived from the few data available in the literature, which can be considered significant but hardly representative of the entire variability of boundary conditions the implant might face over the patients population. This moved the focus of the research on modelling the musculoskeletal system; the ongoing activity is to develop subject-specific musculoskeletal models of the lower limb from medical images. Conclusions. The developed protocol was able to accurately predict known clinical outcomes when applied to a well-established device and, to support the design optimisation phase providing important information on critical characteristics of the patients when applied to a new prosthesis. The presented approach does have a relevant generality that would allow the extension of the protocol to a large set of orthopaedic scenarios with minor changes. Hence, a failure mode analysis criterion can be considered a suitable tool in developing new orthopaedic devices.

Creep and damage models for the service behaviour of structural members

Deformability is often a crucial to the conception of many civil-engineering structural elements. Also, design is all the more burdensome if both long- and short-term deformability has to be considered. In this thesis, long- and short-term deformability has been studied from the material and the structural modelling point of view. Moreover, two materials have been handled: pultruded composites and concrete. A new finite element model for thin-walled beams has been introduced. As a main assumption, cross-sections rigid are considered rigid in their plane; this hypothesis replaces that of the classical beam theory of plane cross-sections in the deformed state. That also allows reducing the total number of degrees of freedom, and therefore making analysis faster compared with twodimensional finite elements. Longitudinal direction warping is left free, allowing describing phenomena such as the shear lag. The new finite-element model has been first applied to concrete thin-walled beams (such as roof high span girders or bridge girders) subject to instantaneous service loadings. Concrete in his cracked state has been considered through a smeared crack model for beams under bending. At a second stage, the FE-model has been extended to the viscoelastic field and applied to pultruded composite beams under sustained loadings. The generalized Maxwell model has been adopted. As far as materials are concerned, long-term creep tests have been carried out on pultruded specimens. Both tension and shear tests have been executed. Some specimen has been strengthened with carbon fibre plies to reduce short- and long- term deformability. Tests have been done in a climate room and specimens kept 2 years under constant load in time. As for concrete, a model for tertiary creep has been proposed. The basic idea is to couple the UMLV linear creep model with a damage model in order to describe nonlinearity. An effective strain tensor, weighting the total and the elasto-damaged strain tensors, controls damage evolution through the damage loading function. Creep strains are related to the effective stresses (defined by damage models) and so associated to the intact material.