Simulation of stromatolite growth using diffusion-limited aggregation and a discrete model of biogeochemical exchanges in microbial mat

Resume : Mieux comprendre les stromatolithes et les tapis microbiens est un sujet important en biogéosciences puisque cela aide à l'étude des premières formes de vie sur Terre, a mieux cerner l'écologie des communautés microbiennes et la contribution des microorganismes a la biominéralisation, et même à poser certains fondements dans les recherches en exobiologie. D'autre part, la modélisation est un outil puissant utilisé dans les sciences naturelles pour appréhender différents phénomènes de façon théorique. Les modèles sont généralement construits sur un système d'équations différentielles et les résultats sont obtenus en résolvant ce système. Les logiciels disponibles pour implémenter les modèles incluent les logiciels mathématiques et les logiciels généraux de simulation. L'objectif principal de cette thèse est de développer des modèles et des logiciels pour aider a comprendre, via la simulation, le fonctionnement des stromatolithes et des tapis microbiens. Ces logiciels ont été développés en C++ en ne partant d'aucun pré-requis de façon a privilégier performance et flexibilité maximales. Cette démarche permet de construire des modèles bien plus spécifiques et plus appropriés aux phénomènes a modéliser. Premièrement, nous avons étudié la croissance et la morphologie des stromatolithes. Nous avons construit un modèle tridimensionnel fondé sur l'agrégation par diffusion limitée. Le modèle a été implémenté en deux applications C++: un moteur de simulation capable d'exécuter un batch de simulations et de produire des fichiers de résultats, et un outil de visualisation qui permet d'analyser les résultats en trois dimensions. Après avoir vérifié que ce modèle peut en effet reproduire la croissance et la morphologie de plusieurs types de stromatolithes, nous avons introduit un processus de sédimentation comme facteur externe. Ceci nous a mené a des résultats intéressants, et permis de soutenir l'hypothèse que la morphologie des stromatolithes pourrait être le résultat de facteurs externes autant que de facteurs internes. Ceci est important car la classification des stromatolithes est généralement fondée sur leur morphologie, imposant que la forme d'un stromatolithe est dépendante de facteurs internes uniquement (c'est-à-dire les tapis microbiens). Les résultats avancés dans ce mémoire contredisent donc ces assertions communément admises. Ensuite, nous avons décidé de mener des recherches plus en profondeur sur les aspects fonctionnels des tapis microbiens. Nous avons construit un modèle bidimensionnel de réaction-diffusion fondé sur la simulation discrète. Ce modèle a été implémenté dans une application C++ qui permet de paramétrer et exécuter des simulations. Nous avons ensuite pu comparer les résultats de simulation avec des données du monde réel et vérifier que le modèle peut en effet imiter le comportement de certains tapis microbiens. Ainsi, nous avons pu émettre et vérifier des hypothèses sur le fonctionnement de certains tapis microbiens pour nous aider à mieux en comprendre certains aspects, comme la dynamique des éléments, en particulier le soufre et l'oxygène. En conclusion, ce travail a abouti à l'écriture de logiciels dédiés à la simulation de tapis microbiens d'un point de vue tant morphologique que fonctionnel, suivant deux approches différentes, l'une holistique, l'autre plus analytique. Ces logiciels sont gratuits et diffusés sous licence GPL (General Public License). Abstract : Better understanding of stromatolites and microbial mats is an important topic in biogeosciences as it helps studying the early forms of life on Earth, provides clues re- garding the ecology of microbial ecosystems and their contribution to biomineralization, and gives basis to a new science, exobiology. On the other hand, modelling is a powerful tool used in natural sciences for the theoretical approach of various phenomena. Models are usually built on a system of differential equations and results are obtained by solving that system. Available software to implement models includes mathematical solvers and general simulation software. The main objective of this thesis is to develop models and software able to help to understand the functioning of stromatolites and microbial mats. Software was developed in C++ from scratch for maximum performance and flexibility. This allows to build models much more specific to a phenomenon rather than general software. First, we studied stromatolite growth and morphology. We built a three-dimensional model based on diffusion-limited aggregation. The model was implemented in two C++ applications: a simulator engine, which can run a batch of simulations and produce result files, and a Visualization tool, which allows results to be analysed in three dimensions. After verifying that our model can indeed reproduce the growth and morphology of several types of stromatolites, we introduced a sedimentation process as an external factor. This lead to interesting results, and allowed to emit the hypothesis that stromatolite morphology may be the result of external factors as much as internal factors. This is important as stromatolite classification is usually based on their morphology, imposing that a stromatolite shape is dependant on internal factors only (i.e. the microbial mat). This statement is contradicted by our findings, Second, we decided to investigate deeper the functioning of microbial mats, We built a two-dimensional reaction-diffusion model based on discrete simulation, The model was implemented in a C++ application that allows setting and running simulations. We could then compare simulation results with real world data and verify that our model can indeed mimic the behaviour of some microbial mats. Thus, we have proposed and verified hypotheses regarding microbial mats functioning in order to help to better understand them, e.g. the cycle of some elements such as oxygen or sulfur. ln conclusion, this PhD provides a simulation software, dealing with two different approaches. This software is free and available under a GPL licence.

An active strain electromechanical model for cardiac tissue

We propose a finite element approximation of a system of partial differential equations describing the coupling between the propagation of electrical potential and large deformations of the cardiac tissue. The underlying mathematical model is based on the active strain assumption, in which it is assumed that a multiplicative decomposition of the deformation tensor into a passive and active part holds, the latter carrying the information of the electrical potential propagation and anisotropy of the cardiac tissue into the equations of either incompressible or compressible nonlinear elasticity, governing the mechanical response of the biological material. In addition, by changing from an Eulerian to a Lagrangian configuration, the bidomain or monodomain equations modeling the evolution of the electrical propagation exhibit a nonlinear diffusion term. Piecewise quadratic finite elements are employed to approximate the displacements field, whereas for pressure, electrical potentials and ionic variables are approximated by piecewise linear elements. Various numerical tests performed with a parallel finite element code illustrate that the proposed model can capture some important features of the electromechanical coupling, and show that our numerical scheme is efficient and accurate.

Diverse set of Turing nanopatterns coat corneae across insect lineages.

Nipple-like nanostructures covering the corneal surfaces of moths, butterflies, and Drosophila have been studied by electron and atomic force microscopy, and their antireflective properties have been described. In contrast, corneal nanostructures of the majority of other insect orders have either been unexamined or examined by methods that did not allow precise morphological characterization. Here we provide a comprehensive analysis of corneal surfaces in 23 insect orders, revealing a rich diversity of insect corneal nanocoatings. These nanocoatings are categorized into four major morphological patterns and various transitions between them, many, to our knowledge, never described before. Remarkably, this unexpectedly diverse range of the corneal nanostructures replicates the complete set of Turing patterns, thus likely being a result of processes similar to those modeled by Alan Turing in his famous reaction-diffusion system. These findings reveal a beautiful diversity of insect corneal nanostructures and shed light on their molecular origin and evolutionary diversification. They may also be the first-ever biological example of Turing nanopatterns.

A fully adaptive numerical approximation for a two-dimensional epidemic model with nonlinear cross-diffusion

An epidemic model is formulated by a reactionâeuro"diffusion system where the spatial pattern formation is driven by cross-diffusion. The reaction terms describe the local dynamics of susceptible and infected species, whereas the diffusion terms account for the spatial distribution dynamics. For both self-diffusion and cross-diffusion, nonlinear constitutive assumptions are suggested. To simulate the pattern formation two finite volume formulations are proposed, which employ a conservative and a non-conservative discretization, respectively. An efficient simulation is obtained by a fully adaptive multiresolution strategy. Numerical examples illustrate the impact of the cross-diffusion on the pattern formation.

CD8 kinetically promotes ligand binding to the T-cell antigen receptor.

The mechanism of CD8 cooperation with the TCR in antigen recognition was studied on live T cells. Fluorescence correlation measurements yielded evidence of the presence of two TCR and CD8 subpopulations with different lateral diffusion rate constants. Independently, evidence for two subpopulations was derived from the experimentally observed two distinct association phases of cognate peptide bound to class I MHC (pMHC) tetramers and the T cells. The fast phase rate constant ((1.7 +/- 0.2) x 10(5) M(-1) s(-1)) was independent of examined cell type or MHC-bound peptides' structure. Its value was much faster than that of the association of soluble pMHC and TCR ((7.0 +/- 0.3) x 10(3) M(-1) s(-1)), and close to that of the association of soluble pMHC with CD8 ((1-2) x 10(5) M(-1) s(-1)). The fast binding phase disappeared when CD8-pMHC interaction was blocked by a CD8-specific mAb. The latter rate constant was slowed down approximately 10-fold after cells treatment with methyl-beta-cyclodextrin. These results suggest that the most efficient pMHC-cell association route corresponds to a fast tetramer binding to a colocalized CD8-TCR subpopulation, which apparently resides within membrane rafts: the reaction starts by pMHC association with the CD8. This markedly faster step significantly increases the probability of pMHC-TCR encounters and thereby promotes pMHC association with CD8-proximal TCR. The slow binding phase is assigned to pMHC association with a noncolocalized CD8-TCR subpopulation. Taken together with results of cytotoxicity assays, our data suggest that the colocalized, raft-associated CD8-TCR subpopulation is the one capable of inducing T-cell activation.

When the Line is Crossed... : Paths to Control and Sanction Behaviour Necessitating a State Reaction

The article presents a special form of a European comparative synopsis. For this case examples have been chosen ranging from administrative or minor (criminal) offences to increasingly serious offences and offenders. In this way it can be comparatively demonstrated how the criminal justice systems studied handle specific cases and whether they do so in a similar or different way.

Alveolar echinococcosis of the liver: Diffusion-weighted MRI findings and potential role in lesion characterisation

PURPOSE: To report the diffusion-weighted MRI findings in alveolar echinococcosis (AE) of the liver and evaluate the potential role of apparent diffusion coefficients (ADCs) in the characterisation of lesions. MATERIALS AND METHODS: We retrospectively included 22 patients with 63 AE liver lesions (≥1cm), examined with 3-T liver MRI, including a free-breathing diffusion-weighted single-shot echo-planar imaging sequence (b-values=50, 300 and 600s/mm(2)). Two radiologists jointly assessed the following lesion features: size, location, presence of cystic and/or solid components (according to Kodama's classification system), relative contrast enhancement, and calcifications (on CT). The ADCtotal, ADCmin and ADCmax were measured in each lesion and the surrounding liver parenchyma. RESULTS: Three type 1, 19 type 2, 17 type 3, three type 4 and 21 type 5 lesions were identified. The mean (±SD) ADCtotal, ADCmin and ADCmax for all lesions were 1.73±0.50, 0.76±0.38 and 2.63±0.76×10(-3)mm(2)/s, respectively. The mean ADCtotal for type 1, type 2, type 3, type 4 and type 5 lesions were 1.97±1.01, 1.76±0.53, 1.73±0.41, 1.15±0.42 and 1.76±0.44×10(-3)mm(2)/s, respectively. No significant differences were found between the five lesion types, except for type 4 (p=0.0363). There was a significant correlation between the presence of a solid component and low ADCmin (r=0.39, p=0.0016), whereas an inverse correlation was found between the relative contrast enhancement and ADCtotal (r=-0.34, p=0.0072). CONCLUSION: The ADCs of AE lesions are relatively low compared to other cystic liver lesions, which may help in the differential diagnosis. Although ADCs are of little use to distinguish between the five lesion types, their low value reflects the underlying solid component.

A behavioral study of alpha-1b adrenergic receptor knockout mice: increased reaction to novelty and selectively reduced learning capacities.

Knockout mice lacking the alpha-1b adrenergic receptor were tested in behavioral experiments. Reaction to novelty was first assessed in a simple test in which the time taken by the knockout mice and their littermate controls to enter a second compartment was compared. Then the mice were tested in an open field to which unknown objects were subsequently added. Special novelty was introduced by moving one of the familiar objects to another location in the open field. Spatial behavior and memory were further studied in a homing board test, and in the water maze. The alpha-1b knockout mice showed an enhanced reactivity to new situations. They were faster to enter the new environment, covered longer paths in the open field, and spent more time exploring the new objects. They reacted like controls to modification inducing spatial novelty. In the homing board test, both the knockout mice and the control mice seemed to use a combination of distant visual and proximal olfactory cues, showing place preference only if the two types of cues were redundant. In the water maze the alpha-1b knockout mice were unable to learn the task, which was confirmed in a probe trial without platform. They were perfectly able, however, to escape in a visible platform procedure. These results confirm previous findings showing that the noradrenergic pathway is important for the modulation of behaviors such as reaction to novelty and exploration, and suggest that this is mediated, at least partly, through the alpha-1b adrenergic receptors. The lack of alpha-1b adrenergic receptors in spatial orientation does not seem important in cue-rich tasks but may interfere with orientation in situations providing distant cues only.

Investigation of field and diffusion time dependence of the diffusion-weighted signal at ultrahigh magnetic fields.

Over the last decade, there has been a significant increase in the number of high-magnetic-field MRI magnets. However, the exact effect of a high magnetic field strength (B0 ) on diffusion-weighted MR signals is not yet fully understood. The goal of this study was to investigate the influence of different high magnetic field strengths (9.4 T and 14.1 T) and diffusion times (9, 11, 13, 15, 17 and 24 ms) on the diffusion-weighted signal in rat brain white matter. At a short diffusion time (9 ms), fractional anisotropy values were found to be lower at 14.1 T than at 9.4 T, but this difference disappeared at longer diffusion times. A simple two-pool model was used to explain these findings. The model describes the white matter as a first hindered compartment (often associated with the extra-axonal space), characterized by a faster orthogonal diffusion and a lower fractional anisotropy, and a second restricted compartment (often associated with the intra-axonal space), characterized by a slower orthogonal diffusion (i.e. orthogonal to the axon direction) and a higher fractional anisotropy. Apparent T2 relaxation time measurements of the hindered and restricted pools were performed. The shortening of the pseudo-T2 value from the restricted compartment with B0 is likely to be more pronounced than the apparent T2 changes in the hindered compartment. This study suggests that the observed differences in diffusion tensor imaging parameters between the two magnetic field strengths at short diffusion time may be related to differences in the apparent T2 values between the pools. Copyright © 2013 John Wiley & Sons, Ltd.

Measuring the diffusion of palliative care in long-term care facilities - a death census

ABSTRACT: BACKGROUND: The dissemination of palliative care for patients presenting complex chronic diseases at various stages has become an important matter of public health. A death census in Swiss long-term care facilities (LTC) was set up with the aim of monitoring the frequency of selected indicators of palliative care. METHODS: The survey covered 150 LTC facilities (105 nursing homes and 45 home health services), each of which was asked to complete a questionnaire for every non-accidental death over a period of six months. The frequency of 4 selected indicators of palliative care (resort to a specialized palliative care service, the administration of opiates, use of any pain measurement scale or other symptom measurement scale) was monitored in respect of the stages of care and analysed based on gender, age, medical condition and place of residence. RESULTS: Overall, 1200 deaths were reported, 29.1% of which were related to cancer. The frequencies of each indicator varied according to the type of LTC, mostly regarding the administration of opiate. It appeared that the access to palliative care remained associated with cancer, terminal care and partly with age, whereas gender and the presence of mental disorders had no effect on the indicators. In addition, the use of drugs was much more frequent than the other indicators. CONCLUSION: The profile of patients with access to palliative care must become more diversified. Among other recommendations, equal access to opiates in nursing homes and in home health services, palliative care at an earlier stage and the systematic use of symptom management scales when resorting to opiates have to become of prime concern.