Control d'un sistema de marcatge làser de caixes

RESUM Avui en dia, l’alta competitivitat que existeix al mercat, fa que les empreses hagin d’esprémer al màxim les seves possibilitats per no quedar-se enrere. Un dels processos en que aquest fet hi és més present és el productiu. L’empresa JCM Technologies també engloba aquest camp i és en un dels seus processos productius on aquest projecte pren part. L’objectiu d’aquest projecte final de carrera ha estat desenvolupar un sistema per poder marcar caixes mitjançant un làser de CO2 i un automatisme manipulador de caixes. D’aquesta manera aquest procés productiu té una durada molt inferior a l’antic procés, que consistia en enganxar una etiqueta al lloc on ara és marcat pel làser. Per satisfer els objectius, s’ha creat una aplicació de Windows que per mitjà d’una interfície gràfica, permet a l’usuari realitzar els passos necessaris per fer el marcatge. Primerament es recullen les dades procedents de la comanda; seguidament es seleccionen les que s’han de marcar a les caixes i s’envien al làser mitjançant comunicació sèrie; una vegada aquesta inicialització ha finalitzat correctament, s’engega la seqüència de marcatge de les caixes, que en marcarà la quantitat indicada. Aquest procés de marcatge consisteix en supervisar l’estat en que es troben certes senyals, procedents de l’automatisme i del làser, i depenent d’aquestes generar-ne unes altres. Aconseguint així realitzar el procés de marcatge de cada caixa. Com a conclusions cal dir, que els objectius s’han complert, ja que s’ha aconseguit un procés de marcatge ràpid i robust. També s’ha aconseguit que les parts de configuració de l’aplicació, i de les caixes siguin de fàcil manipulació. Per tant, amb l’acompliment dels objectius d’aquest projecte, aconseguim completar el sistema de marcatge, i permetre que les caixes siguin marcades de forma correcta, ràpida i eficient.

Contribution of the gap junction proteins Connexin40 and Connexin43 to the control of blood pressure

Summary Cells in tissues and organs coordinate their activities by communicating with each other through intercellular channels named gap junctions. These channels are conduits between the cytoplasmic compartments of adjacent cells, allowing the exchange of small molecules which may be crucial for hormone secretion. Renin is normally secreted in a regulated manner by specific cells of the juxtaglomerular apparatus located within the renal cortex. Gap junctional communication may be requisite to maintain an accurate functioning in coordination of renin-producing cells, more especially as renin is of paramount importance for the control of blood pressure. Connexin43 (Cx43) and Cx40 form gap junctions that link in vivo the cells of the juxtaglomerular apparatus. Cx43 links the endothelial cells, whereas gap junctions made of Cx40 connect the endothelial cells, the renin secreting cells, as well as the endothelial cells of to the renin-secreting cells of the afferent arteriole. The observation that loss of Cx40 results in chronic hypertension associated with altered vasomotion and signal conduction along arterioles, has lead us to suggest that connexins may contribute to control blood pressure by participating to the integration of various mechanical, osmotic and electrochemical stimuli involved in the control of renin secretion and by mediating the adaptive changes of the vascular wall induced by elevated blood pressure and mechanical stress. We therefore postulated that the absence of Cx40 could have deleterious effects on the coordinated functioning of the renin-containing cells, hence accounting for hypertension. In the first part of my thesis, we reported that Cx40-deficient mice (Cx40) are hypertensive due to increased plasma renin levels and numbers of renin-producing cells. Besides, we demonstrated that prostaglandins and nitric oxide, which are possible mediators in the regulation of renin secretion by the macula densa, exert a critical role in the mechanisms controlling blood pressure ín Cx40 knockout hypertensive mice. In view of previous studies that stated avessel-specifc increase in the expression of Cx43 during renin-dependent hypertension, we hypothesized that Cx43 channels are particularly well-matched to integrate the response of cells constituting the vascular wall to hypertensive conditions. Using transgenic mice in which Cx43 was replaced by Cx32, we revealed that the replacement of Cx43 by Cx32 is associated with decreased expression and secretion of renin and prevent the renin-dependent hypertension which is normally induced in the 2K1C model. To gain insights into the regulation of connexins in two separate tissues exposed to the same fluid pressure, the second part of my thesis work was dedicated to the study of the impact of chronic hypertension and related hypertrophy on the expression of the cardiovascular connexins (Cx40, Cx37, Cx43 and Cx45) in mouse aorta and heart. Our results documented that the expression of connexins is differentially regulated in mouse aorta. according to the models of hypertension. Thus, blood pressure induces mechanical forces that differentially alter the expression of vascular connexins in order to respond to an adaptation of the aortic wall observed under pathological conditions. Altogether these data provide the first evidences that intercellular communication mediated by gap junctions is required for a proper renin secretion from the juxtaglomerular apparatus in order to control blood pressure.

Gas identification with tin oxide sensor array and self-organizing maps: adaptive correction of sensor drifts

Low-cost tin oxide gas sensors are inherently nonspecific. In addition, they have several undesirable characteristics such as slow response, nonlinearities, and long-term drifts. This paper shows that the combination of a gas-sensor array together with self-organizing maps (SOM's) permit success in gas classification problems. The system is able to determine the gas present in an atmosphere with error rates lower than 3%. Correction of the sensor's drift with an adaptive SOM has also been investigated

A regulatory network for the efficient control of transgene expression.

BACKGROUND: Expression of heterologous genes in mammalian cells or organisms for therapeutic or experimental purposes often requires tight control of transgene expression. Specifically, the following criteria should be met: no background gene activity in the off-state, high gene expression in the on-state, regulated expression over an extended period, and multiple switching between on- and off-states. METHODS: Here, we describe a genetic switch system for controlled transgene transcription using chimeric repressor and activator proteins functioning in a novel regulatory network. In the off-state, the target transgene is actively silenced by a chimeric protein consisting of multimerized eukaryotic transcriptional repression domains fused to the DNA-binding tetracycline repressor. In the on-state, the inducer drug doxycycline affects both the derepression of the target gene promoter and activation by the GAL4-VP16 transactivator, which in turn is under the control of an autoregulatory feedback loop. RESULTS: The hallmark of this new system is the efficient transgene silencing in the off-state, as demonstrated by the tightly controlled expression of the highly cytotoxic diphtheria toxin A gene. Addition of the inducer drug allows robust activation of transgene expression. In stably transfected cells, this control is still observed after months of repeated cycling between the repressed and activated states of the target genes. CONCLUSIONS: This system permits tight long-term regulation when stably introduced into cell lines. The underlying principles of this network system should have general applications in biotechnology and gene therapy.

Accurate and efficient simulation of multiphase flow in a heterogeneous reservoir by using error estimate and control in the multiscale finite-volume framework

The multiscale finite-volume (MSFV) method is designed to reduce the computational cost of elliptic and parabolic problems with highly heterogeneous anisotropic coefficients. The reduction is achieved by splitting the original global problem into a set of local problems (with approximate local boundary conditions) coupled by a coarse global problem. It has been shown recently that the numerical errors in MSFV results can be reduced systematically with an iterative procedure that provides a conservative velocity field after any iteration step. The iterative MSFV (i-MSFV) method can be obtained with an improved (smoothed) multiscale solution to enhance the localization conditions, with a Krylov subspace method [e.g., the generalized-minimal-residual (GMRES) algorithm] preconditioned by the MSFV system, or with a combination of both. In a multiphase-flow system, a balance between accuracy and computational efficiency should be achieved by finding a minimum number of i-MSFV iterations (on pressure), which is necessary to achieve the desired accuracy in the saturation solution. In this work, we extend the i-MSFV method to sequential implicit simulation of time-dependent problems. To control the error of the coupled saturation/pressure system, we analyze the transport error caused by an approximate velocity field. We then propose an error-control strategy on the basis of the residual of the pressure equation. At the beginning of simulation, the pressure solution is iterated until a specified accuracy is achieved. To minimize the number of iterations in a multiphase-flow problem, the solution at the previous timestep is used to improve the localization assumption at the current timestep. Additional iterations are used only when the residual becomes larger than a specified threshold value. Numerical results show that only a few iterations on average are necessary to improve the MSFV results significantly, even for very challenging problems. Therefore, the proposed adaptive strategy yields efficient and accurate simulation of multiphase flow in heterogeneous porous media.

Adaptive tracking of EEG oscillations.

Neuronal oscillations are an important aspect of EEG recordings. These oscillations are supposed to be involved in several cognitive mechanisms. For instance, oscillatory activity is considered a key component for the top-down control of perception. However, measuring this activity and its influence requires precise extraction of frequency components. This processing is not straightforward. Particularly, difficulties with extracting oscillations arise due to their time-varying characteristics. Moreover, when phase information is needed, it is of the utmost importance to extract narrow-band signals. This paper presents a novel method using adaptive filters for tracking and extracting these time-varying oscillations. This scheme is designed to maximize the oscillatory behavior at the output of the adaptive filter. It is then capable of tracking an oscillation and describing its temporal evolution even during low amplitude time segments. Moreover, this method can be extended in order to track several oscillations simultaneously and to use multiple signals. These two extensions are particularly relevant in the framework of EEG data processing, where oscillations are active at the same time in different frequency bands and signals are recorded with multiple sensors. The presented tracking scheme is first tested with synthetic signals in order to highlight its capabilities. Then it is applied to data recorded during a visual shape discrimination experiment for assessing its usefulness during EEG processing and in detecting functionally relevant changes. This method is an interesting additional processing step for providing alternative information compared to classical time-frequency analyses and for improving the detection and analysis of cross-frequency couplings.

An adaptive multiscale method for density-driven instabilities

Accurate modeling of flow instabilities requires computational tools able to deal with several interacting scales, from the scale at which fingers are triggered up to the scale at which their effects need to be described. The Multiscale Finite Volume (MsFV) method offers a framework to couple fine-and coarse-scale features by solving a set of localized problems which are used both to define a coarse-scale problem and to reconstruct the fine-scale details of the flow. The MsFV method can be seen as an upscaling-downscaling technique, which is computationally more efficient than standard discretization schemes and more accurate than traditional upscaling techniques. We show that, although the method has proven accurate in modeling density-driven flow under stable conditions, the accuracy of the MsFV method deteriorates in case of unstable flow and an iterative scheme is required to control the localization error. To avoid large computational overhead due to the iterative scheme, we suggest several adaptive strategies both for flow and transport. In particular, the concentration gradient is used to identify a front region where instabilities are triggered and an accurate (iteratively improved) solution is required. Outside the front region the problem is upscaled and both flow and transport are solved only at the coarse scale. This adaptive strategy leads to very accurate solutions at roughly the same computational cost as the non-iterative MsFV method. In many circumstances, however, an accurate description of flow instabilities requires a refinement of the computational grid rather than a coarsening. For these problems, we propose a modified iterative MsFV, which can be used as downscaling method (DMsFV). Compared to other grid refinement techniques the DMsFV clearly separates the computational domain into refined and non-refined regions, which can be treated separately and matched later. This gives great flexibility to employ different physical descriptions in different regions, where different equations could be solved, offering an excellent framework to construct hybrid methods.

Feedback control of DnaA-mediated replication initiation by replisome-associated HdaA protein in Caulobacter.

Chromosome replication in Caulobacter crescentus is tightly regulated to ensure that initiation occurs at the right time and only once during the cell cycle. The timing of replication initiation is controlled by both CtrA and DnaA. CtrA binds to and silences the origin. Upon the clearance of CtrA from the cell, the DnaA protein accumulates and allows loading of the replisome at the origin. Here, we identify an additional layer of replication initiation control that is mediated by the HdaA protein. In Escherichia coli, the Hda protein inactivates DnaA after replication initiation. We show that the Caulobacter HdaA homologue is necessary to restrict the initiation of DNA replication to only once per cell cycle and that it dynamically colocalizes with the replisome throughout the cell cycle. Moreover, the transcription of hdaA is directly activated by DnaA, providing a robust feedback regulatory mechanism that adjusts the levels of HdaA to inactivate DnaA.

Robust accelerated failure time regression

Robust estimators for accelerated failure time models with asymmetric (or symmetric) error distribution and censored observations are proposed. It is assumed that the error model belongs to a log-location-scale family of distributions and that the mean response is the parameter of interest. Since scale is a main component of mean, scale is not treated as a nuisance parameter. A three steps procedure is proposed. In the first step, an initial high breakdown point S estimate is computed. In the second step, observations that are unlikely under the estimated model are rejected or down weighted. Finally, a weighted maximum likelihood estimate is computed. To define the estimates, functions of censored residuals are replaced by their estimated conditional expectation given that the response is larger than the observed censored value. The rejection rule in the second step is based on an adaptive cut-off that, asymptotically, does not reject any observation when the data are generat ed according to the model. Therefore, the final estimate attains full efficiency at the model, with respect to the maximum likelihood estimate, while maintaining the breakdown point of the initial estimator. Asymptotic results are provided. The new procedure is evaluated with the help of Monte Carlo simulations. Two examples with real data are discussed.

The adaptive function of melanin-based coloration in the tawny owl (Strix aluco) ; the proximate role of the melanocortin system

Colour polymorphism is common in wild population. One of the main questioning of evolutionary biologists is to understand how different colour variants could have evolved and be maintained in fluctuating environments, a selective process that forces individuals to constantly adapt their strategies in order to survive. This issue is particularly true for traits that are genetically inherited. Natural selection erodes genotypes with lowest fitness (less adapted), reducing in turn global genetic variation within population. In this context, the study of the evolution and maintenance of melanin- based coloration is relevant since inter-individual variation in the deposition of these pigments is common in animal and plant kingdoms and is under strong genetic control. In this thesis, I focus on the specific case of the tawny owl (Strix aluco), a species displaying continuous variation in reddish pheomelanin-based coloration. Interestingly, empirical studies highlighted covariations between melanin-based coloration and important behavioural, physiological and life history traits. Recently, a genetic model pointed out the melanocortin system and their pleiotropic effects as a potential regulator of these covariations. Accordingly, this PhD thesis further investigates colour-specific behavioural, physiological, or life history strategies, while examining the proximate mechanisms underlying these reaction norms. We found that differently coloured tawny owls differently resolve fundamental trade-off between offspring number and quality (Chapter 1), light melanic individuals producing many low- quality offspring and dark, melanic ones producing few high-quality offspring. These reproductive strategies are likely to induce alternative physiological constraints. Indeed, we demonstrated that light melanic individuals produced higher levels of reactive oxygen species (ROS, Chapter 2), but also expressed higher levels of antioxidant (GSH, Chapters 2 & 3). Interestingly, we showed that light melanic breeding females could modulate their POMC prohormone levels according to the environmental conditions, while dark reddish ones produced constant levels of this prohormone {Chapter 4). Finally, we highlighted colour-specific patterns of prohormone convertase 1 (PCI) gene expression (Chapter 5), an enzyme responsible for POMC prohormone processing to ACTH and a- MSH, for instance. Altogether, these results provide strong evidence of colour-specific strategies, light and melanic tawny owls better coping with stressful and relaxed environments, respectively. Variation in melanin-based coloration is likely to be maintained by the heterogeneity of our study area and strong environmental stochasticity within and between years, these process favouring differently coloured tawny owls at different periods of time. From a proximate point of view, this PhD thesis supports the hypothesis that covariations between phenotypic traits and melanin-based coloration stems from the melanocortin system, especially the fundamental role of POMC gene expression and its processing to melanocortin peptides. - Le polymorphisme de couleur est une variation phénotypique très fréquente dans la nature. En biologie évolutive, une des problématiques clés est donc de comprendre comment différent morphes de couleur peuvent être apparus et maintenus au cours du temps dans des environnements aussi variables que les nôtres, surtout que ces fluctuations forcent ces morphes à s'adapter constamment pour assurer leur survie. Cette thématique est particulièrement réelle lorsque les variations phénotypiques sont héréditaires et donc sous forte influence génétique. La sélection naturelle a en effet le pouvoir d'éroder rapidement la variation génétique en éliminant les génotypes mal adaptés. Dans ce sens, l'étude de l'évolution, et de la maintenance de la coloration mélanique est donc tout à fait pertinente car la variation de coloration entre individus est très répandue à travers les règnes animal et végétal et sous forte influence génétique. Dans cette thèse, je me suis concentré sur le cas spécifique de la chouette hulotte (Strix aluco), une espèce présentant une variation continue dans la déposition de pigments pheomélaniques roux. De précédentes études ont déjà montré que cette variation de coloration était associée avec des variations de traits comportementaux, physiologiques ou d'histoire de vie. Récemment, une étude a souligné l'importance du système des mélanocortines et de leurs effets pléiotropes dans la régulation de ces covariations. En conséquence, cette thèse de doctoral a pour but d'étudier un peu plus les stratégies comportementales, physiologiques ou d'histoire de vie spécifiques à chaque morphe de couleur, tout en examinant un peu plus les mécanismes proximaux potentiellement à la base de ces normes de réactions. Nous constatons tout d'abord que les morphes de couleurs étaient associés à différentes stratégies dans la résolution de compromis telle que la production de beaucoup de jeunes ou des jeunes de qualité (Chapitre 1). Les morphes gris (dit peu mélaniques) ont tendance à produire beaucoup de jeunes mains de moindre qualité, alors que les morphes roux (dit fortement mélaniques) produisent moins de jeunes mais de meilleure qualité. Ces stratégies sont susceptibles alors d'induire certaines contraintes physiologiques. Par exemple, nous montrons que les morphes gris produisent plus de dérivés réactifs de l'oxygène (ROS, Chapitre 2), mais aussi plus d'antioxydants (GSH, Chapitres 2 & 3). Nous montrons ensuite que les femelles grises ont une plus grande capacité à moduler leur niveau de POMC prohormone dans le sang en fonction des conditions environnementales, alors que les femelles rousses gardent un niveau constant (Chapitre 4). Finalement, nous démontrons que les patterns d'expression du gène codant pour la prohormone convertase 1 varient chez des jeunes issus de parents gris ou roux (Chapitre 5). Ceci est particulièrement intéressant car cette enzyme permet de scinder la POMC prohormone en plusieurs peptides importants tels que l'ACTH ou l'a-MSH. En conclusion, ces résultats démontrent qu'il y a bel et bien des stratégies évolutives différentes entre les morphes de couleurs, les chouettes hulottes grises et rousses étant respectivement plus adaptés à des environnements stressants ou favorables. L'hétérogénéité de notre zone d'étude et la stochasticité environnementale qui caractérise ses habitats pourraient donc agir comme une source de sélection temporelle, laquelle favoriserait les différents morphes de couleurs à diverses périodes. D'un point de vue plus proximale maintenant, cette thèse de doctorat soutient l'hypothèse que les covariations observées entre la coloration mélanique et des traits phénotypiques importants sont modulées par les effets pléiotropes du système des mélanocortines, et met en avant le rôle prépondérant que pourrait jouer l'expression du gène POMC et sa post traduction en mélanocortines.

Immunization with HIV Gag targeted to dendritic cells followed by recombinant New York vaccinia virus induces robust T-cell immunity in nonhuman primates.

Protein vaccines, if rendered immunogenic, would facilitate vaccine development against HIV and other pathogens. We compared in nonhuman primates (NHPs) immune responses to HIV Gag p24 within 3G9 antibody to DEC205 ("DEC-HIV Gag p24"), an uptake receptor on dendritic cells, to nontargeted protein, with or without poly ICLC, a synthetic double stranded RNA, as adjuvant. Priming s.c. with 60 μg of both HIV Gag p24 vaccines elicited potent CD4(+) T cells secreting IL-2, IFN-γ, and TNF-α, which also proliferated. The responses increased with each of three immunizations and recognized multiple Gag peptides. DEC-HIV Gag p24 showed better cross-priming for CD8(+) T cells, whereas the avidity of anti-Gag antibodies was ∼10-fold higher with nontargeted Gag 24 protein. For both protein vaccines, poly ICLC was essential for T- and B-cell immunity. To determine whether adaptive responses could be further enhanced, animals were boosted with New York vaccinia virus (NYVAC)-HIV Gag/Pol/Nef. Gag-specific CD4(+) and CD8(+) T-cell responses increased markedly after priming with both protein vaccines and poly ICLC. These data reveal qualitative differences in antibody and T-cell responses to DEC-HIV Gag p24 and Gag p24 protein and show that prime boost with protein and adjuvant followed by NYVAC elicits potent cellular immunity.

The evolution, maintenance and adaptive function of genetic colour polymorphism in birds.

The hypothesis that ornaments can honestly signal quality only if their expression is condition-dependent has dominated the study of the evolution and function of colour traits. Much less interest has been devoted to the adaptive function of colour traits for which the expression is not, or is to a low extent, sensitive to body condition and the environment in which individuals live. The aim of the present paper is to review the current theoretical and empirical knowledge of the evolution, maintenance and adaptive function of colour plumage traits for which the expression is mainly under genetic control. The finding that in many bird species the inheritance of colour morphs follows the laws of Mendel indicates that genetic colour polymorphism is frequent. Polymorphism may have evolved or be maintained because each colour morph facilitates the exploitation of alternative ecological niches as suggested by the observation that individuals are not randomly distributed among habitats with respect to coloration. Consistent with the hypothesis that different colour morphs are linked to alternative strategies is the finding that in a majority of species polymorphism is associated with reproductive parameters, and behavioural, life-history and physiological traits. Experimental studies showed that such covariations can have a genetic basis. These observations suggest that colour polymorphism has an adaptive function. Aviary and field experiments demonstrated that colour polymorphism is used as a criterion in mate-choice decisions and dominance interactions confirming the claim that conspecifics assess each other's colour morphs. The factors favouring the evolution and maintenance of genetic variation in coloration are reviewed, but empirical data are virtually lacking to assess their importance. Although current theory predicts that only condition-dependent traits can signal quality, the present review shows that genetically inherited morphs can reveal the same qualities. The study of genetic colour polymorphism will provide important and original insights on the adaptive function of conspicuous traits.

Cardiac raptor deletion impairs adaptive hypertrophy, alters metabolic gene expression and causes heart failure in mice

Background: Mammalian target of rapamycin (mTOR), a central regulator of cell growth, is found in two structurally and functionally distinct multiprotein complexes called mTOR complex (mTORC)1 and mTORC2. The specific roles of each of these branches of mTOR signaling have not been dissected in the adult heart. In the present study, we aimed to bring new insights into the function of cardiac mTORC1-mediated signaling in physiological as well as pathological situations.Methods: We generated mice homozygous for loxP-flanked raptor and positive for the tamoxifen-inducible Cre recombinase (MerCreMer) under control of the α- myosin heavy chain promoter. The raptor gene encodes an essential component of mTORC1. Gene ablation was induced at the age of 10-12 weeks, and two weeks later the raptor cardiac-knockout (raptor-cKO) mice started voluntary cagewheel exercise or were subjected to transverse aortic constriction (TAC) to induce pressure overload.Results: In sedentary raptor-cKO mice, ejection fractions gradually decreased, resulting in significantly reduced values at 38 days (P < 0.001). Raptor-cKO mice started to die during the fifth week after the last tamoxifen injection. At that time, the mortality rate was 36% in sedentary (n = 11) and 64% in exercising (n = 14) mice. TAC-induced pressure overload resulted in severe cardiac dysfunction already at earlier timepoints. Thus, at 7-9 days after surgery, ejection fraction and fractional shortening values were 22.3% vs 43.5% and 10.2% vs 21.5% in raptor-cKO vs wild-type mice, respectively. This was accompanied by significant reductions of ventricular wall and septal thickness as well as an increase in left ventricular internal diameter. Moreover, ventricular weight to tibial length ratios were increased in wild-type, but not in the raptor-cKO TAC mice. Together, this shows that raptor-cKO mice rapidly developed dilated cardiomyopathy without going through a phase of adaptive hypertrophy. Expression of ANP and β-MHC was induced in all raptor-cKO mice irrespective of the cardiac load conditions. Consistent with reduced mTORC1 activity, phosphorylation of ribosomal S6 kinase and 4E-BP1 was blunted, indicating reduced protein synthesis. Moreover, expression of multiple genes involved in the regulation of energy metabolism was altered, and followed by a shift from fatty acid to glucose oxidation.Conclusion: Our study suggests that mTORC1 coordinates protein and energy metabolic pathways in the heart. Moreover, we demonstrate that raptor is essential for the cardiac adaptation to increased workload and importantly, also for normal physiological cardiac function.

The branch-site test of positive selection is surprisingly robust but lacks power under synonymous substitution saturation and variation in GC

Positive selection is widely estimated from protein coding sequence alignments by the nonsynonymous-to-synonymous ratio omega. Increasingly elaborate codon models are used in a likelihood framework for this estimation. Although there is widespread concern about the robustness of the estimation of the omega ratio, more efforts are needed to estimate this robustness, especially in the context of complex models. Here, we focused on the branch-site codon model. We investigated its robustness on a large set of simulated data. First, we investigated the impact of sequence divergence. We found evidence of underestimation of the synonymous substitution rate for values as small as 0.5, with a slight increase in false positives for the branch-site test. When dS increases further, underestimation of dS is worse, but false positives decrease. Interestingly, the detection of true positives follows a similar distribution, with a maximum for intermediary values of dS. Thus, high dS is more of a concern for a loss of power (false negatives) than for false positives of the test. Second, we investigated the impact of GC content. We showed that there is no significant difference of false positives between high GC (up to similar to 80%) and low GC (similar to 30%) genes. Moreover, neither shifts of GC content on a specific branch nor major shifts in GC along the gene sequence generate many false positives. Our results confirm that the branch-site is a very conservative test.

Cross-frequency coupling in EEG revealed by adaptive oscillation tracking.

Introduction: Neuronal oscillations have been the focus of increasing interest in the neuroscientific community, in part because they have been considered as a possible integrating mechanism through which internal states can influence stimulus processing in a top-down way (Engel et al., 2001). Moreover, increasing evidence indicates that oscillations in different frequency bands interact with one other through coupling mechanisms (Jensen and Colgin, 2007). The existence and the importance of these cross-frequency couplings during various tasks have been verified by recent studies (Canolty et al., 2006; Lakatos et al., 2007). In this study, we measure the strength and directionality of two types of couplings - phase-amplitude couplings and phase-phase couplings - between various bands in EEG data recorded during an illusory contour experiment that were identified using a recently-proposed adaptive frequency tracking algorithm (Van Zaen et al., 2010). Methods: The data used in this study have been taken from a previously published study examining the spatiotemporal mechanisms of illusory contour processing (Murray et al., 2002). The EEG in the present study were from a subset of nine subjects. Each stimulus was composed of 'pac-man' inducers presented in two orientations: IC, when an illusory contour was present, and NC, when no contour could be detected. The signals recorded by the electrodes P2, P4, P6, PO4 and PO6 were averaged, and filtered into the following bands: 4-8Hz, 8-12Hz, 15-25Hz, 35-45Hz, 45-55Hz, 55-65Hz and 65-75Hz. An adaptive frequency tracking algorithm (Van Zaen et al., 2010) was then applied in each band in order to extract the main oscillation and estimate its frequency. This additional step ensures that clean phase information is obtained when taking the Hilbert transform. The frequency estimated by the tracker was averaged over sliding windows and then used to compare the two conditions. Two types of cross-frequency couplings were considered: phase-amplitude couplings and phase-phase couplings. Both types were measured with the phase locking value (PLV, Lachaux et al., 1999) over sliding windows. The phase-amplitude couplings were computed with the phase of the low frequency oscillation and the phase of the amplitude of the high frequency one. Different coupling coefficients were used when measuring phase-phase couplings in order to estimate different m:n synchronizations (4:3, 3:2, 2:1, 3:1, 4:1, 5:1, 6:1, 7:1, 8:1 and 9:1) and to take into account the frequency differences across bands. Moreover, the direction of coupling was estimated with a directionality index (Bahraminasab et al., 2008). Finally, the two conditions IC and NC were compared with ANOVAs with 'subject' as a random effect and 'condition' as a fixed effect. Before computing the statistical tests, the PLV values were transformed into approximately normal variables (Penny et al., 2008). Results: When comparing the mean estimated frequency across conditions, a significant difference was found only in the 4-8Hz band, such that the frequency within this band was significantly higher for IC than NC stimuli starting at ~250ms post-stimulus onset (Fig. 1; solid line shows IC and dashed line NC). Significant differences in phase-amplitude couplings were obtained only when the 4-8 Hz band was taken as the low frequency band. Moreover, in all significant situations, the coupling strength is higher for the NC than IC condition. An example of significant difference between conditions is shown in Fig. 2 for the phase-amplitude coupling between the 4-8Hz and 55-65Hz bands (p-value in top panel and mean PLV values in the bottom panel). A decrease in coupling strength was observed shortly after stimulus onset for both conditions and was greater for the condition IC. This phenomenon was observed with all other frequency bands. The results obtained for the phase-phase couplings were more complex. As for the phase-amplitude couplings, all significant differences were obtained when the 4-8Hz band was considered as the low frequency band. The stimulus condition exhibiting the higher coupling strength depended on the ratio of the coupling coefficients. When this ratio was small, the IC condition exhibited the higher phase-phase coupling strength. When this ratio was large, the NC condition exhibited the higher coupling strength. Fig. 3 shows the phase-phase couplings between the 4-8Hz and 35-45Hz bands for the coupling coefficient 6:1, and the coupling strength was significantly higher for the IC than NC condition. By contrast, for the coupling coefficient 9:1 the NC condition gave the higher coupling strength (Fig. 4). Control analyses verified that it is not a consequence of the frequency difference between the two conditions in the 4-8Hz band. The directionality measures indicated a transfer of information from the low frequency components towards the high frequency ones. Conclusions: Adaptive tracking is a feasible method for EEG analyses, revealing information both about stimulus-related differences and coupling patterns across frequencies. Theta oscillations play a central role in illusory shape processing and more generally in visual processing. The presence vs. absence of illusory shapes was paralleled by faster theta oscillations. Phase-amplitude couplings were decreased more for IC than NC and might be due to a resetting mechanism. The complex patterns in phase-phase coupling between theta and beta/gamma suggest that the contribution of these oscillations to visual binding and stimulus processing are not as straightforward as conventionally held. Causality analyses further suggest that theta oscillations drive beta/gamma oscillations (see also Schroeder and Lakatos, 2009). The present findings highlight the need for applying more sophisticated signal analyses in order to establish a fuller understanding of the functional role of neural oscillations.