Geotechnical investigation of soil properties in Hatsalatladi village, Botswana; Insights from aeromagnetic, laboratory soil tests and Rayleigh wave dispersion datasets.

Soil tests and Multichannel Analysis of Surface Waves (MASW) data were conducted in Hatsalatladi village, Botswana, to investigate the occurrence of ground fissures within the village and to identify the likely causes of the fissures and their depth extent. The MASW data were collected to gain insights into the variation of shear wave velocity with depth. The dataset shows that the shear wave velocity ranged from 150 m/s - 500 m/s, with Poisson's ratios ranging from 0.02 to 0.25. A low-velocity zone (LVZ) was observed in the upper 5 m of the subsurface with velocities ranging from 200 m/s to 350 m/s. The soil plasticity was measured through the plastic and liquid Atterberg tests. Atterberg limits measurements obtained from the three survey sites show that the plastic index of the soil samples collected from depths of 1 m fall within the 10-20% range. Specifically, the Filled Crack survey site had a plastic index of 16%, while the Abandoned House and Bridge sites had 18.7% and 13.5%, respectively. Soil samples from Filled Crack and Abandoned House site revealed a linear shrinkage of 6.4%, while the Bridge site soil sample had a linear shrinkage of 2.9%. The sieve analysis test results are also presented.

Thalamic control of sensory processing and spindles in a biophysical somatosensory thalamoreticular circuit model of wakefulness and sleep.

Thalamoreticular circuitry plays a key role in arousal, attention, cognition, and sleep spindles, and is linked to several brain disorders. A detailed computational model of mouse somatosensory thalamus and thalamic reticular nucleus has been developed to capture the properties of over 14,000 neurons connected by 6 million synapses. The model recreates the biological connectivity of these neurons, and simulations of the model reproduce multiple experimental findings in different brain states. The model shows that inhibitory rebound produces frequency-selective enhancement of thalamic responses during wakefulness. We find that thalamic interactions are responsible for the characteristic waxing and waning of spindle oscillations. In addition, we find that changes in thalamic excitability control spindle frequency and their incidence. The model is made openly available to provide a new tool for studying the function and dysfunction of the thalamoreticular circuitry in various brain states.

STEPS 4.0: Fast and memory-efficient molecular simulations of neurons at the nanoscale.

Recent advances in computational neuroscience have demonstrated the usefulness and importance of stochastic, spatial reaction-diffusion simulations. However, ever increasing model complexity renders traditional serial solvers, as well as naive parallel implementations, inadequate. This paper introduces a new generation of the STochastic Engine for Pathway Simulation (STEPS) project (http://steps.sourceforge.net/), denominated STEPS 4.0, and its core components which have been designed for improved scalability, performance, and memory efficiency. STEPS 4.0 aims to enable novel scientific studies of macroscopic systems such as whole cells while capturing their nanoscale details. This class of models is out of reach for serial solvers due to the vast quantity of computation in such detailed models, and also out of reach for naive parallel solvers due to the large memory footprint. Based on a distributed mesh solution, we introduce a new parallel stochastic reaction-diffusion solver and a deterministic membrane potential solver in STEPS 4.0. The distributed mesh, together with improved data layout and algorithm designs, significantly reduces the memory footprint of parallel simulations in STEPS 4.0. This enables massively parallel simulations on modern HPC clusters and overcomes the limitations of the previous parallel STEPS implementation. Current and future improvements to the solver are not sustainable without following proper software engineering principles. For this reason, we also give an overview of how the STEPS codebase and the development environment have been updated to follow modern software development practices. We benchmark performance improvement and memory footprint on three published models with different complexities, from a simple spatial stochastic reaction-diffusion model, to a more complex one that is coupled to a deterministic membrane potential solver to simulate the calcium burst activity of a Purkinje neuron. Simulation results of these models suggest that the new solution dramatically reduces the per-core memory consumption by more than a factor of 30, while maintaining similar or better performance and scalability.

The SONATA data format for efficient description of large-scale network models.

Increasing availability of comprehensive experimental datasets and of high-performance computing resources are driving rapid growth in scale, complexity, and biological realism of computational models in neuroscience. To support construction and simulation, as well as sharing of such large-scale models, a broadly applicable, flexible, and high-performance data format is necessary. To address this need, we have developed the Scalable Open Network Architecture TemplAte (SONATA) data format. It is designed for memory and computational efficiency and works across multiple platforms. The format represents neuronal circuits and simulation inputs and outputs via standardized files and provides much flexibility for adding new conventions or extensions. SONATA is used in multiple modeling and visualization tools, and we also provide reference Application Programming Interfaces and model examples to catalyze further adoption. SONATA format is free and open for the community to use and build upon with the goal of enabling efficient model building, sharing, and reproducibility.

Cortical reliability amid noise and chaos.

Typical responses of cortical neurons to identical sensory stimuli appear highly variable. It has thus been proposed that the cortex primarily uses a rate code. However, other studies have argued for spike-time coding under certain conditions. The potential role of spike-time coding is directly limited by the internally generated variability of cortical circuits, which remains largely unexplored. Here, we quantify this internally generated variability using a biophysical model of rat neocortical microcircuitry with biologically realistic noise sources. We find that stochastic neurotransmitter release is a critical component of internally generated variability, causing rapidly diverging, chaotic recurrent network dynamics. Surprisingly, the same nonlinear recurrent network dynamics can transiently overcome the chaos in response to weak feed-forward thalamocortical inputs, and support reliable spike times with millisecond precision. Our model shows that the noisy and chaotic network dynamics of recurrent cortical microcircuitry are compatible with stimulus-evoked, millisecond spike-time reliability, resolving a long-standing debate.

An algorithm to predict the connectome of neural microcircuits.

Experimentally mapping synaptic connections, in terms of the numbers and locations of their synapses and estimating connection probabilities, is still not a tractable task, even for small volumes of tissue. In fact, the six layers of the neocortex contain thousands of unique types of synaptic connections between the many different types of neurons, of which only a handful have been characterized experimentally. Here we present a theoretical framework and a data-driven algorithmic strategy to digitally reconstruct the complete synaptic connectivity between the different types of neurons in a small well-defined volume of tissue-the micro-scale connectome of a neural microcircuit. By enforcing a set of established principles of synaptic connectivity, and leveraging interdependencies between fundamental properties of neural microcircuits to constrain the reconstructed connectivity, the algorithm yields three parameters per connection type that predict the anatomy of all types of biologically viable synaptic connections. The predictions reproduce a spectrum of experimental data on synaptic connectivity not used by the algorithm. We conclude that an algorithmic approach to the connectome can serve as a tool to accelerate experimental mapping, indicating the minimal dataset required to make useful predictions, identifying the datasets required to improve their accuracy, testing the feasibility of experimental measurements, and making it possible to test hypotheses of synaptic connectivity.

Reconstruction and Simulation of Neocortical Microcircuitry.

We present a first-draft digital reconstruction of the microcircuitry of somatosensory cortex of juvenile rat. The reconstruction uses cellular and synaptic organizing principles to algorithmically reconstruct detailed anatomy and physiology from sparse experimental data. An objective anatomical method defines a neocortical volume of 0.29 ± 0.01 mm(3) containing ~31,000 neurons, and patch-clamp studies identify 55 layer-specific morphological and 207 morpho-electrical neuron subtypes. When digitally reconstructed neurons are positioned in the volume and synapse formation is restricted to biological bouton densities and numbers of synapses per connection, their overlapping arbors form ~8 million connections with ~37 million synapses. Simulations reproduce an array of in vitro and in vivo experiments without parameter tuning. Additionally, we find a spectrum of network states with a sharp transition from synchronous to asynchronous activity, modulated by physiological mechanisms. The spectrum of network states, dynamically reconfigured around this transition, supports diverse information processing strategies. PAPERCLIP: VIDEO ABSTRACT.

Intrinsic morphological diversity of thick-tufted layer 5 pyramidal neurons ensures robust and invariant properties of in silico synaptic connections.

The morphology of neocortical pyramidal neurons is not only highly characteristic but also displays an intrinsic diversity that renders each neuron morphologically unique. We investigated the significance of this intrinsic morphological diversity in in silico networks composed of thick-tufted layer 5 (TTL5) pyramidal neurons, by comparing the in silico and in vitro properties of TTL5 synaptic connections. The synaptic locations of in silico connections were determined by placing 3D reconstructed TTL5 neurons randomly in a volume equivalent to that of layer 5 in the juvenile rat somatosensory cortex and using a 'collision-detection' algorithm to identify the incidental loci of axo-dendritic overlap. The activation time of the modelled synapses and their biophysical properties were characterized based on experimental measurements. We found that the anatomical loci of synapses and the physiological properties of the somatically recorded EPSPs closely matched those recorded experimentally without the need for any fine-tuning. Furthermore, perturbations to both the physiological or anatomical parameters of the model did not alter the average physiological properties of the population of modelled synaptic connections. This microcircuit-level robust behaviour was due to the intrinsic diversity of the morphology of pyramidal neurons in the microcircuit. We conclude that synaptic transmission in a network of TTL5 neurons is highly invariant across microcircuits suggesting that intrinsic diversity is a mechanism to ensure the same average synaptic properties in different animals of the same species. Finally, we show that the average physiological properties of the TTL5 microcircuit are surprisingly robust to anatomical and physiological perturbations also partly due to the intrinsic diversity of pyramidal neuron morphology.

A Component-Based Extension Framework for Large-Scale Parallel Simulations in NEURON.

As neuronal simulations approach larger scales with increasing levels of detail, the neurosimulator software represents only a part of a chain of tools ranging from setup, simulation, interaction with virtual environments to analysis and visualizations. Previously published approaches to abstracting simulator engines have not received wide-spread acceptance, which in part may be to the fact that they tried to address the challenge of solving the model specification problem. Here, we present an approach that uses a neurosimulator, in this case NEURON, to describe and instantiate the network model in the simulator's native model language but then replaces the main integration loop with its own. Existing parallel network models are easily adopted to run in the presented framework. The presented approach is thus an extension to NEURON but uses a component-based architecture to allow for replaceable spike exchange components and pluggable components for monitoring, analysis, or control that can run in this framework alongside with the simulation.