Ready…Go: Amplitude of the fMRI Signal Encodes Expectation of Cue Arrival Time

What happens when the brain awaits a signal of uncertain arrival time, as when a sprinter waits for the starting pistol? And what happens just after the starting pistol fires? Using functional magnetic resonance imaging (fMRI), we have discovered a novel correlate of temporal expectations in several brain regions, most prominently in the supplementary motor area (SMA). Contrary to expectations, we found little fMRI activity during the waiting period; however, a large signal appears after the “go” signal, the amplitude of which reflects learned expectations about the distribution of possible waiting times. Specifically, the amplitude of the fMRI signal appears to encode a cumulative conditional probability, also known as the cumulative hazard function. The fMRI signal loses its dependence on waiting time in a “countdown” condition in which the arrival time of the go cue is known in advance, suggesting that the signal encodes temporal probabilities rather than simply elapsed time. The dependence of the signal on temporal expectation is present in “no-go” conditions, demonstrating that the effect is not a consequence of motor output. Finally, the encoding is not dependent on modality, operating in the same manner with auditory or visual signals. This finding extends our understanding of the relationship between temporal expectancy and measurable neural signals.     

Ready…Go: Amplitude of the fMRI Signal Encodes Expectation of Cue Arrival Time

What happens when the brain awaits a signal of uncertain arrival time, as when a sprinter waits for the starting pistol? And what happens just after the starting pistol fires? Using functional magnetic resonance imaging (fMRI), we have discovered a novel correlate of temporal expectations in several brain regions, most prominently in the supplementary motor area (SMA). Contrary to expectations, we found little fMRI activity during the waiting period; however, a large signal appears after the “go” signal, the amplitude of which reflects learned expectations about the distribution of possible waiting times. Specifically, the amplitude of the fMRI signal appears to encode a cumulative conditional probability, also known as the cumulative hazard function. The fMRI signal loses its dependence on waiting time in a “countdown” condition in which the arrival time of the go cue is known in advance, suggesting that the signal encodes temporal probabilities rather than simply elapsed time. The dependence of the signal on temporal expectation is present in “no-go” conditions, demonstrating that the effect is not a consequence of motor output. Finally, the encoding is not dependent on modality, operating in the same manner with auditory or visual signals. This finding extends our understanding of the relationship between temporal expectancy and measurable neural signals.     

A Neural Model of Auditory Space Compatible with Human Perception under Simulated Echoic Conditions

In a typical auditory scene, sounds from different sources and reflective surfaces summate in the ears, causing spatial cues to fluctuate. Prevailing hypotheses of how spatial locations may be encoded and represented across auditory neurons generally disregard these fluctuations and must therefore invoke additional mechanisms for detecting and representing them. Here, we consider a different hypothesis in which spatial perception corresponds to an intermediate or sub-maximal firing probability across spatially selective neurons within each hemisphere. The precedence or Haas effect presents an ideal opportunity for examining this hypothesis, since the temporal superposition of an acoustical reflection with sounds arriving directly from a source can cause otherwise stable cues to fluctuate. Our findings suggest that subjects’ experiences may simply reflect the spatial cues that momentarily arise under various acoustical conditions and how these cues are represented. We further suggest that auditory objects may acquire “edges” under conditions when interaural time differences are broadly distributed.     

Leveraging Genomic Annotations and Pleiotropic Enrichment for Improved Replication Rates in Schizophrenia GWAS

Most of the genetic architecture of schizophrenia (SCZ) has not yet been identified. Here, we apply a novel statistical algorithm called Covariate-Modulated Mixture Modeling (CM3), which incorporates auxiliary information (heterozygosity, total linkage disequilibrium, genomic annotations, pleiotropy) for each single nucleotide polymorphism (SNP) to enable more accurate estimation of replication probabilities, conditional on the observed test statistic (“z-score”) of the SNP. We use a multiple logistic regression on z-scores to combine information from auxiliary information to derive a “relative enrichment score” for each SNP. For each stratum of these relative enrichment scores, we obtain nonparametric estimates of posterior expected test statistics and replication probabilities as a function of discovery z-scores, using a resampling-based approach that repeatedly and randomly partitions meta-analysis sub-studies into training and replication samples. We fit a scale mixture of two Gaussians model to each stratum, obtaining parameter estimates that minimize the sum of squared differences of the scale-mixture model with the stratified nonparametric estimates. We apply this approach to the recent genome-wide association study (GWAS) of SCZ (n = 82,315), obtaining a good fit between the model-based and observed effect sizes and replication probabilities. We observed that SNPs with low enrichment scores replicate with a lower probability than SNPs with high enrichment scores even when both they are genome-wide significant (p < 5x10^-8). There were 693 and 219 independent loci with model-based replication rates ≥80% and ≥90%, respectively. Compared to analyses not incorporating relative enrichment scores, CM3 increased out-of-sample yield for SNPs that replicate at a given rate. This demonstrates that replication probabilities can be more accurately estimated using prior enrichment information with CM3.     

On the Transmission of Information through Sensory Neurons

Information about muscle length is transmitted to the cerebellum from muscle spindle receptors through the dorsal spinocerebellar tract (DSCT). The “transinformation” about muscle length in single DSCT fibers was calculated from steady-state spike trains by two different methods, assuming that the decoding mechanisms use a frequency code. By the first method, the number of distinguishable muscle lengths (and thus the transiformation) was determined from the rate of convergence of the mean frequency of firing (with increasing number of intervals). The observation time necessary to estimate the mean frequency of the impulse train with a certain accuracy was independent of the stretch level, even though the number of intervals necessary to make this estimate was different at high and low levels of stretch. By the second method an input frequency-output frequency matrix was calculated. The transinformations and the rate of transinformation was then calculated from this matrix. There was an acceptable agreement in the estimates of transinformation by the two methods. The rates of transinformation are significantly increased by the particular time structure of the discharge patterns of the nerve cells. Consequently, the loss of information due to the synaptic coupling is appreciably reduced.     

The Impact of Various Time Intervals on the Supragingival Plaque Dynamic Core Microbiome

ObjectiveThe aim of this study was to examine the influence of various time intervals on the composition of the supragingival plaque microbiome, especially the dynamic core microbiome, and to find a suitable observation interval for further studies on oral microbiota.ResultsA total of 359,565 qualified reads for 64 samples were generated for subsequent analyses, which represents 8,452 operational taxonomic units identified at 3% dissimilarity. The dynamic core microbiome detected in the current study included five phyla, 12 genera and 13 species. At the genus level, the relative abundance of bacterial communities under the “1 day,” “1 month,” and “3 months” intervals was clustered into sub-category. At the species level, the number of overlapping species remained stable between the “1 month” and “3 months” intervals, whereas the number of dynamic core species became stable within only 1 week.ConclusionsThis study emphasized the impact of different time intervals (days, weeks and months) on the composition, commonality and diversity of the supragingival microbiome. The analyses found that for various types of studies, the time interval of a month is more suitable for observing the general composition of the supragingival microbiome, and that a week is better for observing the dynamic core microbiome.     

Auditory Nerve Spike Generator Modeled as a Variable Attenuator Based on a Saddle Node on Invariant Circle Bifurcation

Mammalian inner hair cells transduce the sound waves amplified by the cochlear amplifier (CA) into a graded neurotransmitter release that activates channels on auditory nerve fibers (ANF). These synaptic channels then charge its dendritic spike generator. While the outer hair cells of the CA employ positive feedback, poising on Andronov-Hopf type instabilities which make them extremely sensitive to faint sounds and make CA output strongly nonlinear, the ANF appears to be based on different principles and a different type of dynamical instability. Its spike generator “digitizes” CA output into trains of action potentials and behaves as a linear filter, rate-coding sound intensity across a wide dynamic range. Here we model the spike generator as a 3 dimensional version of a saddle node on invariant circle (SNIC) bifurcation. The generic 2d SNIC increases its spike rate as the square root of the input current above its spiking threshold. We add negative feedback in the form of a low voltage-threshold potassium conductance that slows down the generator’s rate of increase of its spike rate. A Poisson random source simulates an inner hair cell, outputting a series of noisy periodic current pulses to the model ANF whose spikes phase lock to these pulses and have a linear frequency to current relation with a wide dynamic range. Also, the spike generator compartment has a cholinergic feedback connection from the olive and experiments show that such feedback is able to alter the amount of H conductance inside the generator compartment. We show that an olive able to decrease H would be able to shift the spike generator’s dynamic range to higher sound intensities. In a quiet environment by increasing H the olive would be able to make spike trains similar to those caused by synaptic input.     

Directionality of auditory nerve fiber responses to pure tone stimuli in the grassfrog, Rana temporaria. II. Spike timing

We studied the directionality of spike timing in the responses of single auditory nerve fibers of the grass frog, Rana temporaria, to tone burst stimulation. Both the latency of the first spike after stimulus onset and the preferred firing phase during the stimulus were studied. In addition, the directionality of the phase of eardrum vibrations was measured. The response latency showed systematic and statistically significant changes with sound direction at both low and high frequencies. The latency changes were correlated with response strength (spike rate) changes and were probably the result of directional changes in effective stimulus intensity. Systematic changes in the preferred firing phase were seen in all fibers that showed phaselocking (i.e., at frequencies below 500–700 Hz). The mean phase lead for stimulation from the contralateral side was approximately 140° at 200 Hz and decreased to approximately 100° at 700 Hz. These phaseshifts correspond to differences in spike timing of approximately 2 ms and 0.4 ms respectively. The phaseshifts were nearly independent of stimulus intensity. The phase directionality of eardrum vibrations was smaller than that of the nerve fibers. Hence, the strong directional phaseshifts shown by the nerve fibers probably reflect the directional characteristics of extratympanic pathways. Accepted: 23 November 1996     

Properties of synaptic noise in tonically active human motoneurons

The objective of these experiments was to determine the amount of synaptic noise on the cell membrane at various intervals after an action potential in a motoneuron firing at a specified frequency. Sources of noise such as variations in the level of voluntary drive were minimized by selecting only segments of the spike train in which the unit was running within prescribed frequency limits. The level of the membrane potential of the motoneuron during these intervals was determined using two test “pulses” (compound Ia excitatory postsynaptic potentials) of known amplitude. This enabled the probability of the membrane potential falling within a voltage “window” of known size at known times after the preceding spike to be determined. The probability density histograms showed that the fluctuations of membrane potential about a target interspike trajectory (i.e., the membrane noise) increased with time after the preceding spike. These fluctuations in the membrane potential can be accounted for by a one-dimensional “random walk” model of membrane noise. This model explains the salient features of the interval histograms, such as positive skewness at low target frequencies. A quantitative test of the model demonstrated its applicability to the motor pools of tibialis and masseter.     

Proposed Classification of Auriculotemporal Nerve,Based on the Root System

The topography of the auriculotemporal nerve (ATN) root system is the main criterion of this nerve classification. Previous publications indicate that ATN may have between one and five roots. Most common is a one- or two-root variant of the nerve structure. The problem of many publications is the inconsistency of nomenclature which concerns the terms “roots”, “connecting branches”, or “branches” that are used to identify the same structures. This study was performed on 80 specimens (40 adults and 40 fetuses) to propose a classification based on: (i) the number of roots, (ii) way of root division, and (iii) configuration of interradicular fibers that form the ATN trunk. This new classification is a remedy for inconsistency of nomenclature of ATN in the infratemporal fossa. This classification system has proven beneficial when organizing all ATN variants described in previous studies and could become a helpful tool for surgeons and dentists. Examination of ATN from the infratemporal fossa of fetuses (the youngest was at 18 weeks gestational age) showed that, at that stage, the nerve is fully developed.     

The Donnan-dominated resting state of skeletal muscle fibers contributes to resilience and longevity in dystrophic fibers

Duchenne muscular dystrophy (DMD) is an X-linked dystrophin-minus muscle-wasting disease. Ion homeostasis in skeletal muscle fibers underperforms as DMD progresses. But though DMD renders these excitable cells intolerant of exertion, sodium overloaded, depolarized, and spontaneously contractile, they can survive for several decades. We show computationally that underpinning this longevity is a strikingly frugal, robust Pump-Leak/Donnan (P-L/D) ion homeostatic process. Unlike neurons, which operate with a costly “Pump-Leak–dominated” ion homeostatic steady state, skeletal muscle fibers operate with a low-cost “Donnan-dominated” ion homeostatic steady state that combines a large chloride permeability with an exceptionally small sodium permeability. Simultaneously, this combination keeps fiber excitability low and minimizes pump expenditures. As mechanically active, long-lived multinucleate cells, skeletal muscle fibers have evolved to handle overexertion, sarcolemmal tears, ischemic bouts, etc.; the frugality of their Donnan dominated steady state lets them maintain the outsized pump reserves that make them resilient during these inevitable transient emergencies. Here, P-L/D model variants challenged with DMD-type insult/injury (low pump-strength, overstimulation, leaky Nav and cation channels) show how chronic “nonosmotic” sodium overload (observed in DMD patients) develops. Profoundly severe DMD ion homeostatic insult/injury causes spontaneous firing (and, consequently, unwanted excitation–contraction coupling) that elicits cytotoxic swelling. Therefore, boosting operational pump-strength and/or diminishing sodium and cation channel leaks should help extend DMD fiber longevity.     

Categorical Vowel Perception Enhances the Effectiveness and Generalization of Auditory Feedback in Human-Machine-Interfaces

Human-machine interface (HMI) designs offer the possibility of improving quality of life for patient populations as well as augmenting normal user function. Despite pragmatic benefits, utilizing auditory feedback for HMI control remains underutilized, in part due to observed limitations in effectiveness. The goal of this study was to determine the extent to which categorical speech perception could be used to improve an auditory HMI. Using surface electromyography, 24 healthy speakers of American English participated in 4 sessions to learn to control an HMI using auditory feedback (provided via vowel synthesis). Participants trained on 3 targets in sessions 1–3 and were tested on 3 novel targets in session 4. An “established categories with text cues” group of eight participants were trained and tested on auditory targets corresponding to standard American English vowels using auditory and text target cues. An “established categories without text cues” group of eight participants were trained and tested on the same targets using only auditory cuing of target vowel identity. A “new categories” group of eight participants were trained and tested on targets that corresponded to vowel-like sounds not part of American English. Analyses of user performance revealed significant effects of session and group (established categories groups and the new categories group), and a trend for an interaction between session and group. Results suggest that auditory feedback can be effectively used for HMI operation when paired with established categorical (native vowel) targets with an unambiguous cue.     

Auditory Stimuli Mimicking Ambient Sounds Drive Temporal “Delta-Brushes” in Premature Infants

In the premature infant, somatosensory and visual stimuli trigger an immature electroencephalographic (EEG) pattern, “delta-brushes,” in the corresponding sensory cortical areas. Whether auditory stimuli evoke delta-brushes in the premature auditory cortex has not been reported. Here, responses to auditory stimuli were studied in 46 premature infants without neurologic risk aged 31 to 38 postmenstrual weeks (PMW) during routine EEG recording. Stimuli consisted of either low-volume technogenic “clicks” near the background noise level of the neonatal care unit, or a human voice at conversational sound level. Stimuli were administrated pseudo-randomly during quiet and active sleep. In another protocol, the cortical response to a composite stimulus (“click” and voice) was manually triggered during EEG hypoactive periods of quiet sleep. Cortical responses were analyzed by event detection, power frequency analysis and stimulus locked averaging. Before 34 PMW, both voice and “click” stimuli evoked cortical responses with similar frequency-power topographic characteristics, namely a temporal negative slow-wave and rapid oscillations similar to spontaneous delta-brushes. Responses to composite stimuli also showed a maximal frequency-power increase in temporal areas before 35 PMW. From 34 PMW the topography of responses in quiet sleep was different for “click” and voice stimuli: responses to “clicks” became diffuse but responses to voice remained limited to temporal areas. After the age of 35 PMW auditory evoked delta-brushes progressively disappeared and were replaced by a low amplitude response in the same location. Our data show that auditory stimuli mimicking ambient sounds efficiently evoke delta-brushes in temporal areas in the premature infant before 35 PMW. Along with findings in other sensory modalities (visual and somatosensory), these findings suggest that sensory driven delta-brushes represent a ubiquitous feature of the human sensory cortex during fetal stages and provide a potential test of functional cortical maturation during fetal development.     

From decision to action: Detailed modelling of frog tadpoles reveals neuronal mechanisms of decision-making and reproduces unpredictable swimming movements in response to sensory signals

How does the brain process sensory stimuli, and decide whether to initiate locomotor behaviour? To investigate this question we develop two whole body computer models of a tadpole. The “Central Nervous System” (CNS) model uses evidence from whole-cell recording to define 2300 neurons in 12 classes to study how sensory signals from the skin initiate and stop swimming. In response to skin stimulation, it generates realistic sensory pathway spiking and shows how hindbrain sensory memory populations on each side can compete to initiate reticulospinal neuron firing and start swimming. The 3-D “Virtual Tadpole” (VT) biomechanical model with realistic muscle innervation, body flexion, body-water interaction, and movement is then used to evaluate if motor nerve outputs from the CNS model can produce swimming-like movements in a volume of “water”. We find that the whole tadpole VT model generates reliable and realistic swimming. Combining these two models opens new perspectives for experiments.     

Phantoms in the brain: Ambiguous representations of stimulus amplitude and timing in weakly electric fish

In wave-type weakly electric fish, two distinct types of primary afferent fibers are specialized for separately encoding modulations in the amplitude and phase (timing) of electrosensory stimuli. Time-coding afferents phase lock to periodic stimuli and respond to changes in stimulus phase with shifts in spike timing. Amplitude-coding afferents fire sporadically to periodic stimuli. Their probability of firing in a given cycle, and therefore their firing rate, is proportional to stimulus amplitude. However, the spike times of time-coding afferents are also affected by changes in amplitude; similarly, the firing rates of amplitude-coding afferents are also affected by changes in phase. Because identical changes in the activity of an individual primary afferent can be caused by modulations in either the amplitude or phase of stimuli, there is ambiguity regarding the information content of primary afferent responses that can result in ‘phantom’ modulations not present in an actual stimulus. Central electrosensory neurons in the hindbrain and midbrain respond to these phantom modulations. Phantom modulations can also elicit behavioral responses, indicating that ambiguity in the encoding of amplitude and timing information ultimately distorts electrosensory perception. A lack of independence in the encoding of multiple stimulus attributes can therefore result in perceptual illusions. Similar effects may occur in other sensory systems as well. In particular, the vertebrate auditory system is thought to be phylogenetically related to the electrosensory system and it encodes information about amplitude and timing in similar ways. It has been well established that pitch perception and loudness perception are both affected by the frequency and intensity of sounds, raising the intriguing possibility that auditory perception may also be affected by ambiguity in the encoding of sound amplitude and timing.     

Bortezomib Resistance Can Be Reversed by Induced Expression of Plasma Cell Maturation Markers in a Mouse In Vitro Model of Multiple Myeloma

Multiple myeloma (MM), the second most common hematopoietic malignancy, remains an incurable plasma cell (PC) neoplasm. While the proteasome inhibitor, bortezomib (Bz) has increased patient survival, resistance represents a major treatment obstacle as most patients ultimately relapse becoming refractory to additional Bz therapy. Current tests fail to detect emerging resistance; by the time patients acquire resistance, their prognosis is often poor. To establish immunophenotypic signatures that predict Bz sensitivity, we utilized Bz-sensitive and -resistant cell lines derived from tumors of the Bcl-X_L/Myc mouse model of PC malignancy. We identified significantly reduced expression of two markers (CD93, CD69) in “acquired” (Bz-selected) resistant cells. Using this phenotypic signature, we isolated a subpopulation of cells from a drug-naïve, Bz-sensitive culture that displayed “innate” resistance to Bz. Although these genes were identified as biomarkers, they may indicate a mechanism for Bz-resistance through the loss of PC maturation which may be induced and/or selected by Bz. Significantly, induction of PC maturation in both “acquired” and “innate” resistant cells restored Bz sensitivity suggesting a novel therapeutic approach for reversing Bz resistance in refractory MM.     

Audio-Visual Perception of 3D Cinematography: An fMRI Study Using Condition-Based and Computation-Based Analyses

The use of naturalistic stimuli to probe sensory functions in the human brain is gaining increasing interest. Previous imaging studies examined brain activity associated with the processing of cinematographic material using both standard “condition-based” designs, as well as “computational” methods based on the extraction of time-varying features of the stimuli (e.g. motion). Here, we exploited both approaches to investigate the neural correlates of complex visual and auditory spatial signals in cinematography. In the first experiment, the participants watched a piece of a commercial movie presented in four blocked conditions: 3D vision with surround sounds (3D-Surround), 3D with monaural sound (3D-Mono), 2D-Surround, and 2D-Mono. In the second experiment, they watched two different segments of the movie both presented continuously in 3D-Surround. The blocked presentation served for standard condition-based analyses, while all datasets were submitted to computation-based analyses. The latter assessed where activity co-varied with visual disparity signals and the complexity of auditory multi-sources signals. The blocked analyses associated 3D viewing with the activation of the dorsal and lateral occipital cortex and superior parietal lobule, while the surround sounds activated the superior and middle temporal gyri (S/MTG). The computation-based analyses revealed the effects of absolute disparity in dorsal occipital and posterior parietal cortices and of disparity gradients in the posterior middle temporal gyrus plus the inferior frontal gyrus. The complexity of the surround sounds was associated with activity in specific sub-regions of S/MTG, even after accounting for changes of sound intensity. These results demonstrate that the processing of naturalistic audio-visual signals entails an extensive set of visual and auditory areas, and that computation-based analyses can track the contribution of complex spatial aspects characterizing such life-like stimuli.     

The Time Course of Corticospinal Excitability during a Simple Reaction Time Task

The production of movement in a simple reaction time task can be separated into two time periods: the foreperiod, which is thought to include preparatory processes, and the reaction time interval, which includes initiation processes. To better understand these processes, transcranial magnetic stimulation has been used to probe corticospinal excitability at various time points during response preparation and initiation. Previous research has shown that excitability decreases prior to the “go” stimulus and increases following the “go”; however these two time frames have been examined independently. The purpose of this study was to measure changes in CE during both the foreperiod and reaction time interval in a single experiment, relative to a resting baseline level. Participants performed a button press movement in a simple reaction time task and excitability was measured during rest, the foreperiod, and the reaction time interval. Results indicated that during the foreperiod, excitability levels quickly increased from baseline with the presentation of the warning signal, followed by a period of stable excitability leading up to the “go” signal, and finally a rapid increase in excitability during the reaction time interval. This excitability time course is consistent with neural activation models that describe movement preparation and response initiation.     

STRETCHING OF THE SCIATIC NERVE—A Means of Relieving Postoperative Pain Following Removal of Ruptured Lumbar Intervertebral Discs

Stretching the sciatic nerve for the relief of “sciatica” was frequently employed before 1900 and was subsequently abandoned, probably because it was done without sufficient scrutiny of the indications. The procedure has recently been employed in cases in which “sciatica” remains following the operative removal of ruptured intervertebral discs, and it has been instrumental in relieving postoperative “sciatica” when the cause was the formation of adhesions about the lumbar nerve roots. If the nerve root is compressed by recurrent disc protrusion or by adjacent bone, the manipulation usually increases the pain, a phenomenon that has been helpful from a diagnostic standpoint.     

A cross-interval spike train analysis: the correlation between spike generation and temporal integration of doublets

A stochastic spike train analysis technique is introduced to reveal the correlation between the firing of the next spike and the temporal integration period of two consecutive spikes (i.e., a doublet). Statistics of spike firing times between neurons are established to obtain the conditional probability of spike firing in relation to the integration period. The existence of a temporal integration period is deduced from the time interval between two consecutive spikes fired in a reference neuron as a precondition to the generation of the next spike in a compared neuron. This analysis can show whether the coupled spike firing in the compared neuron is correlated with the last or the second-to-last spike in the reference neuron. Analysis of simulated and experimentally recorded biological spike trains shows that the effects of excitatory and inhibitory temporal integration are extracted by this method without relying on any subthreshold potential recordings. The analysis also shows that, with temporal integration, a neuron driven by random firing patterns can produce fairly regular firing patterns under appropriate conditions. This regularity in firing can be enhanced by temporal integration of spikes in a chain of polysynaptically connected neurons. The bandpass filtering of spike firings by temporal integration is discussed. The results also reveal that signal transmission delays may be attributed not just to conduction and synaptic delays, but also to the delay time needed for temporal integration. Received: 3 March 1997 / Accepted in revised form: 6 November 1997