The use of in ovo electroporation for the rapid analysis of neural-specific murine enhancers

The identification and characterization of DNA sequences necessary for proper gene expression have provided insights into gene regulation and generated tools useful for further experimentation. Studies of developmentally regulated genes have demonstrated how transcription factors interact at enhancers to generate restricted patterns of expression during embryogenesis. In vertebrates, the pursuit of such studies has relied on the generation of transgenic mice and thus has been limited by the time and expense required generating and characterizing these mice. The recently developed technique of in ovo electroporation allows the rapid introduction of exogenous DNA into developing chicken embryos. Here we have used this technique to introduce DNA containing murine enhancer/reporter constructs into cells of the chicken neural tube, resulting in appropriate expression of the reporter. This technique has the potential to greatly reduce the effort involved in the study of vertebrate enhancers. Furthermore, we have characterized factors such as timing of electroporation, concentration of DNA, and choice of basal promoters and found that they can influence the degree to which expression of enhancer constructs reflects endogenous gene expression.

Model selection in non-nested hidden Markov models for ion channel gating

An important task in the application of Markov models to the analysis of ion channel data is the determination of the correct gating scheme of the ion channel under investigation. Some prior knowledge from other experiments can reduce significantly the number of possible models. If these models are standard statistical procedures nested like likelihood ratio testing, provide reliable selection methods. In the case of non-nested models, information criteria like AIC, BIC, etc., are used. However, it is not known if any of these criteria provide a reliable selection method and which is the best one in the context of ion channel gating. We provide an alternative approach to model selection in the case of non-nested models with an equal number of open and closed states. The models to choose from are embedded in a properly defined general model. Therefore, we circumvent the problems of model selection in the non-nested case and can apply model selection procedures for nested models.

Tremor-correlated cortical activity in essential tremor

BACKGROUND

In patients with parkinsonian resting tremor, tremor-correlated activity in the contralateral sensorimotor cortex has been studied by both magnetoencephalography (MEG) and electroencephalography (EEG). In essential tremor, MEG failed to detect cortical involvement. The objective of this study was to investigate whether EEG recording can reveal tremor-correlated cortical activity in patients with essential tremor or enhanced physiological tremor.

METHODS

Seven patients with essential tremor and three patients with enhanced physiological tremor participated in the study. Unilateral postural tremor was activated by wrist extension on the right or on the left side. Electromyography (EMG) signals arising from the wrist extensor and flexor muscles, and a high-resolution EEG were recorded simultaneously. Coherences between the time series of the rectified tremor EMG and the EEG were estimated.

FINDINGS

In five of nine arms with essential tremor, we found highly significant coherences at the tremor frequency between the tremor EMG and the EEG. Isocoherence maps illustrating the topography of significant coherences over the scalp showed that the maximum coherences were located over the contralateral sensorimotor cortex. In the patients with enhanced physiological tremor, we were unable to detect consistent significant corticomuscular coherences at the tremor frequency.

INTERPRETATION

Using simultaneous EEG-EMG recordings, we showed that significant corticomuscular coherences at the tremor frequency can be found in essential tremor. This finding contrasts with a recent study based on MEG recordings. The results suggest that the sensorimotor cortex is involved in the generation of essential tremor, in a similar way to that previously shown in parkinsonian resting tremor.

The effects of non-identifiability on testing for detailed balance in aggregated Markov models for ion-channel gating

Aggregated Markov models are a widely used tool to model patch clamp data measured from single ion channels. These channels must obey the principle of detailed balance in thermodynamic equilibrium; otherwise, the channel is driven by an external source of energy. We investigate the power of a likelihood ratio test for detailed balance for a number of data points which is in the order of magnitude of patch clamp experiments. We show that for certain models with nearly equal dwell times, a test for detailed balance suffers from a loss of power to detect violations of detailed balance which is due to the non-identifiability of the transition rates for models with equal dwell times.

Detection of obstructive sleep apnea by analysis of phase angle using the forced oscillation signal

Pharyngeal collapse in patients with obstructive sleep apnea syndrome (OSAS) is linked to decreased upper airway muscle activity. We hypothesised that decreased muscle activity causes decreased stiffness of the upper airway wall and assumed that a decrease in wall stiffness would result in a change point (CP) of the morphology of phase angle time series phi(t) obtained by forced oscillation technique (FOT). We developed an algorithm to detect CPs in φ(t) and investigated phi(t) data measured in parallel to all-night polysomnography in seven patients. A total of 2724 CPs were detected by algorithm. The CPs were marked on the polysomnograms and compared with polysomnogram scoring performed jointly by three sleep expert physicians. A total of 791 (67%) of the 1172 handscored respiratory events showed a CP in phi(t) within a time interval of 8 sec before and 6 sec after the marked onset. A total of 672 (57%) respiratory events were detected at an earlier time by algorithm, and 119 (10%) were algorithmically detected later. The average detection time by the CP algorithm was 4.0+/-1.1 sec earlier than the manual scoring. We conjecture that a CP in phi(t) indicates a change in upper airway collapsibility and that phi(t) might be a potentially useful parameter for detection of impending upper airway obstruction.

The n-PI-method for helical cone-beam CT

A new class of acquisition schemes for helical cone-beam computed tomography (CB-CT) scanning is introduced, and their effect on the reconstruction methods is analyzed. These acquisition schemes are based on a new detector shape that is bounded by the helix. It will be shown that the data acquired with these schemes are compatible with exact reconstruction methods, and the adaptation of exact reconstruction algorithms to the new acquisition geometry is described. At the same time, the so-called PI-sufficiency condition is fulfilled. Moreover, a good fit to the acquisition requirements of the various medical applications of cone-beam CT is achieved. In contrast to other helical cone-beam acquisition and reconstruction methods, the n-PI-method introduced in this publication allows for variable pitches of the acquisition helix. This additional feature will introduce a higher flexibility into the acquisition protocols of future medical cone-beam scanners. An approximative n-PI-filtered backprojection (n-PI-FBP) reconstruction method is presented and verified. It yields convincing image quality.

On identification of Na(+) channel gating schemes using moving-average filtered hidden Markov models

Transitions between distinct kinetic states of an ion channel are described by a Markov process. Hidden Markov models (HMM) have been successfully applied in the analysis of single ion channel recordings with a small signal-to-noise ratio. However, we have recently shown that the anti-aliasing low-pass filter misleads parameter estimation. Here, we show for the case of a Na(+) channel recording that the standard HMM do neither allow parameter estimation nor a correct identification of the gating scheme. In particular, the number of closed and open states is determined incorrectly, whereas a modified HMM considering the anti-aliasing filter (moving-average filtered HMM) is able to reproduce the characteristic properties of the time series and to perform gating scheme identification.

Tremor-correlated cortical activity detected by electroencephalography

OBJECTIVE

In this study we investigated whether cortical activity related to Parkinsonian resting tremor can be detected by electroencephalography (EEG).

METHODS

Seven patients with idiopathic Parkinson's disease suffering from unilateral tremor participated in the study. Electromyography (EMG) signals arising from the wrist extensor and flexor muscles as well as a high resolution EEG were recorded simultaneously. Coherencies between EEG and EMG were calculated.

RESULTS

In all patients, we found highly significant coherencies at the tremor frequency or its first harmonic between the tremor EMG and contralateral EEG channels. There were no significant coherencies between the tremor EMG and ipsilateral EEG channels. Isocoherency maps illustrating the topography of the coherencies over the scalp showed that the maximum coherencies were situated over the cortical motor areas. In one case, a high coherency was also found over the parietal cortex.

CONCLUSIONS

The results show for the first time that tremor-correlated cortical activity can be detected by electroencephalography. The findings underline that motor areas of the cerebral cortex are involved in the neuronal network generating resting tremor in Parkinson's disease.

Pathological tremors: Deterministic chaos or nonlinear stochastic oscillators?

Pathological tremors exhibit a nonlinear oscillation that is not strictly periodic. We investigate whether the deviation from periodicity is due to nonlinear deterministic chaotic dynamics or due to nonlinear stochastic dynamics. To do so, we apply various methods from linear and nonlinear time series analysis to tremor time series. The results of the different methods suggest that the considered types of pathological tremors represent nonlinear stochastic second order processes. Finally, we evaluate whether two earlier proposed features capturing nonlinear effects in the time series allow for a discrimination between two pathological forms of tremor for a much larger sample of time series than previously investigated. (c) 2000 American Institute of Physics.

Linear and nonlinear time series analysis of the black hole candidate cygnus X-1

We analyze the variability in the x-ray lightcurves of the black hole candidate Cygnus X-1 by linear and nonlinear time series analysis methods. While a linear model describes the overall second order properties of the observed data well, surrogate data analysis reveals a significant deviation from linearity. We discuss the relation between shot noise models usually applied to analyze these data and linear stochastic autoregressive models. We debate statistical and interpretational issues of surrogate data testing for the present context. Finally, we suggest a combination of tools from linear and nonlinear time series analysis methods as a procedure to test the predictions of astrophysical models on observed data.

Selecting relevant electrode positions for classification tasks based on the electro-encephalogram

The aim is to describe a general approach to determining important electrode positions when measured electro-encephalogram signals are used for classification. The approach is exemplified in the frame of the brain-computer interface, which crucially depends on the classification of different brain states. To classify two brain states, e.g. planning of movement of right and left index fingers, three different approaches are compared: classification using a physiologically motivated set of four electrodes, a set determined by principal component analysis and electrodes determined by spatial pattern analysis. Spatial pattern analysis enhances the classification rate significantly from 61.3 +/- 1.8% (with four electrodes) to 71.8 +/- 1.4%, whereas the classification rate using principal component analysis is significantly lower (65.2 +/- 1.4%). Most of the 61 electrodes used have no influence on the classification rate, so that, in future experiments, the setup can be simplified drastically to six to eight electrodes without loss of information.

Side-to-side correlation of muscle activity in physiological and pathological human tremors

OBJECTIVE

Many tremors occur always or often bilaterally. The question arises whether this could be explained by a common source or commonly transmitting pathways or by bilaterally represented, independent structures with the same oscillatory properties. A similar tremor frequency does not provide sufficient information to clarify this question.

METHODS

We analyze coherencies between surface electromyographies (EMG) to investigate if bilateral physiologic (PT), essential (ET), Parkinsonian (PD) and orthostatic (OT) tremors originate from a common source for both sides of the body. We show that commonly used techniques to test whether coherencies are significant could lead to false positive results for tremor EMGs. A new estimation procedure is proposed to test EMG tremor time series on their linear independence. We apply this test to bilateral tremors.

RESULTS

All measured EMG-pairs in OT (n = 7) were highly coherent between both sides with reproducible coherency values of up to 0.99. All other investigated tremors, i.e. PT and enhanced physiological tremors (EPT, n = 117), ET (n = 76) and PD resting and postural tremors (n = 70) do not show a significant side-to-side correlation.

CONCLUSIONS

This finding shows that the pathophysiologies of OT and other pathological tremors are definitely different. Either they have different origins or different kinds of transmitting pathways. The proposed method might also be used to investigate other electrophysiological data and is a helpful, easy to use investigation for a daily clinical routine.

Classical conditioning of the electrically elicited blink reflex in humans: a new method of data analysis

The eyeblink conditioning paradigm is a well-established model to study learning processes in humans and animals. Especially results from animal studies have supplied new insight into physiological pathways and brain structures involved in associative motor learning and memory. An important role of the cerebellum and its afferent fiber systems could be shown. Recent studies in humans have given evidence that results of animal experiments can be applied directly to the human condition. A high variation of baseline EMG activity and/or spontaneous blinks may influence the analysis of classical conditioning of the electrically elicited blink reflex in humans. To optimize differentiation between real conditioned responses and stimulus-independent EMG activity, we developed an algorithm which is fully automated and independent of a possible bias of an examiner. In a first step the algorithm decides whether a subject fulfills the criteria of a successful learning process or not. The second step quantifies the learning process. For quantification of the learning process, the following parameters were calculated: number of conditioned responses, onset of conditioning, time and amount of maximal conditioning, speed of conditioning and speed of habituation. According to our criteria, 80% of the healthy volunteers acquired conditioned responses. There is an age-related decline in eyeblink classical conditioning. Analysis of patient groups with different types of lesions will further improve our knowledge and understanding of pathways involved in learning processes in humans. The proposed new algorithm of data analysis takes less than 10 s on a standard computer, is more sensitive and more specific in detecting conditioned responses and, therefore, may further improve the value and reliability of the eyeblink conditioning paradigm in clinical research.

Early detection of upper airway obstructions by analysis of acoustical respiratory input impedance

Repetitive occurrence of partial or total upper airway obstruction characterizes several respiratory dysfunctions such as the obstructive sleep apnea syndrome (OSAS). In OSAS patients, pharyngeal collapses are linked to a decrease in upper airway muscle activity during sleep which causes decreased upper airway wall stiffness. Continuous positive airway pressure (CPAP) is recommended as the treatment of choice. Advancements in CPAP therapy require early detection of respiratory events in real time to adapt the level of the applied pressure to airway collapsibility. The forced oscillation technique (FOT) is a noninvasive method which reflects patients' airway patency by measuring respiratory impedance. The aim of this study was to evaluate by a mathematical model of the respiratory system if FOT can provide an early detection index of total or partial upper airway obstruction. Furthermore, the simulation should suggest which characteristic features are relevant for early apnea detection in measured clinical data. The respiratory system has been treated as a series of cylindrical segments. The oropharynx analog of the model allows simulation of upper airway collapse, mimicking the situation in patients with OSAS. We calculated the input impedance for different degrees of upper airway obstruction ranging from unobstructed airways to total occlusion. Furthermore, we simulated different upper airway wall compliances. We compared the simulation with real data. The results of the study suggest that FOT is a valuable tool for assessing the degree of upper airway obstruction in patients with OSAS. Especially, the phase angle of the impedance seems to be a potentially useful tool for early apnea detection by assessing the upper airway wall collapsibility.

A software for recording and analysis of human tremor

For many diseases various methods for the diagnosis and treatment monitoring are available. Presently, such methods are not established for an investigation of tremor diseases, although the different forms of tremor are common neurological symptoms and occur frequently in various neurological diseases and also other conditions. We developed an easy-to-use application for tremor-analysis and recording, running under MS-Windows, that allows us to investigate different forms of tremor by advanced mathematical methods of time series analysis. The application is also applicable for users who are not familiar with these kind of advanced data analysis methods. It provides tools for the diagnosis and treatment monitoring under laboratory conditions, based on previously developed and established methods of spectral and cross spectral analysis of tremor and electromyographic time series.

Artefacts in spiral-CT images and their relation to pitch and subject morphology

This qualitative study is intended to create awareness of artefacts that are associated with spiral-CT imaging. A simple description of spiral-CT reconstruction is used to explain how these artefacts depend on the pitch and subject morphology, and shows when these artefacts are likely to impair the diagnostic value of the acquired images. We scanned a cone and rod phantom with pitch 2, and used the acquired images to demonstrate how spiral data acquisition and interpolation leads to artefacts in the reconstructed images. We then demonstrated the effects of various pitches in scans of a human cadaver, whereas the slice thickness was kept constant. Some patient studies are presented in order to show the possible clinical consequences. Spiral acquisition may cause geometric distortions and apparent inhomogeneity of homogeneous structures. We were able to link these artefacts to the way in the acquisitions were done, and the reconstructions were performed. We have shown how these artefacts can be anticipated in clinical studies. When areas of low contrast, surrounded by hypo- or hyperdense structures, are scanned with a large pitch and viewed with a narrow window, spiral artefacts may influence the diagnostic quality of the images. These effects should be considered when choosing the pitch.

Involvement of cranial muscles and high intermuscular coherence in orthostatic tremor

Electromyographic recordings were conducted from limb, trunk, and cranial muscles in 6 patients with orthostatic tremor. Spectral analysis revealed a high-frequency tremor not only in the muscles of the limbs and trunk, but also in cranial muscles. The cross spectra were analyzed between various pairs of muscles that displayed a high-frequency tremor pattern. The resulting peak correlations were uniformly very high (near one) suggesting a high level of coherence. The involvement of cranial muscles suggests that supraspinal mechanisms are involved in the generation of orthostatic tremor. The high intermuscular coherence between all muscles indicates the existence of either a unique oscillator that generates tremor in all involved muscles on both sides of the body or a linking mechanism probably at a supraspinal level. The high-frequency tremor was only found when the muscles were contracted isometrically, irrespective of body posture. Thus, tremor generation might be more closely linked to mechanisms responsible for isometric force control than to those involved in stance regulation.

Evaluation of forced oscillation technique for early detection of airway obstruction in sleep apnea: a model study

The forced oscillation technique (FOT) is a non-invasive method which may be suitable for assessing upper airway obstruction in obstructive sleep apnea/hypopnea syndrome (OSAS) patients. The aim of this study was to determine in vitro if FOT can provide an early detection index of total or partial upper airway occlusion. A respiratory system analog was developed which includes an upper airway analog that allows simulation of upper airway collapse (thus mimicking the situation in patients with OSAS). We simulated different degrees of upper airway obstructions ranging from 0 (unobstructed airways) to 1 (total occlusion). Furthermore, we imitated the collapsible upper airway wall by means of elastic membranes with ten different wall compliances ranging from 3.3 x 10(-4) to 1 1/mbar. For the two stiffest rubber membranes (Cwall = 0.01 and 3.3 x 10(-4) l/mbar) the absolute value of the impedance (¿Z¿) showed a marked increase for obstructions greater than 0.6. For the two membranes with the highest wall compliances (Cwall = 0.03 and 1 1/mbar) obstructions with an increase in ¿Z¿ could not be detected before the obstruction reached 0.8. For degrees of obstruction less than 0.6 the phase angle of collapsible membranes with low compliance (stiff airway wall) were about 1.5pi which significantly differed from phase angles of 1.77pi measured in membranes with high compliance (elastic airway wall); p < 0.01. We hypothesized that stiffness of upper airway walls corresponds with their muscle tone, i.e., stiff airway walls are related with high muscle tone and vice versa. Thus, a decrease in upper airway muscle activity would cause an increase of upper airway wall elasticity that enables upper airway collapse. As a consequence the phase angle phi could be expected to change from values characterizing stiff membranes to values characterizing more elastic membranes which could be used as early indicator for obstructive respiratory events. We have frequently observed such changes in morphology of phi(t) data obtained from patients with OSAS.

Cross-spectral analysis of physiological tremor and muscle activity. II. Application to synchronized electromyogram

We investigated the relationship between synchronized muscle activity and tremor time series in (enhanced) physiological tremor by cross-spectral analysis. Special attention was directed to the phase spectrum and its potential to clarify the contribution of reflex mechanisms to physiological tremor. The phase spectra are investigated assuming that the electromyogram (EMG) synchronization was caused by a reflex or a central oscillator. Comparing these results to phase spectra of measured data, we found a significant contribution of reflexes. But reflexes only modify existing peaks in the power spectrum. The main agents of physiological tremor are an efferent pace and the resonant behavior of the biomechanical system.

Cross-spectral analysis of physiological tremor and muscle activity. I. Theory and application to unsynchronized electromyogram

We investigate the relationship between the extensor electromyogram (EMG) and tremor times series in physiological hand tremor by cross-spectral analysis. Special attention is directed to the phase spectrum and the effects of observational noise. We calculate the theoretical phase spectrum for a second-order linear stochastic process and compare the results to measured tremor data recorded from subjects who did not show a synchronized EMG activity in the corresponding extensor muscle. The results show that physiological tremor is well described by the proposed model and that the measured EMG represents a Newtonian force by which the muscle acts on the hand.

Detecting multimodality in saccadic reaction time distributions in gap and overlap tasks

In many cases the distribution of saccadic reaction times (SRT) deviates considerably from a unimodal distribution and may often exhibit several peaks. We present a statistical approach to determining the number and form of the individual peaks. The overall density of the reaction times fi(t), i = 1...M obtained in M different experiments with the same subject is described as the sum of K basis functions xk(t), k = 1...K with different weights and an error term. A change in the experimental conditions is assumed to cause a change in the weights, not in the basis functions. We minimize the square of the difference (measured data minus approximation), divided by the error of the data. Incrementing K step by step we determine the necessary number of basis functions. This method is applied to data of six subjects tested in different saccade tasks. We detect five different modes: two in the range 80-140 ms (express modes), two in the range 145-190 ms (fast-regular mode) and one at about 230 ms (slow-regular mode). These modes are located at about the same positions for different subjects. The method presented here not only proves statistically the existence of several modes in SRT distributions but also allows the distributions to be described by a few characteristic numbers that go beyond the mean values and standard deviations.

Numerical methods to determine calcium release flux from calcium transients in muscle cells

Several methods are currently in use to estimate the rate of depolarization-induced calcium release in muscle cells from measured calcium transients. One approach first characterizes calcium removal of the cell. This is done by determining parameters of a reaction scheme from a fit to the decay of elevated calcium after the depolarizing stimulus. In a second step, the release rate during depolarization is estimated based on the fitted model. Using simulated calcium transients with known underlying release rates, we tested the fidelity of this analysis in determining the time course of calcium release under different conditions. The analysis reproduced in a satisfactory way the characteristics of the input release rate, even when the assumption that release had ended before the start of the fitting interval was severely violated. Equally good reconstructions of the release rate time course could be obtained when the model used for the analysis differed in structure from the one used for simulating the data. We tested the application of a new strategy (multiple shooting) for fitting parameters in nonlinear differential equation systems. This procedure rendered the analysis less sensitive to ill-chosen initial guesses of the parameters and to noise. A locally adaptive kernel estimator for calculating numerical derivatives allowed good reconstructions of the original release rate time course from noisy calcium transients when other methods failed.

Central mechanisms in human enhanced physiological tremor

The sites of the central nervous structures involved in enhanced physiological tremor (EPT) are still unclear. The syndrome of persistent mirror movements (PMM) is characterized by abnormal bilateral corticospinal projections. If a supraspinal mechanism is involved in EPT, the activity of EPT should be coherent between both sides in subjects with this abnormality. We investigated three PMM subjects and three normal controls. Focal transcranial magnetic stimulation (TMS) resulted in contralateral hand muscle responses in the controls. The PMM subjects, in contrast, had bilateral responses. Similarly, long-latency reflexes (LLR) in PMM could be recorded bilaterally, while the control subjects showed responses only on the stimulated side. EPT was evoked by intravenous salbutamol. EMG time series were recorded bilaterally from the wrist extensor muscles and cross spectra were calculated. If there was a significant right-left-coherence, phase analysis was performed. No control subject showed a significant right-left-coherence of tremor activity. In contrast, a significant coherence was found in PMM between 8 and 12 Hz. When the mechanical tremor frequency of one hand was reduced by loading, coherences and phase spectra of the EMGs remained unchanged. By comparing the results from TMS, LLR and cross spectral analysis we come to the conclusion, that the 8 to 12 Hz component of EPT is transmitted transcortically, most likely originating from two separate generators for both sides.

Analysis of multichannel patch clamp recordings by hidden Markov models

Conventional methods of analysis do not allow the kinetics of patch clamped ion channels to be completely determined if more than one channel is present in the patch. This hinders investigations on small ion channels as well as on channel cooperativity and the homogeneity of channel populations. We present a method to extract the rate constants and current amplitudes for each individual channel from multichannel patches by a one-step procedure. For this purpose, the current record is modeled by the superposed Markov processes of the opening and closing of each channel that is contaminated by noise (Hidden Markov Model). Channel parameters are obtained by maximum likelihood methods. Because the parameters can be calculated directly from the unfiltered record, the dwell time and missed event problems are widely diminished. Confidence bounds for the estimated parameters are given. Statistical tests to decide whether channels switch identically and/or independently are introduced. The application of the method is demonstrated with simulated data.

Modeling volatility using state space models

In time series problems, noise can be divided into two categories: dynamic noise which drives the process, and observational noise which is added in the measurement process, but does not influence future values of the system. In this framework, we show that empirical volatilities (the squared relative returns of prices) exhibit a significant amount of observational noise. To model and predict their time evolution adequately, we estimate state space models that explicitly include observational noise. We obtain relaxation times for shocks in the logarithm of volatility ranging from three weeks (for foreign exchange) to three to five months (for stock indices). In most cases, a two-dimensional hidden state is required to yield residuals that are consistent with white noise. We compare these results with ordinary autoregressive models (without a hidden state) and find that autoregressive models underestimate the relaxation times by about two orders of magnitude since they do not distinguish between observational and dynamic noise. This new interpretation of the dynamics of volatility in terms of relaxators in a state space model carries over to stochastic volatility models and to GARCH models, and is useful for several problems in finance, including risk management and the pricing of derivative securities. Data sets used: Olsen & Associates high frequency DEM/USD foreign exchange rates (8 years). Nikkei 225 index (40 years). Dow Jones Industrial Average (25 years).

Saccadic reaction times: a statistical analysis of multimodal distributions

The distributions of saccadic reaction times (SRT) often deviate from unimodal normal distributions. An excess-mass procedure was used to detect peaks in 963 data sets containing 90,927 reaction times from 170 subjects. About 55% showed one, 30% two, 12% three and 3% four peaks. According to their clustering along the reaction time scale the modes could be classified into express (90-120 msec), fast regular (135-170 msec) and slow regular (200-220 msec) modes. Among the unimodal distributions 29% had peaks in the range of the express mode and 46% had peaks in the range of the fast regular mode. Therefore, 87% of the data sets support the notion of saccadic reaction time distributions being the superposition of three modes. All experimental distributions were fitted by as many gamma distributions as determined by the excess-mass test. The significance of the multimodality for saccade generation processes is discussed.

Increased sensitivity of the inositol-phospholipid system in neutrophils from patients with acute major depressive episodes

Evidence from in vitro and in vivo studies suggests that the therapeutic and prophylactic effects of lithium in recurrent affective disorders are due to an attenuation of the inositol-phospholipid (IPL) second messenger system. An increased sensitivity of this signal transduction system might therefore constitute a risk factor for affective illness. The extent of the agonist-induced release of intracellular Ca2+ (Ca2+ response) can be used as an indicator of the sensitivity of the IPL system. Using this paradigm, we have measured the agonist-induced Ca2+ response in neutrophils of 17 unmedicated patients who were experiencing an acute major depressive episode. The neutrophils were stimulated by the chemotactic peptide formylmethionylleucylphenylalanine, which activates the IPL system in the cells. The sensitivity of the IPL system in these patients was significantly greater (dose-response curve shifted to the left) compared with its sensitivity in healthy age- and sex-matched control subjects. The results indicate that acute episodes of major depression are associated with an increased sensitivity of the IPL system.

Quantitative analysis of tremor time series

Spectral analysis is applied to tremor time series in basic research and treatment monitoring. The estimation of the spectra from the data is usually done by averaging the squared modulus of the Fourier transform of segments of the data. We discuss drawbacks of this method and propose an alternative procedure to estimate the spectra adaptively based on the data. Thus, the method can be applied to all types of tremor. Applying the theory of spectral estimation, we propose a method to decide whether a spectrum exhibits multiple significant peaks and discuss different approaches to determine the amplitude of the tremor from the spectrum.

Both adenosine A1- and A2-receptors are required to stimulate microglial proliferation

The neuromodulator adenosine is one of the major endogenous inhibitors of overactive excitatory neurotransmission. Adenosine receptors have been identified on neuronal but also on glial surfaces, indicating a role of glial cells in mediation of adenosine effects. Microglia, the immunocompetent cells of the brain, typically respond with proliferation, migration and production of inflammatory substances to viral or bacterial stimuli or to cell damage and degeneration. Since adenosine is released in large amounts in conditions of, for example, hypoxic or ischemic stress, it might be involved in the activation process of microglia. Proliferation of microglia was determined by incorporation of [3H]thymidine into microglial DNA after stimulation with adenosine A1- and A2-receptor agonists. N6-Cyclopentyl adenosine (CPA) and CGS-21680, a specific adenosine A2-receptor agonist had no effect on microglial proliferation. However, combinations of CPA and CGS-21680 as well as the mixed agonist, N6-ethyl-carboxamido adenosine (NECA) increased incorporation of radiolabel above controls. The effect of NECA was inhibited by the adenosine A1-receptor antagonist 8-cyclopentyl-1,3-dipropylxanthine (DPCPX). From these results, it is concluded that proliferation of microglia can be increased only by simultaneous stimulation of both adenosine A1- and A2-receptors. Targeted interference with the activation of A1-adenosine receptors by specific drugs appears to be sufficient to reduce microglial activation. The findings may have implications for the treatment of neurodegenerative diseases in which microglial activation is supposed to play a causative role.

Clinical neurophysiology of tremor

The neurophysiological analysis of tremor has a long tradition. These attempts were directed to understand the mechanisms underlying tremor, on the one hand, and to develop tools to better diagnose the different types of tremor, on the other. Meanwhile, reasonable criteria are available to distinguish between centrally and peripherally mediated tremors. However, no generally accepted means exist to differentiate the different forms of central tremors. Frequency is a useful classifier for cerebellar tremor, rubral tremor, and orthostatic tremor. Although the highest amplitudes are found in Parkinson's disease, this parameter does not well distinguish between the different tremors. Waveform analysis of tremor is a promising tool to separate between the different tremors. Polymyography is pathognomonic for some rare forms of tremor. New approaches to classify tremors are based on positron emission tomography scanning, analysis of ballistic movement, and reflex testing. The means to separate myoclonias from tremors include EEG/EMG correlation techniques, long-latency reflexes, and polymyography. Provided these techniques are applied in the setting of careful clinical analysis of tremor syndromes, they may prove to be helpful in clinical practice.

Some considerations on estimating event-related brain signals

Understanding the timing of mental acts is one of the prominent questions in information processing research. The analysis of event related potentials (ERP) with their high temporal resolution might make access to cognition related brain activity possible. We consider three major problems which make the application of ERPs questionable and then propose some solutions to these problems. The primary problem concerns the separation of the ERPs from the background EEG which is not related to the stimulus. The most common method used is averaging. We argue that this is not the most appropriate method and suggest an alternative for estimating the signal in single-trial recordings. Artifacts present a second problem. We will first review established methods of dealing with eye-movement artifacts and then propose an alternative. We will also report on current work on the parametrisation of single-trial signal estimates, which constitute the third problem considered.

The gridding method for image reconstruction by Fourier transformation

The authors explore a computational method for reconstructing an n-dimensional signal f from a sampled version of its Fourier transform f;. The method involves a window function w; and proceeds in three steps. First, the convolution g;=w;*f; is computed numerically on a Cartesian grid, using the available samples of f;. Then, g=wf is computed via the inverse discrete Fourier transform, and finally f is obtained as g/w. Due to the smoothing effect of the convolution, evaluating w;*f; is much less error prone than merely interpolating f;. The method was originally devised for image reconstruction in radio astronomy, but is actually applicable to a broad range of reconstructive imaging methods, including magnetic resonance imaging and computed tomography. In particular, it provides a fast and accurate alternative to the filtered backprojection. The basic method has several variants with other applications, such as the equidistant resampling of arbitrarily sampled signals or the fast computation of the Radon (Hough) transform.

Use of sodium hydroxide solution to clear partially occluded vascular access ports

Gradual vascular access port or tunneled central venous catheter occlusion during total parenteral nutrition is a common complication of unknown etiology, usually unresponsive to treatment with urokinase, ethanol, or hydrochloric acid. After establishing safety with regard to catheter materials and efficacy in dissolving the occluding substance in vitro, we treated a total of 13 vascular access ports (10 patients) that had occluded gradually during "three-in-one" total parental nutrition by slow perfusion with 10 to 20 mL of 0.1 mmol/mL sodium hydroxide solution. We found a mean increase in flow (results are flow rates under gravity-dependent vascular access port inlet pressure of 1.5 kPa) from 37 microL/s (95% upper confidence limit, 45 microL/s) to 75 microL/s (95% lower confidence limit, 64 microL/s; p < .001). For normal use, a flow of 50 microL/s is sufficient. There were no side effects. Five vascular access ports remained in excellent function during a follow-up period of 9 months. The others had reocclusions, but the partial reocclusions could be treated successfully.

Cytokine production during sleep and wakefulness and its relationship to cortisol in healthy humans

A growing body of evidence indicates that cytokines, especially interleukin-1 beta, are involved in the regulation of sleep and wakefulness. The aim of the present pilot study was to investigate the relationship between interleukin-1 beta (IL-1 beta) and gamma-interferon (gamma-IFN) production and the regulation of sleep and wakefulness. Four healthy male volunteers were investigated. After one adaptation night, beginning at 8 a.m. in the morning, the EEG was recorded by means of a mobile long-term EEG and blood samples were drawn every 45 min for the analysis of IL-1 beta, gamma-IFN and cortisol for 24 h. For the analysis of cytokines whole blood cultures were established. After 48 h of incubation in the presence of endotoxin Salmonella typhimurium, IL-1 beta and gamma-IFN levels were measured in the culture supernatants using specific immunodetection assays. Methods of stochastic time series analysis were adopted to evaluate the biochemical data. Our results show the capability of cultured blood cells to produce cytokines upon endotoxin challenge to be at a maximum around the time of sleep onset and during the first hours of sleep, declining during the night to a minimum level in the morning hours. The opposite was observed for cortisol. The analysis of autocorrelation functions gives evidence of a 24-hour rhythm of cortisol and cytokines. The results indicate that the cytokines IL-1 beta and gamma-IFN may play a role in sleep regulation.

Characteristics of hand tremor time series

Tremor is classified into physiological, essential, and parkinsonian tremor by means of clinical criteria. The aim of our work was to extract quantitative features from the measurements of the acceleration of human postural hand tremor. Different mathematical methods were adopted and modified in order to separate these three types of tremor. Best discrimination between physiological and pathological tremors has been achieved by methods distinguishing nonlinear from linear behavior. On the other hand, methods separating different forms of nonlinear behavior have been found to be superior in discriminating parkinsonian and essential tremor. By these methods physiological and pathological tremors can be separated with an error rate below 20% and essential and parkinsonian tremor with an error rate below 10%. This may help to classify tremor time series by objective mathematical criteria and may increase the understanding of the pathophysiological differences underlying these kinds of tremor.