Growth curve models of repeated binary response

Experimental designs that include repeated measures of binary response variables over time and under different conditions are common in biology. In such settings, it is often desirable to characterize the response pattern over time. When response variables are continuous, this characterization can be made in terms of a growth model such as the Potthoff-Roy growth curve model. We illustrate how a similar growth curve modeling strategy can be implemented using weighted least squares (WLS) methods for binary response data. The growth models are constructed in terms of polynomial functions across marginal response. However, when growth models are fit to repeated binary response, the nonsignificant higher-order polynomial functions are dropped from the model, rather than used as covariates. Dropping the nonsignificant polynomials from the model will reduce the number of response functions, and help avoid small-sample problems that can occur when the number of correlated response functions is large and sample sizes are small. The reduced set of response functions are then modeled using WLS methods. We illustrate such models with an example of binary fly oviposition response (accept or reject) exhibited by two populations of flies at four ages to two types of fruit.     

Heterogeneous response of microvascular endothelial cells to shear stress

We investigated changes in calcium concentration in cultured bovine aortic endothelial cells (BAECs) and rat adrenomedulary endothelial cells (RAMECs, microvascular) in response to different levels of shear stress. In BAECs, the onset of shear stress elicited a transient increase in intracellular calcium concentration that was spatially uniform, synchronous, and dose dependent. In contrast, the response of RAMECs was heterogeneous in time and space. Shear stress induced calcium waves that originated from one or several cells and propagated to neighboring cells. The number and size of the responding groups of cells did not depend on the magnitude of shear stress or the magnitude of the calcium change in the responding cells. The initiation and the propagation of calcium waves in RAMECs were significantly suppressed under conditions in which either purinergic receptors were blocked by suramin or extracellular ATP was degraded by apyrase. Exogenously applied ATP produced similarly heterogeneous responses. The number of responding cells was dependent on ATP concentration, but the magnitude of the calcium change was not. Our data suggest that shear stress stimulates RAMECs to release ATP, causing the increase in intracellular calcium concentration via purinergic receptors in cells that are heterogeneously sensitive to ATP. The propagation of the calcium signal is also mediated by ATP, and the spatial pattern suggests a locally elevated ATP concentration in the vicinity of the initially responding cells.     

Separating spatially distributed response to stimulation from background. I. Optical imaging

We consider the problem of estimating a small stimulus-induced response to stimulation that is masked by a fluctuating background when measurements of the background in the absence of stimulation are available, as is common in optical imaging of the cortex and in many other experimental situations. Two related methods based on the Karhunen-Loève procedure are discussed. One seeks the function, an indicator function, that is most parallel to the response data and most orthogonal to the background data. The second removes the subspace spanned by the background from the response. Numerical investigations on simulated optical imaging data show that the first method is generally superior. Connections between the two methods and techniques for assessing the quality of the result are discussed. Received: 17 January 1997 / Accepted in revised form: 4 August 1997     

Memory Maintenance in Synapses with Calcium-Based Plasticity in the Presence of Background Activity

Most models of learning and memory assume that memories are maintained in neuronal circuits by persistent synaptic modifications induced by specific patterns of pre- and postsynaptic activity. For this scenario to be viable, synaptic modifications must survive the ubiquitous ongoing activity present in neural circuits in vivo. In this paper, we investigate the time scales of memory maintenance in a calcium-based synaptic plasticity model that has been shown recently to be able to fit different experimental data-sets from hippocampal and neocortical preparations. We find that in the presence of background activity on the order of 1 Hz parameters that fit pyramidal layer 5 neocortical data lead to a very fast decay of synaptic efficacy, with time scales of minutes. We then identify two ways in which this memory time scale can be extended: (i) the extracellular calcium concentration in the experiments used to fit the model are larger than estimated concentrations in vivo. Lowering extracellular calcium concentration to in vivo levels leads to an increase in memory time scales of several orders of magnitude; (ii) adding a bistability mechanism so that each synapse has two stable states at sufficiently low background activity leads to a further boost in memory time scale, since memory decay is no longer described by an exponential decay from an initial state, but by an escape from a potential well. We argue that both features are expected to be present in synapses in vivo. These results are obtained first in a single synapse connecting two independent Poisson neurons, and then in simulations of a large network of excitatory and inhibitory integrate-and-fire neurons. Our results emphasise the need for studying plasticity at physiological extracellular calcium concentration, and highlight the role of synaptic bi- or multistability in the stability of learned synaptic structures.     

A dual receptor crosstalk model of G-protein-coupled signal transduction

Macrophage cells that are stimulated by two different ligands that bind to G-protein-coupled receptors (GPCRs) usually respond as if the stimulus effects are additive, but for a minority of ligand combinations the response is synergistic. The G-protein-coupled receptor system integrates signaling cues from the environment to actuate cell morphology, gene expression, ion homeostasis, and other physiological states. We analyze the effects of the two signaling molecules complement factors 5a (C5a) and uridine diphosphate (UDP) on the intracellular second messenger calcium to elucidate the principles that govern the processing of multiple signals by GPCRs. We have developed a formal hypothesis, in the form of a kinetic model, for the mechanism of action of this GPCR signal transduction system using data obtained from RAW264.7 macrophage cells. Bayesian statistical methods are employed to represent uncertainty in both data and model parameters and formally tie the model to experimental data. When the model is also used as a tool in the design of experiments, it predicts a synergistic region in the calcium peak height dose response that results when cells are simultaneously stimulated by C5a and UDP. An analysis of the model reveals a potential mechanism for crosstalk between the Galphai-coupled C5a receptor and the Galphaq-coupled UDP receptor signaling systems that results in synergistic calcium release.     

Time-dependent restoration of the trigger pool of calcium after termination of angiotensin II action in adrenal glomerulosa cells

In adrenal glomerulosa cells, angiotensin II causes an immediate release of calcium from an intracellular trigger pool (Kojima, I., Kojima, K., and Rasmussen, H. (1985) Am. J. Physiol. 247, E36-E43). The present study was conducted to determine how the trigger pool of calcium is restored after cessation of the agonist action. Upon termination of angiotensin II action, calcium influx rate decreased immediately while total cell calcium increased rapidly. The increase in total cell calcium is not affected by 1 microM nitrendipine, which blocks angiotensin II-stimulated calcium influx without inhibiting basal influx of calcium. In contrast, total cell calcium did not increase in medium containing 1 microM calcium, in which basal calcium influx is negligible. A rapid increase in total cell calcium after an addition of the antagonist was not accompanied by changes in cytoplasmic free calcium concentration. A second stimulation of cells with either angiotensin II or carbachol did not cause calcium release when the interval of two stimulations was shorter than 20 min. The longer the interval, the greater the magnitude of calcium release in response to the second stimulator. The maximum response was obtained when the interval was 40 min or more. When exogenous arachidonic acid, which mobilized calcium by acting directly on the inositol trisphosphate-sensitive pool, was employed as a second stimulator, the magnitude of the decrease in total cell calcium was also dependent on the interval. These results suggest that, upon termination of angiotensin II action, calcium is rapidly accumulated first in an intracellular pool which is insensitive to either inositol 1,4,5-trisphosphate or arachidonic acid and that the trigger pool is restored gradually thereafter.     

Mixtures of tight-binding enzyme inhibitors. Kinetic analysis by a recursive rate equation.

When two or more tight-binding inhibitors are present in an enzyme assay, the equation that relates the initial velocity v to the concentration of reactants cannot be written in an algebraically explicit form. Rather, for n inhibitors it is an implicit polynomial equation of degree n + 1 with respect to v. The complexity of the polynomial coefficients dramatically increases with each added inhibitor. Solving the transcendental rate equation by traditional methods of numerical mathematics has proven tedious because of the sensitivity of these methods to initial estimates and because of the existence of multiple roots. However, the equation can be rearranged into a convenient recursive form, one in which the velocity appears on both sides and the solution is found iteratively. The algebraic form of the recursive rate equation is remarkably simple and differs from the rate equation for classical rather than tight-binding inhibition only by an added term. The numerical stability and the speed of convergence were tested on the case of two competitive inhibitors. Initial estimates of velocity that spanned 12 orders of magnitude converged within five iterations. The velocities computed with the recursive method for a single tight-binding inhibitor were identical with the values predicted by the Morrison equation. The method is used to analyze experimental data for the inhibition of rat liver dihydrofolate reductase by mixtures of the anticancer drug methotrexate and its metabolic precursor form, methotrexate-alpha-aspartate (a prodrug).     

A polynomial based model for cell fate prediction in human diseases

Background

Cell fate regulation directly affects tissue homeostasis and human health. Research on cell fate decision sheds light on key regulators, facilitates understanding the mechanisms, and suggests novel strategies to treat human diseases that are related to abnormal cell development.

Results

In this study, we proposed a polynomial based model to predict cell fate. This model was derived from Taylor series. As a case study, gene expression data of pancreatic cells were adopted to test and verify the model. As numerous features (genes) are available, we employed two kinds of feature selection methods, i.e. correlation based and apoptosis pathway based. Then polynomials of different degrees were used to refine the cell fate prediction function. 10-fold cross-validation was carried out to evaluate the performance of our model. In addition, we analyzed the stability of the resultant cell fate prediction model by evaluating the ranges of the parameters, as well as assessing the variances of the predicted values at randomly selected points. Results show that, within both the two considered gene selection methods, the prediction accuracies of polynomials of different degrees show little differences. Interestingly, the linear polynomial (degree 1 polynomial) is more stable than others. When comparing the linear polynomials based on the two gene selection methods, it shows that although the accuracy of the linear polynomial that uses correlation analysis outcomes is a little higher (achieves 86.62%), the one within genes of the apoptosis pathway is much more stable.

Conclusions

Considering both the prediction accuracy and the stability of polynomial models of different degrees, the linear model is a preferred choice for cell fate prediction with gene expression data of pancreatic cells. The presented cell fate prediction model can be extended to other cells, which may be important for basic research as well as clinical study of cell development related diseases.

     

Transferable species distribution modelling: Comparative performance of Generalised Functional Response models

Predictive species distribution models (SDMs) are becoming increasingly important in ecology, in the light of rapid environmental change. However, the predictions of most current SDMs are specific to the habitat composition of the environments in which they were fitted. This may limit SDM predictive power because species may respond differently to a given habitat depending on the availability of all habitats in their environment, a phenomenon known as a functional response in resource selection. The Generalised Functional Response (GFR) framework captures this dependence by formulating the SDM coefficients as functions of habitat availability. The original GFR implementation used global polynomial functions of habitat availability to describe the functional responses. In this study, we develop several refinements of this approach and compare their predictive performance using two simulated and two real datasets. We first use local radial basis functions (RBF), a more flexible approach than global polynomials, to represent the habitat selection coefficients, and balance bias with precision via regularization to prevent overfitting. Second, we use the RBF-GFR and GFR models in combination with the classification and regression tree CART, which has more flexibility and better predictive powers for non-linear modelling. As further extensions, we use random forests (RFs) and extreme gradient boosting (XGBoost), ensemble approaches that consistently lead to variance reduction in generalization error. We find that the different methods are ranked consistently across the datasets for out-of-data prediction. The traditional stationary approach to SDMs and the GFR model consistently perform at the bottom of the ranking (simple SDMs underfit, and polynomial GFRs overfit the data). The best methods in our list provide non-negligible improvements in predictive performance, in some cases taking the out-of-sample R2 from 0.3 up to 0.7 across datasets. At times of rapid environmental change and spatial non-stationarity ignoring the effects of functional responses on SDMs, results in two different types of prediction bias (under-prediction or mis-positioning of distribution hotspots). However, not all functional response models perform equally well. The more volatile polynomial GFR models can generate biases through over-prediction. Our results indicate that there are consistently robust GFR approaches that achieve impressive gains in transferability across very different datasets.     

Timing of Ca2+ response in pancreatic beta-cells is related to mitochondrial mass

The timing and magnitude of calcium response are cell-specific in individual beta-cells. This may indicate that the cells have different roles in the intact islet. It is unknown what mechanisms determine these characteristics. We previously found that the mechanisms setting cell-specific response timing are disturbed in beta-cells from hyperglycemic mice and one of the causes is likely to be an altered mitochondrial metabolism. Mitochondria play a key role in the control of nutrient-induced insulin secretion. Here, we used confocal microscopy with the fluorescent probe MitoTracker Red CMXRos and Fluo-3 to study how the amount of active mitochondria is related to the lag-time and the magnitude of calcium response to 20mM glucose in isolated beta-cells and in cells within intact lean and ob/ob mouse islets. Results show that the mitochondrial mass is inversely correlated with the lag-times for calcium response both in lean and ob/ob mouse beta-cells (r=-0.73 and r=-0.43, respectively, P<0.05). Thus, the state of mitochondria may determine the timing of calcium response.     

Fura 2 analysis of cytosolic calcium regulation in elutriated rat gastric parietal cells

The calcium probe, Fura 2, is used to establish and partially characterize histamine-, carbachol-, and forskolin-induced calcium transients in enriched parietal cell populations prepared by centrifugal elutriation of dispersed rat fundic mucosa cell isolates. The magnitude of the maximal carbachol response, which is blocked by atropine but not cimetidine, is nearly five times that of histamine or forskolin. Time to peak responses for carbachol, forskolin, and histamine are approximately 7, 17, and 28 sec, respectively. Carbachol-, histamine-, and forskolin-induced increases in Fura 2 fluorescence appear dependent upon extracellular calcium, since these responses are attenuated in low calcium media and blocked by EGTA in low-calcium media or by lanthanum in high- or low-calcium medium. Trifluoperazine and fenoctimine, at concentrations that inhibit secretion, have no effect on either carbachol- or histamine-induced increases in cytosolic calcium. Seven major calcium/EGTA-sensitive phosphoproteins are identified by SDS-PAGE electrophoresis of ATP 32P-labeled cell sonicates. We conclude that cytosolic calcium in enriched rat gastric parietal cell populations is regulated by secretagogue receptor-controlled calcium channels. We postulate that these channels may be controlled by cyclic AMP-dependent phosphorylation, since neither changes in cyclic AMP nor calcium alone mediate the effects of secretagogues entirely, but the interplay between these two second-messenger systems potentiates the actions of these agents. The role of cytosolic calcium as a second messenger in secretagogue action appears similar to that of cyclic AMP in that a specific cellular concentration must be reached to initiate acid secretion.     

Osmometric analysis of quasi-indefinite mixed associations of the first kind

Mixed associations of the type A + B----AB, A + AB----A2B, ..., A + Ai-1 B----AiB, ... are readily analyzed by osmometric methods. The equilibrium molar concentration of A, mA, is obtained very simply from mA = meq-m0B; here meq = c/Meqn is the equilibrium molar concentration of all associating species and m0B denotes the stoichiometric or original molar concentration of B. The quantity mB can then be obtained from methods developed by Steiner. The value of the binding polynomial lambda is given by lambda = m0B/mB; lambda is a function of mA only. In principle, one can evaluate the equilibrium constants (kA,B,etc.) by fitting lambda to the appropriate polynomial in mA of degree n (n = 2, 3, ...). The binding polynomial lambda is analogous to polynomials encountered in the analysis of self-associations. By making some simple assumptions one can develop four analogs of two sequential, equal equilibrium constant (SEK) or two attenuated equilibrium constant (AK) models. With the aid of r (the number average degree of binding), g (the osmotic coefficient), lambda, as well as mA and mB, one can evaluate the equilibrium constant or constants. The methods developed here can be extended to the nonideal case.     

Thapsigargin potentiates histamine-stimulated HCl secretion in gastric parietal cells but does not mimic cholinergic responses.

The role of calcium in control of HCl secretion by the gastric parietal cell was examined using a recently available intracellular calcium-releasing agent, thapsigargin, which has been shown, in some cell types, to induce sustained elevation of intracellular calcium ([Ca2+]i), an action that appears to be independent of inositol lipid breakdown and protein kinase C activation and to be mediated, at least partially, by selective inhibition of endoplasmic reticulum Ca2(+)-ATPase. Using the calcium-sensitive fluorescent probe, fura-2, in combination with digitized video image analysis of single cells as well as standard fluorimetric techniques, we found that thapsigargin induced sustained elevation of [Ca2+]i in single parietal cells and in parietal cells populations. Chelation of medium calcium led to a transient rise and fall in [Ca2+]i, indicating that the sustained elevation in [Ca2+]i in response to thapsigargin was due to both intracellular calcium release and influx. Although thapsigargin appeared to affect the same calcium pool(s) regulated by the cholinergic agonist, carbachol, and the pattern of thapsigargin-induced increases in [Ca2+]i were similar to the plateau phase of the cholinergic response, thapsigargin did not induce acid secretory responses of the same magnitude as those initiated by carbachol (28 vs 600% of basal). The protein kinase C activator, 12-O-tetradecanoyl phorbol-13-acetate (TPA) potentiated the secretory response to thapsigargin but this combined response also did not attain the same magnitude as the maximal cholinergic response. In the presence but not the absence of medium calcium, thapsigargin potentiated acid secretory responses to histamine, which elevate both cyclic AMP (cAMP) and [Ca2+]i in parietal cells, as well as forskolin and cAMP analogues but had no effect on submaximal and an inhibitory effect on maximal cholinergic stimulation. Furthermore, thapsigargin did not fully mimic potentiating interactions between histamine and carbachol, either in magnitude or in the pattern of temporal response. Assuming that the action of thapsigargin is specific for intracellular calcium release mechanisms, these data suggest that 1) sustained influx of calcium is necessary but not sufficient for cholinergic activation of parietal cell HCl secretion and for potentiating interactions between cAMP-dependent agonists and carbachol; 2) mechanisms in addition to elevated [Ca2+]i and protein kinase C activation may be involved in cholinergic regulation; and 3) increases in [Ca2+]i in response to histamine are not directly involved in the mechanism of histamine-stimulated secretion.     

Changes in gravity rapidly alter the magnitude and direction of a cellular calcium current

In single-celled spores of the fern Ceratopteris richardii, gravity directs polarity of development and induces a directional, trans-cellular calcium (Ca²⁺) current. To clarify how gravity polarizes this electrophysiological process, we measured the kinetics of the cellular response to changes in the gravity vector, which we initially estimated using the self-referencing calcium microsensor. In order to generate more precise and detailed data, we developed a silicon microfabricated sensor array which facilitated a lab-on-a-chip approach to simultaneously measure calcium currents from multiple cells in real time. These experiments revealed that the direction of the gravity-dependent polar calcium current is reversed in less than 25 s when the cells are inverted, and that changes in the magnitude of the calcium current parallel rapidly changing g-forces during parabolic flight on the NASA C-9 aircraft. The data also revealed a hysteresis in the response of cells in the transition from 2g to micro-g in comparison to cells in the micro-g to 2-g transition, a result consistent with a role for mechanosensitive ion channels in the gravity response. The calcium current is suppressed by either nifedipine (calcium-channel blocker) or eosin yellow (plasma membrane calcium pump inhibitor). Nifedipine disrupts gravity-directed cell polarity, but not spore germination. These results indicate that gravity perception in single plant cells may be mediated by mechanosensitive calcium channels, an idea consistent with some previously proposed models of plant gravity perception.     

Effects of light intensity and restraint on dark-pulse-induced circadian phase shifting during subjective night in Syrian hamsters

Dark pulses presented on a background of constant light (LL) result in phase advances during midsubjective day and early subjective night, and phase delays during late subjective night, as shown in the dark-pulse phase response curve. In hamsters, the phase-shifting effects of dark pulses are thought to be mediated by increased activity, as previous studies have shown that restraining animals during dark pulses blocks the phase shifts observed in midsubjective day and late subjective night. This study focuses on dark-pulse-induced phase shifting during early subjective night, examining the influence of both LL intensity and restraint on the magnitude of these phase shifts. Syrian hamsters were maintained in LL of four different illumination levels (1, 10, 100, or 600 lux) and periodically presented with 6-h pulses (dark pulse alone, restraint alone, or dark pulse plus restraint) beginning at circadian time 11. Phase advances were observed in response to dark pulses alone, and the magnitude of these shifts was dependent on background illumination, with significantly larger advances seen under higher intensities. No relationship was found between the amount of activity displayed during dark pulses and phase shift magnitude. Six-hour periods of restraint resulted in phase delays, the magnitude of which was also dependent on background illumination. Restraining hamsters during dark pulses reduced the magnitude of phase advances, but the extent of this reduction could be predicted from the additive effects of the dark-pulse-alone and restraint-alone conditions. These results indicate that the phase-shifting effects of dark pulses during early subjective night are not mediated by behavioral activation and may instead reflect a mirror image of the phase-delaying effects of light pulses at this phase.     

A solution to the problems of cytolysis assays with additional applications to other immunological and biochemical assays

After cellular immunoassays are compared with classical bioassays, conventional methods and consequent problems of data analysis for cytolysis assays are reviewed and a new solution is proposed. This solution incorporates new methods, called dose-response surface assays and analysis (DRSA), which estimate cytolytic activity coefficients on a surface in a three-dimensional space with two dose variables (killers and targets) and one response variable (counts). These new methods based on dose-response surfaces are demonstrated to be more informative and reliable than classical methods based on dose-response curves. In a test of the methods' robustness (sensitivity of parameter estimates to changes in the dose levels of the assay design), cytolytic activity coefficients estimated by DRSA varied by less than or equal to 30% over a reduction of three to four orders of magnitude in the dose levels. This remarkable robustness should be compared with the corresponding figures of as much as 500% over less than 1 order of magnitude for previously published results of coefficients estimated by conventional methods. DRSA is distinguished from replot-of-plots methods such as those used for enzyme inhibition assays in biochemistry, and is recommended as a more efficient method that should replace replot-of-plot methods now antiquated by the advent of microcomputers. DRSA can be applied to any experimental system that requires an activity coefficient to be estimated on a dose-response surface in a space of greater than or equal to 3 dimensions (greater than or equal to 2 dose variables and one response variable), regardless of the mathematical model and statistical estimators used to analyze the dose-response interaction. Finally, DRSA is compared with the methods known as response surface methodology (RSM), and is described as a new class of methods to be added to those that constitute RSM.     

Corticosterone Enhances Kainic Acid-Induced Calcium Elevation in Cultured Hippocampal Neurons

Abstract: Corticosterone, a steroid secreted during stress, increases hippocampal neuronal vulnerability to excitotoxins, hypoxia-ischemia, and antimetabolites. Energy supplementation and N -methyl-d-aspartate receptor antagonists prevent this corticosterone-enhanced neurotoxicity. Because neuronal calcium regulation is energy dependent and a large calcium influx accompanies N -methyl-d-aspartate receptor activation, we investigated whether corticosterone exacerbates the elevation of hippocampal neuronal calcium induced by the glutamatergic excitotoxin kainic acid. Corticosterone caused a 23-fold increase in the magnitude of the calcium response to kainic acid, a sevenfold increase in the peak magnitude of the calcium response, and a twofold increase in calcium recovery time. This corticosterone effect may be energetic in nature as corticosterone decreases hippocampal neuronal glucose transport. Glucose supplementation reduced the corticosterone effect on the magnitude and peak magnitude of the calcium response to kainic acid. Glucose reduction, by the approximate magnitude by which corticosterone inhibits glucose transport, mimicked the corticosterone effect on the peak magnitude of the calcium response to kainic acid. Thus, corticosterone increases calcium after kainic acid exposure in hippocampal neurons in an energy-dependent manner. Elevated calcium is strongly implicated in stimulating neurotoxic cascades during other energetic insults and may be the mechanism for the corticosterone-induced hippocampal neuronal vulnerability and toxicity.     

Agents that raise cAMP in human T lymphocytes release an intracellular pool of calcium in the absence of inositol phosphate production

Prostaglandins (PGs) of the E series are recognized by specific receptors on T lymphocytes which lead to an increase in cAMP. The role of cAMP in modulation of T lymphocyte function is unknown. Here, we demonstrate that agents which increase cAMP in human T cells raise the intracellular free calcium concentration ([Ca2+]i). This increase in [Ca2+]i occurred following receptor stimulation with PGEs or by bypassing the receptor with the cell-permeant analog 8-(4-chlorophenylthio)-cAMP or forskolin, a direct activator of adenylyl cyclase. The calcium response to a submaximally stimulatory concentration of PGE2 was potentiated by the cAMP phosphodiesterase inhibitor isobutylmethylxanthine. A time course of cAMP production in response to PGE2 stimulation closely resembled the calcium response and suggested that the two events were coincident. The PGE2 concentrations required to achieve 50% maximum effect of cAMP production and increases in [Ca2+]i were similar, 0.07 and 0.15 microM respectively. Chelation of extracellular Ca2+ did not abolish the PGE2-stimulated Ca2+ response, suggesting that an intracellular source of calcium was sensitive to cAMP. Significant inositol phosphate production was not detected in response to PGE2 over a wide concentration range. The PGE2-induced calcium response curves were of lesser magnitude with shorter times to peak than those of a known inositol 1,4,5 trisphosphate-producing agonist, anti-CD3, suggesting distinct Ca2+ release mechanisms. However, the cAMP-releasable store appeared to be contained within the inositol trisphosphate-releasable store since no response could be seen with cAMP-elevating agents following emptying of the inositol trisphosphate-sensitive pool of Ca2+.     

Integrated mismatch negativity (MMNi): a noise-free representation of evoked responses allowing single-point distribution-free statistical tests

If the repeated presentation of a single (standard) auditory stimulus is randomly interspersed with a second acoustically different (deviant) stimulus, the cortical activity evoked by the deviant stimulus can contain a negative component known as the mismatch negativity (MMN). The MMN is derived by subtracting the averaged response evoked by the standard stimulus from that evoked by the deviant stimulus. When the magnitude of the response is small or the signal-to-noise ratio is poor, it is difficult to judge the presence or absence of the MMN simply by visual inspection, and statistical detection techniques become necessary. A method of analysis is proposed to quantify the magnitude and statistically evaluate the presence of the MMN based on time-integrated evoked responses. This paper demonstrates the use of this integrated mismatch negativity (MMN_i) analysis to detect the MMN evoked by stimulus contrasts near the perceptual threshold of two subjects. The MMN_i, by virtue of being equivalent to a low-pass filtered response, presents an almost noise-free estimate of MMN magnitude. A single measure of the integrated evoked response at a fixed time point is used in a distribution-free statistic that compares the magnitude of the averaged response evoked by the deviant stimulus with a magnitude distribution derived from 200 subaveraged responses to the standard stimulus (with the number of sweeps per average equal to that of the deviant stimulus). This allows a calculation of the exact probability for the null hypothesis that the negative magnitude of the response evoked by the deviant stimulus is drawn from the magnitude distribution of responses evoked by the standard stimulus. Rejection of this hypothesis provides objective evidence of the presence of the MMN.