Fast numerical algorithms for the computation of invariant tori in Hamiltonian systems

In this paper, we develop numerical algorithms that use small requirements of storage and operations for the computation of invariant tori in Hamiltonian systems (exact symplectic maps and Hamiltonian vector fields). The algorithms are based on the parameterization method and follow closely the proof of the KAM theorem given in [LGJV05] and [FLS07]. They essentially consist in solving a functional equation satisfied by the invariant tori by using a Newton method. Using some geometric identities, it is possible to perform a Newton step using little storage and few operations. In this paper we focus on the numerical issues of the algorithms (speed, storage and stability) and we refer to the mentioned papers for the rigorous results. We show how to compute efficiently both maximal invariant tori and whiskered tori, together with the associated invariant stable and unstable manifolds of whiskered tori. Moreover, we present fast algorithms for the iteration of the quasi-periodic cocycles and the computation of the invariant bundles, which is a preliminary step for the computation of invariant whiskered tori. Since quasi-periodic cocycles appear in other contexts, this section may be of independent interest. The numerical methods presented here allow to compute in a unified way primary and secondary invariant KAM tori. Secondary tori are invariant tori which can be contracted to a periodic orbit. We present some preliminary results that ensure that the methods are indeed implementable and fast. We postpone to a future paper optimized implementations and results on the breakdown of invariant tori.

The Argentine ant (Linepithema humile) invasion effects on Mediterranean forests’ trophic web: blue tit (Parus caeruleus) reproduction and development

Les invasions biològiques representen una greu amenaça per al funcionament dels ecosistemes i per a la preservació de la biodiversitat.. La formiga argentina (Linepithema humile) està considerada com una de les 100 espècies invasores més nocives. Prospera en extenses àrees de clima mediterrani de regions temperades i subtropicals de tots els continents amb l’excepció de l’Antàrtida. És una formiga dominant i una competidora agressiva que mitjançant múltiples mecanismes, des de predació directe a competència, produeix efectes negatius en una amplia varietat de taxons, principalment formigues i altres artròpodes, però també vertebrats. S’ha investigat, per primera vegada, els efectes de la formiga invasiva sobre les comunitats d’artròpodes de fullatge i com aquestes pertorbacions es transmeten en la xarxa tròfica del bosc esclerofil•le mediterrani. En les suredes estudiades la invasió de formiga argentina és causa directe de la extinció local de la gran majoria de poblacions de formigues natives. En el període mostrejat s’han constatat també impactes negatius en la diversitat i en l’abundància d’artròpodes natius en les capçades dels arbres, particularment d’erugues. Una avaluació preliminar basada únicament amb dades del 2005 indica que, reduint la disponibilitat d’erugues, la formiga argentina empobreix l’hàbitat reproductiu de la mallerenga blava (Parus caeruleus). La mallerenga blava basa la dieta insectívora estricte de la seva pollada fonamentalment en les erugues. No hem detectat impactes en l’èxit reproductiu de les mallerengues blaves en zones envaïdes. Els polls crescuts en àrees envaïdes assoleixen una condició física similar als de les zones no envaïdes, però la reducció en la disponibilitat d’erugues associada a la invasió de formiga argentina es tradueix en un creixement descompassat i en una menor mida estructural del polls volanders. Així, les pertorbacions en la comunitat d’artròpodes associades a la invasió de la formiga argentina promouen efectes bottom-up que acaben perjudicant el desenvolupament dels polls de mallerenga blava.

The effects of feeding Karlodinium veneficum (PLY # 103; Gymnodinium veneficum Ballantine) to the blue mussel Mytilus edulis

The effects of exposure to the type species for Karlodinium veneficum (PLY # 103) on immune function and histopathology in the blue mussel Mytilus edulis were investigated. Mussels from Whitsand Bay, Cornwall (UK) were exposed to K. veneficum (PLY # 103) for 3 and 6 days. Assays for immune function included total and differential cells counts, phagocytosis and release of extra cellular reactive oxygen species. Histology was carried out on digestive gland and mantle tissues. The toxin cell quota for K. veneficum (PLY #103) was measured by liquid chromatography-mass spectrometry detecting two separable toxins KvTx1 (11.6 ± 5.4 ng/ml) and KvTx2 (47.7 ± 4.2 ng/ml). There were significant effects of K. veneficum exposure with increasing phagocytosis and release of reactive oxygen species following 6 days exposure. There were no significant effects on total cell counts. However, differential cell counts did show significant effects after 3 days exposure to the toxic alga. All mussels produced faeces but not pseudofaeces indicating that algae were not rejected prior to ingestion. Digestive glands showed ingestion of the algae and hemocyte infiltration after 3 days of exposure, whereas mantle tissue did not show differences between treatments. As the effects of K. veneficum were not observed in the mantle tissue it can be hypothesized that the algal concentration was not high enough, or exposure long enough, to affect all the tissues. Despite being in culture for more than 50 years the original K. veneficum isolate obtained by Mary Parke still showed toxic effects on mussels.

Pathology and immune response of the blue mussel (Mytilus edulis L.) after an exposure to the harmful dinoflagellate Prorocentrum minimum

The harmful dinoflagellate Prorocentrum minimum has different effects upon various species of grazing bivalves, and these effects also vary with life-history stage. Possible effects of this dinoflagellate upon mussels have not been reported; therefore, experiments exposing adult blue mussels, Mytilus edulis, to P. minimum were conducted. Mussels were exposed to cultures of toxic P. minimum or benign Rhodomonas sp. in glass aquaria. After a short period of acclimation, samples were collected on day 0 (before the exposure) and after 3, 6, and 9 days of continuous-exposure experiment. Hemolymph was extracted for flow-cytometric analyses of hemocyte, immune-response functions, and soft tissues were excised for histopathology. Mussels responded to P. minimum exposure with diapedesis of hemocytes into the intestine, presumably to isolate P. minimum cells within the gut, thereby minimizing damage to other tissues. This immune response appeared to have been sustained throughout the 9-day exposure period, as circulating hemocytes retained hematological and functional properties. Bacteria proliferated in the intestines of the P. minimum-exposed mussels. Hemocytes within the intestine appeared to be either overwhelmed by the large number of bacteria or fully occupied in the encapsulating response to P. minimum cells; when hemocytes reached the intestine lumina, they underwent apoptosis and bacterial degradation. This experiment demonstrated that M. edulis is affected by ingestion of toxic P. minimum; however, the specific responses observed in the blue mussel differed from those reported for other bivalve species. This finding highlights the need to study effects of HABs on different bivalve species, rather than inferring that results from one species reflect the exposure responses of all bivalves.

Experimental exposure of the blue mussel (Mytilus edulis, L.) to the toxic dinoflagellate Alexandrium fundyense: Histopathology, immune responses, and recovery

Mussels (Mytilus edulis) were exposed to cultures of the toxic dinoflagellate Alexandrium fundyense or the non-toxic alga Rhodomonas sp. to evaluate the effects of the harmful alga on the mussels and to study recovery after discontinuation of the A. fundyense exposure. Mussels were exposed for 9 days to the different algae and then all were fed Rhodomonas sp. for 6 more days. Samples of hemolymph for hemocyte analyses and tissues for histology were collected before the exposure and periodically during exposure and recovery periods. Mussels filtered and ingested both microalgal cultures, producing fecal pellets containing degraded, partially degraded, and intact cells of both algae. Mussels exposed to A. fundyense had an inflammatory response consisting of degranulation and diapedesis of hemocytes into the alimentary canal and, as the exposure continued, hemocyte migration into the connective tissue between the gonadal follicles. Evidence of lipid peroxidation, similar to the detoxification pathway described for various xenobiotics, was found; insoluble lipofuchsin granules formed (ceroidosis), and hemocytes carried the granules to the alimentary canal, thus eliminating putative dinoflagellate toxins in feces. As the number of circulating hemocytes in A. fundyense-exposed mussels became depleted, mussels were immunocompromised, and pathological changes followed, i.e., increased prevalences of ceroidosis and trematodes after 9 days of exposure. Moreover, the total number of pathological changes increased from the beginning of the exposure until the last day (day 9). After 6 days of the exposure, mussels in one of the three tanks exposed to A. fundyense mass spawned; these mussels showed more severe effects of the toxic algae than non-spawning mussels exposed to A. fundyense. No significant differences were found between the two treatments during the recovery period, indicating rapid homeostatic processes in tissues and circulating hemocytes.

Rigorous derivation of a nonlinear diffusion equation as fast-reaction limit of a continuous coagulation-fragmentation model with diffusion

Weak solutions of the spatially inhomogeneous (diffusive) Aizenmann-Bak model of coagulation-breakup within a bounded domain with homogeneous Neumann boundary conditions are shown to converge, in the fast reaction limit, towards local equilibria determined by their mass. Moreover, this mass is the solution of a nonlinear diffusion equation whose nonlinearity depends on the (size-dependent) diffusion coefficient. Initial data are assumed to have integrable zero order moment and square integrable first order moment in size, and finite entropy. In contrast to our previous result [CDF2], we are able to show the convergence without assuming uniform bounds from above and below on the number density of clusters.

Fast screening and confirmation of doping agents by UHPLC-QTOF-MS/MS

Introduction: The general strategy to perform anti-doping analysis starts with a screening followed by a confirmatory step when a sample is suspected to be positive. The screening step should be fast, generic and able to highlight any sample that may contain a prohibited substance by avoiding false negative and reducing false positive results. The confirmatory step is a dedicated procedure comprising a selective sample preparation and detection mode. Aim: The purpose of the study is to develop rapid screening and selective confirmatory strategies to detect and identify 103 doping agents in urine. Methods: For the screening, urine samples were simply diluted by a factor 2 with ultra-pure water and directly injected ("dilute and shoot") in the ultrahigh- pressure liquid chromatography (UHPLC). The UHPLC separation was performed in two gradients (ESI positive and negative) from 5/95 to 95/5% of MeCN/Water containing 0.1% formic acid. The gradient analysis time is 9 min including 3 min reequilibration. Analytes detection was performed in full scan mode on a quadrupole time-of-flight (QTOF) mass spectrometer by acquiring the exact mass of the protonated (ESI positive) or deprotonated (ESI negative) molecular ion. For the confirmatory analysis, urine samples were extracted on SPE 96-well plate with mixed-mode cation (MCX) for basic and neutral compounds or anion exchange (MAX) sorbents for acidic molecules. The analytes were eluted in 3 min (including 1.5 min reequilibration) with a S1-25 Ann Toxicol Anal. 2009; 21(S1) Abstracts gradient from 5/95 to 95/5% of MeCN/Water containing 0.1% formic acid. Analytes confirmation was performed in MS and MS/MS mode on a QTOF mass spectrometer. Results: In the screening and confirmatory analysis, basic and neutral analytes were analysed in the positive ESI mode, whereas acidic compounds were analysed in the negative mode. The analyte identification was based on retention time (tR) and exact mass measurement. "Dilute and shoot" was used as a generic sample treatment in the screening procedure, but matrix effect (e.g., ion suppression) cannot be avoided. However, the sensitivity was sufficient for all analytes to reach the minimal required performance limit (MRPL) required by the World Anti Doping Agency (WADA). To avoid time-consuming confirmatory analysis of false positive samples, a pre-confirmatory step was added. It consists of the sample re-injection, the acquisition of MS/MS spectra and the comparison to reference material. For the confirmatory analysis, urine samples were extracted by SPE allowing a pre-concentration of the analyte. A fast chromatographic separation was developed as a single analyte has to be confirmed. A dedicated QTOF-MS and MS/MS acquisition was performed to acquire within the same run a parallel scanning of two functions. Low collision energy was applied in the first channel to obtain the protonated molecular ion (QTOF-MS), while dedicated collision energy was set in the second channel to obtain fragmented ions (QTOF-MS/MS). Enough identification points were obtained to compare the spectra with reference material and negative urine sample. Finally, the entire process was validated and matrix effects quantified. Conclusion: Thanks to the coupling of UHPLC with the QTOF mass spectrometer, high tR repeatability, sensitivity, mass accuracy and mass resolution over a broad mass range were obtained. The method was sensitive, robust and reliable enough to detect and identify doping agents in urine. Keywords: screening, confirmatory analysis, UHPLC, QTOF, doping agents

Central and Cortical Gray Mater Segmentation of Magnetic Resonance Images of the Fetal Brain

Motivation. The study of human brain development in itsearly stage is today possible thanks to in vivo fetalmagnetic resonance imaging (MRI) techniques. Aquantitative analysis of fetal cortical surfacerepresents a new approach which can be used as a markerof the cerebral maturation (as gyration) and also forstudying central nervous system pathologies [1]. However,this quantitative approach is a major challenge forseveral reasons. First, movement of the fetus inside theamniotic cavity requires very fast MRI sequences tominimize motion artifacts, resulting in a poor spatialresolution and/or lower SNR. Second, due to the ongoingmyelination and cortical maturation, the appearance ofthe developing brain differs very much from thehomogenous tissue types found in adults. Third, due tolow resolution, fetal MR images considerably suffer ofpartial volume (PV) effect, sometimes in large areas.Today extensive efforts are made to deal with thereconstruction of high resolution 3D fetal volumes[2,3,4] to cope with intra-volume motion and low SNR.However, few studies exist related to the automatedsegmentation of MR fetal imaging. [5] and [6] work on thesegmentation of specific areas of the fetal brain such asposterior fossa, brainstem or germinal matrix. Firstattempt for automated brain tissue segmentation has beenpresented in [7] and in our previous work [8]. Bothmethods apply the Expectation-Maximization Markov RandomField (EM-MRF) framework but contrary to [7] we do notneed from any anatomical atlas prior. Data set &Methods. Prenatal MR imaging was performed with a 1-Tsystem (GE Medical Systems, Milwaukee) using single shotfast spin echo (ssFSE) sequences (TR 7000 ms, TE 180 ms,FOV 40 x 40 cm, slice thickness 5.4mm, in plane spatialresolution 1.09mm). Each fetus has 6 axial volumes(around 15 slices per volume), each of them acquired inabout 1 min. Each volume is shifted by 1 mm with respectto the previous one. Gestational age (GA) ranges from 29to 32 weeks. Mother is under sedation. Each volume ismanually segmented to extract fetal brain fromsurrounding maternal tissues. Then, in-homogeneityintensity correction is performed using [9] and linearintensity normalization is performed to have intensityvalues that range from 0 to 255. Note that due tointra-tissue variability of developing brain someintensity variability still remains. For each fetus, ahigh spatial resolution image of isotropic voxel size of1.09 mm is created applying [2] and using B-splines forthe scattered data interpolation [10] (see Fig. 1). Then,basal ganglia (BS) segmentation is performed on thissuper reconstructed volume. Active contour framework witha Level Set (LS) implementation is used. Our LS follows aslightly different formulation from well-known Chan-Vese[11] formulation. In our case, the LS evolves forcing themean of the inside of the curve to be the mean intensityof basal ganglia. Moreover, we add local spatial priorthrough a probabilistic map created by fitting anellipsoid onto the basal ganglia region. Some userinteraction is needed to set the mean intensity of BG(green dots in Fig. 2) and the initial fitting points forthe probabilistic prior map (blue points in Fig. 2). Oncebasal ganglia are removed from the image, brain tissuesegmentation is performed as described in [8]. Results.The case study presented here has 29 weeks of GA. Thehigh resolution reconstructed volume is presented in Fig.1. The steps of BG segmentation are shown in Fig. 2.Overlap in comparison with manual segmentation isquantified by the Dice similarity index (DSI) equal to0.829 (values above 0.7 are considered a very goodagreement). Such BG segmentation has been applied on 3other subjects ranging for 29 to 32 GA and the DSI hasbeen of 0.856, 0.794 and 0.785. Our segmentation of theinner (red and blue contours) and outer cortical surface(green contour) is presented in Fig. 3. Finally, torefine the results we include our WM segmentation in theFreesurfer software [12] and some manual corrections toobtain Fig.4. Discussion. Precise cortical surfaceextraction of fetal brain is needed for quantitativestudies of early human brain development. Our workcombines the well known statistical classificationframework with the active contour segmentation forcentral gray mater extraction. A main advantage of thepresented procedure for fetal brain surface extraction isthat we do not include any spatial prior coming fromanatomical atlases. The results presented here arepreliminary but promising. Our efforts are now in testingsuch approach on a wider range of gestational ages thatwe will include in the final version of this work andstudying as well its generalization to different scannersand different type of MRI sequences. References. [1]Guibaud, Prenatal Diagnosis 29(4) (2009). [2] Rousseau,Acad. Rad. 13(9), 2006, [3] Jiang, IEEE TMI 2007. [4]Warfield IADB, MICCAI 2009. [5] Claude, IEEE Trans. Bio.Eng. 51(4) (2004). [6] Habas, MICCAI (Pt. 1) 2008. [7]Bertelsen, ISMRM 2009 [8] Bach Cuadra, IADB, MICCAI 2009.[9] Styner, IEEE TMI 19(39 (2000). [10] Lee, IEEE Trans.Visual. And Comp. Graph. 3(3), 1997, [11] Chan, IEEETrans. Img. Proc, 10(2), 2001 [12] Freesurfer,http://surfer.nmr.mgh.harvard.edu.

Real-time imaging of regional myocardial function using fast-SENC.

A technique for fast imaging of regional myocardial function using a spiral acquisition in combination with strain-encoded (SENC) magnetic resonance imaging (MRI) is presented in this paper. This technique, which is termed fast-SENC, enables scan durations as short as a single heartbeat. A reduced field of view (FOV) without foldover artifacts was achieved by localized SENC, which selectively excited the region around the heart. The two images required for SENC imaging (low- and high-tuning) were acquired in an interleaved fashion throughout the cardiac cycle to further shorten the scan time. Regional circumferential contraction and longitudinal shortening of both the left ventricle (LV) and right ventricle (RV) were examined in long- and short-axis views, respectively. The in vivo results obtained from five human subjects and five infarcted dogs are presented. The results of the fast-SENC technique in a single heartbeat acquisition were comparable to those obtained by conventional SENC in a long acquisition time. Therefore, fast-SENC may prove useful for imaging during stress or arrhythmia.

Apport du bleu de toluidine en cancérologie bucco-pharyngo-oesophagienne [Contribution of Toluidine blue to bucco-pharyngo-oesophageal cancerology]

En muqueuse malpighienne lisse, le Bleu de Toluidine présente un intérêt diagnostic considérable, non seulement lors de l'endoscopie, mais lors de la chirurgie d'exérèse des carcinomes intra-épithéliaux et micro-invasifs de la voie digestive supérieure (bouche, pharynx, oesophage). L'utilisation systématique de ce colorant vital nous a permis le diagnostic et l'étude morphologique endoscopique de 18 carcinomes « précoces » de l'oesophage (9 in situ, 3 dysplasies sévères, 1 microinvasif, 5 limités à la sous-muqueuse). La multicentricité et la multifocalité de ces lésions est fréquente (80 à 90 % des cas), ce qui impose leur exérèse chirurgicale sous Bleu de Toluidine, afin de détecter les foyers de carcinome in situ à distance de la lésion primaire. La connaissance exacte des faux positifs et négatifs permet de limiter les erreurs d'interprétation lors de la coloration.

Survival, development and growth of larvae of the blue swimmer crab, Portunus pelagicus, cultured under different photoperiod conditions

The blue swimmer crab is a commercially important species of the tropical Indo-Pacific regions that shows substantial potential as a candidate species for aquaculture. Optimization of larval rearing conditions, including photoperiod, is therefore important to establish a method for the intensive hatchery culture of this species. Newly hatched larvae of Portunuspelagicus in first zoeal stage (ZI) were reared under five photoperiod regimes 0L: 24D, 6L: 18D, 12L: 12D, 18L: 6D, and 24L: 0D (5 replicates per treatment) till they metamorphosed to megalopae (ranged from 8.5 ± 0.3 days (18L: 6D) to 10.8 ± 1.8 days (0L: 24D) at 29 ± 1 °C). Daily, larvae of each treatment were fed an identical diet of mixed rotifer and Artemia nauplii, and the survival and molt to successive stages was monitored. Newly hatched ZI larvae of P. pelagicus could successfully develop to the megalopal stage under all tested photoperiod conditions, but we detected significant differences in survival among treatments (p & 0.05). The constant darkness treatment (0L: 24D) had the lowest (19.2 ± 7.2%, mean ± S.E.) cumulative survival from ZI to the megalopal stage, while the 18L: 6D treatment achieved the highest survival (51.2 ± 23.6%). Similarly, the photoperiod significantly affected zoeal development. Constant darkness led to the longest cumulative zoeal duration (10.8 ± 1.8 days), whereas the 18L: 6D treatment rendered the shortest larval development (8.5 ± 0.3 days). In addition, larvae reared under constant darkness resulted in the smallest megalopae (carapace length = 1.44 ± 0.09 mm) and the lowest dry weight (0.536 ± 0.188 mg). In conclusion, photoperiod significantly affected the survival, development, and growth of P. pelagicus zoeal larvae. Constant darkness led to the lowest larval survival and developmental rate, while a photoperiod regime of 18L: 6D appeared to be the most suitable condition for the rearing of zoeal larvae of P. pelagicus.

Characterization of pupil responses to blue and red light stimuli in autosomal dominant retinitis pigmentosa due to NR2E3 mutation.

PURPOSE: We characterized the pupil responses that reflect rod, cone, and melanopsin function in a genetically homogeneous cohort of patients with autosomal dominant retinitis pigmentosa (adRP). METHODS: Nine patients with Gly56Arg mutation of the NR2E3 gene and 12 control subjects were studied. Pupil and subjective visual responses to red and blue light flashes over a 7 log-unit range of intensities were recorded under dark and light adaptation. The pupil responses were plotted against stimulus intensity to obtain red-light and blue-light response curves. RESULTS: In the dark-adapted blue-light stimulus condition, patients showed significantly higher threshold intensities for visual perception and for a pupil response compared to controls (P = 0.02 and P = 0.006, respectively). The rod-dependent, blue-light pupil responses decreased with disease progression. In contrast, the cone-dependent pupil responses (light-adapted red-light stimulus condition) did not differ between patients and controls. The difference in the retinal sensitivity to blue and red stimuli was the most sensitive parameter to detect photoreceptor dysfunction. Unexpectedly, the melanopsin-mediated pupil response was decreased in patients (P = 0.02). CONCLUSIONS: Pupil responses of patients with NR2E3-associated adRP demonstrated reduced retinal sensitivity to dim blue light under dark adaptation, presumably reflecting decreased rod function. Rod-dependent pupil responses were quantifiable in all patients, including those with non-recordable scotopic electroretinogram, and correlated with the extent of clinical disease. Thus, the chromatic pupil light reflex can be used to monitor photoreceptor degeneration over a larger range of disease progression compared to standard electrophysiology.

Predictive value of the acid fast smear for detection of Mycobacterium tuberculosis in respiratory specimens in a Reference Center of HIV/Aids in Rio de Janeiro, Brazil

In order to evaluate the predictive value of acid fast bacilii (AFB) smear for the diagnosis of Mycobacterium tuberculosis in respiratory specimens in a setting with a high prevalence of Aids and an unknown prevalence of nontuberculous mycobacteria (NTM), we retrospectively examined specimens cultured for mycobacteria between 1 September 1993 and 30 September 1994 and medical records of patients with positive culture in a General Hospital, Aids reference in Rio de Janeiro, Brazil. Seventy three per cent (1517/2077) of samples were respiratory specimens and mycobacteria were recovered from 20.6% (313/1517) of these. M. tuberculosis was identified in 94.2% (295/313) and NTM in 5.8% (18/313). The yield of positive AFB smear and of positive culture was 6.1% (93/1517) and 20.6% (313/1517), respectively. The positive predictive value (PPV) of AFB for M. tuberculosis was 98.4% in expectorated sputum and 96.4% in bronchoalveolar lavage. Forty four percent (130/295) of specimens with positive culture for M. tuberculosis and 66.7% (12/18) for NTM were from patients HIV positive. The conclusion was that in our study population, the PPV of AFB for M. tuberculosis in respiratory specimens was high and the prevalence of NTM was low despite the high prevalence of HIV positive.