|Home | About | Journals | Submit | Contact Us | Français|
The study of brain functioning is a major challenge in neuroscience fields as human brain has a dynamic and ever changing information processing. Case is worsened with conditions where brain undergoes major changes in so-called different conscious states. Even though the exact definition of consciousness is a hard one, there are certain conditions where the descriptions have reached a consensus. The sleep and the anesthesia are different conditions which are separable from each other and also from wakefulness. The aim of our group has been to tackle the issue of brain functioning with setting up similar research conditions for these three conscious states.
In order to achieve this goal we have designed an auditory stimulation battery with changing conditions to be recorded during a 40 channel EEG polygraph (Nuamps) session. The stimuli (modified mismatch, auditory evoked etc.) have been administered both in the operation room and the sleep lab via Embedded Interactive Stimulus Unit which was developed in our lab. The overall study has provided some results for three domains of consciousness. In order to be able to monitor the changes we have incorporated Bispectral Index Monitoring to both sleep and anesthesia conditions.
The first stage results have provided a basic understanding in these altered states such that auditory stimuli have been successfully processed in both light and deep sleep stages. The anesthesia provides a sudden change in brain responsiveness; therefore a dosage dependent anesthetic administration has proved to be useful. The auditory processing was exemplified targeting N1 wave, with a thorough analysis from spectrogram to sLORETA. The frequency components were observed to be shifting throughout the stages. The propofol administration and the deeper sleep stages both resulted in the decreasing of N1 component. The sLORETA revealed similar activity at BA7 in sleep (BIS 70) and target propofol concentration of 1.2 µg/mL.
The current study utilized similar stimulation and recording system and incorporated BIS dependent values to validate a common approach to sleep and anesthesia. Accordingly the brain has a complex behavior pattern, dynamically changing its responsiveness in accordance with stimulations and states.
The human brain represents one of the most complex biological systems in the world. When it comes to its computational power, the estimates are overwhelming and are not comparable to man-made systems . This complex organ also has a very dynamic nature, during which it shifts from one state to another, almost entirely changing its functional properties. These major states can be exemplified in the case of conscious states. Hence the consciousness can also be defined with this concept of altered state of brain cognitive functioning. A perfect definition of consciousness is still not available; however the existence of differentiations between some conscious states is out of debate. The sleep, anesthesia and wakefulness are three of such separable states with distinct features (Figure (Figure11).
The changes in the dynamic features of the brain can be obtained through its responsiveness. The conscious states are usually determined with the degree of lack of responses to outer world. Therefore a gradually changing response monitoring system can be useful in this type of approach. The auditory domain could be a proper stimulation modality, as it has been a common choice in neuroscience [2,3]. The auditory evoked potentials have been described with the acoustic features and stimulation patterns. Out of these sub features, the frequency, intensity and location are among the most well known properties [2,4]. Typically, tones (sinusoidal) or clicks are applied as the stimuli and the sensory/cognitive responses are recorded via electrodes attached to the scalp. The basic sensory features (bottom-up) are inherently accompanied with other features such as attention, memory (top-down) etc. Thus, the top-down and bottom-up properties all paint a complex landscape with very short latencies such as Auditory Brain-stem Responses (ABR), Middle Latency Responses (MLRs) to long latency wave forms . A recently revisited auditory stimulation pattern has been the dichotic one, with a certain scope on brain asymmetry (syllables, tones etc.) [6-8]. From tones and clicks to syllables, these simple and semi-complex stimulations form a plausible domain for assessment of brain responsiveness. The need to obtain different features has forced the experimental design to include acoustical properties, attention-free assessment with mismatch, attention-bound oddball, dichotic and diotic features. Therefore these different properties have become a bundle of blocks constituting Brain biophysics battery (BBB).
On the other hand, the conditions of performing ideal experiments may not be an easy task as the operation room and the sleep environment are far from perfect neuropsychological setups. One of the major problems in surgery room comes from being electromagnetically hostile environment [9,10], therefore effecting the outcome of electrophysiological recording. During sleep, the ever changing body position as well as depth of sleep becomes problematic in conducting sensory/cognitive auditory tests (i.e. presenting perfect symmetrical and ideal auditory stimulations). Thirdly the stimulations especially in sleep may cause some alterations in the state of responsiveness of the brain .
Recently, the state of consciousness and brain functioning during sleep and anesthesia have become a topic of interest [5,12-16]. Furthermore Tung and Medelson  reviewed the studies and issues related to both sleep and anesthesia addressing components from neuromodulators to nuclei. Having noted above mentioned points, the present paper will address issues related to brain responsiveness in anesthesiology and sleep. The presentation of auditory stimulations with different physical and design related properties would enable the acquisition of some of the dynamic and complex response patterns of the brain.
The sleep experiments were performed in the Sleep Dynamics Research Laboratory of the Biophysics Department. The participants slept for one night in the laboratory and therefore, the current study is based on their first night sleep recordings. Overnight sleep data were collected from 12 healthy volunteer individuals (mean age: 24.5, range: 18-32, 10 males). Except one, all participants were right handed. The exclusion criteria included a past and/or present history of any neurological, psychiatric or chronic disorders and use of any drugs that effect cognitive functions. The participants were advised not to take any caffeine on the night of the recording day.
Participants arrived at the laboratory approximately two hours before their normal sleeping time and recordings were realized approximately between 23.00-08.00 o’clock. Over the entire course of the recording session, researchers stayed in the laboratory in which the recording systems were present. The participants spent a whole night lying on the bed in the isolated room which limited electromagnetic interference. Furthermore, the walls of the laboratory have been enveloped with an acoustic material to isolate room from the external auditory noise. The sleeping room was dimly lighted and the subject was monitored.
Bispectral index (BIS) recording was performed using a bispectral index monitor (Aspect-A2000) with a sensor (BIS Quatro). Sensors were checked for signal quality (SQI), assuring impedance below 5 kOhms. In every five seconds BIS was recorded via the RS232 cable using a HyperTerminal protocol. The BIS SQI values had been planned to be rejected below 50. However obtained SQI values were over 50; thus no rejection was required. Sleep scoring was performed according to the Rechtschaffen and Kales (R&K)  criteria in 30-s time windows. For every R&K score there were six BIS values (a BIS value obtained in every five seconds). Sleep data were analyzed in five BIS clusters, which have BIS values 90, 80, 70, 60, 50 respectively. For each stimulus, sweeps with 500 ms prestim and 1500 ms poststim data segments were formed in these five BIS clusters. The sweeps which had amplitudes higher than ±100 µV in Electrooculography (EOG) channel were automatically rejected. The recordings were corrected based on horizontal axis (baseline corrected) and digitally filtered with the 0.5-48 Hz band pass filter (12 dB/oct and zero phase shift, Neuroscan 4.3). Henceforth, every single subject’s data were averaged and all of the subjects’ grand average was prepared.
The anesthesia related experiments were conducted in the surgery room. 12 subjects (9 females) participated with an age range of 25-46 (mean 37.75) years. These were patients with a similar surgery protocol undergoing for lumbar disc hernia. The ethical approval has been received from the local Ethics Committee. All patients have signed a consent form prior to procedures.
The propofol administration was controlled with a target-controlled infusion mechanism (Fresenius Vial Orchestra Base Primea, Le Grand Chemin, France) where precalculated brain dosage (effect-site) would be obtained. Accordingly, the concentration levels of 0 µg/mL to 1.6 µg/mL of Propofol were obtained. Besides the various monitoring devices, the bispectral index monitor was used. The BIS provides a dimensionless number from 0 to 100, which denotes to anesthesia depth. It utilizes electrophysiological (EEG) parameters including bispectral index. The wakefulness values range around >95 and with the propofol administration sudden drops in the index values occur. At around 1.0 µg/mL of Propofol administration, the activity would be in the range around 80. Commonly levels would be kept around 50-60 range to perform adequate surgery.
The brain biophysics battery is an amalgamation of auditory blocks with different features:
Auditory evoked potential (AEP): 1500 Hz, 80 dBSPL, 500 ms, 2–3s inter stimulus interval (ISI). Unless mentioned otherwise throughout this manuscript AEP was used as the main example of stimulus condition.
Auditory event related potential (AERP): Auditory oddball paradigm applied via headphones (binaural). Auditory stimuli were 500 ms, 80 dBSPL, sinusoidal (target: 1600 Hz, nontarget:1500 Hz) sounds. Targets occurred in 20-25% of the cases.
Dichotic linguistic (DL): 36 different combinations of /ba/, /da/, /ga/, /ka/, /pa/, /ta/ syllables were represented by headphones, one from each ear. ISI was varied randomly between 2.5 and 3 seconds.
Dichotic tone (DT): Two different tones applied via headphones binaurally. Tones were 291 Hz and 392 Hz.
Mismatch negativity (MMN): A modified version of optimized MMN paradigm developed by Pakarinen et al.  was used. Standard tone was 75 ms in duration, at an intensity of 70 dB, composed of three sinusoidal partials (523 Hz as fundamental frequency) and presented simultaneously via headphones (resulting in perception of the sound source as localized in the centre). There were four type of deviants (duration, intensity, frequency, location), each with three levels. Location deviance was achieved by delaying the tone monaurally, right or left with equal probability.
These abovementioned stimuli were given in varying blocks from the onset to the end of the experimental sessions.
A special stimulus unit developed in our lab Embedded Microcontroller Interactive Stimulus Unit (EMISU)  was used for the generation of stimuli and the recording of electrophysiological data both in sleep and anesthesia experiments. The EEG recording was achieved in the operation room by means of 40 channel EEG (Nuamps), headphones with noise cancellation feature (Creative HN-700, Republic of Singapore). MATLAB, EEGLAB, and EMISU software were used. The auditory stimuli were administered through noise cancelling head phones at a level of 85 dB SPL. The stimuli were consisted of tones 1500 Hz, 500 ms duration 30 ms rise/fall time, with randomly varying intervals ranging from 2.5 to 5 sec.
The continuous EEG data files (*.cnt) were epoched for the evoked condition or segmented for the spontaneous condition to equal data lengths. Each segment contained 601 data points (sampling rate 1000 Hz) which equaled to -200 to 400 ms for the epoched sections. Spontaneous segment would be selected from an artifact free section with no stimulation. In a section with stimulations, then the segment would be temporally separated by 2-3 seconds pre and post stimuli.
By means of manual observation the artifacts were removed. These artifact free sections were baseline corrected. The obtained sweeps were than averaged. Due to hostile conditions in the surgery room not all the channels were artifact free. For this reason 18 channels [21,22] with equal distribution across anterior/posterior and hemispheric locations were fixed across all subjects. These data segments were then transferred into standardized low resolution electromagnetic tomography (sLORETA) protocol. The data were transformed with proper electrode locations to obtain X, Y, Z coordinates, which were then used to create transformation matrix file in sLORETA application. “Global field power” (GFP) was used to estimate the target peaks. Therefore at a given sLORETA image the current density value (CDV) of maximal GFP peaks are presented.
sLORETA was applied to estimate the neural sources of event related scalp potentials [23,24] in the LORETA 2D potential distribution to 3D brain volume. Montreal Neurological Institute (MNI) was referred for standard solution space which was subdivided into a 3D grid consisting of 2394 volume elements with a regular cube size of 7x7x7 mm. An equivalent current dipole is positioned, for which a current density (µA/mm2) is computed on each of these voxels [23,24].
The spectrograms for the presentations of EEG data segments were obtained by means of a MATLAB routine . The outline of the method is given below:
Wavelet transform of a signal gives the time-scale information found by the Wavelet transform of a signal s(t) gives the time-scale information found by the following formula
where the mother wavelet function was chosen as Complex Morlet.
* denotes the complex conjugate, b is the shift parameter, the corresponding scale is given by dilation parameter a. Since most of the signal processing applications focus on the time-frequency analysis, it is essential to convert the scale information into frequency. The formula given below provides corresponding frequency Fa for a specified scale value a
with the parameters Fc denoting center frequency of the wavelet in Hz, Δ denoting the sampling period. For obtaining the spectrogram the Central (Cz) electrode data were used.
The Shapiro-Wilk test was used for determining distribution of the groups and deciding statistical methods. When the groups were in normal distribution paired t-test was used to compare the above mentioned two conditions. The correlation of BIS and propofol dosage was determined with nominal logistic regression.
The electrophysiological recordings, the auditory stimulations, and the BIS monitoring have been successfully performed in all subjects of the two experimental groups. Additionally in anesthesia group, drug infusion system was used. The conditions for the stimulations were kept similar.
All subjects have undergone the step-wise propofol administration. One subject’s data were not available for both BIS and perfusion control system, therefore it was dismissed. The subjects had been instructed to stay awake. Around 0.8-1.2 µg/mL (propofol) they have lost consciousness (LOC) at which point they ceased to respond to comments. The electrophysiological recordings were continued regardless of LOC till reaching level of 1.6 µg/mL. Therefore all subjects received a similar dosage regime of propofol as well as auditory stimulations.
The figure figure22 represents the distribution of BIS values across the propofol dosage.
Nominal logistic regression: propofol concentrations versus BIS Scores
In table table11 nominal logistic regression analysis has been applied to model the relationship between propofol concentrations and BIS values. Positive coefficient value for logit 1 indicates that as BIS score increases propofol concentration tends to be in propofol concentration 1.2 µg/mL compared to propofol concentration 1.6 µg/mL. Since all coefficients are positive, the conclusions are all same. p-value = 0.000 for all coefficients show that BIS score affects being in propofol concentrations 1.2, 0.8, 0.4 and 0 µg/mL compared to concentration 1.6 µg/mL. The results reveal that the odds increase when we compare lower concentrations with concentration 1.6 µg/mL. Nominal logistic regression has been chosen over ordinal logistic regression since effect of BIS are not same for each concentration of propofol which is obvious from the larger regression coefficients and odds ratios for logit 2, 3 and 4.
Accordingly, the increments of propofol resulted in decreasing of the BIS index (p<0.0001) (Table (Table2).2). The auditory stimulations were presented throughout the session. The evoked waveforms are presented in figure figure3.3. Here, the upper waveforms represent the averaged response with a distinct N1 component (yellow vertical box points out). The lower lines represent the increasing level of propofol. Accordingly, the highest level of propofol displays no clear N1 waveform. In order to highlight the differences between the conditions, the lower dosages (0.4 and 0.8 µg/mL) and the higher dosage levels (1.2 and 1.6 µg/mL) were grouped. The comparison of higher dosage group with the lower one revealed a significant difference (p=0.003).
A sample sLORETA of group average is provided in figure figure4.4. On the left hand side, the spectrogram of the evoked response for N1 is given. The 0 denotes the stimulus onset time. The major part of the activity lies from delta to alpha ranges (3.5-11 Hz). The sLORETA reveals highest power to be located at Broadman Area (BA) 6. Additionally, with the administration of propofol, especially reaching around 0.8 µg/mL, the distinctive spindle oscillations at around 11-14 Hz were observed. These waveforms lasted approximately for 1 second and were predominantly present in the central line electrodes (Figure (Figure55).
Lower propofol dosage (0.4 and 0.8 µg/mL) level revealed sLORETA locations of maximal activity areas to include, BA6 (5 subjects), BA7 (7 subjects), BA20 (3 subjects), BA3 (1 subject), BA11 (1 subjects), BA18 (3 subject), BA19 (1 subject), BA38 (1 subject), BA47 (2 subjects).
For the higher propofol dosage (1.2 and 1.6 µg/mL) the locations from the maximal activity areas included, BA6 (6 subjects), BA7 (7 subjects), BA3 (1 subject), BA8 (2 subjects), BA11 (2 subjects), BA19 (2 subjects), BA20 (1 subject), BA21 (1 subject), BA22 (1 subject), BA37 (1 subject).
The Bispectral Index (BIS) measurements have been utilized both for anesthesia and sleep conditions. The BIS values follow the depth of sleep closely. The figure figure66 presents the BIS (blue) and classical R&K stages (red) in vertical axes.
According to Spearman analysis, the correlation rate between BIS values and sleep stages in this sleep period was r =0.82.
During the administration of auditory stimulations, upon changing of the type of stimulation of BBB, we have observed changes in depth of sleep of subjects. These changes have been noted e.g. figure figure77 especially with the presentation of syllables (DL) after simple tones (DT). The second application of syllables resulted in a similar increase in BIS values (arrow in Fig. Fig.7).7). As every subject would have a unique sleep pattern-and it is not predictable- the causal relationship of external stimulation changes on the sleep patterns remains a further study topic.
The auditory evoked potentials have been succesfully obtained during different stages of sleep. The figure figure88 and Table Table33 represent the N1 potential across different BIS values. “The lower the BIS value deeper the sleep” was already obtained above (figure (figure6)6) . Therefore the lower BIS values in figure figure88 represent deeper sleep stages through which the N1 wave starts to dissappear.
Similar to the anesthesia, the higher BIS group (80 and 70) and the lower BIS group (60 and 50) were grouped. The comparison of higher BIS group with the lower one revealed a significant difference (p=0.005).
A sample sLORETA of group average is provided in figure figure9.9. On the left hand side, the spectrogram of the evoked response for N1 during sleep is given. The major part of the activity lies around delta with some activity towards theta ranges (3.5-7 Hz). The sLORETA reveals highest power to be located at Broadman Area 7.
The current manuscript presents data from different experimental setups and from different conscious states. Therefore a confined approach was applied. Accordingly, the effect of simple auditory stimulations in the range of 100-150 ms was focused. The anesthesia procedure by means of step-wise administration of propofol, has displayed a gradual loss of N1. Likewise with deepening of sleep -which could be monitored also with BIS- the N1 component diminished.
It has been reported that approximately after the 100-150 ms of the auditory stimulation a waveform which has a negative peak was observed in awake subjects [6,27]. This waveform is called (N100) “N1”. The neuronal generators for this waveform are suggested to be planum temporale in the secondary cortical area .
The current results revealed activities located at around BA3, BA6, BA7, BA8, BA11, BA19, BA20, BA21, BA22, BA37, BA38 and BA47 regions during the auditory stimulations under propofol anesthesia. These locations were found across the subjects using the GFP maxima matching the N1 peak (113 ms). The presence of BA6 and BA7 were more frequent than the other locations. The BA7 area had been located in other studies highlighting the functional prospects of precuneus [29-31]. Additionally, the vegetative state has been linked to impaired activity at this area using PET imaging .
Furthermore the current sleep experiment findings revealed locations around BA4, BA7, and BA9. Likewise, in the study of Anderer , the sleep spindles which are more frequent in Stage 2 (and thus around BIS level 70) were found to be related to BA9 and BA10. In , BA7, BA9, and BA10 were among the common activation areas. The activation of these areas has been linked to sleep spindles. Accordingly the areas from the current study are plausible target areas of sleep and anaesthesia domains. Interestingly, the sleep and anaesthesia share similar areas to a degree for the functional changes. The short term spectral changes like in the case of spindles, need to be further addressed in parallel experimental designs.
With the administration of propofol, especially reaching around 0.8 µg/mL the distinctive spindle oscillations at around 11-14 Hz were observed. At this level also the subjects started losing their consciousness and secondly the N1 waveform began to disappear. While the formation of LOC, the disappearance of N1 and the appearance of spindle oscillations can be coincidental, there have been various reports supporting the case [33-35]. The presence of this type oscillatory activity has been attributed to inductive versus disruptive properties [35-38]. Therefore the dual nature of these activities needs to be unveiled via controlled studies.
The sleep related EEG segments and the spectrograms also revealed spindle activity at around BIS levels of 80-70. The causal mechanisms and relationship of these spindles to anesthesia are beyond the scope of current manuscript, however there have been various studies indicating functional aspects of the spindles [22,33,35,39,40].
The spatiotemporal properties of brain functioning undergo some degree of changes during alterations in the conscious states. Some of these mechanisms may be linked to minor to large scale networks in the brain [41-43]. These changes may be necessary for brain to shift from one state to another (i.e. sleep), or they can be direct effects of medications (i.e. propofol anesthesia).
The disappearance of basic electrophysiological patterns does not necessarily point to cessation of auditory cognitive processing. This could be analyzed with implicit memory paradigms as well as other special designs. Henceforth, caution is necessary while commenting on electrophysiological data in relationship to cognition.
The first stage results have provided a basic understanding in these altered states such that various auditory stimuli are successfully processed in both light and deep sleep stages. In addition to information processing during sleep, a large group of researchers is concerned with sleep disorders [44-46]. However; all these studies are still far away from elucidating the sleep field. In order to shed light on sleep processes, the brain functioning in wakefulness shouldn’t be completely neglected. It has been indicated that sleep has a marked effect on N1 waveform . Researchers have stated that at the beginning of sleep there is a prolongation in latency and decrease in amplitude of N100 component [29,47]. It has been also reported that N100 waveform is highly sensitive to physical properties of the stimulus and with the repetition of the same stimuli N100 amplitude decreases . It was stated that during NREM sleep amplitude of N100 decreases and the latency lengthens . In a research it has been found that the amplitude of N100 reduces in slow wave sleep and slightly increases in REM sleep . Näätänen has pronounced that in alert wakefulness the N100 component is correlated with task performance . As can be seen from the above mentioned studies, many researchers have pointed that N100 waveform is susceptible to the state of the brain. A common feature of stage 2, the sleep spindles, is thought to be related with suppression of information processing . It has been associated that decrease in N100 amplitude in N2 sleep to the aforementioned function of the sleep spindles . Atienza et al. have interpreted the decrease of N100 amplitude as reduced afferent sensorial processing at the subcortical level  In the light of these results one could explain the decrease in N100 amplitude in deeper sleep stages in our study to be due to the suppression of auditory stimulation in order to maintain sleeping.
From a higher scope, the stages of sleep are artificial separations, whereas during natural sleep the brain may shift abruptly from one state to another. The external and internal conditions may all play a role in this sleep equilibrium. These shorter term stages may be defined as microstates. The external stimulations do have effect on immediate neurocognitive properties. Additionally, these stimulations may also play a causative role for the brain change its response state. Therefore, the observations that we have stated in our study will remain as a study area for our group. The brain biophysics battery is made up of different auditory stimulation blocks. Hence the different properties of brain responsiveness can be addressed while the block changes might serve to highlight the microstates level.
The field of studying brain functioning as in the present approach, would benefit from further experiments or research methods. These can be enriched by incorporating “Directed information transfer”, “Entropy and coherence”, “Causal relationships”, and “Brain asymmetry (i.e. dichotomy) under different conditions”. Additionally, an extended inclusion of various parameters such as body temperature and other physiological data might be beneficial to explain further mechanisms of the brain that interact with the body.
Both sleep and anesthesia are dynamic conditions which require a uniform approach. Auditory stimulations with distinct features may provide a thorough insight into the brain responsiveness in different conscious states. Additionally spectro-temporal properties of the dynamic states and the analysis of microstates can constitute to better understanding of underlying mechanisms. The other line of plausible applications includes the neuropsychiatric pathologies, coma, and other major states of the brain.
The Embedded Interactive Stimulus Unit (EMISU) system has a “patent pending” status.
MO conceived of the study, and constructed its design and coordination. OB participated in the anesthesia experiments and supported the sLORETA assessments, SK participated in the sleep experiments and supported with the analysis of the data, NG participated in the anesthesia experiments and AO participated in the design of the study and manuscript. All authors read and approved the final manuscript.
ABR: Auditory Brain-stem Responses; AEP: Auditory Evoked Potentials; AERP: Auditory Event Related Potential; BA: Broadman Area; BBB: Brain Biophysics Battery; BIS: Bispectral Index; CDV: Current Density Value; DL: Dichotic Listening; DT: Dichotic Tone; ECoG: Electrocorticogram; EEG: Electroencephalography; EMISU: Embedded Microcontroller Interactive Stimulus Unit; EOG: Electrooculography; fMRI: functional Magnetic Resonance Imaging; GFP: Global Field Power; ISI: Inter Stimulus Interval; LOC: Loss Of Consciousness; MEG: Magneto Encephalography; MLRs: Middle Latency Responses; MNI: Montreal Neurological Institute; MMN: Mismatch Negativity; PET: Positron Emission Tomography; R&K: Rechtschaffen and Kales; ROI: Region Of Interest; sLORETA: standardized Low Resolution Electromagnetic Tomography; SPL: Sound Pressure Level; SQI: Signal Quality Index)
The authors are grateful to Emre Cek for his support on wavelet based spectral analysis tool and to Cagdas Guducu for his support on the graphics and data sets; and to Ugras Erdogan for participation in the design of the EMISU system and to Dr. Ozgul Vupa for her support on statistical analysis. We thank to the neurosurgery and anesthesia teams. The study was supported by means of projects DEU.2006.KB.SAG.017, DEU.2007.KB.SAG.061/65 and TUBITAK 108S113, COST BM601 (Neuromath). The study was partially presented at the conference of “Consciousness and its Measures”, November 29 - December 1, 2009 Limassol, Cyprus (http://www.humanbraindynamics.com/neuroconsci2009).
This article has been published as part of Nonlinear Biomedical Physics Volume 4 Supplement 1, 2010: Consciousness and its Measures: Joint Workshop for COST Actions Neuromath and Consciousness. The full contents of the supplement are available online at http://www.nonlinearbiomedphys.com/supplements/4/S1.