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.

Magnetoencephalography demonstrates multiple asynchronous generators during human sleep spindles.

Sleep spindles are approximately 1 s bursts of 10-16 Hz activity that occur during stage 2 sleep. Spindles are highly synchronous across the cortex and thalamus in animals, and across the scalp in humans, implying correspondingly widespread and synchronized cortical generators. However, prior studies have noted occasional dissociations of the magnetoencephalogram (MEG) from the EEG during spindles, although detailed studies of this phenomenon have been lacking. We systematically compared high-density MEG and EEG recordings during naturally occurring spindles in healthy humans. As expected, EEG was highly coherent across the scalp, with consistent topography across spindles. In contrast, the simultaneously recorded MEG was not synchronous, but varied strongly in amplitude and phase across locations and spindles. Overall, average coherence between pairs of EEG sensors was approximately 0.7, whereas MEG coherence was approximately 0.3 during spindles. Whereas 2 principle components explained approximately 50% of EEG spindle variance, >15 were required for MEG. Each PCA component for MEG typically involved several widely distributed locations, which were relatively coherent with each other. These results show that, in contrast to current models based on animal experiments, multiple asynchronous neural generators are active during normal human sleep spindles and are visible to MEG. It is possible that these multiple sources may overlap sufficiently in different EEG sensors to appear synchronous. Alternatively, EEG recordings may reflect diffusely distributed synchronous generators that are less visible to MEG. An intriguing possibility is that MEG preferentially records from the focal core thalamocortical system during spindles, and EEG from the distributed matrix system.

The role of PLK-1 in asynchronous cell division of two-cell stage "C. elegans" embryos

Summary Acquisition of lineage-specific cell cycle duration is an important feature of metazoan development. In Caenorhabditis a/egans, differences in cell cycle duration are already apparent in two-cell stage embryos, when the larger anterior blastomere AB divides before the smaller posterior blastomere P1. This time difference is under the control of anterior-posterior (A-P) polarity cues set by the PAR proteins. The mechanism by which these cues regulate the cell cycle machinery differentially in AB and P1 are incompletely understood. Previous work established that retardation of P1 cell division is due in part to preferential activation of an ATL1/CHK-1 dependent checkpoint in P1 but how the remaining time difference is controlled was not known at the onset of my work. The principal line of work in this thesis established that differential timing relies also on a mechanism that promotes mitosis onset preferentially in AB. The polo-like kinase PLK-1, a positive regulator of mitotic entry, is distributed in an asymmetric manner in two-cell stage embryos, with more protein present in AB than in P1. We find that PLK-1 asymmetry is regulated by anterior-posterior (A-P) polarity cues through preferential protein retention in the embryo anterior. Importantly, mild inactivation of plk-1 by RNAi delays entry into mitosis in P1 but not in AB, in a manner that is independent of ATL-1/CHK-1. Together, these findings favor a model in which differential timing of mitotic entry in C. elegans embryos relies on two complementary mechanisms: ATL-1/CHK-1 dependent preferential retardation in P1 and PLK-1 dependent preferential promotion in AB, which together couple polarity cues and cell cycle progression during early development. Besides analyzing PLK-1 asymmetry and its role in differential timing of two-cells stage embryos, we also characterized t2190, a mutant that exhibits reduced differential timing between AB and P1. We found this mutant to be a new allele of par-1. Additionally, we analyzed the role of NMY-2 in regulating the asynchrony of two-cell stage embryos, which may be uncoupled from its role in A-P polarity establishment and carried out a preliminary analysis of the mechanism underlying CDC-25 asymmetry between AB and P,. Overall, our works bring new insights into the mechanism controlling cell cycle progression in early C. elegans embryos. As most of the players important in C. elegans are conserved in other organisms, analogous mechanisms may be utilized in polarized cells of other species. Résumé Au cours du développement, les processus de division cellulaire sont régulés dans l'espace et le temps afin d'aboutir à la formation d'un organisme fonctionnel. Chez les Métazoaires, l'un des mécanismes de contrôle s'effectue au niveau de la durée du cycle cellulaire, celle-ci étant specifiée selon la lignée cellulaire. L'embryon du nématode Caenorhabditis elegans apparaît comme un excellent modèle d'étude de la régulation temporelle du cycle cellulaire. En effet, suite à la première division du zygote, l'embryon est alors composé de deux cellules de taille et d'identité différentes, appelées blastomères AB et P1. Ces deux cellules vont ensuite se diviser de manière asynchrone, le grand blastomère antérieur AB se divisant plus rapidement que le petit blastomère postérieur P1. Cette asynchronie de division est sous le contrôle des protéines PAR qui sont impliquées dans l'établissement de l'axe antéro-postérieur de l'embryon. A ce jour, les mécanismes moléculaires gouvernant ce processus d'asynchronie ne sont que partiellement compris. Des études menées précédemment ont établit que le retard de division observé dans le petit blastomère postérieur P1 était dû, en partie, à l'activation préférentielle dans cette cellule de ATL-1/CHK-1, protéines contrôlant la réponse à des erreurs dans le processus de réplication de l'ADN. L'analyse des autres mécanismes responsables de la différence temporelle d'entrée en mitose des deux cellules a été entreprise au cours de cette thèse. Nous avons considéré la possibilité que l'asynchronie de division était du à l'entrée préférentielle en mitose du grand blastomère AB. Nous avons établi que la protéine kinase PLK-1 (polo-like kinase 1), impliquée dans la régulation positive de la mitose, était distribuée de manière asymétrique dans l'embryon deux cellules. PLK-1 est en effet enrichi dans le blastomère AB. Cette localisation asymétrique de PLK-1 est sous le contrôle des protéines PAR et semble établie via une rétention de PLK-1 dans la cellule AB. Par ailleurs, nous avons démontré que l'inactivation partielle de plk-7 par interférence à ARN (RNAi) conduit à un délai de l'entrée en mitose de la cellule P1 spécifiquement, indépendamment des protéines régulatrices ATL-1/CHK-1. En conclusion, nous proposons un modèle de régulation temporelle de l'entrée en mitose dans l'embryon deux cellules de C. elegans basé sur deux mécanismes complémentaires. Le premier implique l'activation préférentielle des protéines ATL-1/CHK-1, et conduit à un retard d'entrée en mitose spécifiquement dans la cellule P1. Le second est basé sur la localisation asymétrique de la protéine kinase PLK-1 dans la cellule AB et induit une entrée précoce en mitose de cette cellule. Par ailleurs, nous avons étudié un mutant appelé t2190 qui réduit la différence temporelle d'entrée en mitose entre les cellules AB et P1. Nous avons démontré que ce mutant correspondait à un nouvel allèle du Bene par-1. De plus, nous avons analysé le rôle de NMY-2, une protéine myosine qui agit comme moteur moléculaire sur les filaments d'active; dans la régulation de l'asynchronie de division des blastomères AB et P1, indépendamment de sa fonction dans l'établissement de l'axe antéro-postérieur. Par ailleurs, nous avons commencé l'étude du mécanisme moléculaire régulant la localisation asymétrique entre les cellules AB et P1 de la protéine phosphatase CDC25, qui est également un important régulateur de l'entrée en mitose. En conclusion, ce travail de thèse a permis une meilleure compréhension des mécanismes gouvernant la progression du cycle cellulaire dans l'embryon précoce de C. elegans. Etant donné que la plupart des protéines impliquées dans ces processus sont conservées chez d'autres organismes multicellulaires, il apparaît probable que les mécanismes moléculaires révélés dans cette étude soit aussi utilisés chez ceux-ci.

Pathologic ventricular hypertrophy in the offspring of diabetic mothers : a retrospective study

L'hypertrophie ventriculaire pathologique chez les nouveau-nés des mères diabétiques une étude rétrospective RESUME Objectif L'incidence du diabète chez les femmes enceintes ne cesse de croître, de même que les complications chez leurs nouveau-nés. C'est pourquoi, nous avons étudié la population de mères diabétiques suivies dans notre établissement entre les années 2003-2005 dans le but d'analyser spécifiquement le problème d'hypertrophie ventriculaire pathologique (HVP) chez les nouveau-nés de cette population. Méthode et résultats Dans notre étude rétrospective comprenant 87 grossesses de femmes diabétiques (92 nouveau-nés), 16 présentaient un diabète de type 1, 17 de type 2 et 54 ont développé un diabète gestationnel (DG). Le médian des hémoglobines glycquées (HbAlc) pour cette population est de 5.8% (5.3-6.5) : 17 avaient une HbAlc au-dessus de la norme, dont 2 souffrant d'une cardiomyopathie congénitale (CMC) et six d'une HVP. Un total de 75 nouveaux-nés étaient normaux, cinq avaient une CMC et 12 une HVP (1/12 décédé post-natalement, 1/12 mort-né, 2/12 nécessitant un accouchement prématuré, 8/12 normaux). Les 16 mères avec un diabète de type 1 accouchèrent de trois nouveau-nés avec une CMC et de 50% avec une HVP, comprenant un enfant décédé et un prématuré né par césarienne à cause d'une HVP. Dans le groupe des 17 nouveau-nés issus d'une mère connue pour un diabète de type 2, un cas présentait une CMC et 25% des cas une HVP. Parmi les 54 grossesses avec un DG, on dénombre un cas de CMC et un cas de HVP. Conclusion Les grossesses de mères souffrant d'un diabète de type 1 et de type 2 comportent toutes deux un risque augmenté de développement d'une HVP comparées à celles de mères ayant développé un diabète gestationnel. Les contrôles glycémiques sont insuffisants pour éviter la survenue d'une HVP. Comme aucun autre paramètre prédictif n'a pu été défini jusqu'alors, nous concluons qu'un suivi échographique rapproché de ces grossesses peut prévenir des complications périnatales sévères.

The L-Type Ca+ and KATP channels may contribute to pacing-induced protection against anoxia-reoxygenation in the embryonic heart model.

L-Type Ca(2+) and K(ATP) Channels in Pacing-Induced Cardioprotection. AIMS: The L-type Ca(2+) channel, the sarcolemmal (sarcK(ATP)), and mitochondrial K(ATP) (mitoK(ATP)) channels are involved in myocardial preconditioning. We aimed at determining to what extent these channels can also participate in pacing-induced cardioprotection. METHODS: Hearts of 4-day-old chick embryos were paced in ovo during 12 hour using asynchronous intermittent ventricular stimulation at 110% of the intrinsic rate. Sham operated and paced hearts were then submitted in vitro to anoxia (30 minutes) and reoxygenation (60 minutes). These hearts were exposed to L-type Ca(2+) channel agonist Bay-K-8644 (BAY-K) or blocker verapamil, nonselective K(ATP) channel antagonist glibenclamide (GLIB), mitoK(ATP) channel agonist diazoxide (DIAZO), or antagonist 5-hydroxydecanoate. Electrocardiogram, electromechanical delay (EMD) reflecting excitation-contraction (E-C) coupling, and contractility were determined. RESULTS: Under normoxia, heart rate, QT duration, conduction, EMD, and ventricular shortening were similar in sham and paced hearts. During reoxygenation, arrhythmias ceased earlier and ventricular EMD recovered faster in paced hearts than in sham hearts. In sham hearts, BAY-K (but not verapamil), DIAZO (but not 5-hydroxydecanoate) or GLIB accelerated recovery of ventricular EMD, reproducing the pacing-induced protection. By contrast, none of these agents further ameliorated recovery of the paced hearts. CONCLUSION: The protective effect of chronic asynchronous pacing at near physiological rate on ventricular E-C coupling appears to be associated with subtle activation of L-type Ca(2+) channel, inhibition of sarcK(ATP) channel, and/or opening of mitoK(ATP) channel.

Dynamical modeling and analysis of large cellular regulatory networks.

The dynamical analysis of large biological regulatory networks requires the development of scalable methods for mathematical modeling. Following the approach initially introduced by Thomas, we formalize the interactions between the components of a network in terms of discrete variables, functions, and parameters. Model simulations result in directed graphs, called state transition graphs. We are particularly interested in reachability properties and asymptotic behaviors, which correspond to terminal strongly connected components (or "attractors") in the state transition graph. A well-known problem is the exponential increase of the size of state transition graphs with the number of network components, in particular when using the biologically realistic asynchronous updating assumption. To address this problem, we have developed several complementary methods enabling the analysis of the behavior of large and complex logical models: (i) the definition of transition priority classes to simplify the dynamics; (ii) a model reduction method preserving essential dynamical properties, (iii) a novel algorithm to compact state transition graphs and directly generate compressed representations, emphasizing relevant transient and asymptotic dynamical properties. The power of an approach combining these different methods is demonstrated by applying them to a recent multilevel logical model for the network controlling CD4+ T helper cell response to antigen presentation and to a dozen cytokines. This model accounts for the differentiation of canonical Th1 and Th2 lymphocytes, as well as of inflammatory Th17 and regulatory T cells, along with many hybrid subtypes. All these methods have been implemented into the software GINsim, which enables the definition, the analysis, and the simulation of logical regulatory graphs.

Imaging EEG synchronization in schizophrenia patients with S-estimator

Recently a new measure of the cooperative behavior of simultaneous time series was introduced (Carmeli et al. NeuroImage 2005). This measure called S-estimator is defined from the embedding dimension in a state space. S-estimator quantifies the amount of synchronization within a data set by comparing the actual dimensionality of the set with the expected full dimensionality of the asynchronous set. It has the advantage of being a multivariate measure over traditionally used in systems neuroscience bivariate measures of synchronization. Multivariate measures of synchronization are of particular interest for applications in the field of modern multichannel EEG research, since they easily allow mapping of local and/or regional synchronization and are compatible with other imaging techniques. We applied Sestimator to the analysis of EEG synchronization in schizophrenia patients vs. matched controls. The whole-head mapping with S-estimator revealed a specific pattern of local synchronization in schizophrenia patients. The differences in the landscape of synchronization included decreased local synchronization in the territories over occipital and midline areas and increased synchronization over temporal areas. In frontal areas, the S-estimator revealed a tendency for an asymmetry: decreased S-values over the left hemisphere were adjacent to increased values over the right hemisphere. Separate calculations showed reproducibility of this pattern across the main EEG frequency bands. The maintenance of the same synchronization landscape across EEG frequencies probably implies the structural changes in the cortical circuitry of schizophrenia patients. These changes are regionally specific and suggest that schizophrenia is a misconnectivity rather than hypo- or hyper-connectivity disorder.

Ambulatory measurement of ankle kinetics for clinical applications.

This study aimed to design and validate the measurement of ankle kinetics (force, moment, and power) during consecutive gait cycles and in the field using an ambulatory system. An ambulatory system consisting of plantar pressure insole and inertial sensors (3D gyroscopes and 3D accelerometers) on foot and shank was used. To test this system, 12 patients and 10 healthy elderly subjects wore shoes embedding this system and walked many times across a gait lab including a force-plate surrounded by seven cameras considered as the reference system. Then, the participants walked two 50-meter trials where only the ambulatory system was used. Ankle force components and sagittal moment of ankle measured by ambulatory system showed correlation coefficient (R) and normalized RMS error (NRMSE) of more than 0.94 and less than 13% in comparison with the references system for both patients and healthy subjects. Transverse moment of ankle and ankle power showed R>0.85 and NRMSE<23%. These parameters also showed high repeatability (CMC>0.7). In contrast, the ankle coronal moment of ankle demonstrated high error and lower repeatability. Except for ankle coronal moment, the kinetic features obtained by the ambulatory system could distinguish the patients with ankle osteoarthritis from healthy subjects when measured in 50-meter trials. The proposed ambulatory system can be easily accessible in most clinics and could assess main ankle kinetics quantities with acceptable error and repeatability for clinical evaluations. This system is therefore suggested for field measurement in clinical applications.

How local excitation-inhibition ratio impacts the whole brain dynamics.

The spontaneous activity of the brain shows different features at different scales. On one hand, neuroimaging studies show that long-range correlations are highly structured in spatiotemporal patterns, known as resting-state networks, on the other hand, neurophysiological reports show that short-range correlations between neighboring neurons are low, despite a large amount of shared presynaptic inputs. Different dynamical mechanisms of local decorrelation have been proposed, among which is feedback inhibition. Here, we investigated the effect of locally regulating the feedback inhibition on the global dynamics of a large-scale brain model, in which the long-range connections are given by diffusion imaging data of human subjects. We used simulations and analytical methods to show that locally constraining the feedback inhibition to compensate for the excess of long-range excitatory connectivity, to preserve the asynchronous state, crucially changes the characteristics of the emergent resting and evoked activity. First, it significantly improves the model's prediction of the empirical human functional connectivity. Second, relaxing this constraint leads to an unrealistic network evoked activity, with systematic coactivation of cortical areas which are components of the default-mode network, whereas regulation of feedback inhibition prevents this. Finally, information theoretic analysis shows that regulation of the local feedback inhibition increases both the entropy and the Fisher information of the network evoked responses. Hence, it enhances the information capacity and the discrimination accuracy of the global network. In conclusion, the local excitation-inhibition ratio impacts the structure of the spontaneous activity and the information transmission at the large-scale brain level.

Vaccination with a Melan-A peptide selects an oligoclonal T cell population with increased functional avidity and tumor reactivity.

Both the underlying molecular mechanisms and the kinetics of TCR repertoire selection following vaccination against tumor Ags in humans have remained largely unexplored. To gain insight into these questions, we performed a functional and structural longitudinal analysis of the TCR of circulating CD8(+) T cells specific for the HLA-A2-restricted immunodominant epitope from the melanocyte differentiation Ag Melan-A in a melanoma patient who developed a vigorous and sustained Ag-specific T cell response following vaccination with the corresponding synthetic peptide. We observed an increase in functional avidity of Ag recognition and in tumor reactivity in the postimmune Melan-A-specific populations as compared with the preimmune blood sample. Improved Ag recognition correlated with an increase in the t(1/2) of peptide/MHC interaction with the TCR as assessed by kinetic analysis of A2/Melan-A peptide multimer staining decay. Ex vivo analysis of the clonal composition of Melan-A-specific CD8(+) T cells at different time points during vaccination revealed that the response was the result of asynchronous expansion of several distinct T cell clones. Some of these T cell clones were also identified at a metastatic tumor site. Collectively, these data show that tumor peptide-driven immune stimulation leads to the selection of high-avidity T cell clones of increased tumor reactivity that independently evolve within oligoclonal populations.

Functional calibration procedure for 3D knee joint angle description using inertial sensors.

Measurement of three-dimensional (3D) knee joint angle outside a laboratory is of benefit in clinical examination and therapeutic treatment comparison. Although several motion capture devices exist, there is a need for an ambulatory system that could be used in routine practice. Up-to-date, inertial measurement units (IMUs) have proven to be suitable for unconstrained measurement of knee joint differential orientation. Nevertheless, this differential orientation should be converted into three reliable and clinically interpretable angles. Thus, the aim of this study was to propose a new calibration procedure adapted for the joint coordinate system (JCS), which required only IMUs data. The repeatability of the calibration procedure, as well as the errors in the measurement of 3D knee angle during gait in comparison to a reference system were assessed on eight healthy subjects. The new procedure relying on active and passive movements reported a high repeatability of the mean values (offset<1 degrees) and angular patterns (SD<0.3 degrees and CMC>0.9). In comparison to the reference system, this functional procedure showed high precision (SD<2 degrees and CC>0.75) and moderate accuracy (between 4.0 degrees and 8.1 degrees) for the three knee angle. The combination of the inertial-based system with the functional calibration procedure proposed here resulted in a promising tool for the measurement of 3D knee joint angle. Moreover, this method could be adapted to measure other complex joint, such as ankle or elbow.

Safety, pharmacokinetics, and preliminary clinical activity of inotuzumab ozogamicin, a novel immunoconjugate for the treatment of B-cell non-Hodgkin's lymphoma: results of a phase I study.

PURPOSE Inotuzumab ozogamicin (CMC-544) is an antibody-targeted chemotherapy agent composed of a humanized anti-CD22 antibody conjugated to calicheamicin, a potent cytotoxic agent. This was a phase I study to determine the maximum-tolerated dose (MTD), safety, and preliminary efficacy of inotuzumab ozogamicin in an expanded MTD cohort of patients with relapsed or refractory CD22(+) B-cell non-Hodgkin's lymphoma (NHL). PATIENTS AND METHODS Inotuzumab ozogamicin was administered intravenously as a single agent once every 3 or 4 weeks at doses ranging from 0.4 to 2.4 mg/m(2). Outcomes included MTD, safety, pharmacokinetics, response, progression-free survival (PFS), and overall survival. Results Seventy-nine patients were enrolled. The MTD was determined to be 1.8 mg/m(2). Common adverse events at the MTD were thrombocytopenia (90%), asthenia (67%), and nausea and neutropenia (51% each). The objective response rate at the end of treatment was 39% for the 79 enrolled patients, 68% for all patients with follicular NHL treated at the MTD, and 15% for all patients with diffuse large B-cell lymphoma treated at the MTD. Median PFS was 317 days (approximately 10.4 months) and 49 days for patients with follicular NHL and diffuse large B-cell lymphoma, respectively. CONCLUSION Inotuzumab ozogamicin has demonstrated efficacy against CD22(+) B-cell NHL, with reversible thrombocytopenia as the main toxicity.

Pacing redistributes glycogen within the developing myocardium.

Electrical pacing at physiological rate induces myocardial remodeling associated with regional changes in workload, blood flow and oxygen consumption. However, to what extent energy-producing pathways are also modified within the paced heart remains to be investigated. Pacing could particularly affect glycogen metabolism since hypertrophy stimulates glycolysis and increased workload favors glucose over fat oxidation. In order to test this hypothesis, we used the embryonic chick heart model in which ventricular pacing rapidly resulted in thinning of the ventricle wall and thickening of the atrial wall. Hearts of stage 22HH chick embryos were submitted in ovo to asynchronous and intermittent ventricular pacing delivered at physiological rate during 24 h. The resulting alterations of glycogen content were determined in atrium, ventricle and conotruncus of paced and sham-operated hearts. Hemodynamic parameters of the paced and spontaneously beating hearts were derived from computerized image analysis of video recordings. With respect to sham, paced hearts showed a significant decrease in glycogen content (nmoles glucose units/microg protein; mean+/-S.D.) only in atrium (1.48+/-0.40 v 0.84+/-0.34, n=8) and conotruncus (0.75+/-0.28 v 0.42+/-0.23, n=8). Pacing decreased the end diastolic and stroke volumes by 34 and 44%, respectively. Thus, the rapid glycogen depletion in regions remote from the stimulation site appears to be associated with regional changes in workload and remodeling. These findings underscore the importance of the coupling mechanisms between metabolic pathways and myocardial remodeling in the ectopically paced heart.

Temporo-Parietal Junction encodes the Embodied Self: Neuro-Robotics, fMRI, and Lesion mapping.

Introduction: Neuroimaging of the self focused on high-level mechanisms such as language, memory or imagery of the self. Recent evidence suggests that low-level mechanisms of multisensory and sensorimotor integration may play a fundamental role in encoding self-location and the first-person perspective (Blanke and Metzinger, 2009). Neurological patients with out-of body experiences (OBE) suffer from abnormal self-location and the first-person perspective due to a damage in the temporo-parietal junction (Blanke et al., 2004). Although self-location and the first-person perspective can be studied experimentally (Lenggenhager et al., 2009), the neural underpinnings of self-location have yet to be investigated. To investigate the brain network involved in self-location and first-person perspective we used visuo-tactile multisensory conflict, magnetic resonance (MR)-compatible robotics, and fMRI in study 1, and lesion analysis in a sample of 9 patients with OBE due to focal brain damage in study 2. Methods: Twenty-two participants saw a video showing either a person's back or an empty room being stroked (visual stimuli) while the MR-compatible robotic device stroked their back (tactile stimulation). Direction and speed of the seen stroking could either correspond (synchronous) or not (asynchronous) to those of the seen stroking. Each run comprised the four conditions according to a 2x2 factorial design with Object (Body, No-Body) and Synchrony (Synchronous, Asynchronous) as main factors. Self-location was estimated using the mental ball dropping (MBD; Lenggenhager et al., 2009). After the fMRI session participants completed a 6-item adapted from the original questionnaire created by Botvinick and Cohen (1998) and based on questions and data obtained by Lenggenhager et al. (2007, 2009). They were also asked to complete a questionnaire to disclose the perspective they adopted during the illusion. Response times (RTs) for the MBD and fMRI data were analyzed with a 3-way mixed model ANOVA with the in-between factor Perspective (up, down) and the two with-in factors Object (body, no-body) and Stroking (synchronous, asynchronous). Quantitative lesion analysis was performed using MRIcron (Rorden et al., 2007). We compared the distributions of brain lesions confirmed by multimodality imaging (Knowlton, 2004) in patients with OBE with those showing complex visual hallucinations involving people or faces, but without any disturbance of self-location and first person perspective. Nine patients with OBE were investigated. The control group comprised 8 patients. Structural imaging data were available for normalization and co-registration in all the patients. Normalization of each patient's lesion into the common MNI (Montreal Neurological Institute) reference space permitted simple, voxel-wise, algebraic comparisons to be made. Results: Even if in the scanner all participants were lying on their back and were facing upwards, analysis of perspective showed that half of the participants had the impression to be looking down at the virtual human body below them, despite any cues about their body position (Down-group). The other participants had the impression to be looking up at the virtual body above them (Up-group). Analysis of Q3 ("How strong was the feeling that the body you saw was you?") indicated stronger self-identification with the virtual body during the synchronous stroking. RTs in the MBD task confirmed these subjective data (significant 3-way interaction between perspective, object and stroking). fMRI results showed eight cortical regions where the BOLD signal was significantly different during at least one of the conditions resulting from the combination of Object and Stroking, relative to baseline: right and left temporo-parietal junction, right EBA, left middle occipito-temporal gyrus, left postcentral gyrus, right medial parietal lobe, bilateral medial occipital lobe (Fig 1). The activation patterns in right and left temporo-parietal junction and right EBA reflected changes in self-location and perspective as revealed by statistical analysis that was performed on the percentage of BOLD change with respect to the baseline. Statistical lesion overlap comparison (using nonparametric voxel based lesion symptom mapping) with respect to the control group revealed the right temporo-parietal junction, centered at the angular gyrus (Talairach coordinates x = 54, y =-52, z = 26; p>0.05, FDR corrected). Conclusions: The present questionnaire and behavioural results show that - despite the noisy and constraining MR environment) our participants had predictable changes in self-location, self-identification, and first-person perspective when robotic tactile stroking was applied synchronously with the robotic visual stroking. fMRI data in healthy participants and lesion data in patients with abnormal self-location and first-person perspective jointly revealed that the temporo-parietal cortex especially in the right hemisphere encodes these conscious experiences. We argue that temporo-parietal activity reflects the experience of the conscious "I" as embodied and localized within bodily space.