Adaptive filtering methods for identifying cross-frequency couplings in human EEG.

Oscillations have been increasingly recognized as a core property of neural responses that contribute to spontaneous, induced, and evoked activities within and between individual neurons and neural ensembles. They are considered as a prominent mechanism for information processing within and communication between brain areas. More recently, it has been proposed that interactions between periodic components at different frequencies, known as cross-frequency couplings, may support the integration of neuronal oscillations at different temporal and spatial scales. The present study details methods based on an adaptive frequency tracking approach that improve the quantification and statistical analysis of oscillatory components and cross-frequency couplings. This approach allows for time-varying instantaneous frequency, which is particularly important when measuring phase interactions between components. We compared this adaptive approach to traditional band-pass filters in their measurement of phase-amplitude and phase-phase cross-frequency couplings. Evaluations were performed with synthetic signals and EEG data recorded from healthy humans performing an illusory contour discrimination task. First, the synthetic signals in conjunction with Monte Carlo simulations highlighted two desirable features of the proposed algorithm vs. classical filter-bank approaches: resilience to broad-band noise and oscillatory interference. Second, the analyses with real EEG signals revealed statistically more robust effects (i.e. improved sensitivity) when using an adaptive frequency tracking framework, particularly when identifying phase-amplitude couplings. This was further confirmed after generating surrogate signals from the real EEG data. Adaptive frequency tracking appears to improve the measurements of cross-frequency couplings through precise extraction of neuronal oscillations.

Amino acid coevolution reveals three-dimensional structure and functional domains of insect odorant receptors.

Insect odorant receptors (ORs) comprise an enormous protein family that translates environmental chemical signals into neuronal electrical activity. These heptahelical receptors are proposed to function as ligand-gated ion channels and/or to act metabotropically as G protein-coupled receptors (GPCRs). Resolving their signalling mechanism has been hampered by the lack of tertiary structural information and primary sequence similarity to other proteins. We use amino acid evolutionary covariation across these ORs to define restraints on structural proximity of residue pairs, which permit de novo generation of three-dimensional models. The validity of our analysis is supported by the location of functionally important residues in highly constrained regions of the protein. Importantly, insect OR models exhibit a distinct transmembrane domain packing arrangement to that of canonical GPCRs, establishing the structural unrelatedness of these receptor families. The evolutionary couplings and models predict odour binding and ion conduction domains, and provide a template for rationale structure-activity dissection.

Anti-tumour effect of monoclonal antibodies against selected antigens expressed by B-cells in Chronic Lymphocytic Leukaemia : strategies to enhance the efficacy of immunotherapy

Résumé : Les anticorps monoclonaux ont une place de plus en plus prépondérante dans le traitement des lymphomes et leucémies. Dans cette étude, trois anticorps monoclonaux murins, dirigés contre les antigènes CDS, CD71 et HLA-DR exprimés à la surface des cellules de leucémies lymphoïdes chroniques (LLC), ont été évalués. In vitro, les anticorps radiomarqués ont montrés des bonnes liaisons spécifiques sur les différentes cellules cibles. L'anti-CD71 inhibait la prolifération de la plupart des lignées cellulaires testées avec une accumulation des cellules en phase S précoce du cycle cellulaire. L'anti-HLA-DR inhibait aussi la prolifération des lignées leucémique JOK1-5.3 et lymphoïde Daudi. Cette inhibition était associée à une agrégation des cellules. Aucune induction d'apoptose n'a pu être clairement observée avec ces anticorps. L'anti-CD5 n'a montré aucun effet d'inhibition de croissance in vitro. In vivo, l'injection des anticorps individuellement augmentait significativement la survie médiane de souris SCID greffées avec des cellules JOK1-5.3 en i.p. De plus, l'anticorps antiCD5 combiné à l'anti-HLA-DR ou l'anti-CD71, sous certaines conditions, inhibait complètement le développement tumoral dans la quasi totalité des souris traitées avec une augmentation significative de l'efficacité comparée aux anticorps seuls. L'augmentation de l'efficacité thérapeutique des anticorps monoclonaux par les cytokines, dont l'IL-2, a déjà été montrée dans la littérature. Au regard du meilleur comportement de l'IL-2 sous la forme complexée à un anticorps anti-IL-2, nous avons évalué l'efficacité de l'IL-2/anti-IL-2 seul ou combinés au rituximab chez différents modèles tumoraux s.c. (BL60.2, Daudi, Ramos) ou i.p. (JOK15.3) de souris SCID. Le complexe IL-2/anti-IL-2 a montré un effet anti-tumoral dans les souris greffées avec BL60.2 et Daudi. Le traitement IL-2/anti-IL-2 combiné au rituximab a montré une efficacité accrue chez des souris avec BL60.2 par rapport au rituximab seul. En revanche, nous n'avons pas observé de différence avec IL-2/anti-IL-2 seul.Aussi, nous avons évalué l'utilisation de l'agent couplant tri-fonctionnel TMEA pour produire des anticorps bispecifiques. Les expériences préliminaires avec les anticorps rituximab et herceptine, ont mis en évidence sur gel SDS-Page la formation de dimers (~100kDa) et de trimers (~150kDa). Les anticorps bispecifiques sont composés d'un fragment Fab' d'une spécificité et de un ou deux fragments Fab' de l'autre spécificité permettant de moduler la capacité de liaison. Nous avons enfin montré qu'une construction anti-CD5/anti-CD20 était capable de se lier indépendamment ou simultanément à ses antigènes cibles. En conclusion, ce travail a montré l'efficacité thérapeutique des trois anticorps monoclonaux étudiés dans un model de LLC in vivo, et plus particulièrement l'intérêt de certaines combinaisons. D'autre part, nous avons montré l'efficacité anti-tumorale du complexe IL-2/anti-IL-2 in vivo. Des études futures devront permettre de définir un régime favorable pour augmenter l'efficacité de la thérapie avec les anticorps monoclonaux. Enfin, nous avons montré la faisabilité d'utiliser l'agent couplant TMEA pour produire des anticorps bispécifiques fonctionnels.Abstract : Monoclonal antibody (mAb) therapy has become an integral part in different treatments of lymphomas and leukaemias. In this study, we describe three murine mAbs directed against the CD5, CD71 and HLA-DR antigens expressed on chronic lymphocytic leukaemia cells (CLL). In vitro, radiolabeled purified mAbs showed good specific binding on live target cells. Anti-CD71 mAb inhibited proliferation of most cell lines with an accumulation of responding cells in early S-phase of the cell cycle, but without induction of apoptosis. Anti-HLA-DR mAb showed proliferation inhibition of leukaemia JOK1-5.3 and lymphoid Daudi cells, associated with cell aggregation, but again no specific sign of apoptosis was observed. Anti-CD5 mAb did not show any growth inhibitory effect in vitro. In vivo, in a model of SCID mice grafted i.p. with JOK1-5.3 cells, injection of individual mAbs induced significant prolongation of median survival, up to complete inhibition of tumour growth in some mice. Antibody combination of anti-CD5 with anti-HLA-DR or anti-CD71, evaluated in an early treatment, completely inhibited tumour growth in most mice, with a significant efficacy enhancement as compared to mAb used as single agents. Previous reports described the improved efficacy of mAb therapy when combined with cytokines such as IL-2. Relying further on the improved efficacy of IL-2 when administered as an immune complex with anti-IL-2 mAb, we evaluated the anti-tumour effect of the IL-2/anti-IL-2 complex alone or combined with rituximab in subcutaneous (BL60.2, Daudi, Ramos) or i.p. (JOK1-5.3) tumour models in SCID mice. The IL-2/anti-IL-2 complex demonstrated an anti-tumour effect in BL60.2 and Daudi grafted SCID mice. Combination of IL-2/anti-IL-2 treatment with rituximab showed increased efficacy as compared to rituximab alone in BL60.2 grafted mice. However, no difference was observed with IL-2/anti-IL-2 complex alone in these experiments. Finally, we evaluated the feasibility of producing bispecific antibodies (bsAbs) using a trifunctional coupling agent, called TMEA. In preliminary experiments coupling rituximab with herceptine Fab' fragments we obtained the formation of dimers (~100kDa) and trimers (~150kDa) as observed on SDS-Page gel. This method allowed us to produce bsAb with one Fab' fragments of one specificity and one or two Fab' fragments of the second specificity. An anti-CD5/anti-CD20 bsAb was shown to bind targeted antigen either independently or simultaneously. In conclusion, these data show that the three mAbs were all able to induce significant growth inhibition of the JOK1-5.3 cell line in vivo, and efficacy was enhanced when used in combination. IL2/anti-IL-2 complex displayed anti-tumour efficacy in vivo. Further evaluation is necessary to define the most favourable combination to improve mAb therapy. BsAb were produced using the tri-functional agent allowing antibody fragments with relatively good binding. The poor yield obtained with such chemical couplings limited the use of these constructs in preclinical experiments.

Resting-brain functional connectivity predicted by analytic measures of network communication.

The complex relationship between structural and functional connectivity, as measured by noninvasive imaging of the human brain, poses many unresolved challenges and open questions. Here, we apply analytic measures of network communication to the structural connectivity of the human brain and explore the capacity of these measures to predict resting-state functional connectivity across three independently acquired datasets. We focus on the layout of shortest paths across the network and on two communication measures-search information and path transitivity-which account for how these paths are embedded in the rest of the network. Search information is an existing measure of information needed to access or trace shortest paths; we introduce path transitivity to measure the density of local detours along the shortest path. We find that both search information and path transitivity predict the strength of functional connectivity among both connected and unconnected node pairs. They do so at levels that match or significantly exceed path length measures, Euclidean distance, as well as computational models of neural dynamics. This capacity suggests that dynamic couplings due to interactions among neural elements in brain networks are substantially influenced by the broader network context adjacent to the shortest communication pathways.

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.

Stimulus statistics shape oscillations in nonlinear recurrent neural networks.

Rhythmic activity plays a central role in neural computations and brain functions ranging from homeostasis to attention, as well as in neurological and neuropsychiatric disorders. Despite this pervasiveness, little is known about the mechanisms whereby the frequency and power of oscillatory activity are modulated, and how they reflect the inputs received by neurons. Numerous studies have reported input-dependent fluctuations in peak frequency and power (as well as couplings across these features). However, it remains unresolved what mediates these spectral shifts among neural populations. Extending previous findings regarding stochastic nonlinear systems and experimental observations, we provide analytical insights regarding oscillatory responses of neural populations to stimulation from either endogenous or exogenous origins. Using a deceptively simple yet sparse and randomly connected network of neurons, we show how spiking inputs can reliably modulate the peak frequency and power expressed by synchronous neural populations without any changes in circuitry. Our results reveal that a generic, non-nonlinear and input-induced mechanism can robustly mediate these spectral fluctuations, and thus provide a framework in which inputs to the neurons bidirectionally regulate both the frequency and power expressed by synchronous populations. Theoretical and computational analysis of the ensuing spectral fluctuations was found to reflect the underlying dynamics of the input stimuli driving the neurons. Our results provide insights regarding a generic mechanism supporting spectral transitions observed across cortical networks and spanning multiple frequency bands.

A critical analysis of three psychological research programs of doping behaviour

Objectives To consider the various specific substances-taking activities in sport an examination of three psychological models of doping behaviour utilised by researchers is presented in order to evaluate their real and potential impact, and to improve the relevance and efficiency of anti-doping campaigns. Design Adopting the notion of a "research program" (Lakatos, 1978) from the philosophy of science, a range of studies into the psychology of doping behaviour are classified and critically analysed. Method Theoretical and practical parameters of three research programs are critically evaluated (i) cognitive; (ii) drive; and (iii) situated-dynamic. Results The analysis reveals the diversity of theoretical commitments of the research programs and their practical consequences. The «cognitive program» assumes that athletes are accountable for their acts that reflect the endeavour to attain sporting and non-sporting goals. Attitudes, knowledge and rational decisions are understood to be the basis of doping behaviour. The «drive program» characterises the variety of traces and consequences on psychological and somatic states coming from athlete's experience with sport. Doping behaviour here is conceived of as a solution to reduce unconscious psychological and somatic distress. The «situated-dynamic program» considers a broader context of athletes' doping activity and its evolution during a sport career. Doping is considered as emergent and self-organized behaviour, grounded on temporally critical couplings between athletes' actions and situations and the specific dynamics of their development during the sporting life course. Conclusions These hypothetical, theoretical and methodological considerations offer a more nuanced understanding of doping behaviours, making an effective contribution to anti-doping education and research by enabling researchers and policy personnel to become more critically reflective about their explicit and implicit assumptions regarding models of explanations for doping behaviour.

Proton NMR of (15)N-choline metabolites enhanced by dynamic nuclear polarization.

Chemical shifts of protons can report on metabolic transformations such as the conversion of choline to phosphocholine. To follow such processes in vivo, magnetization can be enhanced by dynamic nuclear polarization (DNP). We have hyperpolarized in this manner nitrogen-15 spins in (15)N-labeled choline up to 3.3% by irradiating the 94 GHz electron spin resonance of admixed TEMPO nitroxide radicals in a magnetic field of 3.35 T during ca. 3 h at 1.2 K. The sample was subsequently transferred to a high-resolution magnet, and the enhanced polarization was converted from (15)N to methyl- and methylene protons, using the small (2,3)J((1)H,(15)N) couplings in choline. The room-temperature lifetime of nitrogen polarization in choline, T(1)((15)N) approximately 200 s, could be considerably increased by partial deuteration of the molecule. This procedure enables studies of choline metabolites in vitro and in vivo using DNP-enhanced proton NMR.

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.

Predicting human resting-state functional connectivity from structural connectivity

In the cerebral cortex, the activity levels of neuronal populations are continuously fluctuating. When neuronal activity, as measured using functional MRI (fMRI), is temporally coherent across 2 populations, those populations are said to be functionally connected. Functional connectivity has previously been shown to correlate with structural (anatomical) connectivity patterns at an aggregate level. In the present study we investigate, with the aid of computational modeling, whether systems-level properties of functional networks-including their spatial statistics and their persistence across time-can be accounted for by properties of the underlying anatomical network. We measured resting state functional connectivity (using fMRI) and structural connectivity (using diffusion spectrum imaging tractography) in the same individuals at high resolution. Structural connectivity then provided the couplings for a model of macroscopic cortical dynamics. In both model and data, we observed (i) that strong functional connections commonly exist between regions with no direct structural connection, rendering the inference of structural connectivity from functional connectivity impractical; (ii) that indirect connections and interregional distance accounted for some of the variance in functional connectivity that was unexplained by direct structural connectivity; and (iii) that resting-state functional connectivity exhibits variability within and across both scanning sessions and model runs. These empirical and modeling results demonstrate that although resting state functional connectivity is variable and is frequently present between regions without direct structural linkage, its strength, persistence, and spatial statistics are nevertheless constrained by the large-scale anatomical structure of the human cerebral cortex.