High-sensitivity rod photoreceptor input to the blue-yellow color opponent pathway in macaque retina.

Small bistratified cells (SBCs) in the primate retina carry a major blue-yellow opponent signal to the brain. We found that SBCs also carry signals from rod photoreceptors, with the same sign as S cone input. SBCs exhibited robust responses under low scotopic conditions. Physiological and anatomical experiments indicated that this rod input arose from the AII amacrine cell-mediated rod pathway. Rod and cone signals were both present in SBCs at mesopic light levels. These findings have three implications. First, more retinal circuits may multiplex rod and cone signals than were previously thought to, efficiently exploiting the limited number of optic nerve fibers. Second, signals from AII amacrine cells may diverge to most or all of the approximately 20 retinal ganglion cell types in the peripheral primate retina. Third, rod input to SBCs may be the substrate for behavioral biases toward perception of blue at mesopic light levels.

Rainfall input, throughfall and stemflow of nutrients in a central African rain forest dominated by ectomycorrhizal trees

Incident rainfall is a major source of nutrient input to a forest ecosystem and the consequent throughfall and stemflow contribute to nutrient cycling. These rain-based fluxes were measured over 12 mo in two forest types in Korup National Park, Cameroon, one with low (LEM) and one with high (HEM) ectomycorrhizal abundances of trees. Throughfall was 96.6 and 92.4% of the incident annual rainfall (5370 mm) in LEM and HEM forests respectively; stemflow was correspondingly 1.5 and 2.2%. Architectural analysis showed that ln(funneling ratio) declined linearly with increasing ln(basal area) of trees. Mean annual inputs of N, P, K, Mg and Ca in incident rainfall were 1.50, 1.07, 7.77, 5.25 and 9.27 kg ha(-1), and total rain-based inputs to the forest floor were 5.0, 3.2, 123.4, 14.4 and 37.7 kg ha-1 respectively. The value for K is high for tropical forests and that for N is low. Nitrogen showed a significantly lower loading of throughfall and stemflow in HEM than in LEM forest, this being associated in the HEM forest with a greater abundance of epiphytic bryophytes which may absorb more N. Incident rainfall provided c. 35% of the gross input of P to the forest floor (i. e., rain-based plus small litter inputs), a surprisingly high contribution given the sandy P-poor soils. At the start of the wet season leaching of K from the canopy was particularly high. Calcium in the rain was also highest at this time, most likely due to washing off of dry-deposited Harmattan dusts. It is proposed that throughfall has an important `priming' function in the rapid decomposition of litter and mineralization of P at the start of the wet season. The contribution of P inputted from the atmosphere appears to be significant when compared to the rates of P mineralization from leaf litter.

Helium transfer from water into quartz crystals: A new approach for porewater dating

Several important fundamental and applied problems require a quantification of slow rates of groundwater flow. To resolve these problems helium appears to be a promising tracer. In this contribution we discuss a new approach, which gives the helium inventory in a rock – pore water system by using the relevant mineral record, i.e., without extraction and investigation of the porewater samples. Some U- and Th-poor minerals such as quartz (quartz separates from Permo-Carboniferous Formation, sandstone–shale interlayering, Molasses Basin, Northern Switzerland, hereafter PCF, are used in this study) contain excessive helium having migrated into their internal helium-accessible volume (HAV) from the surrounding porewater [I.N. Tolstikhin, B.E. Lehmann, H.H. Loosli, A. Gautschi, Helium and argon isotopes in rocks, minerals and related groundwaters: a case study in Northern Switzerland, Geochim. Cosmochim. Acta 60 (1996) 1497–1514]. These volumes are estimated by using helium as a nano-size penetrating tool, i.e., by saturation of the minerals with helium under controlled pressure–temperature conditions and subsequent measurements of the helium-saturated concentrations. In the quartz separates HAV/total volume ratios vary from 0.017% to 0.16%; along with the measured initial (unsaturated) He concentration the HAV gives the internal helium pressure, the mean value obtained for 7 samples (25 sample aliquots) is P=0.45F0.15 atm (1 r). The product of helium pressure and solubility (7.35_10_3 cc STP He/cc H2O for the temperature and salinity of PCF aquifers reported in [F.J. Pearson, W. Balderer, H.H. Loosli, B.E. Lehmann, A. Matter, T. Peters, H. Schmassmann, A. Gautschi, Applied Isotope Hydrogeology–A Case Study in Northern Switzerland, Elsevier Amsterdam, 1991, 439 pp.]) is the mineral-derived He concentration in the respective porewater, CPW=0.0035F0.0017 cc He/cc H2O. This value is in full accord with measured He concentrations in PCF aquifers, CPCF, varying from 0.0045 to 0.0016 cc He/cc H2O. This agreement validates the proposed approach and also shows that the mineral–porewater helium–concentration equilibrium has been established. Indeed, estimates of the He-migration rates through our quartz samples show that in ~6000 years the internal pressure should equilibrate with He-concentration in related porewater of PCF, and this time interval is short compared to independent estimates [I.N. Tolstikhin, B.E. Lehmann, H.H. Loosli, A. Gautschi, Helium and argon isotopes in rocks, minerals and related groundwaters: a case study in Northern Switzerland, Geochim. Cosmochim. Acta 60 (1996) 1497–1514]. The helium inventory in the rock–porewater assemblage shows that helium abundance in pore waters is indeed important. In shale samples (with ~15% porosity) porewaters contain more helium than the host minerals altogether. Porewater heliumconcentration profiles, available from the mineral record, along with helium production rates are input parameters allowing model(s) of helium migration through a hydrological structure to be developed. Quite high helium concentrations in PCF porewaters imply slow removal mechanisms, which will be discussed elsewhere.

Local environmental knowledge: a key input for the sustainable development of Andean mountain communities

Local knowledge is crucial to both human development and environmental conservation. This is especially the case in mountain regions, where a combination of remoteness, harsh climatic conditions, rich cultural heritage, and high biological diversity has led to the development of complex local environmental knowledge systems. In the Andes for instance, rural populations mainly rely on their own environmental knowledge to ensure their food security and health. Recent studies conducted within Quechua communities in Peru and Bolivia showed that this knowledge was both persistent and dynamic, and that it responded to socio-economic and environmental changes through cultural resistance and adaptation. As this paper argues, combining local knowledge and so-called scientific knowledge – especially in development projects – can lead to innovative solutions to the socio-environmental challenges facing mountain communities in our globalized world. Based on experiences from the Andes, this paper will provide concrete recommendations to policymakers and practitioners for integrating local knowledge into development and natural resource management initiatives.

Comparison of (31)P saturation and inversion magnetization transfer in human liver and skeletal muscle using a clinical MR system and surface coils.

(31)P MRS magnetization transfer ((31)P-MT) experiments allow the estimation of exchange rates of biochemical reactions, such as the creatine kinase equilibrium and adenosine triphosphate (ATP) synthesis. Although various (31)P-MT methods have been successfully used on isolated organs or animals, their application on humans in clinical scanners poses specific challenges. This study compared two major (31)P-MT methods on a clinical MR system using heteronuclear surface coils. Although saturation transfer (ST) is the most commonly used (31)P-MT method, sequences such as inversion transfer (IT) with short pulses might be better suited for the specific hardware and software limitations of a clinical scanner. In addition, small NMR-undetectable metabolite pools can transfer MT to NMR-visible pools during long saturation pulses, which is prevented with short pulses. (31)P-MT sequences were adapted for limited pulse length, for heteronuclear transmit-receive surface coils with inhomogeneous B1 , for the need for volume selection and for the inherently low signal-to-noise ratio (SNR) on a clinical 3-T MR system. The ST and IT sequences were applied to skeletal muscle and liver in 10 healthy volunteers. Monte-Carlo simulations were used to evaluate the behavior of the IT measurements with increasing imperfections. In skeletal muscle of the thigh, ATP synthesis resulted in forward reaction constants (k) of 0.074 ± 0.022 s(-1) (ST) and 0.137 ± 0.042 s(-1) (IT), whereas the creatine kinase reaction yielded 0.459 ± 0.089 s(-1) (IT). In the liver, ATP synthesis resulted in k = 0.267 ± 0.106 s(-1) (ST), whereas the IT experiment yielded no consistent results. ST results were close to literature values; however, the IT results were either much larger than the corresponding ST values and/or were widely scattered. To summarize, ST and IT experiments can both be implemented on a clinical body scanner with heteronuclear transmit-receive surface coils; however, ST results are much more robust against experimental imperfections than the current implementation of IT.

A Baseline Wander Tracking System for Artifact Rejection in Long-Term Electrocardiography

Long-term electrocardiogram (ECG) signals might suffer from relevant baseline disturbances during physical activity. Motion artifacts in particular are more pronounced with dry surface or esophageal electrodes which are dedicated to prolonged ECG recording. In this paper we present a method called baseline wander tracking (BWT) that tracks and rejects strong baseline disturbances and avoids concurrent saturation of the analog front-end. The proposed algorithm shifts the baseline level of the ECG signal to the middle of the dynamic input range. Due to the fast offset shifts, that produce much steeper signal portions than the normal ECG waves, the true ECG signal can be reconstructed offline and filtered using computationally intensive algorithms. Based on Monte Carlo simulations we observed reconstruction errors mainly caused by the non-linearity inaccuracies of the DAC. However, the signal to error ratio of the BWT is higher compared to an analog front-end featuring a dynamic input ranges above 15 mV if a synthetic ECG signal was used. The BWT is additionally able to suppress (electrode) offset potentials without introducing long transients. Due to its structural simplicity, memory efficiency and the DC coupling capability, the BWT is dedicated to high integration required in long-term and low-power ECG recording systems.

Focal Epilepsy: MR Imaging of Nonhemodynamic Field Effects by Using a Phase-cycled Stimulus-induced Rotary Saturation Approach with Spin-Lock Preparation.

Purpose To investigate whether nonhemodynamic resonant saturation effects can be detected in patients with focal epilepsy by using a phase-cycled stimulus-induced rotary saturation (PC-SIRS) approach with spin-lock (SL) preparation and whether they colocalize with the seizure onset zone and surface interictal epileptiform discharges (IED). Materials and Methods The study was approved by the local ethics committee, and all subjects gave written informed consent. Eight patients with focal epilepsy undergoing presurgical surface and intracranial electroencephalography (EEG) underwent magnetic resonance (MR) imaging at 3 T with a whole-brain PC-SIRS imaging sequence with alternating SL-on and SL-off and two-dimensional echo-planar readout. The power of the SL radiofrequency pulse was set to 120 Hz to sensitize the sequence to high gamma oscillations present in epileptogenic tissue. Phase cycling was applied to capture distributed current orientations. Voxel-wise subtraction of SL-off from SL-on images enabled the separation of T2* effects from rotary saturation effects. The topography of PC-SIRS effects was compared with the seizure onset zone at intracranial EEG and with surface IED-related potentials. Bayesian statistics were used to test whether prior PC-SIRS information could improve IED source reconstruction. Results Nonhemodynamic resonant saturation effects ipsilateral to the seizure onset zone were detected in six of eight patients (concordance rate, 0.75; 95% confidence interval: 0.40, 0.94) by means of the PC-SIRS technique. They were concordant with IED surface negativity in seven of eight patients (0.88; 95% confidence interval: 0.51, 1.00). Including PC-SIRS as prior information improved the evidence of the standard EEG source models compared with the use of uninformed reconstructions (exceedance probability, 0.77 vs 0.12; Wilcoxon test of model evidence, P < .05). Nonhemodynamic resonant saturation effects resolved in patients with favorable postsurgical outcomes, but persisted in patients with postsurgical seizure recurrence. Conclusion Nonhemodynamic resonant saturation effects are detectable during interictal periods with the PC-SIRS approach in patients with epilepsy. The method may be useful for MR imaging-based detection of neuronal currents in a clinical environment. (©) RSNA, 2016 Online supplemental material is available for this article.

CHOLINERGIC MODULATION OF THALAMIC INPUT INTO THE CINGULATE CORTEX

This research demonstrates cholinergic modulation of thalamic input into the limbic cortex. A projection from the mediodorsal thalamus (MD) to the anterior cingulate cortex was defined anatomically and physiologically. Injections of horse-radish peroxidase into the anterior cingulate cortex labels neurons in the lateral, parvocellular, region of MD. Electrical Stimulation of this area produces a complex field potential in the anterior cingulate cortex which was further characterized by current density analysis and single cell recordings.^ The monsynaptic component of the response was identified as a large negative field which is maximal in layer IV of the anterior cingulate cortex. This response shows remarkable tetanic potentiation of frequencies near 7 Hz. During a train of 50 or more stimuli, the response would grow quickly and remain at a fairly stable potentiated level throughout the train.^ Cholinergic modulation of this thalamic response was demonstrated by iontophoretic application of the cholinergic agonist carbachol decreased the effectiveness of the thalamic imput by rapidly attenuation the response during a train of stimuli. The effect was apparently mediated by muscarinic receptors since the effect of carbachol was blocked by atropine but not by hexamethonium.^ To determine the source of the cingulate cortex cholinergic innervation, lesions were made in the anterior and medial thalamus and in the nucleus of the diagonal band of Broca. The effects of these lesions on choline acetyltranferase activity in the cingulate cortex were determined by a micro-radio-enzymatical assay. Only the lesions of the nucleus of the diagonal band significantly decreased the choline acetyltransferase activity in the cingulate cortex regions. Therefore, the diagonal band appears to be a major source of sensory cholinergic innervation and may be involved in gating of sensory information from the thalamus into the limbic cortex. Attempts to modulate the cingulate response to MD stimulation with electrical stimulation of the diagonal band, however were not successful.^

Measuring Scale Efficiency from the Translog Multi-Input, Multi-Output Distance Function

Ray (1998) developed measures of input- and output-oriented scale efficiency that can be directly computed from an estimated Translog frontier production function. This note extends the earlier results from Ray (1998) to the multiple-output multiple input case.

Comparing Input- and Output-Oriented Measures of Technical Efficiency to Determine Local Returns to Scale in DEA Models

This paper shows how one can infer the nature of local returns to scale at the input- or output-oriented efficient projection of a technically inefficient input-output bundle, when the input- and output-oriented measures of efficiency differ.

Input Aggregation in Models of Data Envelopment Analysis: A Statistical Test with an Application to Indian Manufacturing

A problem frequently encountered in Data Envelopment Analysis (DEA) is that the total number of inputs and outputs included tend to be too many relative to the sample size. One way to counter this problem is to combine several inputs (or outputs) into (meaningful) aggregate variables reducing thereby the dimension of the input (or output) vector. A direct effect of input aggregation is to reduce the number of constraints. This, in its turn, alters the optimal value of the objective function. In this paper, we show how a statistical test proposed by Banker (1993) may be applied to test the validity of a specific way of aggregating several inputs. An empirical application using data from Indian manufacturing for the year 2002-03 is included as an example of the proposed test.

Input Price Variation Across Locations and a Generalized Measure of Cost Efficiency

We propose a nonparametric model for global cost minimization as a framework for optimal allocation of a firm's output target across multiple locations, taking account of differences in input prices and technologies across locations. This should be useful for firms planning production sites within a country and for foreign direct investment decisions by multi-national firms. Two illustrative examples are included. The first example considers the production location decision of a manufacturing firm across a number of adjacent states of the US. In the other example, we consider the optimal allocation of US and Canadian automobile manufacturers across the two countries.

TIME SERIES ANALYSIS AS INPUT FOR PREDICTIVE MODELING: PREDICTING CARDIAC ARREST IN A PEDIATRIC INTENSIVE CARE UNIT

The first manuscript, entitled "Time-Series Analysis as Input for Clinical Predictive Modeling: Modeling Cardiac Arrest in a Pediatric ICU" lays out the theoretical background for the project. There are several core concepts presented in this paper. First, traditional multivariate models (where each variable is represented by only one value) provide single point-in-time snapshots of patient status: they are incapable of characterizing deterioration. Since deterioration is consistently identified as a precursor to cardiac arrests, we maintain that the traditional multivariate paradigm is insufficient for predicting arrests. We identify time series analysis as a method capable of characterizing deterioration in an objective, mathematical fashion, and describe how to build a general foundation for predictive modeling using time series analysis results as latent variables. Building a solid foundation for any given modeling task involves addressing a number of issues during the design phase. These include selecting the proper candidate features on which to base the model, and selecting the most appropriate tool to measure them. We also identified several unique design issues that are introduced when time series data elements are added to the set of candidate features. One such issue is in defining the duration and resolution of time series elements required to sufficiently characterize the time series phenomena being considered as candidate features for the predictive model. Once the duration and resolution are established, there must also be explicit mathematical or statistical operations that produce the time series analysis result to be used as a latent candidate feature. In synthesizing the comprehensive framework for building a predictive model based on time series data elements, we identified at least four classes of data that can be used in the model design. The first two classes are shared with traditional multivariate models: multivariate data and clinical latent features. Multivariate data is represented by the standard one value per variable paradigm and is widely employed in a host of clinical models and tools. These are often represented by a number present in a given cell of a table. Clinical latent features derived, rather than directly measured, data elements that more accurately represent a particular clinical phenomenon than any of the directly measured data elements in isolation. The second two classes are unique to the time series data elements. The first of these is the raw data elements. These are represented by multiple values per variable, and constitute the measured observations that are typically available to end users when they review time series data. These are often represented as dots on a graph. The final class of data results from performing time series analysis. This class of data represents the fundamental concept on which our hypothesis is based. The specific statistical or mathematical operations are up to the modeler to determine, but we generally recommend that a variety of analyses be performed in order to maximize the likelihood that a representation of the time series data elements is produced that is able to distinguish between two or more classes of outcomes. The second manuscript, entitled "Building Clinical Prediction Models Using Time Series Data: Modeling Cardiac Arrest in a Pediatric ICU" provides a detailed description, start to finish, of the methods required to prepare the data, build, and validate a predictive model that uses the time series data elements determined in the first paper. One of the fundamental tenets of the second paper is that manual implementations of time series based models are unfeasible due to the relatively large number of data elements and the complexity of preprocessing that must occur before data can be presented to the model. Each of the seventeen steps is analyzed from the perspective of how it may be automated, when necessary. We identify the general objectives and available strategies of each of the steps, and we present our rationale for choosing a specific strategy for each step in the case of predicting cardiac arrest in a pediatric intensive care unit. Another issue brought to light by the second paper is that the individual steps required to use time series data for predictive modeling are more numerous and more complex than those used for modeling with traditional multivariate data. Even after complexities attributable to the design phase (addressed in our first paper) have been accounted for, the management and manipulation of the time series elements (the preprocessing steps in particular) are issues that are not present in a traditional multivariate modeling paradigm. In our methods, we present the issues that arise from the time series data elements: defining a reference time; imputing and reducing time series data in order to conform to a predefined structure that was specified during the design phase; and normalizing variable families rather than individual variable instances. The final manuscript, entitled: "Using Time-Series Analysis to Predict Cardiac Arrest in a Pediatric Intensive Care Unit" presents the results that were obtained by applying the theoretical construct and its associated methods (detailed in the first two papers) to the case of cardiac arrest prediction in a pediatric intensive care unit. Our results showed that utilizing the trend analysis from the time series data elements reduced the number of classification errors by 73%. The area under the Receiver Operating Characteristic curve increased from a baseline of 87% to 98% by including the trend analysis. In addition to the performance measures, we were also able to demonstrate that adding raw time series data elements without their associated trend analyses improved classification accuracy as compared to the baseline multivariate model, but diminished classification accuracy as compared to when just the trend analysis features were added (ie, without adding the raw time series data elements). We believe this phenomenon was largely attributable to overfitting, which is known to increase as the ratio of candidate features to class examples rises. Furthermore, although we employed several feature reduction strategies to counteract the overfitting problem, they failed to improve the performance beyond that which was achieved by exclusion of the raw time series elements. Finally, our data demonstrated that pulse oximetry and systolic blood pressure readings tend to start diminishing about 10-20 minutes before an arrest, whereas heart rates tend to diminish rapidly less than 5 minutes before an arrest.