Generalizing sample tree information with semiparametric and parametric models

Summary

Role of the inflammasome in murine experimental models of arthritis

Summary : The purpose of this study was to investigate the role of the inflammasome in human and experimental murine models (such as ΑΙΑ and K/BxN) of rheumatoid arthritis (RA)RA, affecting 1% of the population is the most frequent inflammatory disease characterized by synovial hyperplasia and cartilage and bone erosion, leading to joint destruction. In general, women are 3 times more affected by RA suggesting a role of estrogen in this disease. The inflammasome is a multiproteic complex triggering the activation of caspase-1 leading to the activation of IL-1 β, an important pro-inflammatory cytokine implicated in arthritis. The inflammasome has been implicated in several inflammatory diseases and particularly in gout. To highlight a possible role of the inflammasome in murine arthritis, we obtained ASC, caspase-1 and NALP3 +/+ and -/- littermate mice to perform ΑΙΑ and K/BxN arthritis. NALP3 -/- and caspase-1 -/- mice were as arthritic as wild type littermate mice in both ΑΙΑ and K/BxN models implicating that the NALP3 inflammasome is not involved in experimental arthritis. By contrast, ΑΙΑ severity was significantly diminished in ASC- deficient male and female mice, and in the K/BxN model, in ASC-deficient female mice. These results were supported by histological scoring and acute phase protein serum amyloid A (SAA) levels that were equivalent between NALP+/+ and NALP3-/- mice and diminished in ASC -/- mice. In ΑΙΑ and K/BxN murine experimental models, we observed a sexdependent phenotype. We studied the role of estradiol in both the ALA and the K/BxN models. Castrated female or male ASC -/- mice that received estradiol had a decreased arthritis severity. This implies a protective role of estrogen in the absence of ASC. In the ΑΙΑ model, proliferation assay were performed using splenocytes from mBSA- immunized ASC +/+ and -/- mice. The mBSA-induced proliferation was significantly lower in ASC-/- splenocytes. Moreover the CD3-specific proliferation of purified splenic Τ cells was significantly lower in ASC-/- cells. Finally, Τ cells from ASC-/- mice produced significantly decreased levels of IFN-gamma associated with increased levels of IL-10. These results imply a possible role of ASC in the TCR-signaling pathway and Τ cell cytokine production. In parallel the expression of the different inflammasome components were analyzed in biopsies from rheumatoid arthritis (RA) and osteoarthritis (OA) patiens. The expression of the 14 different NALPs, their effector protein ASC, and caspase-1 and -5 was readily measurable by RT-PCR in a similar proportion in RA and OA synovial samples, with the exception of NALP-5 and NALP-13, which weren't found in samples from either disease. The corresponding NALP1, -3, -12 and ASC proteins were expressed at similar levels in both OA and RA biopsies, as determined by immunohistochemistry and Western-blot analysis. By contrast, caspase-1 levels were significantly enhanced in RA synovial tissues compared to those from OA patients. NALP-1, -2, -3, -10, -12 and -14, as well as ASC, caspase-1, and -5 were detected in RNA from unstimulated and stimulated RA synoviocytes. In FLS, only ASC and caspase-1 were expressed at the protein level. NALP1, 3 and 12 were not detected. However, upon stimulation, no secreted IL-Ιβ was detectable in either RA or in OA synoviocytes culture medium. Résumé : Le but de ce projet était d'étudier le rôle de l'inflammasome dans des modèles expérimentaux d'arthrite tels que les modèles ΑΙΑ et K/BxN ainsi que dans la polyarthrite humaine (RA). La polyarthrite est une maladie inflammatoire très fréquente avec 1 % de la population affectée et touche 3 fois plus les femmes que les hommes, suggérant un rôle des hormones sexuelles dans cette pathologie. L'inflammasome est un complexe multiprotéique qui permet l'activation de la caspase-1, une cystéine protéase qui va ensuite cliver et activer rinterleukine-ΐβ (IL-Ιβ). L'inflammasome a été impliqué ces dernières années dans de nombreuses maladies inflammatoires notamment dans la goutte. Pour mettre en évidence un éventuel rôle de l'inflammasome dans l'arthrite expérimentale nous avons obtenu des souris déficientes pour certains des composants de l'inflammasome tels que ASC, NALP3 et caspase-1. Les souris NALP3 déficientes et caspase-1 déficientes sont aussi arthritiques que les souris wild type correspondantes que ce soit dans le modèle ΑΙΑ ou K/BxN. Par contre les souris mâles et femelles ASC-déficientes sont moins arthritiques que les souris +/+ correspondantes dans le modèle ΑΙΑ. Dans le modèle KRN, le même phénotype (diminution de la sévérité de l'arthrite) est observé uniquement chez les femelles ASC-/- Ce phénotype est corrélé avec l'histologie ainsi qu'avec le dosage du serum amyloid A (SAA) qui reflète l'inflammation systémique et qui est diminué chez les souris ASC-déficientes. Nous avons ensuite étudié le rôle de Γ estradiol (une des formes active des estrogènes) dans les modèles K/BxN et ΑΙΑ. Les souris castrées maies ou femelles déficientes pour ASC ayant reçu de l'estradiol ont une arthrite moins sévère ce qui implique que les estradiol ont un effet protecteur en l'absence de ASC. Dans le modèle ΑΙΑ, nous nous sommes aussi intéressés à la réponse immune. Des tests de prolifération ont été effectués sur des splénocytes en présence de mBSA (qui est l'antigène utilisé dans le modèle ΑΙΑ). Les splénocytes ASC -/- ont une proliferation qui est diminuée en présence de l'antigène. De plus la proliferation de cellules Τ spléniques purifiées en présence d'anti-CD3 est diminuée chez les cellules Τ ASC-/-. Ces résultats nous indiquent une éventuelle implication de ASC dans la signalisation par le récépteur des cellules T. En parallèle l'expression des différents composants de l'inflammasome a été analysée dans des biopsies de patients atteints de polyarthrite rhumatoide (RA) et d'arthrose (OA). L'expression des 14 différents NALPs, de l'adaptateur ASC, ainsi que des caspase-1 et -5 était similaires dans les échantillons RA et OA, à l'exception de NALP5 et 13 qui n'étaient pas détéctables. L'expression protéique de NALP1, 3, 12 et ASC effectuée par Western blot et immunohistochimie était similaire dans les biopsies RA et OA. Par contre la quantité de la caspase-1 mesurée par ELISA était augmentée de façon significative dans les extraits protéiques de biopsies RA. NALP-1, -2. -3, -10, -12, and -14 ainsi que ASC, caspase-1 et -5 étaient exprimés de façon similaire par les synoviocytes RA non stimulés et stimulés. Dans les synoviocytes seuls ASC et caspase-1 étaient détéctable au niveau protéique. NALP-1, -3 et -12 n'était pas détéctables. Cependant après stimulation il n'y avait d'IL-Ιβ sécrété que ce soit dans les surnageants de cultures de synoviocytes RA ou OA.

SBML qualitative models: a model representation format and infrastructure to foster interactions between qualitative modelling formalisms and tools.

BACKGROUND: Qualitative frameworks, especially those based on the logical discrete formalism, are increasingly used to model regulatory and signalling networks. A major advantage of these frameworks is that they do not require precise quantitative data, and that they are well-suited for studies of large networks. While numerous groups have developed specific computational tools that provide original methods to analyse qualitative models, a standard format to exchange qualitative models has been missing. RESULTS: We present the Systems Biology Markup Language (SBML) Qualitative Models Package ("qual"), an extension of the SBML Level 3 standard designed for computer representation of qualitative models of biological networks. We demonstrate the interoperability of models via SBML qual through the analysis of a specific signalling network by three independent software tools. Furthermore, the collective effort to define the SBML qual format paved the way for the development of LogicalModel, an open-source model library, which will facilitate the adoption of the format as well as the collaborative development of algorithms to analyse qualitative models. CONCLUSIONS: SBML qual allows the exchange of qualitative models among a number of complementary software tools. SBML qual has the potential to promote collaborative work on the development of novel computational approaches, as well as on the specification and the analysis of comprehensive qualitative models of regulatory and signalling networks.

Pattern recognition ecological niche models fit to presence-only and presence-absence data

1. Identifying the boundary of a species' niche from observational and environmental data is a common problem in ecology and conservation biology and a variety of techniques have been developed or applied to model niches and predict distributions. Here, we examine the performance of some pattern-recognition methods as ecological niche models (ENMs). Particularly, one-class pattern recognition is a flexible and seldom used methodology for modelling ecological niches and distributions from presence-only data. The development of one-class methods that perform comparably to two-class methods (for presence/absence data) would remove modelling decisions about sampling pseudo-absences or background data points when absence points are unavailable. 2. We studied nine methods for one-class classification and seven methods for two-class classification (five common to both), all primarily used in pattern recognition and therefore not common in species distribution and ecological niche modelling, across a set of 106 mountain plant species for which presence-absence data was available. We assessed accuracy using standard metrics and compared trade-offs in omission and commission errors between classification groups as well as effects of prevalence and spatial autocorrelation on accuracy. 3. One-class models fit to presence-only data were comparable to two-class models fit to presence-absence data when performance was evaluated with a measure weighting omission and commission errors equally. One-class models were superior for reducing omission errors (i.e. yielding higher sensitivity), and two-classes models were superior for reducing commission errors (i.e. yielding higher specificity). For these methods, spatial autocorrelation was only influential when prevalence was low. 4. These results differ from previous efforts to evaluate alternative modelling approaches to build ENM and are particularly noteworthy because data are from exhaustively sampled populations minimizing false absence records. Accurate, transferable models of species' ecological niches and distributions are needed to advance ecological research and are crucial for effective environmental planning and conservation; the pattern-recognition approaches studied here show good potential for future modelling studies. This study also provides an introduction to promising methods for ecological modelling inherited from the pattern-recognition discipline.

Selecting statistical models to study the relationship between soybean yield and soil physical properties

Statistical models allow the representation of data sets and the estimation and/or prediction of the behavior of a given variable through its interaction with the other variables involved in a phenomenon. Among other different statistical models, are the autoregressive state-space models (ARSS) and the linear regression models (LR), which allow the quantification of the relationships among soil-plant-atmosphere system variables. To compare the quality of the ARSS and LR models for the modeling of the relationships between soybean yield and soil physical properties, Akaike's Information Criterion, which provides a coefficient for the selection of the best model, was used in this study. The data sets were sampled in a Rhodic Acrudox soil, along a spatial transect with 84 points spaced 3 m apart. At each sampling point, soybean samples were collected for yield quantification. At the same site, soil penetration resistance was also measured and soil samples were collected to measure soil bulk density in the 0-0.10 m and 0.10-0.20 m layers. Results showed autocorrelation and a cross correlation structure of soybean yield and soil penetration resistance data. Soil bulk density data, however, were only autocorrelated in the 0-0.10 m layer and not cross correlated with soybean yield. The results showed the higher efficiency of the autoregressive space-state models in relation to the equivalent simple and multiple linear regression models using Akaike's Information Criterion. The resulting values were comparatively lower than the values obtained by the regression models, for all combinations of explanatory variables.

Predicting the growth response to thinning for Scots pine stands using individual-tree growth models

Summary

Calibration of the century, apsim and ndicea models of decomposition and n mineralization of plant residues in the humid tropics

The aim of this study was to calibrate the CENTURY, APSIM and NDICEA simulation models for estimating decomposition and N mineralization rates of plant organic materials (Arachis pintoi, Calopogonium mucunoides, Stizolobium aterrimum, Stylosanthes guyanensis) for 360 days in the Atlantic rainforest bioma of Brazil. The models´ default settings overestimated the decomposition and N-mineralization of plant residues, underlining the fact that the models must be calibrated for use under tropical conditions. For example, the APSIM model simulated the decomposition of the Stizolobium aterrimum and Calopogonium mucunoides residues with an error rate of 37.62 and 48.23 %, respectively, by comparison with the observed data, and was the least accurate model in the absence of calibration. At the default settings, the NDICEA model produced an error rate of 10.46 and 14.46 % and the CENTURY model, 21.42 and 31.84 %, respectively, for Stizolobium aterrimum and Calopogonium mucunoides residue decomposition. After calibration, the models showed a high level of accuracy in estimating decomposition and N- mineralization, with an error rate of less than 20 %. The calibrated NDICEA model showed the highest level of accuracy, followed by the APSIM and CENTURY. All models performed poorly in the first few months of decomposition and N-mineralization, indicating the need of an additional parameter for initial microorganism growth on the residues that would take the effect of leaching due to rainfall into account.

Nuclear curvature energy in relativistic models

The difficulties arising in the calculation of the nuclear curvature energy are analyzed in detail, especially with reference to relativistic models. It is underlined that the implicit dependence on curvature of the quantal wave functions is directly accessible only in a semiclassical framework. It is shown that also in the relativistic models quantal and semiclassical calculations of the curvature energy are in good agreement.

Nonideal particle distributions from kinetic freeze-out models

In fluid dynamical models the freeze-out of particles across a three-dimensional space-time hypersurface is discussed. The calculation of final momentum distribution of emitted particles is described for freeze-out surfaces, with both spacelike and timelike normals, taking into account conservation laws across the freeze-out discontinuity.

Freeze-out in hydrodynamical models

We study the effects of strict conservation laws and the problem of negative contributions to final momentum distribution during the freeze-out through 3-dimensional hypersurfaces with spacelike normal. We study some suggested solutions for this problem, and demonstrate it in one example.

Classical motions from pseudoclassical spin-1/2 particle models

The classical trajectory and spin precessions of Bargmann, Michel, and Telegdi are deduced from a pseudoclassical model of a relativistic spin-(1/2) particle. The corresponding deduction from a non- relativistic model is also given.

Do habitat suitability models reliably predict the recovery areas of threatened species?

1. Identifying those areas suitable for recolonization by threatened species is essential to support efficient conservation policies. Habitat suitability models (HSM) predict species' potential distributions, but the quality of their predictions should be carefully assessed when the species-environment equilibrium assumption is violated.2. We studied the Eurasian otter Lutra lutra, whose numbers are recovering in southern Italy. To produce widely applicable results, we chose standard HSM procedures and looked for the models' capacities in predicting the suitability of a recolonization area. We used two fieldwork datasets: presence-only data, used in the Ecological Niche Factor Analyses (ENFA), and presence-absence data, used in a Generalized Linear Model (GLM). In addition to cross-validation, we independently evaluated the models with data from a recolonization event, providing presences on a previously unoccupied river.3. Three of the models successfully predicted the suitability of the recolonization area, but the GLM built with data before the recolonization disagreed with these predictions, missing the recolonized river's suitability and badly describing the otter's niche. Our results highlighted three points of relevance to modelling practices: (1) absences may prevent the models from correctly identifying areas suitable for a species spread; (2) the selection of variables may lead to randomness in the predictions; and (3) the Area Under Curve (AUC), a commonly used validation index, was not well suited to the evaluation of model quality, whereas the Boyce Index (CBI), based on presence data only, better highlighted the models' fit to the recolonization observations.4. For species with unstable spatial distributions, presence-only models may work better than presence-absence methods in making reliable predictions of suitable areas for expansion. An iterative modelling process, using new occurrences from each step of the species spread, may also help in progressively reducing errors.5. Synthesis and applications. Conservation plans depend on reliable models of the species' suitable habitats. In non-equilibrium situations, such as the case for threatened or invasive species, models could be affected negatively by the inclusion of absence data when predicting the areas of potential expansion. Presence-only methods will here provide a better basis for productive conservation management practices.