Evaluating the influence of spatial uncertainty in locality points for species distributional modeling

1. Species distribution modelling is used increasingly in both applied and theoretical research to predict how species are distributed and to understand attributes of species' environmental requirements. In species distribution modelling, various statistical methods are used that combine species occurrence data with environmental spatial data layers to predict the suitability of any site for that species. While the number of data sharing initiatives involving species' occurrences in the scientific community has increased dramatically over the past few years, various data quality and methodological concerns related to using these data for species distribution modelling have not been addressed adequately. 2. We evaluated how uncertainty in georeferences and associated locational error in occurrences influence species distribution modelling using two treatments: (1) a control treatment where models were calibrated with original, accurate data and (2) an error treatment where data were first degraded spatially to simulate locational error. To incorporate error into the coordinates, we moved each coordinate with a random number drawn from the normal distribution with a mean of zero and a standard deviation of 5 km. We evaluated the influence of error on the performance of 10 commonly used distributional modelling techniques applied to 40 species in four distinct geographical regions. 3. Locational error in occurrences reduced model performance in three of these regions; relatively accurate predictions of species distributions were possible for most species, even with degraded occurrences. Two species distribution modelling techniques, boosted regression trees and maximum entropy, were the best performing models in the face of locational errors. The results obtained with boosted regression trees were only slightly degraded by errors in location, and the results obtained with the maximum entropy approach were not affected by such errors. 4. Synthesis and applications. To use the vast array of occurrence data that exists currently for research and management relating to the geographical ranges of species, modellers need to know the influence of locational error on model quality and whether some modelling techniques are particularly robust to error. We show that certain modelling techniques are particularly robust to a moderate level of locational error and that useful predictions of species distributions can be made even when occurrence data include some error.

Dynamical modeling and analysis of large cellular regulatory networks.

The dynamical analysis of large biological regulatory networks requires the development of scalable methods for mathematical modeling. Following the approach initially introduced by Thomas, we formalize the interactions between the components of a network in terms of discrete variables, functions, and parameters. Model simulations result in directed graphs, called state transition graphs. We are particularly interested in reachability properties and asymptotic behaviors, which correspond to terminal strongly connected components (or "attractors") in the state transition graph. A well-known problem is the exponential increase of the size of state transition graphs with the number of network components, in particular when using the biologically realistic asynchronous updating assumption. To address this problem, we have developed several complementary methods enabling the analysis of the behavior of large and complex logical models: (i) the definition of transition priority classes to simplify the dynamics; (ii) a model reduction method preserving essential dynamical properties, (iii) a novel algorithm to compact state transition graphs and directly generate compressed representations, emphasizing relevant transient and asymptotic dynamical properties. The power of an approach combining these different methods is demonstrated by applying them to a recent multilevel logical model for the network controlling CD4+ T helper cell response to antigen presentation and to a dozen cytokines. This model accounts for the differentiation of canonical Th1 and Th2 lymphocytes, as well as of inflammatory Th17 and regulatory T cells, along with many hybrid subtypes. All these methods have been implemented into the software GINsim, which enables the definition, the analysis, and the simulation of logical regulatory graphs.

DETERMINATION OF CELL-SPECIFIC NEUROTOXICITY OF MALONATE, METHYLMALONATE AND PROPIONATE IN A 3D RAT BRAIN CELL AGGREGATE SYSTEM

Malonate, methylmalonate and propionate are potentially neurotoxic metabolites in branched-chain organic acidurias. Their effects were tested on cultured 3D rat brain cell aggregates, using dosages of 0.1, 1.0 and 10.0 mM with a short but intense (twice a day over 3 days) and a longer but less intense treatment (every 3 rdday over 9 days). CNS cell-specific immunohistochemical stainings allowed the follow-up of neurons (axons, phosphorylated medium-weight neurofilament), astrocytes (glial fibrillary acidic protein) and oligodendrocytes (myelin basic protein). Methylmalonate and malonate were quantified by tandem mass spectrometry. Tandem mass spectrometry analysis of harvested brain cell aggregates revealed clear intracellular accumulation of methylmalonate and malonate. In immunohistochemical stainings oligodendrocytes appeared the most affected brain cells. The MBP signal disappeared already at 0.1 mM treatment with each metabolite. Mature astrocytes were not affected by propionate, while immature astrocytes on intense treatment with propionate developed cell swelling. 1 mM methylmalonate induced cell swelling of both immature and mature astrocytes , while 1 mM malonate only affected mature astrocytes. Neurons were not affected by methylmalonate, but 10.0 mM malonate on less intense treatment and 0.1, 1.0 and 10.0 mM propionate on intense treatment affected axonal growth. Our study shows significant uptake and deleterious effects of these metabolites on brain cells, principally on astrocytes and oligodendrocytes. This may be explained by the absence of the pathway in glial cells, which thus are not able to degrade these metabolites. Further studies are ongoing to elucidate the underlying mechanisms of the observed neurotoxic effects.

3D visualization software to analyze topological outcomes of topoisomerase reactions.

The action of various DNA topoisomerases frequently results in characteristic changes in DNA topology. Important information for understanding mechanistic details of action of these topoisomerases can be provided by investigating the knot types resulting from topoisomerase action on circular DNA forming a particular knot type. Depending on the topological bias of a given topoisomerase reaction, one observes different subsets of knotted products. To establish the character of topological bias, one needs to be aware of all possible topological outcomes of intersegmental passages occurring within a given knot type. However, it is not trivial to systematically enumerate topological outcomes of strand passage from a given knot type. We present here a 3D visualization software (TopoICE-X in KnotPlot) that incorporates topological analysis methods in order to visualize, for example, knots that can be obtained from a given knot by one intersegmental passage. The software has several other options for the topological analysis of mechanisms of action of various topoisomerases.

Multimodel inference and multimodel averaging in empirical modeling of occupational exposure levels.

Empirical modeling of exposure levels has been popular for identifying exposure determinants in occupational hygiene. Traditional data-driven methods used to choose a model on which to base inferences have typically not accounted for the uncertainty linked to the process of selecting the final model. Several new approaches propose making statistical inferences from a set of plausible models rather than from a single model regarded as 'best'. This paper introduces the multimodel averaging approach described in the monograph by Burnham and Anderson. In their approach, a set of plausible models are defined a priori by taking into account the sample size and previous knowledge of variables influent on exposure levels. The Akaike information criterion is then calculated to evaluate the relative support of the data for each model, expressed as Akaike weight, to be interpreted as the probability of the model being the best approximating model given the model set. The model weights can then be used to rank models, quantify the evidence favoring one over another, perform multimodel prediction, estimate the relative influence of the potential predictors and estimate multimodel-averaged effects of determinants. The whole approach is illustrated with the analysis of a data set of 1500 volatile organic compound exposure levels collected by the Institute for work and health (Lausanne, Switzerland) over 20 years, each concentration having been divided by the relevant Swiss occupational exposure limit and log-transformed before analysis. Multimodel inference represents a promising procedure for modeling exposure levels that incorporates the notion that several models can be supported by the data and permits to evaluate to a certain extent model selection uncertainty, which is seldom mentioned in current practice.

3D gait assessment in young and elderly subjects using foot-worn inertial sensors.

This study describes the validation of a new wearable system for assessment of 3D spatial parameters of gait. The new method is based on the detection of temporal parameters, coupled to optimized fusion and de-drifted integration of inertial signals. Composed of two wirelesses inertial modules attached on feet, the system provides stride length, stride velocity, foot clearance, and turning angle parameters at each gait cycle, based on the computation of 3D foot kinematics. Accuracy and precision of the proposed system were compared to an optical motion capture system as reference. Its repeatability across measurements (test-retest reliability) was also evaluated. Measurements were performed in 10 young (mean age 26.1±2.8 years) and 10 elderly volunteers (mean age 71.6±4.6 years) who were asked to perform U-shaped and 8-shaped walking trials, and then a 6-min walking test (6MWT). A total of 974 gait cycles were used to compare gait parameters with the reference system. Mean accuracy±precision was 1.5±6.8cm for stride length, 1.4±5.6cm/s for stride velocity, 1.9±2.0cm for foot clearance, and 1.6±6.1° for turning angle. Difference in gait performance was observed between young and elderly volunteers during the 6MWT particularly in foot clearance. The proposed method allows to analyze various aspects of gait, including turns, gait initiation and termination, or inter-cycle variability. The system is lightweight, easy to wear and use, and suitable for clinical application requiring objective evaluation of gait outside of the lab environment.

Coltop3D: A new software for structural analysis with high resolution 3D point clouds and DEM

Coltop3D is a software that performs structural analysis by using digital elevation model (DEM) and 3D point clouds acquired with terrestrial laser scanners. A color representation merging slope aspect and slope angle is used in order to obtain a unique code of color for each orientation of a local slope. Thus a continuous planar structure appears in a unique color. Several tools are included to create stereonets, to draw traces of discontinuities, or to compute automatically density stereonet. Examples are shown to demonstrate the efficiency of the method.

Cell Adhesion Peptide RGDSP and Laminin Support Proliferation and Migration of Postnatal Retinal Stem Cells in a 3D Hydrogel Culture System

Purpose: Retinal stem cells (RSCs) can be isolated from radial glia population of the newborn mouse retina (Angénieux et al., 2006). These RSCs have great capacity to renew and generate neurons including cells differentiated towards the photoreceptor lineage (Mehri-Soussi et al., 2006). However, our published results showed poor integration and survival rate after cell grafting into the retina. The uncontrollable environment of retina seems to be the problem. To bypass this, we are trying to generate hemi-retinal tissue in vitro that can be used for transplantation. Methods: Expanded RSCs were seeded in a mixture of poly-ethylene-glycol (PEG)-polymer-based hydrogels crosslinked by peptides that also serve as substrates for matrix metalloproteinases. Different doses of crosslinker peptides were tested. Several growth factors were studied to stimulate cell proliferation and differentiation. Results: Cells were trapped in hydrogels and cultured in the presence of FGF2 and EGF. Spherical cell clusters indicating proliferation appeared within several days, but there was no cell migration within the gel. We then added cell adhesion molecules integrin ligand RGDSP, or laminin, or a combination of both, into the gel. Cells grown with laminin showed the best proliferation. Cells grown with RGDSP proliferated a few times and then started to spread out. Cells grown with the combination of RGDSP and laminin showed better proliferation than with RGDSP alone and larger spread-outs than with laminin alone. After stimulations with first FGF2 and EGF, and then only FGF2, some cells showed neuronal morphology after 2 weeks. The neuronal population was assessed by the presence of neuronal marker b-tubulin-III. Glial cells were also present. Further characterizations are undergoing. Conclusions: RSC can grow and migrate in 3D hydrogel with the addition of FGF2, EGF, RGDSP and laminin. Further developments are necessary to form a homogenous tissue containing retinal cells.

The borehole-fluid effect in electrical resistivity imaging

Fluid that fills boreholes in crosswell electrical resistivity investigations provides the necessary electrical contact between the electrodes and the rock formation but it is also the source of image artifacts in standard inversions that do not account for the effects of the boreholes. The image distortions can be severe for large resistivity contrasts between the rock formation and borehole fluid and for large borehole diameters. We have carried out 3D finite-element modeling using an unstructured-grid approach to quantify the magnitude of borehole effects for different resistivity contrasts, borehole diameters, and electrode configurations. Relatively common resistivity contrasts of 100:1 and borehole diameters of 10 and 20 cm yielded, for a bipole length of 5 m, apparent resistivity underestimates of approximately 12% and 32% when using AB-MN configurations and apparent resistivity overestimates of approximately 24% and 95% when using AM-BN configurations. Effects are generally more severe at shorter bipole spacings. We report the results obtained by either including or ignoring the boreholes in inversions of 3D field data from a test site in Switzerland, where approximately 10,000 crosswell resistivity-tomography measurements were made across six acquisition planes among four boreholes. Inversions of raw data that ignored the boreholes filled with low-resistivity fluid paradoxically produced high-resistivity artifacts around the boreholes. Including correction factors based on the modeling results fora ID model with and without the boreholes did not markedly improve the images. The only satisfactory approach was to use a 3D inversion code that explicitly incorporated the boreholes in the actual inversion. This new approach yielded an electrical resistivity image that was devoid of artifacts around the boreholes and that correlated well with coincident crosswell radar images.

Evaluation of the 3D symmetry of the knees when walking using na ambulatory inertia system: results 14 months after ACL Plasty.

Purpose of the study: Reconstruction of the anterior cruciate ligament (ACL) controls laxity but does not enable restoration of strictly normal 3D kinematics. The purpose of this study was to compare the kinematics of the pathological knee with that of the healthy knee after ACL plasty. This study applied a new ambulatory system using miniature captors. Material and method: Five patients with an isolated injury of the ACL participated in this study. The patients were assessed after injury (T1), at five months (T2), and at 14 months (T3) after surgery. The assessment included laxity (KT-1000), the IKDC score and the Lysholm score. The 3D angles of the knees were measured when walking 30 m on flat ground using a system composed of to small inertia units (3D accelerometer and 3D gyroscope) and a portable recorder. Functional settings were optimised and validating to ensure easy precise measurement of the 3D angles. Symmetry of the two knees was quantified using a symmetry index (SI) (difference in amplitude normalised in relation to mean amplitude) and the correlation coefficient CC. Results: Clinical indicators improved during the follow-up (IKDC T1: 3C, 2C; T2: 5B; T3: 2A, 3B; subjective IKD: 53-95; Lysholm 67-96). Mean laxity improved from 8.6m to 2.5 mm. The gait analysis showed increased symmetry in terms of amplitude for flexion-extension (SI: −17% at T1, −1% at T2, 1% at T3), and an increase in symmetry in terms of the rotation signature (CC: 0.16 at T1, 0.99 at T2, 0.99 at T3). There was no trend to varus-valgus. Discussion: This study demonstrates the clinical application of the new ambulatory system for measuring 3D angles of the knee joint. Joint symmetry increased after ACL plasty but still showed some perturbation at 14 months. The results observed here are in agreement with the literature. Other patients and other types of gait are being analysed. Conclusion: This portable system allows gait analysis outside the laboratory, before and after ACL injury. It is very useful for follow-up after surgery.

Simulation des caractéristiques ECG du syndrome de Brugada à l'aide d'un programme de modélisation 3D

Introduction : Le syndrome de Brugada, décrit en 1992 par Pedro et Josep Brugada, est un syndrome cardiaque caractérisé par un sus-décalage particulier du segment ST associé à un bloc de branche droit atypique au niveau des dérivations ECG V1 à V3. Les altérations ECG du syndrome de Brugada sont classifiées en 3 types dont seul le type 1 est diagnostique. Les mécanismes physiopathologiques exacts de ce syndrome sont pour le moment encore controversés. Plusieurs hypothèses sont proposées dans la littérature dont deux principales retiennent l'attention : 1) le modèle du trouble de repolarisation stipule des potentiels d'action réduits en durée et en amplitude liés à un changement de répartition de canaux potassiques 2) le modèle du trouble de dépolarisation spécifie un retard de conduction se traduisant par une dépolarisation retardée. Dans le STEMI, un sus-décalage ST ressemblant à celui du syndrome de Brugada est expliqué par deux théories : 1) le courant de lésion diastolique suggère une élévation du potentiel diastolique transformé artificiellement en sus-décalage ST par les filtres utilisés dans tous les appareils ECG.¦Objectif : Recréer les manifestations ECG du syndrome de Brugada en appliquant les modifications du potentiel d'action des cardiomyocytes rapportées dans la littérature.¦Méthode : Pour ce travail, nous avons utilisé "ECGsim", un simulateur informatique réaliste d'ECG disponible gratuitement sur www.ecgsim.org. Ce programme est basé sur une reconstruction de l'ECG de surface à l'aide de 1500 noeuds représentant chacun les potentiels d'action des ventricules droit et gauche, épicardiques et endocardiques. L'ECG simulé peut être donc vu comme l'intégration de l'ensemble de ces potentiels d'action en tenant compte des propriétés de conductivité des tissus s'interposant entre les électrodes de surface et le coeur. Dans ce programme, nous avons définit trois zones, de taille différente, comprenant la chambre de chasse du ventricule droit. Pour chaque zone, nous avons reproduit les modifications des potentiels d'action citées dans les modèles du trouble de repolarisation et de dépolarisation et des théories de courant de lésion systolique et diastolique. Nous avons utilisé, en plus des douze dérivations habituelles, une électrode positionnée en V2IC3 (i.e. 3ème espace intercostal) sur le thorax virtuel du programme ECGsim.¦Résultats : Pour des raisons techniques, le modèle du trouble de repolarisation n'a pas pu être entièrement réalisée dans ce travail. Le modèle du trouble de dépolarisation ne reproduit pas d'altération de type Brugada mais un bloc de branche droit plus ou moins complet. Le courant de lésion diastolique permet d'obtenir un sus-décalage ST en augmentant le potentiel diastolique épicardique des cardiomyocytes de la chambre de chasse du ventricule droit. Une inversion de l'onde T apparaît lorsque la durée du potentiel d'action est prolongée. L'amplitude du sus-décalage ST dépend de la valeur du potentiel diastolique, de la taille de la lésion et de sa localisation épicardique ou transmurale. Le courant de lésion systolique n'entraîne pas de sus-décalage ST mais accentue l'amplitude de l'onde T.¦Discussion et conclusion : Dans ce travail, l'élévation du potentiel diastolique avec un prolongement de la durée du potentiel d'action est la combinaison qui reproduit le mieux les altérations ECG du Brugada. Une persistance de cellules de type nodal au niveau de la chambre de chasse du ventricule droit pourrait être une explication à ces modifications particulières du potentiel d'action. Le risque d'arythmie dans la Brugada pourrait également être expliqué par une automaticité anormale des cellules de type nodal. Ainsi, des altérations des mécanismes cellulaires impliqués dans le maintien du potentiel diastolique pourraient être présentes dans le syndrome de Brugada, ce qui, à notre connaissance, n'a jamais été rapporté dans la littérature.