Human electromagnetic and haemodynamic networks systematically converge in unimodal cortex and diverge in transmodal cortex

Altered long-range functional connectivity in PTSD: Role of the infraslow oscillations of cortical activity amplitude envelopes

OBJECTIVE

Coupling between the amplitude envelopes (AEs) of regional cortical activity reflects mechanisms that coordinate the excitability of large-scale cortical networks. We used resting-state MEG recordings to investigate the association between alterations in the coupling of cortical AEs and symptoms of post-traumatic stress disorder (PTSD).

METHODS

Participants (n = 96) were service members with combat exposure and various levels of post-traumatic stress severity (PTSS). We assessed the correlation between PTSS and (1) coupling of broadband cortical AEs of beta band activity, (2) coupling of the low- (<0.5 Hz) and high-frequency (>0.5 Hz) components of the AEs, and (3) their time-varying patterns.

RESULTS

PTSS was associated with widespread hypoconnectivity assessed from the broadband AE fluctuations, which correlated with subscores for the negative thoughts and feelings/emotional numbing (NTF/EN) and hyperarousal clusters of symptoms. Higher NTF/EN scores were also associated with smaller increases in resting-state functional connectivity (rsFC) with time during the recordings. The distinct patterns of rsFC in PTSD were primarily due to differences in the coupling of low-frequency (infraslow) fluctuations of the AEs of beta band activity.

CONCLUSIONS

Our findings implicate the mechanisms underlying the regulation/coupling of infraslow oscillations in the alterations of rsFC assessed from broadband AEs and in PTSD symptomatology.

SIGNIFICANCE

Altered coordination of infraslow amplitude fluctuations across large-scale cortical networks can contribute to network dysfunction and may provide a target for treatment in PTSD.

The neurophysiological brain-fingerprint of Parkinson's disease

BACKGROUND

Research in healthy young adults shows that characteristic patterns of brain activity define individual "brain-fingerprints" that are unique to each person. However, variability in these brain-fingerprints increases in individuals with neurological conditions, challenging the clinical relevance and potential impact of the approach. Our study shows that brain-fingerprints derived from neurophysiological brain activity are associated with pathophysiological and clinical traits of individual patients with Parkinson's disease (PD).

METHODS

We created brain-fingerprints from task-free brain activity recorded through magnetoencephalography in 79 PD patients and compared them with those from two independent samples of age-matched healthy controls (N = 424 total). We decomposed brain activity into arrhythmic and rhythmic components, defining distinct brain-fingerprints for each type from recording durations of up to 4 min and as short as 30 s.

FINDINGS

The arrhythmic spectral components of cortical activity in patients with Parkinson's disease are more variable over short periods, challenging the definition of a reliable brain-fingerprint. However, by isolating the rhythmic components of cortical activity, we derived brain-fingerprints that distinguished between patients and healthy controls with about 90% accuracy. The most prominent cortical features of the resulting Parkinson's brain-fingerprint are mapped to polyrhythmic activity in unimodal sensorimotor regions. Leveraging these features, we also demonstrate that Parkinson's symptom laterality can be decoded directly from cortical neurophysiological activity. Furthermore, our study reveals that the cortical topography of the Parkinson's brain-fingerprint aligns with that of neurotransmitter systems affected by the disease's pathophysiology.

INTERPRETATION

The increased moment-to-moment variability of arrhythmic brain-fingerprints challenges patient differentiation and explains previously published results. We outline patient-specific rhythmic brain signaling features that provide insights into both the neurophysiological signature and symptom laterality of Parkinson's disease. Thus, the proposed definition of a rhythmic brain-fingerprint of Parkinson's disease may contribute to novel, refined approaches to patient stratification. Symmetrically, we discuss how rhythmic brain-fingerprints may contribute to the improved identification and testing of therapeutic neurostimulation targets.

FUNDING

Data collection and sharing for this project was provided by the Quebec Parkinson Network (QPN), the Pre-symptomatic Evaluation of Novel or Experimental Treatments for Alzheimer's Disease (PREVENT-AD; release 6.0) program, the Cambridge Centre for Aging Neuroscience (Cam-CAN), and the Open MEG Archives (OMEGA). The QPN is funded by a grant from Fonds de Recherche du Québec - Santé (FRQS). PREVENT-AD was launched in 2011 as a $13.5 million, 7-year public-private partnership using funds provided by McGill University, the FRQS, an unrestricted research grant from Pfizer Canada, the Levesque Foundation, the Douglas Hospital Research Centre and Foundation, the Government of Canada, and the Canada Fund for Innovation. The Brainstorm project is supported by funding to SB from the NIH (R01-EB026299-05). Further funding to SB for this study included a Discovery grant from the Natural Sciences and Engineering Research Council of Canada of Canada (436355-13), and the CIHR Canada research Chair in Neural Dynamics of Brain Systems (CRC-2017-00311).

Deciphering the functional specialization of whole-brain spatiomolecular gradients in the adult brain

Cortical arealization arises during neurodevelopment from the confluence of molecular gradients representing patterned expression of morphogens and transcription factors. However, whether similar gradients are maintained in the adult brain remains unknown. Here, we uncover three axes of topographic variation in gene expression in the adult human brain that specifically capture previously identified rostral-caudal, dorsal-ventral, and medial-lateral axes of early developmental patterning. The interaction of these spatiomolecular gradients i) accurately reconstructs the position of brain tissue samples, ii) delineates known functional territories, and iii) can model the topographical variation of diverse cortical features. The spatiomolecular gradients are distinct from canonical cortical axes differentiating the primary sensory cortex from the association cortex, but radiate in parallel with the axes traversed by local field potentials along the cortex. We replicate all three molecular gradients in three independent human datasets as well as two nonhuman primate datasets and find that each gradient shows a distinct developmental trajectory across the lifespan. The gradients are composed of several well-known transcription factors (e.g., PAX6 and SIX3), and a small set of genes shared across gradients are strongly enriched for multiple diseases. Together, these results provide insight into the developmental sculpting of functionally distinct brain regions, governed by three robust transcriptomic axes embedded within brain parenchyma.

Frequency-specific segregation and integration of human cerebral cortex: An intrinsic functional atlas

The cognitive and behavioral functions of the human brain are supported by its frequency multiplexing mechanism. However, there is limited understanding of the dynamics of the functional network topology. This study aims to investigate the frequency-specific topology of the functional human brain using 7T rs-fMRI data. Frequency-specific parcellations were first performed, revealing frequency-dependent dynamics within the frontoparietal control, parietal memory, and visual networks. An intrinsic functional atlas containing 456 parcels was proposed and validated using stereo-EEG. Graph theory analysis suggested that, in addition to the task-positive vs. task-negative organization observed in static networks, there was a cognitive control system additionally from a frequency perspective. The reproducibility and plausibility of the identified hub sets were confirmed through 3T fMRI analysis, and their artificial removal had distinct effects on network topology. These results indicate a more intricate and subtle dynamics of the functional human brain and emphasize the significance of accurate topography.

Controlling the human connectome with spatially diffuse input signals

The human brain is never at "rest"; its activity is constantly fluctuating over time, transitioning from one brain state-a whole-brain pattern of activity-to another. Network control theory offers a framework for understanding the effort - energy - associated with these transitions. One branch of control theory that is especially useful in this context is "optimal control", in which input signals are used to selectively drive the brain into a target state. Typically, these inputs are introduced independently to the nodes of the network (each input signal is associated with exactly one node). Though convenient, this input strategy ignores the continuity of cerebral cortex - geometrically, each region is connected to its spatial neighbors, allowing control signals, both exogenous and endogenous, to spread from their foci to nearby regions. Additionally, the spatial specificity of brain stimulation techniques is limited, such that the effects of a perturbation are measurable in tissue surrounding the stimulation site. Here, we adapt the network control model so that input signals have a spatial extent that decays exponentially from the input site. We show that this more realistic strategy takes advantage of spatial dependencies in structural connectivity and activity to reduce the energy (effort) associated with brain state transitions. We further leverage these dependencies to explore near-optimal control strategies such that, on a per-transition basis, the number of input signals required for a given control task is reduced, in some cases by two orders of magnitude. This approximation yields network-wide maps of input site density, which we compare to an existing database of functional, metabolic, genetic, and neurochemical maps, finding a close correspondence. Ultimately, not only do we propose a more efficient framework that is also more adherent to well-established brain organizational principles, but we also posit neurobiologically grounded bases for optimal control.

Individual word representations dissociate from linguistic context along a cortical unimodal to heteromodal gradient

Language comprehension involves multiple hierarchical processing stages across time, space, and levels of representation. When processing a word, the sensory input is transformed into increasingly abstract representations that need to be integrated with the linguistic context. Thus, language comprehension involves both input-driven as well as context-dependent processes. While neuroimaging research has traditionally focused on mapping individual brain regions to the distinct underlying processes, recent studies indicate that whole-brain distributed patterns of cortical activation might be highly relevant for cognitive functions, including language. One such pattern, based on resting-state connectivity, is the 'principal cortical gradient', which dissociates sensory from heteromodal brain regions. The present study investigated the extent to which this gradient provides an organizational principle underlying language function, using a multimodal neuroimaging dataset of functional magnetic resonance imaging (fMRI) and magnetoencephalography (MEG) recordings from 102 participants during sentence reading. We found that the brain response to individual representations of a word (word length, orthographic distance, and word frequency), which reflect visual; orthographic; and lexical properties, gradually increases towards the sensory end of the gradient. Although these properties showed opposite effect directions in fMRI and MEG, their association with the sensory end of the gradient was consistent across both neuroimaging modalities. In contrast, MEG revealed that properties reflecting a word's relation to its linguistic context (semantic similarity and position within the sentence) involve the heteromodal end of the gradient to a stronger extent. This dissociation between individual word and contextual properties was stable across earlier and later time windows during word presentation, indicating interactive processing of word representations and linguistic context at opposing ends of the principal gradient. To conclude, our findings indicate that the principal gradient underlies the organization of a range of linguistic representations while supporting a gradual distinction between context-independent and context-dependent representations. Furthermore, the gradient reveals convergent patterns across neuroimaging modalities (similar location along the gradient) in the presence of divergent responses (opposite effect directions).

The neurophysiological brain-fingerprint of Parkinson's disease

In this study, we investigate the clinical potential of brain-fingerprints derived from electrophysiological brain activity for diagnostics and progression monitoring of Parkinson's disease (PD). We obtained brain-fingerprints from PD patients and age-matched healthy controls using short, task-free magnetoencephalographic recordings. The rhythmic components of the individual brain-fingerprint distinguished between patients and healthy participants with approximately 90% accuracy. The most prominent cortical features of the Parkinson's brain-fingerprint mapped to polyrhythmic activity in unimodal sensorimotor regions. Leveraging these features, we also show that Parkinson's disease stages can be decoded directly from cortical neurophysiological activity. Additionally, our study reveals that the cortical topography of the Parkinson's brain-fingerprint aligns with that of neurotransmitter systems affected by the disease's pathophysiology. We further demonstrate that the arrhythmic components of cortical activity are more variable over short periods of time in patients with Parkinson's disease than in healthy controls, making individual differentiation between patients based on these features more challenging and explaining previous negative published results. Overall, we outline patient-specific rhythmic brain signaling features that provide insights into both the neurophysiological signature and clinical staging of Parkinson's disease. For this reason, the proposed definition of a rhythmic brain-fingerprint of Parkinson's disease may contribute to novel, refined approaches to patient stratification and to the improved identification and testing of therapeutic neurostimulation targets.

Integrating brainstem and cortical functional architectures

The brainstem is a fundamental component of the central nervous system yet it is typically excluded from in vivo human brain mapping efforts, precluding a complete understanding of how the brainstem influences cortical function. Here we use high-resolution 7 Tesla fMRI to derive a functional connectome encompassing cortex as well as 58 brainstem nuclei spanning the midbrain, pons and medulla. We identify a compact set of integrative hubs in the brainstem with widespread connectivity with cerebral cortex. Patterns of connectivity between brainstem and cerebral cortex manifest as multiple emergent phenomena including neurophysiological oscillatory rhythms, patterns of cognitive functional specialization, and the unimodal-transmodal functional hierarchy. This persistent alignment between cortical functional topographies and brainstem nuclei is shaped by the spatial arrangement of multiple neurotransmitter receptors and transporters. We replicate all findings using 3 Tesla data from the same participants. Collectively, we find that multiple organizational features of cortical activity can be traced back to the brainstem.

Living on the edge: network neuroscience beyond nodes

Network neuroscience has emphasized the connectional properties of neural elements - cells, populations, and regions. This has come at the expense of the anatomical and functional connections that link these elements to one another. A new perspective - namely one that emphasizes 'edges' - may prove fruitful in addressing outstanding questions in network neuroscience. We highlight one recently proposed 'edge-centric' method and review its current applications, merits, and limitations. We also seek to establish conceptual and mathematical links between this method and previously proposed approaches in the network science and neuroimaging literature. We conclude by presenting several avenues for future work to extend and refine existing edge-centric analysis.

Integrating brainstem and cortical functional architectures

The brainstem is a fundamental component of the central nervous system yet it is typically excluded from in vivo human brain mapping efforts, precluding a complete understanding of how the brainstem influences cortical function. Here we use high-resolution 7 Tesla fMRI to derive a functional connectome encompassing cortex as well as 58 brainstem nuclei spanning the midbrain, pons and medulla. We identify a compact set of integrative hubs in the brainstem with widespread connectivity with cerebral cortex. Patterns of connectivity between brainstem and cerebral cortex manifest as multiple emergent phenomena including neurophysiological oscillatory rhythms, patterns of cognitive functional specialization, and the unimodal-transmodal functional hierarchy. This persistent alignment between cortical functional topographies and brainstem nuclei is shaped by the spatial arrangement of multiple neurotransmitter receptors and transporters. We replicate all findings using 3 Tesla data from the same participants. Collectively, we find that multiple organizational features of cortical activity can be traced back to the brainstem.

Network connectivity and local transcriptomic vulnerability underpin cortical atrophy progression in Parkinson's disease

Parkinson's disease pathology is hypothesized to spread through the brain via axonal connections between regions and is further modulated by local vulnerabilities within those regions. The resulting changes to brain morphology have previously been demonstrated in both prodromal and de novo Parkinson's disease patients. However, it remains unclear whether the pattern of atrophy progression in Parkinson's disease over time is similarly explained by network-based spreading and local vulnerability. We address this gap by mapping the trajectory of cortical atrophy rates in a large, multi-centre cohort of Parkinson's disease patients and relate this atrophy progression pattern to network architecture and gene expression profiles. Across 4-year follow-up visits, increased atrophy rates were observed in posterior, temporal, and superior frontal cortices. We demonstrated that this progression pattern was shaped by network connectivity. Regional atrophy rates were strongly related to atrophy rates across structurally and functionally connected regions. We also found that atrophy progression was associated with specific gene expression profiles. The genes whose spatial distribution in the brain was most related to atrophy rate were those enriched for mitochondrial and metabolic function. Taken together, our findings demonstrate that both global and local brain features influence vulnerability to neurodegeneration in Parkinson's disease.

Resolving inter-regional communication capacity in the human connectome

Applications of graph theory to the connectome have inspired several models of how neural signaling unfolds atop its structure. Analytic measures derived from these communication models have mainly been used to extract global characteristics of brain networks, obscuring potentially informative inter-regional relationships. Here we develop a simple standardization method to investigate polysynaptic communication pathways between pairs of cortical regions. This procedure allows us to determine which pairs of nodes are topologically closer and which are further than expected on the basis of their degree. We find that communication pathways delineate canonical functional systems. Relating nodal communication capacity to meta-analytic probabilistic patterns of functional specialization, we also show that areas that are most closely integrated within the network are associated with higher order cognitive functions. We find that these regions' proclivity towards functional integration could naturally arise from the brain's anatomical configuration through evenly distributed connections among multiple specialized communities. Throughout, we consider two increasingly constrained null models to disentangle the effects of the network's topology from those passively endowed by spatial embedding. Altogether, the present findings uncover relationships between polysynaptic communication pathways and the brain's functional organization across multiple topological levels of analysis and demonstrate that network integration facilitates cognitive integration.

Integrating multimodal and multiscale connectivity blueprints of the human cerebral cortex in health and disease

The brain is composed of disparate neural populations that communicate and interact with one another. Although fiber bundles, similarities in molecular architecture, and synchronized neural activity all reflect how brain regions potentially interact with one another, a comprehensive study of how all these interregional relationships jointly reflect brain structure and function remains missing. Here, we systematically integrate 7 multimodal, multiscale types of interregional similarity ("connectivity modes") derived from gene expression, neurotransmitter receptor density, cellular morphology, glucose metabolism, haemodynamic activity, and electrophysiology in humans. We first show that for all connectivity modes, feature similarity decreases with distance and increases when regions are structurally connected. Next, we show that connectivity modes exhibit unique and diverse connection patterns, hub profiles, spatial gradients, and modular organization. Throughout, we observe a consistent primacy of molecular connectivity modes-namely correlated gene expression and receptor similarity-that map onto multiple phenomena, including the rich club and patterns of abnormal cortical thickness across 13 neurological, psychiatric, and neurodevelopmental disorders. Finally, to construct a single multimodal wiring map of the human cortex, we fuse all 7 connectivity modes and show that the fused network maps onto major organizational features of the cortex including structural connectivity, intrinsic functional networks, and cytoarchitectonic classes. Altogether, this work contributes to the integrative study of interregional relationships in the human cerebral cortex.

Spatially heterogeneous structure-function coupling in haemodynamic and electromagnetic brain networks

The relationship between structural and functional connectivity in the brain is a key question in connectomics. Here we quantify patterns of structure-function coupling across the neocortex, by comparing structural connectivity estimated using diffusion MRI with functional connectivity estimated using both neurophysiological (MEG-based) and haemodynamic (fMRI-based) recordings. We find that structure-function coupling is heterogeneous across brain regions and frequency bands. The link between structural and functional connectivity is generally stronger in multiple MEG frequency bands compared to resting state fMRI. Structure-function coupling is greater in slower and intermediate frequency bands compared to faster frequency bands. We also find that structure-function coupling systematically follows the archetypal sensorimotor-association hierarchy, as well as patterns of laminar differentiation, peaking in granular layer IV. Finally, structure-function coupling is better explained using structure-informed inter-regional communication metrics than using structural connectivity alone. Collectively, these results place neurophysiological and haemodynamic structure-function relationships in a common frame of reference and provide a starting point for a multi-modal understanding of structure-function coupling in the brain.

Adult lifespan trajectories of neuromagnetic signals and interrelations with cortical thickness

Oscillatory power and phase synchronization map neuronal dynamics and are commonly studied to differentiate the healthy and diseased brain. Yet, little is known about the course and spatial variability of these features from early adulthood into old age. Leveraging magnetoencephalography (MEG) resting-state data in a cross-sectional adult sample (n = 350), we probed lifespan differences (18-88 years) in connectivity and power and interaction effects with sex. Building upon recent attempts to link brain structure and function, we tested the spatial correspondence between age effects on cortical thickness and those on functional networks. We further probed a direct structure-function relationship at the level of the study sample. We found MEG frequency-specific patterns with age and divergence between sexes in low frequencies. Connectivity and power exhibited distinct linear trajectories or turning points at midlife that might reflect different physiological processes. In the delta and beta bands, these age effects corresponded to those on cortical thickness, pointing to co-variation between the modalities across the lifespan. Structure-function coupling was frequency-dependent and observed in unimodal or multimodal regions. Altogether, we provide a comprehensive overview of the topographic functional profile of adulthood that can form a basis for neurocognitive and clinical investigations. This study further sheds new light on how the brain's structural architecture relates to fast oscillatory activity.

Multi-modal and multi-model interrogation of large-scale functional brain networks

Existing whole-brain models are generally tailored to the modelling of a particular data modality (e.g., fMRI or MEG/EEG). We propose that despite the differing aspects of neural activity each modality captures, they originate from shared network dynamics. Building on the universal principles of self-organising delay-coupled nonlinear systems, we aim to link distinct features of brain activity - captured across modalities - to the dynamics unfolding on a macroscopic structural connectome. To jointly predict connectivity, spatiotemporal and transient features of distinct signal modalities, we consider two large-scale models - the Stuart Landau and Wilson and Cowan models - which generate short-lived 40 Hz oscillations with varying levels of realism. To this end, we measure features of functional connectivity and metastable oscillatory modes (MOMs) in fMRI and MEG signals - and compare them against simulated data. We show that both models can represent MEG functional connectivity (FC), functional connectivity dynamics (FCD) and generate MOMs to a comparable degree. This is achieved by adjusting the global coupling and mean conduction time delay and, in the WC model, through the inclusion of balance between excitation and inhibition. For both models, the omission of delays dramatically decreased the performance. For fMRI, the SL model performed worse for FCD and MOMs, highlighting the importance of balanced dynamics for the emergence of spatiotemporal and transient patterns of ultra-slow dynamics. Notably, optimal working points varied across modalities and no model was able to achieve a correlation with empirical FC higher than 0.4 across modalities for the same set of parameters. Nonetheless, both displayed the emergence of FC patterns that extended beyond the constraints of the anatomical structure. Finally, we show that both models can generate MOMs with empirical-like properties such as size (number of brain regions engaging in a mode) and duration (continuous time interval during which a mode appears). Our results demonstrate the emergence of static and dynamic properties of neural activity at different timescales from networks of delay-coupled oscillators at 40 Hz. Given the higher dependence of simulated FC on the underlying structural connectivity, we suggest that mesoscale heterogeneities in neural circuitry may be critical for the emergence of parallel cross-modal functional networks and should be accounted for in future modelling endeavours.

A rapid theta network mechanism for flexible information encoding

Flexible behavior requires gating mechanisms that encode only task-relevant information in working memory. Extant literature supports a theoretical division of labor whereby lateral frontoparietal interactions underlie information maintenance and the striatum enacts the gate. Here, we reveal neocortical gating mechanisms in intracranial EEG patients by identifying rapid, within-trial changes in regional and inter-regional activities that predict subsequent behavioral outputs. Results first demonstrate information accumulation mechanisms that extend prior fMRI (i.e., regional high-frequency activity) and EEG evidence (inter-regional theta synchrony) of distributed neocortical networks in working memory. Second, results demonstrate that rapid changes in theta synchrony, reflected in changing patterns of default mode network connectivity, support filtering. Graph theoretical analyses further linked filtering in task-relevant information and filtering out irrelevant information to dorsal and ventral attention networks, respectively. Results establish a rapid neocortical theta network mechanism for flexible information encoding, a role previously attributed to the striatum.

Modules in connectomes of phase-synchronization comprise anatomically contiguous, functionally related regions

Modules in brain functional connectomes are essential to balancing segregation and integration of neuronal activity. Connectomes are the complete set of pairwise connections between brain regions. Non-invasive Electroencephalography (EEG) and Magnetoencephalography (MEG) have been used to identify modules in connectomes of phase-synchronization. However, their resolution is suboptimal because of spurious phase-synchronization due to EEG volume conduction or MEG field spread. Here, we used invasive, intracerebral recordings from stereo-electroencephalography (SEEG, N = 67), to identify modules in connectomes of phase-synchronization. To generate SEEG-based group-level connectomes affected only minimally by volume conduction, we used submillimeter accurate localization of SEEG contacts and referenced electrode contacts in cortical gray matter to their closest contacts in white matter. Combining community detection methods with consensus clustering, we found that the connectomes of phase-synchronization were characterized by distinct and stable modules at multiple spatial scales, across frequencies from 3 to 320 Hz. These modules were highly similar within canonical frequency bands. Unlike the distributed brain systems identified with functional Magnetic Resonance Imaging (fMRI), modules up to the high-gamma frequency band comprised only anatomically contiguous regions. Notably, the identified modules comprised cortical regions involved in shared repertoires of sensorimotor and cognitive functions including memory, language and attention. These results suggest that the identified modules represent functionally specialised brain systems, which only partially overlap with the brain systems reported with fMRI. Hence, these modules might regulate the balance between functional segregation and functional integration through phase-synchronization.

Task matters: Individual MEG signatures from naturalistic and neurophysiological brain states

The discovery that human brain connectivity data can be used as a "fingerprint" to identify a given individual from a population, has become a burgeoning research area in the neuroscience field. Recent studies have identified the possibility to extract these brain signatures from the temporal rich dynamics of resting-state magneto encephalography (MEG) recordings. Nevertheless, it is still uncertain to what extent MEG signatures can serve as an indicator of human identifiability during task-related conduct. Here, using MEG data from naturalistic and neurophysiological tasks, we show that identification improves in tasks relative to resting-state, providing compelling evidence for a task dependent axis of MEG signatures. Notably, improvements in identifiability were more prominent in strictly controlled tasks. Lastly, the brain regions contributing most towards individual identification were also modified when engaged in task activities. We hope that this investigation advances our understanding of the driving factors behind brain identification from MEG signals.

Connectional Hierarchy in Human Brain Revealed by Individual Variability of Functional Network Edges

The human cerebral cortex is connected by intricate inter-areal wiring at the macroscale. The cortical hierarchy from primary sensorimotor to higher-order association areas is a unifying organizational principle across various neurobiological properties; however, previous studies have not clarified whether the connections between cortical regions exhibit a similar hierarchical pattern. Here, we identify a connectional hierarchy indexed by inter-individual variability of functional connectivity edges, which continuously progresses along a hierarchical gradient from within-network connections to between-network edges connecting sensorimotor and association networks. We found that this connectional hierarchy of variability aligns with both hemodynamic and electromagnetic connectivity strength and is constrained by structural connectivity strength. Moreover, the patterning of connectional hierarchy is related to inter-regional similarity in transcriptional and neurotransmitter receptor profiles. Using the Neurosynth cognitive atlas and cortical vulnerability maps in 13 brain disorders, we found that the connectional hierarchy of variability is associated with similarity networks of cognitive relevance and that of disorder vulnerability. Finally, we found that the prominence of this hierarchical gradient of connectivity variability declines during youth. Together, our results reveal a novel hierarchal organizational principle at the connectional level that links multimodal and multiscale human connectomes to individual variability in functional connectivity.

Correspondence between gene expression and neurotransmitter receptor and transporter density in the human brain

Neurotransmitter receptors modulate signaling between neurons. Thus, neurotransmitter receptors and transporters play a key role in shaping brain function. Due to the lack of comprehensive neurotransmitter receptor/transporter density datasets, microarray gene expression measuring mRNA transcripts is often used as a proxy for receptor densities. In the present report, we comprehensively test the spatial correlation between gene expression and protein density for a total of 27 neurotransmitter receptors, receptor binding-sites, and transporters across 9 different neurotransmitter systems, using both PET and autoradiography radioligand-based imaging modalities. We find poor spatial correspondences between gene expression and density for all neurotransmitter receptors and transporters except four single-protein metabotropic receptors (5-HT_1A, CB₁, D₂, and MOR). These expression-density associations are related to gene differential stability and can vary between cortical and subcortical structures. Altogether, we recommend using direct measures of receptor and transporter density when relating neurotransmitter systems to brain structure and function.

Mapping neurotransmitter systems to the structural and functional organization of the human neocortex

Neurotransmitter receptors support the propagation of signals in the human brain. How receptor systems are situated within macro-scale neuroanatomy and how they shape emergent function remain poorly understood, and there exists no comprehensive atlas of receptors. Here we collate positron emission tomography data from more than 1,200 healthy individuals to construct a whole-brain three-dimensional normative atlas of 19 receptors and transporters across nine different neurotransmitter systems. We found that receptor profiles align with structural connectivity and mediate function, including neurophysiological oscillatory dynamics and resting-state hemodynamic functional connectivity. Using the Neurosynth cognitive atlas, we uncovered a topographic gradient of overlapping receptor distributions that separates extrinsic and intrinsic psychological processes. Finally, we found both expected and novel associations between receptor distributions and cortical abnormality patterns across 13 disorders. We replicated all findings in an independently collected autoradiography dataset. This work demonstrates how chemoarchitecture shapes brain structure and function, providing a new direction for studying multi-scale brain organization.