The physiological variability of channel density in hippocampal CA1 pyramidal cells and interneurons explored using a unified data-driven modeling workflow.

Every neuron is part of a network, exerting its function by transforming multiple spatiotemporal synaptic input patterns into a single spiking output. This function is specified by the particular shape and passive electrical properties of the neuronal membrane, and the composition and spatial distribution of ion channels across its processes. For a variety of physiological or pathological reasons, the intrinsic input/output function may change during a neuron's lifetime. This process results in high variability in the peak specific conductance of ion channels in individual neurons. The mechanisms responsible for this variability are not well understood, although there are clear indications from experiments and modeling that degeneracy and correlation among multiple channels may be involved. Here, we studied this issue in biophysical models of hippocampal CA1 pyramidal neurons and interneurons. Using a unified data-driven simulation workflow and starting from a set of experimental recordings and morphological reconstructions obtained from rats, we built and analyzed several ensembles of morphologically and biophysically accurate single cell models with intrinsic electrophysiological properties consistent with experimental findings. The results suggest that the set of conductances expressed in any given hippocampal neuron may be considered as belonging to two groups: one subset is responsible for the major characteristics of the firing behavior in each population and the other is responsible for a robust degeneracy. Analysis of the model neurons suggests several experimentally testable predictions related to the combination and relative proportion of the different conductances that should be expressed on the membrane of different types of neurons for them to fulfill their role in the hippocampus circuitry.

Enhanced Sensitivity to Rapid Input Fluctuations by Nonlinear Threshold Dynamics in Neocortical Pyramidal Neurons.

The way in which single neurons transform input into output spike trains has fundamental consequences for network coding. Theories and modeling studies based on standard Integrate-and-Fire models implicitly assume that, in response to increasingly strong inputs, neurons modify their coding strategy by progressively reducing their selective sensitivity to rapid input fluctuations. Combining mathematical modeling with in vitro experiments, we demonstrate that, in L5 pyramidal neurons, the firing threshold dynamics adaptively adjust the effective timescale of somatic integration in order to preserve sensitivity to rapid signals over a broad range of input statistics. For that, a new Generalized Integrate-and-Fire model featuring nonlinear firing threshold dynamics and conductance-based adaptation is introduced that outperforms state-of-the-art neuron models in predicting the spiking activity of neurons responding to a variety of in vivo-like fluctuating currents. Our model allows for efficient parameter extraction and can be analytically mapped to a Generalized Linear Model in which both the input filter--describing somatic integration--and the spike-history filter--accounting for spike-frequency adaptation--dynamically adapt to the input statistics, as experimentally observed. Overall, our results provide new insights on the computational role of different biophysical processes known to underlie adaptive coding in single neurons and support previous theoretical findings indicating that the nonlinear dynamics of the firing threshold due to Na+-channel inactivation regulate the sensitivity to rapid input fluctuations.

Automated High-Throughput Characterization of Single Neurons by Means of Simplified Spiking Models.

Single-neuron models are useful not only for studying the emergent properties of neural circuits in large-scale simulations, but also for extracting and summarizing in a principled way the information contained in electrophysiological recordings. Here we demonstrate that, using a convex optimization procedure we previously introduced, a Generalized Integrate-and-Fire model can be accurately fitted with a limited amount of data. The model is capable of predicting both the spiking activity and the subthreshold dynamics of different cell types, and can be used for online characterization of neuronal properties. A protocol is proposed that, combined with emergent technologies for automatic patch-clamp recordings, permits automated, in vitro high-throughput characterization of single neurons.

Mutations in the polyglutamine binding protein 1 gene cause X-linked mental retardation.

We found mutations in the gene PQBP1 in 5 of 29 families with nonsyndromic (MRX) and syndromic (MRXS) forms of X-linked mental retardation (XLMR). Clinical features in affected males include mental retardation, microcephaly, short stature, spastic paraplegia and midline defects. PQBP1 has previously been implicated in the pathogenesis of polyglutamine expansion diseases. Our findings link this gene to XLMR and shed more light on the pathogenesis of this common disorder.

DEVILS: a tool for the visualization of large datasets with a high dynamic range.

The number of grey values that can be displayed on monitors and be processed by the human eye is smaller than the dynamic range of image-based sensors. This makes the visualization of such data a challenge, especially with specimens where small dim structures are equally important as large bright ones, or whenever variations in intensity, such as non-homogeneous staining efficiencies or light depth penetration, becomes an issue. While simple intensity display mappings are easily possible, these fail to provide a one-shot observation that can display objects of varying intensities. In order to facilitate the visualization-based analysis of large volumetric datasets, we developed an easy-to-use ImageJ plugin enabling the compressed display of features within several magnitudes of intensities. The Display Enhancement for Visual Inspection of Large Stacks plugin (DEVILS) homogenizes the intensities by using a combination of local and global pixel operations to allow for high and low intensities to be visible simultaneously to the human eye. The plugin is based on a single, intuitively understandable parameter, features a preview mode, and uses parallelization to process multiple image planes. As output, the plugin is capable of producing a BigDataViewer-compatible dataset for fast visualization. We demonstrate the utility of the plugin for large volumetric image data.

Inferring and validating mechanistic models of neural microcircuits based on spike-train data.

The interpretation of neuronal spike train recordings often relies on abstract statistical models that allow for principled parameter estimation and model selection but provide only limited insights into underlying microcircuits. In contrast, mechanistic models are useful to interpret microcircuit dynamics, but are rarely quantitatively matched to experimental data due to methodological challenges. Here we present analytical methods to efficiently fit spiking circuit models to single-trial spike trains. Using derived likelihood functions, we statistically infer the mean and variance of hidden inputs, neuronal adaptation properties and connectivity for coupled integrate-and-fire neurons. Comprehensive evaluations on synthetic data, validations using ground truth in-vitro and in-vivo recordings, and comparisons with existing techniques demonstrate that parameter estimation is very accurate and efficient, even for highly subsampled networks. Our methods bridge statistical, data-driven and theoretical, model-based neurosciences at the level of spiking circuits, for the purpose of a quantitative, mechanistic interpretation of recorded neuronal population activity.

Characterization of FBX25, encoding a novel brain-expressed F-box protein.

F-box proteins (FBPs) confer substrate specificity to the SCF-type (Skp1/Cul1/FBP) of ubiquitin ligase complexes through their F-box. Multiple FBPs have been predicted, but experimental evidence is lagging. We report on the predicted human FBP hFBX25 which we found to be disrupted in a mentally retarded translocation carrier suffering from epileptic seizures. We investigated hFBX25's genomic organization and established hFBX25 as an FBP by verifying its interaction with Skp1 and Cul1. In the process, we identified an atypical serine residue in the F-box which is crucial for the hFBX25-Skp1 binding. We determined hFBX25's subcellular localization. We found strong transcription in human brain. In mouse embryonic sections, mFbx25 shows predominantly neuronal expression and in adult mouse brain, expression is confined to the hippocampus, the cerebral cortex and the Purkinje cell layer. Interestingly, aberrations in the ubiquitin pathway have been linked to neurological conditions.

A new standard nomenclature for proteins related to Apx and Shroom.

Shroom is a recently-described regulator of cell shape changes in the developing nervous system. This protein is a member of a small family of related proteins that are defined by sequence similarity and in most cases by some link to the actin cytoskeleton. At present these proteins are named Shroom, APX, APXL, and KIAA1202. In light of the growing interest in this family of proteins, we propose here a new standard nomenclature.

Exonic microdeletions in the X-linked PQBP1 gene in mentally retarded patients: a pathogenic mutation and in-frame deletions of uncertain effect.

Mutations in PQBP1 were recently identified in families with syndromic and non-syndromic X-linked mental retardation (XLMR). Clinical features frequently associated with MR were microcephaly and/or short stature. The predominant mutations detected so far affect a stretch of six AG dinucleotides in the polar-amino-acid-rich domain (PRD), causing frameshifts in the fourth coding exon. We searched for PQBP1 exon 4 frameshifts in 57 mentally retarded males in whom initial referral description indicated at least one of the following criteria: microcephaly, short stature, spastic paraplegia or family history compatible with XLMR, and in 772 mentally retarded males not selected for specific clinical features or family history. We identified a novel frameshift mutation (23 bp deletion) in two half-brothers with specific clinical features, and performed prenatal diagnosis in this family. We also found two different 21 bp in-frame deletions (c.334-354del(21 bp) and c.393-413del(21 bp)) in four unrelated probands from various ethnic origins, each deleting one of five copies of an imperfect seven amino-acid repeat. Although such deletions have not been detected in 1180 X chromosomes from European controls, the c. 334-354del(21 bp) was subsequently found in two of 477 Xs from Indian controls. We conclude that pathogenic frameshift mutations in PQBP1 are rare in mentally retarded patients lacking specific associated signs and that the 21 bp in-frame deletions may be non-pathogenic, or alternatively could act subtly on PQBP1 function. This touches upon a common dilemma in XLMR, that is, how to distinguish between mutations and variants that may be non-pathogenic or represent risk factors for cognitive impairment.

Structural model of the OPA1 GTPase domain may explain the molecular consequences of a novel mutation in a family with autosomal dominant optic atrophy.

Autosomal dominant optic atrophy (ADOA) is the most frequent hereditary optic neuropathy. Three loci have been reported for ADOA: a major locus, harboring all identified mutations to date, maps to 3q28 (OPA1), a second locus is linked to 18q12.2-q12.3 (OPA4) and a third locus on 22q12.1-q13.1 (OPA5) has been reported recently. We describe a six-generation Iranian family in which optic atrophy runs as an autosomal dominant trait with an age of onset at 14-15years. We performed linkage analysis with markers mapping to 3q28 and 18q12.2-q12.3 and found linkage to 3q28. Subsequent sequencing of OPA1 identified a novel heterozygous missense mutation (c.1313A>G) replacing aspartic acid by glycine (p.D438G) in the GTPase domain of OPA1. Interestingly, another missense mutation at the same position (c.1313A>T, D438V) has been reported before in two unrelated German families, indicating a possible mutation hot spot. Further evidence supporting the importance of D438 is its conservation from human to acoelomata. OPA1 is believed to be the human orthologue of yeast MGM1, a dynamin-related protein required for the integrity of mitochondrial DNA. Homology modeling of the OPA1 GTPase domain revealed extensive structural similarity to the Dictyostelium dynamin A GTPase domain and showed that D438 may interact with residues of the G1 and the G4 motifs, which are crucial in coordinating GTP. Based on this analysis, we propose a mechanism which explains the gradual decline of vision in ADOA patients with OPA1 mutations at position 438.

Truncation of the CNS-expressed JNK3 in a patient with a severe developmental epileptic encephalopathy.

We have investigated the breakpoints in a male child with pharmacoresistant epileptic encephalopathy and a de novo balanced translocation t(Y;4)(q11.2;q21). By fluorescence in situ hybridisation, we have identified genomic clones from both chromosome 4 and chromosome Y that span the breakpoints. Precise mapping of the chromosome 4 breakpoint indicated that the c-Jun N-terminal kinase 3 (JNK3) gene is disrupted in the patient. This gene is predominantly expressed in the central nervous system, and it plays an established role in both neuronal differentiation and apoptosis. Expression studies in the patient lymphoblastoid cell line show that the truncated JNK3 protein is expressed, i.e. the disrupted transcript is not immediately subject to nonsense-mediated mRNA decay, as is often the case for truncated mRNAs or those harbouring premature termination codons. Over-expression studies with the mutant protein in various cell lines, including neural cells, indicate that both its solubility and cellular localisation differ from that of the wild-type JNK3. It is plausible, therefore, that the presence of the truncated JNK3 disrupts normal JNK3 signal transduction in neuronal cells. JNK3 is one of the downstream effectors of the GTPase-regulated MAP kinase cascade, several members of which have been implicated in cognitive function. In addition, two known JNK3-interacting proteins, beta-arrestin 2 and JIP3, play established roles in neurite outgrowth and neurological development. These interactions are likely affected by a truncated JNK3 protein, and thereby provide an explanation for the link between alterations in MAP kinase signal transduction and brain disorders.

Disruptions of the novel KIAA1202 gene are associated with X-linked mental retardation.

The extensive heterogeneity underlying the genetic component of mental retardation (MR) is the main cause for our limited understanding of the aetiology of this highly prevalent condition. Hence we set out to identify genes involved in MR. We investigated the breakpoints of two balanced X;autosome translocations in two unrelated female patients with mild/moderate MR and found that the Xp11.2 breakpoints disrupt the novel human KIAA1202 (hKIAA1202) gene in both cases. We also identified a missense exchange in this gene, segregating with the Stocco dos Santos XLMR syndrome in a large four-generation pedigree but absent in >1,000 control X-chromosomes. Among other phenotypic characteristics, the affected males in this family present with severe MR, delayed or no speech, seizures and hyperactivity. Molecular studies of hKIAA1202 determined its genomic organisation, its expression throughout the brain and the regulation of expression of its mouse homologue during development. Transient expression of the wild-type KIAA1202 protein in HeLa cells showed partial colocalisation with the F-actin based cytoskeleton. On the basis of its domain structure, we argue that hKIAA1202 is a new member of the APX/Shroom protein family. Members of this family contain a PDZ and two ASD domains of unknown function and have been shown to localise at the cytoskeleton, and play a role in neurulation, cellular architecture, actin remodelling and ion channel function. Our results suggest that hKIAA1202 may be important in cognitive function and/or development.