Topographic time-frequency decomposition of the EEG

Frequency-transformed EEG resting data has been widely used to describe normal and abnormal brain functional states as function of the spectral power in different frequency bands. This has yielded a series of clinically relevant findings. However, by transforming the EEG into the frequency domain, the initially excellent time resolution of time-domain EEG is lost. The topographic time-frequency decomposition is a novel computerized EEG analysis method that combines previously available techniques from time-domain spatial EEG analysis and time-frequency decomposition of single-channel time series. It yields a new, physiologically and statistically plausible topographic time-frequency representation of human multichannel EEG. The original EEG is accounted by the coefficients of a large set of user defined EEG like time-series, which are optimized for maximal spatial smoothness and minimal norm. These coefficients are then reduced to a small number of model scalp field configurations, which vary in intensity as a function of time and frequency. The result is thus a small number of EEG field configurations, each with a corresponding time-frequency (Wigner) plot. The method has several advantages: It does not assume that the data is composed of orthogonal elements, it does not assume stationarity, it produces topographical maps and it allows to include user-defined, specific EEG elements, such as spike and wave patterns. After a formal introduction of the method, several examples are given, which include artificial data and multichannel EEG during different physiological and pathological conditions.

Code-specific learning rules improve action selection by populations of spiking neurons

Population coding is widely regarded as a key mechanism for achieving reliable behavioral decisions. We previously introduced reinforcement learning for population-based decision making by spiking neurons. Here we generalize population reinforcement learning to spike-based plasticity rules that take account of the postsynaptic neural code. We consider spike/no-spike, spike count and spike latency codes. The multi-valued and continuous-valued features in the postsynaptic code allow for a generalization of binary decision making to multi-valued decision making and continuous-valued action selection. We show that code-specific learning rules speed up learning both for the discrete classification and the continuous regression tasks. The suggested learning rules also speed up with increasing population size as opposed to standard reinforcement learning rules. Continuous action selection is further shown to explain realistic learning speeds in the Morris water maze. Finally, we introduce the concept of action perturbation as opposed to the classical weight- or node-perturbation as an exploration mechanism underlying reinforcement learning. Exploration in the action space greatly increases the speed of learning as compared to exploration in the neuron or weight space.

Learning by the Dendritic Prediction of Somatic Spiking

Recent modeling of spike-timing-dependent plasticity indicates that plasticity involves as a third factor a local dendritic potential, besides pre- and postsynaptic firing times. We present a simple compartmental neuron model together with a non-Hebbian, biologically plausible learning rule for dendritic synapses where plasticity is modulated by these three factors. In functional terms, the rule seeks to minimize discrepancies between somatic firings and a local dendritic potential. Such prediction errors can arise in our model from stochastic fluctuations as well as from synaptic input, which directly targets the soma. Depending on the nature of this direct input, our plasticity rule subserves supervised or unsupervised learning. When a reward signal modulates the learning rate, reinforcement learning results. Hence a single plasticity rule supports diverse learning paradigms.

Matching Recall and Storage in Sequence Learning with Spiking Neural Networks

Storing and recalling spiking sequences is a general problem the brain needs to solve. It is, however, unclear what type of biologically plausible learning rule is suited to learn a wide class of spatiotemporal activity patterns in a robust way. Here we consider a recurrent network of stochastic spiking neurons composed of both visible and hidden neurons. We derive a generic learning rule that is matched to the neural dynamics by minimizing an upper bound on the Kullback–Leibler divergence from the target distribution to the model distribution. The derived learning rule is consistent with spike-timing dependent plasticity in that a presynaptic spike preceding a postsynaptic spike elicits potentiation while otherwise depression emerges. Furthermore, the learning rule for synapses that target visible neurons can be matched to the recently proposed voltage-triplet rule. The learning rule for synapses that target hidden neurons is modulated by a global factor, which shares properties with astrocytes and gives rise to testable predictions.

Development of a bead-based multiplex assay for the simultaneous detection of porcine inflammation markers using xMAP technology

Commercially available assays for the simultaneous detection of multiple inflammatory and cardiac markers in porcine blood samples are currently lacking. Therefore, this study was aimed at developing a bead-based, multiplexed flow cytometric assay to simultaneously detect porcine cytokines [interleukin (IL)-1β, IL-6, IL-10, and tumor necrosis factor alpha], chemokines (IL-8 and monocyte chemotactic protein 1), growth factors [basic fibroblast growth factor (bFGF), vascular endothelial growth factor, and platelet-derived growth factor-bb], and injury markers (cardiac troponin-I) as well as complement activation markers (C5a and sC5b-9). The method was based on the Luminex xMAP technology, resulting in the assembly of a 6- and 11-plex from the respective individual singleplex situation. The assay was evaluated for dynamic range, sensitivity, cross-reactivity, intra-assay and interassay variance, spike recovery, and correlation between multiplex and commercially available enzyme-linked immunosorbent assay as well as the respective singleplex. The limit of detection ranged from 2.5 to 30,000 pg/ml for all analytes (6- and 11-plex assays), except for soluble C5b-9 with a detection range of 2-10,000 ng/ml (11-plex). Typically, very low cross-reactivity (<3% and <1.4% by 11- and 6-plex, respectively) between analytes was found. Intra-assay variances ranged from 4.9 to 7.4% (6-plex) and 5.3 to 12.9% (11-plex). Interassay variances for cytokines were between 8.1 and 28.8% (6-plex) and 10.1 and 26.4% (11-plex). Correlation coefficients with singleplex assays for 6-plex as well as for 11-plex were high, ranging from 0.988 to 0.997 and 0.913 to 0.999, respectively. In this study, a bead-based porcine 11-plex and 6-plex assay with a good assay sensitivity, broad dynamic range, and low intra-assay variance and cross-reactivity was established. These assays therefore represent a new, useful tool for the analysis of samples generated from experiments with pigs.

Petrogenesis of rift-related tephrites, phonolites and trachytes (Central European Volcanic Province, Rhon, FRG): Constraints from Sr, Nd, Pb and O isotopes

The volcanic rocks of the Rhön area (Central European Volcanic Province, Germany) belong to a moderately alkali basaltic suite that is associated with minor tephriphonolites, phonotephrites, tephrites, phonolites and trachytes. Based on isotope sytematics (87Sr/86Sr: 0.7033–0.7042; 143Nd/144Nd: 0.51279–0.51287; 206Pb/204Pb: 19.1–19.5), the inferred parental magmas formed by variable degrees of partial melting of a common asthenospheric mantle source (EAR: European Asthenospheric Reservoir of Cebriá and Wilson, 1995). Tephrites, tephriphonolites, phonotephrites, phonolites and trachytes show depletions and enrichments in some trace elements (Sr, Ba, Nb, Zr, Y) indicating that they were generated by broadly similar differentiation processes that were dominated by fractionation of olivine, clinopyroxene, amphibole, apatite and titaniferous magnetite ± plagioclase ± alkalifeldspar. The fractionated samples seem to have evolved by two distinct processes. One is characterized by pure fractional crystallization indicated by increasing Nb (and other incompatible trace element) concentrations at virtually constant 143Nd/144Nd ~ 0.51280 and 87Sr/86Sr ~ 0.7035. The other process involved an assimilation–fractional crystallization (AFC) process where moderate assimilation to crystallization rates produced evolved magmas characterized by higher Nb concentrations at slightly lower 143Nd/144Nd down to 0.51275. Literature data for some of the evolved rocks show more variable 87Sr/86Sr ranging from 0.7037 to 0.7089 at constant 143Nd/144Nd ~ 0.51280. These features may result from assimilation of upper crustal rocks by highly differentiated low-Sr (< 100 ppm Sr) lavas. However, based on the displacement of the differentiated rocks from this study towards lower 143Nd/144Nd ratios and modeled AFC processes in 143Nd/144Nd vs. 87Sr/86Sr and 207Pb/204Pb vs. 143Nd/144Nd space assimilation of lower crustal rocks seems more likely. The view that assimilation of lower crustal rocks played a role is confirmed by high-precision double-spike Pb isotope data that reveal higher 207Pb/204Pb ratios (15.62–15.63) in the differentiated rocks than in the primitive basanites (15.58–15.61). This is compatible with incorporation of radiogenic Pb from lower crustal xenoliths (207Pb/204Pb: 15.63–15.69) into the melt. However, 206Pb/204Pb ratios are similar for the differentiated rocks (19.13–19.35) and the primitive basanites (19.12–19.55) implying that assimilation involved an ancient crustal end member with a higher U/Pb ratio than the mantle source of the basanites. In addition, alteration-corrected δ18O values of the differentiated rocks range from c. 5 to 7‰ which is the same range as observed in the primitive alkaline rocks. This study confirms previous interpretations that highlighted the role of AFC processes in the evolution of alkaline volcanic rocks in the Rhön area of the Central European Volcanic Province.

Inhibition of dendritic Ca2+ spikes by GABAB receptors in cortical pyramidal neurons is mediated by a direct Gi/o- -subunit interaction with Cav1 channels

Voltage-dependent calcium channels (VDCCs) serve a wide range of physiological functions and their activity is modulated by different neurotransmitter systems. GABAergic inhibition of VDCCs in neurons has an important impact in controlling transmitter release, neuronal plasticity, gene expression and neuronal excitability. We investigated the molecular signalling mechanisms by which GABAB receptors inhibit calcium-mediated electrogenesis (Ca2+ spikes) in the distal apical dendrite of cortical layer 5 pyramidal neurons. Ca2+ spikes are the basis of coincidence detection and signal amplification of distal tuft synaptic inputs characteristic for the computational function of cortical pyramidal neurons. By combining dendritic whole-cell recordings with two-photon fluorescence Ca2+ imaging we found that all subtypes of VDCCs were present in the Ca2+ spike initiation zone, but that they contribute differently to the initiation and sustaining of dendritic Ca2+ spikes. Particularly, Cav1 VDCCs are the most abundant VDCC present in this dendritic compartment and they generated the sustained plateau potential characteristic for the Ca2+ spike. Activation of GABAB receptors specifically inhibited Cav1 channels. This inhibition of L-type Ca2+ currents was transiently relieved by strong depolarization but did not depend on protein kinase activity. Therefore, our findings suggest a novel membrane-delimited interaction of the Gi/o-βγ-subunit with Cav1 channels identifying this mechanism as the general pathway of GABAB receptor-mediated inhibition of VDCCs. Furthermore, the characterization of the contribution of the different VDCCs to the generation of the Ca2+ spike provides new insights into the molecular mechanism of dendritic computation.

Development of dendritic tonic GABAergic inhibition regulates excitability and plasticity in CA1 pyramidal neurons

Synaptic plasticity rules change during development: while hippocampal synapses can be potentiated by a single action potential pairing protocol in young neurons, mature neurons require burst firing to induce synaptic potentiation. An essential component for spike timing-dependent plasticity is the backpropagating action potential (BAP). BAP along the dendrites can be modulated by morphology and ion channel composition, both of which change during late postnatal development. However it is unclear whether these dendritic changes can explain the developmental changes in synaptic plasticity induction rules. Here, we show that tonic GABAergic inhibition regulates dendritic action potential backpropagation in adolescent but not pre-adolescent CA1 pyramidal neurons. These developmental changes in tonic inhibition also altered the induction threshold for spike timing-dependent plasticity in adolescent neurons. This GABAergic regulatory effect upon backpropagation is restricted to distal regions of apical dendrites (>200 μm) and mediated by α5-containing GABA(A) receptors. Direct dendritic recordings demonstrate α5-mediated tonic GABA(A) currents in adolescent neurons which can modulate backpropagating action potentials. These developmental modulations in dendritic excitability could not be explained by concurrent changes in dendritic morphology. To explain our data, model simulations propose a distally-increasing or localized distal expression of dendritic α5 tonic inhibition in mature neurons. Overall, our results demonstrate that dendritic integration and plasticity in more mature dendrites are significantly altered by tonic α5 inhibition in a dendritic region-specific and developmentally-regulated manner.

Dynamic directed interictal connectivity in left and right temporal lobe epilepsy.

OBJECTIVE There is increasing evidence that epileptic activity involves widespread brain networks rather than single sources and that these networks contribute to interictal brain dysfunction. We investigated the fast-varying behavior of epileptic networks during interictal spikes in right and left temporal lobe epilepsy (RTLE and LTLE) at a whole-brain scale using directed connectivity. METHODS In 16 patients, 8 with LTLE and 8 with RTLE, we estimated the electrical source activity in 82 cortical regions of interest (ROIs) using high-density electroencephalography (EEG), individual head models, and a distributed linear inverse solution. A multivariate, time-varying, and frequency-resolved Granger-causal modeling (weighted Partial Directed Coherence) was applied to the source signal of all ROIs. A nonparametric statistical test assessed differences between spike and baseline epochs. Connectivity results between RTLE and LTLE were compared between RTLE and LTLE and with neuropsychological impairments. RESULTS Ipsilateral anterior temporal structures were identified as key drivers for both groups, concordant with the epileptogenic zone estimated invasively. We observed an increase in outflow from the key driver already before the spike. There were also important temporal and extratemporal ipsilateral drivers in both conditions, and contralateral only in RTLE. A different network pattern between LTLE and RTLE was found: in RTLE there was a much more prominent ipsilateral to contralateral pattern than in LTLE. Half of the RTLE patients but none of the LTLE patients had neuropsychological deficits consistent with contralateral temporal lobe dysfunction, suggesting a relationship between connectivity changes and cognitive deficits. SIGNIFICANCE The different patterns of time-varying connectivity in LTLE and RTLE suggest that they are not symmetrical entities, in line with our neuropsychological results. The highest outflow region was concordant with invasive validation of the epileptogenic zone. This enhanced characterization of dynamic connectivity patterns could better explain cognitive deficits and help the management of epilepsy surgery candidates.

Random Sampling with Interspike-Intervals of the Exponential Integrate and Fire Neuron: A Computational Interpretation of UP-States

Oscillations between high and low values of the membrane potential (UP and DOWN states respectively) are an ubiquitous feature of cortical neurons during slow wave sleep and anesthesia. Nevertheless, a surprisingly small number of quantitative studies have been conducted only that deal with this phenomenon’s implications for computation. Here we present a novel theory that explains on a detailed mathematical level the computational benefits of UP states. The theory is based on random sampling by means of interspike intervals (ISIs) of the exponential integrate and fire (EIF) model neuron, such that each spike is considered a sample, whose analog value corresponds to the spike’s preceding ISI. As we show, the EIF’s exponential sodium current, that kicks in when balancing a noisy membrane potential around values close to the firing threshold, leads to a particularly simple, approximative relationship between the neuron’s ISI distribution and input current. Approximation quality depends on the frequency spectrum of the current and is improved upon increasing the voltage baseline towards threshold. Thus, the conceptually simpler leaky integrate and fire neuron that is missing such an additional current boost performs consistently worse than the EIF and does not improve when voltage baseline is increased. For the EIF in contrast, the presented mechanism is particularly effective in the high-conductance regime, which is a hallmark feature of UP-states. Our theoretical results are confirmed by accompanying simulations, which were conducted for input currents of varying spectral composition. Moreover, we provide analytical estimations of the range of ISI distributions the EIF neuron can sample from at a given approximation level. Such samples may be considered by any algorithmic procedure that is based on random sampling, such as Markov Chain Monte Carlo or message-passing methods. Finally, we explain how spike-based random sampling relates to existing computational theories about UP states during slow wave sleep and present possible extensions of the model in the context of spike-frequency adaptation.

Komatiites constrain molybdenum isotope composition of the Earth's mantle

In order to estimate the Mo isotope composition and Mo abundance in the Bulk Silicate Earth (BSE), a total of thirty komatiite samples from five localities on three continents were analyzed using an isotope dilution double spike technique. Calculated Mo concentrations of the emplaced komatiite lavas range from 25±325±3 to 66±22 ng/g66±22 ng/g, and the inferred Mo concentrations in the deep mantle sources of the komatiites range between 17±417±4 and 30±12 ng/g30±12 ng/g, with an average value of 23±7 ng/g23±7 ng/g (2SE). This average value represents our best estimate for the Mo concentration in the BSE; it is identical, within the uncertainty, to published previous estimates of 39±16 ng/g39±16 ng/g, but is at least a factor of 2 more precise. The Mo isotope compositions of the komatiite mantle sources overlap within uncertainty and range from View the MathML sourceδMo98=−0.04±0.28 to 0.11±0.10‰0.11±0.10‰, with an average of 0.04±0.06‰0.04±0.06‰ (2SE). This value is analytically indistinguishable from published Mo isotope compositions of ordinary and enstatite chondrites and represents the best estimate for the Mo isotope composition of the BSE. The inferred δ98Mo for the BSE is therefore lighter than the suggested average of the upper continental crust (0.3 to 0.4‰). Thus, from the mass balance standpoint, a reservoir with lighter Mo isotope composition should exist in the Earth's mantle; this reservoir can potentially be found in subducted oceanic crust. The similarity of δ98Mo between chondritic meteorites and estimates for the BSE from this study indicates that during the last major equilibration between Earth's core and mantle, i.e., the one that occurred during the giant impact that produced the Moon, chemical and isotopic equilibrium of Mo between Fe metal of the core and the silicate mantle was largely achieved.

Evidence for an Ancestral Association of Human Coronavirus 229E with Bats.

UNLABELLED We previously showed that close relatives of human coronavirus 229E (HCoV-229E) exist in African bats. The small sample and limited genomic characterizations have prevented further analyses so far. Here, we tested 2,087 fecal specimens from 11 bat species sampled in Ghana for HCoV-229E-related viruses by reverse transcription-PCR (RT-PCR). Only hipposiderid bats tested positive. To compare the genetic diversity of bat viruses and HCoV-229E, we tested historical isolates and diagnostic specimens sampled globally over 10 years. Bat viruses were 5- and 6-fold more diversified than HCoV-229E in the RNA-dependent RNA polymerase (RdRp) and spike genes. In phylogenetic analyses, HCoV-229E strains were monophyletic and not intermixed with animal viruses. Bat viruses formed three large clades in close and more distant sister relationships. A recently described 229E-related alpaca virus occupied an intermediate phylogenetic position between bat and human viruses. According to taxonomic criteria, human, alpaca, and bat viruses form a single CoV species showing evidence for multiple recombination events. HCoV-229E and the alpaca virus showed a major deletion in the spike S1 region compared to all bat viruses. Analyses of four full genomes from 229E-related bat CoVs revealed an eighth open reading frame (ORF8) located at the genomic 3' end. ORF8 also existed in the 229E-related alpaca virus. Reanalysis of HCoV-229E sequences showed a conserved transcription regulatory sequence preceding remnants of this ORF, suggesting its loss after acquisition of a 229E-related CoV by humans. These data suggested an evolutionary origin of 229E-related CoVs in hipposiderid bats, hypothetically with camelids as intermediate hosts preceding the establishment of HCoV-229E. IMPORTANCE The ancestral origins of major human coronaviruses (HCoVs) likely involve bat hosts. Here, we provide conclusive genetic evidence for an evolutionary origin of the common cold virus HCoV-229E in hipposiderid bats by analyzing a large sample of African bats and characterizing several bat viruses on a full-genome level. Our evolutionary analyses show that animal and human viruses are genetically closely related, can exchange genetic material, and form a single viral species. We show that the putative host switches leading to the formation of HCoV-229E were accompanied by major genomic changes, including deletions in the viral spike glycoprotein gene and loss of an open reading frame. We reanalyze a previously described genetically related alpaca virus and discuss the role of camelids as potential intermediate hosts between bat and human viruses. The evolutionary history of HCoV-229E likely shares important characteristics with that of the recently emerged highly pathogenic Middle East respiratory syndrome (MERS) coronavirus.

Selenium isotope analysis by N-TIMS: Potential and challenges

The isotope composition of selenium (Se) can provide important constraints on biological, geochemical, and cosmochemical processes taking place in different reservoirs on Earth and during planet formation. To provide precise qualitative and quantitative information on these processes, accurate and highly precise isotope data need to be obtained. The currently applied ICP-MS methods for Se isotope measurements are compromised by the necessity to perform a large number of interference corrections. Differences in these correction methods can lead to discrepancies in published Se isotope values of rock standards which are significantly higher than the acclaimed precision. An independent analytical approach applying a double spike (DS) and state-of-the-art TIMS may yield better precision due to its smaller number of interferences and could test the accuracy of data obtained by ICP-MS approaches. This study shows that the precision of Se isotope measurements performed with two different Thermo Scientific™ Triton™ Plus TIMS is distinctly deteriorated by about ±1‰ (2 s.d.) due to δ80/78Se by a memory Se signal of up to several millivolts and additional minor residual mass bias which could not be corrected for with the common isotope fractionation laws. This memory Se has a variable isotope composition with a DS fraction of up to 20% and accumulates with increasing number of measurements. Thus it represents an accumulation of Se from previous Se measurements with a potential addition from a sample or machine blank. Several cleaning techniques of the MS parts were tried to decrease the memory signal, but were not sufficient to perform precise Se isotope analysis. If these serious memory problems can be overcome in the future, the precision and accuracy of Se isotope analysis with TIMS should be significantly better than those of the current ICP-MS approaches.

Influence of sample matrix on the alkaline extraction of Cr(VI) in soils and industrial materials

An accurate and efficient determination of the highly toxic Cr(VI) in solid materials is important to determine the total Cr(VI) inventory of contaminated sites and the Cr(VI) release potential from such sites into the environment. Most commonly, total Cr(VI) is extracted from solid materials following a hot alkaline extraction procedure (US EPA method 3060A) where a complete release of water-extractable and sparingly soluble Cr(VI) phase is achieved. This work presents an evaluation of matrix effects that may occur during the hot alkaline extraction and in the determination of the total Cr(VI) inventory of variably composed contaminated soils and industrial materials (cement, fly ash) and is compared to water-extractable Cr(VI) results. Method validation including multiple extractions and matrix spiking along with chemical and mineralogical characterization showed satisfying results for total Cr(VI) contents for most of the tested materials. However, unreliable results were obtained by applying method 3060A to anoxic soils due to the degradation of organic material and/or reactions with Fe2+-bearing mineral phases. In addition, in certain samples discrepant spike recoveries have to be also attributed to sample heterogeneity. Separation of possible extracted Cr(III) by applying cation-exchange cartridges prior to solution analysis further shows that under the hot alkaline extraction conditions only Cr(VI) is present in solution in measurable amounts, whereas Cr(III) gets precipitated as amorphous Cr(OH)3(am). It is concluded that prior to routine application of method 3060A to a new material type, spiking tests are recommended for the identification of matrix effects. In addition, the mass of extracted solid material should to be well adjusted to the heterogeneity of the Cr(VI) distribution in the material in question.