The neocortical microcircuit collaboration portal: a resource for rat somatosensory cortex

We have established a multi-constraint, data-driven process to digitally reconstruct, and simulate prototypical neocortical microcircuitry, using sparse experimental data. We applied this process to reconstruct the microcircuitry of the somatosensory cortex in juvenile rat at the cellular and synaptic levels. The resulting reconstruction is broadly consistent with current knowledge about the neocortical microcircuit and provides an array of predictions on its structure and function. To engage the community in exploring, challenging, and refining the reconstruction, we have developed a collaborative, internet-accessible facility—the Neocortical Microcircuit Collaboration portal (NMC portal; https://bbp.epfl.ch/nmc-portal). The NMC portal allows users to access the experimental data used in the reconstruction process, download cellular and synaptic models, and analyze the predicted properties of the microcircuit: six layers, ~31,000 neurons, 55 morphological types, 11 electrical types, 207 morpho-electrical types, 1941 unique synaptic connection types between neurons of specific morphological types, predicted properties for the anatomy and physiology of ~40 million intrinsic synapses. It also provides data supporting comparison of the anatomy and physiology of the reconstructed microcircuit against results in the literature. The portal aims to catalyze consensus on the cellular and synaptic organization of neocortical microcircuitry (ion channel, neuron and synapse types and distributions, connectivity, etc.). Community feedback will contribute to refined versions of the reconstruction to be released periodically. We consider that the reconstructions and the simulations they enable represent a major step in the development of in silico neuroscience.

Introduction

The Blue Brain Project has developed a unifying process for the pragmatic integration of available data and knowledge into in silico reconstructions of neuronal microcircuits. We have used this process to reconstruct the microcircuitry of developing rat [postnatal days (P) 13–16] somatosensory cortex from sparse data on its cellular and synaptic organization (Markram et al., 2015). The resulting reconstruction is made up of component data covering multiple levels of detail: ion channels, synapses, neurons, and the entire microcircuit. The reconstruction evolves in successive cycles of reconstruction-validation-experiment-reconstruction, and will continuously integrate new properties and principles as they become available.

The anatomy and physiology of the reconstruction can be studied in the same way as those of a block of neocortical tissue. To enable such studies, we have developed an internet-accessible public resource—The NMC Portal. The portal provides access to the data, literature, and models used in the reconstruction of P13-16 rat somatosensory cortex, together with interactive tools to browse and query its detailed anatomy and physiology. The resource also provides means for the community to explore, analyze, and annotate the properties of the reconstruction, to compare them against results from the literature, and to provide feedback for future releases.

The Reconstructed Neocortical Microcircuit

The reconstructed microcircuit represents six cortical layers with an overall thickness of ~2.1 mm, and a volume of ~0.29 mm³. It is constituted by ~31,000 neurons belonging to 55 morphological types (m-types), 11 electrical types (e-types), and 207 morpho-electrical types (me-types); ~7.5 million connections belonging to 1941 unique m-type specific connection types; and ~40 million synapses belonging to six synapse types (s-types) (Markram et al., 2015). The overall dimensions, individual layer thicknesses, and the ratio of excitatory to inhibitory neurons in the reconstruction are consistent with results from previous studies of rat barrel cortex (Meyer et al., 2010, 2011; Wimmer et al., 2010). However, the total number of neurons in the reconstruction based on measurements in juvenile (P13-16) rat somatosensory cortex is around 1.5 times higher than in adult (P27) rat barrel cortex (Meyer et al., 2010).

The NMC Portal—an Overview

The portal allows users to access the experimental data used in the reconstruction process, download cellular and synaptic models, and analyze the predicted properties of the microcircuit. The portal is structured into the following sections (see Supplementary Figure 1):

Microcircuit

An interactive browser across three levels: layers, neurons, and synapses. The layer level provides data on thickness, neuron, and synapse densities, and the distributions of neuron and synapse types for each of the six layers of the reconstruction. The neuronal level describes the anatomical and physiological properties of morphological, electrical and morpho-electrical neuron types (m-, e-, and me-types, respectively). The synaptic level represents the anatomical and physiological properties of synaptic connections between specific pre-post combinations of m-types, and the complete map of intrinsic input and output synapses (three excitatory, and three inhibitory synaptic types), from and to neurons of different types.

Literature Consistency

Comparisons between the morphological, molecular, electrical, synaptic, and physiology properties of the reconstruction and its overall circuit anatomy against the published literature.

Experimental Data

Experimental datasets used in the reconstruction process.

Videos

Animations of simulated microcircuit activity under a variety of experimental conditions.

Images

Images illustrating key steps in the reconstruction process.

Tools

Tools for analyzing, and simulating the reconstructed microcircuit.

Downloads

Models of neurons, ion channels, and synapses for the NEURON simulation environment.

The Microcircuit

The microcircuit section of the portal contains data computed during the reconstruction process, using methods described in (Markram et al., 2015). The data is presented in fact sheets, each representing a different level of biological organization.

The fact sheet for the whole reconstructed microcircuit provides an integrated view of its dimensions, layer-wise distributions and densities of neurons, total number of morphological types (m-types), electrical types (e-types), morpho-electrical types (me-types), and synapse types (s-types), numbers of intrinsic, and extrinsic synaptic connections and associated synapses, and the number of unique synaptic connection types between neurons of specific source and target m-types (see Figure 1).

FIGURE 1

Figure 1. An overview of the reconstructed microcircuit—facts and figures. Top left: facts and figures on the neuronal anatomy of the reconstructed microcircuit. Top right: overview of neuronal physiology. Bottom left: facts and figures on synaptic anatomy. Bottom right: overview of synaptic physiology.

Layers

The layer level fact sheets (see Figure 2) provide a unified picture of the anatomical and physiological properties of each layer in the reconstructed microcircuit as outlined below.

1. The total number of neurons in each designated layer and for each morphological type (m-type) in a layer are given (see Figure 2B), where the naming of m-types was based on the most common names used in previous studies (Larkman, 1991; Kawaguchi and Kubota, 1997; DeFelipe et al., 2002; Wang et al., 2006; Lübke, 2003; Ascoli et al., 2008; Romand et al., 2011).

2. Total axonal length and volume, total dendritic length and volume, synapse density, total number of morphology-specific synaptic connections are provided, and categorized into: geometrically possible and viable pathways, excitatory and inhibitory pathways, and intra- and inter-laminar pathways (Thomson and Deuchars, 1997; Somogyi et al., 1998; Feldmeyer et al., 1999, 2002; Gupta et al., 2000; Wang et al., 2002; Silberberg and Markram, 2007).

3. The proportion of different electrical types (e-types) for each layer is given (see Figure 2C for an example of e-type proportions in layer 5).

4. The mapping of inhibitory m-types to the main calcium-binding proteins [parvalbumin (PV), calbindin (CB), and calretinin (CR)], and neuropeptides [somatostatin (SOM), vasoactive intestinal polypeptide (VIP), neuropeptide Y (NPY), and cholecystokinin (CCK)] they express is given (see Figure 2D) as described previously (DeFelipe, 1993; Gonchar and Burkhalter, 1997; Kawaguchi and Kubota, 1997; Kawaguchi and Kondo, 2002; Toledo-Rodriguez et al., 2005; Gonchar et al., 2008; O'Connor et al., 2009; Meyer et al., 2011; Santana et al., 2013).

FIGURE 2

Figure 2. Layer fact sheet. (A) The pathway navigator predicting the post-synaptic target layers of all m-types in layer 5. The thickness of a connecting ribbon indicates the number of possible pathways (connections or synapses, based on the user selection in the top left) of a given layer 5 m-type to the recipient layer. (B) Overview of the m-type composition in layer 5, morphometrics, and synaptic anatomy. (C) Proportional composition of different e-types in layer 5. (D) Map of commonly expressed calcium binding proteins and neuropeptides within inhibitory neurons of different m-types in layer 5. The thickness of the connecting ribbon indicates the prevalence of specific markers.

Morphological Types (m-types)

The m-type fact sheets (see Figure 3) provide an anatomical overview of total axonal and dendritic lengths, and volumes for individual m-types (see Figure 3B1). Previous studies have established that neurons of particular m-types display diverse electrical behavior (Kawaguchi and Kubota, 1997; Gupta et al., 2000; Markram et al., 2004; Toledo-Rodriguez et al., 2004). Therefore, the m-type fact sheets also list the different e-types associated with the m-type (see Figure 3B2). In addition, the fact sheets provide detailed anatomical characterizations of the axonal and dendritic properties of neurons belonging to the m-type. These include section bifurcation angles (Figures 3C,G), total lengths (see Figures 3D,H), individual section lengths (see Figures 3E,I), and volumes (Figures 3F,J).

FIGURE 3

Figure 3. m-type fact sheet. (A) Predicted post-synaptic targets for neurons belonging to the L5_SBC m-type; (B1) anatomical fact sheet for all neurons of the L5_SBC m-type. (B2) e-types expressed in neurons belonging to the L5_SBC m-type (C) Histogram of axonal section bifurcation angles compiled from all reconstructed L5_SBC morphologies (N = 106). (D) Histogram of total axonal lengths. (E) Histogram of individual axonal section lengths. (F) Histogram of axonal volumes. (G) Histogram of dendritic section bifurcation angles. (H) Histogram of total dendritic lengths. (I) Histogram of individual axonal section lengths. (J) Histogram of dendritic volumes.

Morpho-Electrical Types (me-types)

The me-type fact sheets provide an in-depth description of the anatomical and physiological properties of each of the 207 individual me-types in the reconstruction, and detailed data for five exemplar neurons of each me-type (Figure 4; see also Figure 4E for a single exemplar): total length and volume, maximal section length, maximal branch order, and soma diameter (see Figure 4A). The physiological characterization includes standard metrics for intrinsic membrane properties (Connors and Gutnick, 1990; Kasper et al., 1994; Zhu, 2000): resting potential, input resistance, and membrane time constant (see Figure 4B). The fact sheets also provide a list of relevant ion channels models for specific me-types, and data on the distribution of ion channels along neuronal arbors, identified from the literature (see Figure 4C) (Stuart and Sakmann, 1994; Korngreen and Sakmann, 2000; Kole et al., 2006). These data make it possible to model the electrical firing patterns of specific me-types, as shown in Figure 4D, which illustrates modeled electrical responses to somatic step current injections of varying intensities. The fact sheets link to the relevant model packages, which are designed for use in the NEURON simulation environment.

FIGURE 4

Figure 4. me-type fact sheet. (A) Anatomical properties of all neuron models belonging to the L5_SBC-dNAC me-type. (B) Physiological properties of all neuron models belonging to the L5_SBC-dNAC me-type. (C) Density distributions of ion channels used in neuron models. (D) Representative response traces in a L5_SBC_dNAC me-type model to different intensities of somatically injected step currents. (E) 3D reconstructed morphology of an exemplar L5_SBC neuron. (F) Predicted dendritic attenuation of EPSPs. (G) Predicted attenuation of the back-propagating AP.

Neuronal dendrites not only funnel synaptic input toward the axon, but also sustain action potentials (AP) back-propagating from the axon (Stuart and Sakmann, 1994). However, experimental data characterizing the active dendritic properties of neocortical neurons is only available for a small subset of me-types (Larkum et al., 2007, 2009; Nevian et al., 2007; Ledergerber and Larkum, 2010). The reconstruction predicts these values. The me-type fact sheets provide predictions (movies) of dendritic back-propagating action potentials (bAP), and the attenuation of post-synaptic potentials (PSP) for all 207 me-types (see Figures 4F,G).

Connection Types

The 55 m-types constituting the microcircuit could theoretically form 3025 morphology-specific connection types. Consideration of the overlap of axons and dendrites, or lack thereof, in the reconstruction predicts that 1941 of these connection types are geometrically viable. Although, only around 1% of such connection types have been characterized experimentally (Thomson et al., 1993; Markram et al., 1997; Reyes et al., 1998; Feldmeyer, 2006), the reconstruction makes it possible to predict their properties, which are presented in connection type fact sheets. Each fact sheet shows the predicted map of synaptic connections from neurons of one m-type to neurons of other m-types (Supplementary Figure 2A), together with data for their predicted anatomy, and physiology, (Supplementary Figure 2B).

The anatomical data includes total synapse count, average connection probability, number of common neighbors (see Perin et al., 2011), number of convergent and divergent neurons, and the mean number of synapses per connection (synapses/connection) for all connections between pairs of neurons belonging to a specific connection type (Supplementary Figure 2B). It also provides more detailed data, including a graphical representation of an exemplar pair of synaptically connected neurons (see Figure 5A), overlapping axo-dendritic density clouds for all pre- post-synaptic neuron pairs that belong to the connection type (Figure 5B), the overlaid axonal and dendritic density clouds for individual pre-synaptic and post-synaptic morphologies (Figures 5C,E, respectively), the axogram for the morphology of the pre-synaptic neuron in the exemplar pair (as shown in Figure 5A), predicting the axonal location of synaptic contacts (Figure 5D), and the dendrogram for the morphology of the post-synaptic neuron in the exemplar pair, predicting the dendritic location of synaptic contacts (Figure 5F).

FIGURE 5

Figure 5. Overview of in silico synaptic anatomy. (A) An exemplar in silico pair of synaptically connected L23_SBC (right, yellow) and L23_PC (left, black) neurons. The synaptic contacts mediating the connection are shown as yellow circles on the L23_PC. (B) Overlap of axo-dendritic density clouds of all pre-synaptic L23_SBCs and post-synaptic L23_PCs. (C) Overlap of axonal density clouds of all pre-synaptic L23_SBCs. (D) Axogram of the pre-synaptic L23_SBC showing the location of synaptic contacts (red circles) along the axon (in blue). (E) Overlap of dendritic density clouds of all pre-synaptic L23_PCs. (F) Dendrogram of the post-synaptic L23_PC, showing the location of synaptic contacts (blue circles) along the dendrites (in red).

Additionally, the fact sheets provide detailed statistical distribution for the number of synapses per connection (Figure 6A), the total number of synapses from all pre-synaptic neurons of a given type to a post-synaptic neuron of a given type (Figure 6B), the number of post-synaptic neurons of a given type targeted by individual neurons of a given type (i.e., neuronal divergence; Figure 6C, left), the number of pre-synaptic neurons of a given type, targeting individual post-synaptic neurons of a given type (i.e., neuronal convergence; Figure 6C, right), the fraction of efferent synapses from a given neuron type targeting neurons of a specific type (i.e., synaptic divergence; Figure 6D, left), the fraction of afferent synapses from neurons of a specific type targeting neurons of a specific type (i.e., synaptic convergence; Figure 6D, right), patterns of axonal and dendritic synaptic innervation based on branch order (Figure 6E; left and right, respectively), and path distance from the soma (Figure 6F; left and right, respectively).

FIGURE 6

Figure 6. Characterization of in silico synaptic anatomy. (A) Distribution of number of synapses per connection for pairs of synaptically connected L23_SBC to L23_PC neurons. (B) Distribution of the total number of synapses from all L23_SBCs to L23_PCs. (C) Left: neuronal divergence; number of L23_PCs targeted by individual L23_SBCs. Right: neuronal convergence - number of L23_SBCs targeting individual L23_PCs. (D) Left: synaptic divergence - fraction of all synapses formed by L23_SBCs that target L23_PCs. Right: synaptic convergence; fraction of all synapses formed onto L23_PCs that originate from L23_SBCs. (E) Axonal (left) and dendritic (right) innervation patterns, in terms of branch order of synaptic contacts. (F) Same as (E), but in terms of geometrical distance of synaptic contacts.

The physiological data provides predictions at the level of single synapses and connections (Supplementary Figure 2B). Each synapse is assigned to one of six s-types (Gupta et al., 2000; Wang et al., 2006). Individual synapses are characterized by their peak conductance (g_syn; in nS), neurotransmitter release probability (U), time constant for recovery from depression (D; in ms), time constant for recovery from facilitation (F; in ms), and the predicted correction factor necessary to account for dendritic filtering (space-clamp artifacts). Physiological data at the connection level summarizes the kinetics and time course of post-synaptic potentials (PSPs). The data provided include onset latencies, peak amplitudes, 20–80% rise times, decay time constants, transmission failures, and the coefficient of variation of PSP amplitudes (c.v.; defined as the ratio of the standard deviation to the mean of PSP amplitude) (see Ramaswamy and Markram, 2015).

Additionally, the fact sheets provide data on the trial-to-trial response variability of unitary PSPs for given types of post-synaptic neurons (Figure 7A for the exemplar pair shown in Figure 5A), synaptic response to pre-synaptic stimulation at 10, 30, 50, and 70 Hz (Figure 7C, clockwise from top left), and the relationship between the frequency of pre-synaptic stimulation and PSP amplitude. These data illustrate the 1/f law for synapses found in previous experimental studies [(Figure 7B); (Tsodyks and Markram, 1997); also see (Ramaswamy and Markram, 2015)]. The fact sheets also provide data from virtual recordings of local dendritic PSPs (Figure 7D). It is currently not possible to obtain comparable data through experiments alone.

FIGURE 7

Figure 7. Overview of in silico synaptic physiology. (A) Average time course and amplitude of in silico unitary PSPs measured in a post-synaptic L23_PC by evoked unitary APs in a pre-synaptic L23_SBC (for the same pair of neurons shown in Figure 5A); blue line to the left indicates the pre-synaptic AP; pink traces, individual post-synaptic responses across 30 trials; red trace, average of 30 trials. (B) The 1/f rule demonstrating an inverse relationship between stationary PSPs and the frequency of pre-synaptic stimulation; the solid line in red shows the inverse relationship with frequency, and the cross shows the limiting frequency. (C) Frequency dependence of synaptic transmission. Clockwise from top left: in silico PSPs evoked in the post-synaptic L23_PC upon stimulating the pre-synaptic L23_SBC with pulse trains at frequencies of 10, 30, 70, and 50 Hz, respectively (for the same pair of neurons shown in Figure 5A). (D) The digital microcircuit reconstruction enables direct in silico recordings of single synaptic contacts on dendrites. For the synaptically connected pair (shown in Figure 5A), post-synaptic responses were recorded directly at synaptic locations on dendrites of the L23_PC (dendrites in black; soma in red); the connection was mediated by nine synaptic contacts (blue circles); dendritic recording sites are shown as light blue pipettes; corresponding PSPs are shown in blue.

The fact sheets also provide detailed statistical distributions for the kinetics and average time course of PSPs, including mean PSP amplitude (Figure 8A), 20–80% PSP rise time (Figure 8B), PSP onset latency (Figure 8C), PSP decay time constant (Figure 8D), transmission failures (Figure 8E), c.v. of PSP amplitude (Figure 8F), and the inverse relationship between transmission failures and c.v. of PSP amplitude against mean PSP amplitude [Figures 8G,H, respectively; also see (Markram et al., 1997; Feldmeyer et al., 1999)].

FIGURE 8

Figure 8. Characterization of in silico synaptic physiology. (A) Histogram of in silico PSP amplitudes for synaptically connected pairs of L23_SBC and L23_PC neurons (n = 100 pairs, sampled at inter-somatic distances ≤ 100 μm); the dashed line in black indicates the mean of the distribution; the error bar shows the standard deviation (SD). (B) Histogram of in silico 20–80% PSP rise times. (C) Histogram of in silico PSP onset latencies. (D) Histogram of in silico PSP decay time constants. (E) Histogram of transmission failures. (F) Histogram of the c.v. of PSP amplitudes. (G) Predicted inverse relationship between the rate of transmission failures and PSP amplitude; the reconstructed microcircuit predicted a decrease in failure rates with increasing PSP amplitudes, consistent with a binomial model of transmitter release. (H) Same as in (G), but for c.v. of PSP amplitudes.

The reconstruction makes it possible to predict the full complement of synaptic inputs and outputs for any given neuron. For each connection-type, a connection-type fact sheet provides predictions on the number and locations of afferent synaptic input to and efferent synaptic output from neurons belonging to specific m-types (see Figure 9A). These data categorize inputs to and outputs from all m-, e-, and s-types in the reconstructed microcircuit (see Figures 9B,C). Web links are provided to download NEURON model packages for exemplar neurons for each of the 55 m-types.

FIGURE 9

Figure 9. Predicted map of afferent synaptic input to, and efferent synaptic output from an exemplar L23_PC neuron. (A) Rendering of a L23_PC, predicting the map of afferent synapses color coded by m-type. (B,C) The intrinsic synaptic input and output, respectively from and to all m-, e-, and s-types across six layers for all L23_PCs in the reconstructed microcircuit.

Literature Consistency

The literature consistency contains papers considered in the validation process. For each paper we identified one possible evaluation criterion. Each paper is annotated to show the evaluation criterion, and the result of the validation (consistent, possibly consistent, inconsistent). Where a paper was not actually used, it is marked as not implemented, or not applicable. Additional annotations describe the species and age of the animals used in the study, the brain region, cortical layer, and cell type concerned, the key finding, and other relevant results. The PubMed ID of the publication (PubMed is a free resource that is developed and maintained by the National Center for Biotechnology Information, at the U.S. National Library of Medicine, located at the National Institutes of Health), a web link to the publication, the name of the contact author, and the names of the persons who undertook the literature search are also provided.

The evaluation criteria are based on the morphological, molecular, electrical, synaptic, circuit anatomical, and physiological properties of the reconstructed microcircuit. For ease of presentation, papers are grouped according to the class of properties considered in their evaluation criteria.

Morphological Properties

This sub-section presents 64 papers used for the validation of experimentally reconstructed neuron morphologies. The papers provide data for Sholl statistics, dendritic and axonal lengths, dendritic and axonal volumes, soma diameters, and dendritic and axonal segment branch order and path distance statistics. Overall, the properties of the morphologies used in the reconstruction were consistent with the results reported in 17 papers, possibly consistent with 10 papers and inconsistent with nine papers. Six papers whose data were not comparable with the reconstruction (different species, different brain regions etc.) were excluded from the validation process.

Molecular Properties

During the validation process, gene expression data associated with specific neuron types was used to perform in silico staining. This sub-section presents six papers providing the necessary data. The papers present data on the expression of the main calcium binding proteins and neuropeptides expressed by neurons of different types.

Electrical Properties

This sub-section presents 77 papers describing the electrical properties of different cortical neuron types. Data are provided for ion channel kinetics, electrical firing patterns, and passive membrane properties, including the resting membrane potential, membrane time constant, and input resistance. Other data describe electrophysiological features including AP amplitude, AP half-width, firing frequency, inter-spike interval, and dendritic properties including bAP and EPSP/IPSP attenuation.

Synaptic Properties

This sub-section contains 81 publications used to validate the anatomical and physiological properties of synaptic connection types in the reconstruction. The data provided includes spine and bouton densities for different neuron types, as well as data for pairs of neurons belonging to specific m-types: synaptic innervation patterns, number of synaptic contacts per connection, connection probabilities, transmitter release probabilities, AMPA, NMDA, kainate, metabotropic glutamate, GABA_A, and GABA_B receptor peak conductances, EPSC/IPSC properties (onset latency, amplitude, rise time, and decay time constant), EPSP/IPSP properties (onset latency, amplitude, rise time, and decay time constant), transmission failure rates, c.v. of EPSP/IPSP amplitudes, short term depression and facilitation, and tonic excitation and inhibition.

Microcircuit Anatomical Properties

This sub-section presents 43 papers used to validate the anatomical properties of the reconstructed microcircuit. The properties considered include microcircuit thickness, individual layer thicknesses, layer-wise distributions and densities of neurons, total axon and dendritic lengths, total number of synapses and their densities, and layer-wise distributions of bouton densities from thalamocortical axons.

Microcircuit Physiological Properties

This sub-section contains 47 papers used to assess the physiological properties of the reconstructed microcircuit. Relevant properties include the frequency of network oscillations, patterns of propagation of electrical activity across different cortical layers, the physiological properties of thalamocortical synapses, and the balance of excitatory and inhibitory synaptic conductances in different neuron types due to network activity.

Experimental Data

Over the past two decades, Markram and colleagues have experimentally characterized the cellular, and synaptic microcircuitry of developing (P13-16) rat neocortex, where the anatomical, physiological and synaptic properties are not completely mature (Markram et al., 1997; Gupta et al., 2000; Wang et al., 2002, 2004, 2006; Toledo-Rodriguez et al., 2005; Silberberg and Markram, 2007; Berger et al., 2009). The portal provides access to the data sets used in the reconstruction process, together with descriptions of the experimental protocols used to generate the data. They include:

1. Measurements of individual layer heights.

• Data are provided in the form of microscopy images of NeuN (neuron-specific nuclear protein) stained coronal slices with annotations of individual layer extents, and spreadsheets summarizing measurements of layer heights.

2. Layer-wise distributions and densities of neurons.

• Data are given as microscopy images of NeuN stained slices with annotations of individual layer extents, and spreadsheets summarizing measurements of neuron counts across different layers.

3. Molecular characterization of single neurons based on the expression of calcium binding proteins [parvalbumin (PV), calbindin (CB), and calretinin (CR)], and neuropeptides [somatostatin (SOM), vasoactive intestinal polypeptide (VIP), neuropeptide Y (NPY), and cholecystokinin (CCK)].

• Data can be obtained as text files containing the expression patterns of calcium binding proteins, and neuropeptides for morphologically identified neurons.

4. More than 1000 morphological reconstructions of neurons of different m-types.

• Data are provided as text files containing 3D representations of neuronal morphologies reconstructed using Neurolucida (neuron tracing, reconstruction, and analysis software).

5. More than 1500 electrical recordings from neurons of different e-types.

• Data are given for current stimuli (in pA), and voltage responses (in mV) from whole-cell patch clamp recordings in single neurons acquired using IGOR Pro (scientific graphing, and data analysis software; WaveMetrics, Inc., Portland, OR, USA).

6. More than 5000 experiments from synaptically connected neurons.

• Data are provided for post-synaptic responses obtained from whole-cell patch clamp recordings in synaptically connected pairs of neurons.

Videos

The videos section provides movies of simulated spontaneous and evoked activity in virtual cortical slices, obtained from the reconstructed microcircuit under a variety of simulated experimental conditions. Movies of spontaneous activity show electrical activity in the reconstructed microcircuit at different levels of depolarization and extracellular Ca²⁺. Movies of evoked activity show electrical activity in response to the stimulation of thalamocortical fibers at different levels of extracellular Ca²⁺. Additional movies visualize the distribution of afferent synapses for a single exemplar neuron for each of the 55 different m-types, and predicted dendritic attenuation of EPSPs and bAPs for 5 exemplars for each of the 207 me-types.

Tools

The tools section contains three tools used for the analysis, repair and cloning of neuronal morphologies. NeuroM provides a range of morphometric analyses that allow a user to quantify properties of the axonal and dendritic morphologies of neurons, and provides features that are used in the classification of different neurons into one of 55 m-types. NeuroR utilizes the results of NeuroM to identify and repair arbors that were severed during slice preparation. NeuroC produces clones of the neurons repaired by NeuroR, introducing statistical variations in the arbors of each cloned neuron. This procedure makes it possible to generate a limitless number of unique instances of neurons belonging to a given m-type.

Downloads

The downloads section provides models of single neurons, synapses, and predicted maps of input-output synapses. Models are available either in the native format used in the reconstruction process, based on the NEURON simulation environment, or in the emerging NeuroML 2.0 (a model description language for computational neuroscience based on the extensible markup language–XML), and LEMS (Low Entropy Model Specification language, which aims to provide a compact, minimally redundant way of expressing models of biological systems). Models of individual neurons (in the NEURON simulation environment) can also be obtained from the me-type fact sheets [see Morpho-Electrical Types (me-types)]. Model packages for individual neurons contain a 3D reconstructed morphological model, models of ion channels and synapses, synaptic model parameter descriptions, and a template model of the electrical type. Additional helper scripts for the NEURON simulation environment are provided to instantiate a morphoelectrical neuron model, distribute ion channel mechanisms on axonal and dendritic arbors, and simulate electrophysiological and synaptic experiments.

Discussion

The NMC portal was developed to serve as an online resource for experimentalists and theoreticians, allowing them to explore, analyze, challenge, and refine the reconstructed microcircuit and to derive predictions to be validated against experiments. The reconstructed microcircuit is a first draft, which will be refined in successive cycles of reconstruction-validation-experiment-reconstruction. The interactive environment provided by the NMC portal will facilitate the integration of new datasets, contributing to the refinement of the reconstruction. For example, the total length of axons in the reconstruction falls short of experimental estimates, possibly because of the difficulty of reconstructing the thin terminal axonal segments of biocytin filled neurons in vitro. The shortfall can be overcome by in silico synthesis of axons, utilizing data from in vivo filled neurons (see Oberlaender et al., 2011, 2012; Narayanan et al., 2015). Other examples of data that can be used to refine the reconstruction are profiles of the density of ion channels along neuronal arbors, which provide constraints for the optimization of single neuron models. Additional datasets will help to fill in biological details that are not included in the current first draft reconstruction (e.g., data on gap junctions, cholinergic modulatory mechanisms, rules for activity dependent plasticity, extrasynaptic glutamate, and GABA receptors, kinetics of metabotropic glutamate receptors, non-synaptic transmission, autaptic connections etc.). We are currently developing further data integration techniques and quality control measures, which will ensure strict compatibility of different datasets with our own experimental data in terms of species, age, cortical area and region, slice orientation and thickness, composition of extracellular and intracellular solutions, recording temperature, and technique, liquid junction potential correction etc.

The NMC Portal is part of a broad trend toward the development of large repositories of anatomical and physiological data. Other examples include NeuroMorpho, ModelDB, NeuroElectro, Open Source Brain, and the recent Allen Brain Cell Types Database. The portal contributes to this trend by integrating experimental data, models of single neurons, ion channels, synapses and detailed reconstructions of the microcircuit, and validation data in a single interactive environment. The ultimate goal is to catalyze consensus on the cellular and synaptic organization of neocortical microcircuitry (ion channel, neuron and synapse types and distributions, connectivity, etc.) and to enable community-driven refinement of the reconstruction. Newer versions of the reconstruction will be released periodically, consisting of data structures consistent with the current version. We consider that these reconstructions and the simulations they enable represent a major step in the development of in silico neuroscience.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This work was supported by funding from the EPFL to the Laboratory of Neural Microcircuitry (LNMC), funding from the ETH Domain for the Blue Brain Project (BBP) and by funding to the Human Brain Project from the European Union Seventh Framework Program (FP7/2007–2013) under grant agreement no. 604102 (HBP). Additional funding came from The Gatsby Charitable Foundation, The Cajal Blue Brain project, Ministerio de Economia y Competitividad Spanish Ministry of Education and Science, and an EPFL-Hebrew University Collaborative Grant. In the years 2005–2009, the Blue Gene/L system was funded by the EPFL. Financial support for the subsequent CADMOS Blue Gene/P and Blue Gene/Q systems was provided by the Canton of Geneva, Canton of Vaud, Hans Wilsdorf Foundation, Louis-Jeantet Foundation, University of Geneva, University of Lausanne and École Polytechnique Fédérale de Lausanne. The BlueBrain IV BlueGene/Q system is financed by ETH Board Funding to the Blue Brain Project as a National Research Infrastructure and hosted at the Swiss National Supercomputing Center (CSCS). PG was supported by the Wellcome Trust (grant no. 101445). We thank Oren Amsalem, and Guy Eyal from the Hebrew University of Jerusalem, Lea Goetz, and Dr. Arnd Roth from University College London, Dr. Michele Giugliano from the University of Antwerp, and Drs. Costas Anastassiou, Anton Arkhipov, Stefan Mihalas, Susan Sunkin, and Prof. Christof Koch from the Allen Institute for Brain Science for insightful discussions and suggestions in improving the presentation of the NMC portal. We are grateful to Prof. Alex Thomson for many years of collaboration, support, and discussions.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fncir.2015.00044

References

Ascoli, G. A., Alonso-Nanclares, L., Anderson, S. A., Barrionuevo, G., Benavides-Piccione, R., Burkhalter, A., et al. (2008). Petilla terminology: nomenclature of features of GABAergic interneurons of the cerebral cortex. Nat. Rev. Neurosci. 9, 557–568. doi: 10.1038/nrn2402

PubMed Abstract | CrossRef Full Text | Google Scholar

Berger, T. K., Perin, R., Silberberg, G., and Markram, H. (2009). Frequency-dependent disynaptic inhibition in the pyramidal network: a ubiquitous pathway in the developing rat neocortex. J. Physiol. 587, 5411–5425. doi: 10.1113/jphysiol.2009.176552

PubMed Abstract | CrossRef Full Text | Google Scholar

Connors, B. W., and Gutnick, M. J. (1990). Intrinsic firing patterns of diverse neocortical neurons. Trends Neurosci. 13, 99–104. doi: 10.1016/0166-2236(90)90185-D

PubMed Abstract | CrossRef Full Text | Google Scholar

DeFelipe, J., Alonso-Nanclares, L., and Arellano, J. I. (2002). Microstructure of the neocortex: comparative aspects. J. Neurocytol. 31, 299–316. doi: 10.1023/A:1024130211265

PubMed Abstract | CrossRef Full Text | Google Scholar

DeFelipe, J. (1993). Neocortical neuronal diversity: chemical heterogeneity revealed by colocalization studies of classic neurotransmitters, neuropeptides, calcium-binding proteins, and cell surface molecules. Cereb. Cortex 3, 273–289. doi: 10.1093/cercor/3.4.273

PubMed Abstract | CrossRef Full Text | Google Scholar

Feldmeyer, D., Egger, V., Lübke, J., and Sakmann, B. (1999). Reliable synaptic connections between pairs of excitatory layer 4 neurones within a single “barrel” of developing rat somatosensory cortex. J. Physiol. 521(Pt 1), 169–190. doi: 10.1111/j.1469-7793.1999.00169.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Feldmeyer, D., Lübke, J., Silver, R. A., and Sakmann, B. (2002). Synaptic connections between layer 4 spiny neurone- layer 2/3 pyramidal cell pairs in juvenile rat barrel cortex: physiology and anatomy of interlaminar signalling within a cortical column. J. Physiol. 538, 803–822. doi: 10.1113/jphysiol.2001.012959

PubMed Abstract | CrossRef Full Text | Google Scholar

Feldmeyer, D. (2006). Efficacy and connectivity of intracolumnar pairs of layer 2/3 pyramidal cells in the barrel cortex of juvenile rats. J. Physiol. 575, 583–602. doi: 10.1113/jphysiol.2006.105106

PubMed Abstract | CrossRef Full Text | Google Scholar

Gonchar, Y., and Burkhalter, A. (1997). Three distinct families of GABAergic neurons in rat visual cortex. Cereb. Cortex 7, 347–358. doi: 10.1093/cercor/7.4.347

PubMed Abstract | CrossRef Full Text | Google Scholar

Gonchar, Y., Wang, Q., Burkhalter, A. H., Gonchar, Y., Wang, Q., and Burkhalter, A. (2008). Multiple distinct subtypes of GABAergic neurons in mouse visual cortex identified by triple immunostaining. Front. Neuroanat. 1:3. doi: 10.3389/neuro.05.003.2007

PubMed Abstract | CrossRef Full Text | Google Scholar

Gupta, A., Wang, Y., and Markram, H. (2000). Organizing principles for a diversity of GABAergic interneurons and synapses in the neocortex. Science 287, 273–278. doi: 10.1126/science.287.5451.273

PubMed Abstract | CrossRef Full Text | Google Scholar

Kasper, E. M., Larkman, A. U., Lübke, J., and Blakemore, C. (1994). Pyramidal neurons in layer 5 of the rat visual cortex. II. Development of electrophysiological properties. J. Comp. Neurol. 339, 475–494. doi: 10.1002/cne.903390403

PubMed Abstract | CrossRef Full Text | Google Scholar

Kawaguchi, Y., and Kondo, S. (2002). Parvalbumin, somatostatin and cholecystokinin as chemical markers for specific GABAergic interneuron types in the rat frontal cortex. J. Neurocytol. 31, 277–287. doi: 10.1023/A:1024126110356

PubMed Abstract | CrossRef Full Text | Google Scholar

Kawaguchi, Y., and Kubota, Y. (1997). GABAergic cell subtypes and their synaptic connections in rat frontal cortex. Cereb. Cortex 7, 476–486. doi: 10.1093/cercor/7.6.476

PubMed Abstract | CrossRef Full Text | Google Scholar

Kole, M. H. P., Hallermann, S., and Stuart, G. J. (2006). Single Ih channels in pyramidal neuron dendrites: properties, distribution, and impact on action potential output. J. Neurosci. 26, 1677–1687. doi: 10.1523/JNEUROSCI.3664-05.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Korngreen, A., and Sakmann, B. (2000). Voltage-gated K+ channels in layer 5 neocortical pyramidal neurones from young rats: subtypes and gradients. J. Physiol. 525, 621–639. doi: 10.1111/j.1469-7793.2000.00621.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Larkman, A. U. (1991). Dendritic morphology of pyramidal neurones of the visual cortex of the rat: I. Branching patterns. J. Comp. Neurol. 306, 307–319. doi: 10.1002/cne.903060207

PubMed Abstract | CrossRef Full Text | Google Scholar

Larkum, M. E., Nevian, T., Sandler, M., Polsky, A., and Schiller, J. (2009). Synaptic integration in tuft dendrites of layer 5 pyramidal neurons: a new unifying principle. Science 325, 756–760. doi: 10.1126/science.1171958

PubMed Abstract | CrossRef Full Text | Google Scholar

Larkum, M. E., Waters, J., Sakmann, B., and Helmchen, F. (2007). Dendritic spikes in apical dendrites of neocortical layer 2/3 pyramidal neurons. J. Neurosci. 27, 8999–9008. doi: 10.1523/JNEUROSCI.1717-07.2007

PubMed Abstract | CrossRef Full Text | Google Scholar

Ledergerber, D., and Larkum, M. E. (2010). Properties of layer 6 pyramidal neuron apical dendrites. J. Neurosci. 30, 13031–13044. doi: 10.1523/JNEUROSCI.2254-10.2010

PubMed Abstract | CrossRef Full Text | Google Scholar

Lübke, J. (2003). Morphometric analysis of the columnar innervation domain of neurons connecting layer 4 and layer 2/3 of Juvenile rat barrel cortex. Cereb. Cortex 13, 1051–1063. doi: 10.1093/cercor/13.10.1051

PubMed Abstract | CrossRef Full Text | Google Scholar

Markram, H., Lübke, J., Frotscher, M., Roth, A., and Sakmann, B. (1997). Physiology and anatomy of synaptic connections between thick tufted pyramidal neurones in the developing rat neocortex. J. Physiol. 500, 409–440. doi: 10.1113/jphysiol.1997.sp022031

PubMed Abstract | CrossRef Full Text | Google Scholar

Markram, H., Muller, E., Ramaswamy, S., Reimann, M. W., Abdellah, M., Aguado Sanchez, C., et al. (2015). Reconstruction and simulation of neocortical microcircuitry. Cell 163, 456–492. doi: 10.1016/j.cell.2015.09.029

CrossRef Full Text

Markram, H., Toledo-Rodriguez, M., Wang, Y., Gupta, A., Silberberg, G., and Wu, C. (2004). Interneurons of the neocortical inhibitory system. Nat. Rev. Neurosci. 5, 793–807. doi: 10.1038/nrn1519

PubMed Abstract | CrossRef Full Text | Google Scholar

Meyer, H. S., Schwarz, D., Wimmer, V. C., Schmitt, A. C., Kerr, J. N. D., Sakmann, B., et al. (2011). Inhibitory interneurons in a cortical column form hot zones of inhibition in layers 2 and 5A. Proc. Natl. Acad. Sci. U.S.A. 108, 16807–16812. doi: 10.1073/pnas.1113648108

PubMed Abstract | CrossRef Full Text | Google Scholar

Meyer, H. S., Wimmer, V. C., Oberlaender, M., de Kock, C. P. J., Sakmann, B., and Helmstaedter, M. (2010). Number and laminar distribution of neurons in a thalamocortical projection column of rat vibrissal cortex. Cereb. Cortex 20, 2277–2286. doi: 10.1093/cercor/bhq067

PubMed Abstract | CrossRef Full Text | Google Scholar

Narayanan, R. T., Egger, R., Johnson, A. S., Mansvelder, H. D., Sakmann, B., Kock, C. P. J., et al. (2015). Beyond columnar organization: cell type- and target layer-specific principles of horizontal axon projection patterns in rat vibrissal cortex. Cereb. Cortex. doi: 10.1093/cercor/bhv053. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Nevian, T., Larkum, M. E., Polsky, A., and Schiller, J. (2007). Properties of basal dendrites of layer 5 pyramidal neurons: a direct patch-clamp recording study. Nat. Neurosci. 10, 206–214. doi: 10.1038/nn1826

PubMed Abstract | CrossRef Full Text | Google Scholar

O'Connor, D. H., Huber, D., and Svoboda, K. (2009). Reverse engineering the mouse brain. Nature 461, 923–929. doi: 10.1038/nature08539

PubMed Abstract | CrossRef Full Text | Google Scholar

Oberlaender, M., Boudewijns, Z. S. R. M., Kleele, T., Mansvelder, H. D., Sakmann, B., and de Kock, C. P. J. (2011). Three-dimensional axon morphologies of individual layer 5 neurons indicate cell type-specific intracortical pathways for whisker motion and touch. Proc. Natl. Acad. Sci. U.S.A. 108, 4188–4193. doi: 10.1073/pnas.1100647108

PubMed Abstract | CrossRef Full Text | Google Scholar

Oberlaender, M., de Kock, C. P. J., de Bruno, R. M., Ramirez, A., Meyer, H. S., Dercksen, V. J., et al. (2012). Cell type–specific three-dimensional structure of thalamocortical circuits in a column of rat vibrissal cortex. Cereb. Cortex 22, 2375–2391. doi: 10.1093/cercor/bhr317

PubMed Abstract | CrossRef Full Text | Google Scholar

Perin, R., Berger, T. K., and Markram, H. (2011). A synaptic organizing principle for cortical neuronal groups. Proc. Natl. Acad. Sci. U.S.A. 108, 5419–5424. doi: 10.1073/pnas.1016051108

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramaswamy, S., and Markram, H. (2015). Anatomy and physiology of the thick-tufted layer 5 pyramidal neuron. Front. Cell Neurosci. 9:233. doi: 10.3389/fncel.2015.00233

PubMed Abstract | CrossRef Full Text | Google Scholar

Reyes, A., Lujan, R., Rozov, A., Burnashev, N., Somogyi, P., and Sakmann, B. (1998). Target-cell-specific facilitation and depression in neocortical circuits. Nat. Neurosci. 1, 279–285. doi: 10.1038/1092

PubMed Abstract | CrossRef Full Text | Google Scholar

Romand, S., Wang, Y., Toledo-Rodriguez, M., and Markram, H. (2011). Morphological development of thick-tufted layer v pyramidal cells in the rat somatosensory cortex. Front. Neuroanat. 5:5. doi: 10.3389/fnana.2011.00005

PubMed Abstract | CrossRef Full Text | Google Scholar

Santana, R., McGarry, L. M., Bielza, C., Larrañaga, P., and Yuste, R. (2013). Classification of neocortical interneurons using affinity propagation. Front. Neural Circuits 7:185. doi: 10.3389/fncir.2013.00185

PubMed Abstract | CrossRef Full Text | Google Scholar

Silberberg, G., and Markram, H. (2007). Disynaptic inhibition between neocortical pyramidal cells mediated by martinotti cells. Neuron 53, 735–746. doi: 10.1016/j.neuron.2007.02.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Somogyi, P., Tamás, G., Lujan, R., and Buhl, E. H. (1998). Salient features of synaptic organisation in the cerebral cortex. Brain Res. Rev. 26, 113–135. doi: 10.1016/S0165-0173(97)00061-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Stuart, G. J., and Sakmann, B. (1994). Active propagation of somatic action potentials into neocortical pyramidal cell dendrites. Nature 367, 69–72. doi: 10.1038/367069a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomson, A. M., and Deuchars, J. (1997). Synaptic interactions in neocortical local circuits: dual intracellular recordings in vitro. Cereb. Cortex 7, 510–522. doi: 10.1093/cercor/7.6.510

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomson, A. M., Deuchars, J., and West, D. C. (1993). Large, deep layer pyramid-pyramid single axon EPSPs in slices of rat motor cortex display paired pulse and frequency-dependent depression, mediated presynaptically and self-facilitation, mediated postsynaptically. J. Neurophysiol. 70, 2354–2369.

PubMed Abstract | Google Scholar

Toledo-Rodriguez, M., Blumenfeld, B., Wu, C., Luo, J., Attali, B., Goodman, P., et al. (2004). Correlation maps allow neuronal electrical properties to be predicted from single-cell gene expression profiles in rat neocortex. Cereb. Cortex 14, 1310–1327. doi: 10.1093/cercor/bhh092

PubMed Abstract | CrossRef Full Text | Google Scholar

Toledo-Rodriguez, M., Goodman, P., Illic, M., Wu, C., and Markram, H. (2005). Neuropeptide and calcium-binding protein gene expression profiles predict neuronal anatomical type in the juvenile rat. J. Physiol. 567, 401–413. doi: 10.1113/jphysiol.2005.089250

PubMed Abstract | CrossRef Full Text

Tsodyks, M. V., and Markram, H. (1997). The neural code between neocortical pyramidal neurons depends on neurotransmitter release probability. Proc. Natl. Acad. Sci. U.S.A. 94, 5495. doi: 10.1073/pnas.94.2.719

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Gupta, A., Toledo-Rodriguez, M., Wu, C. Z., and Markram, H. (2002). Anatomical, physiological, molecular and circuit properties of nest basket cells in the developing somatosensory cortex. Cereb. Cortex 12, 395–410. doi: 10.1093/cercor/12.4.395

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Markram, H., Goodman, P. H., Berger, T. K., Ma, J., and Goldman-Rakic, P. S. (2006). Heterogeneity in the pyramidal network of the medial prefrontal cortex. Nat. Neurosci. 9, 534–542. doi: 10.1038/nn1670

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Toledo-Rodriguez, M., Gupta, A., Wu, C., Silberberg, G., Luo, J., et al. (2004). Anatomical, physiological and molecular properties of Martinotti cells in the somatosensory cortex of the juvenile rat. J. Physiol. 561, 65–90. doi: 10.1113/jphysiol.2004.073353

PubMed Abstract | CrossRef Full Text | Google Scholar

Wimmer, V. C., Bruno, R. M., de Kock, C. P., Kuner, T., and Sakmann, B. (2010). Dimensions of a projection column and architecture of VPM and POm axons in rat vibrissal cortex. Cereb. Cortex 20, 2265–2276. doi: 10.1093/cercor/bhq068

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, J. J. (2000). Maturation of layer 5 neocortical pyramidal neurons: amplifying salient layer 1 and layer 4 inputs by Ca2+ action potentials in adult rat tuft dendrites. J. Physiol. 526, 571–587. doi: 10.1111/j.1469-7793.2000.00571.x

PubMed Abstract | CrossRef Full Text | Google Scholar