5-HT and GABA modulate intrinsic excitability of type I interneurons in Hermissenda.

The sensory neurons (photoreceptors) in the visual system of Hermissenda are one site of plasticity produced by Pavlovian conditioning. A second site of plasticity produced by conditioning is the type I interneurons in the cerebropleural ganglia. Both photoreceptors and statocyst hair cells of the graviceptive system form monosynaptic connections with identified type I interneurons. Two proposed neurotransmitters in the graviceptive system, serotonin (5-HT) and gamma-aminobutyric acid (GABA), have been shown to modify synaptic strength and intrinsic neuronal excitability in identified photoreceptors. However, the potential role of 5-HT and GABA in plasticity of type I interneurons has not been investigated. Here we show that 5-HT increased the peak amplitude of light-evoked complex excitatory postsynaptic potentials (EPSPs), enhanced intrinsic excitability, and increased spike activity of identified type I(e(A)) interneurons. In contrast, 5-HT decreased spike activity and intrinsic excitability of type I(e(B)) interneurons. The classification of two categories of type I(e) interneurons was also supported by the observation that 5-HT produced opposite effects on whole cell steady-state outward currents in type I(e) interneurons. Serotonin produced a reduction in the amplitude of light-evoked complex inhibitory PSPs (IPSPs), increased spontaneous spike activity, decreased intrinsic excitability, and depolarized the resting membrane potential of identified type I(i) interneurons. In contrast to the effects of 5-HT, GABA produced inhibition in both types of I(e) interneurons and type I(i) interneurons. These results show that 5-HT and GABA can modulate the intrinsic excitability of type I interneurons independent of the presynaptic effects of the same transmitters on excitability and synaptic efficacy of photoreceptors.

Spike timing dependent plasticity: a consequence of more fundamental learning rules.

Spike timing dependent plasticity (STDP) is a phenomenon in which the precise timing of spikes affects the sign and magnitude of changes in synaptic strength. STDP is often interpreted as the comprehensive learning rule for a synapse - the "first law" of synaptic plasticity. This interpretation is made explicit in theoretical models in which the total plasticity produced by complex spike patterns results from a superposition of the effects of all spike pairs. Although such models are appealing for their simplicity, they can fail dramatically. For example, the measured single-spike learning rule between hippocampal CA3 and CA1 pyramidal neurons does not predict the existence of long-term potentiation one of the best-known forms of synaptic plasticity. Layers of complexity have been added to the basic STDP model to repair predictive failures, but they have been outstripped by experimental data. We propose an alternate first law: neural activity triggers changes in key biochemical intermediates, which act as a more direct trigger of plasticity mechanisms. One particularly successful model uses intracellular calcium as the intermediate and can account for many observed properties of bidirectional plasticity. In this formulation, STDP is not itself the basis for explaining other forms of plasticity, but is instead a consequence of changes in the biochemical intermediate, calcium. Eventually a mechanism-based framework for learning rules should include other messengers, discrete change at individual synapses, spread of plasticity among neighboring synapses, and priming of hidden processes that change a synapse's susceptibility to future change. Mechanism-based models provide a rich framework for the computational representation of synaptic plasticity.

A biophysical model of synaptic plasticity and metaplasticity can account for the dynamics of the backward shift of hippocampal place fields.

Hippocampal place cells in the rat undergo experience-dependent changes when the rat runs stereotyped routes. One such change, the backward shift of the place field center of mass, has been linked by previous modeling efforts to spike-timing-dependent plasticity (STDP). However, these models did not account for the termination of the place field shift and they were based on an abstract implementation of STDP that ignores many of the features found in cortical plasticity. Here, instead of the abstract STDP model, we use a calcium-dependent plasticity (CaDP) learning rule that can account for many of the observed properties of cortical plasticity. We use the CaDP learning rule in combination with a model of metaplasticity to simulate place field dynamics. Without any major changes to the parameters of the original model, the present simulations account both for the initial rapid place field shift and for the subsequent slowing down of this shift. These results suggest that the CaDP model captures the essence of a general cortical mechanism of synaptic plasticity, which may underlie numerous forms of synaptic plasticity observed both in vivo and in vitro.

Convergent and divergent modulatory pathways in {\it Aplysia\/}: Cellular and network studies

Neuromodulation is essential to many functions of the nervous system. In the simple gastropod mollusk Aplysia californica, neuromodulation of the circuits for the defensive withdrawal reflexes has been associated with several forms of learning. In the present work, the neurotransmitters and neural circuitry which contribute to the modulation of the tail-siphon withdrawal reflex were examined.^ A recently-identified neuropeptide transmitter, buccalin A was found to modulate the biophysical properties of the sensory neurons that mediate the reflex. The actions of buccalin A on the sensory neurons were compared with those of the well-characterized modulatory transmitter serotonin, and convergence and divergence in the actions of these two transmitters were evaluated. Buccalin A dramatically increased the excitability of sensory neurons and occluded further enhancement of excitability by serotonin. Buccalin A produced no significant change in spike duration, and it did not block serotonin-induced spike broadening. Voltage-clamp analysis revealed the currents that may be involved in the effects on spike duration and excitability. Buccalin A decreased an outward current similar to the S-K$\sp+$ current (I$\sb{\rm K,S}$). Buccalin A appeared to occlude further modulation of I$\sb{\rm K,S}$ by serotonin, but did not block serotonin-induced modulation of the voltage-dependent delayed rectifier K$\sp+$ current (I$\sb{\rm K,V}$). These results suggest that buccalin A converges on some, but not all, of the same subcellular modulatory pathways as serotonin.^ In order to begin to understand neuromodulation in a more physiological context for the tail-siphon withdrawal reflex, the modulatory circuitry for the tail-withdrawal circuit was examined. Mechanoafferent neurons in the J cluster of the cerebral ganglion were identified as elements of a modulatory circuit for the reflex. Excitatory and inhibitory connections were observed between the J cells and the pleural sensory neurons, the tail motor neurons, and several classes of interneurons for the tail-siphon withdrawal circuit. The J cells produced both fast and slow PSPs in these neurons. Of particular interest was the ability of the J cells to produce slow EPSPs in the pleural sensory neurons. These slow EPSPs were associated with an increase in the excitability of the sensory neurons. The J cells appear to mediate both sensory and modulatory inputs to the circuit for the tail-siphon withdrawal reflex from the anterior part of the animal. ^

THE SITE AND MECHANISM OF ACTION OF BENZODIAZEPINE IN RABBIT HIPPOCAMPUS

The purpose of this study was to determine the effects of benzodiazepine in area CA1 of the hippocampus, and to explore possible mechanisms of action for these agents in this brain area. Two distinctly different benzodiazepine-induced changes in hippocampal physiology have been identified. First, benzodiazepine depresses the population spike recorded in stratum pyramidale, indicating a decrease in action potential generation. Second, benzodiazepine decreases the magnitude of post-tetanic potentiation of the population EPSP recorded in stratum radiatum, and shortens the duration. The effect of benzodiazepine on pyramidal cell excitation was reversed by the GABA antagonis bicuculline, and mimicked by GABA itself. Thus the available evidence is consistent with the hypothesis that benzodiazepine acts by enhancing the effect of GABA in this area. In stratum radiatum, on the other hand, the effect of benzodiazepine on post tetanic potentiation of the population EPSP was not altered by bicuculline although bicuculline did antagonize GABA in this area. In addition, application of GABA, while it caused profound changes in the population EPSP,p, did not cause the same changes that were induced by benzodiazepine. Thus the evidence does not support the hypothesis that benzodiazepine is acting in stratum radiatum by enhancing the effects of GABA. ^

Bayesian joint modeling of longitudinal and survival data

The joint modeling of longitudinal and survival data is a new approach to many applications such as HIV, cancer vaccine trials and quality of life studies. There are recent developments of the methodologies with respect to each of the components of the joint model as well as statistical processes that link them together. Among these, second order polynomial random effect models and linear mixed effects models are the most commonly used for the longitudinal trajectory function. In this study, we first relax the parametric constraints for polynomial random effect models by using Dirichlet process priors, then three longitudinal markers rather than only one marker are considered in one joint model. Second, we use a linear mixed effect model for the longitudinal process in a joint model analyzing the three markers. In this research these methods were applied to the Primary Biliary Cirrhosis sequential data, which were collected from a clinical trial of primary biliary cirrhosis (PBC) of the liver. This trial was conducted between 1974 and 1984 at the Mayo Clinic. The effects of three longitudinal markers (1) Total Serum Bilirubin, (2) Serum Albumin and (3) Serum Glutamic-Oxaloacetic transaminase (SGOT) on patients' survival were investigated. Proportion of treatment effect will also be studied using the proposed joint modeling approaches. ^ Based on the results, we conclude that the proposed modeling approaches yield better fit to the data and give less biased parameter estimates for these trajectory functions than previous methods. Model fit is also improved after considering three longitudinal markers instead of one marker only. The results from analysis of proportion of treatment effects from these joint models indicate same conclusion as that from the final model of Fleming and Harrington (1991), which is Bilirubin and Albumin together has stronger impact in predicting patients' survival and as a surrogate endpoints for treatment. ^

RNA-SEQUENCING APPLICATIONS: GENE EXPRESSION QUANTIFICATION AND METHYLATOR PHENOTYPE IDENTIFICATION

My dissertation focuses on two aspects of RNA sequencing technology. The first is the methodology for modeling the overdispersion inherent in RNA-seq data for differential expression analysis. This aspect is addressed in three sections. The second aspect is the application of RNA-seq data to identify the CpG island methylator phenotype (CIMP) by integrating datasets of mRNA expression level and DNA methylation status. Section 1: The cost of DNA sequencing has reduced dramatically in the past decade. Consequently, genomic research increasingly depends on sequencing technology. However it remains elusive how the sequencing capacity influences the accuracy of mRNA expression measurement. We observe that accuracy improves along with the increasing sequencing depth. To model the overdispersion, we use the beta-binomial distribution with a new parameter indicating the dependency between overdispersion and sequencing depth. Our modified beta-binomial model performs better than the binomial or the pure beta-binomial model with a lower false discovery rate. Section 2: Although a number of methods have been proposed in order to accurately analyze differential RNA expression on the gene level, modeling on the base pair level is required. Here, we find that the overdispersion rate decreases as the sequencing depth increases on the base pair level. Also, we propose four models and compare them with each other. As expected, our beta binomial model with a dynamic overdispersion rate is shown to be superior. Section 3: We investigate biases in RNA-seq by exploring the measurement of the external control, spike-in RNA. This study is based on two datasets with spike-in controls obtained from a recent study. We observe an undiscovered bias in the measurement of the spike-in transcripts that arises from the influence of the sample transcripts in RNA-seq. Also, we find that this influence is related to the local sequence of the random hexamer that is used in priming. We suggest a model of the inequality between samples and to correct this type of bias. Section 4: The expression of a gene can be turned off when its promoter is highly methylated. Several studies have reported that a clear threshold effect exists in gene silencing that is mediated by DNA methylation. It is reasonable to assume the thresholds are specific for each gene. It is also intriguing to investigate genes that are largely controlled by DNA methylation. These genes are called “L-shaped” genes. We develop a method to determine the DNA methylation threshold and identify a new CIMP of BRCA. In conclusion, we provide a detailed understanding of the relationship between the overdispersion rate and sequencing depth. And we reveal a new bias in RNA-seq and provide a detailed understanding of the relationship between this new bias and the local sequence. Also we develop a powerful method to dichotomize methylation status and consequently we identify a new CIMP of breast cancer with a distinct classification of molecular characteristics and clinical features.

Neocortical hyperexcitability defect in a mutant mouse model of spike-wave epilepsy, {\it stargazer}

Single-locus mutations in mice can express epileptic phenotypes and provide critical insights into the naturally occurring defects that alter excitability and mediate synchronization in the central nervous system (CNS). One such recessive mutation (on chromosome (Chr) 15), stargazer(stg/stg) expresses frequent bilateral 6-7 cycles per second (c/sec) spike-wave seizures associated with behavioral arrest, and provides a valuable opportunity to examine the inherited lesion associated with spike-wave synchronization.^ The existence of distinct and heterogeneous defects mediating spike-wave discharge (SWD) generation has been demonstrated by the presence of multiple genetic loci expressing generalized spike-wave activity and the differential effects of pharmacological agents on SWDs in different spike-wave epilepsy models. Attempts at understanding the different basic mechanisms underlying spike-wave synchronization have focused on $\gamma$-aminobutyric acid (GABA) receptor-, low threshold T-type Ca$\sp{2+}$ channel-, and N-methyl-D-aspartate receptor (NMDA-R)-mediated transmission. It is believed that defects in these modes of transmission can mediate the conversion of normal oscillations in a trisynaptic circuit, which includes the neocortex, reticular nucleus and thalamus, into spike-wave activity. However, the underlying lesions involved in spike-wave synchronization have not been clearly identified.^ The purpose of this research project was to locate and characterize a distinct neuronal hyperexcitability defect favoring spike-wave synchronization in the stargazer brain. One experimental approach for anatomically locating areas of synchronization and hyperexcitability involved an attempt to map patterns of hypersynchronous activity with antibodies to activity-induced proteins.^ A second approach to characterizing the neuronal defect involved examining the neuronal responses in the mutant following application of pharmacological agents with well known sites of action.^ In order to test the hypothesis that an NMDA receptor mediated hyperexcitability defect exists in stargazer neocortex, extracellular field recordings were used to examine the effects of CPP and MK-801 on coronal neocortical brain slices of stargazer and wild type perfused with 0 Mg$\sp{2+}$ artificial cerebral spinal fluid (aCSF).^ To study how NMDA receptor antagonists might promote increased excitability in stargazer neocortex, two basic hypotheses were tested: (1) NMDA receptor antagonists directly activate deep layer principal pyramidal cells in the neocortex of stargazer, presumably by opening NMDA receptor channels altered by the stg mutation; and (2) NMDA receptor antagonists disinhibit the neocortical network by blocking recurrent excitatory synaptic inputs onto inhibitory interneurons in the deep layers of stargazer neocortex.^ In order to test whether CPP might disinhibit the 0 Mg$\sp{2+}$ bursting network in the mutant by acting on inhibitory interneurons, the inhibitory inputs were pharmacologically removed by application of GABA receptor antagonists to the cortical network, and the effects of CPP under 0 Mg$\sp{2+}$aCSF perfusion in layer V of stg/stg were then compared with those found in +/+ neocortex using in vitro extracellular field recordings. (Abstract shortened by UMI.) ^