Normal forms, stability and splitting of invariant manifolds II. Finitely differentiable Hamiltonians

This paper is a sequel to ``Normal forms, stability and splitting of invariant manifolds I. Gevrey Hamiltonians", in which we gave a new construction of resonant normal forms with an exponentially small remainder for near-integrable Gevrey Hamiltonians at a quasi-periodic frequency, using a method of periodic approximations. In this second part we focus on finitely differentiable Hamiltonians, and we derive normal forms with a polynomially small remainder. As applications, we obtain a polynomially large upper bound on the stability time for the evolution of the action variables and a polynomially small upper bound on the splitting of invariant manifolds for hyperbolic tori.

Assessment of liver tumor response by high-field (3T) MRI after radiofrequency ablation: Short- and mid-term evolution of diffusion parameters within the ablation zone.

PURPOSE: To compare the apparent diffusion coefficient (ADC) values of malignant liver lesions on diffusion-weighted MRI (DWI) before and after successful radiofrequency ablation (RF ablation). MATERIALS AND METHODS: Thirty-two patients with 43 malignant liver lesions (23/20: metastases/hepatocellular carcinomas (HCC)) underwent liver MRI (3.0T) before (<1month) and after RF ablation (at 1, 3 and 6months) using T2-, gadolinium-enhanced T1- and DWI-weighted MR sequences. Jointly, two radiologists prospectively measured ADCs for each lesion by means of two different regions of interest (ROIs), first including the whole lesion and secondly the area with the visibly most restricted diffusion (MRDA) on ADC map. Changes of ADCs were evaluated with ANOVA and Dunnett tests. RESULTS: Thirty-one patients were successfully treated, while one patient was excluded due to focal recurrence. In metastases (n=22), the ADC in the whole lesion and in MRDA showed an up-and-down evolution. In HCC (n=20), the evolution of ADC was more complex, but with significantly higher values (p=0.013) at 1 and 6months after RF ablation. CONCLUSION: The ADC values of malignant liver lesions successfully treated by RF ablation show a predictable evolution and may help radiologists to monitor tumor response after treatment.

A New Imaging Software To Evaluate The Centration Of Capsulorhexis And Intraocular Lens (IOL) After Cataract Surgery

Purpose: IOL centration and stability after cataract surgery is of high interest for cataract surgeons and IOL-producing companies. We present a new imaging software to evaluate the centration of the rhexis and the centration of the IOL after cataract surgery.Methods: We developed, in collaboration with the Biomedical Imaging Group (BIG), EPFL, Lausanne, a new working tool in order to assess precisely outcomes after IOL-implantation, such as ideal capsulorhexis and IOL-centration. The software is a plug-in of ImageJ, a general-purpose image processing and image-analysis package. The specifications of this software are: evaluation of the rhexis-centration and evaluation the position of the IOL in the posterior chamber. The end points are to analyze the quality of the centration of a rhexis after cataract surgery, the deformation of the rhexis with capsular bag retraction and the centration of the IOL after implantation.Results: This software delivers tools to interactively measure the distances between limbus, IOL and capsulorhexis and its changes over time. The user is invited to adjust nodes of three radial curves for the limbus, rhexis and the optic of the IOL. The radial distances of the curves are computed to evaluate the IOL implantation. The user is also able to define patterns for ideal capsulorhexis and optimal IOL-centration. We are going to present examples of calculations after cataract surgery.Conclusions: Evaluation of the centration of the rhexis and of the IOL after cataract surgery is an important end point for optimal IOL implantation after cataract surgery. Especially multifocal or accommodative lenses need a precise position in the bag with a good stability over time. This software is able to evaluate these parameters just after the surgery but also its changes over time. The results of these evaluations can lead to an optimizing of surgical procedures and materials.

Cleaning wrasse species vary with respect to dependency on the mutualism and behavioural adaptations in interactions

Interspecific mutualisms are an essential feature of life on earth, yet we know little about their evolution and stability. In many mutualisms several species are available as partners, raising questions about the similarity in function and behavioural repertoire depending on the partner species. Furthermore, variation between species in the quantity and quality of interactions resulting in variation in payoffs may allow us to infer the potential evolutionary origin of a multispecies mutualism complex. We addressed these issues in the marine cleaning mutualism, in which so-called 'cleaners' remove ectoparasites from so-called 'client' reef fish. We measured several parameters concerning the quantity and quality of cleaning interactions in six sympatric cleaner wrasse species. We found significant variation between cleaner species with respect to client diversity, the number of interactions with predatory clients, the duration of interactions, the frequency of client jolts as a correlate of 'cheating' by cleaners, and behaviours used for manipulation of client decisions. Exploratory correlations between cleaner species' dependency and our variables of interest suggest that cleaning originated as a conflict-free by-product mutualism and evolved towards more sophisticated behaviours, including strategic behaviours for interactions with predators, cheating and manipulation specifically adapted to the client type.

Sex-biased dispersal in a migratory bat: a characterization using sex-specific demographic parameters.

We studied the noctule bat (Nyctalus noctula), in which the mitochondrial F(ST) is about 10 times that revealed by nuclear markers, to address two questions. We first verified whether random dispersal of one sex is compatible with highly contrasted mitochondrial and nuclear population structures. Using computer simulations, we then assessed the power of multilocus population differentiation tests when the expected population structure departs only slightly from panmixia. Using an island model with sex-specific demographic parameters, we found that random male dispersal is consistent with the population structure observed in the noctule. However, other parameter combinations are also compatible with the data. We computed the minimum sex bias in dispersal (at least 69% of the dispersing individuals are males), a result that would not be available if we had used more classical population genetic models. The power of multilocus population differentiation tests was unexpectedly high, the tests being significant in almost 100% of the replicates, although the observed population structure infered from nuclear markers was extremely low (F(ST) = 0.6%).

Stochastic stability and the evolution of coordination in spatially structured populations.

Animals can often coordinate their actions to achieve mutually beneficial outcomes. However, this can result in a social dilemma when uncertainty about the behavior of partners creates multiple fitness peaks. Strategies that minimize risk ("risk dominant") instead of maximizing reward ("payoff dominant") are favored in economic models when individuals learn behaviors that increase their payoffs. Specifically, such strategies are shown to be "stochastically stable" (a refinement of evolutionary stability). Here, we extend the notion of stochastic stability to biological models of continuous phenotypes at a mutation-selection-drift balance. This allows us to make a unique prediction for long-term evolution in games with multiple equilibria. We show how genetic relatedness due to limited dispersal and scaled to account for local competition can crucially affect the stochastically-stable outcome of coordination games. We find that positive relatedness (weak local competition) increases the chance the payoff dominant strategy is stochastically stable, even when it is not risk dominant. Conversely, negative relatedness (strong local competition) increases the chance that strategies evolve that are neither payoff nor risk dominant. Extending our results to large multiplayer coordination games we find that negative relatedness can create competition so extreme that the game effectively changes to a hawk-dove game and a stochastically stable polymorphism between the alternative strategies evolves. These results demonstrate the usefulness of stochastic stability in characterizing long-term evolution of continuous phenotypes: the outcomes of multiplayer games can be reduced to the generic equilibria of two-player games and the effect of spatial structure can be analyzed readily.

Effect of oxidized cholesterol on age-associated changes to immune parameters in spleen lymphocytes and peritoneal exudate cells derived from rats

The effects of oxidized cholesterol on immune parameters were examined by using spleen lymphocytes and peritoneal exudate cells (PEC) derived from 5-week- (Young) and 9-month-old (Adult) rats. The immunoglobulin (Ig) G and IgM production was inhibited by oxidized cholesterol in the rats of both ages when lymphocytes were exposed to 30 micrograms/ml of oxidized cholesterol for 24 hr. The intracellular IgA level was also lowered by 30 micrograms/ml of oxidized cholesterol, irrespective of age. In contrast, IgE production was significantly increased by the addition of 30 micrograms/ml of oxidized cholesterol in only young lymphocytes. Moreover, oxidized cholesterol enhanced the intracellular histamine accumulation in only adult PEC, although the total histamine level produced by PEC was similar in the rats of both ages. These results thus suggest the possibility that oxidized cholesterol can have different effects on the age-related modulation of immune functions such as Igs production and histamine release.

Dose reconstruction in the context of tomotherapy

In vivo dosimetry is a way to verify the radiation dose delivered to the patient in measuring the dose generally during the first fraction of the treatment. It is the only dose delivery control based on a measurement performed during the treatment. In today's radiotherapy practice, the dose delivered to the patient is planned using 3D dose calculation algorithms and volumetric images representing the patient. Due to the high accuracy and precision necessary in radiation treatments, national and international organisations like ICRU and AAPM recommend the use of in vivo dosimetry. It is also mandatory in some countries like France. Various in vivo dosimetry methods have been developed during the past years. These methods are point-, line-, plane- or 3D dose controls. A 3D in vivo dosimetry provides the most information about the dose delivered to the patient, with respect to ID and 2D methods. However, to our knowledge, it is generally not routinely applied to patient treatments yet. The aim of this PhD thesis was to determine whether it is possible to reconstruct the 3D delivered dose using transmitted beam measurements in the context of narrow beams. An iterative dose reconstruction method has been described and implemented. The iterative algorithm includes a simple 3D dose calculation algorithm based on the convolution/superposition principle. The methodology was applied to narrow beams produced by a conventional 6 MV linac. The transmitted dose was measured using an array of ion chambers, as to simulate the linear nature of a tomotherapy detector. We showed that the iterative algorithm converges quickly and reconstructs the dose within a good agreement (at least 3% / 3 mm locally), which is inside the 5% recommended by the ICRU. Moreover it was demonstrated on phantom measurements that the proposed method allows us detecting some set-up errors and interfraction geometry modifications. We also have discussed the limitations of the 3D dose reconstruction for dose delivery error detection. Afterwards, stability tests of the tomotherapy MVCT built-in onboard detector was performed in order to evaluate if such a detector is suitable for 3D in-vivo dosimetry. The detector showed stability on short and long terms comparable to other imaging devices as the EPIDs, also used for in vivo dosimetry. Subsequently, a methodology for the dose reconstruction using the tomotherapy MVCT detector is proposed in the context of static irradiations. This manuscript is composed of two articles and a script providing further information related to this work. In the latter, the first chapter introduces the state-of-the-art of in vivo dosimetry and adaptive radiotherapy, and explains why we are interested in performing 3D dose reconstructions. In chapter 2 a dose calculation algorithm implemented for this work is reviewed with a detailed description of the physical parameters needed for calculating 3D absorbed dose distributions. The tomotherapy MVCT detector used for transit measurements and its characteristics are described in chapter 3. Chapter 4 contains a first article entitled '3D dose reconstruction for narrow beams using ion chamber array measurements', which describes the dose reconstruction method and presents tests of the methodology on phantoms irradiated with 6 MV narrow photon beams. Chapter 5 contains a second article 'Stability of the Helical TomoTherapy HiArt II detector for treatment beam irradiations. A dose reconstruction process specific to the use of the tomotherapy MVCT detector is presented in chapter 6. A discussion and perspectives of the PhD thesis are presented in chapter 7, followed by a conclusion in chapter 8. The tomotherapy treatment device is described in appendix 1 and an overview of 3D conformai- and intensity modulated radiotherapy is presented in appendix 2. - La dosimétrie in vivo est une technique utilisée pour vérifier la dose délivrée au patient en faisant une mesure, généralement pendant la première séance du traitement. Il s'agit de la seule technique de contrôle de la dose délivrée basée sur une mesure réalisée durant l'irradiation du patient. La dose au patient est calculée au moyen d'algorithmes 3D utilisant des images volumétriques du patient. En raison de la haute précision nécessaire lors des traitements de radiothérapie, des organismes nationaux et internationaux tels que l'ICRU et l'AAPM recommandent l'utilisation de la dosimétrie in vivo, qui est devenue obligatoire dans certains pays dont la France. Diverses méthodes de dosimétrie in vivo existent. Elles peuvent être classées en dosimétrie ponctuelle, planaire ou tridimensionnelle. La dosimétrie 3D est celle qui fournit le plus d'information sur la dose délivrée. Cependant, à notre connaissance, elle n'est généralement pas appliquée dans la routine clinique. Le but de cette recherche était de déterminer s'il est possible de reconstruire la dose 3D délivrée en se basant sur des mesures de la dose transmise, dans le contexte des faisceaux étroits. Une méthode itérative de reconstruction de la dose a été décrite et implémentée. L'algorithme itératif contient un algorithme simple basé sur le principe de convolution/superposition pour le calcul de la dose. La dose transmise a été mesurée à l'aide d'une série de chambres à ionisations alignées afin de simuler la nature linéaire du détecteur de la tomothérapie. Nous avons montré que l'algorithme itératif converge rapidement et qu'il permet de reconstruire la dose délivrée avec une bonne précision (au moins 3 % localement / 3 mm). De plus, nous avons démontré que cette méthode permet de détecter certaines erreurs de positionnement du patient, ainsi que des modifications géométriques qui peuvent subvenir entre les séances de traitement. Nous avons discuté les limites de cette méthode pour la détection de certaines erreurs d'irradiation. Par la suite, des tests de stabilité du détecteur MVCT intégré à la tomothérapie ont été effectués, dans le but de déterminer si ce dernier peut être utilisé pour la dosimétrie in vivo. Ce détecteur a démontré une stabilité à court et à long terme comparable à d'autres détecteurs tels que les EPIDs également utilisés pour l'imagerie et la dosimétrie in vivo. Pour finir, une adaptation de la méthode de reconstruction de la dose a été proposée afin de pouvoir l'implémenter sur une installation de tomothérapie. Ce manuscrit est composé de deux articles et d'un script contenant des informations supplémentaires sur ce travail. Dans ce dernier, le premier chapitre introduit l'état de l'art de la dosimétrie in vivo et de la radiothérapie adaptative, et explique pourquoi nous nous intéressons à la reconstruction 3D de la dose délivrée. Dans le chapitre 2, l'algorithme 3D de calcul de dose implémenté pour ce travail est décrit, ainsi que les paramètres physiques principaux nécessaires pour le calcul de dose. Les caractéristiques du détecteur MVCT de la tomothérapie utilisé pour les mesures de transit sont décrites dans le chapitre 3. Le chapitre 4 contient un premier article intitulé '3D dose reconstruction for narrow beams using ion chamber array measurements', qui décrit la méthode de reconstruction et présente des tests de la méthodologie sur des fantômes irradiés avec des faisceaux étroits. Le chapitre 5 contient un second article intitulé 'Stability of the Helical TomoTherapy HiArt II detector for treatment beam irradiations'. Un procédé de reconstruction de la dose spécifique pour l'utilisation du détecteur MVCT de la tomothérapie est présenté au chapitre 6. Une discussion et les perspectives de la thèse de doctorat sont présentées au chapitre 7, suivies par une conclusion au chapitre 8. Le concept de la tomothérapie est exposé dans l'annexe 1. Pour finir, la radiothérapie «informationnelle 3D et la radiothérapie par modulation d'intensité sont présentées dans l'annexe 2.

A repeated GGA motif is critical for the activity and stability of the riboregulator RsmY of Pseudomonas fluorescens.

The riboregulator RsmY of Pseudomonas fluorescens strain CHA0 is an example of small regulatory RNAs belonging to the global Rsm/Csr regulatory systems controlling diverse cellular processes such as glycogen accumulation, motility, or formation of extracellular products in various bacteria. By binding multiple molecules of the small regulatory protein RsmA, RsmY relieves the negative effect of RsmA on the translation of several target genes involved in the biocontrol properties of strain CHA0. RsmY and functionally related riboregulators have repeated GGA motifs predicted to be exposed in single-stranded regions, notably in the loops of hairpins. The secondary structure of RsmY was corroborated by in vivo cleavage with lead acetate. RsmY mutants lacking three or five (out of six) of the GGA motifs showed reduced ability to derepress the expression of target genes in vivo and failed to bind the RsmA protein efficiently in vitro. The absence of GGA motifs in RsmY mutants resulted in reduced abundance of these transcripts and in a shorter half-life (< or = 6 min as compared with 27 min for wild type RsmY). These results suggest that both the interaction of RsmY with RsmA and the stability of RsmY strongly depend on the GGA repeats and that the ability of RsmY to interact with small regulatory proteins such as RsmA may protect this RNA from degradation.

Simulation methods with extended stability for stiff biochemical kinetics

Background: With increasing computer power, simulating the dynamics of complex systems in chemistry and biology is becoming increasingly routine. The modelling of individual reactions in (bio)chemical systems involves a large number of random events that can be simulated by the stochastic simulation algorithm (SSA). The key quantity is the step size, or waiting time, τ, whose value inversely depends on the size of the propensities of the different channel reactions and which needs to be re-evaluated after every firing event. Such a discrete event simulation may be extremely expensive, in particular for stiff systems where τ can be very short due to the fast kinetics of some of the channel reactions. Several alternative methods have been put forward to increase the integration step size. The so-called τ-leap approach takes a larger step size by allowing all the reactions to fire, from a Poisson or Binomial distribution, within that step. Although the expected value for the different species in the reactive system is maintained with respect to more precise methods, the variance at steady state can suffer from large errors as τ grows. Results: In this paper we extend Poisson τ-leap methods to a general class of Runge-Kutta (RK) τ-leap methods. We show that with the proper selection of the coefficients, the variance of the extended τ-leap can be well-behaved, leading to significantly larger step sizes.Conclusions: The benefit of adapting the extended method to the use of RK frameworks is clear in terms of speed of calculation, as the number of evaluations of the Poisson distribution is still one set per time step, as in the original τ-leap method. The approach paves the way to explore new multiscale methods to simulate (bio)chemical systems.

The control of microtubule stability in vitro and in transfected cells by MAP1B and SCG10.

In neurons, the regulation of microtubules plays an important role for neurite outgrowth, axonal elongation, and growth cone steering. SCG10 family proteins are the only known neuronal proteins that have a strong destabilizing effect, are highly enriched in growth cones and are thought to play an important role during axonal elongation. MAP1B, a microtubule-stabilizing protein, is found in growth cones as well, therefore it was important to test their effect on microtubules in the presence of both proteins. We used recombinant proteins in microtubule assembly assays and in transfected COS-7 cells to analyze their combined effects in vitro and in living cells, respectively. Individually, both proteins showed their expected activities in microtubule stabilization and destruction respectively. In MAP1B/SCG10 double-transfected cells, MAP1B could not protect microtubules from SCG10-induced disassembly in most cells, in particular not in cells that contained high levels of SCG10. This suggests that SCG10 is more potent to destabilize microtubules than MAP1B to rescue them. In microtubule assembly assays, MAP1B promoted microtubule formation at a ratio of 1 MAP1B per 70 tubulin dimers while a ratio of 1 SCG10 per two tubulin dimers was needed to destroy microtubules. In addition to its known binding to tubulin dimers, SCG10 binds also to purified microtubules in growth cones of dorsal root ganglion neurons in culture. In conclusion, neuronal microtubules are regulated by antagonistic effects of MAP1B and SCG10 and a fine tuning of the balance of these proteins may be critical for the regulation of microtubule dynamics in growth cones.

Review of Stability Berm Alternatives for Environmentally Senstitive Areas, 2005

Stability berms are commonly constructed where roadway embankments cross soft or unstable ground conditions. Under certain circumstances, the construction of stability berms cause unfavorable environmental impacts, either directly or indirectly, through their effect on wetlands, endangered species habitat, stream channelization, longer culvert lengths, larger right-of-way purchases, and construction access limits. Due to an ever more restrictive regulatory environment, these impacts are problematic. The result is the loss of valuable natural resources to the public, lengthy permitting review processes for the department of transportation and permitting agencies, and the additional expenditures of time and money for all parties. The purpose of this project was to review existing stability berm alternatives for potential use in environmentally sensitive areas. The project also evaluates how stabilization technologies are made feasible, desirable, and cost-effective for transportation projects and determines which alternatives afford practical solutions for avoiding and minimizing impacts to environmentally sensitive areas. An online survey of engineers at state departments of transportation was also conducted to assess the frequency and cost effectiveness of the various stabilization technologies. Geotechnical engineers that responded to the survey overwhelmingly use geosynthetic reinforcement as a suitable and cost-effective solution for stabilizing embankments and cut slopes. Alternatively, chemical stabilization and installation of lime/cement columns is rarely a remediation measure employed by state departments of transportation.

Do pornographic websites influence boys' sexual behavior?

Purpose: To compare the sexual behavior of adolescent males who do and do not watch pornographic websites. Methods: This study was presented as a school survey. Data were drawn from the 2002 Swiss Multicenter Adolescent Survey on Health (SMASH02) database, a survey including 7,548 adolescents age 16-20. The setting was post-mandatory schools in Switzerland. A total of 2,891 male students who connected to the internet in the last 30 days were enrolled and distributed into two groups: boys who deliberately watched pornographic websites in the last 30 days (n ¼ 942; 33%) and boys who did not (n ¼ 1,949; 67%). Socio-demographic characteristics; frequency of connection to the internet; sexual behavior parameters (having a girlfriend and if yes, for more or less than 6 months; having had sexual intercourse; age at first sexual intercourse; use of a condom at last sexual intercourse; number of sexual partners; having made a partner pregnant). Results: A logistic regression was performed using STATA 9.2. The only significant socio-demographic variable was having a low socioeconomic status (adjusted odds ratio [AOR] 1.66); no difference was found for age and academic track between the two groups. Boys who watch pornographic websites were also significantly more likely to connect frequently to the internet (one day a week: AOR 1.75; several days a week: AOR 2.36; every day: AOR 3.11), to have had sexual intercourse (AOR 2.06), and to have had their first sexual intercourse before age 15 (AOR 1.48). The stability of the relationship with their girlfriend did not appear to have any influence on the search for pornography on the internet. Conclusions: About one third of boys in our sample report having accessed pornographic websites in the last 30 days, a proportion similar to other studies. Watching such websites increases with the frequency of connection to the internet and seems to be correlated with an earlier sexual activity debut among adolescent males. However, having had first sexual intercourse before age 15 is the only sexual risk behavior that seems to be increased when watching pornographic websites among boys. Further studies should address the causality of this correlation and the factors influencing the search for pornography on the web among boys, in order to explore some new ways of prevention about sexual risk behaviors. Sources of Support: The SMASH02 survey was carried out with the financial support of the Swiss Federal Office of Public Health and the participating cantons.