Simple models of quantitative firing phenotypes in hippocampal neurons: Comprehensive coverage of intrinsic diversity.

Patterns of periodic voltage spikes elicited by a neuron help define its dynamical identity. Experimentally recorded spike trains from various neurons show qualitatively distinguishable features such as delayed spiking, spiking with or without frequency adaptation, and intrinsic bursting. Moreover, the input-dependent responses of a neuron not only show different quantitative features, such as higher spike frequency for a stronger input current injection, but can also exhibit qualitatively different responses, such as spiking and bursting under different input conditions, thus forming a complex phenotype of responses. In previous work, the comprehensive knowledge base of hippocampal neuron types Hippocampome.org systematically characterized various spike pattern phenotypes experimentally identified from 120 neuron types/subtypes. In this paper, we present a complete set of simple phenomenological models that quantitatively reproduce the diverse and complex phenotypes of hippocampal neurons. In addition to point-neuron models, we created compact multi-compartment models with up to four compartments, which will allow spatial segregation of synaptic integration in network simulations. Electrotonic compartmentalization observed in our compact multi-compartment models is qualitatively consistent with experimental observations. The models were created using an automated pipeline based on evolutionary algorithms. This work maps 120 neuron types/subtypes in the rodent hippocampus to a low-dimensional model space and adds another dimension to the knowledge accumulated in Hippocampome.org. Computationally efficient representations of intrinsic dynamics, along with other pieces of knowledge available in Hippocampome.org, provide a biologically realistic platform to explore the large-scale interactions of various neuron types at the mesoscopic level.

Itinerant complexity in networks of intrinsically bursting neurons.

Active neurons can be broadly classified by their intrinsic oscillation patterns into two classes characterized by spiking or bursting. Here, we show that networks of identical bursting neurons with inhibitory pulsatory coupling exhibit itinerant dynamics. Using the relative phases of bursts between neurons, we numerically demonstrate that the network exhibits endogenous transitions between multiple modes of transient synchrony. This is true even for bursts consisting of two spikes. In contrast, our simulations reveal that networks of identical singlet-spiking neurons do not exhibit such complexity. These results suggest a role for bursting dynamics in realizing itinerant complexity in neural circuits.

Biophysical modulation and robustness of itinerant complexity in neuronal networks.

Transient synchronization of bursting activity in neuronal networks, which occurs in patterns of metastable itinerant phase relationships between neurons, is a notable feature of network dynamics observed in vivo. However, the mechanisms that contribute to this dynamical complexity in neuronal circuits are not well understood. Local circuits in cortical regions consist of populations of neurons with diverse intrinsic oscillatory features. In this study, we numerically show that the phenomenon of transient synchronization, also referred to as metastability, can emerge in an inhibitory neuronal population when the neurons' intrinsic fast-spiking dynamics are appropriately modulated by slower inputs from an excitatory neuronal population. Using a compact model of a mesoscopic-scale network consisting of excitatory pyramidal and inhibitory fast-spiking neurons, our work demonstrates a relationship between the frequency of pyramidal population oscillations and the features of emergent metastability in the inhibitory population. In addition, we introduce a method to characterize collective transitions in metastable networks. Finally, we discuss potential applications of this study in mechanistically understanding cortical network dynamics.

Increased perivascular space volume in white matter and basal ganglia is associated with cognition in Parkinson's Disease.

Perivascular spaces (PVS), fluid-filled compartments surrounding brain vasculature, are an essential component of the glymphatic system responsible for transport of waste and nutrients. Glymphatic system impairment may underlie cognitive deficits in Parkinson's disease (PD). Studies have focused on the role of basal ganglia PVS with cognition in PD, but the role of white matter PVS is unknown. This study examined the relationship of white matter and basal ganglia PVS with domain-specific and global cognition in individuals with PD. Fifty individuals with PD underwent 3T T1w magnetic resonance imaging (MRI) to determine PVS volume fraction, defined as PVS volume normalized to total regional volume, within (i) centrum semiovale, (ii) prefrontal white matter (medial orbitofrontal, rostral middle frontal, superior frontal), and (iii) basal ganglia. A neuropsychological battery included assessment of global cognitive function (Montreal Cognitive Assessment, and global cognitive composite score), and cognitive-specific domains (executive function, memory, visuospatial function, attention, and language). Higher white matter rostral middle frontal PVS was associated with lower scores in both global cognitive and visuospatial function. In the basal ganglia higher PVS was associated with lower scores for memory with a trend towards lower global cognitive composite score. While previous reports have shown that greater amount of PVS in the basal ganglia is associated with decline in global cognition in PD, our findings suggest that increased white matter PVS volume may also underlie changes in cognition.

Conduction Velocity, G-ratio, and Extracellular Water as Microstructural Characteristics of Autism Spectrum Disorder.

The neuronal differences contributing to the etiology of autism spectrum disorder (ASD) are still not well defined. Previous studies have suggested that myelin and axons are disrupted during development in ASD. By combining structural and diffusion MRI techniques, myelin and axons can be assessed using extracellular water, aggregate g-ratio, and a novel metric termed aggregate conduction velocity, which is related to the capacity of the axon to carry information. In this study, several innovative cellular microstructural methods, as measured from magnetic resonance imaging (MRI), are combined to characterize differences between ASD and typically developing adolescent participants in a large cohort. We first examine the relationship between each metric, including microstructural measurements of axonal and intracellular diffusion and the T1w/T2w ratio. We then demonstrate the sensitivity of these metrics by characterizing differences between ASD and neurotypical participants, finding widespread increases in extracellular water in the cortex and decreases in aggregate g-ratio and aggregate conduction velocity throughout the cortex, subcortex, and white matter skeleton. We finally provide evidence that these microstructural differences are associated with higher scores on the Social Communication Questionnaire (SCQ) a commonly used diagnostic tool to assess ASD. This study is the first to reveal that ASD involves MRI-measurable in vivo differences of myelin and axonal development with implications for neuronal and behavioral function. We also introduce a novel neuroimaging metric, aggregate conduction velocity, that is highly sensitive to these changes. We conclude that ASD may be characterized by otherwise intact structural connectivity but that functional connectivity may be attenuated by network properties affecting neural transmission speed. This effect may explain the putative reliance on local connectivity in contrast to more distal connectivity observed in ASD.

Conduction velocity, G-ratio, and extracellular water as microstructural characteristics of autism spectrum disorder.

The neuronal differences contributing to the etiology of autism spectrum disorder (ASD) are still not well defined. Previous studies have suggested that myelin and axons are disrupted during development in ASD. By combining structural and diffusion MRI techniques, myelin and axons can be assessed using extracellular water, aggregate g-ratio, and a new approach to calculating axonal conduction velocity termed aggregate conduction velocity, which is related to the capacity of the axon to carry information. In this study, several innovative cellular microstructural methods, as measured from magnetic resonance imaging (MRI), are combined to characterize differences between ASD and typically developing adolescent participants in a large cohort. We first examine the relationship between each metric, including microstructural measurements of axonal and intracellular diffusion and the T1w/T2w ratio. We then demonstrate the sensitivity of these metrics by characterizing differences between ASD and neurotypical participants, finding widespread increases in extracellular water in the cortex and decreases in aggregate g-ratio and aggregate conduction velocity throughout the cortex, subcortex, and white matter skeleton. We finally provide evidence that these microstructural differences are associated with higher scores on the Social Communication Questionnaire (SCQ) a commonly used diagnostic tool to assess ASD. This study is the first to reveal that ASD involves MRI-measurable in vivo differences of myelin and axonal development with implications for neuronal and behavioral function. We also introduce a novel formulation for calculating aggregate conduction velocity, that is highly sensitive to these changes. We conclude that ASD may be characterized by otherwise intact structural connectivity but that functional connectivity may be attenuated by network properties affecting neural transmission speed. This effect may explain the putative reliance on local connectivity in contrast to more distal connectivity observed in ASD.

Combinatorial quantification of distinct neural projections from retrograde tracing.

Comprehensive quantification of neuronal architectures underlying anatomical brain connectivity remains challenging. We introduce a method to identify the distinct axonal projection patterns from a source to a set of target regions and the count of neurons with each pattern. For a source region projecting to n targets, there are 2n - 1 theoretically possible projection types, although only a subset of these types typically exists. By injecting uniquely labeled retrograde tracers in k regions (k < n), one can experimentally count the cells expressing different combinations of colors in the source region1,2. Such an experiment can be performed for n choose k combinations. The counts of cells with different color combinations from all experiments provide constraints for a system of equations that include 2n - 1 unknown variables, each corresponding to the count of neurons for a projection pattern. Evolutionary algorithms prove to be effective at solving the resultant system of equations, thus allowing the determination of the counts of neurons with each of the possible projection patterns. Numerical analysis of simulated 4 choose 3 retrograde injection experiments using surrogate data demonstrates reliable and precise count estimates for all projection neuron types. We illustrate the experimental application of this framework by quantifying the projections of mouse primary motor cortex to four prominent targets: the primary and secondary somatosensory and motor cortices.

Combinatorial quantification of distinct neural projections from retrograde tracing.

Comprehensive quantification of neuronal architectures underlying anatomical brain connectivity remains challenging. We introduce a method to identify distinct axonal projection patterns from a source to a set of target regions and the count of neurons with each pattern. A source region projecting to n targets could have 2n-1 theoretically possible projection types, although only a subset of these types typically exists. By injecting uniquely labeled retrograde tracers in k target regions (k < n), one can experimentally count the cells expressing different color combinations in the source region. The neuronal counts for different color combinations from n-choose-k experiments provide constraints for a model that is robustly solvable using evolutionary algorithms. Here, we demonstrate this method's reliability for 4 targets using simulated triple injection experiments. Furthermore, we illustrate the experimental application of this framework by quantifying the projections of male mouse primary motor cortex to the primary and secondary somatosensory and motor cortices.

Magnetic resonance spectroscopy shows associations between neurometabolite levels and perivascular space volume in Parkinson's disease: a pilot and feasibility study.

OBJECTIVE: Higher volume fraction of perivascular space (PVS) has recently been reported in Parkinson's disease (PD) and related disorders. Both elevated PVS and altered levels of neurometabolites, assayed by proton magnetic resonance spectroscopy (MRS), are suspected indicators of neuroinflammation, but no published reports have concurrently examined PVS and MRS neurometabolites. METHODS: In an exploratory pilot study, we acquired multivoxel 3-T MRS using a semi-Localization by Adiabatic SElective Refocusing (sLASER) pulse-sequence (repetition time/echo time = 2810/60 ms, voxels 10 × 10 × 10 mm3) from a 2D slab sampling bilateral frontal white matter (FWM) and anterior middle cingulate cortex (aMCC). PVS maps obtained from high-resolution (0.8 × 0.8 × 0.8 mm3) T1-weighted MRI were co-registered with MRS. In each MRS voxel, PVS volume and neurometabolite levels were measured. RESULTS: Linear regression accounting for age, sex, and BMI found greater PVS volume for higher levels of choline-containing compounds (Cho; P = 0.047) in FWM and lower PVS volume for higher levels of N-acetyl compounds (NAA; P = 0.012) in aMCC. Since (putatively) higher Cho is associated with inflammation while NAA has anti-inflammatory properties, these observations add to evidence that higher PVS load is a sign of inflammation. Additionally, lower Montreal Cognitive Assessment scores were associated with lower NAA in aMCC (P = 0.002), suggesting that local neuronal dysfunction and inflammation contribute to cognitive impairment in PD. CONCLUSION: These exploratory findings indicate that co-analysis of PVS and MRS is feasible and may help elucidate the cellular and metabolic substrates of glymphatic and inflammatory processes in PD.

Physical activity intensity is associated with cognition and functional connectivity in Parkinson's disease.

BACKGROUND: Cognitive impairment is common in Parkinson's disease (PD) and often leads to dementia, with no effective treatment. Aging studies suggest that physical activity (PA) intensity has a positive impact on cognition and enhanced functional connectivity may underlie these benefits. However, less is known in PD. This cross-sectional study examined the relationship between PA intensity, cognitive performance, and resting state functional connectivity in PD and whether PA intensity influences the relationship between functional connectivity and cognitive performance. METHODS: 96 individuals with mild-moderate PD completed a comprehensive neuropsychological battery. Intensity of PA was objectively captured over a seven-day period using a wearable device (ActiGraph). Time spent in light and moderate intensity PA was determined based on standardized actigraphy cut points. Resting-state fMRI was assessed in a subset of 50 individuals to examine brain-wide functional connectivity. RESULTS: Moderate intensity PA (MIPA), but not light PA, was associated with better global cognition, visuospatial function, memory, and executive function. Individuals who met the WHO recommendation of ≥150 min/week of MIPA demonstrated better global cognition, executive function, and visuospatial function. Resting-state functional connectivity associated with MIPA included a combination of brainstem, hippocampus, and regions in the frontal, cingulate, and parietal cortices, which showed higher connectivity across the brain in those achieving the WHO MIPA recommendation. Meeting this recommendation positively moderated the associations between identified functional connectivity and global cognition, visuospatial function, and language. CONCLUSION: Encouraging MIPA, particularly the WHO recommendation of ≥150 min of MIPA/week, may represent an important prescription for PD cognition.

Integrative Models of Brain Structure and Dynamics: Concepts, Challenges, and Methods.

The anatomical architecture of the brain constrains the dynamics of interactions between various regions. On a microscopic scale, neural plasticity regulates the connections between individual neurons. This microstructural adaptation facilitates coordinated dynamics of populations of neurons (mesoscopic scale) and brain regions (macroscopic scale). However, the mechanisms acting on multiple timescales that govern the reciprocal relationship between neural network structure and its intrinsic dynamics are not well understood. Studies empirically investigating such relationships on the whole-brain level rely on macroscopic measurements of structural and functional connectivity estimated from various neuroimaging modalities such as Diffusion-weighted Magnetic Resonance Imaging (dMRI), Electroencephalography (EEG), Magnetoencephalography (MEG), and functional Magnetic Resonance Imaging (fMRI). dMRI measures the anisotropy of water diffusion along axonal fibers, from which structural connections are estimated. EEG and MEG signals measure electrical activity and magnetic fields induced by the electrical activity, respectively, from various brain regions with a high temporal resolution (but limited spatial coverage), whereas fMRI measures regional activations indirectly via blood oxygen level-dependent (BOLD) signals with a high spatial resolution (but limited temporal resolution). There are several studies in the neuroimaging literature reporting statistical associations between macroscopic structural and functional connectivity. On the other hand, models of large-scale oscillatory dynamics conditioned on network structure (such as the one estimated from dMRI connectivity) provide a platform to probe into the structure-dynamics relationship at the mesoscopic level. Such investigations promise to uncover the theoretical underpinnings of the interplay between network structure and dynamics and could be complementary to the macroscopic level inquiries. In this article, we review theoretical and empirical studies that attempt to elucidate the coupling between brain structure and dynamics. Special attention is given to various clinically relevant dimensions of brain connectivity such as the topological features and neural synchronization, and their applicability for a given modality, spatial or temporal scale of analysis is discussed. Our review provides a summary of the progress made along this line of research and identifies challenges and promising future directions for multi-modal neuroimaging analyses.

Quantitative firing pattern phenotyping of hippocampal neuron types.

Systematically organizing the anatomical, molecular, and physiological properties of cortical neurons is important for understanding their computational functions. Hippocampome.org defines 122 neuron types in the rodent hippocampal formation based on their somatic, axonal, and dendritic locations, putative excitatory/inhibitory outputs, molecular marker expression, and biophysical properties. We augmented the electrophysiological data of this knowledge base by collecting, quantifying, and analyzing the firing responses to depolarizing current injections for every hippocampal neuron type from published experiments. We designed and implemented objective protocols to classify firing patterns based on 5 transients (delay, adapting spiking, rapidly adapting spiking, transient stuttering, and transient slow-wave bursting) and 4 steady states (non-adapting spiking, persistent stuttering, persistent slow-wave bursting, and silence). This automated approach revealed 9 unique (plus one spurious) families of firing pattern phenotypes while distinguishing potential new neuronal subtypes. Novel statistical associations emerged between firing responses and other electrophysiological properties, morphological features, and molecular marker expression. The firing pattern parameters, experimental conditions, spike times, references to the original empirical evidences, and analysis scripts are released open-source through Hippocampome.org for all neuron types, greatly enhancing the existing search and browse capabilities. This information, collated online in human- and machine-accessible form, will help design and interpret both experiments and model simulations.

Evolving Simple Models of Diverse Intrinsic Dynamics in Hippocampal Neuron Types.

The diversity of intrinsic dynamics observed in neurons may enhance the computations implemented in the circuit by enriching network-level emergent properties such as synchronization and phase locking. Large-scale spiking network models of entire brain regions offer a platform to test theories of neural computation and cognitive function, providing useful insights on information processing in the nervous system. However, a systematic in-depth investigation requires network simulations to capture the biological intrinsic diversity of individual neurons at a sufficient level of accuracy. The computationally efficient Izhikevich model can reproduce a wide range of neuronal behaviors qualitatively. Previous studies using optimization techniques, however, were less successful in quantitatively matching experimentally recorded voltage traces. In this article, we present an automated pipeline based on evolutionary algorithms to quantitatively reproduce features of various classes of neuronal spike patterns using the Izhikevich model. Employing experimental data from Hippocampome.org, a comprehensive knowledgebase of neuron types in the rodent hippocampus, we demonstrate that our approach reliably fit Izhikevich models to nine distinct classes of experimentally recorded spike patterns, including delayed spiking, spiking with adaptation, stuttering, and bursting. Importantly, by leveraging the parameter-exploration capabilities of evolutionary algorithms, and by representing qualitative spike pattern class definitions in the error landscape, our approach creates several suitable models for each neuron type, exhibiting appropriate feature variabilities among neurons. Moreover, we demonstrate the flexibility of our methodology by creating multi-compartment Izhikevich models for each neuron type in addition to single-point versions. Although the results presented here focus on hippocampal neuron types, the same strategy is broadly applicable to any neural systems.