High constant incidence of second primary colorectal cancer.

Patients who had a colorectal cancer have a 1.5- to 2-fold excess risk of a second colorectal cancer as compared to the general population, the excess being higher at younger age at diagnosis. To further investigate the risk and the age-relation of the incidence of second primary colorectal cancer, we considered 9,389 first colon and rectal cancers registered in the Vaud Cancer Registry, Switzerland, between 1974 and 2008, and followed-up to the end of 2008 for a total of 44,113 person-years. There were 136 second colorectal cancers versus 90.5 expected, corresponding to a standardized incidence ratio (SIR) of 1.5 (95% confidence interval, CI, 1.3-1.8). The SIRs were not heterogeneous between men and women, and in strata of calendar year at diagnosis, duration of follow-up, and subsite. However, the SIR was 7.5 (95% CI 4.2-12.4) for subjects diagnosed below age 50 and declined thereafter to reach 1.0 (95% CI 0.6-1.6) at age 80 or over. Consequently, the incidence of second primary colorectal cancer was stable, and exceedingly high, around 300-400/100,000 between age 30-39 and 70 or over. This age pattern is consistent with the existence of a single mutational event in a population of highly susceptible individuals.

Quantitative comparison between simulations of seismic wave propagation in heterogeneous poro-elastic media and equivalent visco-elastic solids for marine-type environments

There is increasing evidence to suggest that the presence of mesoscopic heterogeneities constitutes an important seismic attenuation mechanism in porous rocks. As a consequence, centimetre-scale perturbations of the rock physical properties should be taken into account for seismic modelling whenever detailed and accurate responses of specific target structures are desired, which is, however, computationally prohibitive. A convenient way to circumvent this problem is to use an upscaling procedure to replace each of the heterogeneous porous media composing the geological model by corresponding equivalent visco-elastic solids and to solve the visco-elastic equations of motion for the inferred equivalent model. While the overall qualitative validity of this procedure is well established, there are as of yet no quantitative analyses regarding the equivalence of the seismograms resulting from the original poro-elastic and the corresponding upscaled visco-elastic models. To address this issue, we compare poro-elastic and visco-elastic solutions for a range of marine-type models of increasing complexity. We found that despite the identical dispersion and attenuation behaviour of the heterogeneous poro-elastic and the equivalent visco-elastic media, the seismograms may differ substantially due to diverging boundary conditions, where there exist additional options for the poro-elastic case. In particular, we observe that at the fluid/porous-solid interface, the poro- and visco-elastic seismograms agree for closed-pore boundary conditions, but differ significantly for open-pore boundary conditions. This is an important result which has potentially far-reaching implications for wave-equation-based algorithms in exploration geophysics involving fluid/porous-solid interfaces, such as, for example, wavefield decomposition.

MRS glucose mapping and PET joining forces: re-evaluation of the lumped constant in the rat brain under isoflurane anaesthesia.

Although numerous positron emission tomography (PET) studies with (18) F-fluoro-deoxyglucose (FDG) have reported quantitative results on cerebral glucose kinetics and consumption, there is a large variation between the absolute values found in the literature. One of the underlying causes is the inconsistent use of the lumped constants (LCs), the derivation of which is often based on multiple assumptions that render absolute numbers imprecise and errors hard to quantify. We combined a kinetic FDG-PET study with magnetic resonance spectroscopic imaging (MRSI) of glucose dynamics in Sprague-Dawley rats to obtain a more comprehensive view of brain glucose kinetics and determine a reliable value for the LC under isoflurane anaesthesia. Maps of Tmax /CMRglc derived from MRSI data and Tmax determined from PET kinetic modelling allowed to obtain an LC-independent CMRglc . The LC was estimated to range from 0.33 ± 0.07 in retrosplenial cortex to 0.44 ± 0.05 in hippocampus, yielding CMRglc between 62 ± 14 and 54 ± 11 μmol/min/100 g, respectively. These newly determined LCs for four distinct areas in the rat brain under isoflurane anaesthesia provide means of comparing the growing amount of FDG-PET data available from translational studies.

Estimating parameters for generalized mass action models using constraint propagation.

As modern molecular biology moves towards the analysis of biological systems as opposed to their individual components, the need for appropriate mathematical and computational techniques for understanding the dynamics and structure of such systems is becoming more pressing. For example, the modeling of biochemical systems using ordinary differential equations (ODEs) based on high-throughput, time-dense profiles is becoming more common-place, which is necessitating the development of improved techniques to estimate model parameters from such data. Due to the high dimensionality of this estimation problem, straight-forward optimization strategies rarely produce correct parameter values, and hence current methods tend to utilize genetic/evolutionary algorithms to perform non-linear parameter fitting. Here, we describe a completely deterministic approach, which is based on interval analysis. This allows us to examine entire sets of parameters, and thus to exhaust the global search within a finite number of steps. In particular, we show how our method may be applied to a generic class of ODEs used for modeling biochemical systems called Generalized Mass Action Models (GMAs). In addition, we show that for GMAs our method is amenable to the technique in interval arithmetic called constraint propagation, which allows great improvement of its efficiency. To illustrate the applicability of our method we apply it to some networks of biochemical reactions appearing in the literature, showing in particular that, in addition to estimating system parameters in the absence of noise, our method may also be used to recover the topology of these networks.

Growth mechanisms and timing of emplacement of a vertically layered mafic intrusion in the feeder zone of an ocean island volcano (Fuerteventura, Canary islands)

Résumé pour le grand public L'île de Fuerteventura (Canaries) offre l'occasion rare d'observer les racines d'un volcan océanique édifié il y a 25 à 30 millions d'années et complètement érodé. On y voit de nombreux petits plutons de forme et composition variées, témoignant d'autant d'épisodes de l'activité magmatique. L'un de ces plutons, appelé PX1, présente une structure inhabituelle formée d'une alternance de bandes verticales d'épaisseur métrique à hectométrique de roches sombres de composition pyroxénilique ou gabbroïque. Les pyroxénites résultent clairement de l'accumulation de cristaux de pyroxènes et non de la simple solidification d'un magma? Se pose dès lors la question de la nature du processus qui a conduit à l'accumulation verticale de niveaux concentrés en pyroxènes. En effet, les litages pyroxénitiques classiques sont subhorizontaux, car ils résultent de l'accumulation gravitaire des cristaux séparés du magma dont ils cristalli¬sent par sédimentation. Cette étude vise à identifier et comprendre les mécanismes qui ont engendré ce Iitage minéralogique vertical et l'im¬portant volume de ces faciès cumulatifs. Nous nous sommes également intéressés aux conditions de pression et de température régnant au moment de la mise en place du pluton, ainsi qu'à sa durée de vie et à sa vitesse de refroidis¬sement. Enfin une approche géochimique nous a permis de préciser la nature de la source mantellique des magmas liés à cette activité magmatique. PX1 est en réalité un complexe filonien formé à des conditions de pression et de température de 1-2 kbar et 1050- 1100°C; sa construction a nécessité au moins 150 km3 de magma. L'alternance d'horizons gabbroïques et pyroxéniti¬ques représente des injections successives de magma sous la forme de filons verticaux, mis en place dans un contexte régional en extension. L'étude des orientations des minéraux dans ces faciès révèle que les horizons gabbroïques enregistrent l'extension régionale, alors que les pyroxénites sont générées par une compaction au sein du pluton. Ceci suggère que le régime des contraintes, qui était extensif lors de l'initiation de la mise en place de PX1, est pério¬diquement devenu compressif au sein même du pluton. Cette compression serait liée à des cycles de mise en place où la vitesse de croissance du pluton dépassait celle de l'extension régionale. La différenciation observée au sein de chaque horizon, depuis des pyroxénites riches en olivine jusqu'à des pyroxé¬nites à plagioclase interstitiel et des gabbros, ainsi que la composition géochimique des minéraux qui les constituent suggèrent que chaque filon vertical s'est mis en place à partir d'un magma de composition identique, puis a évolué indépendamment des autres en fonction du régime thermique et du régime des contraintes local. Lorsque le magma en train de cristalliser s'est trouvé en compression, le liquide résiduel a été séparé des cristaux déjà formés et extrait du système, laissant derrière lui une accumulation de cristaux dont la nature et les proportions dépendaient du stade de cristallisation atteint par le magma au moment de l'extraction. Ainsi, les niveaux de pyroxénites à olivine (premier minéral à cristalliser) ont été formés lorsque le magma correspondant était encore peu cristallisé; à l'inverse, les py¬roxénites riches en plagioclase (minéral plus tardif dans la séquence de cristallisation) et certains gabbros à caractère cumulatif résultent d'une compression tardive dans le processus de cristallisation du filon concerné. Les liquides résiduels extraits des niveaux pyroxénitiques sont rarement observés dans PX1, certaines poches et filonets de com¬position anorthositique pourraient en être les témoins. L'essentiel de ces liquides a probablement gagné des niveaux supérieurs du pluton, voire la surface du volcan. L'origine du régime compressif périodique affectant les filons en voie de cristallisation est attribuée aux injections suivantes de magma au sein du pluton, qui se sont succédées à un rythme plus rapide que la vitesse de consolidation des filons. Des datations U/Pb de haute précision sur des cristaux de zircon et de baddeleyite ainsi que40Ar/39Ar sur des cris¬taux d'amphibole révèlent une initiation de la mise en place de PX1 il y a 22.1 ± 0,7 Ma; celle-ci a duré quelque 0,48 ± 0,22 à 0,52 ± 0,29 Ma. Ce laps de temps est compatible avec celui nécessaire à la cristallisation des filons individuels, qui va de moins d'une année lors de l'initiation du magmatisme à 5 ans lors du maximum d'activité de PX1. La présence de cristaux résorbés enregistrant une cristallisation complexe suggère l'existence d'une chambre mag¬matique convective sous-jacente à PX1 et périodiquement rechargée. Les compositions isotopiques des roches étu¬diées révèlent une source mantellique profonde de type point chaud avec une contribution du manteau lithosphéri- que métasomatisé présent sous les îles Canaries. Résumé L'intrusion mafique Miocène PX1 fait partie du soubassement superficiel (0.15-0.2 GPa, 1100 °Q d'un volcan d'île océanique. La particularité de ce pluton est l'existence d'alternances d'unités de gabbros et de pyroxénites qui met¬tent en évidence un litage magmatique vertical (NNE-SSW). Les horizons gabbroiques et pyroxénitiques sont constitués d'unités de différenciation métriques qui suggèrent tine mise en place par injections périodiques de filons verticaux de magma formant un complexe filonien. Chaque filon vertical a subi une différenciation parallèle à un front de solidification sub-vertical parallèle aux bords du filon. Les pyroxénites résultent du fractionnement et de l'accumulation d'olivine ± clinopyroxene ± plagioclase à partir d'un magma basaltique faiblement alcalin et sont interprétées comme étant des imités de différenciation tronquées dont le liquide interstitiel a été extrait par compaction. L'orientation préférentielle des clinopyroxènes dans ces pyroxe- nites (obtenues par analyse EBSD et micro-tomographique) révèle une composante de cisaillement simple dans la genèse de ces roches, ce qui confirme cette interprétation. La compaction des pyroxénites est probablement causée par a mise en place de filons de magma suivants. Le liquide interstitiel expulsé est probablement par ces derniers. Les clinopyroxènes des gabbros, montrent une composante de cisaillement pure suggérant qu'ils sont affectés par une déformation syn-magmatique parallèle aux zones de cisaillement NNE-SSW observées autour de PX1 et liées au contexte tectonique Miocène d'extension régionale. Ceci suggère que les gabbros sont liés à des taux de mise en place faibles à la fin de cycles d'activité magmatique et sont peu ou pas affectés par la compaction. L'initiation et la géométrie de PX1 sont donc contrôlées par le contexte tectonique régional d'extension alors que les taux et les volumes de magma dépendent de facteurs liés à la source. Des taux d'injection élevés résultent probable¬ment en une croissance du pluton supérieure à la place crée par cette extension. Dans ce cas de figure, la propagation des nouveaux dykes et l'inaptitude du magma à circuler à travers les anciens dykes cristallisés pourrait causer une augmentation de la pression non-lithostatique sur ces derniers, exprimée par un cisaillement simple et l'expulsion du liquide interstitiel qu'ils contiennent (documenté par les zones de collecte anorthositiques). Les compositions en éléments majeurs et traces des gabbros et pyroxenites de PX1 sont globalement homogènes et dépendent de la nature cumulative des échantillons. Cependant, de petites variations des concentrations en éléments traces ainsi que les teneurs en éléments traces des bordures de clinopyroxenes suggèrent que ces derniers ont subi un processus de rééquilibrage et de cristallisation in situ. L'homogénéité des compositions chimiques des échantillons, ainsi que la présence de grains de clinopyroxene résorbés suggère que le complexe filonien PX1 s'est mis en place au dessus d'une chambre magmatique périodiquement rechargée dans laquelle la convection est efficace. Chaque filon est donc issu d'un même magma, mais a subi une différenciation par cristallisation in situ (jusqu'à 70% de fraction¬nement) indépendamment des autres. Dans ces filons cristallisés, les minéraux cumulatifs subissent un rééquilibrage partiel avec les liquide interstitiel avant que ce dernier ne soit expulsé lors de la compaction (mettant ainsi un terme à la différenciation). Ce modèle de mise en place signifie qu'un minimum de 150Km3 de magma est nécessaire à la genèse de PX1, une partie de ce volume ayant été émis par le 'Central Volcanic Complex' de Fuerteventura. Les rapports isotopiques radiogéniques mesurés révèlent la contribution de trois pôles mantelliques dans la genèse du magma formant PX1. Le mélange de ces pôles HIMU, DMM et EM1 refléterai l'interaction du point chaud Cana¬rien avec un manteau lithosphérique hétérogène métasomatisé. Les petites variations de ces rapports et des teneurs en éléments traces au sein des faciès pourrait refléter des taux de fusion partielle variable de la source, résultant en un échantillonnage variable du manteau lithosphérique métasomatisé lors de son interaction avec le point chaud. Des datations U/Pb de haute précision (TIMS) sur des cristaux de zircon et de baddeleyite extraits de gabbros de PX1 révèlent que l'initiation de la cristallisation du magma a eu lieu il y a 22.10±0.07 Ma et que l'activité magmatique a duré un minimum de 0.48 à 0.52 Ma. Des âges 40Ar/39Ar obtenus sur amphibole sont de 21.9 ± 0.6 à 21.8 ± 0.3 Ma, identiques aux âges U/Pb. La combinaison de ces méthodes de datations, suggère que le temps maximum nécessaire à PX1 pour se refroidir en dessous de la température de fermeture de l'amphibole est de 0.8Ma. Ceci signifie que la durée de vie de PX1 est de 520 000 à 800 000 ans. La coexistence de cristaux de baddeleyite et de zircon dans un gabbro est attribuée à son interaction avec un fluide riche en C02 relâché par les carbonatites encaissantes lors du métamorphisme de contact généré par la mise en place de PX1 environ 160 000 ans après le début de sa mise en place. Les durées de vie obtenue sont en accord avec le modèle de mise en place suggérant une durée de cristallisation poux chaque filon allant de 1 an à 5 ans. Abstract The Miocene PX1 gabbro-pyroxenite intrusion (Fuerteventura, Canary Islands), is interpreted as the shallow-level feeder-zone (0.15-0.2 GPa and 1100-1120°C), to an ocean island volcano. The particularity of PX1 is that it displays a NNE-SSW trending vertical magmatic banding expressed by alternating gabbro and pyroxeriite sequences. The gabbro and pyroxenite sequences consist of metre-thick differentiation units, which suggest emplacement by pe¬riodic injection of magma pulses as vertical dykes that amalgamated, similarly to a sub-volcanic sheeted dyke com¬plex. Individual dykes underwent internal differentiation following a solidification front (favoured by a significant lateral/horizontal thermal gradient) parallel to the dyke edges. Pyroxenitic layers result from the fractionation and accumulation of clinopyroxene ± olivine ± plagioclase crystals from a mildly alkaline basaltic liquid and are interpre¬ted as truncated differentiation sequences, from which residual melts were extracted by compaction. Clinopyroxene mineral orientation in pyroxenites (evidenced by EBSD and micro X-ray tomography analysis) display a marked pure shear component, supporting this interpretation. Compaction and squeezing of the crystal mush is ascribed to the incoming and inflating magma pulses. The resulting expelled interstitial liquid was likely collected and erupted along with the magma flowing through the newly injected dykes. Gabbro sequences represent crystallised coalesced magma batches, emplaced at lower rates at the end of eruptive cycles, and underwent minor melt extraction as evi¬denced by clinopyroxene orientations that record a simple shear component suggesting syn-magmatic deformation parallel to observed NNF.-SSW trending shear-zones induced by the regional tensional Miocene stress-field. The initiation and geometry of PX1 is controlled by the regional extensional tectonic regime whereas rates and vo¬lumes of magma depend on source-related factors. High injection rates are likely to induce intrusion growth rates larger than could be accommodated by the regional extension. In this case, dyke tip geometry and the inability of magma to circulate through previously emplaced and crystallised dykes could result in an increase of non-lithostatic pressure on previously emplaced mushy dyke walls; generating strong pure-shear compaction and interstitial melt expulsion within the feeder-zone as recorded by the cumulitic pyroxenite bands and anorthositic collection zones. The whole-rock major and trace-element chemistry of PX1 gabbros and pyroxenites is globally homogeneous and controlled by the cumulate nature of the samples (i.e. on the modal proportions of olivine, pyroxene, plagioclase and oxides). However, small variations of whole-rock trace-element contents as well as trace-element contents of clinopyroxene rims suggest that in-situ re-equilibration and crystallisation has occurred. Additionally, the global homogeneity and presence of complex zoning of rare resorbed clinopyroxene crystals suggest that the PX1 feeder- zone overlies a periodically replenished and efficiently mixed magma chamber. Each individual dyke of magma thus originated from a compositionally constant mildly alkaline magma and differentiated independently from the others reaching up to 70% fractionation. Following dyke arrest these are affected by interaction with the trapped interstitial liquid prior to its compaction-linked expulsion (thus stopping the differentiation process). This emplacement model implies that minimum amount of approximately 150 km3 of magma is needed to generate PX1, part of it having been erupted through the overlying Central Volcanic Complex of Fuerteventura. The radiogenic isotope ratios of PX1 samples reveal the contribution on three end-members during magma genesis. This mixing of the H1MU, EMI and DMM end-members could reflect the interaction of the deep-seated Canarian mantle plume with a heterogeneous metasomatic and sepentininsed lithospheric mantle. Additionally, the observed trace-element and isotopic variations within the same fades groups could reflect varying degrees of partial melting of the source region, thus tapping more or less large areas of the metasomatised lithospheric mantle during interac¬tion with the plume. High precision ID-TIMS U/Pb zircon and baddeleyite ages from the PX1 gabbro samples, indicate initiation of magma crystallisation at 22.10 ± 0.07 Ma. The magmatic activity lasted a minimum of 0.48 to 0.52 Ma. 40Ar/39Ar amphibole ages are of 21.9 ± 0.6 to 21.8 ± 0.3, identical within errors to the U/Pb ages. The combination of the 40Ar/39Ar and U/Pb datasets imply that the maximum amount of time PX1 took to cool below amphibole Tc is 0.8 Ma, suggesting PX1 lifetime of 520 000 to 800 000 years. On top of this, the coexistence of baddeleyite and zircon in a single sample is ascribed to the interaction of PX1 with C02-rich carbonatite-derived fluids released from the host-rock carbonatites during contact metamorphism 160 000 years after PX1 initiation. These ages are in agreement with the emplacement model, implying a crystallisation time of less than 1 to 5 years for individual dykes.