Detection and localization of multiple rate changes in Poisson spike trains

Poster presentation from Twentieth Annual Computational Neuroscience Meeting: CNS*2011 Stockholm, Sweden. 23-28 July 2011. In statistical spike train analysis, stochastic point process models usually assume stationarity, in particular that the underlying spike train shows a constant firing rate (e.g. [1]). However, such models can lead to misinterpretation of the associated tests if the assumption of rate stationarity is not met (e.g. [2]). Therefore, the analysis of nonstationary data requires that rate changes can be located as precisely as possible. However, present statistical methods focus on rejecting the null hypothesis of stationarity without explicitly locating the change point(s) (e.g. [3]). We propose a test for stationarity of a given spike train that can also be used to estimate the change points in the firing rate. Assuming a Poisson process with piecewise constant firing rate, we propose a Step-Filter-Test (SFT) which can work simultaneously in different time scales, accounting for the high variety of firing patterns in experimental spike trains. Formally, we compare the numbers N1=N1(t,h) and N2=N2(t,h) of spikes in the time intervals (t-h,t] and (h,t+h]. By varying t within a fine time lattice and simultaneously varying the interval length h, we obtain a multivariate statistic D(h,t):=(N1-N2)/V(N1+N2), for which we prove asymptotic multivariate normality under homogeneity. From this a practical, graphical device to spot changes of the firing rate is constructed. Our graphical representation of D(h,t) (Figure 1A) visualizes the changes in the firing rate. For the statistical test, a threshold K is chosen such that under homogeneity, |D(h,t)|<K holds for all investigated h and t with probability 0.95. This threshold can indicate potential change points in order to estimate the inhomogeneous rate profile (Figure 1B). The SFT is applied to a sample data set of spontaneous single unit activity recorded from the substantia nigra of anesthetized mice. In this data set, multiple rate changes are identified which agree closely with visual inspection. In contrast to approaches choosing one fixed kernel width [4], our method has advantages in the flexibility of h

Statistical Properties of Charmonium Spectrum and a New Mechanism of J/\psi Suppression

The statistical properties of Charmonium energy spectrum determined by the Bethe-Salpeter equation are investigated. It is found that the regular motion of the $c\bar{c}$ system can be expected at a small value of color screening mass but the chaotic motion at a large one. It is shown that the level mixing due to color screening serves as a new mechanism resulting in the J/$\psi$ suppression. Moreover, this kind of suppression can occur before the color screening mass reaches its critical value for J/$\psi$ dissociation. It implies that a strong J/$\psi$ suppression is possible in the absence of dissociation of J/$\psi$.Comment: 13 latex pages, 2 figures. Phys. Rev. C in pres

Aptamer-based multiplexed proteomic technology for biomarker discovery

Interrogation of the human proteome in a highly multiplexed and efficient manner remains a coveted and challenging goal in biology. We present a new aptamer-based proteomic technology for biomarker discovery capable of simultaneously measuring thousands of proteins from small sample volumes (15 [mu]L of serum or plasma). Our current assay allows us to measure ~800 proteins with very low limits of detection (1 pM average), 7 logs of overall dynamic range, and 5% average coefficient of variation. This technology is enabled by a new generation of aptamers that contain chemically modified nucleotides, which greatly expand the physicochemical diversity of the large randomized nucleic acid libraries from which the aptamers are selected. Proteins in complex matrices such as plasma are measured with a process that transforms a signature of protein concentrations into a corresponding DNA aptamer concentration signature, which is then quantified with a DNA microarray. In essence, our assay takes advantage of the dual nature of aptamers as both folded binding entities with defined shapes and unique sequences recognizable by specific hybridization probes. To demonstrate the utility of our proteomics biomarker discovery technology, we applied it to a clinical study of chronic kidney disease (CKD). We identified two well known CKD biomarkers as well as an additional 58 potential CKD biomarkers. These results demonstrate the potential utility of our technology to discover unique protein signatures characteristic of various disease states. More generally, we describe a versatile and powerful tool that allows large-scale comparison of proteome profiles among discrete populations. This unbiased and highly multiplexed search engine will enable the discovery of novel biomarkers in a manner that is unencumbered by our incomplete knowledge of biology, thereby helping to advance the next generation of evidence-based medicine

Validation of a blood protein signature for non-small cell lung cancer

Background: CT screening for lung cancer is effective in reducing mortality, but there are areas of concern, including a positive predictive value of 4% and development of interval cancers. A blood test that could manage these limitations would be useful, but development of such tests has been impaired by variations in blood collection that may lead to poor reproducibility across populations. Results: Blood-based proteomic profiles were generated with SOMAscan technology, which measured 1033 proteins. First, preanalytic variability was evaluated with Sample Mapping Vectors (SMV), which are panels of proteins that detect confounders in protein levels related to sample collection. A subset of well collected serum samples not influenced by preanalytic variability was selected for discovery of lung cancer biomarkers. The impact of sample collection variation on these candidate markers was tested in the subset of samples with higher SMV scores so that the most robust markers could be used to create disease classifiers. The discovery sample set (n = 363) was from a multi-center study of 94 non-small cell lung cancer (NSCLC) cases and 269 long-term smokers and benign pulmonary nodule controls. The analysis resulted in a 7-marker panel with an AUC of 0.85 for all cases (68% adenocarcinoma, 32% squamous) and an AUC of 0.93 for squamous cell carcinoma in particular. This panel was validated by making blinded predictions in two independent cohorts (n = 138 in the first validation and n = 135 in the second). The model was recalibrated for a panel format prior to unblinding the second cohort. The AUCs overall were 0.81 and 0.77, and for squamous cell tumors alone were 0.89 and 0.87. The estimated negative predictive value for a 15% disease prevalence was 93% overall and 99% for squamous lung tumors. The proteins in the classifier function in destruction of the extracellular matrix, metabolic homeostasis and inflammation. Conclusions: Selecting biomarkers resistant to sample processing variation led to robust lung cancer biomarkers that performed consistently in independent validations. They form a sensitive signature for detection of lung cancer, especially squamous cell histology. This non-invasive test could be used to improve the positive predictive value of CT screening, with the potential to avoid invasive evaluation of nonmalignant pulmonary nodules

‘There is a Time to be Born and a Time to Die’ (Ecclesiastes 3:2a): Jewish Perspectives on Euthanasia

Reviewing the publications of prominent American rabbis who have (extensively) published on Jewish biomedical ethics, this article highlights Orthodox, Conservative and Reform opinions on a most pressing contemporary bioethical issue: euthanasia. Reviewing their opinions against the background of the halachic character of Jewish (biomedical) ethics, this article shows how from one traditional Jewish textual source diverse, even contradictory, opinions emerge through different interpretations. In this way, in the Jewish debate on euthanasia the specific methodology of Jewish (bio)ethical reasoning comes forward as well as a diversity of opinion within Judaism and its branches

A Naturally Occurring Polymorphism at Drosophila melanogaster Lim3 Locus, a Homolog of Human LHX3/4, Affects Lim3 Transcription and Fly Lifespan

Lim3 encodes an RNA polymerase II transcription factor with a key role in neuron specification. It was also identified as a candidate gene that affects lifespan. These pleiotropic effects indicate the fundamental significance of the potential interplay between neural development and lifespan control. The goal of this study was to analyze the causal relationships between Lim3 structural variations, and gene expression and lifespan changes, and to provide insights into regulatory pathways controlling lifespan. Fifty substitution lines containing second chromosomes from a Drosophila natural population were used to analyze the association between lifespan and sequence variation in the 5′-regulatory region, and first exon and intron of Lim3A, in which we discovered multiple transcription start sites (TSS). The core and proximal promoter organization for Lim3A and a previously unknown mRNA named Lim3C were described. A haplotype of two markers in the Lim3A regulatory region was significantly associated with variation in lifespan. We propose that polymorphisms in the regulatory region affect gene transcription, and consequently lifespan. Indeed, five polymorphic markers located within 380 to 680 bp of the Lim3A major TSS, including two markers associated with lifespan variation, were significantly associated with the level of Lim3A transcript, as evaluated by real time RT-PCR in embryos, adult heads, and testes. A naturally occurring polymorphism caused a six-fold change in gene transcription and a 25% change in lifespan. Markers associated with long lifespan and intermediate Lim3A transcription were present in the population at high frequencies. We hypothesize that polymorphic markers associated with Lim3A expression are located within the binding sites for proteins that regulate gene function, and provide general rather than tissue-specific regulation of transcription, and that intermediate levels of Lim3A expression confer a selective advantage and longer lifespan

State-Space Analysis of Time-Varying Higher-Order Spike Correlation for Multiple Neural Spike Train Data

Precise spike coordination between the spiking activities of multiple neurons is suggested as an indication of coordinated network activity in active cell assemblies. Spike correlation analysis aims to identify such cooperative network activity by detecting excess spike synchrony in simultaneously recorded multiple neural spike sequences. Cooperative activity is expected to organize dynamically during behavior and cognition; therefore currently available analysis techniques must be extended to enable the estimation of multiple time-varying spike interactions between neurons simultaneously. In particular, new methods must take advantage of the simultaneous observations of multiple neurons by addressing their higher-order dependencies, which cannot be revealed by pairwise analyses alone. In this paper, we develop a method for estimating time-varying spike interactions by means of a state-space analysis. Discretized parallel spike sequences are modeled as multi-variate binary processes using a log-linear model that provides a well-defined measure of higher-order spike correlation in an information geometry framework. We construct a recursive Bayesian filter/smoother for the extraction of spike interaction parameters. This method can simultaneously estimate the dynamic pairwise spike interactions of multiple single neurons, thereby extending the Ising/spin-glass model analysis of multiple neural spike train data to a nonstationary analysis. Furthermore, the method can estimate dynamic higher-order spike interactions. To validate the inclusion of the higher-order terms in the model, we construct an approximation method to assess the goodness-of-fit to spike data. In addition, we formulate a test method for the presence of higher-order spike correlation even in nonstationary spike data, e.g., data from awake behaving animals. The utility of the proposed methods is tested using simulated spike data with known underlying correlation dynamics. Finally, we apply the methods to neural spike data simultaneously recorded from the motor cortex of an awake monkey and demonstrate that the higher-order spike correlation organizes dynamically in relation to a behavioral demand

Relationship Contexts as Sources of Socialization: An Exploration of Intimate Partner Violence Experiences of Economically Disadvantaged African American Adolescents

Intimate partner violence (IPV) among African Americans is a serious public health concern. Research suggest that African Americans adolescents, particularly those from economically disadvantaged communities, are at heightened risk for experiencing and perpetrating dating violence compared to youth from other racial and ethnic groups. In the present study, we examined different relationship contexts that are sources of IPV socialization. Semi-structured interviews were conducted with 22 economically disadvantaged African American adolescents. Content analysis yielded five relationship contexts through which the participants witnessed, experienced, and perpetrated IPV: (a) adolescents’ own dating relationships (64%), (b) siblings and extended family members (e.g., cousins, aunts, uncles) (59%), (c) parent-partners (27%), (d) friends (23%), and (e) neighbors (18%). Adolescents also frequently described IPV in their own dating relationships and in parent-partner relationships as mutual. Moreover, they appeared to minimize the experience of IPV in their own relationships. Efforts to reduce rates of IPV among economically disadvantaged African American adolescents should consider these relational contexts through which adolescents are socialized with regards to IPV and adolescents’ beliefs about mutual violence in relationships. Results highlight the importance of culturally relevant prevention and intervention programs that consider these relationship contexts

HEATR2 Plays a Conserved Role in Assembly of the Ciliary Motile Apparatus

Cilia are highly conserved microtubule-based structures that perform a variety of sensory and motility functions during development and adult homeostasis. In humans, defects specifically affecting motile cilia lead to chronic airway infections, infertility and laterality defects in the genetically heterogeneous disorder Primary Ciliary Dyskinesia (PCD). Using the comparatively simple Drosophila system, in which mechanosensory neurons possess modified motile cilia, we employed a recently elucidated cilia transcriptional RFX-FOX code to identify novel PCD candidate genes. Here, we report characterization of CG31320/HEATR2, which plays a conserved critical role in forming the axonemal dynein arms required for ciliary motility in both flies and humans. Inner and outer arm dyneins are absent from axonemes of CG31320 mutant flies and from PCD individuals with a novel splice-acceptor HEATR2 mutation. Functional conservation of closely arranged RFX-FOX binding sites upstream of HEATR2 orthologues may drive higher cytoplasmic expression of HEATR2 during early motile ciliogenesis. Immunoprecipitation reveals HEATR2 interacts with DNAI2, but not HSP70 or HSP90, distinguishing it from the client/chaperone functions described for other cytoplasmic proteins required for dynein arm assembly such as DNAAF1-4. These data implicate CG31320/HEATR2 in a growing intracellular pre-assembly and transport network that is necessary to deliver functional dynein machinery to the ciliary compartment for integration into the motile axoneme