Frontostriatal salience network expansion in individuals in depression.

Decades of neuroimaging studies have shown modest differences in brain structure and connectivity in depression, hindering mechanistic insights or the identification of risk factors for disease onset1. Furthermore, whereas depression is episodic, few longitudinal neuroimaging studies exist, limiting understanding of mechanisms that drive mood-state transitions. The emerging field of precision functional mapping has used densely sampled longitudinal neuroimaging data to show behaviourally meaningful differences in brain network topography and connectivity between and in healthy individuals2-4, but this approach has not been applied in depression. Here, using precision functional mapping and several samples of deeply sampled individuals, we found that the frontostriatal salience network is expanded nearly twofold in the cortex of most individuals with depression. This effect was replicable in several samples and caused primarily by network border shifts, with three distinct modes of encroachment occurring in different individuals. Salience network expansion was stable over time, unaffected by mood state and detectable in children before the onset of depression later in adolescence. Longitudinal analyses of individuals scanned up to 62 times over 1.5 years identified connectivity changes in frontostriatal circuits that tracked fluctuations in specific symptoms and predicted future anhedonia symptoms. Together, these findings identify a trait-like brain network topology that may confer risk for depression and mood-state-dependent connectivity changes in frontostriatal circuits that predict the emergence and remission of depressive symptoms over time.

Psilocybin desynchronizes the human brain.

A single dose of psilocybin, a psychedelic that acutely causes distortions of space-time perception and ego dissolution, produces rapid and persistent therapeutic effects in human clinical trials1-4. In animal models, psilocybin induces neuroplasticity in cortex and hippocampus5-8. It remains unclear how human brain network changes relate to subjective and lasting effects of psychedelics. Here we tracked individual-specific brain changes with longitudinal precision functional mapping (roughly 18 magnetic resonance imaging visits per participant). Healthy adults were tracked before, during and for 3 weeks after high-dose psilocybin (25 mg) and methylphenidate (40 mg), and brought back for an additional psilocybin dose 6-12 months later. Psilocybin massively disrupted functional connectivity (FC) in cortex and subcortex, acutely causing more than threefold greater change than methylphenidate. These FC changes were driven by brain desynchronization across spatial scales (areal, global), which dissolved network distinctions by reducing correlations within and anticorrelations between networks. Psilocybin-driven FC changes were strongest in the default mode network, which is connected to the anterior hippocampus and is thought to create our sense of space, time and self. Individual differences in FC changes were strongly linked to the subjective psychedelic experience. Performing a perceptual task reduced psilocybin-driven FC changes. Psilocybin caused persistent decrease in FC between the anterior hippocampus and default mode network, lasting for weeks. Persistent reduction of hippocampal-default mode network connectivity may represent a neuroanatomical and mechanistic correlate of the proplasticity and therapeutic effects of psychedelics.

A somato-cognitive action network alternates with effector regions in motor cortex.

Motor cortex (M1) has been thought to form a continuous somatotopic homunculus extending down the precentral gyrus from foot to face representations1,2, despite evidence for concentric functional zones3 and maps of complex actions4. Here, using precision functional magnetic resonance imaging (fMRI) methods, we find that the classic homunculus is interrupted by regions with distinct connectivity, structure and function, alternating with effector-specific (foot, hand and mouth) areas. These inter-effector regions exhibit decreased cortical thickness and strong functional connectivity to each other, as well as to the cingulo-opercular network (CON), critical for action5 and physiological control6, arousal7, errors8 and pain9. This interdigitation of action control-linked and motor effector regions was verified in the three largest fMRI datasets. Macaque and pediatric (newborn, infant and child) precision fMRI suggested cross-species homologues and developmental precursors of the inter-effector system. A battery of motor and action fMRI tasks documented concentric effector somatotopies, separated by the CON-linked inter-effector regions. The inter-effectors lacked movement specificity and co-activated during action planning (coordination of hands and feet) and axial body movement (such as of the abdomen or eyebrows). These results, together with previous studies demonstrating stimulation-evoked complex actions4 and connectivity to internal organs10 such as the adrenal medulla, suggest that M1 is punctuated by a system for whole-body action planning, the somato-cognitive action network (SCAN). In M1, two parallel systems intertwine, forming an integrate-isolate pattern: effector-specific regions (foot, hand and mouth) for isolating fine motor control and the SCAN for integrating goals, physiology and body movement.

Structure-function coupling in highly sampled individual brains.

Structural connectivity (SC) between distant regions of the brain support synchronized function known as functional connectivity (FC) and give rise to the large-scale brain networks that enable cognition and behavior. Understanding how SC enables FC is important to understand how injuries to SC may alter brain function and cognition. Previous work evaluating whole-brain SC-FC relationships showed that SC explained FC well in unimodal visual and motor areas, but only weakly in association areas, suggesting a unimodal-heteromodal gradient organization of SC-FC coupling. However, this work was conducted in group-averaged SC/FC data. Thus, it could not account for inter-individual variability in the locations of cortical areas and white matter tracts. We evaluated the correspondence of SC and FC within three highly sampled healthy participants. For each participant, we collected 78 min of diffusion-weighted MRI for SC and 360 min of resting state fMRI for FC. We found that FC was best explained by SC in visual and motor systems, as well as in anterior and posterior cingulate regions. A unimodal-to-heteromodal gradient could not fully explain SC-FC coupling. We conclude that the SC-FC coupling of the anterior-posterior cingulate circuit is more similar to unimodal areas than to heteromodal areas.

Functional parcellation of the neonatal cortical surface.

The cerebral cortex is organized into distinct but interconnected cortical areas, which can be defined by abrupt differences in patterns of resting state functional connectivity (FC) across the cortical surface. Such parcellations of the cortex have been derived in adults and older infants, but there is no widely used surface parcellation available for the neonatal brain. Here, we first demonstrate that existing parcellations, including surface-based parcels derived from older samples as well as volume-based neonatal parcels, are a poor fit for neonatal surface data. We next derive a set of 283 cortical surface parcels from a sample of n = 261 neonates. These parcels have highly homogenous FC patterns and are validated using three external neonatal datasets. The Infomap algorithm is used to assign functional network identities to each parcel, and derived networks are consistent with prior work in neonates. The proposed parcellation may represent neonatal cortical areas and provides a powerful tool for neonatal neuroimaging studies.

Thalamo-cortical and cerebello-cortical functional connectivity in development.

The thalamus is a critical relay center for neural pathways involving sensory, motor, and cognitive functions, including cortico-striato-thalamo-cortical and cortico-ponto-cerebello-thalamo-cortical loops. Despite the importance of these circuits, their development has been understudied. One way to investigate these pathways in human development in vivo is with functional connectivity MRI, yet few studies have examined thalamo-cortical and cerebello-cortical functional connectivity in development. Here, we used resting-state functional connectivity to measure functional connectivity in the thalamus and cerebellum with previously defined cortical functional networks in 2 separate data sets of children (7-12 years old) and adults (19-40 years old). In both data sets, we found stronger functional connectivity between the ventral thalamus and the somatomotor face cortical functional network in children compared with adults, extending previous cortico-striatal functional connectivity findings. In addition, there was more cortical network integration (i.e. strongest functional connectivity with multiple networks) in the thalamus in children than in adults. We found no developmental differences in cerebello-cortical functional connectivity. Together, these results suggest different maturation patterns in cortico-striato-thalamo-cortical and cortico-ponto-cerebellar-thalamo-cortical pathways.

Network-specific selectivity of functional connections in the neonatal brain.

The adult human brain is organized into functional brain networks, groups of functionally connected segregated brain regions. A key feature of adult functional networks is long-range selectivity, the property that spatially distant regions from the same network have higher functional connectivity than spatially distant regions from different networks. Although it is critical to establish the status of functional networks and long-range selectivity during the neonatal period as a foundation for typical and atypical brain development, prior work in this area has been mixed. Although some studies report distributed adult-like networks, other studies suggest that neonatal networks are immature and consist primarily of spatially isolated regions. Using a large sample of neonates (n = 262), we demonstrate that neonates have long-range selective functional connections for the default mode, fronto-parietal, and dorsal attention networks. An adult-like pattern of functional brain networks is evident in neonates when network-detection algorithms are tuned to these long-range connections, when using surface-based registration (versus volume-based registration), and as per-subject data quantity increases. These results help clarify factors that have led to prior mixed results, establish that key adult-like functional network features are evident in neonates, and provide a foundation for studies of typical and atypical brain development.

Cingulo-opercular control network and disused motor circuits joined in standby mode.

Whole-brain resting-state functional MRI (rs-fMRI) during 2 wk of upper-limb casting revealed that disused motor regions became more strongly connected to the cingulo-opercular network (CON), an executive control network that includes regions of the dorsal anterior cingulate cortex (dACC) and insula. Disuse-driven increases in functional connectivity (FC) were specific to the CON and somatomotor networks and did not involve any other networks, such as the salience, frontoparietal, or default mode networks. Censoring and modeling analyses showed that FC increases during casting were mediated by large, spontaneous activity pulses that appeared in the disused motor regions and CON control regions. During limb constraint, disused motor circuits appear to enter a standby mode characterized by spontaneous activity pulses and strengthened connectivity to CON executive control regions.

Parallel hippocampal-parietal circuits for self- and goal-oriented processing.

The hippocampus is critically important for a diverse range of cognitive processes, such as episodic memory, prospective memory, affective processing, and spatial navigation. Using individual-specific precision functional mapping of resting-state functional MRI data, we found the anterior hippocampus (head and body) to be preferentially functionally connected to the default mode network (DMN), as expected. The hippocampal tail, however, was strongly preferentially functionally connected to the parietal memory network (PMN), which supports goal-oriented cognition and stimulus recognition. This anterior-posterior dichotomy of resting-state functional connectivity was well-matched by differences in task deactivations and anatomical segmentations of the hippocampus. Task deactivations were localized to the hippocampal head and body (DMN), relatively sparing the tail (PMN). The functional dichotomization of the hippocampus into anterior DMN-connected and posterior PMN-connected parcels suggests parallel but distinct circuits between the hippocampus and medial parietal cortex for self- versus goal-oriented processing.

Brain hubs defined in the group do not overlap with regions of high inter-individual variability.

Connector 'hubs' are brain regions with links to multiple networks. These regions are hypothesized to play a critical role in brain function. While hubs are often identified based on group-average functional magnetic resonance imaging (fMRI) data, there is considerable inter-subject variation in the functional connectivity profiles of the brain, especially in association regions where hubs tend to be located. Here we investigated how group hubs are related to locations of inter-individual variability. To answer this question, we examined inter-individual variation at group-level hubs in both the Midnight Scan Club and Human Connectome Project datasets. The top group hubs defined based on the participation coefficient did not overlap strongly with the most prominent regions of inter-individual variation (termed 'variants' in prior work). These hubs have relatively strong similarity across participants and consistent cross-network profiles, similar to what was seen for many other areas of cortex. Consistency across participants was further improved when these hubs were allowed to shift slightly in local position. Thus, our results demonstrate that the top group hubs defined with the participation coefficient are generally consistent across people, suggesting they may represent conserved cross-network bridges. More caution is warranted with alternative hub measures, such as community density (which are based on spatial proximity to network borders) and intermediate hub regions which show higher correspondence to locations of individual variability.

Neuroimaging-based classification of PTSD using data-driven computational approaches: A multisite big data study from the ENIGMA-PGC PTSD consortium.

BACKGROUND: Recent advances in data-driven computational approaches have been helpful in devising tools to objectively diagnose psychiatric disorders. However, current machine learning studies limited to small homogeneous samples, different methodologies, and different imaging collection protocols, limit the ability to directly compare and generalize their results. Here we aimed to classify individuals with PTSD versus controls and assess the generalizability using a large heterogeneous brain datasets from the ENIGMA-PGC PTSD Working group. METHODS: We analyzed brain MRI data from 3,477 structural-MRI; 2,495 resting state-fMRI; and 1,952 diffusion-MRI. First, we identified the brain features that best distinguish individuals with PTSD from controls using traditional machine learning methods. Second, we assessed the utility of the denoising variational autoencoder (DVAE) and evaluated its classification performance. Third, we assessed the generalizability and reproducibility of both models using leave-one-site-out cross-validation procedure for each modality. RESULTS: We found lower performance in classifying PTSD vs. controls with data from over 20 sites (60 % test AUC for s-MRI, 59 % for rs-fMRI and 56 % for d-MRI), as compared to other studies run on single-site data. The performance increased when classifying PTSD from HC without trauma history in each modality (75 % AUC). The classification performance remained intact when applying the DVAE framework, which reduced the number of features. Finally, we found that the DVAE framework achieved better generalization to unseen datasets compared with the traditional machine learning frameworks, albeit performance was slightly above chance. CONCLUSION: These results have the potential to provide a baseline classification performance for PTSD when using large scale neuroimaging datasets. Our findings show that the control group used can heavily affect classification performance. The DVAE framework provided better generalizability for the multi-site data. This may be more significant in clinical practice since the neuroimaging-based diagnostic DVAE classification models are much less site-specific, rendering them more generalizable.

Individualized Functional Subnetworks Connect Human Striatum and Frontal Cortex.

The striatum and cerebral cortex are interconnected via multiple recurrent loops that play a major role in many neuropsychiatric conditions. Primate corticostriatal connections can be precisely mapped using invasive tract-tracing. However, noninvasive human research has not mapped these connections with anatomical precision, limited in part by the practice of averaging neuroimaging data across individuals. Here we utilized highly sampled resting-state functional connectivity MRI for individual-specific precision functional mapping (PFM) of corticostriatal connections. We identified ten individual-specific subnetworks linking cortex-predominately frontal cortex-to striatum, most of which converged with nonhuman primate tract-tracing work. These included separable connections between nucleus accumbens core/shell and orbitofrontal/medial frontal gyrus; between anterior striatum and dorsomedial prefrontal cortex; between dorsal caudate and lateral prefrontal cortex; and between middle/posterior putamen and supplementary motor/primary motor cortex. Two subnetworks that did not converge with nonhuman primates were connected to cortical regions associated with human language function. Thus, precision subnetworks identify detailed, individual-specific, neurobiologically plausible corticostriatal connectivity that includes human-specific language networks.

Default-mode network streams for coupling to language and control systems.

The human brain is organized into large-scale networks identifiable using resting-state functional connectivity (RSFC). These functional networks correspond with broad cognitive domains; for example, the Default-mode network (DMN) is engaged during internally oriented cognition. However, functional networks may contain hierarchical substructures corresponding with more specific cognitive functions. Here, we used individual-specific precision RSFC to test whether network substructures could be identified in 10 healthy human brains. Across all subjects and networks, individualized network subdivisions were more valid-more internally homogeneous and better matching spatial patterns of task activation-than canonical networks. These measures of validity were maximized at a hierarchical scale that contained ∼83 subnetworks across the brain. At this scale, nine DMN subnetworks exhibited topographical similarity across subjects, suggesting that this approach identifies homologous neurobiological circuits across individuals. Some DMN subnetworks matched known features of brain organization corresponding with cognitive functions. Other subnetworks represented separate streams by which DMN couples with other canonical large-scale networks, including language and control networks. Together, this work provides a detailed organizational framework for studying the DMN in individual humans.

Individual-specific functional connectivity of the amygdala: A substrate for precision psychiatry.

The amygdala is central to the pathophysiology of many psychiatric illnesses. An imprecise understanding of how the amygdala fits into the larger network organization of the human brain, however, limits our ability to create models of dysfunction in individual patients to guide personalized treatment. Therefore, we investigated the position of the amygdala and its functional subdivisions within the network organization of the brain in 10 highly sampled individuals (5 h of fMRI data per person). We characterized three functional subdivisions within the amygdala of each individual. We discovered that one subdivision is preferentially correlated with the default mode network; a second is preferentially correlated with the dorsal attention and fronto-parietal networks; and third subdivision does not have any networks to which it is preferentially correlated relative to the other two subdivisions. All three subdivisions are positively correlated with ventral attention and somatomotor networks and negatively correlated with salience and cingulo-opercular networks. These observations were replicated in an independent group dataset of 120 individuals. We also found substantial across-subject variation in the distribution and magnitude of amygdala functional connectivity with the cerebral cortex that related to individual differences in the stereotactic locations both of amygdala subdivisions and of cortical functional brain networks. Finally, using lag analyses, we found consistent temporal ordering of fMRI signals in the cortex relative to amygdala subdivisions. Altogether, this work provides a detailed framework of amygdala-cortical interactions that can be used as a foundation for models relating aberrations in amygdala connectivity to psychiatric symptoms in individual patients.

Accuracy and reliability of diffusion imaging models.

Diffusion imaging aims to non-invasively characterize the anatomy and integrity of the brain's white matter fibers. We evaluated the accuracy and reliability of commonly used diffusion imaging methods as a function of data quantity and analysis method, using both simulations and highly sampled individual-specific data (927-1442 diffusion weighted images [DWIs] per individual). Diffusion imaging methods that allow for crossing fibers (FSL's BedpostX [BPX], DSI Studio's Constant Solid Angle Q-Ball Imaging [CSA-QBI], MRtrix3's Constrained Spherical Deconvolution [CSD]) estimated excess fibers when insufficient data were present and/or when the data did not match the model priors. To reduce such overfitting, we developed a novel Bayesian Multi-tensor Model-selection (BaMM) method and applied it to the popular ball-and-stick model used in BedpostX within the FSL software package. BaMM was robust to overfitting and showed high reliability and the relatively best crossing-fiber accuracy with increasing amounts of diffusion data. Thus, sufficient data and an overfitting resistant analysis method enhance precision diffusion imaging. For potential clinical applications of diffusion imaging, such as neurosurgical planning and deep brain stimulation (DBS), the quantities of data required to achieve diffusion imaging reliability are lower than those needed for functional MRI.

A comparison of methods to harmonize cortical thickness measurements across scanners and sites.

Results of neuroimaging datasets aggregated from multiple sites may be biased by site-specific profiles in participants' demographic and clinical characteristics, as well as MRI acquisition protocols and scanning platforms. We compared the impact of four different harmonization methods on results obtained from analyses of cortical thickness data: (1) linear mixed-effects model (LME) that models site-specific random intercepts (LMEINT), (2) LME that models both site-specific random intercepts and age-related random slopes (LMEINT+SLP), (3) ComBat, and (4) ComBat with a generalized additive model (ComBat-GAM). Our test case for comparing harmonization methods was cortical thickness data aggregated from 29 sites, which included 1,340 cases with posttraumatic stress disorder (PTSD) (6.2-81.8 years old) and 2,057 trauma-exposed controls without PTSD (6.3-85.2 years old). We found that, compared to the other data harmonization methods, data processed with ComBat-GAM was more sensitive to the detection of significant case-control differences (Χ2(3) = 63.704, p < 0.001) as well as case-control differences in age-related cortical thinning (Χ2(3) = 12.082, p = 0.007). Both ComBat and ComBat-GAM outperformed LME methods in detecting sex differences (Χ2(3) = 9.114, p = 0.028) in regional cortical thickness. ComBat-GAM also led to stronger estimates of age-related declines in cortical thickness (corrected p-values < 0.001), stronger estimates of case-related cortical thickness reduction (corrected p-values < 0.001), weaker estimates of age-related declines in cortical thickness in cases than controls (corrected p-values < 0.001), stronger estimates of cortical thickness reduction in females than males (corrected p-values < 0.001), and stronger estimates of cortical thickness reduction in females relative to males in cases than controls (corrected p-values < 0.001). Our results support the use of ComBat-GAM to minimize confounds and increase statistical power when harmonizing data with non-linear effects, and the use of either ComBat or ComBat-GAM for harmonizing data with linear effects.

Cortical volume abnormalities in posttraumatic stress disorder: an ENIGMA-psychiatric genomics consortium PTSD workgroup mega-analysis.

Studies of posttraumatic stress disorder (PTSD) report volume abnormalities in multiple regions of the cerebral cortex. However, findings for many regions, particularly regions outside commonly studied emotion-related prefrontal, insular, and limbic regions, are inconsistent and tentative. Also, few studies address the possibility that PTSD abnormalities may be confounded by comorbid depression. A mega-analysis investigating all cortical regions in a large sample of PTSD and control subjects can potentially provide new insight into these issues. Given this perspective, our group aggregated regional volumes data of 68 cortical regions across both hemispheres from 1379 PTSD patients to 2192 controls without PTSD after data were processed by 32 international laboratories using ENIGMA standardized procedures. We examined whether regional cortical volumes were different in PTSD vs. controls, were associated with posttraumatic stress symptom (PTSS) severity, or were affected by comorbid depression. Volumes of left and right lateral orbitofrontal gyri (LOFG), left superior temporal gyrus, and right insular, lingual and superior parietal gyri were significantly smaller, on average, in PTSD patients than controls (standardized coefficients = -0.111 to -0.068, FDR corrected P values < 0.039) and were significantly negatively correlated with PTSS severity. After adjusting for depression symptoms, the PTSD findings in left and right LOFG remained significant. These findings indicate that cortical volumes in PTSD patients are smaller in prefrontal regulatory regions, as well as in broader emotion and sensory processing cortical regions.

Individual-Specific Areal-Level Parcellations Improve Functional Connectivity Prediction of Behavior.

Resting-state functional magnetic resonance imaging (rs-fMRI) allows estimation of individual-specific cortical parcellations. We have previously developed a multi-session hierarchical Bayesian model (MS-HBM) for estimating high-quality individual-specific network-level parcellations. Here, we extend the model to estimate individual-specific areal-level parcellations. While network-level parcellations comprise spatially distributed networks spanning the cortex, the consensus is that areal-level parcels should be spatially localized, that is, should not span multiple lobes. There is disagreement about whether areal-level parcels should be strictly contiguous or comprise multiple noncontiguous components; therefore, we considered three areal-level MS-HBM variants spanning these range of possibilities. Individual-specific MS-HBM parcellations estimated using 10 min of data generalized better than other approaches using 150 min of data to out-of-sample rs-fMRI and task-fMRI from the same individuals. Resting-state functional connectivity derived from MS-HBM parcellations also achieved the best behavioral prediction performance. Among the three MS-HBM variants, the strictly contiguous MS-HBM exhibited the best resting-state homogeneity and most uniform within-parcel task activation. In terms of behavioral prediction, the gradient-infused MS-HBM was numerically the best, but differences among MS-HBM variants were not statistically significant. Overall, these results suggest that areal-level MS-HBMs can capture behaviorally meaningful individual-specific parcellation features beyond group-level parcellations. Multi-resolution trained models and parcellations are publicly available (https://github.com/ThomasYeoLab/CBIG/tree/master/stable_projects/brain_parcellation/Kong2022_ArealMSHBM).

Trait-like variants in human functional brain networks.

Resting-state functional magnetic resonance imaging (fMRI) has provided converging descriptions of group-level functional brain organization. Recent work has revealed that functional networks identified in individuals contain local features that differ from the group-level description. We define these features as network variants. Building on these studies, we ask whether distributions of network variants reflect stable, trait-like differences in brain organization. Across several datasets of highly-sampled individuals we show that 1) variants are highly stable within individuals, 2) variants are found in characteristic locations and associate with characteristic functional networks across large groups, 3) task-evoked signals in variants demonstrate a link to functional variation, and 4) individuals cluster into subgroups on the basis of variant characteristics that are related to differences in behavior. These results suggest that distributions of network variants may reflect stable, trait-like, functionally relevant individual differences in functional brain organization.

Precision Inhibitory Stimulation of Individual-Specific Cortical Hubs Disrupts Information Processing in Humans.

Noninvasive brain stimulation (NIBS) is a promising treatment for psychiatric and neurologic conditions, but outcomes are variable across treated individuals. In principle, precise targeting of individual-specific features of functional brain networks could improve the efficacy of NIBS interventions. Network theory predicts that the role of a node in a network can be inferred from its connections; as such, we hypothesized that targeting individual-specific "hub" brain areas with NIBS should impact cognition more than nonhub brain areas. Here, we first demonstrate that the spatial positioning of hubs is variable across individuals but reproducible within individuals upon repeated imaging. We then tested our hypothesis in healthy individuals using a prospective, within-subject, double-blind design. Inhibition of a hub with continuous theta burst stimulation disrupted information processing during working-memory more than inhibition of a nonhub area, despite targets being separated by only a few centimeters on the right middle frontal gyrus of each subject. Based upon these findings, we conclude that individual-specific brain network features are functionally relevant and could leveraged as stimulation sites in future NIBS interventions.