PyRates-A code-generation tool for modeling dynamical systems in biology and beyond.

The mathematical study of real-world dynamical systems relies on models composed of differential equations. Numerical methods for solving and analyzing differential equation systems are essential when complex biological problems have to be studied, such as the spreading of a virus, the evolution of competing species in an ecosystem, or the dynamics of neurons in the brain. Here we present PyRates, a Python-based software for modeling and analyzing differential equation systems via numerical methods. PyRates is specifically designed to account for the inherent complexity of biological systems. It provides a new language for defining models that mirrors the modular organization of real-world dynamical systems and thus simplifies the implementation of complex networks of interacting dynamic entities. Furthermore, PyRates provides extensive support for the various forms of interaction delays that can be observed in biological systems. The core of PyRates is a versatile code-generation system that translates user-defined models into "backend" implementations in various languages, including Python, Fortran, Matlab, and Julia. This allows users to apply a wide range of analysis methods for dynamical systems, eliminating the need for manual translation between code bases. PyRates may also be used as a model definition interface for the creation of custom dynamical systems tools. To demonstrate this, we developed two extensions of PyRates for common analyses of dynamic models of biological systems: PyCoBi for bifurcation analysis and RectiPy for parameter fitting. We demonstrate in a series of example models how PyRates can be used in combination with PyCoBi and RectiPy for model analysis and fitting. Together, these tools offer a versatile framework for applying computational modeling and numerical analysis methods to dynamical systems in biology and beyond.

Laminar neural dynamics of auditory evoked responses: Computational modeling of local field potentials in auditory cortex of non-human primates.

Evoked neural responses to sensory stimuli have been extensively investigated in humans and animal models both to enhance our understanding of brain function and to aid in clinical diagnosis of neurological and neuropsychiatric conditions. Recording and imaging techniques such as electroencephalography (EEG), magnetoencephalography (MEG), local field potentials (LFPs), and calcium imaging provide complementary information about different aspects of brain activity at different spatial and temporal scales. Modeling and simulations provide a way to integrate these different types of information to clarify underlying neural mechanisms. In this study, we aimed to shed light on the neural dynamics underlying auditory evoked responses by fitting a rate-based model to LFPs recorded via multi-contact electrodes which simultaneously sampled neural activity across cortical laminae. Recordings included neural population responses to best-frequency (BF) and non-BF tones at four representative sites in primary auditory cortex (A1) of awake monkeys. The model considered major neural populations of excitatory, parvalbumin-expressing (PV), and somatostatin-expressing (SOM) neurons across layers 2/3, 4, and 5/6. Unknown parameters, including the connection strength between the populations, were fitted to the data. Our results revealed similar population dynamics, fitted model parameters, predicted equivalent current dipoles (ECD), tuning curves, and lateral inhibition profiles across recording sites and animals, in spite of quite different extracellular current distributions. We found that PV firing rates were higher in BF than in non-BF responses, mainly due to different strengths of tonotopic thalamic input, whereas SOM firing rates were higher in non-BF than in BF responses due to lateral inhibition. In conclusion, we demonstrate the feasibility of the model-fitting approach in identifying the contributions of cell-type specific population activity to stimulus-evoked LFPs across cortical laminae, providing a foundation for further investigations into the dynamics of neural circuits underlying cortical sensory processing.

White matter brain structure predicts language performance and learning success.

Individual differences in the ability to process language have long been discussed. Much of the neural basis of these, however, is yet unknown. Here we investigated the relationship between long-range white matter connectivity of the brain, as revealed by diffusion tractography, and the ability to process syntactically complex sentences in the participants' native language as well as the improvement thereof by multiday training. We identified specific network motifs by singular value decomposition that indeed related white matter structural connectivity to individual language processing performance. First, for two such motifs, one in the left and one in the right hemisphere, their individual prevalence significantly predicted the individual language performance, suggesting an anatomical predisposition for the individual ability to process syntactically complex sentences. Both motifs comprise a number of cortical regions, but seem to be dominated by areas known for the involvement in working memory rather than the classical language network itself. Second, we identified another left hemispheric network motif, whose change of prevalence over the training period significantly correlated with the individual change in performance, thus reflecting training induced white matter plasticity. This motif comprises diverse cortical areas including regions known for their involvement in language processing, working memory and motor functions. The present findings suggest that individual differences in language processing and learning can be explained, in part, by individual differences in the brain's white matter structure. Brain structure may be a crucial factor to be considered when discussing variations in human cognitive performance, more generally.

Resting-State fMRI Functional Connectivity Strength Predicts Local Activity Change in the Dorsal Cingulate Cortex: A Multi-Target Focused rTMS Study.

Previous resting state functional magnetic resonance imaging (RS-fMRI) studies suggested that repetitive transcranial magnetic stimulation (rTMS) can modulate local activity in distant areas via functional connectivity (FC). A brain region has more than one connection with the superficial cortical areas. The current study proposed a multi-target focused rTMS protocol for indirectly stimulating a deep region, and to investigate 1) whether FC strength between stimulation targets (right middle frontal gyrus [rMFG] and right inferior parietal lobule [rIPL]) and effective region (dorsal anterior cingulate cortex [dACC]) can predict local activity changes of dACC and 2) whether multiple stimulation targets can focus on the dACC via FC. A total of 24 healthy participants received rTMS with two stimulation targets, both showing strong FC with the dACC. There were four rTMS conditions (>1 week apart, 10 Hz, 1800 pulses for each): rMFG-target, rIPL-target, Double-targets (900 pulses for each target), and Sham. The results failed to validate the multi-target focused rTMS hypothesis. But rMFG-target significantly decreased the local activity in the dACC. In addition, stronger dACC-rMFG FC was associated with a greater local activity change in the dACC. Future studies should use stronger FC to focus stimulation effects on the deep region.

On the Role of Arkypallidal and Prototypical Neurons for Phase Transitions in the External Pallidum.

The external pallidum (globus pallidus pars externa [GPe]) plays a central role for basal ganglia functions and dynamics and, consequently, has been included in most computational studies of the basal ganglia. These studies considered the GPe as a homogeneous neural population. However, experimental studies have shown that the GPe contains at least two distinct cell types (prototypical and arkypallidal cells). In this work, we provide in silico insight into how pallidal heterogeneity modulates dynamic regimes inside the GPe and how they affect the GPe response to oscillatory input. We derive a mean-field model of the GPe system from a microscopic spiking neural network of recurrently coupled prototypical and arkypallidal neurons. Using bifurcation analysis, we examine the influence of dopamine-dependent changes of intrapallidal connectivity on the GPe dynamics. We find that increased self-inhibition of prototypical cells can induce oscillations, whereas increased inhibition of prototypical cells by arkypallidal cells leads to the emergence of a bistable regime. Furthermore, we show that oscillatory input to the GPe, arriving from striatum, leads to characteristic patterns of cross-frequency coupling observed at the GPe. Based on these findings, we propose two different hypotheses of how dopamine depletion at the GPe may lead to phase-amplitude coupling between the parkinsonian beta rhythm and a GPe-intrinsic γ rhythm. Finally, we show that these findings generalize to realistic spiking neural networks of sparsely coupled Type I excitable GPe neurons.SIGNIFICANCE STATEMENT Our work provides (1) insight into the theoretical implications of a dichotomous globus pallidus pars externa (GPe) organization, and (2) an exact mean-field model that allows for future investigations of the relationship between GPe spiking activity and local field potential fluctuations. We identify the major phase transitions that the GPe can undergo when subject to static or periodic input and link these phase transitions to the emergence of synchronized oscillations and cross-frequency coupling in the basal ganglia. Because of the close links between our model and experimental findings on the structure and dynamics of prototypical and arkypallidal cells, our results can be used to guide both experimental and computational studies on the role of the GPe for basal ganglia dynamics in health and disease.

Cross-frequency phase-amplitude coupling in repetitive movements in patients with Parkinson's disease.

Bradykinesia is a cardinal motor symptom in Parkinson's disease (PD), the pathophysiology of which is not fully understood. We analyzed the role of cross-frequency coupling of oscillatory cortical activity in motor impairment in patients with PD and healthy controls. High-density EEG signals were recorded during various motor activities and at rest. Patients performed a repetitive finger-pressing task normally, but were slower than controls during tapping. Phase-amplitude coupling (PAC) between ß (13-30 Hz) and broadband γ (50-150 Hz) was computed from individual EEG source signals in the premotor, primary motor, and primary somatosensory cortices, and the primary somatosensory complex. In all four regions, averaging the entire movement period resulted in higher PAC in patients than in controls for the resting condition and the pressing task (similar performance between groups). However, this was not the case for the tapping tasks where patients performed slower. This suggests the strength of state-related ß-γ PAC does not determine Parkinsonian bradykinesia. Examination of the dynamics of oscillatory EEG signals during motor transitions revealed a distinctive motif of PAC rise and decay around press onset. This pattern was also present at press offset and slow tapping onset, linking such idiosyncratic PAC changes to transitions between different movement states. The transition-related PAC modulation in patients was similar to controls in the pressing task but flattened during slow tapping, which related to normal and abnormal performance, respectively. These findings suggest that the dysfunctional evolution of neuronal population dynamics during movement execution is an important component of the pathophysiology of Parkinsonian bradykinesia.NEW & NOTEWORTHY Our findings using noninvasive EEG recordings provide evidence that PAC dynamics might play a role in the physiological cortical control of movement execution and may encode transitions between movement states. Results in patients with Parkinson's disease suggest that bradykinesia is related to a deficit of the dynamic regulation of PAC during movement execution rather than its absolute strength. Our findings may contribute to the development of a new concept of the pathophysiology of bradykinesia.

Ephaptic coupling in white matter fibre bundles modulates axonal transmission delays.

Axonal connections are widely regarded as faithful transmitters of neuronal signals with fixed delays. The reasoning behind this is that extracellular potentials caused by spikes travelling along axons are too small to have an effect on other axons. Here we devise a computational framework that allows us to study the effect of extracellular potentials generated by spike volleys in axonal fibre bundles on axonal transmission delays. We demonstrate that, although the extracellular potentials generated by single spikes are of the order of microvolts, the collective extracellular potential generated by spike volleys can reach several millivolts. As a consequence, the resulting depolarisation of the axonal membranes increases the velocity of spikes, and therefore reduces axonal delays between brain areas. Driving a neural mass model with such spike volleys, we further demonstrate that only ephaptic coupling can explain the reduction of stimulus latencies with increased stimulus intensities, as observed in many psychological experiments.

Spatiotemporal features of ß-γ phase-amplitude coupling in Parkinson's disease derived from scalp EEG.

Abnormal phase-amplitude coupling between ß and broadband-γ activities has been identified in recordings from the cortex or scalp of patients with Parkinson's disease. While enhanced phase-amplitude coupling has been proposed as a biomarker of Parkinson's disease, the neuronal mechanisms underlying the abnormal coupling and its relationship to motor impairments in Parkinson's disease remain unclear. To address these issues, we performed an in-depth analysis of high-density EEG recordings at rest in 19 patients with Parkinson's disease and 20 age- and sex-matched healthy control subjects. EEG signals were projected onto the individual cortical surfaces using source reconstruction techniques and separated into spatiotemporal components using independent component analysis. Compared to healthy controls, phase-amplitude coupling of Parkinson's disease patients was enhanced in dorsolateral prefrontal cortex, premotor cortex, primary motor cortex and somatosensory cortex, the difference being statistically significant in the hemisphere contralateral to the clinically more affected side. ß and γ signals involved in generating abnormal phase-amplitude coupling were not strictly phase-phase coupled, ruling out that phase-amplitude coupling merely reflects the abnormal activity of a single oscillator in a recurrent network. We found important differences for couplings between the ß and γ signals from identical components as opposed to those from different components (originating from distinct spatial locations). While both couplings were abnormally enhanced in patients, only the latter were correlated with clinical motor severity as indexed by part III of the Movement Disorder Society Unified Parkinson's Disease Rating Scale. Correlations with parkinsonian motor symptoms of such inter-component couplings were found in premotor, primary motor and somatosensory cortex, but not in dorsolateral prefrontal cortex, suggesting motor domain specificity. The topography of phase-amplitude coupling demonstrated profound differences in patients compared to controls. These findings suggest, first, that enhanced phase-amplitude coupling in Parkinson's disease patients originates from the coupling between distinct neural networks in several brain regions involved in motor control. Because these regions included the somatosensory cortex, abnormal phase-amplitude coupling is not exclusively tied to the hyperdirect tract connecting cortical regions monosynaptically with the subthalamic nucleus. Second, only the coupling between ß and γ signals from different components appears to have pathophysiological significance, suggesting that therapeutic approaches breaking the abnormal lateral coupling between neuronal circuits may be more promising than targeting phase-amplitude coupling per se.

Efficient high-resolution TMS mapping of the human motor cortex by nonlinear regression.

Transcranial magnetic stimulation (TMS) is a powerful tool to investigate causal structure-function relationships in the human brain. However, a precise delineation of the effectively stimulated neuronal populations is notoriously impeded by the widespread and complex distribution of the induced electric field. Here, we propose a method that allows rapid and feasible cortical localization at the individual subject level. The functional relationship between electric field and behavioral effect is quantified by combining experimental data with numerically modeled fields to identify the cortical origin of the modulated effect. Motor evoked potentials (MEPs) from three finger muscles were recorded for a set of random stimulations around the primary motor area. All induced electric fields were nonlinearly regressed against the elicited MEPs to identify their cortical origin. We could distinguish cortical muscle representation with high spatial resolution and localized them primarily on the crowns and rims of the precentral gyrus. A post-hoc analysis revealed exponential convergence of the method with the number of stimulations, yielding a minimum of about 180 random stimulations to obtain stable results. Establishing a functional link between the modulated effect and the underlying mode of action, the induced electric field, is a fundamental step to fully exploit the potential of TMS. In contrast to previous approaches, the presented protocol is particularly easy to implement, fast to apply, and very robust due to the random coil positioning and therefore is suitable for practical and clinical applications.

Functional brain plasticity during L1 training on complex sentences: Changes in gamma-band oscillatory activity.

The adult human brain remains plastic even after puberty. However, whether first language (L1) training in adults can alter the language network is yet largely unknown. Thus, we conducted a longitudinal training experiment on syntactically complex German sentence comprehension. Sentence complexity was varied by the depth of the center embedded relative clauses (i.e., single or double embedded). Comprehension was tested after each sentence with a question on the thematic role assignment. Thirty adult, native German speakers were recruited for 4 days of training. Magnetoencephalography (MEG) data were recorded and subjected to spectral power analysis covering the classical frequency bands (i.e., theta, alpha, beta, low gamma, and gamma). Normalized spectral power, time-locked to the final closure of the relative clause, was subjected to a two-factor analysis ("sentence complexity" and "training days"). Results showed that for the more complex sentences, the interaction of sentence complexity and training days was observed in Brodmann area 44 (BA 44) as a decrease of gamma power with training. Moreover, in the gamma band (55-95 Hz) functional connectivity between BA 44 and other brain regions such as the inferior frontal sulcus and the inferior parietal cortex were correlated with behavioral performance increase due to training. These results show that even for native speakers, complex L1 sentence training improves language performance and alters neural activities of the left hemispheric language network. Training strengthens the use of the dorsal processing stream with working-memory-related brain regions for syntactically complex sentences, thereby demonstrating the brain's functional plasticity for L1 training.

A novel approach to localize cortical TMS effects.

Despite the widespread use of transcranial magnetic stimulation (TMS), the precise cortical locations underlying the resulting physiological and behavioral effects are still only coarsely known. To date, mapping strategies have relied on projection approaches (often termed "center of gravity" approaches) or maximum electric field value evaluation, and therefore localize the stimulated cortical site only approximately and indirectly. Focusing on the motor cortex, we present and validate a novel method to reliably determine the effectively stimulated cortical site at the individual subject level. The approach combines measurements of motor evoked potentials (MEPs) at different coil positions and orientations with numerical modeling of induced electric fields. We identify sharply bounded cortical areas, around the gyral crowns and rims of the motor hand area, as the origin of MEPs and show that the magnitude of the tangential component and the overall magnitude of the electric field are most relevant for the observed effect. To validate our approach, we identified the coil location and orientation that produces the maximal electric field at the predicted stimulation site, and then experimentally show that this location produces MEPs more efficiently than other tested locations/orientations. Moreover, we used extensive uncertainty and sensitivity analyses to verify the robustness of the method and identify the most critical model parameters. Our generic approach improves the localization of the cortical area effectively stimulated by TMS and may be transferred to other modalities such as language mapping.

A Mean-Field Description of Bursting Dynamics in Spiking Neural Networks with Short-Term Adaptation.

Bursting plays an important role in neural communication. At the population level, macroscopic bursting has been identified in populations of neurons that do not express intrinsic bursting mechanisms. For the analysis of phase transitions between bursting and non-bursting states, mean-field descriptions of macroscopic bursting behavior are a valuable tool. In this article, we derive mean-field descriptions of populations of spiking neurons and examine whether states of collective bursting behavior can arise from short-term adaptation mechanisms. Specifically, we consider synaptic depression and spike-frequency adaptation in networks of quadratic integrate-and-fire neurons. Analyzing the mean-field model via bifurcation analysis, we find that bursting behavior emerges for both types of short-term adaptation. This bursting behavior can coexist with steady-state behavior, providing a bistable regime that allows for transient switches between synchronized and nonsynchronized states of population dynamics. For all of these findings, we demonstrate a close correspondence between the spiking neural network and the mean-field model. Although the mean-field model has been derived under the assumptions of an infinite population size and all-to-all coupling inside the population, we show that this correspondence holds even for small, sparsely coupled networks. In summary, we provide mechanistic descriptions of phase transitions between bursting and steady-state population dynamics, which play important roles in both healthy neural communication and neurological disorders.

Action potential propagation and synchronisation in myelinated axons.

With the advent of advanced MRI techniques it has become possible to study axonal white matter non-invasively and in great detail. Measuring the various parameters of the long-range connections of the brain opens up the possibility to build and refine detailed models of large-scale neuronal activity. One particular challenge is to find a mathematical description of action potential propagation that is sufficiently simple, yet still biologically plausible to model signal transmission across entire axonal fibre bundles. We develop a mathematical framework in which we replace the Hodgkin-Huxley dynamics by a spike-diffuse-spike model with passive sub-threshold dynamics and explicit, threshold-activated ion channel currents. This allows us to study in detail the influence of the various model parameters on the action potential velocity and on the entrainment of action potentials between ephaptically coupled fibres without having to recur to numerical simulations. Specifically, we recover known results regarding the influence of axon diameter, node of Ranvier length and internode length on the velocity of action potentials. Additionally, we find that the velocity depends more strongly on the thickness of the myelin sheath than was suggested by previous theoretical studies. We further explain the slowing down and synchronisation of action potentials in ephaptically coupled fibres by their dynamic interaction. In summary, this study presents a solution to incorporate detailed axonal parameters into a whole-brain modelling framework.

A principled approach to conductivity uncertainty analysis in electric field calculations.

Uncertainty surrounding ohmic tissue conductivity impedes accurate calculation of the electric fields generated by non-invasive brain stimulation. We present an efficient and generic technique for uncertainty and sensitivity analyses, which quantifies the reliability of field estimates and identifies the most influential parameters. For this purpose, we employ a non-intrusive generalized polynomial chaos expansion to compactly approximate the multidimensional dependency of the field on the conductivities. We demonstrate that the proposed pipeline yields detailed insight into the uncertainty of field estimates for transcranial magnetic stimulation (TMS) and transcranial direct current stimulation (tDCS), identifies the most relevant tissue conductivities, and highlights characteristic differences between stimulation methods. Specifically, we test the influence of conductivity variations on (i) the magnitude of the electric field generated at each gray matter location, (ii) its normal component relative to the cortical sheet, (iii) its overall magnitude (indexed by the 98th percentile), and (iv) its overall spatial distribution. We show that TMS fields are generally less affected by conductivity variations than tDCS fields. For both TMS and tDCS, conductivity uncertainty causes much higher uncertainty in the magnitude as compared to the direction and overall spatial distribution of the electric field. Whereas the TMS fields were predominantly influenced by gray and white matter conductivity, the tDCS fields were additionally dependent on skull and scalp conductivities. Comprehensive uncertainty analyses of complex systems achieved by the proposed technique are not possible with classical methods, such as Monte Carlo sampling, without extreme computational effort. In addition, our method has the advantages of directly yielding interpretable and intuitive output metrics and of being easily adaptable to new problems.

Emergence of cognitive priming and structure building from the hierarchical interaction of canonical microcircuit models.

The concept of connectionism states that higher cognitive functions emerge from the interaction of many simple elements. Accordingly, research on canonical microcircuits conceptualizes findings on fundamental neuroanatomical circuits as well as recurrent organizational principles of the cerebral cortex and examines the link between architectures and their associated functionality. In this study, we establish minimal canonical microcircuit models as elements of hierarchical processing networks. Based on a combination of descriptive time simulations and explanatory state-space mappings, we show that minimal canonical microcircuits effectively segregate feedforward and feedback information flows and that feedback information conditions basic processing operations in minimal canonical microcircuits. Further, we derive and examine two prototypical meta-circuits of cooperating minimal canonical microcircuits for the neurocognitive problems of priming and structure building. Through the application of these findings to a language network of syntax parsing, this study embodies neurocognitive research on hierarchical communication in light of canonical microcircuits, cell assembly theory, and predictive coding.

A generic deviance detection principle for cortical On/Off responses, omission response, and mismatch negativity.

Neural responses to sudden changes can be observed in many parts of the sensory pathways at different organizational levels. For example, deviants that violate regularity at various levels of abstraction can be observed as simple On/Off responses of individual neurons or as cumulative responses of neural populations. The cortical deviance-related responses supporting different functionalities (e.g., gap detection, chunking, etc.) seem unlikely to arise from different function-specific neural circuits, given the relatively uniform and self-similar wiring patterns across cortical areas and spatial scales. Additionally, reciprocal wiring patterns (with heterogeneous combinations of excitatory and inhibitory connections) in the cortex naturally speak in favor of a generic deviance detection principle. Based on this concept, we propose a network model consisting of reciprocally coupled neural masses as a blueprint of a universal change detector. Simulation examples reproduce properties of cortical deviance-related responses including the On/Off responses, the omitted-stimulus response (OSR), and the mismatch negativity (MMN). We propose that the emergence of change detectors relies on the involvement of disinhibition. An analysis of network connection settings further suggests a supportive effect of synaptic adaptation and a destructive effect of N-methyl-D-aspartate receptor (NMDA-r) antagonists on change detection. We conclude that the nature of cortical reciprocal wiring gives rise to a whole range of local change detectors supporting the notion of a generic deviance detection principle. Several testable predictions are provided based on the network model. Notably, we predict that the NMDA-r antagonists would generally dampen the cortical Off response, the cortical OSR, and the MMN.

Cortical parcellation based on structural connectivity: A case for generative models.

One of the major challenges in systems neuroscience is to identify brain networks and unravel their significance for brain function -this has led to the concept of the 'connectome'. Connectomes are currently extensively studied in large-scale international efforts at multiple scales, and follow different definitions with respect to their connections as well as their elements. Perhaps the most promising avenue for defining the elements of connectomes originates from the notion that individual brain areas maintain distinct (long-range) connection profiles. These connectivity patterns determine the areas' functional properties and also allow for their anatomical delineation and mapping. This rationale has motivated the concept of connectivity-based cortex parcellation. In the past ten years, non-invasive mapping of human brain connectivity has led to immense advances in the development of parcellation techniques and their applications. Unfortunately, many of these approaches primarily aim for confirmation of well-known, existing architectonic maps and, to that end, unsuitably incorporate prior knowledge and frequently build on circular argumentation. Often, current approaches also tend to disregard the specific apertures of connectivity measurements, as well as the anatomical specificities of cortical areas, such as spatial compactness, regional heterogeneity, inter-subject variability, the multi-scaling nature of connectivity information, and potential hierarchical organisation. From a methodological perspective, however, a useful framework that regards all of these aspects in an unbiased way is technically demanding. In this commentary, we first outline the concept of connectivity-based cortex parcellation and discuss its prospects and limitations in particular with respect to structural connectivity. To improve reliability and efficiency, we then strongly advocate for connectivity-based cortex parcellation as a modelling approach; that is, an approximation of the data based on (model) parameter inference. As such, a parcellation algorithm can be formally tested for robustness -the precision of its predictions can be quantified and statistics about potential generalization of the results can be derived. Such a framework also allows the question of model constraints to be reformulated in terms of hypothesis testing through model selection and offers a formative way to integrate anatomical knowledge in terms of prior distributions.

Resting state functional connectivity of the ventral auditory pathway in musicians with absolute pitch.

Absolute pitch (AP) is the ability to recognize pitch chroma of tonal sound without external references, providing a unique model of the human auditory system (Zatorre: Nat Neurosci 6 () 692-695). In a previous study (Kim and Knösche: Hum Brain Mapp () 3486-3501), we identified enhanced intracortical myelination in the right planum polare (PP) in musicians with AP, which could be a potential site for perceptional processing of pitch chroma information. We speculated that this area, which initiates the ventral auditory pathway, might be crucially involved in the perceptual stage of the AP process in the context of the "dual pathway hypothesis" that suggests the role of the ventral pathway in processing nonspatial information related to the identity of an auditory object (Rauschecker: Eur J Neurosci 41 () 579-585). To test our conjecture on the ventral pathway, we investigated resting state functional connectivity (RSFC) using functional magnetic resonance imaging (fMRI) from musicians with varying degrees of AP. Should our hypothesis be correct, RSFC via the ventral pathway is expected to be stronger in musicians with AP, whereas such group effect is not predicted in the RSFC via the dorsal pathway. In the current data, we found greater RSFC between the right PP and bilateral anteroventral auditory cortices in musicians with AP. In contrast, we did not find any group difference in the RSFC of the planum temporale (PT) between musicians with and without AP. We believe that these findings support our conjecture on the critical role of the ventral pathway in AP recognition. Hum Brain Mapp 38:3899-3916, 2017. © 2017 Wiley Periodicals, Inc.

Intracortical myelination in musicians with absolute pitch: Quantitative morphometry using 7-T MRI.

Absolute pitch (AP) is known as the ability to recognize and label the pitch chroma of a given tone without external reference. Known brain structures and functions related to AP are mainly of macroscopic aspects. To shed light on the underlying neural mechanism of AP, we investigated the intracortical myeloarchitecture in musicians with and without AP using the quantitative mapping of the longitudinal relaxation rates with ultra-high-field magnetic resonance imaging at 7 T. We found greater intracortical myelination for AP musicians in the anterior region of the supratemporal plane, particularly the medial region of the right planum polare (PP). In the same region of the right PP, we also found a positive correlation with a behavioral index of AP performance. In addition, we found a positive correlation with a frequency discrimination threshold in the anterolateral Heschl's gyrus in the right hemisphere, demonstrating distinctive neural processes of absolute recognition and relative discrimination of pitch. Regarding possible effects of local myelination in the cortex and the known importance of the anterior superior temporal gyrus/sulcus for the identification of auditory objects, we argue that pitch chroma may be processed as an identifiable object property in AP musicians. Hum Brain Mapp 37:3486-3501, 2016. © 2016 Wiley Periodicals, Inc.

MEG source reconstruction based on identification of directed source interactions on whole-brain anatomical networks.

We present an MEG source reconstruction method that simultaneously reconstructs source amplitudes and identifies source interactions across the whole brain. In the proposed method, a full multivariate autoregressive (MAR) model formulates directed interactions (i.e., effective connectivity) between sources. The MAR coefficients (the entries of the MAR matrix) are constrained by the prior knowledge of whole-brain anatomical networks inferred from diffusion MRI. Moreover, to increase the accuracy and robustness of our method, we apply an fMRI prior on the spatial activity patterns and a sparse prior on the MAR coefficients. The observation process of MEG data, the source dynamics, and a series of the priors are combined into a Bayesian framework using a state-space representation. The parameters, such as the source amplitudes and the MAR coefficients, are jointly estimated from a variational Bayesian learning algorithm. By formulating the source dynamics in the context of MEG source reconstruction, and unifying the estimations of source amplitudes and interactions, we can identify the effective connectivity without requiring the selection of regions of interest. Our method is quantitatively and qualitatively evaluated on simulated and experimental data, respectively. Compared with non-dynamic methods, in which the interactions are estimated after source reconstruction with no dynamic constraints, the proposed dynamic method improves most of the performance measures in simulations, and provides better physiological interpretation and inter-subject consistency in real data applications.