Synchronous versus asynchronous modeling of gene regulatory networks.

MOTIVATION: In silico modeling of gene regulatory networks has gained some momentum recently due to increased interest in analyzing the dynamics of biological systems. This has been further facilitated by the increasing availability of experimental data on gene-gene, protein-protein and gene-protein interactions. The two dynamical properties that are often experimentally testable are perturbations and stable steady states. Although a lot of work has been done on the identification of steady states, not much work has been reported on in silico modeling of cellular differentiation processes. RESULTS: In this manuscript, we provide algorithms based on reduced ordered binary decision diagrams (ROBDDs) for Boolean modeling of gene regulatory networks. Algorithms for synchronous and asynchronous transition models have been proposed and their corresponding computational properties have been analyzed. These algorithms allow users to compute cyclic attractors of large networks that are currently not feasible using existing software. Hereby we provide a framework to analyze the effect of multiple gene perturbation protocols, and their effect on cell differentiation processes. These algorithms were validated on the T-helper model showing the correct steady state identification and Th1-Th2 cellular differentiation process. AVAILABILITY: The software binaries for Windows and Linux platforms can be downloaded from http://si2.epfl.ch/~garg/genysis.html.

The Early Toarcian anoxia, a synchronous event in the Western Tethys? An approach by quantitative biochronology (Unitary Associations), applied on calcareous nannofossils

During the Early Toarcian, major paleoenvironnemental and paleoceanographical changes occurred, leading to an oceanic anoxic event (OAE) and to a perturbation of the carbon isotope cycle. Although the standard biochronology of the Lower Jurassic is essentially based upon ammonites, in recent years biostratigraphy based on calcareous nannofossils and dinoflagellate cysts is increasingly used to date Jurassic rocks. However, the precise dating and correlation of the Early Toarcian OAE, and of the associated delta C-13 anomaly in different settings of the western Tethys, are still partly problematic, and it is still unclear whether these events are synchronous or not. In order to allow more accurate correlations of the organic rich levels recorded in the Lower Toarcian OAE, this account proposes a new biozonation based on a quantitative biochronology approach, the Unitary Associations (UA), applied to calcareous nannofossils. This study represents the first attempt to apply the UA method to Jurassic nannofossils. The study incorporates eighteen sections distributed across western Tethys and ranging from the Pliensbachian to Aalenian, comprising 1220 samples and 72 calcareous nannofossil taxa. The BioGraph [Savary, J., Guex, J., 1999. Discrete biochronological scales and unitary associations: description of the Biograph Computer program. Memoires de Geologie de Lausanne 34, 282 pp] and UA-Graph (Copyright Hammer O., Guex and Savary, 2002) softwares provide a discrete biochronological framework based upon multi-taxa concurrent range zones in the different sections. The optimized dataset generates nine UAs using the co-occurrences of 56 taxa. These UAs are grouped into six Unitary Association Zones (UA-Z), which constitute a robust biostratigraphic synthesis of all the observed or deduced biostratigraphic relationships between the analysed taxa. The UA zonation proposed here is compared to ``classic'' calcareous nannofossil biozonations, which are commonly used for the southern and the northern sides of Tethys. The biostratigraphic resolution of the UA-Zones varies from one nannofossil subzone or part of it to several subzones, and can be related to the pattern of calcareous nannoplankton originations and extinctions during the studied time interval. The Late Pliensbachian - Early Toarcian interval (corresponding to the UA-Z II) represents a major step in the Jurassic nannoplankton radiation. The recognized UA-Zones are also compared to the carbon isotopic negative excursion and TOC maximum in five sections of central Italy, Germany and England, with the aim of providing a more reliable correlation tool for the Early Toarcian OAE, and of the associated isotopic anomaly, between the southern and northern part of western Tethys. The results of this work show that the TOC maximum and delta C-13 negative excursion correspond to the upper part of the UA-Z II (i.e., UA 3) in the sections analysed. This suggests that the Early Toarcian OAE was a synchronous event within the western Tethys. (c) 2006 Elsevier B.V. All rights reserved.

Synchronous and baroceptor-sensitive oscillations in skin microcirculation: evidence for central autonomic control.

To determine whether skin blood flow is local or takes part in general regulatory mechanisms, we recorded laser-Doppler flowmetry (LDF; left and right index fingers), blood pressure, muscle sympathetic nerve activity (MSNA), R-R interval, and respiration in 10 healthy volunteers and 3 subjects after sympathectomy. We evaluated 1) the synchronism of LDF fluctuations in two index fingers, 2) the relationship with autonomically mediated fluctuations in other signals, and 3) the LDF ability to respond to arterial baroreflex stimulation (by neck suction at frequencies from 0.02 to 0.20 Hz), using spectral analysis (autoregressive uni- and bivariate, time-variant algorithms). Synchronous LDF fluctuations were observed in the index fingers of healthy subjects but not in sympathectomized patients. LDF fluctuations were coherent with those obtained for blood pressure, MSNA, and R-R interval. LDF fluctuations were leading blood pressure in the low-frequency (LF; 0.1 Hz) band and lagging in the respiratory, high-frequency (HF; approximately 0.25 Hz) band, suggesting passive "downstream" transmission only for HF and "upstream" transmission for LF from the microvessels. LDF fluctuations were responsive to sinusoidal neck suction up to 0.1 Hz, indicating response to sympathetic modulation. Skin blood flow thus reflects modifications determined by autonomic activity, detectable by frequency analysis of spontaneous fluctuations.

The "reversed" treatment approach for synchronous liver metastases and colorectal cancer. Feasability and first results

Background: The « reversed treatment» approach inverts the treatment sequence of¦advanced synchronous colorectal liver metastases - i.e. the liver metastasis is¦treated first, followed by resection of the primary tumor. Chemotherapy is performed¦before and after liver surgery. We recently started to use a reversed treatment¦approach in selected patients. The aim of this study is to critically assess this new¦treatment modality.¦Methods: Nine patients (7 male, 2 female, mean age 62 years) benefited from this¦new treatment between November 2008 and May 2010. The data were collected¦retrospectively.¦Results: All patients responded to the neoadjuvant chemotherapy. The median¦number of liver metastases was 6 (range 1 - 22). The median size of the largest liver¦metastases was 4.3 cm (range 2.6 - 13 cm). Three patients had portal vein¦embolization prior to liver surgery. Two patients could not complete the treatment.¦One had to undergo emergency surgery for occluding colonic tumor. The second one¦showed liver recurrence before starting the adjuvant chemotherapy. The seven¦patients who completed the treatment are still alive after a median time of 27 months¦(range 17 - 37 months). Seven of them had recurrence (1 rectal, 6 liver). The median¦disease-free survival was 9 months (range 0 - 17 months).¦Conclusion: Based on our preliminary experiences, the reversed strategy shows¦encouraging results for the treatment of advanced synchronous colorectal liver¦metastases in well selected patients. The treatment was generally well tolerated and¦long term survival seems to be prolonged.

Simulations in Evolution. III. Randomness as a Generator of Opportunities.

In Neo-Darwinism, variation and natural selection are the two evolutionary mechanisms which propel biological evolution. Our previous reports presented a histogram model to simulate the evolution of populations of individuals classified into bins according to an unspecified, quantifiable phenotypic character, and whose number in each bin changed generation after generation under the influence of fitness, while the total population was maintained constant. The histogram model also allowed Shannon entropy (SE) to be monitored continuously as the information content of the total population decreased or increased. Here, a simple Perl (Practical Extraction and Reporting Language) application was developed to carry out these computations, with the critical feature of an added random factor in the percent of individuals whose offspring moved to a vicinal bin. The results of the simulations demonstrate that the random factor mimicking variation increased considerably the range of values covered by Shannon entropy, especially when the percentage of changed offspring was high. This increase in information content is interpreted as facilitated adaptability of the population.

Whole-body vibration training: metabolic cost of synchronous, side-alternating or no vibrations.

Whole-body vibration training improves strength and can increase maximal oxygen consumption ([·V]O(2max)). No study has compared the metabolic demand of synchronous and side-alternating whole-body vibration. We measured [·V]O₂ and heart rate during a typical synchronous or side-alternating whole-body vibration session in 10 young female sedentary participants. The 20-min session consisted of three sets of six 45-s exercises, with 15 s recovery between exercises. Three conditions were randomly tested on separate days: synchronous at 35 Hz and 4 mm amplitude, side-alternating at 26 Hz and 7.5 mm amplitude (peak acceleration matched at 20 g in both vibration conditions), and no vibrations. Mean [·V]O₂ (expressed as %[·V]O(2max)) did not differ between conditions: 29.7 ± 4.2%, 32.4 ± 6.5%, and 28.7 ± 6.7% for synchronous, side-alternating, and no vibrations respectively (P = 0.103). Mean heart rate (% maximal heart rate) was 65.6 ± 7.3%, 69.8 ± 7.9%, and 64.7 ± 5.6% for synchronous, side-alternating, and no vibrations respectively, with the side-alternating vibrations being significantly higher (P = 0.019). When analysing changes over exercise sessions, mean [·V]O₂ was higher for side-alternating (P < 0.001) than for synchronous and no vibrations. In conclusion, side-alternating whole-body vibration elicits higher heart rate responses than synchronous or no vibrations, and could elevate [·V]O₂, provided the session lasts more than 20 min.

A synchronous, middle early Aptian age for the demise of the Helvetic Urgonian platform related to the unfolding oceanic anoxic event 1a ("Selli event")

Un âge synchrone (partie moyenne de l'Aptien inférieur) de l'ennoiement de la plate-forme Urgonienne helvétique en relation avec l'événement océanique anoxique 1a ("événement Selli"). - La fin de la plate-forme urgonienne, calibrée par analyse des isotopes stables du carbone sur roche totale et par biostratigraphie basée sur les ammonites, est datée du milieu de l'Aptien inférieur (Près de la limite des zones weissi et deshayesi). Cet arrêt, synchrone dans des coupes représentatives du domaine helvétique alpin, est un événement environemental majeur renregistré en France, en Espagne, au Protugal, en Oman, au Mexique et dans le domaine Pacifique. En tenant compte des limites de résolution de la biostatrigraphie et des autres techniques de datation, cet épisode semble également être synchrone à l'échelle globale. Pour beaucoup d'auteurs, la disparition de récifs de coraux et de rudistes corrélée à la fin de la sédimentation urgonienne correspond à la mise en place de conditions anoxiques à l'Aptien inférieur. Celles-ci caractérisent un événement d'importance global: l'événement anoxique OAE 1a.

Influence of the implanted pulse generator as reference electrode in finite element model of monopolar deep brain stimulation.

Electrical deep brain stimulation (DBS) is an efficient method to treat movement disorders. Many models of DBS, based mostly on finite elements, have recently been proposed to better understand the interaction between the electrical stimulation and the brain tissues. In monopolar DBS, clinically widely used, the implanted pulse generator (IPG) is used as reference electrode (RE). In this paper, the influence of the RE model of monopolar DBS is investigated. For that purpose, a finite element model of the full electric loop including the head, the neck and the superior chest is used. Head, neck and superior chest are made of simple structures such as parallelepipeds and cylinders. The tissues surrounding the electrode are accurately modelled from data provided by the diffusion tensor magnetic resonance imaging (DT-MRI). Three different configurations of RE are compared with a commonly used model of reduced size. The electrical impedance seen by the DBS system and the potential distribution are computed for each model. Moreover, axons are modelled to compute the area of tissue activated by stimulation. Results show that these indicators are influenced by the surface and position of the RE. The use of a RE model corresponding to the implanted device rather than the usually simplified model leads to an increase of the system impedance (+48%) and a reduction of the area of activated tissue (-15%).

Tracking the source of cerebellar epilepsy: hemifacial seizures associated with cerebellar cortical dysplasia.

Traditionally, subcortical structures such as the cerebellum are supposed to exert a modulatory effect on epileptic seizures, rather than being the primary seizure generator. We report a 14-month old girl presenting, since birth, with seizures symptomatic of a right cerebellar dysplasia, manifested as paroxystic contralateral hemifacial spasm and ipsilateral facial weakness. Multimodal imaging was used to investigate both anatomical landmarks related to the cerebellar lesion and mechanisms underlying seizure generation. Electric source imaging (ESI) supported the hypothesis of a right cerebellar epileptogenic generator in concordance with nuclear imaging findings; subsequently validated by intra-operative intralesional recordings. Diffusion spectrum imaging-related tractography (DSI) showed severe cerebellar structural abnormalities confirmed by histological examination. We suggest that hemispheric cerebellar lesions in cases like this are likely to cause epilepsy via an effect on the facial nuclei through ipsilateral and contralateral aberrant connections.

Distinct neurophysiological mechanisms mediate mixing costs and switch costs.

Using event-related potentials (ERPs), we investigated the neural response associated with preparing to switch from one task to another. We used a cued task-switching paradigm in which the interval between the cue and the imperative stimulus was varied. The difference between response time (RT) to trials on which the task switched and trials on which the task repeated (switch cost) decreased as the interval between cue and target (CTI) was increased, demonstrating that subjects used the CTI to prepare for the forthcoming task. However, the RT on repeated-task trials in blocks during which the task could switch (mixed-task blocks) were never as short as RTs during single-task blocks (mixing cost). This replicates previous research. The ERPs in response to the cue were compared across three conditions: single-task trials, switch trials, and repeat trials. ERP topographic differences were found between single-task trials and mixed-task (switch and repeat) trials at approximately 160 and approximately 310 msec after the cue, indicative of changes in the underlying neural generator configuration as a basis for the mixing cost. In contrast, there were no topographic differences evident between switch and repeat trials during the CTI. Rather, the response of statistically indistinguishable generator configurations was stronger at approximately 310 msec on switch than on repeat trials. By separating differences in ERP topography from differences in response strength, these results suggest that a reappraisal of previous research is appropriate.

Age-specific incidence of all neoplasms after colorectal cancer.

PURPOSE: Patients diagnosed with a specific neoplasm tend to have a subsequent excess risk of the same neoplasm. The age incidence of a second neoplasm at the same site is approximately constant with age, and consequently the relative risk is greater at younger age. It is unclear whether such a line of reasoning can be extended from a specific neoplasm to the incidence of all neoplasms in subjects diagnosed with a defined neoplasm. METHODS: We considered the age-specific incidence of all non-hormone-related epithelial neoplasms after a first primary colorectal cancer (n = 9542) in the Vaud Cancer Registry data set. RESULTS: In subjects with a previous colorectal cancer, the incidence rate of all other epithelial non-hormone-related cancers was stable around 800 per 100,000 between age 30 and 60 years, and rose only about twofold to reach 1685 at age 70 to 79 years and 1826 per 100,000 at age 80 years or older. After excluding synchronous cancers, the rise was only about 1.5-fold, that is, from about 700 to 1000. In the general population, the incidence rate of all epithelial non-hormone-related cancers was 29 per 100,000 at age 30 to 39 years, and rose 30-fold to 883 per 100,000 at age 70 to 79 years. Excluding colorectal cancers, the rise of all non-hormone-related cancers was from 360 per 100,000 at age 40 to 49 years to 940 at age 70 to 79 years after colorectal cancer, and from 90 to 636 per 100,000 in the general population (i.e., 2.6- vs. 7.1-fold). CONCLUSIONS: The rise of incidence with age of all epithelial non-hormone-related second cancers after colorectal cancer is much smaller than in the general population. This can possibly be related to the occurrence of a single mutational event in a population of susceptible individuals, although alternative models are plausible within the complexity of the process of carcinogenesis.

Food availability dictates the timing of parturition in insectivorous mouse-eared bats

Which of these two confounding factors, weather or food availability - that largely correlate and interact - controls the timing of parturition in insectivorous bats? To answer this question. we took advantage of a predator-prey system that offers a unique opportunity to perform natural experiments. The phenology of reproduction of two sibling bat species that inhabit the same colonial roosts, but exploit different feeding niches. was investigated. Myotis myotis feeds mainly on carabid beetles, a food source available from the end of hibernation onwards, whereas bush crickets, the main prey of M. blythii, are not available early in the season due to their successive instars; cockchafers are actually the sole possible alternative prey for M. blythii at that time of the year, but they occur every third year only, independently of local weather conditions. By comparing the species responses to the presence/absence of cockchafers, we could test the hypothesis that food availability, rather than climate. influences the timing of bat parturition. Our data show that Nt. blythii gave birth, on average. 10 d later than M. myotis in years without cockchafers, whilst parturition (1) was synchronous during cockchafer years, and (2) did not show much among-year time variation in M. myotis. This suggests that food availability is the chief factor regulating the timing of parturition in mouse-eared bats.

Systematic review and meta-analysis of risk factors for survival after lung metastasectomy in colorectal cancer patients

Objective: Resection of lung metastases (LM) from colorectal cancer (CRC)¦is increasingly performed with a curative intent.Most series report small groups¦of patients, and it is currently not possible to identify those CRC patients who¦may benefit the most of surgical management. It is clinically relevant to assess¦risk factors for prolonged survival after this type of procedures.¦Methods: A meta analysis of 24 series published between 2000 and 2011¦which focused on surgical management of LM from CRC and included more¦than 40 patients each, with or without prior resection of in transit liver¦metastases. Random effects were calculated for five variables considered as¦potential prognostic factors.¦Results: A total of 2815 patients who underwent surgery with a curative¦intent were considered in this analysis. Four parameters were associated with¦a decreased survival: 1) a short disease-free interval between primary tumor¦resection and development of LM (HR = 1·59, 95% CI 1·27-1·98); 2) multiple¦LM (HR = 2·04, 95%CI 1·72-2·41); 3) positive hilar/mediastinal lymph nodes¦(HR = 1·65, 95% CI 1·35-2·02); and 4) a high prethoracotomy CEA value (HR¦=1·91, 95% CI 1·57-2·32). By comparison, a history of resected liver metastases¦(HR = 1·36, 95% CI 0·92-2·03) did not achieve statistical significance.¦Conclusion: Risk factors for poor clinical outcome after surgery for lung¦metastases in CRC patients include: 1) synchronous lung metastases; 2) high¦pre-thoracotomy CEA; 3) hilar nodes involvement; and 4) multiple pulmonary¦lesions.

Garnet growth in the metasediments of the Zermatt-Saas Fee zone (western alps)

The Zermatt-Saas Fee Zone (ZSZ) in the Western Alps consists of multiple slices of ultramafic, mafic and metasedimentary rocks. They represent the remnants of the Mesozoic Piemonte-Ligurian oceanic basin which was subducted to eclogite facies conditions with peak pressures and temperatures of up to 20-28 kbar and 550-630 °C, followed by a greenschist overprint during exhumation. Previous studies, emphasizing on isotopie geochronology and modeling of REE-behavior in garnets from mafic eclogites, suggest that the ZSZ is buildup of tectonic slices which underwent a protracted diachronous subduction followed by a rapid synchronous exhumation. In this study Rb/Sr geochronology is applied to phengite included in garnets from metasediments of two different slices of the ZSZ to date garnet growth. Inclusion ages for 2 metapelitic samples from the same locality from the first slice are 44.25 ± 0.48 Ma and 43.19 ± 0.32 Ma. Those are about 4 Ma older than the corresponding matrix mica ages of respectively 40.02 ± 0.13 Ma and 39.55 ± 0.25 Ma. The inclusion age for a third calcschist sample, collected from a second slice, is 40.58 ± 0.24 Ma and the matrix age is 39.8 ± 1.5 Ma. The results show that garnet effectively functioned as a shield, preventing a reset of the Rb/Sr isotopie clock in the included phengites to temperatures well above the closure of Sr in mica. The results are consistent with the results of former studies on the ZSZ using both Lu/Hf and Sm/Nd geochronology on mafic eclogites. They confirm that at least parts of the ZSZ underwent close to peak metamorphic HP conditions younger than 43 m.y. ago before being rapidly exhumed about 40 m.y. ago. Fluid infiltration in rocks of the second slice occurred likely close to the peak metamorphic conditions, resulting in rapid growth of garnets. Similar calcschists from the same slice contain two distinct types of porphyroblast garnets with indications of multiple growth pulses and resorption indicated by truncated chemical zoning patterns. In-situ oxygen isotope Sensitive High Resolution Ion Microprobe (SHRIMP) analyses along profiles on central sections of the garnets reveal variations of up to 5 %o in individual garnets. The complex compositional zoning and graphite inclusion patterns as well as the variations in oxygen isotopes correspond to growing under changing fluid composition conditions caused by external infiltrated fluids. The ultramafic and mafic rocks, which were subducted along with the sediments and form the volumetrically most important part of the ZSZ, are the likely source of those mainly aqueous fluids. - La Zone de Zermatt-Saas Fee (ZZS) est constituée de multiples écailles de roches ultramafiques, mafiques et méta-sédimentaires. Cette zone, qui affleure dans les Alpes occidentales, représente les restes du basin océanique Piémontais-Ligurien d'âge mésozoïque. Lors de la subduction de ce basin océanique à l'Eocène, les différentes roches composant le planché océanique ont atteint les conditions du faciès éclogitique avec des pressions et des températures maximales estimées entre 20 - 28 kbar et 550 - 630 °C respectivement, avant de subir une rétrogression au faciès schiste vert pendant l'exhumation. Différentes études antérieures combinant la géochronologie isotopique et la modélisation des mécanismes gouvernant l'incorporation des terres rares dans les grenats des éclogites mafiques, suggèrent que la ZZS ne correspond pas à une seule unité, mais est constituée de différentes écailles tectoniques qui ont subi une subduction prolongée et diachrone suivie d'une exhumation rapide et synchrone. Afin de tester cette hypothèse, j'ai daté, dans cette étude, des phengites incluses dans les grenats des méta-sédiments de deux différentes écailles tectoniques de la ZZS, afin de dater la croissance relative de ces grenats. Pour cela j'ai utilisé la méthode géochronologique basée sur la décroissance du Rb87 en Sr87. J'ai daté trois échantillons de deux différentes écailles. Les premiers deux échantillons proviennent de Triftji, au nord du Breithorn, d'une première écaille dont les méta-sédiments sont caractérisés par des bandes méta-pélitiques à grenat et des calcschistes. Le troisième échantillon a été collectionné au Riffelberg, dans une écaille dont les méta-sédiments sont essentiellement des calcschistes qui sont mélangés avec des roches mafiques et des serpentinites. Ce mélange se trouve au-dessus de la grande masse de serpentinites qui forment le Riffelhorn, le Trockenersteg et le Breithorn, et qui est connu sous le nom de la Zone de mélange de Riffelberg (Bearth, 1953). Les inclusions dans les grenats de deux échantillons méta-pélitiques de la première écaille sont datées à 44.25 ± 0.48 Ma et à 43.19 ± 0.32 Ma. Ces âges sont à peu près 4 Ma plus vieux que les âges obtenus sur les phengites provenant de la matrice de ces mêmes échantillons qui donnent des âges de 40.02 ± 0.13 Ma et 39.55 ± 0.25 Ma respectivement. Les inclusions de phengite dans les grenats appartenant à un calcschiste de la deuxième écaille ont un âge de 40.58 ± 0.24 Ma alors que les phengites de la matrice ont un âge de 39.8 ± 1.5 Ma. Pour expliquer ces différences d'âge entre les phengites incluses dans le grenat et les phengites provenant de la matrice, nous suggérons que la cristallisation de grenat ait permis d'isoler ces phengites et de les préserver de tous rééquilibrage lors de la suite du chemin métamorphique prograde, puis rétrograde. Ceci est particulièrement important pour expliquer l'absence de rééquilibrage des phengites dans des conditions de températures supérieures à la température de fermeture du système Rb/Sr pour les phengites. Les phengites en inclusions n'ayant pas pu être datées individuellement, nous interprétons l'âge de 44 Ma pour les inclusions de phengite comme un âge moyen pour l'incorporation de ces phengites dans le grenat. Ces résultats sont cohérents avec les résultats des études antérieures de la ZZS utilisant les systèmes isotopiques de Sm/Nd et Lu/Hf sur des eclogites mafiques. ils confirment qu'aux moins une partie de la ZZS a subi des conditions de pression et de température maximale il y a moins de 44 à 42 Ma avant d'être rapidement exhumée à des conditions métamorphiques du faciès schiste vert supérieur autour de 40 Ma. Cette étude détaillée des grenats a permis, également, de mettre en évidence le rôle des fluides durant le métamorphisme prograde. En effet, si tous les grenats montrent des puises de croissance et de résorption, on peut distinguer, dans différents calcschists provenant de la deuxième écaille, deux types distincts de porphyroblast de grenat en fonction de la présence ou non d'inclusions de graphite. Nous lions ces puises de croissances/résorptions ainsi que la présence ou l'absence de graphite en inclusion dans les grenats à l'infiltration de fluides dans le système, et ceci durant tous le chemin prograde mais plus particulièrement proche et éventuellement peu après du pic du métamorphisme comme le suggère l'âge de 40 Ma mesuré dans les inclusions de phengites de l'échantillon du Riffelberg. Des analyses in-situ d'isotopes d'oxygène réalisé à l'aide de la SHRIMP (Sensitive High Resolution Ion Microprobe) dans des coupes centrales des grenats indiquent des variations jusqu'à 5 %o au sein même d'un grenat. Les motifs de zonations chimiques et d'inclusions de graphite complexes, ainsi que les variations du δ180 correspondent à une croissance de grenat sous des conditions de fluides changeantes dues aux infiltrations de fluides externes. Nous lions l'origine de ces fluides aqueux aux unités ultramafiques et mafiques qui ont été subductés avec les méta-sédiments ; unités ultramafiques et mafiques qui forment la partie volumétrique la plus importante de la ZZS.

Assessment of left coronary artery beat-to-beat repositioning in order to propose an optimal acquisition window for coronary MRA

Introduction: Coronary magnetic resonance angiography (MRA) is a medical imaging technique that involves collecting data from consecutive heartbeats, always at the same time in the cardiac cycle, in order to minimize heart motion artifacts. This technique relies on the assumption that coronary arteries always follow the same trajectory from heartbeat to heartbeat. Until now, choosing the acquisition window in the cardiac cycle was based exclusively on the position of minimal coronary motion. The goal of this study was to test the hypothesis that there are time intervals during the cardiac cycle when coronary beat-to-beat repositioning is optimal. The repositioning uncertainty values in these time intervals were then compared with the intervals of low coronary motion in order to propose an optimal acquisition window for coronary MRA. Methods: Cine breath-hold x-ray angiograms with synchronous ECG were collected from 11 patients who underwent elective routine diagnostic coronarography. Twenty-three bifurcations of the left coronary artery were selected as markers to evaluate repositioning uncertainty and velocity during cardiac cycle. Each bifurcation was tracked by two observers, with the help of a user-assisted algorithm implemented in Matlab (The Mathworks, Natick, MA, USA) that compared the trajectories of the markers coming from consecutive heartbeats and computed the coronary repositioning uncertainty with steps of 50ms until 650ms after the R-wave. Repositioning uncertainty was defined as the diameter of the smallest circle encompassing the points to be compared at the same time after the R-wave. Student's t-tests with a false discovery rate (FDR, q=0.1) correction for multiple comparison were applied to see whether coronary repositioning and velocity vary statistically during cardiac cycle. Bland-Altman plots and linear regression were used to assess intra- and inter-observer agreement. Results: The analysis of left coronary artery beat-to-beat repositioning uncertainty shows a tendency to have better repositioning in mid systole (less than 0.84±0.58mm) and mid diastole (less than 0.89±0.6mm) than in the rest of the cardiac cycle (highest value at 50ms=1.35±0.64mm). According to Student's t-tests with FDR correction for multiple comparison (q=0.1), two intervals, in mid systole (150-200ms) and mid diastole (550-600ms), provide statistically better repositioning in comparison with the early systole and the early diastole. Coronary velocity analysis reveals that left coronary artery moves more slowly in end systole (14.35±11.35mm/s at 225ms) and mid diastole (11.78±11.62mm/s at 625ms) than in the rest of the cardiac cycle (highest value at 25ms: 55.96±22.34mm/s). This was confirmed by Student's t-tests with FDR correction for multiple comparison (q=0.1, FDR-corrected p-value=0.054): coronary velocity values at 225, 575 and 625ms are not much different between them but they are statistically inferior to all others. Bland-Altman plots and linear regression show that intra-observer agreement (y=0.97x+0.02 with R²=0.93 at 150ms) is better than inter-observer (y=0.8x+0.11 with R²=0.67 at 150ms). Discussion: The present study has demonstrated that there are two time intervals in the cardiac cycle, one in mid systole and one in mid diastole, where left coronary artery repositioning uncertainty reaches points of local minima. It has also been calculated that the velocity is the lowest in end systole and mid diastole. Since systole is less influenced by heart rate variability than diastole, it was finally proposed to test an acquisition window between 150 and 200ms after the R-wave.