Linking macroscopic with microscopic neuroanatomy using synthetic neuronal populations

A realistic morpho-anatomical connection strategy for modelling full-scale point-neuron microcircuits

The modeling of extended microcircuits is emerging as an effective tool to simulate the neurophysiological correlates of brain activity and to investigate brain dysfunctions. However, for specific networks, a realistic modeling approach based on the combination of available physiological, morphological and anatomical data is still an open issue. One of the main problems in the generation of realistic networks lies in the strategy adopted to build network connectivity. Here we propose a method to implement a neuronal network at single cell resolution by using the geometrical probability volumes associated with pre- and postsynaptic neurites. This allows us to build a network with plausible connectivity properties without the explicit use of computationally intensive touch detection algorithms using full 3D neuron reconstructions. The method has been benchmarked for the mouse hippocampus CA1 area, and the results show that this approach is able to generate full-scale brain networks at single cell resolution that are in good agreement with experimental findings. This geometric reconstruction of axonal and dendritic occupancy, by effectively reflecting morphological and anatomical constraints, could be integrated into structured simulators generating entire circuits of different brain areas facilitating the simulation of different brain regions with realistic models.

Reconstruction of the Hippocampus

The hippocampus is a widely studied brain region thought to play an important role in higher cognitive functions such as learning, memory, and navigation. The amount of data on this region increases every day and delineates a complex and fragmented picture, but an integrated understanding of hippocampal function remains elusive. Computational methods can help to move the research forward, and reconstructing a full-scale model of the hippocampus is a challenging yet feasible task that the research community should undertake.

In this chapter, we present strategies for reconstructing a large-scale model of the hippocampus. Based on a previously published approach to reconstruct and simulate brain tissue, which is also explained in Chap. 10, we discuss the characteristics of the hippocampus in the light of its special anatomical and physiological features, data availability, and existing large-scale hippocampus models. A large-scale model of the hippocampus is a compound model of several building blocks: ion channels, morphologies, single cell models, connections, synapses. We discuss each of those building blocks separately and discuss how to merge them back and simulate the resulting network model.

Generation of Granule Cell Dendritic Morphologies by Estimating the Spatial Heterogeneity of Dendritic Branching

Biological realism of dendritic morphologies is important for simulating electrical stimulation of brain tissue. By adding point process modeling and conditional sampling to existing generation strategies, we provide a novel means of reproducing the nuanced branching behavior that occurs in different layers of granule cell dendritic morphologies. In this study, a heterogeneous Poisson point process was used to simulate branching events. Conditional distributions were then used to select branch angles depending on the orthogonal distance to the somatic plane. The proposed method was compared to an existing generation tool and a control version of the proposed method that used a homogeneous Poisson point process. Morphologies were generated with each method and then compared to a set of digitally reconstructed neurons. The introduction of a conditionally dependent branching rate resulted in the generation of morphologies that more accurately reproduced the emergent properties of dendritic material per layer, Sholl intersections, and proximal passive current flow. Conditional dependence was critically important for the generation of realistic granule cell dendritic morphologies.

Heterogeneities in intrinsic excitability and frequency-dependent response properties of granule cells across the blades of the rat dentate gyrus

The dentate gyrus (DG), the input gate to the hippocampus proper, is anatomically segregated into three different sectors, namely, the suprapyramidal blade, the crest region, and the infrapyramidal blade. Although there are well-established differences between these sectors in terms of neuronal morphology, connectivity patterns, and activity levels, differences in electrophysiological properties of granule cells within these sectors have remained unexplored. Here, employing somatic whole cell patch-clamp recordings from the rat DG, we demonstrate that granule cells in these sectors manifest considerable heterogeneities in their intrinsic excitability, temporal summation, action potential characteristics, and frequency-dependent response properties. Across sectors, these neurons showed positive temporal summation of their responses to inputs mimicking excitatory postsynaptic currents and showed little to no sag in their voltage responses to pulse currents. Consistently, the impedance amplitude profile manifested low-pass characteristics and the impedance phase profile lacked positive phase values at all measured frequencies and voltages and for all sectors. Granule cells in all sectors exhibited class I excitability, with broadly linear firing rate profiles, and granule cells in the crest region fired significantly fewer action potentials compared with those in the infrapyramidal blade. Finally, we found weak pairwise correlations across the 18 different measurements obtained individually from each of the three sectors, providing evidence that these measurements are indeed reporting distinct aspects of neuronal physiology. Together, our analyses show that granule cells act as integrators of afferent information and emphasize the need to account for the considerable physiological heterogeneities in assessing their roles in information encoding and processing.NEW & NOTEWORTHY We employed whole cell patch-clamp recordings from granule cells in the three subregions of the rat dentate gyrus to demonstrate considerable heterogeneities in their intrinsic excitability, temporal summation, action potential characteristics, and frequency-dependent response properties. Across sectors, granule cells did not express membrane potential resonance, and their impedance profiles lacked inductive phase leads at all measured frequencies. Our analyses also show that granule cells manifest class I excitability characteristics, categorizing them as integrators of afferent information.

Spatial registration of neuron morphologies based on maximization of volume overlap

BACKGROUND

Morphological features are widely used in the study of neuronal function and pathology. Invertebrate neurons are often structurally stereotypical, showing little variance in gross spatial features but larger variance in their fine features. Such variability can be quantified using detailed spatial analysis, which however requires the morphologies to be registered to a common frame of reference.

RESULTS

We outline here new algorithms - Reg-MaxS and Reg-MaxS-N - for co-registering pairs and groups of morphologies, respectively. Reg-MaxS applies a sequence of translation, rotation and scaling transformations, estimating at each step the transformation parameters that maximize spatial overlap between the volumes occupied by the morphologies. We test this algorithm with synthetic morphologies, showing that it can account for a wide range of transformation differences and is robust to noise. Reg-MaxS-N co-registers groups of more than two morphologies by iteratively calculating an average volume and registering all morphologies to this average using Reg-MaxS. We test Reg-MaxS-N using five groups of morphologies from the Droshophila melanogaster brain and identify the cases for which it outperforms existing algorithms and produce morphologies very similar to those obtained from registration to a standard brain atlas.

CONCLUSIONS

We have described and tested algorithms for co-registering pairs and groups of neuron morphologies. We have demonstrated their application to spatial comparison of stereotypic morphologies and calculation of dendritic density profiles, showing how our algorithms for registering neuron morphologies can enable new approaches in comparative morphological analyses and visualization.

A Density-Driven Method for the Placement of Biological Cells Over Two-Dimensional Manifolds

We introduce a graphical method originating from the computer graphics domain that is used for the arbitrary placement of cells over a two-dimensional manifold. Using a bitmap image whose luminance provides cell density, this method guarantees a discrete distribution of the positions of the cells respecting the local density. This method scales to any number of cells, allows one to specify arbitrary enclosing shapes and provides a scalable and versatile alternative to the more classical assumption of a uniform spatial distribution. The method is illustrated on a discrete homogeneous neural field, on the distribution of cones and rods in the retina and on the neural density of a flattened piece of cortex.

T2N as a new tool for robust electrophysiological modeling demonstrated for mature and adult-born dentate granule cells

Compartmental models are the theoretical tool of choice for understanding single neuron computations. However, many models are incomplete, built ad hoc and require tuning for each novel condition rendering them of limited usability. Here, we present T2N, a powerful interface to control NEURON with Matlab and TREES toolbox, which supports generating models stable over a broad range of reconstructed and synthetic morphologies. We illustrate this for a novel, highly detailed active model of dentate granule cells (GCs) replicating a wide palette of experiments from various labs. By implementing known differences in ion channel composition and morphology, our model reproduces data from mouse or rat, mature or adult-born GCs as well as pharmacological interventions and epileptic conditions. This work sets a new benchmark for detailed compartmental modeling. T2N is suitable for creating robust models useful for large-scale networks that could lead to novel predictions. We discuss possible T2N application in degeneracy studies.

Semi-automatic 3D morphological reconstruction of neurons with densely branching morphology: Application to retinal AII amacrine cells imaged with multi-photon excitation microscopy

BACKGROUND

Accurate reconstruction of the morphology of single neurons is important for morphometric studies and for developing compartmental models. However, manual morphological reconstruction can be extremely time-consuming and error-prone and algorithms for automatic reconstruction can be challenged when applied to neurons with a high density of extensively branching processes.

NEW METHOD

We present a procedure for semi-automatic reconstruction specifically adapted for densely branching neurons such as the AII amacrine cell found in mammalian retinas. We used whole-cell recording to fill AII amacrine cells in rat retinal slices with fluorescent dyes and acquired digital image stacks with multi-photon excitation microscopy. Our reconstruction algorithm combines elements of existing procedures, with segmentation based on adaptive thresholding and reconstruction based on a minimal spanning tree. We improved this workflow with an algorithm that reconnects neuron segments that are disconnected after adaptive thresholding, using paths extracted from the image stacks with the Fast Marching method.

RESULTS

By reducing the likelihood that disconnected segments were incorrectly connected to neighboring segments, our procedure generated excellent morphological reconstructions of AII amacrine cells.

COMPARISON WITH EXISTING METHODS

Reconstructing an AII amacrine cell required about 2h computing time, compared to 2-4days for manual reconstruction. To evaluate the performance of our method relative to manual reconstruction, we performed detailed analysis using a measure of tree structure similarity (DIADEM score), the degree of projection area overlap (Dice coefficient), and branch statistics.

CONCLUSIONS

We expect our procedure to be generally useful for morphological reconstruction of neurons filled with fluorescent dyes.

AII amacrine cells: quantitative reconstruction and morphometric analysis of electrophysiologically identified cells in live rat retinal slices imaged with multi-photon excitation microscopy

AII amacrine cells have been found in all mammalian retinas examined and play an important role for visual processing under both scotopic and photopic conditions. Whereas ultrastructural investigations have provided a detailed understanding of synaptic connectivity, there is little information available with respect to quantitative properties and variation of cellular morphology. Here, we performed whole-cell recordings from AII amacrine cells in rat retinal slices and filled the cells with fluorescent dyes. Multi-photon excitation microscopy was used to acquire image stacks and after deconvolution, we performed quantitative morphological reconstruction by computer-aided manual tracing. We reconstructed and performed morphometric analysis on 43 AII amacrine cells, with a focus on branching pattern, dendritic lengths and diameters, surface area, and number and distribution of dendritic varicosities. Compared to previous descriptions, the most surprising result was the considerable extent of branching, with the maximum branch order ranging from approximately 10-40. We found that AII amacrine cells conform to a recently described general structural design principle for neural arbors, where arbor density decreases proportionally to increasing territory size. We confirmed and quantified the bi-stratified morphology of AII amacrine cells by analyzing the arborizations as a function of retinal localization or with Sholl spheres. Principal component and cluster analysis revealed no evidence for morphological subtypes of AII amacrines. These results establish a database of morphometric properties important for studies of development, regeneration, degeneration, and disease processes, as well as a workflow compatible with compartmental modeling.

Interneuronal mechanisms of hippocampal theta oscillations in a full-scale model of the rodent CA1 circuit

The hippocampal theta rhythm plays important roles in information processing; however, the mechanisms of its generation are not well understood. We developed a data-driven, supercomputer-based, full-scale (1:1) model of the rodent CA1 area and studied its interneurons during theta oscillations. Theta rhythm with phase-locked gamma oscillations and phase-preferential discharges of distinct interneuronal types spontaneously emerged from the isolated CA1 circuit without rhythmic inputs. Perturbation experiments identified parvalbumin-expressing interneurons and neurogliaform cells, as well as interneuronal diversity itself, as important factors in theta generation. These simulations reveal new insights into the spatiotemporal organization of the CA1 circuit during theta oscillations.

Interneuronal mechanisms of hippocampal theta oscillations in a full-scale model of the rodent CA1 circuit

The hippocampal theta rhythm plays important roles in information processing; however, the mechanisms of its generation are not well understood. We developed a data-driven, supercomputer-based, full-scale (1:1) model of the rodent CA1 area and studied its interneurons during theta oscillations. Theta rhythm with phase-locked gamma oscillations and phase-preferential discharges of distinct interneuronal types spontaneously emerged from the isolated CA1 circuit without rhythmic inputs. Perturbation experiments identified parvalbumin-expressing interneurons and neurogliaform cells, as well as interneuronal diversity itself, as important factors in theta generation. These simulations reveal new insights into the spatiotemporal organization of the CA1 circuit during theta oscillations.

DOI:http://dx.doi.org/10.7554/eLife.18566.001

Adult-born dentate granule cells show a critical period of dendritic reorganization and are distinct from developmentally born cells

Adult-born dentate granule cells (abGCs) exhibit a critical developmental phase during function integration. The time window of this phase is debated and whether abGCs become indistinguishable from developmentally born mature granule cells (mGCs) is uncertain. We analyzed complete dendritic reconstructions from abGCs and mGCs using viral labeling. AbGCs from 21-77 days post intrahippocampal injection (dpi) exhibited comparable dendritic arbors, suggesting that structural maturation precedes functional integration. In contrast, significant structural differences were found compared to mGCs: AbGCs had more curved dendrites, more short terminal segments, a different branching pattern, and more proximal terminal branches. Morphological modeling attributed these differences to developmental dendritic pruning and postnatal growth of the dentate gyrus. We further correlated GC morphologies with the responsiveness to unilateral medial perforant path stimulation using the immediate-early gene Arc as a marker of synaptic activation. Only abGCs at 28 and 35 dpi but neither old abGCs nor mGCs responded to stimulation with a remodeling of their dendritic arbor. Summarized, abGCs stay distinct from mGCs and their dendritic arbor can be shaped by afferent activity during a narrow critical time window.

Dendritic and Axonal Wiring Optimization of Cortical GABAergic Interneurons

The way in which a neuronal tree expands plays an important role in its functional and computational characteristics. We aimed to study the existence of an optimal neuronal design for different types of cortical GABAergic neurons. To do this, we hypothesized that both the axonal and dendritic trees of individual neurons optimize brain connectivity in terms of wiring length. We took the branching points of real three-dimensional neuronal reconstructions of the axonal and dendritic trees of different types of cortical interneurons and searched for the minimal wiring arborization structure that respects the branching points. We compared the minimal wiring arborization with real axonal and dendritic trees. We tested this optimization problem using a new approach based on graph theory and evolutionary computation techniques. We concluded that neuronal wiring is near-optimal in most of the tested neurons, although the wiring length of dendritic trees is generally nearer to the optimum. Therefore, wiring economy is related to the way in which neuronal arborizations grow irrespective of the marked differences in the morphology of the examined interneurons.

Microcircuits in Epilepsy: Heterogeneity and Hub Cells in Network Synchronization

Epilepsy is a complex disorder involving neurological alterations that lead to the pathological development of spontaneous, recurrent seizures. For decades, seizures were thought to be largely repetitive, and had been examined at the macrocircuit level using electrophysiological recordings. However, research mapping the dynamics of large neuronal populations has revealed that seizures are not simply recurrent bursts of hypersynchrony. Instead, it is becoming clear that seizures involve a complex interplay of different neurons and circuits. Herein, we will review studies examining microcircuit changes that may underlie network hyperexcitability, discussing observations from network theory, computational modeling, and optogenetics. We will delve into the idea of hub cells as pathological centers for seizure activity, and will explore optogenetics as a novel avenue to target and treat pathological circuits. Finally, we will conclude with a discussion on future directions in the field.