Cell-based computational modeling of vascular morphogenesis using Tissue Simulation Toolkit.

Computational modeling has become a widely used tool for unraveling the mechanisms of higher level cooperative cell behavior during vascular morphogenesis. However, experimenting with published simulation models or adding new assumptions to those models can be daunting for novice and even for experienced computational scientists. Here, we present a step-by-step, practical tutorial for building cell-based simulations of vascular morphogenesis using the Tissue Simulation Toolkit (TST). The TST is a freely available, open-source C++ library for developing simulations with the two-dimensional cellular Potts model, a stochastic, agent-based framework to simulate collective cell behavior. We will show the basic use of the TST to simulate and experiment with published simulations of vascular network formation. Then, we will present step-by-step instructions and explanations for building a recent simulation model of tumor angiogenesis. Demonstrated mechanisms include cell-cell adhesion, chemotaxis, cell elongation, haptotaxis, and haptokinesis.

Cutting edge: apoptosis of superantigen-activated T cells occurs preferentially after a discrete number of cell divisions in vivo.

Staphylococcal enterotoxins are bacterial products that display superantigen activity in vitro as well as in vivo. For instance, staphylococcal enterotoxin B (SEB) polyclonally activates T cells that bear the Vbeta8 gene segment of the TCR. SEB-activated T cells undergo a burst of proliferation that is followed by apoptosis. Using an in vivo adaptation of a fluorescent cell division monitoring technique, we show here that SEB-activated T cells divide asynchronously, and that apoptosis of superantigen-activated T cells is preferentially restricted to cells which have undergone a discrete number of cell divisions. Collectively, our data suggest that superantigen-activated T cells are programmed to undergo a fixed number of cell divisions before undergoing apoptosis. A delayed death program may provide a mechanistic compromise between effector functions and homeostasis of activated T cells.

Regularity and mass conservation for discrete coagulation-fragmentation equations with diffusion

We present a new a-priori estimate for discrete coagulation fragmentation systems with size-dependent diffusion within a bounded, regular domain confined by homogeneous Neumann boundary conditions. Following from a duality argument, this a-priori estimate provides a global L2 bound on the mass density and was previously used, for instance, in the context of reaction-diffusion equations. In this paper we demonstrate two lines of applications for such an estimate: On the one hand, it enables to simplify parts of the known existence theory and allows to show existence of solutions for generalised models involving collision-induced, quadratic fragmentation terms for which the previous existence theory seems difficult to apply. On the other hand and most prominently, it proves mass conservation (and thus the absence of gelation) for almost all the coagulation coefficients for which mass conservation is known to hold true in the space homogeneous case.

Two weight inequalities for discrete positive operators

"Vegeu el resum a l'inici del document del fitxer adjunt."

The impact of past direct-personal traumatic events on 12-month outcome in first episode psychotic mania: trauma and early psychotic mania.

OBJECTIVE: Past traumatic events have been associated with poorer clinical outcomes in people with bipolar disorder. However, the impact of these events in the early stages of the illness remains unclear. The aim of this study was to investigate whether prior traumatic events were related to poorer outcomes 12 months following a first episode of psychotic mania. METHODS: Traumatic events were retrospectively evaluated from patient files in a sample of 65 participants who had experienced first episode psychotic mania. Participants were aged between 15 and 28 years and were treated at a specialised early psychosis service. Clinical outcomes were measured by a variety of symptomatic and functioning scales at the 12-month time-point. RESULTS: Direct-personal traumatic experiences prior to the onset of psychotic mania were reported by 48% of the sample. Participants with past direct-personal trauma had significantly higher symptoms of mania (p=0.02), depression (p=0.03) and psychopathology (p=0.01) 12 months following their first episode compared to participants without past direct-personal trauma, with medium to large effects observed. After adjusting for baseline scores, differences in global functioning (as measured by the Global Assessment of Functioning scale) were non-significant (p=0.05); however, participants with past direct-personal trauma had significantly poorer social and occupational functioning (p=0.04) at the 12-month assessment with medium effect. CONCLUSIONS: Past direct-personal trauma may predict poorer symptomatic and functional outcomes after first episode psychotic mania. Limitations include that the findings represent individuals treated at a specialist early intervention centre for youth and the retrospective assessment of traumatic events may have been underestimated.

Nash equilibrium strategies in discrete-time finite-horizon dynamic games with risk-and effort-averse players

The objective of this paper is to re-examine the risk-and effort attitude in the context of strategic dynamic interactions stated as a discrete-time finite-horizon Nash game. The analysis is based on the assumption that players are endogenously risk-and effort-averse. Each player is characterized by distinct risk-and effort-aversion types that are unknown to his opponent. The goal of the game is the optimal risk-and effort-sharing between the players. It generally depends on the individual strategies adopted and, implicitly, on the the players' types or characteristics.

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.

Agent based simulation to optimise emergency departments

Nowadays, many of the health care systems are large and complex environments and quite dynamic, specifically Emergency Departments, EDs. It is opened and working 24 hours per day throughout the year with limited resources, whereas it is overcrowded. Thus, is mandatory to simulate EDs to improve qualitatively and quantitatively their performance. This improvement can be achieved modelling and simulating EDs using Agent-Based Model, ABM and optimising many different staff scenarios. This work optimises the staff configuration of an ED. In order to do optimisation, objective functions to minimise or maximise have to be set. One of those objective functions is to find the best or optimum staff configuration that minimise patient waiting time. The staff configuration comprises: doctors, triage nurses, and admissions, the amount and sort of them. Staff configuration is a combinatorial problem, that can take a lot of time to be solved. HPC is used to run the experiments, and encouraging results were obtained. However, even with the basic ED used in this work the search space is very large, thus, when the problem size increases, it is going to need more resources of processing in order to obtain results in an acceptable time.

Evolution of anaesthesia care and related events between 1996 and 2010 in a developed country

Background: Anaesthesia Databank Switzerland (ADS) is a voluntary data registry introduced in 1996. The goal was to promote quality in anaesthesiology. Methods: Analysis of routinely recorded adverse events. Internal and external benchmark comparisons between anaesthesia departments. Results: In 2010, the database included 2'158'735 anaesthetic procedures. Forty-four anaesthesia departments were participating to the data collection in 2010. Over time, the number of patients in older age groups increased, the largest group being patients aged 50 to 64 years. Over time, the percentage of patients with ASA physical status score 1 decreased while the number of ASA 2 or 3 patients increased. The most frequent co-morbidities were hypertension (21%), smoking (16%), allergy (15%), and obesity (12%). Between 1996 and 2010, 146'459 adverse events were recorded, of which 34% were cardiovascular, 7% respiratory, 39% specific to anaesthesia and 17% nonspecific. The overall proportion of adverse events decreased over time, whatever their severity. Conclusion: The ADS routine data collection contributes to monitoring the trends of anaesthesia care in Switzerland.

Trace-element and phosphorus contents in sediment associated with Cretaceous oceanic anoxic events

ABSTRACT : During my SNSF-funded Ph.D. thesis project, I studied the evolution of redox conditions and organic-carbon preservation in the western Tethyan realm during three major positive excursions in the Cretaceous δ13C record, corresponding to the Valanginian, Early Aptian and Late Cenomanian. These periods were characterized by important global environmental and climate change, which was associated with perturbations in the carbon cycle. For the period of the Valanginian δ13C excursion, total organic carbon (TOC) contents and the quality of preserved organic matter are typical of oxic pelagic settings in the western Tethys. This is confirmed by the absence of major excursions in the stratigraphic distribution of RSTE during the δ13C shift. Published TOC data from other parts of the Valanginian oceans indicate that dys- to anaerobic zones were restricted to marginal seas within the Atlantic and Southern Ocean, and to the Pacific. Phosphorus (P) and mineralogical contents suggest a stepwise climatic evolution during the Valanginian, with a humid and warm climate prior to the δ13C shift leading to an increase in continental runoff. During the δ13C shift, a decrease in detrital input and P contents suggests a change in the climate towards more and conditions. During the early Aptian oceanic anoxic event (OAE 1a), a general increase followed by a rapid decrease in P contents suggests enhanced nutrient input at the beginning of OAE 1a. The return to lower values during OAE 1 a, associated with an increase in RSTE contents, may have been related to the weakened capacity to retain P in the sedimentary reservoir due to bottom-water oxygen depletion. In basinal settings, the RSTE distribution indicates well-developed anoxic conditions during OAE la, whereas in the shallower-water environments, conditions were oxic to suboxic, rather than anoxic. Furthermore, in the deeper part of the Tethys, two distinct enrichments have been observed, indicating fluctuations in the intensity of water column anoxia during the δ73C excursion. We also studied the effect of the end-Cenomanian oceanic anoxic event (OAE 2) on an expanded section in the Chrummflueschlucht (E of Euthal, Ct Switzerland). The goal here was to identify paleoceanographic and paleoenvironmental conditions during OAE 2 in this part of the northern Tethyan margin. The results show that this section is one of the most complete sections for the Cenomanian-Turonian boundary interval known from the Helvetic realm, despite a small hiatus between sediments corresponding to peaks 1 and 2 in the δ13C record. The evolution of P contents points to an increase in the input of this nutrient at the onset of OAE 2. The trends in RSTE contents show, however, that this part of the Helvetic realm was not affected by a strong depletion in oxygen conditions during OAE 2, despite its hemipelagic position. A further goal of this project was to submit the samples to a total extraction method (a combined HF/HNO3/HCI acid digestion) and compare the results obtained by the partial HNO3 acid extraction in order to standardize the analytical prócedures in the extraction of RSTE. The obtained results for samples of OAE 1 a suggest that RSTE trends using the partial HNO3 digestion are very comparable to those obtained by the total digestion method and subsequently normalized with regards to AI contents. RÉSUMÉ : Durant ce projet de thèse, financé par le Swiss National Science Funding (SNSF), j'ai étudié l'évolution des conditions redox et de la préservation de carbone organique dans le domnaine ouesttéthysien pendant trois excursions majeures du δ13C au Crétacé correspondant au Valanginien, à l'Aptien inférieur et à la limite Cénomanien-Turonien. Ces périodes sont caractérisées par des changements climatiques et environnementaux globaux associés à des perturbations dans le cylce du carbone. Pour L'excursion positive en δ13C du Valanginien, les analyses du carbone organique total (COT) et les observations palynologiques du domaine téthysien ont présenté des indications d'environnement pélagique relativementbienoxygéné. L'absence d'enrichissements en éléments traces sensibles aux conditions redox (TE) pendant l'excursion positive en δ13C confirme ces interprétations. Les données publiées de COT dans d'autres partie du globe indiquent cependant l'existence de conditions dys- à anaérobiques dans certains bassins restreints de l'Atlantique, l'Océan Austral et du Pacifique. L'évolution du phosphore (P) et la composition minéralogique des sédiments semblent indiquer un climat relativement chaud et humide avant l'excursion en δ13C entraînant une augmentation de l'altération continentale. Pendant le shift isotopique, une diminution des apports détritiques et du P suggèrent une transition vers des conditions plus arides. À l'Aptien Inférieur, le début de l'événement anoxique (OAE 1a) est marqué par une augmentation générale du P dans les sédiments indiquant une augmentation du niveau trophique à la base de l'excursion isotopique. Durant l'événement anoxique, les sédiments sont relativement appauvris en P. Cette diminution rapide associée à des enrichissements en TE est probablement liée à une remobilisation plus importante du P lors de la mise en place de conditions anoxiques dans les eaux de fond. Dans les environnements de bassin, le comportement des TE (enrichissements bien marqués) attestent de conditions réductrices bien marquées alors que dans les environnements moins profonds, les conditions semblent plutôt oxiques à dysoxiques. De plus, deux niveaux d'enrichissement en TE ont été observés dans la partie plus profonde de la Téthys, indiquant des fluctuations assez rapides dans l'intensité de l'anoxie de la colonne d'eau. Nous avons ensuite étudié les effets de l'événement anoxique de la fin du Cenomanien (OAE 2) dans un basin marginal de la marge nord de la Téthys avec la coupe de Chrummflueschlucht (à l'est de Euthal, Ct Schwyz). Les résultats ont montré que cette coupe présente un des enregistrements sédimentaires des plus complets de l'OAE 2 dans le domaine helvétique malgré un hiatus entre le pic 1 et 2 de l'excursion en δ13C. L'évolution du P montre une augmentation au début de l'OAE 2. Cependant, la distribution des TE indique que cette région n'a pas été affectée par des conditions réductrices trop importantes. Un second aspect de ce travail a été l'étude des différentes méthodes sur l'analyse de la distribution des TE. Des échantillons de l'OAE 1a ont été soumis à deux types d'extractions, l'une dite «totale » (attaque combinée d'acides HF/HNO3/HCI) et l'autre dite partielle » (HNO3). Les résultats obtenus suggèrent que les courbes de tendances des TE acquises par extraction partielle sont semblables à celle obtenues par extraction totale et normalisées par l'AI.

Power laws and scaling of rain events and dry spells in the Catalonia region

We analyze the statistics of rain-event sizes, rain-event durations, and dry-spell durations in a network of 20 rain gauges scattered in an area situated close to the NW Mediterranean coast. Power-law distributions emerge clearly for the dryspell durations, with an exponent around 1.50 ± 0.05, although for event sizes and durations the power-law ranges are rather limited, in some cases. Deviations from power-law behavior are attributed to finite-size effects. A scaling analysis helps to elucidate the situation, providing support for the existence of scale invariance in these distributions. It is remarkable that rain data of not very high resolution yield findings in agreement with self-organized critical phenomena.

Artificial Snow Optimization in Winter Sport Destinations Using a Multi-agent Simulation

This paper presents the Juste-Neige system for predicting the snow height on the ski runs of a resort using a multi-agent simulation software. Its aim is to facilitate snow cover management in order to i) reduce the production cost of artificial snow and to improve the profit margin for the companies managing the ski resorts; and ii) to reduce the water and energy consumption, and thus to reduce the environmental impact, by producing only the snow needed for a good skiing experience. The software provides maps with the predicted snow heights for up to 13 days. On these maps, the areas most exposed to snow erosion are highlighted. The software proceeds in three steps: i) interpolation of snow height measurements with a neural network; ii) local meteorological forecasts for every ski resort; iii) simulation of the impact caused by skiers using a multi-agent system. The software has been evaluated in the Swiss ski resort of Verbier and provides useful predictions.

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.