Bringing Anatomical Information into Neuronal Network Models.

For constructing neuronal network models computational neuroscientists have access to wide-ranging anatomical data that nevertheless tend to cover only a fraction of the parameters to be determined. Finding and interpreting the most relevant data, estimating missing values, and combining the data and estimates from various sources into a coherent whole is a daunting task. With this chapter we aim to provide guidance to modelers by describing the main types of anatomical data that may be useful for informing neuronal network models. We further discuss aspects of the underlying experimental techniques relevant to the interpretation of the data, list particularly comprehensive data sets, and describe methods for filling in the gaps in the experimental data. Such methods of "predictive connectomics" estimate connectivity where the data are lacking based on statistical relationships with known quantities. Exploiting organizational principles that link the plethora of data in a unifying framework can be useful for informing computational models. Besides overarching principles, we touch upon the most prominent features of brain organization that are likely to influence predicted neuronal network dynamics, with a focus on the mammalian cerebral cortex. Given the still existing need for modelers to navigate a complex data landscape full of holes and stumbling blocks, it is vital that the field of neuroanatomy is moving toward increasingly systematic data collection, representation, and publication.

Unified neural field theory of brain dynamics underlying oscillations in Parkinson's disease and generalized epilepsies.

The mechanisms underlying pathologically synchronized neural oscillations in Parkinson's disease (PD) and generalized epilepsies are explored in parallel via a physiologically-based neural field model of the corticothalamic-basal ganglia (CTBG) system. The basal ganglia (BG) are approximated as a single effective population and their roles in the modulation of oscillatory dynamics of the corticothalamic (CT) system and vice versa are analyzed. In addition to normal EEG rhythms, enhanced activity around 4 Hz and 20 Hz exists in the model, consistent with the characteristic frequencies observed in PD. These rhythms result from resonances in loops formed between the BG and CT populations, analogous to those that underlie epileptic oscillations in a previous CT model, and which are still present in the combined CTBG system. Dopamine depletion is argued to weaken the dampening of these loop resonances in PD, and network connections then explain the significant coherence observed between BG, thalamic, and cortical population activity around 4-8 Hz and 20 Hz. Parallels between the afferent and efferent connection sites of the thalamic reticular nucleus (TRN) and BG predict low dopamine to correspond to a reduced likelihood of tonic-clonic (grand mal) seizures, which agrees with experimental findings. Furthermore, the model predicts an increased likelihood of absence (petit mal) seizure resulting from pathologically low dopamine levels in accordance with experimental observations. Suppression of absence seizure activity is demonstrated when afferent and efferent BG connections to the CT system are strengthened, which is consistent with other CTBG modeling studies. The BG are demonstrated to have a suppressive effect on activity of the CTBG system near tonic-clonic seizure states, which provides insight into the reported efficacy of current treatments in BG circuits. Sleep states of the TRN are also found to suppress pathological PD activity in accordance with observations. Overall, the findings demonstrate strong parallels between coherent oscillations in generalized epilepsies and PD, and provide insights into possible comorbidities.

Relationships between Electroencephalographic Spectral Peaks Across Frequency Bands.

The degree to which electroencephalographic spectral peaks are independent, and the relationships between their frequencies have been debated. A novel fitting method was used to determine peak parameters in the range 2-35 Hz from a large sample of eyes-closed spectra, and their interrelationships were investigated. Findings were compared with a mean-field model of thalamocortical activity, which predicts near-harmonic relationships between peaks. The subject set consisted of 1424 healthy subjects from the Brain Resource International Database. Peaks in the theta range occurred on average near half the alpha peak frequency, while peaks in the beta range tended to occur near twice and three times the alpha peak frequency on an individual-subject basis. Moreover, for the majority of subjects, alpha peak frequencies were significantly positively correlated with frequencies of peaks in the theta and low and high beta ranges. Such a harmonic progression agrees semiquantitatively with theoretical predictions from the mean-field model. These findings indicate a common or analogous source for different rhythms, and help to define appropriate individual frequency bands for peak identification.

Age trends and sex differences of alpha rhythms including split alpha peaks.

OBJECTIVE: To investigate age trends, sex differences, and splitting of alpha peaks of the EEG spectrum in the healthy population. METHODS: An automated multi-site algorithm was used to parametrize the alpha rhythm in 1498 healthy subjects aged 6-86 years. Alpha peaks identified from multiple electrode sites were organized into clusters of similar frequencies whose sex differences and age trends were investigated. RESULTS: Significant age-related trends were observed for frequency, position, and amplitude of dominant alpha peaks. Occipital sites had alpha clusters of higher average frequency, higher power, and greater presence across the scalp. Frequency and power differences were found between the sexes. CONCLUSION: Observed increases in alpha frequency in children and decreases in the elderly were consistent with those from earlier studies. A large fraction of participants (≈ 44%) showed multiple distinct alpha rhythm thus investigations which only examine the alpha frequency with the highest peak power can produce misleading results. The strong dependence of alpha frequency on age and anterior-posterior position indicates use of a fixed alpha frequency band is insufficient to capture the full characteristics of the alpha rhythm. SIGNIFICANCE: This study establishes alpha rhythm parameter ranges (including power and frequency) in the healthy population, and quantifies the variation in alpha frequency across the scalp. The automated characterization enables objective evaluations of alpha band activities over large samples. These findings are potentially useful in testing theories of alpha generation, where splitting of the alpha rhythm has been theoretically predicted to occur in individuals with large differences in axon length between anterior and posterior corticothalamic loops.

Neurophysiological changes with age probed by inverse modeling of EEG spectra.

OBJECTIVE: To investigate age-associated changes in physiologically-based EEG spectral parameters in the healthy population. METHODS: Eyes-closed EEG spectra of 1498 healthy subjects aged 6-86 years were fitted to a mean-field model of thalamocortical dynamics in a cross-sectional study. Parameters were synaptodendritic rates, cortical wave decay rates, connection strengths (gains), axonal delays for thalamocortical loops, and power normalizations. Age trends were approximated using smooth asymptotically linear functions with a single turning point. We also considered sex differences and relationships between model parameters and traditional quantitative EEG measures. RESULTS: The cross-sectional data suggest that changes tend to be most rapid in childhood, generally leveling off at age 15-20 years. Most gains decrease in magnitude with age, as does power normalization. Axonal and dendritic delays decrease in childhood and then increase. Axonal delays and gains show small but significant sex differences. CONCLUSIONS: Mean-field brain modeling allows interpretation of age-associated EEG trends in terms of physiological processes, including the growth and regression of white matter, influencing axonal delays, and the establishment and pruning of synaptic connections, influencing gains. SIGNIFICANCE: This study demonstrates the feasibility of inverse modeling of EEG spectra as a noninvasive method for investigating large-scale corticothalamic dynamics, and provides a basis for future comparisons.

Mean-field modeling of the basal ganglia-thalamocortical system. I Firing rates in healthy and parkinsonian states.

Parkinsonism leads to various electrophysiological changes in the basal ganglia-thalamocortical system (BGTCS), often including elevated discharge rates of the subthalamic nucleus (STN) and the output nuclei, and reduced activity of the globus pallidus external (GPe) segment. These rate changes have been explained qualitatively in terms of the direct/indirect pathway model, involving projections of distinct striatal populations to the output nuclei and GPe. Although these populations partly overlap, evidence suggests dopamine depletion differentially affects cortico-striato-pallidal connection strengths to the two pallidal segments. Dopamine loss may also decrease the striatal signal-to-noise ratio, reducing both corticostriatal coupling and striatal firing thresholds. Additionally, nigrostriatal degeneration may cause secondary changes including weakened lateral inhibition in the GPe, and mesocortical dopamine loss may decrease intracortical excitation and especially inhibition. Here a mean-field model of the BGTCS is presented with structure and parameter estimates closely based on physiology and anatomy. Changes in model rates due to the possible effects of dopamine loss listed above are compared with experiment. Our results suggest that a stronger indirect pathway, possibly combined with a weakened direct pathway, is compatible with empirical evidence. However, altered corticostriatal connection strengths are probably not solely responsible for substantially increased STN activity often found. A lower STN firing threshold, weaker intracortical inhibition, and stronger striato-GPe inhibition help explain the relatively large increase in STN rate. Reduced GPe-GPe inhibition and a lower GPe firing threshold can account for the comparatively small decrease in GPe rate frequently observed. Changes in cortex, GPe, and STN help normalize the cortical rate, also in accord with experiments. The model integrates the basal ganglia into a unified framework along with an existing thalamocortical model that already accounts for a wide range of electrophysiological phenomena. A companion paper discusses the dynamics and oscillations of this combined system.

Mean-field modeling of the basal ganglia-thalamocortical system. II Dynamics of parkinsonian oscillations.

Neuronal correlates of Parkinson's disease (PD) include a shift to lower frequencies in the electroencephalogram (EEG) and enhanced synchronized oscillations at 3-7 and 7-30 Hz in the basal ganglia, thalamus, and cortex. This study describes the dynamics of a recent physiologically based mean-field model of the basal ganglia-thalamocortical system, and shows how it accounts for many key electrophysiological correlates of PD. Its detailed functional connectivity comprises partially segregated direct and indirect pathways through two populations of striatal neurons, a hyperdirect pathway involving a corticosubthalamic projection, thalamostriatal feedback, and local inhibition in striatum and external pallidum (GPe). In a companion paper, realistic steady-state firing rates were obtained for the healthy state, and after dopamine loss modeled by weaker direct and stronger indirect pathways, reduced intrapallidal inhibition, lower firing thresholds of the GPe and subthalamic nucleus (STN), a stronger projection from striatum to GPe, and weaker cortical interactions. Here it is shown that oscillations around 5 and 20 Hz can arise with a strong indirect pathway, which also causes increased synchronization throughout the basal ganglia. Furthermore, increased theta power with progressive nigrostriatal degeneration is correlated with reduced alpha power and peak frequency, in agreement with empirical results. Unlike the hyperdirect pathway, the indirect pathway sustains oscillations with phase relationships that coincide with those found experimentally. Alterations in the responses of basal ganglia to transient stimuli accord with experimental observations. Reduced cortical gains due to both nigrostriatal and mesocortical dopamine loss lead to slower changes in cortical activity and may be related to bradykinesia. Finally, increased EEG power found in some studies may be partly explained by a lower effective GPe firing threshold, reduced GPe-GPe inhibition, and/or weaker intracortical connections in parkinsonian patients. Strict separation of the direct and indirect pathways is not necessary to obtain these results.

Transformation of arbitrary distributions to the normal distribution with application to EEG test-retest reliability.

Many variables in the social, physical, and biosciences, including neuroscience, are non-normally distributed. To improve the statistical properties of such data, or to allow parametric testing, logarithmic or logit transformations are often used. Box-Cox transformations or ad hoc methods are sometimes used for parameters for which no transformation is known to approximate normality. However, these methods do not always give good agreement with the Gaussian. A transformation is discussed that maps probability distributions as closely as possible to the normal distribution, with exact agreement for continuous distributions. To illustrate, the transformation is applied to a theoretical distribution, and to quantitative electroencephalographic (qEEG) measures from repeat recordings of 32 subjects which are highly non-normal. Agreement with the Gaussian was better than using logarithmic, logit, or Box-Cox transformations. Since normal data have previously been shown to have better test-retest reliability than non-normal data under fairly general circumstances, the implications of our transformation for the test-retest reliability of parameters were investigated. Reliability was shown to improve with the transformation, where the improvement was comparable to that using Box-Cox. An advantage of the general transformation is that it does not require laborious optimization over a range of parameters or a case-specific choice of form.