Soil fertility determines whether ectomycorrhizal fungi accelerate or decelerate decomposition in a temperate forest

Ectomycorrhizal (ECM) fungi can both accelerate and decelerate decomposition of organic matter in forest soils, but a mechanistic understanding of this differential influence is limited. Here, we tested how ECM fungi affect decomposition along a natural fertility gradient in a temperate forest of European beech. Trees were girdled to reduce belowground carbon supply to the soil. Girdling shifted soil fungal community composition and decreased hyphal biomass production and soil CO₂ efflux, indicating a reduced ECM fungal activity. Girdling also affected decomposition processes, but the effects depended on fertility. Our results indicate that ECM fungi decelerate decomposition under conditions of low fertility while under conditions of high fertility ECM fungi and their host roots have an accelerating effect. We conclude that both acceleration and deceleration of decomposition of organic matter by ECM fungi can occur within a forest, with soil fertility determining the direction and magnitude of these effects. We suggest a positive feedback between fertility, stand productivity and soil carbon and nitrogen dynamics that is mediated to a large extent by ECM fungi.

Soil fertility relates to fungal-mediated decomposition and organic matter turnover in a temperate mountain forest

Fungi are known to exert a significant influence over soil organic matter (SOM) turnover, however understanding of the effects of fungal community structure on SOM dynamics and its consequences for ecosystem fertility is fragmentary. Here we studied soil fungal guilds and SOM decomposition processes along a fertility gradient in a temperate mountain beech forest. High-throughput sequencing was used to investigate fungal communities. Carbon and nitrogen stocks, enzymatic activity and microbial respiration were measured. While ectomycorrhizal fungal abundance was not related to fertility, saprotrophic ascomycetes showed higher relative abundances under more fertile conditions. The activity of oxidising enzymes and respiration rates in mineral soil were related positively to fertility and saprotrophic fungi. In addition, organic layer carbon and nitrogen stocks were lower on the more fertile plots, although tree biomass and litter input were higher. Together, the results indicated a faster SOM turnover at the fertile end of the gradient. We suggest that there is a positive feedback mechanism between SOM turnover and fertility that is mediated by soil fungi to a significant extent. By underlining the importance of fungi for soil fertility and plant growth, these findings furthermore emphasise the dependency of carbon cycling on fungal communities below ground.

Modeling the biogeochemistry of sulfur in beech (Fagus sylvatica L.) stands of the Vienna Woods

In this article, we describe the setup and the application of a novel Sulfur Dynamics Model (SDM), with the aim to identify and quantify processes, which are crucial for the understanding of the biogeochemical sulfur (S) cycle of forest ecosystems. The simulator takes into account abiotic processes as well as biotic interactions between atmosphere, plant and soil. We applied the model to two Austrian beech stands where deposition of S and soil solution chemistry were monitored intensively over a two-year period. Under consideration of high historic loads and the more recent recession of atmospheric S deposition, we found a suitable model configuration where it was possible to assign both intra-annual fluctuations of the SO4-S in soil solution and long-time trends in the stream discharge to specific S transformation processes. We identified the interplay of microbial immobilization (the microbial conversion of solute sulfate to organic soil S) and mineralization as key driver of short-term fluctuations in the soil solution. In the long term, the delayed release of historically accumulated S is driven mainly by the slow mineralization of S rich plant biomass, recalcitrant to decomposition. Adsorption and desorption processes seem to play only a negligible role on our investigated stands. We conclude that our proposed model which is based on the current understanding of S biogeochemistry is sufficient to describe S dynamics on the investigated forest stands. The code file (SAS) for all model functions will be provided by the authors after request.

The impact of rising temperatures on water balance and phenology of European beech (Fagus sylvatica L.) stands

In this article, we outline the set-up and the application of an eco-hydrological box model, with the aim to describe the water balance of deciduous (Fagus Sylvatica L.) forest stands. The water balance model (WBM) uses standard meteorological parameters as input variables and runs on a daily time step. It consists of two modules. The aboveground module (1) comprises routines for fog precipitation generation, precipitation interception and snowfall/snowmelt dynamics. Covered belowground processes (2) are bypass flow, percolation, soil evaporation and transpiration, where the latter two processes are considered separately. Preceding to the WBM, a routine is introduced, specifying the intra-annual foliage dynamics of beech. Emphasis is also laid on the inter-annual variation of beech phenology. Leaf sprouting and leaf senescence are calculated as functions of day-length and air temperature. The WBM was applied to four European beech dominated forest stands in the northeastern part of Austria. They are located on a gradient of declining annual precipitation (from west to east). The two easterly sites are located close to the (dry) limit of the natural distribution of beech. Records of soil moisture were used for the adjustment of 26 parameters. On all sites the calibration process (simulated annealing) delivered good predictions of soil moisture (Nash–Sutcliffe efficiency≥ 0.925). Then, the obtained parameterization was used to apply different scenarios of global warming. The temperature was increased step-wisely up to 4 °C. All scenarios were run (1) with present phenological conditions and (2) with phenology responding to higher temperatures. This way, we wanted to assign the effect of higher temperatures and longer growing seasons on the water dynamics of the forest stands. A warming of 1 °C corresponded roughly to an elongation of the growing season of 4.5 days, where the start of the growing season was affected more strongly than the end. Apparently, higher temperatures led to drier soils. The strongest change was observed in early summer, also amplified by an earlier start of the growing season. Rising temperatures led to lower export fluxes of liquid water, simultaneously increasing evapotranspiration (ET). The gain in ET was almost entirely assignable to increased soil evaporation. Drier soils led to a sharp depression of transpiration during summer months. This decline was compensated by the effect of elongated growing seasons. The risk of severe drought was increased by higher temperatures, but here the contribution of growing season length was negligible. Drier soils seem to hamper the stands’ productivity. For all warming scenarios, the estimated increase of the gross primary production, caused by longer periods of assimilation, is nullified by the effect of soil water deficit in mid-summer.

Reconstructing Soil Recovery from Acid Rain in Beech (Fagus sylvatica) Stands of the Vienna Woods as Indicated by Removal of Stemflow and Dendrochemistry

Our goal was to reconstruct soil recovery from Acid Rain based upon removal of stemflow at beech (Fagus sylvatica) stands of known historic and recent soil status. Fourteen beech stands in the Vienna Woods were selected in 1984 and again in 2012 to study changes in soil and foliar chemistry over time. A part of those stands had been strip cut, and to assess reversibility of soil acidification, we analyzed soils around beech stumps from different years of felling, representing the years when acidic stemflow ceased to affect the soil. Furthermore, it was hypothesized that changes of soil chemistry are reflected in the stemwood of beech. Half-decadal samples of tree cores were analyzed for Ca, Mg, K, Mn, Fe, and Al. Soil analyses indicated recovery in the top soil of the stemflow area but recovery was delayed in the between trees areas and deeper soil horizons. Differences in soil pH between proximal and distal area from beech stumps were still detectable after 30 years indicating that soils may not recover fully from acidification or do so at a rather slow rate. Stemwood contents indicated mobilization of base cations during the early 80s followed by a steady decrease thereafter. Backward reconstructions of soil pH and soil nutrients, building on regressions between recent stemwood and soil chemistry, could not be verified by measured soil data in 1984, but matched with declining cation foliar contents from 1984 to 2012. Dendrochemical reconstructions showed highest values in the 1980s, but measured soil exchangeable cation contents were clearly lower in 1984. Hence, we conclude that our reconstructions mimicked soil solution rather than soil exchanger chemistry.

Declining atmospheric deposition of heavy metals over the last three decades is reflected in soil and foliage of 97 beech (Fagus sylvatica) stands in the Vienna Woods

Rigorous studies on long-term changes of heavy metal distribution in forest soils since the implementation of emission controls are rare. Hence, we resampled 97 old-growth beech stands in the Vienna Woods. This study exploits an extensive data set of soil (infiltration zone of stemflow and between trees area) and foliar chemistry from three decades ago. It was hypothesized that declining deposition of heavy metals is reflected in soil and foliar total contents of Pb, Cu, Zn, Ni, Mn and Fe. Mean soil contents of Pb in the stemflow area declined at the highest rate from 223 to 50 mg kg^-1 within the last three decades. Soil contents of Pb and Ni decreased significantly both in the stemflow area and the between trees area down to 80-90 cm soil depth from 1984 to 2012. Top soil (0-5 cm) accumulation and simultaneous loss in the lower soil over time for the plant micro nutrients Cu and Zn are suggested to be caused by plant uptake from deep horizons. Reduced soil leaching, due to a mean soil pH (H₂O) increase from 4.3 to 4.9, and increased plant cycling are put forward to explain the significant increase of total Mn contents in the infiltration zone of beech stemflow. Top soil Pb contents in the stemflow area presently exceed the critical value at which toxicity symptoms may occur at numerous sites. Mean foliar contents of all six studied heavy metals decreased within the last three decades, but plant supply with the micro nutrients Cu, Zn, Mn and Fe is still in the optimum range for beech trees. It is concluded that heavy metal pollution is not critical for the studied beech stands any longer.

Fractionation of sulfur (S) in beech (Fagus sylvatica) forest soils in relation to distance from the stem base as useful tool for modeling S biogeochemistry

The investigation of the fractionation of S compounds in forest soils is a powerful tool for interpreting S dynamics and S biogeochemistry in forest ecosystems. Beech stands on high pH (nutrient-rich) sites on Flysch and on low pH (nutrient-poor) sites on Molasse were selected for testing the influence of stemflow, which represents a high input of water and dissolved elements to the soil, on spatial patterns of sulfur (S) fractions. Soil cores were taken at six distances from a beech stem per site at 55 cm uphill and at 27, 55, 100, 150 and 300 cm downhill from the stem. The cores were divided into the mineral soil horizons 0-3, 3-10, 10-20, 20-30 and 30-50 cm. Soil samples were characterized for pH, C_org, pedogenic Al and Fe oxides and S fractions. Sequential extraction by NH₄Cl, NH₄H₂PO₄ and HCl yielded readily available sulfate-S (RAS), adsorbed sulfate-S (AS) and HCl-soluble sulfate-S (HCS). Organic sulfur (OS) was estimated as the difference between total sulfur (ToS) and inorganic sulfur (RAS + AS + HCS). Organic sulfur was further divided into ester sulfate-S (ES, HI-reduction) and carbon bonded sulfur (CS). On Flysch, RAS represented 3-6%, AS 2-12%, HCS 0-8% and OS 81-95% of ToS. On Molasse, RAS amounted 1-6%, AS 1-60%, HCS 0-8% and OS 37-95% of ToS. Spatial S distribution patterns with respect to the distance from the tree stem base could be clearly observed at all investigated sites. The presented data is a contribution to current reports on negative input-output S budgets of forest watersheds, suggesting that mineralization of OS on nutrient rich soils and desorption of historic AS on nutrient-poor soils are the dominant S sources, which have to be considered in future modeling of sulfur.

Sources of errors and uncertainties in the assessment of forest soil carbon stocks at different scales-review and recommendations

Spatially explicit knowledge of recent and past soil organic carbon (SOC) stocks in forests will improve our understanding of the effect of human- and non-human-induced changes on forest C fluxes. For SOC accounting, a minimum detectable difference must be defined in order to adequately determine temporal changes and spatial differences in SOC. This requires sufficiently detailed data to predict SOC stocks at appropriate scales within the required accuracy so that only significant changes are accounted for. When designing sampling campaigns, taking into account factors influencing SOC spatial and temporal distribution (such as soil type, topography, climate and vegetation) are needed to optimise sampling depths and numbers of samples, thereby ensuring that samples accurately reflect the distribution of SOC at a site. Furthermore, the appropriate scales related to the research question need to be defined: profile, plot, forests, catchment, national or wider. Scaling up SOC stocks from point sample to landscape unit is challenging, and thus requires reliable baseline data. Knowledge of the associated uncertainties related to SOC measures at each particular scale and how to reduce them is crucial for assessing SOC stocks with the highest possible accuracy at each scale. This review identifies where potential sources of errors and uncertainties related to forest SOC stock estimation occur at five different scales-sample, profile, plot, landscape/regional and European. Recommendations are also provided on how to reduce forest SOC uncertainties and increase efficiency of SOC assessment at each scale.

A slight recovery of soils from Acid Rain over the last three decades is not reflected in the macro nutrition of beech (Fagus sylvatica) at 97 forest stands of the Vienna Woods

Rigorous studies of recovery from soil acidification are rare. Hence, we resampled 97 old-growth beech stands in the Vienna Woods. This study exploits an extensive data set of soil (infiltration zone of stemflow and between trees area at different soil depths) and foliar chemistry from three decades ago. It was hypothesized that declining acidic deposition is reflected in soil and foliar chemistry. Top soil pH within the stemflow area increased significantly by 0.6 units in both H2O and KCl extracts from 1984 to 2012. Exchangeable Ca and Mg increased markedly in the stemflow area and to a lower extent in the top soil of the between trees area. Trends of declining base cations in the lower top soil were probably caused by mobilization of organic S and associated leaching with high amounts of sulfate. Contents of C, N and S decreased markedly in the stemflow area from 1984 to 2012, suggesting that mineralization rates of organic matter increased due to more favorable soil conditions. It is concluded that the top soil will continue to recover from acidic deposition. However, in the between trees areas and especially in deeper soil horizons recovery may be highly delayed. The beech trees of the Vienna Woods showed no sign of recovery from acidification although S deposition levels decreased. Release of historic S even increased foliar S contents. Base cation levels in the foliage declined but are still adequate for beech trees. Increasing N/nutrient ratios over time were considered not the result of marginally higher N foliar contents in 2012 but of diminishing nutrient uptake due to the decrease in ion concentration in soil solution. The mean foliar N/P ratio already increased to the alarming value of 31. Further nutritional imbalances will predispose trees to vitality loss.

Novel diffusive gradients in thin films technique to assess labile sulfate in soil

A novel diffusive gradients in thin films (DGT) technique for sampling labile soil sulfate was developed, based on a strong basic anion exchange resin (Amberlite IRA-400) for sulfate immobilization on the binding gel. For reducing the sulfate background on the resin gels, photopolymerization was applied instead of ammonium persulfate-induced polymerization. Agarose cross-linked polyacrylamide (APA) hydrogels were used as diffusive layer. The sulfate diffusion coefficient in APA gel was determined as 9.83 × 10⁻⁶ ± 0.35 × 10^-6 cm² s⁻¹ at 25 °C. The accumulated sulfate was eluted in 1 mol L⁻¹ HNO₃ with a recovery of 90.9 ± 1.6 %. The developed method was tested against two standard extraction methods for soil sulfate measurement. The obtained low correlation coefficients indicate that DGT and conventional soil test methods assess differential soil sulfate pools, rendering DGT a potentially important tool for measuring labile soil sulfate.

A new approach to predict soil temperature under vegetated surfaces

In this article, the setup and the application of an empirical model, based on Newton’s law of cooling, capable to predict daily mean soil temperature (T_soil) under vegetated surfaces, is described. The only input variable, necessary to run the model, is a time series of daily mean air temperature. The simulator employs 9 empirical parameters, which were estimated by inverse modeling. The model, which primarily addresses forested sites, incorporates the effect of snow cover and soil freezing on soil temperature. The model was applied to several temperate forest sites, managing the split between Central Europe (Austria) and the United States (Harvard Forest, Massachusetts; Hubbard Brook, New Hampshire), aiming to cover a broad range of site characteristics. Investigated stands differ fundamentally in stand composition, elevation, exposition, annual mean temperature, precipitation regime, as well as in the duration of winter snow cover. At last, to explore the limits of the formulation, the simulator was applied to non-forest sites (Illinois), where soil temperature was recorded under short cut grass. The model was parameterized, specifically to site and measurement depth. After calibration of the model, an evaluation was performed, using ~50 % of the available data. In each case, the simulator was capable to deliver a feasible prediction of soil temperature in the validation time interval. To evaluate the practical suitability of the simulator, the minimum amount of soil temperature point measurements, necessary to yield expedient model performance was determined. In the investigated case 13–20 point observations, uniformly distributed within an 11-year timeframe, have been proven sufficient to yield sound model performance (root mean square error <0.9 °C, Nash–Sutcliffe efficiency >0.97). This makes the model suitable for the application on sites, where the information on soil temperature is discontinuous or scarce.

MC ICP-MS δ(34)S(VCDT) measurement of dissolved sulfate in environmental aqueous samples after matrix separation by means of an anion exchange membrane

Analysis of ³⁴S/³²S of sulfate in rainwater and soil solutions can be seen as a powerful tool for the study of the sulfur cycle. Therefore, it is considered as a useful means, e.g., for amelioration and calibration of ecological or biogeochemical models. Due to several analytical limitations, mainly caused by low sulfate concentration in rainwater, complex matrix of soil solutions, limited sample volume, and high number of samples in ecosystem studies, a straightforward analytical protocol is required to provide accurate S isotopic data on a large set of diverse samples. Therefore, sulfate separation by anion exchange membrane was combined with precise isotopic measurement by multicollector inductively coupled plasma mass spectrometry (MC ICP-MS). The separation method proved to be able to remove quantitatively sulfate from matrix cations (Ca, K, Na, or Li) which is a precondition in order to avoid a matrix-induced analytical bias in the mass spectrometer. Moreover, sulfate exchange on the resin is capable of preconcentrating sulfate from low concentrated solutions (to factor 3 in our protocol). No significant sulfur isotope fractionation was observed during separation and preconcentration. MC ICP-MS operated at edge mass resolution has enabled the direct ³⁴S/³²S analysis of sulfate eluted from the membrane, with an expanded uncertainty U (k = 2) down to 0.3 ‰ (a single measurement). The protocol was optimized and validated using different sulfate solutions and different matrix compositions. The optimized method was applied in a study on solute samples retrieved in a beech (Fagus sylvatica) forest in the Vienna Woods. Both rainwater (precipitation and tree throughfall) and soil solution δ³⁴S_VCDT ranged between 4 and 6 ‰, the ratio in soil solution being slightly lower. The lower ratio indicates that a considerable portion of the atmospherically deposited sulfate is cycled through the organic S pool before being released to the soil solution. Nearly the same trends and variations were observed in soil solution and rainwater δ³⁴S_VCDT values showing that sulfate adsorption/desorption are not important processes in the studied soil.

Graphical abstract

Decomposition of beech (Fagus sylvatica) and pine (Pinus nigra) litter along an Alpine elevation gradient: Decay and nutrient release

Litter decomposition is an important process for cycling of nutrients in terrestrial ecosystems. The objective of this study was to evaluate direct and indirect effects of climate on litter decomposition along an altitudinal gradient in a temperate Alpine region. Foliar litter of European beech (Fagus sylvatica) and Black pine (Pinus nigra) was incubated in litterbags during two years in the Hochschwab massif of the Northern Limestone Alps of Austria. Eight incubation sites were selected following an altitudinal/climatic transect from 1900 to 900 m asl. The average remaining mass after two years of decomposition amounted to 54% (beech) and 50% (pine). Net release of N, P, Na, Al, Fe and Mn was higher in pine than in beech litter due to high immobilization (retention) rates of beech litter. However, pine litter retained more Ca than beech litter. Altitude retarded decay (mass loss and associated C release) in beech litter during the first year only but had a longer lasting effect on decaying pine litter. Altitude comprises a suite of highly auto-correlated characteristics (climate, vegetation, litter, soil chemistry, soil microbiology, snow cover) that influence litter decomposition. Hence, decay and nutrient release of incubated litter is difficult to predict by altitude, except during the early stage of decomposition, which seemed to be controlled by climate. Reciprocal litter transplant along the elevation gradient yielded even relatively higher decay of pine litter on beech forest sites after a two-year adaptation period of the microbial community.

Does mixing of beech (Fagus sylvatica) and spruce (Picea abies) litter hasten decomposition?

BACKGROUND AND AIMS

It is of practical relevance to know how much beech must be admixed to pure spruce stands in order to increase litter decomposition and associated nutrient cycling, since the formation of thick organic layers is commonly ascribed to the recalcitrance of spruce needles. We addressed the impact of tree species mixture within forest stands and within litter on mass loss and nutritional release from litter.

METHODS

Litter decomposition was measured in three adjacent stands of pure spruce (Picea abies), mixed beech-spruce and pure beech (Fagus sylvatica) on a nutrient-rich site and a nutrient-poor site over a 2-year period using litterbags which were filled with five different mixtures of beech and spruce litter.

RESULTS

Mass loss of beech litter was not higher than mass loss of spruce litter. Decay was primarily affected by tree species composition of the incubation stand and was faster in (mixed) beech forests stands than in spruce forests, while the influence of litter species and their mixtures on decay rates was small. Net transfers of nutrients between the two litter species (direct effects) in the mixed bags were minimal, since initial beech and spruce litter did not have different litter quality. However, in a few cases indirect effects (e.g., changing decomposer abundance and activity) caused non-additive patterns for the totals within the mixed bags, hastening decomposition within the first year.

CONCLUSIONS

Greater accumulation of litter in spruce compared to beech stands is not a consequence of the inherent recalcitrance of needles. Adverse environmental conditions in spruce stands retard decomposition. Indirect effects on decomposition caused by stand mixture are not mimicked by litter mixtures within mesh bags.

On parsing the neural code in the prefrontal cortex of primates using principal dynamic modes

Nonlinear modeling of multi-input multi-output (MIMO) neuronal systems using Principal Dynamic Modes (PDMs) provides a novel method for analyzing the functional connectivity between neuronal groups. This paper presents the PDM-based modeling methodology and initial results from actual multi-unit recordings in the prefrontal cortex of non-human primates. We used the PDMs to analyze the dynamic transformations of spike train activity from Layer 2 (input) to Layer 5 (output) of the prefrontal cortex in primates performing a Delayed-Match-to-Sample task. The PDM-based models reduce the complexity of representing large-scale neural MIMO systems that involve large numbers of neurons, and also offer the prospect of improved biological/physiological interpretation of the obtained models. PDM analysis of neuronal connectivity in this system revealed "input-output channels of communication" corresponding to specific bands of neural rhythms that quantify the relative importance of these frequency-specific PDMs across a variety of different tasks. We found that behavioral performance during the Delayed-Match-to-Sample task (correct vs. incorrect outcome) was associated with differential activation of frequency-specific PDMs in the prefrontal cortex.

Design of optimal stimulation patterns for neuronal ensembles based on Volterra-type hierarchical modeling

This paper presents a general methodology for the optimal design of stimulation patterns applied to neuronal ensembles in order to elicit a desired effect. The methodology follows a variant of the hierarchical Volterra modeling approach that utilizes input-output data to construct predictive models that describe the effects of interactions among multiple input events in an ascending order of interaction complexity. The illustrative example presented in this paper concerns the multi-unit activity of CA1 neurons in the hippocampus of a rodent performing a learned delayed-nonmatch-to-sample (DNMS) task. The multi-unit activity of the hippocampal CA1 neurons is recorded via chronically implanted multi-electrode arrays during this task. The obtained model quantifies the likelihood of having correct performance of the specific task for a given multi-unit (spatiotemporal) activity pattern of a CA1 neuronal ensemble during the 'sample presentation' phase of the DNMS task. The model can be used to determine computationally (off-line) the 'optimal' multi-unit stimulation pattern that maximizes the likelihood of inducing the correct performance of the DNMS task. Our working hypothesis is that application of this optimal stimulation pattern will enhance performance of the DNMS task due to enhancement of memory formation and storage during the 'sample presentation' phase of the task.

Modeling neuron-glia interactions: from parametric model to neuromorphic hardware

Recent experimental evidence suggests that glial cells are more than just supporting cells to neurons - they play an active role in signal transmission in the brain. We herein propose to investigate the importance of these mechanisms and model neuron-glia interactions at synapses using three approaches: A parametric model that takes into account the underlying mechanisms of the physiological system, a non-parametric model that extracts its input-output properties, and an ultra-low power, fast processing, neuromorphic hardware model. We use the EONS (Elementary Objects of the Nervous System) platform, a highly elaborate synaptic modeling platform to investigate the influence of astrocytic glutamate transporters on postsynaptic responses in the detailed micro-environment of a tri-partite synapse. The simulation results obtained using EONS are then used to build a non-parametric model that captures the essential features of glutamate dynamics. The structure of the non-parametric model we use is specifically designed for efficient hardware implementation using ultra-low power subthreshold CMOS building blocks. The utilization of the approach described allows us to build large-scale models of neuron/glial interaction and consequently provide useful insights on glial modulation during normal and pathological neural function.

Restorative encoding memory integrative neural device: "REMIND"

Construction and application of a neural prosthesis device that enhances existing and replaces lost memory capacity in humans is the focus of research described here in rodents. A unique approach for the analysis and application of neural population firing has been developed to decipher the pattern in which information is successfully encoded by the hippocampus where mnemonic accuracy is critical. A nonlinear dynamic multi-input multi-output (MIMO) model is utilized to extract memory relevant firing patterns in CA3 and CA1 and to predict online what the consequences of the encoded firing patterns reflect for subsequent information retrieval for successful performance of delayed-nonmatch-to-sample (DNMS) memory task in rodents. The MIMO model has been tested successfully in a number of different contexts, each of which produced improved performance by a) utilizing online predicted codes to regulate task difficulty, b) employing electrical stimulation of CA1 output areas in the same pattern as successful cell firing, c) employing electrical stimulation to recover cell firing compromised by pharmacological agents and d) transferring and improving performance in naïve animals using the same stimulation patterns that are effective in fully trained animals. The results in rodents formed the basis for extension of the MIMO model to nonhuman primates in the same type of memory task that is now being tested in the last step prior to its application in humans.

Greater accumulation of litter in spruce (Picea abies) compared to beech (Fagus sylvatica) stands is not a consequence of the inherent recalcitrance of needles

BACKGROUND AND AIMS

Replacement of beech by spruce is associated with changes in soil acidity, soil structure and humus form, which are commonly ascribed to the recalcitrance of spruce needles. It is of practical relevance to know how much beech must be admixed to pure spruce stands in order to increase litter decomposition and associated nutrient cycling. We addressed the impact of tree species mixture within forest stands and within litter on mass loss and nutritional release from litter.

METHODS

Litter decomposition was measured in three adjacent stands of pure spruce (Picea abies), mixed beech-spruce and pure beech (Fagus sylvatica) on three nutrient-rich sites and three nutrient-poor sites over a three-year period using the litterbag method (single species and mixed species bags).

RESULTS

Mass loss of beech litter was not higher than mass loss of spruce litter. Mass loss and nutrient release were not affected by litter mixing. Litter decay indicated non-additive patterns, since similar remaining masses under pure beech (47%) and mixed beech-spruce (48%) were significantly lower than under pure spruce stands (67%). Release of the main components of the organic substance (C_org, N_tot, P, S, lignin) and associated K were related to mass loss, while release of other nutrients was not related to mass loss.

CONCLUSIONS

In contradiction to the widely held assumption of slow decomposition of spruce needles, we conclude that accumulation of litter in spruce stands is not caused by recalcitrance of spruce needles to decay; rather adverse environmental conditions in spruce stands retard decomposition. Mixed beech-spruce stands appear to be as effective as pure beech stands in counteracting these adverse conditions.

Dynamic nonlinear modeling of interactions between neuronal ensembles using principal dynamic modes

We present a novel methodology for modeling the interactions between neuronal ensembles that utilizes the concept of Principal Dynamic Modes (PDM) and their associated nonlinear functions (ANF). This new approach seeks to reduce the complexity of the multi-input/multi-output (MIMO) model of the interactions between neuronal ensembles--an issue of critical practical importance in scaling up the MIMO models to incorporate hundreds (or even thousands) of input-output neurons. Global PDMs were extracted from the data using estimated first-order and second-order kernels and singular value decomposition (SVD). These global PDMs represent an efficient "coordinate system" for the representation of the MIMO model. The ANFs of the PDMs are estimated from the histograms of the combinations of PDM output values that lead to output spikes. For initial testing and validation of this approach, we applied it to a set of data collected at the pre-frontal cortex of a non-human primate during a behavioral task (Delayed Match-to-Sample). Recorded spike trains from Layer-2 neurons were viewed as the "inputs" and from Layer-5 neurons as the outputs. Model prediction performance was evaluated by means of computed Receiver Operating Characteristic (ROC) curves. The results indicate that this methodology may greatly reduce the complexity of the MIMO model without significant degradation of performance.

Modeling of nonlinear nonstationary dynamic systems with a novel class of artificial neural networks

This paper introduces a novel neural-network architecture that can be used to model time-varying Volterra systems from input-output data. The Volterra systems constitute a very broad class of stable nonlinear dynamic systems that can be extended to cover nonstationary (time-varying) cases. This novel architecture is composed of parallel subnets of three-layer perceptrons with polynomial activation functions, with the output of each subnet modulated by an appropriate time function that gives the summative output its time-varying characteristics. The paper shows the equivalence between this network architecture and the class of time-varying Volterra systems, and demonstrates the range of applicability of this approach with computer-simulated examples and real data. Although certain types of nonstationarities may not be amenable to this approach, it is hoped that this methodology will provide the practical tools for modeling some broad classes of nonlinear, nonstationary systems from input-output data, thus advancing the state of the art in a problem area that is widely viewed as a daunting challenge.

Improving bioassay sensitivity for neurotoxins detection using volterra based third order nonlinear analysis

Based on a novel analytical method for analyzing short-term plasticity (STP) of the CA1 hippocampal region in vitro, a screening tool for the detection and classification of unknown chemical compounds affecting the nervous system was recently introduced [1], [2]. The recorded signal consisted of evoked population spike in response to Poisson distributed random train impulse stimuli. The developed analytical approach used the first order Volterra kernel and the Laguerre coefficients of the second order Volterra model as classification features [3]. The biosensor showed encouraging results, and was able to classify out of sample compounds correctly [2]. We have taken an exploratory step to investigate the advantage of introducing a third order model [4]. DAP5, an NMDA channel blocker, did not show major changes in the second order kernel and in its corresponding Laguerre coefficients. Data were reanalyzed using a third order model. DAP5 showed discernable changes in the third order kernel as well as in the some of the corresponding Laguerre coefficients. Hence, the third order Volterra based model has the potential to improve the sensitivity and the discriminatory power of the proposed bioassay.

A modeling paradigm incorporating parametric and non-parametric methods

A novel parametric/non-parametric modeling paradigm was defined and used in characterization of synaptic transmission. In this paradigm, parametric and nonparametric techniques were incorporated in a complementary manner. Non-parametric method was used to generalize experimental data and extract system input/output properties. It provided a quantitative and intuitive way to validate a parametric model with respect to general, complete input patterns. Biological processes or mechanisms missed by the conventional parametric modeling approach were revealed and subsequently included into the modified parametric model.

17beta-Estradiol potentiates field excitatory postsynaptic potentials within each subfield of the hippocampus with greatest potentiation of the associational/commissural afferents of CA3

We sought to determine the impact of 17beta-estradiol throughout the hippocampal trisynaptic pathway and to investigate the afferent fiber systems within CA1 and CA3 in detail. To achieve this objective, we utilized multielectrode arrays to simultaneously record the field excitatory postsynaptic potentials from the CA1, dentate gyrus, and CA3 of rat hippocampal slices in the presence or absence of 100 pM 17beta-estradiol. We confirmed our earlier findings in CA1, where 17beta-estradiol significantly increased field excitatory postsynaptic potentials amplitude (20%+/-3%) and slope (22%+/-7%). 17beta-Estradiol significantly potentiated the field excitatory postsynaptic potentials in dentate gyrus, amplitude (15%+/-4%) and slope (17%+/-5), and in CA3, amplitude (15%+/-4%) and slope (19%+/-5%). Using a high-density multielectrode array, we sought to determine the source of potentiation in CA1 and CA3 by determining the impact of 17beta-estradiol on the apical afferents and the basal afferents within CA1 and on the mossy fibers and the associational/commissural fibers within CA3. In CA1, 17beta-estradiol induced a modest increase in the amplitude (7%+/-2%) and slope (9%+/-3%) following apical stimulation with similar magnitude of increase following basal stimulation amplitude (10%+/-2%) and slope (12%+/-3%). In CA3, 17beta-estradiol augmented the mossy fiber amplitude (15%+/-3%) and slope (18%+/-6%) and the associational/commissural fiber amplitude (31%+/-13%) and slope (40%+/-15%). These results indicate that 17beta-estradiol potentiated synaptic transmission in each subfield of the hippocampal slice, with the greatest magnitude of potentiation at the associational/commissural fibers in CA3. 17beta-Estradiol regulation of CA3 responses provides a novel site of 17beta-estradiol action that corresponds to the density of estrogen receptors within the hippocampus. The implications of 17beta-estradiol potentiation of the field potential in each of the hippocampal subfields and in particular CA3 associational/commissural fibers for memory function and clinical assessment are discussed.

General methodology for nonlinear modeling of neural systems with Poisson point-process inputs

This paper presents a general methodological framework for the practical modeling of neural systems with point-process inputs (sequences of action potentials or, more broadly, identical events) based on the Volterra and Wiener theories of functional expansions and system identification. The paper clarifies the distinctions between Volterra and Wiener kernels obtained from Poisson point-process inputs. It shows that only the Wiener kernels can be estimated via cross-correlation, but must be defined as zero along the diagonals. The Volterra kernels can be estimated far more accurately (and from shorter data-records) by use of the Laguerre expansion technique adapted to point-process inputs, and they are independent of the mean rate of stimulation (unlike their P-W counterparts that depend on it). The Volterra kernels can also be estimated for broadband point-process inputs that are not Poisson. Useful applications of this modeling approach include cases where we seek to determine (model) the transfer characteristics between one neuronal axon (a point-process 'input') and another axon (a point-process 'output') or some other measure of neuronal activity (a continuous 'output', such as population activity) with which a causal link exists.

Hierarchical model of the population dynamics of hippocampal dentate granule cells

A hierarchical modeling approach is used as the basis for a mathematical representation of the population activity of hippocampal dentate granule cells. Using neural field equations, the variation in time and space of dentate granule cell activity is derived from the summed synaptic potential and summed action potential responses of a population of granule cells evoked by monosynaptic excitatory input from entorhinal cortical afferents. In this formulation of the problem, we have considered a two-level hierarchy: the synapses of entorhinal cortical axons define the first level of organization, and dentate granule cells, which include these synapses, define the second, higher level of organization. The model is specified by two state field variables, for membrane potential and for synaptic efficacy, respectively, with both evolving according to different time scales. The two state field variables introduce new parameters, physiological and anatomical, which characterize the dentate from the point of view of neuronal and synaptic populations: (1) a set of geometrical constraints corresponding to the morphological properties of granule cells and anatomical characteristics of entorhinal-dentate connections; and (2) a set of neuronal parameters corresponding to physiological mechanisms. Assuming no interaction between granule cells, i.e., neither ephaptic nor synaptic coupling, the model is shown to be mathematically tractable and allows solution of the field equations leading to the determination of activity. This treatment leads to the definition of two state variables, volume of stimulated synapses and firing time, which describe observed activity. Numerical simulations are used to investigate the populational characterization of the dentate by individual parameters: (1) the relationship between the conditions of stimulation of active perforant path fibers, e.g., stimulating intensity, and activity in the granule cell layer; and (2) the influence of geometry on the generation of activity, i.e., the influence of neuron density and synaptic density-connectivity. As an example application of the model, the granule cell population spike is reconstructed and compared with experimental data.

Projection of the magnocellular red nucleus to the region of the accessory abducens nucleus in the rabbit

The projection of the magnocellular red nucleus (RNm) to the region of the accessory abducens nucleus (AABD) was traced in rabbit using the bidirectional tracer wheat germ agglutinin-horseradish peroxidase (WGA-HRP). In one set of animals, recordings of antidromic responses from RNm neurons elicited by electrical stimulation of the rubrospinal tract were used to localize injections of WGA-HRP for orthograde labeling of RNm terminals. In a different set of animals, horseradish peroxidase was injected into the retractor bulbi muscle to retrogradely label motoneurons of the AABD. The positions of RNm fibers and terminals were examined and compared to the locations and distribution of AABD cell bodies and labeled dendrites. Analyses revealed that along the entire rostrocaudal extent of the AABD, RNm efferents terminate primarily lateral to, or in the lateral aspects of, labeled motoneurons. For the rostral AABD, RNm efferents terminate only lateral to the nucleus. Although the terminals are not positioned to contact cell bodies of the AABD, they could overlap with dendrites that extend in the lateral direction. RNm efferents terminate more extensively within the posterior AABD, overlapping within both dendritic and cell body regions of the nucleus. Even in this posterior region, however, RNm efferents were distributed primarily over the lateral half of the nucleus. These data show that RNm can monosynaptically influence the AABD, through primarily its lateral and posterior aspects. Our findings also show that a major target of RNm efferents is the reticular cell population located lateral to the AABD, suggesting that the RNm also may affect AABD motoneuronal output indirectly through its projection to reticular cells.

Differential effect of TEA on long-term synaptic modification in hippocampal CA1 and dentate gyrus in vitro

The effectiveness of tetraethylammonium (TEA) and high-frequency stimulation (HFS) in inducing long-term synaptic modification is compared in CA1 and dentate gyrus (DG) in vitro. High-frequency stimulation induces long-term potentiation (LTP) at synapses of both perforant path-DG granule cell and Schaffer collateral-CA1 pyramidal cell pathways. By contrast, TEA (25 mM) induces long-term depression in DG while inducing LTP in CA1. The mechanisms underlying the differential effect of TEA in CA1 and DG were investigated. It was observed that T-type voltage-dependent calcium channel (VDCC) blocker, Ni2+ (50 microM), partially blocked TEA-induced LTP in CA1. A complete blockade of the TEA-induced LTP occurred when Ni2+ was applied together with the NMDA receptor antagonist, D-APV. The L-type VDCC blocker, nifidipine (20 microM), had no effect on CA1 TEA-induced LTP. In DG of the same slice, TEA actually induced long-term depression (LTD) instead of LTP, an effect that was blocked by D-APV. Neither T-type nor L-type VDCC blockade could prevent this LTD. When the calcium concentration in the perfusion medium was increased, TEA induced a weak LTP in DG that was blocked by Ni2+. During exposure to TEA, the magnitude of field EPSPs was increased in both CA1 and DG, but the increase was substantially greater in CA1. Tetraethylammonium application also was associated with a large, late EPSP component in CA1 that persisted even after severing the connections between CA3 and CA1. All of the TEA effects in CA1, however, were dramatically reduced by Ni2+. The results of this study indicate that TEA indirectly acts via both T-type VDCCs and NMDA receptors in CA1 and, as a consequence, induces LTP. By contrast, TEA indirectly acts via only NMDA receptors in DG and results in LTD. The results raise the possibility of a major synaptic difference in the density and/or distribution of T-type VDCCs and NMDA receptors in CA1 and DG of the rat hippocampus.

A novel network for nonlinear modeling of neural systems with arbitrary point-process inputs

This paper address the issue of nonlinear model estimation for neural systems with arbitrary point-process inputs using a novel network that is composed of a pre-processing stage of a Laguerre filter bank followed by a single hidden layer with polynomial activation functions. The nonlinear modeling problem for neural systems has been attempted thus far only with Poisson point-process inputs and using cross-correlation methods to estimate low-order nonlinearities. The specific contribution of this paper is the use of the described novel network to achieve practical estimation of the requisite nonlinear model in the case of arbitrary (i.e. non-Poisson) point-process inputs and high-order nonlinearities. The success of this approach has critical implications for the study of neuronal ensembles, for which nonlinear modeling has been hindered by the requirement of Poisson process inputs and by the presence of high-order nonlinearities. The proposed methodology yields accurate models even for short input-output data records and in the presence of considerable noise. The efficacy of this approach is demonstrated with computer-simulated examples having continuous output and point-process output, and with real data from the dentate gyrus of the hippocampus.

Intraocular pharmacokinetics and safety of a humanized monoclonal antibody in rabbits after intravitreal administration of a solution or a PLGA microsphere formulation

Poly(lactic-co-glycolic) acid (PLGA) bioresorbable microspheres are used for controlled-release drug delivery and are particularly promising for ocular indications. The objective of the current study was to evaluate the pharmacokinetics and safety of a recombinant human monoclonal antibody (rhuMAb HER2) in rabbits after bolus intravitreal administration of a solution or a PLGA-microsphere formulation. On Day 0, forty-eight male New Zealand white rabbits (2.3-2.6 kg) were immobilized with intramuscular ketamine/xylazine, and the test materials were injected directly into the vitreous compartment. Group 1 animals received rhuMAb HER2 in 50:50 lactide: glycolide PLGA microspheres; Group 2 animals received rhuMAb HER2 in solution (n = 24/group). The dose for each eye was 25 microg (50 microl). After dosing, animals were sacrificed at 2 min, and on 1, 2, 4, 7, 14, 23, 29, 37, 44, 50, and 56 days (n = 2/timepoint/group). Safety assessment included direct ophthalmoscopy, clinical observations, body weight, and hematology and clinical chemistry panels. At necropsy, vitreous and plasma were collected for pharmacokinetics and analysis for antibodies to rhuMAb HER2, and the vitreal pellet (Group 1) was prepared for histologic evaluation. All animals completed the study per protocol-both treatments were well tolerated, and no suppurative or mixed inflammatory cell reaction was observed in the vitreal samples (Group 1) at any of the time points examined. Antibodies to rhuMAb HER2 were detected in plasma samples by Day 7 in both treatment groups, but infrequently in vitreous samples. There were no safety implications associated with this immune response. The in vitro characterization of the PLGA microspheres provided reasonable projections of the in vivo rhuMAb HER2 release kinetics (Group 1). The total amount of antibody that was released was similar in vitro (25.9%) and in vivo (32.4%). RhuMAb HER2 (Group 2) was cleared slowly from the vitreous compartment, with initial and terminal half-lives of 0.9 and 5.6 days, respectively. The volume of distribution approximated the vitreous volume in a rabbit eye.

Application of a novel modeling method to the nonstationary properties of potentiation in the rabbit hippocampus

This paper presents the first application of a novel methodology for nonstationary nonlinear modeling to neurobiological data consisting of extracellular population field potentials recorded from the dendritic layer of the dentate gyrus of the rabbit hippocampus under conditions of stimulus-induced potentiation. The experimental stimulus was a Poisson random sequence with a mean rate of 5 impulses/s applied to the perforant path, which was sufficient to induce a progressive potentiation of perforant path-evoked granule cell response. The modeling method utilizes a novel artificial neural network architecture, which is based on the general time-varying Volterra model. The artificial neural network is composed of parallel subnets of three-layer perceptrons with polynomial activation functions, with the output of each subnet modulated by an appropriate time function that models the system nonstationarities and gives the summative output its time-varying characteristics. For the specific application presented herein these time functions are sigmoidal functions with trainable slopes and inflection points. A possible mapping between the nonstationary components of the model and the mechanisms underlying potentiation changes in the hippocampus is discussed.

17beta-estradiol enhances NMDA receptor-mediated EPSPs and long-term potentiation

Gonadal steroid hormones influence CNS functioning through a variety of different mechanisms. To test the hypothesis that estrogen modulates synaptic plasticity in the hippocampus, in vitro hippocampal slices from 2-mo-old Sprague-Dawley male rats were used to determine the effect of 17beta-estradiol on both N-methyl-D-aspartate (NMDA) receptor-mediated excitatory postsynaptic potentials (EPSPs) through intracellular recordings and long-term potentiation (LTP) through extracellular recordings. Intracellular EPSPs and extracellular field EPSPs (fEPSPs) were recorded from CA1 pyramidal cells by stimulating Schaffer collateral fibers. In intracellular experiments, slices were perfused with medium containing bicuculline (5 microM) and low Mg2+ (0.1 mM) to enhance the NMDA receptor-mediated currents and 6, 7-dinitroquinoxaline-2,3-dione (DNQX) (10 microM) to block the alpha-amino-3-hydroxy-5-methyl-4-isoxazoleproprianate (AMPA) receptor-mediated component. The effects of 17beta-estradiol on NMDA receptor-mediated activity were excitatory; concentrations >10 nM induced seizure activity, and lower concentrations (1 nM) markedly increased the amplitude of NMDA-mediated EPSPs (both the first and second responses increased during paired pulse stimulation by 180 and 197%, respectively). In extracellular experiments, slices perfused with 17beta-estradiol (100 pM) exhibited a pronounced, persisting, and significant enhancement of LTP of both the fEPSP slope (192%) and fEPSP amplitude (177%) compared with control slices (fEPSP slope = 155%; fEPSP amplitude = 156%) 30 min after high-frequency stimulation. These data demonstrate that estrogen enhances NMDA receptor-mediated currents and promotes an enhancement of LTP magnitude.

Spatial distribution of potentiated synapses in hippocampus: dependence on cellular mechanisms and network properties

Long-term potentiation (LTP) of synaptic transmission, studied intensively in reduced brain preparations such as hippocampal brain slices, is the leading candidate for the cellular/molecular basis of learning and memory. Serious consideration of LTP as underlying information storage in the intact brain, however, requires understanding how LTP can be induced selectively at specific synaptic sites in a neural system when the mechanisms underlying LTP are regulated by other structural and functional properties of the same neural system. In the studies reported here, we tested the hypothesis that different patterns of activity within the same population of entorhinal cortical afferents could lead to a selective potentiation of spatially distinct populations of synapses across different regions of the hippocampus, including those activated multisynaptically. We focused specifically on potentiation of direct, monosynaptic entorhinal input to dentate granule cells, which expresses an NMDA receptor-dependent LTP, and on potentiation of indirect, disynaptic entorhinal input to CA3 pyramidal cells, which is transmitted by the mossy fiber projection of dentate granule cells and expresses an NMDA receptor-independent LTP. The principal findings of these experiments show that lower stimulation frequencies (10-20 Hz) of entorhinal cortical axons selectively induce LTP of mossy fiber input to CA3 transsynaptically via excitation of dentate granule cells, and that patterns of stimulation of that mimic neuronal firing in the entorhinal cortex during endogenous theta rhythm (five-impulse bursts at 200 Hz, interburst intervals of 200 msec) induce LTP both monosynaptically for input to dentate granule cells and transsynaptically for mossy fiber input to CA3.

Novel expression mechanism for synaptic potentiation: alignment of presynaptic release site and postsynaptic receptor

A combination of experimental and modeling approaches was used to study cellular-molecular mechanisms underlying the expression of short-term potentiation (STP) and long-term potentiation (LTP) of glutamatergic synaptic transmission in the hippocampal slice. Electrophysiological recordings from dentate granule cells revealed that high-frequency stimulation of perforant path afferents induced a robust STP and LTP of both (+/-)-alpha-amino-3-hydroxy-5-methylisoxazole-4-propionic acid (AMPA) and N-methyl-D-aspartic acid (NMDA) receptor-mediated synaptic responses. However, the decay time constant for STP of the AMPA receptor-mediated excitatory postsynaptic potential was approximately 6 min, whereas the decay time constant for STP of the NMDA receptor-mediated excitatory postsynaptic potential was only 1 min. In addition, focal application of agonists during the expression of STP revealed that the magnitude of conductance change elicited by NMDA application was significantly enhanced, whereas the magnitude of conductance change elicited by application of AMPA remained constant. These findings are most consistent with a postsynaptic mechanism of STP and LTP. Different putative mechanisms were evaluated formally using a computational model that included diffusion of glutamate within the synaptic cleft, different kinetic properties of AMPA and NMDA receptor/channels, and geometric relations between presynaptic release sites and postsynaptic receptor/channels. Simulation results revealed that the only hypothesis consistent with experimental data is that STP and LTP reflect a relocation of AMPA receptor/channels in the postsynaptic membrane such that they become more closely "aligned" with presynaptic release sites. The same mechanism cannot account for STP or LTP of NMDA receptor-mediated responses; instead, potentiation of the NMDA receptor subtype is most consistent with an increase in receptor sensitivity or number.

Transport and fate of trifluoroacetate in upland forest and wetland ecosystems

Although trifluoroacetate (TFA), a breakdown product of chlorofluorocarbon replacements, is being dispersed widely within the biosphere, its ecological fate is largely unknown. TFA was added experimentally to an upland, northern hardwood forest and to a small forest wetland ecosystem within the Hubbard Brook Experimental Forest in New Hampshire. Inputs of TFA were not transported conservatively through these ecosystems; instead, significant amounts of TFA were retained within the vegetation and soil compartments. More TFA was retained by the wetland ecosystem than by the upland forest ecosystem. Using simulation modeling, TFA concentrations were predicted for soil and drainage water until the year 2040.

Differential expression of short-term potentiation by AMPA and NMDA receptors in dentate gyrus

Both alpha-amino-3-hydroxy-5-methyl-4-isoxazole propionate (AMPA) and N-methyl-D-aspartate (NMDA) glutamatergic receptor subtypes in hippocampus have been shown to express long-term potentiation (LTP), a form of synaptic modification believed to be involved in memory formation. Because of their postsynaptic localization, any differential expression of LTP by the two receptor subtypes would strongly support the existence of a postsynaptic mechanism of LTP expression. In this study, electrophysiological recordings from dentate granule cells were used to compare the potentiation of AMPA and NMDA receptor-mediated responses occurring during the initial phase of LTP, typically identified as STP. Results revealed that high-frequency stimulation (HFS) of perforant path afferents induces a robust STP of both AMPA and NMDA receptor-mediated components of granule cell EPSPs (referred to as AMPA STP and NMDA STP, respectively). Although STP for both receptor subtypes decayed to an asymptotic, steady-state level of LTP and could be induced repetitively, there were substantial differences in several aspects of AMPA and NMDA STP dynamics. STP of the AMPA receptor reached its peak magnitude approximately 30 sec after HFS and decayed with a time constant of approximately 6 min. In contrast, peak magnitude of NMDA STP always appeared immediately after HFS and decayed with a time constant of only 1 min. Single-pulse stimulation of perforant path afferents paired with postsynaptic depolarization also induced LTP of both AMPA and NMDA components. When this induction paradigm was used, however, only the AMPA component showed significant STP. Our results demonstrate that AMPA and NMDA receptors exhibit markedly different degrees of activity-dependent, short-term modifiability, with the possibility that STP of the NMDA receptor reflects primarily post-tetanic potentiation (PTP). In addition, our results strongly suggest that the mechanisms underlying STP of the AMPA receptor are postsynaptic in origin.

Dynamic synapse: a new concept of neural representation and computation

Presynaptic mechanisms influencing the probability of neurotransmitter release from an axon terminal, such as facilitation, augmentation, and presynaptic feedback inhibition, are fundamental features of biological neurons and are cardinal physiological properties of synaptic connections in the hippocampus. The consequence of these presynaptic mechanisms is that the probability of release becomes a function of the temporal pattern of action potential occurrence, and hence, the strength of a given synapse varies upon the arrival of each action potential invading the terminal region. From the perspective of neural information processing, the capability of dynamically tuning the synaptic strength as a function of the level of neuronal activation gives rise to a significant representational and processing power of temporal spike patterns at the synaptic level. Furthermore, there is an exponential growth in such computational power when the specific dynamics of presynaptic mechanisms varies quantitatively across axon terminals of a single neuron, a recently established characteristic of hippocampal synapses. During learning, alterations in the presynaptic mechanisms lead to different pattern transformation functions, whereas changes in the postsynaptic mechanisms determine how the synaptic signals are to be combined. We demonstrate the computational capability of dynamic synapses by performing speech recognition from unprocessed, noisy raw waveforms of words spoken by multiple speakers with a simple neural network consisting of a small number of neurons connected with synapses incorporating dynamically determined probability of release. The dynamics included in the model are consistent with available experimental data on hippocampal neurons in that parameter values were chosen so as to be consistent with time constants of facilitative and inhibitory processes governing the dynamics of hippocampal synaptic transmission studied using nonlinear systems analytic procedures.

NMDA receptor-dependent LTD in different subfields of hippocampus in vivo and in vitro

In simulations with artificial neural networks, efficient information processing and storage has been shown to require that the strength of connections between network elements has the capacity to both increase and decrease in a use-dependent manner. In contrast to long-term potentiation (LTP) of excitatory synaptic transmission, activity-dependent long-term depression (LTD) has been difficult to demonstrate in forebrain in vivo. Theoretical arguments indicate that coincidence of presynaptic excitation and low-magnitude postsynaptic activation are the necessary prerequisites for LTD induction. Here we report that stimulation paradigms which cause 1) sufficient excitation to result in NMDA receptor activation and simultaneously 2) attenuate the level of postsynaptic activation by recruitment of GABAA receptor-mediated inhibition consistently produce LTD of commissural input to area CA1 in the hippocampus of anesthetized adult rats, and of the perforant path input to the dentate gyrus in the hippocampus of anesthetized and unanesthetized adult rabbits. A functionally similar pre- and postsynaptic activation pattern applied to the hippocampal slice preparation by injecting hyperpolarizing current into the postsynaptic cell during NMDA receptor-mediated excitation also was effective in consistently inducing LTD. Results of studies in vitro show that Ca2+ influx through the NMDA channel is necessary for the induction of LTD, and moreover, that NMDA receptors also participate in the expression of LTD. Our findings demonstrate a general mechanism for the implementation of a theoretically derived learning rule in adult forebrain in vivo and in vitro and provide justification for the inclusion of use-dependent decreases of connection weights in formal models of cognitive processing.

Excitatory stimulation during postsynaptic inhibition induces long-term depression in hippocampus in vivo

1. As part of an effort to evaluate the biological plausibility of theoretically derived principles of synaptic modification, we studied activity-dependent long-term depression (LTD) of glutamatergic transmission in the hippocampus of anesthetized adult rats. Field potentials of CA1 pyramidal cells evoked by single-pulse stimulation (0.1 Hz) of the commissural afferents were recorded before and after paired-pulse stimulation (0.5 Hz) of the same pathway. A train of 150 or 200 paired pulses produced robust LTD of the commissural input to the CA1 pyramidal neurons when the interstimulus interval (ISI) of the pairs was short (25 ms) but not when the ISI was long (1,000 ms). 2. Paired-pulse stimulation with the short but not with the long ISI also was associated with pronounced inhibition of pyramidal cell firing upon the second pulse of a pair, despite the fact that the excitatory input was facilitated with the short-ISI paradigm. The inhibition of pyramidal cell activity was mediated by input to the pyramidal cells from local gamma-aminobutyric acid (GABA)-releasing interneurons activated by commissural fibers and/or CA1 recurrent collaterals, because the inhibition was eliminated by local administration of the selective GABAA receptor antagonist, bicuculline (50 microM), near the recording site. 3. Postsynaptic input from GABAergic interneurons was necessary for the induction of LTD, because short-ISI paired-pulse stimulation failed to produce LTD in the presence of bicuculline. 4. N-methyl-D-aspartate (NMDA) receptor-mediated excitation also was necessary for the induction of LTD, because administration of the selective NMDA receptor antagonist, D-2-amino-5-phosphonvaleric acid (100 microM), near the recording site prevented the development of LTD.(ABSTRACT TRUNCATED AT 250 WORDS)

Identification and characterization of striatal cell subtypes using in vivo intracellular recording and dye-labeling in rats: III. Morphological correlates and compartmental localization

In the first two reports of this series, in vivo intracellular recording techniques were used to characterize the electrophysiological properties of two types of striatal neurons that had been identified by their distinct response patterns to stimulation of corticostriatal afferents. In this paper, we examined whether cells showing Type I or Type II response patterns also differed with respect to their morphology or compartmental localization by combining intracellular recording and Lucifer yellow staining with immunocytochemical localization of calbindin 28 kd immunoreactivity. In the majority of cases, both Type I and Type II neurons exhibited similar morphological characteristics, with 80% of the Type I cells (13/16) and all of the Type II cells (n = 40) being small or medium spiny neurons. In each case where the morphological class of the cell was different than the spiny cell class, the cell exhibited a Type I response pattern. These spiny neurons had somata that averaged 17.1 +/- 1.3 microns in diameter and gave rise to between four and eight primary dendrites. The axons typically arose from cell bodies (7/13 for Type I and 25/40 for Type II cells) and emitted extensive local axonal collaterals. However, the axons of Type I cells more frequently originated from the dorsal surface of the somata (9/13; 69%), whereas Type II axons more frequently arose from the ventral surface of the somata (25/35; 71%), which may account for their different extracellular waveforms. In coronally sectioned tissue (n = 18), the axons always projected laterally when the somata were located in the medial striatum and projected medially when the somata were in the lateral striatal region. In a subset of experiments (N = 22), Lucifer yellow-stained neurons were localized with respect to their position within the patch and matrix compartments of the striatum using subsequent staining for calbindin 28 kd immunoreactivity. Of the 20 labeled medium spiny neurons examined (Type II: N = 13; Type I: N = 7), 19 were located in the calbindin-positive matrix compartment. The only neuron localized to the patch compartment was a medium spiny cell that exhibited a Type II paired impulse response pattern. In addition, of the two aspiny neurons from this group with beaded dendrites, one was localized to the border between adjacent patch and matrix compartments, whereas the other was located completely within the matrix compartment. Therefore, despite their distinct paired impulse response patterns, the majority of both Type I and Type II neurons were medium spiny cells located in the matrix compartment of the striatum.(ABSTRACT TRUNCATED AT 400 WORDS)

Identification and characterization of striatal cell subtypes using in vivo intracellular recording in rats: I. Basic physiology and response to corticostriatal fiber stimulation

The electrophysiological characteristics of two subtypes of striatal neurons, identified by their distinct patterns of response to paired impulse stimulation of corticostriatal afferents, were compared using in vivo intracellular recordings in rats. As observed in previous extracellular recording studies, the majority of neurons (73%) were found to be of the Type II class, with the remaining cells exhibiting the Type I response pattern. For all cells, cortical stimulation elicited 5-30 mV EPSPs at latencies ranging from 2.0-5.3 msec. Increasing the stimulating current intensity caused a progressive increase in the amplitude of the evoked EPSPs without altering their latencies, suggesting that the EPSPs are monosynaptically mediated. Both the average amplitude and duration of the evoked EPSPs at spike threshold in Type I neurons (9.8 +/- 1.7 mV, 11.8 +/- 2.8 msec; mean +/- SEM) were significantly smaller than those of Type II cells (20.3 +/- 1.4 mV, 22.7 +/- 2.1 msec). Although the average latency to the onset of the EPSP was similar for both cell classes (Type I cells: 2.3 +/- 0.3 msec; Type II cells: 2.2 +/- 0.2 msec), the EPSPs in Type I cells reached peak amplitude and the spikes were triggered at significantly longer latencies than in the Type II cells (peak I: 11.2 +/- 2.5 msec vs. II: 7.6 +/- 0.7 msec; spike I: 8.0 +/- 1.2 msec vs. II: 5.7 +/- 0.4 msec). Striatal neurons had a comparatively hyperpolarized resting membrane potential (-70.2 +/- 2.1 mV) and had an average input resistance of 35.4 +/- 7.6 M omega. Overall, striatal neurons exhibited low levels of spontaneous activity (0.6 +/- 0.7 Hz) with 50% of the neurons being quiescent. Type I cells exhibited significantly higher firing rates (3.2 +/- 0.8 Hz) than Type II cells (0.8 +/- 0.3 Hz), although their resting membrane potentials were not significantly different. Spontaneously occurring spikes had an average amplitude of 72.7 +/- 3.4 mV and spike thresholds of -50.1 +/- 1.5 mV. Irregularly occurring depolarizing plateau potentials, which typically gave rise to spike discharge, were frequently observed in both spontaneously firing and quiescent neurons. A small proportion of the cells recorded (3/55) exhibited a Type I response pattern but had unique physiological characteristics that were similar to those described by others as arising from large, aspiny striatal neurons.(ABSTRACT TRUNCATED AT 400 WORDS)

Identification and characterization of striatal cell subtypes using in vivo intracellular recording in rats: II. Membrane factors underlying paired-pulse response profiles

Two subtypes of neurons in the striatum have been defined on the basis of their different response patterns to paired-impulse stimulation of corticostriatal afferents, with type I cells showing a longer spike latency, facilitation at short intervals, and inhibition at long intervals, and type II cells defined by the facilitation occurring at long interstimulus intervals. Nevertheless, the companion report has shown that this distinction of cell types cannot be accounted for by differences in the basic physiological properties of these cells, but instead is likely to be due to differences in their synaptic connectivity. The experiments performed in this study were directed at examining in detail the membrane factors and synaptic responses that may contribute to these distinct response patterns. When pairs of stimuli were delivered to the corticostriatal fibers at 10-30 ms interstimulus intervals, the EPSPs elicited in type I neurons exhibited a temporal summation, resulting in a facilitation of spike firing to the second stimulus relative to the first. In contrast, type II cells showed decreased EPSP amplitude at short intervals, and in cells showing a short-interval inhibition, there was a significant increase in spike threshold (+5.3 +/- 1.4 mV) during the second response. All type II neurons recorded with KCl-filled microelectrodes showed short-interval facilitation with little or no change in spike threshold. Although the use of KCl electrodes did not alter the facilitation at short intervals in type I neurons it did increase the rate of rise of the EPSP, causing spikes to be triggered at a latency similar to that of type II cells. Paired stimuli delivered at 75-150 ms interstimulus intervals showed inverse effects on type I and type II cells with respect to the probability of spike firing. In type I cells, the evoked EPSP was followed by a long-latency membrane hyperpolarization that prevented the second EPSP from reaching spike threshold. In contrast, the smaller-amplitude hyperpolarization evoked in type II cells enabled the second stimulus to activate an EPSP at the same membrane potential as the first stimulus, resulting in a facilitation of spiking.(ABSTRACT TRUNCATED AT 400 WORDS)

Na+ influx through Ca2+ channels can promote striatal GABA efflux in Ca(2+)-deficient conditions in response to electrical field depolarization

Electrical field depolarization releases gamma-aminobutyric acid (GABA) in rat striatal slices in the absence of external Ca2+. omega-Conotoxin GVIA (omega-CgTx; 1-50 nM), a neuronal Ca2+ channel blocker, inhibits electrically evoked efflux of newly taken up [3H]GABA in a concentration-dependent manner in either normal or Ca(2+)-free medium. This suggests that ion influx occurs through Ca2+ channels in the absence of external Ca2+ and contributes to the efflux of GABA. Reducing external Na+ concentration to 27.25 mM (low [Na+]o medium) by equimolarly substituting choline chloride for sodium chloride has differential effects on electrically evoked GABA efflux depending on the external Ca2+ concentrations. In normal Ca2+ medium, electrically evoked GABA efflux increases whereas, in Ca(2+)-free medium, it is greatly inhibited when [Na+]o is reduced to 27.25 mM. In low [Na+]o medium, GABA efflux is largely tetrodotoxin (TTX)-sensitive, however, spike firing evoked by antidromic stimulation of striatal cells is inhibited. In Na(+)-free medium, resting GABA efflux increases 17-fold whereas evoked GABA efflux diminishes. In Ca(2+)-free medium, 70 min of incubation with 1-2-bis-(1-aminophenoxy)ethane-N,N,N',N' tetraacetoxy methyl ester (BATPA-AM, 1 microM), an intracellular calcium chelator, increases both resting GABA efflux and electrically evoked GABA overflow by approximately 100%. These results suggest that: (1) in Ca(2+)-free conditions, Na+ permeability of cells increases via Ca2+ channels and this profoundly affects GABA efflux. (2) Electrical field depolarization is likely to release GABA by directly depolarizing axon terminals. (3) Ca(2+)-independent GABA efflux is not promoted by an increase in intracellular free Ca2+ concentration via Na+/Ca2+ exchange processes from internal pools.

Depression of glutamatergic and GABAergic synaptic responses in striatal spiny neurons by stimulation of presynaptic GABAB receptors

The influence of gamma-aminobutyric acidB (GABAB) receptor stimulation on the excitatory and inhibitory synaptic potentials and membrane properties of identified striatal spiny neurons was examined in a corticostriatal slice preparation. Stimulation of the subcortical white matter evoked a monosynaptic, excitatory postsynaptic potential (EPSP) and a polysynaptic, inhibitory postsynaptic potential (IPSP) in spiny neurons. The EPSP had two components: a large amplitude response which could be blocked by the kainate/quisqualate receptor antagonist, 6-cyano-7-nitroquinoxaline-2,3-dione (CNQX, 10 microM), and a small amplitude, long-duration depolarization which could be blocked by the N-methyl-D-aspartate receptor antagonist, d-(-)-2-amino-5-phosphonovaleric acid (APV, 100 microM). The IPSP was observed as a membrane depolarization when recorded from neurons at resting membrane potential. However, when neurons were injected with the Na(+)-channel blocker, QX-314, allowing cells to be depolarized above their spike thresholds, a prominent hyperpolarizing IPSP was readily observed which could be blocked by the GABAA antagonist, bicuculline (10-50 microM). This bicuculline-sensitive IPSP was responsible for the inhibition of EPSP amplitude and probability of spike discharge revealed using paired stimulation of the subcortical white matter. The amplitude of both the EPSP and the IPSP were depressed by the GABAB receptor agonist, p-chlophenyl-GABA (baclofen, 0.5-100 microM) in a concentration-dependent manner. Baclofen also blocked paired stimulus inhibition of spike discharge. These effects of baclofen persisted in slices in which the cortex was removed and were reversed by the GABAB receptor antagonist, 3-amino-3-hydroxy-2-(4-chlorophenyl)-propanesulphonic acid (saclofen, 100-500 microM). In contrast to its profound influence on synaptic input, baclofen did not alter resting membrane potential, input resistance, membrane current-voltage relationship, or spike threshold of the cells recorded, and therefore did not appear to exert direct postsynaptic effects on the striatal spiny neurons. Taken together, these data indicate that the depressant effects of baclofen on EPSPs are mediated through GABAB receptors located on the terminals of glutamatergic afferents within the striatum. Moreover, the results suggest that the actions of baclofen on IPSPs and paired stimulus inhibition are produced by activation of GABAB receptors within the striatum at a site presynaptic to spiny neurons, either on the terminals of GABAergic afferents or on an interposed non-spiny GABAergic cell. Thus, GABAB hetero- and auto-receptors have the capacity to provide a negative feedback mechanism through which the major excitatory and inhibitory inputs to striatal spiny neurons are regulated.

Recovery of hippocampal dentate granule cell responsiveness to entorhinal cortical input following norepinephrine depletion

Hippocampal dentate granule cell responsivity to excitatory input from entorhinal perforant path fibers was examined in the chronic rabbit preparation following norepinephrine (NE) depletion induced with the neurotoxin DSP4. To examine granule cell responsivity as a function of perforant path activation, constant low frequency stimulation (0.1 Hz) was applied to the perforant path using an ascending intensity series. To examine granule cell responsivity to more complex patterns of stimulation, a train of impulses, with a random interstimulus interval (Poisson distribution; mean frequency of 2 Hz), was applied to the perforant path. Both single impulse and random interval impulse stimulation revealed that NE depletion increased the average amplitude of the perforant path-granule cell population spike. The random interval impulse stimulation revealed that NE depletion also increased the magnitude and duration of second order inhibitory interactions. These changes were transient, however, and recovered over the 21 day test period. Hippocampal NE levels were reduced an average of 80% between 23 and 38 days post-DSP4. The activity of the rate-limiting enzyme for NE synthesis, tyrosine hydroxylase (TH), was reduced an average of 60%. That NE levels were reduced to a greater extent than was TH activity is suggestive of increased NE synthesis within the remaining nerve terminals. Such an increase in NE synthesis may reflect a compensatory response underlying the functional recovery of electrophysiological responsiveness following partial NE depletion.

Characterizationin vivo of the NMDA receptor-mediated component of dentate granule cell population synaptic responses to perforant path input

The NMDA receptor-mediated component of the hippocampal granule cell population excitatory postsynaptic potential response to low frequency (< 0.2 Hz) stimulation of the medial perforant path was characterized in vivo. Extracellular recordings were obtained from the dentate molecular layer in anesthetized rabbits, and glutamatergic and GABAergic antagonists were applied locally by pressure ejection. To measure the NMDA-mediated component, the NMDA receptor antagonist D-5-aminophosphonovalerate (APV) was applied during the constant ejection of physiological saline, 6-cyano-7-nitroquinoxaline-2,3-dione (CNQX), and/or bicuculline methiodide. In general agreement with the results of attempts by other investigators to identify NMDA responses in vivo, APV did not significantly reduce the response to a single stimulus impulse in the presence of saline. However, an NMDA-mediated response was revealed when alpha-amino-3-hydroxy-5-methyl-4-isoxazoleproprianate receptor-mediated current flow was eliminated by applying the non-NMDA receptor antagonist CNQX. The NMDA component was negative-going as predicted, but its duration was considerably less than indicated in other studies of the dentate in vitro. The relative magnitudes of the NMDA and non-NMDA components of the EPSP were found to vary as a function of stimulus intensity or frequency. The NMDA receptor-mediated component represented 12% of the control response and increased to over 25% in response to higher stimulus intensities. A brief, high-frequency burst of impulses evoked a larger NMDA component in the presence of CNQX and was able to evoke an NMDA component in the presence of saline. Surprisingly, short trains of stimulation at lower frequencies typically produced suppression of the NMDA component. In a final series of experiments, it was found that many characteristics of the NMDA component were substantially altered by GABAergic inhibition. In the presence of the GABAA antagonist bicuculline, the magnitude of NMDA receptor-mediated responses was increased and their duration was greatly extended. Additionally, in the presence of bicuculline, the NMDA component facilitated markedly in response to frequencies of stimulus input > 20 Hz. These results indicate in vivo that the initiation and duration of NMDA current flow depend strongly upon the intensity and frequency of perforant path stimulation. In addition, the NMDA response to a single impulse appears to be reduced and truncated by input from GABAA receptor-mediated feedback and/or feedforward inhibition, and this inhibition affects temporal summation of NMDA receptor-mediated responses over a wide range of input frequencies. It is suggested that such inhibition results from the activation of GABAA receptors located on granule cell dendritic shafts.(ABSTRACT TRUNCATED AT 400 WORDS)

Functionally distinct subpopulations of striatal neurons are differentially regulated by GABAergic and dopaminergic inputs--II. In vitro analysis

In the companion report [Nisenbaum and Berger (1992) Neuroscience 48, 561-578] the contrasting paired impulse responses to stimulation of the corticostriatal pathway which define the Type I and Type II subpopulations of striatal neurons were shown to reflect differential regulation by GABAergic and dopaminergic inputs. More specifically, the decreased probability of spike discharge (inhibition) to long interstimulus intervals (60-260 ms) characteristic of Type I neurons was found to be dependent on dopaminergic input via D1 receptor activation, whereas the inhibition to short interstimulus intervals (10-20 ms) distinctive of Type II neurons was found to be mediated by GABAergic input acting through GABAA receptor stimulation. The present experiments have further investigated the contribution of GABAergic and dopaminergic feedforward and/or feedback circuits to the functional identities of Type I and Type II neurons using an in vitro corticostriatal slice preparation. In this preparation, the cortical afferents to the striatum are preserved, allowing for activation of striatal cells in a manner similar to that used in vivo; however, all axons arising from midbrain and brainstem structures including the substantia nigra are transected, and intrastriatal GABAergic pathways are reduced. Consistent with the predicted effect of disrupting these two neurotransmitter pathways, the paired impulse responses of striatal neurons recorded in vitro were not similar to the responses of either Type I or Type II neurons recorded in vivo. Indeed, the paired impulse profiles of striatal neurons recorded in vitro were relatively homogeneous in that virtually all cells displayed an increased probability of spike discharge (facilitation) to the second impulse of all interstimulus intervals (10-500ms) tested. Low concentrations of allosteric agonists for the GABAA receptor, pregnanolone (5 microM) and pentobarbital (50 microM), selectively inhibited spike discharge in response to short interstimulus intervals (10-20 ms) for approximately 40% of the neurons sampled, but produced no change in facilitation to longer interstimulus intervals (30-500 ms). The agonist-induced inhibition to short interstimulus intervals was blocked by bicuculline (10-20 microM), and was not mimicked by the GABAB receptor agonist, baclofen (1-5 microM). In addition, application of dopamine (5-10 microM) or the D1 receptor agonist, SKF38393 (2,3,4,5-tetrahydro-7,8-dihydroxy-1-phenyl-1H-3-benzazepine; 5 microM), inhibited spike discharge to longer interstimulus intervals (40-500 ms) for approximately 10% of striatal cells recorded. The inhibition to longer interstimulus intervals was blocked by the D1 receptor antagonist, SCH23390 [R-(+)-8-chloro-2,3,4,5-tetrahydro-3-methyl-5-phenyl-1H-3-benzazepin+ ++-7-ol], but not the D2 antagonist, sulpiride.(ABSTRACT TRUNCATED AT 400 WORDS)