Precision Functional Mapping of Individual Human Brains

TS-AI: A deep learning pipeline for multimodal subject-specific parcellation with task contrasts synthesis

Accurate mapping of brain functional subregions at an individual level is crucial. Task-based functional MRI (tfMRI) captures subject-specific activation patterns during various functions and behaviors, facilitating the individual localization of functionally distinct subregions. However, acquiring high-quality tfMRI is time-consuming and resource-intensive in both scientific and clinical settings. The present study proposes a two-stage network model, TS-AI, to individualize an atlas on cortical surfaces through the prediction of tfMRI data. TS-AI first synthesizes a battery of task contrast maps for each individual by leveraging tract-wise anatomical connectivity and resting-state networks. These synthesized maps, along with feature maps of tract-wise anatomical connectivity and resting-state networks, are then fed into an end-to-end deep neural network to individualize an atlas. TS-AI enables the synthesized task contrast maps to be used in individual parcellation without the acquisition of actual task fMRI scans. In addition, a novel feature consistency loss is designed to assign vertices with similar features to the same parcel, which increases individual specificity and mitigates overfitting risks caused by the absence of individual parcellation ground truth. The individualized parcellations were validated by assessing test-retest reliability, homogeneity, and cognitive behavior prediction using diverse reference atlases and datasets, demonstrating the superior performance and generalizability of TS-AI. Sensitivity analysis yielded insights into region-specific features influencing individual variation in functional regionalization. Additionally, TS-AI identified accelerated shrinkage in the medial temporal and cingulate parcels during the progression of Alzheimer's disease, suggesting its potential in clinical research and applications.

Integration of estimated regional gene expression with neuroimaging and clinical phenotypes at biobank scale

An understanding of human brain individuality requires the integration of data on brain organization across people and brain regions, molecular and systems scales, as well as healthy and clinical states. Here, we help advance this understanding by leveraging methods from computational genomics to integrate large-scale genomic, transcriptomic, neuroimaging, and electronic-health record data sets. We estimated genetically regulated gene expression (gr-expression) of 18,647 genes, across 10 cortical and subcortical regions of 45,549 people from the UK Biobank. First, we showed that patterns of estimated gr-expression reflect known genetic-ancestry relationships, regional identities, as well as inter-regional correlation structure of directly assayed gene expression. Second, we performed transcriptome-wide association studies (TWAS) to discover 1,065 associations between individual variation in gr-expression and gray-matter volumes across people and brain regions. We benchmarked these associations against results from genome-wide association studies (GWAS) of the same sample and found hundreds of novel associations relative to these GWAS. Third, we integrated our results with clinical associations of gr-expression from the Vanderbilt Biobank. This integration allowed us to link genes, via gr-expression, to neuroimaging and clinical phenotypes. Fourth, we identified associations of polygenic gr-expression with structural and functional MRI phenotypes in the Human Connectome Project (HCP), a small neuroimaging-genomic data set with high-quality functional imaging data. Finally, we showed that estimates of gr-expression and magnitudes of TWAS were generally replicable and that the p-values of TWAS were replicable in large samples. Collectively, our results provide a powerful new resource for integrating gr-expression with population genetics of brain organization and disease.

Inferotemporal face patches are histo-architectonically distinct

An interconnected group of cortical regions distributed across the primate inferotemporal cortex forms a network critical for face perception. Understanding the microarchitecture of this face network can refine mechanistic accounts of how individual areas function and interact to support visual perception. To address this, we acquire a unique dataset in macaque monkeys combining fMRI to localize face patches in vivo and then ex vivo histology to resolve their histo-architecture across cortical depths in the same individuals. Our findings reveal that face patches differ based on cytochrome oxidase (CO) and, to a lesser extent, myelin staining, with the middle lateral (ML) face patch exhibiting pronounced CO staining. Histo-architectonic differences are less pronounced when using probabilistic definitions of face patches, underscoring the importance of precision mapping integrating in vivo and ex vivo measurements in the same individuals. This study indicates that the macaque face patch network is composed of architectonically distinct components.

Leveraging ultra-high field (7T) MRI in psychiatric research

Non-invasive brain imaging has played a critical role in establishing our understanding of the neural properties that contribute to the emergence of psychiatric disorders. However, characterizing core neurobiological mechanisms of psychiatric symptomatology requires greater structural, functional, and neurochemical specificity than is typically obtainable with standard field strength MRI acquisitions (e.g., 3T). Ultra-high field (UHF) imaging at 7 Tesla (7T) provides the opportunity to identify neurobiological systems that confer risk, determine etiology, and characterize disease progression and treatment outcomes of major mental illnesses. Increases in scanner availability, regulatory approval, and sequence availability have made the application of UHF to clinical cohorts more feasible than ever before, yet the application of UHF approaches to the study of mental health remains nascent. In this technical review, we describe core neuroimaging methodologies which benefit from UHF acquisition, including high resolution structural and functional imaging, single (1H) and multi-nuclear (e.g., 31P) MR spectroscopy, and quantitative MR techniques for assessing brain tissue iron and myelin. We discuss advantages provided by 7T MRI, including higher signal- and contrast-to-noise ratio, enhanced spatial resolution, increased test-retest reliability, and molecular and neurochemical specificity, and how these have begun to uncover mechanisms of psychiatric disorders. Finally, we consider current limitations of UHF in its application to clinical cohorts, and point to ongoing work that aims to overcome technical hurdles through the continued development of UHF hardware, software, and protocols.

Bringing Imaging Biomarkers Into Clinical Reality in Psychiatry

Importance

Advancing precision psychiatry, where treatments are based on an individual's biology rather than solely their clinical presentation, requires attention to several key attributes for any candidate biomarker. These include test-retest reliability, sensitivity to relevant neurophysiology, cost-effectiveness, and scalability. Unfortunately, these issues have not been systematically addressed by biomarker development efforts that use common neuroimaging tools like magnetic resonance imaging (MRI) and electroencephalography (EEG). Here, the critical barriers that neuroimaging methods will need to overcome to achieve clinical relevance in the near to intermediate term are examined.

Observations

Reliability is often overlooked, which together with sensitivity to key aspects of neurophysiology and replicated predictive utility, favors EEG-based methods. The principal barrier for EEG has been the lack of large-scale data collection among multisite psychiatric consortia. By contrast, despite its high reliability, structural MRI has not demonstrated clinical utility in psychiatry, which may be due to its limited sensitivity to psychiatry-relevant neurophysiology. Given the prevalence of structural MRIs, establishment of a compelling clinical use case remains its principal barrier. By contrast, low reliability and difficulty in standardizing collection are the principal barriers for functional MRI, along with the need for demonstration that its superior spatial resolution over EEG and ability to directly image subcortical regions in fact provide unique clinical value. Often missing, moreover, is consideration of how these various scientific issues can be balanced against practical economic realities of psychiatric health care delivery today, for which embedding economic modeling into biomarker development efforts may help direct research efforts.

Conclusions and Relevance

EEG seems most ripe for near- to intermediate-term clinical impact, especially considering its scalability and cost-effectiveness. Recent efforts to broaden its collection, as well as development of low-cost turnkey systems, suggest a promising pathway by which neuroimaging can impact clinical care. Continued MRI research focused on its key barriers may hold promise for longer-horizon utility.

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.

Principles of cortical areas and their implications for neuroimaging

Cortical organization should constrain the study of how the brain performs behavior and cognition. A fundamental concept in cortical organization is that of arealization: that the cortex is parceled into discrete areas. In part one of this report, we review how non-human animal studies have illuminated principles of cortical arealization by revealing: (1) what defines a cortical area, (2) how cortical areas are formed, (3) how cortical areas interact with one another, and (4) what "computations" or "functions" areas perform. In part two, we discuss how these principles apply to neuroimaging research. In doing so, we highlight several examples where the commonly accepted interpretation of neuroimaging observations requires assumptions that violate the principles of arealization, including nonstationary areas that move on short time scales, large-scale gradients as organizing features, and cortical areas with singular functionality that perfectly map psychological constructs. Our belief is that principles of neurobiology should strongly guide the nature of computational explanations.

Neuroimaging Biomarkers in Addiction

As a neurobiological process, addiction involves pathological patterns of engagement with substances and a range of behaviors with a chronic and relapsing course. Neuroimaging technologies assess brain activity, structure, physiology, and metabolism at scales ranging from neurotransmitter receptors to large-scale brain networks, providing unique windows into the core neural processes implicated in substance use disorders. Identified aberrations in the neural substrates of reward and salience processing, response inhibition, interoception, and executive functions with neuroimaging can inform the development of pharmacological, neuromodulatory, and psychotherapeutic interventions to modulate the disordered neurobiology. Based on our systematic search, 409 protocols registered on ClinicalTrials.gov include the use of one or more neuroimaging paradigms as an outcome measure in addiction, with the majority (N=268) employing functional magnetic resonance imaging (fMRI), followed by positron emission tomography (PET) (N=71), electroencephalography (EEG) (N=50), structural magnetic resonance imaging (MRI) (N=35) and magnetic resonance spectroscopy (MRS) (N=35). Furthermore, in a PubMed systematic review, we identified 61 meta-analyses including 30 fMRI, 22 structural MRI, 8 EEG, 7 PET, and 3 MRS meta-analyses suggesting potential biomarkers in addictions. These studies can facilitate the development of a range of biomarkers that may prove useful in the arsenal of addiction treatments in the coming years. There is evidence that these markers of large-scale brain structure and activity may indicate vulnerability or separate disease subtypes, predict response to treatment, or provide objective measures of treatment response or recovery. Neuroimaging biomarkers can also suggest novel targets for interventions. Closed or open loop interventions can integrate these biomarkers with neuromodulation in real-time or offline to personalize stimulation parameters and deliver the precise intervention. This review provides an overview of neuroimaging modalities in addiction, potential neuroimaging biomarkers, and their physiologic and clinical relevance. Future directions and challenges in bringing these putative biomarkers from the bench to the bedside are also discussed.

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.

The Promise and Challenges of Integrating Biological and Prevention Sciences: A Community-Engaged Model for the Next Generation of Translational Research

Beginning with the successful sequencing of the human genome two decades ago, the possibility of developing personalized health interventions based on one's biology has captured the imagination of researchers, medical providers, and individuals seeking health care services. However, the application of a personalized medicine approach to emotional and behavioral health has lagged behind the development of personalized approaches for physical health conditions. There is potential value in developing improved methods for integrating biological science with prevention science to identify risk and protective mechanisms that have biological underpinnings, and then applying that knowledge to inform prevention and intervention services for emotional and behavioral health. This report represents the work of a task force appointed by the Board of the Society for Prevention Research to explore challenges and recommendations for the integration of biological and prevention sciences. We present the state of the science and barriers to progress in integrating the two approaches, followed by recommended strategies that would promote the responsible integration of biological and prevention sciences. Recommendations are grounded in Community-Based Participatory Research approaches, with the goal of centering equity in future research aimed at integrating the two disciplines to ultimately improve the well-being of those who have disproportionately experienced or are at risk for experiencing emotional and behavioral problems.

White matter microstructure and functional connectivity in the brains of infants with Turner syndrome

Turner syndrome, caused by complete or partial loss of an X-chromosome, is often accompanied by specific cognitive challenges. Magnetic resonance imaging studies of adults and children with Turner syndrome suggest these deficits reflect differences in anatomical and functional connectivity. However, no imaging studies have explored connectivity in infants with Turner syndrome. Consequently, it is unclear when in development connectivity differences emerge. To address this gap, we compared functional connectivity and white matter microstructure of 1-year-old infants with Turner syndrome to typically developing 1-year-old boys and girls. We examined functional connectivity between the right precentral gyrus and five regions that show reduced volume in 1-year old infants with Turner syndrome compared to controls and found no differences. However, exploratory analyses suggested infants with Turner syndrome have altered connectivity between right supramarginal gyrus and left insula and right putamen. To assess anatomical connectivity, we examined diffusivity indices along the superior longitudinal fasciculus and found no differences. However, an exploratory analysis of 46 additional white matter tracts revealed significant group differences in nine tracts. Results suggest that the first year of life is a window in which interventions might prevent connectivity differences observed at later ages, and by extension, some of the cognitive challenges associated with Turner syndrome.

Multimodal Connectivity-Based Individual Parcellation and Analysis for Humans and Rhesus Monkeys

Individual brains vary greatly in morphology, connectivity and organization. Individualized brain parcellation is capable of precisely localizing subject-specific functional regions. However, most individualization approaches have examined single modalities of data and have not generalized to nonhuman primates. The present study proposed a novel multimodal connectivity-based individual parcellation (MCIP) method, which optimizes within-region homogeneity, spatial continuity and similarity to a reference atlas with the fusion of personal functional and anatomical connectivity. Comprehensive evaluation demonstrated that MCIP outperformed state-of-the-art multimodal individualization methods in terms of functional and anatomical homogeneity, predictability of cognitive measures, heritability, reproducibility and generalizability across species. Comparative investigation showed a higher topographic variability in humans than that in macaques. Therefore, MCIP provides improved accurate and reliable mapping of brain functional regions over existing methods at an individual level across species, and could facilitate comparative and translational neuroscience research.

Replicability and generalizability in population psychiatric neuroimaging

Studies linking mental health with brain function in cross-sectional population-based association studies have historically relied on small, underpowered samples. Given the small effect sizes typical of such brain-wide associations, studies require samples into the thousands to achieve the statistical power necessary for replicability. Here, we detail how small sample sizes have hampered replicability and provide sample size targets given established association strength benchmarks. Critically, while replicability will improve with larger samples, it is not guaranteed that observed effects will meaningfully apply to target populations of interest (i.e., be generalizable). We discuss important considerations related to generalizability in psychiatric neuroimaging and provide an example of generalizability failure due to "shortcut learning" in brain-based predictions of mental health phenotypes. Shortcut learning is a phenomenon whereby machine learning models learn an association between the brain and an unmeasured construct (the shortcut), rather than the intended target of mental health. Given the complex nature of brain-behavior interactions, the future of epidemiological approaches to brain-based studies of mental health will require large, diverse samples with comprehensive assessment.

Ventral attention network connectivity is linked to cortical maturation and cognitive ability in childhood

The human brain experiences functional changes through childhood and adolescence, shifting from an organizational framework anchored within sensorimotor and visual regions into one that is balanced through interactions with later-maturing aspects of association cortex. Here, we link this profile of functional reorganization to the development of ventral attention network connectivity across independent datasets. We demonstrate that maturational changes in cortical organization link preferentially to within-network connectivity and heightened degree centrality in the ventral attention network, whereas connectivity within network-linked vertices predicts cognitive ability. This connectivity is associated closely with maturational refinement of cortical organization. Children with low ventral attention network connectivity exhibit adolescent-like topographical profiles, suggesting that attentional systems may be relevant in understanding how brain functions are refined across development. These data suggest a role for attention networks in supporting age-dependent shifts in cortical organization and cognition across childhood and adolescence.

Precision Network Modeling of TMS Across Individuals Suggests Therapeutic Targets and Potential for Improvement

Higher-order cognitive and affective functions are supported by large-scale networks in the brain. Dysfunction in different networks is proposed to associate with distinct symptoms in neuro-psychiatric disorders. However, the specific networks targeted by current clinical transcranial magnetic stimulation (TMS) approaches are unclear. While standard-of-care TMS relies on scalp-based landmarks, recent FDA-approved TMS protocols use individualized functional connectivity with the subgenual anterior cingulate cortex (sgACC) to optimize TMS targeting. Leveraging previous work on precision network estimation and recent advances in network-level TMS targeting, we demonstrate that clinical TMS approaches target different functional networks between individuals. Homotopic scalp positions (left F3 and right F4) target different networks within and across individuals, and right F4 generally favors a right-lateralized control network. We also modeled the impact of targeting the dorsolateral prefrontal cortex (dlPFC) zone anticorrelated with the sgACC and found that the individual-specific anticorrelated region variably targets a network coupled to reward circuitry. Combining individualized, precision network mapping and electric field (E-field) modeling, we further illustrate how modeling can be deployed to prospectively target distinct closely localized association networks in the dlPFC with meaningful spatial selectivity and E-field intensity. Lastly, our findings emphasize differences between selectivity and maximal intensity, highlighting the need to consider both metrics in precision TMS efforts.

Connectome-based fingerprinting: reproducibility, precision, and behavioral prediction

Functional magnetic resonance imaging-based functional connectivity enables the non-invasive mapping of individual differences in brain functional organization to individual differences in a vast array of behavioral phenotypes. This flexibility has renewed the search for neuroimaging-based biomarkers that exhibit reproducibility, prediction, and precision. Functional connectivity-based measures that share these three characteristics are key to achieving this goal. Here, we review the functional connectome fingerprinting approach and discuss its value, not only as a simple and intuitive conceptualization of the "functional connectome" that provides new insights into how the connectome is altered in association with psychiatric symptoms, but also as a straightforward and interpretable method for indexing the reproducibility of functional connectivity-based measures. We discuss how these advantages provide new avenues for strengthening reproducibility, precision, and behavioral prediction for functional connectomics and we consider new directions toward discovering better biomarkers for neuropsychiatric conditions.

Complexity organization of resting-state functional-MRI networks

Entropy measures are increasingly being used to analyze the structure of neural activity observed by functional magnetic resonance imaging (fMRI), with resting-state networks (RSNs) being of interest for their reproducible descriptions of the brain's functional architecture. Temporal correlations have shown a dichotomy among these networks: those that engage with the environment, known as extrinsic, which include the visual and sensorimotor networks; and those associated with executive control and self-referencing, known as intrinsic, which include the default mode network and the frontoparietal control network. While these inter-voxel temporal correlations enable the assessment of synchrony among the components of individual networks, entropic measures introduce an intra-voxel assessment that quantifies signal features encoded within each blood oxygen level-dependent (BOLD) time series. As a result, this framework offers insights into comprehending the representation and processing of information within fMRI signals. Multiscale entropy (MSE) has been proposed as a useful measure for characterizing the entropy of neural activity across different temporal scales. This measure of temporal entropy in BOLD data is dependent on the length of the time series; thus, high-quality data with fine-grained temporal resolution and a sufficient number of time frames is needed to improve entropy precision. We apply MSE to the Midnight Scan Club, a highly sampled and well-characterized publicly available dataset, to analyze the entropy distribution of RSNs and evaluate its ability to distinguish between different functional networks. Entropy profiles are compared across temporal scales and RSNs. Our results have shown that the spatial distribution of entropy at infra-slow frequencies (0.005-0.1 Hz) reproduces known parcellations of RSNs. We found a complexity hierarchy between intrinsic and extrinsic RSNs, with intrinsic networks robustly exhibiting higher entropy than extrinsic networks. Finally, we found new evidence that the topography of entropy in the posterior cerebellum exhibits high levels of entropy comparable to that of intrinsic RSNs.

Psychiatric neuroimaging designs for individualised, cohort, and population studies

Psychiatric neuroimaging faces challenges to rigour and reproducibility that prompt reconsideration of the relative strengths and limitations of study designs. Owing to high resource demands and varying inferential goals, current designs differentially emphasise sample size, measurement breadth, and longitudinal assessments. In this overview and perspective, we provide a guide to the current landscape of psychiatric neuroimaging study designs with respect to this balance of scientific goals and resource constraints. Through a heuristic data cube contrasting key design features, we discuss a resulting trade-off among small sample, precision longitudinal studies (e.g., individualised studies and cohorts) and large sample, minimally longitudinal, population studies. Precision studies support tests of within-person mechanisms, via intervention and tracking of longitudinal course. Population studies support tests of generalisation across multifaceted individual differences. A proposed reciprocal validation model (RVM) aims to recursively leverage these complementary designs in sequence to accumulate evidence, optimise relative strengths, and build towards improved long-term clinical utility.

Current best practices and future opportunities for reproducible findings using large-scale neuroimaging in psychiatry

Research into the brain basis of psychopathology is challenging due to the heterogeneity of psychiatric disorders, extensive comorbidities, underdiagnosis or overdiagnosis, multifaceted interactions with genetics and life experiences, and the highly multivariate nature of neural correlates. Therefore, increasingly larger datasets that measure more variables in larger cohorts are needed to gain insights. In this review, we present current "best practice" approaches for using existing databases, collecting and sharing new repositories for big data analyses, and future directions for big data in neuroimaging and psychiatry with an emphasis on contributing to collaborative efforts and the challenges of multi-study data analysis.

Individual patterns of functional connectivity in neonates as revealed by surface-based Bayesian modeling

Resting-state functional connectivity is a widely used approach to study the functional brain network organization during early brain development. However, the estimation of functional connectivity networks in individual infants has been rather elusive due to the unique challenges involved with functional magnetic resonance imaging (fMRI) data from young populations. Here, we use fMRI data from the developing Human Connectome Project (dHCP) database to characterize individual variability in a large cohort of term-born infants (N = 289) using a novel data-driven Bayesian framework. To enhance alignment across individuals, the analysis was conducted exclusively on the cortical surface, employing surface-based registration guided by age-matched neonatal atlases. Using 10 minutes of resting-state fMRI data, we successfully estimated subject-level maps for fourteen brain networks/subnetworks along with individual functional parcellation maps that revealed differences between subjects. We also found a significant relationship between age and mean connectivity strength in all brain regions, including previously unreported findings in higher-order networks. These results illustrate the advantages of surface-based methods and Bayesian statistical approaches in uncovering individual variability within very young populations.

Measuring neuroplasticity in human development: the potential to inform the type and timing of mental health interventions

Neuroplasticity during sensitive periods, the molecular and cellular process of enduring neural change in response to external stimuli during windows of high environmental sensitivity, is crucial for adaptation to expected environments and has implications for psychiatry. Animal research has characterized the developmental sequence and neurobiological mechanisms that govern neuroplasticity, yet gaps in our ability to measure neuroplasticity in humans limit the clinical translation of these principles. Here, we present a roadmap for the development and validation of neuroimaging and electrophysiology measures that index neuroplasticity to begin to address these gaps. We argue that validation of measures to track neuroplasticity in humans will elucidate the etiology of mental illness and inform the type and timing of mental health interventions to optimize effectiveness. We outline criteria for evaluating putative neuroimaging measures of plasticity in humans including links to neurobiological mechanisms shown to govern plasticity in animal models, developmental change that reflects heightened early life plasticity, and prediction of neural and/or behavior change. These criteria are applied to three putative measures of neuroplasticity using electroencephalography (gamma oscillations, aperiodic exponent of power/frequency) or functional magnetic resonance imaging (amplitude of low frequency fluctuations). We discuss the use of these markers in psychiatry, envision future uses for clinical and developmental translation, and suggest steps to address the limitations of the current putative neuroimaging measures of plasticity. With additional work, we expect these markers will significantly impact mental health and be used to characterize mechanisms, devise new interventions, and optimize developmental trajectories to reduce psychopathology risk.

Moving beyond processing- and analysis-related variation in resting-state functional brain imaging

When fields lack consensus standard methods and accessible ground truths, reproducibility can be more of an ideal than a reality. Such has been the case for functional neuroimaging, where there exists a sprawling space of tools and processing pipelines. We provide a critical evaluation of the impact of differences across five independently developed minimal preprocessing pipelines for functional magnetic resonance imaging. We show that, even when handling identical data, interpipeline agreement was only moderate, critically shedding light on a factor that limits cross-study reproducibility. We show that low interpipeline agreement can go unrecognized until the reliability of the underlying data is high, which is increasingly the case as the field progresses. Crucially we show that, when interpipeline agreement is compromised, so too is the consistency of insights from brain-wide association studies. We highlight the importance of comparing analytic configurations, because both widely discussed and commonly overlooked decisions can lead to marked variation.

Reliability of energy landscape analysis of resting-state functional MRI data

Energy landscape analysis is a data-driven method to analyse multidimensional time series, including functional magnetic resonance imaging (fMRI) data. It has been shown to be a useful characterization of fMRI data in health and disease. It fits an Ising model to the data and captures the dynamics of the data as movement of a noisy ball constrained on the energy landscape derived from the estimated Ising model. In the present study, we examine test-retest reliability of the energy landscape analysis. To this end, we construct a permutation test that assesses whether or not indices characterizing the energy landscape are more consistent across different sets of scanning sessions from the same participant (i.e. within-participant reliability) than across different sets of sessions from different participants (i.e. between-participant reliability). We show that the energy landscape analysis has significantly higher within-participant than between-participant test-retest reliability with respect to four commonly used indices. We also show that a variational Bayesian method, which enables us to estimate energy landscapes tailored to each participant, displays comparable test-retest reliability to that using the conventional likelihood maximization method. The proposed methodology paves the way to perform individual-level energy landscape analysis for given data sets with a statistically controlled reliability.

Longitudinal associations between language network characteristics in the infant brain and school-age reading abilities are mediated by early-developing phonological skills

Reading acquisition is a prolonged learning process relying on language development starting in utero. Behavioral longitudinal studies reveal prospective associations between infant language abilities and preschool/kindergarten phonological development that relates to subsequent reading performance. While recent pediatric neuroimaging work has begun to characterize the neural network underlying language development in infants, how this neural network scaffolds long-term language and reading acquisition remains unknown. We addressed this question in a 7-year longitudinal study from infancy to school-age. Seventy-six infants completed resting-state fMRI scanning, and underwent standardized language assessments in kindergarten. Of this larger cohort, forty-one were further assessed on their emergent word reading abilities after receiving formal reading instructions. Hierarchical clustering analyses identified a modular infant language network in which functional connectivity (FC) of the inferior frontal module prospectively correlated with kindergarten-age phonological skills and emergent word reading abilities. These correlations were obtained when controlling for infant age at scan, nonverbal IQ and parental education. Furthermore, kindergarten-age phonological skills mediated the relationship between infant FC and school-age reading abilities, implying a critical mid-way milestone for long-term reading development from infancy. Overall, our findings illuminate the neurobiological mechanisms by which infant language capacities could scaffold long-term reading acquisition.

Brain-wide functional connectivity artifactually inflates throughout functional magnetic resonance imaging scans

Functional magnetic resonance imaging (fMRI) is a central tool for investigating human brain function, organization and disease. Here, we show that fMRI-based estimates of functional brain connectivity artifactually inflate at spatially heterogeneous rates during resting-state and task-based scans. This produces false positive connection strength changes and spatial distortion of brain connectivity maps. We demonstrate that this artefact is driven by temporal inflation of the non-neuronal, systemic low-frequency oscillation (sLFO) blood flow signal during fMRI scanning and is not addressed by standard denoising procedures. We provide evidence that sLFO inflation reflects perturbations in cerebral blood flow by respiration and heart rate changes that accompany diminishing arousal during scanning, although the mechanisms of this pathway are uncertain. Finally, we show that adding a specialized sLFO denoising procedure to fMRI processing pipelines mitigates the artifactual inflation of functional connectivity, enhancing the validity and within-scan reliability of fMRI findings.

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.

Macroscale brain states support the control of semantic cognition

A crucial aim in neuroscience is to understand how the human brain adapts to varying cognitive demands. This study investigates network reconfiguration during controlled semantic retrieval in differing contexts. We analyze brain responses to two semantic tasks of varying difficulty - global association and feature matching judgments - which are contrasted with non-semantic tasks on the cortical surface and within a whole-brain state space. Demanding semantic association tasks elicit activation in anterior prefrontal and temporal regions, while challenging semantic feature matching and non-semantic tasks predominantly activate posterior regions. Task difficulty also modulates activation along different dimensions of functional organization, suggesting different mechanisms of cognitive control. More demanding semantic association judgments engage cognitive control and default mode networks together, while feature matching and non-semantic tasks are skewed towards cognitive control networks. These findings highlight the brain's dynamic ability to tailor its networks to support diverse neurocognitive states, enriching our understanding of controlled cognition.

What ergodicity means for you

This paper explores the relation between within-person and between-person research designs using the concept of ergodicity from statistical mechanics in physics. We demonstrate the consequences of ergodicity using several real data examples from previously published studies. We then create several simulated examples that illustrate the independence of within-person processes from between-person differences, and pair these examples with analytic results that reinforce our conclusions. Finally, we discuss the plausibility of ergodicity being the general rule rather than the exception for social and behavioral processes, address common arguments against heeding the implications of ergodicity for behavioral research, and offer several possible solutions.

The promise of precision functional mapping for neuroimaging in psychiatry

Precision functional mapping (PFM) is a neuroimaging approach to reliably estimate metrics of brain function from individual people via the collection of large amounts of fMRI data (hours per person). This method has revealed much about the inter-individual variation of functional brain networks. While standard group-level studies, in which we average brain measures across groups of people, are important in understanding the generalizable neural underpinnings of neuropsychiatric disorders, many disorders are heterogeneous in nature. This heterogeneity often complicates clinical care, leading to patient uncertainty when considering prognosis or treatment options. We posit that PFM methods may help streamline clinical care in the future, fast-tracking the choice of personalized treatment that is most compatible with the individual. In this review, we provide a history of PFM studies, foundational results highlighting the benefits of PFM methods in the pursuit of an advanced understanding of individual differences in functional network organization, and possible avenues where PFM can contribute to clinical translation of neuroimaging research results in the way of personalized treatment in psychiatry.

A subset of brain regions within adult functional connectivity networks demonstrate high reliability across early development

The human cerebral cortex contains groups of areas that support sensory, motor, cognitive, and affective functions, often categorized as functional networks. These areas show stronger internal and weaker external functional connectivity (FC) and exhibit similar FC profiles within rather than between networks. Previous studies have demonstrated the development of these networks from nascent forms present before birth to their mature, adult-like topography in childhood. However, analyses often still use definitions based on adult functional networks. We aim to assess how this might lead to the misidentification of functional networks and explore potential consequences and solutions. Our findings suggest that even though adult networks provide only a marginally better than-chance description of the infant FC organization, misidentification was largely driven by specific areas. By restricting functional networks to areas showing adult-like network clustering, we observed consistent within-network FC both within and across scans and throughout development. Additionally, these areas were spatially closer to locations with low variability in network identity among adults. Our analysis aids in understanding the potential consequences of using adult networks "as is" and provides guidance for future research on selecting and utilizing functional network models based on the research question and scenario.

Measuring functional connectivity in frequency-domain helps to better characterize brain function

Resting-state functional connectivity (FC) is widely used in multivariate pattern analysis of functional magnetic resonance imaging (fMRI), including identifying the locations of putative brain functional borders, predicting individual phenotypes, and diagnosing clinical mental diseases. However, limited attention has been paid to the analysis of functional interactions from a frequency perspective. In this study, by contrasting coherence-based and correlation-based FC with two machine learning tasks, we observed that measuring FC in the frequency domain helped to identify finer functional subregions and achieve better pattern discrimination capability relative to the temporal correlation. This study has proven the feasibility of coherence in the analysis of fMRI, and the results indicate that modeling functional interactions in the frequency domain may provide richer information than that in the time domain, which may provide a new perspective on the analysis of functional neuroimaging.

Utilizing connectome fingerprinting functional MRI models for motor activity prediction in presurgical planning: A feasibility study

Presurgical planning prior to brain tumor resection is critical for the preservation of neurologic function post-operatively. Neurosurgeons increasingly use advanced brain mapping techniques pre- and intra-operatively to delineate brain regions which are "eloquent" and should be spared during resection. Functional MRI (fMRI) has emerged as a commonly used non-invasive modality for individual patient mapping of critical cortical regions such as motor, language, and visual cortices. To map motor function, patients are scanned using fMRI while they perform various motor tasks to identify brain networks critical for motor performance, but it may be difficult for some patients to perform tasks in the scanner due to pre-existing deficits. Connectome fingerprinting (CF) is a machine-learning approach that learns associations between resting-state functional networks of a brain region and the activations in the region for specific tasks; once a CF model is constructed, individualized predictions of task activation can be generated from resting-state data. Here we utilized CF to train models on high-quality data from 208 subjects in the Human Connectome Project (HCP) and used this to predict task activations in our cohort of healthy control subjects (n = 15) and presurgical patients (n = 16) using resting-state fMRI (rs-fMRI) data. The prediction quality was validated with task fMRI data in the healthy controls and patients. We found that the task predictions for motor areas are on par with actual task activations in most healthy subjects (model accuracy around 90%-100% of task stability) and some patients suggesting the CF models can be reliably substituted where task data is either not possible to collect or hard for subjects to perform. We were also able to make robust predictions in cases in which there were no task-related activations elicited. The findings demonstrate the utility of the CF approach for predicting activations in out-of-sample subjects, across sites and scanners, and in patient populations. This work supports the feasibility of the application of CF models to presurgical planning, while also revealing challenges to be addressed in future developments. PRACTITIONER POINTS: Precision motor network prediction using connectome fingerprinting. Carefully trained models' performance limited by stability of task-fMRI data. Successful cross-scanner predictions and motor network mapping in patients with tumor.

Estimating the total variance explained by whole-brain imaging for zero-inflated outcomes

There is a dearth of statistical models that adequately capture the total signal attributed to whole-brain imaging features. The total signal is often widely distributed across the brain, with individual imaging features exhibiting small effect sizes for predicting neurobehavioral phenotypes. The challenge of capturing the total signal is compounded by the distribution of neurobehavioral data, particularly responses to psychological questionnaires, which often feature zero-inflated, highly skewed outcomes. To close this gap, we have developed a novel Variational Bayes algorithm that characterizes the total signal captured by whole-brain imaging features for zero-inflated outcomes. Our zero-inflated variance (ZIV) estimator estimates the fraction of variance explained (FVE) and the proportion of non-null effects (PNN) from large-scale imaging data. In simulations, ZIV demonstrates superior performance over other linear models. When applied to data from the Adolescent Brain Cognitive DevelopmentSM (ABCD) Study, we found that whole-brain imaging features contribute to a larger FVE for externalizing behaviors compared to internalizing behaviors. Moreover, focusing on features contributing to the PNN, ZIV estimator localized key neurocircuitry associated with neurobehavioral traits. To the best of our knowledge, the ZIV estimator is the first specialized method for analyzing zero-inflated neuroimaging data, enhancing future studies on brain-behavior relationships and improving the understanding of neurobehavioral disorders.

Transdiagnostic Neurobiology of Social Cognition and Individual Variability as Measured by Fractional Amplitude of Low-Frequency Fluctuation in Schizophrenia and Autism Spectrum Disorders

Fractional amplitude of low-frequency fluctuation (fALFF) is a validated measure of resting-state spontaneous brain activity. Previous fALFF findings in autism and schizophrenia spectrum disorders (ASDs and SSDs) have been highly heterogeneous. We aimed to use fALFF in a large sample of typically developing control (TDC), ASD and SSD participants to explore group differences and relationships with inter-individual variability of fALFF maps and social cognition. fALFF from 495 participants (185 TDC, 68 ASD, and 242 SSD) was computed using functional magnetic resonance imaging as signal power within two frequency bands (i.e., slow-4 and slow-5), normalized by the power in the remaining frequency spectrum. Permutation analysis of linear models was employed to investigate the relationship of fALFF with diagnostic groups, higher-level social cognition, and lower-level social cognition. Each participant's average distance of fALFF map to all others was defined as a variability score, with higher scores indicating less typical maps. Lower fALFF in the visual and higher fALFF in the frontal regions were found in both SSD and ASD participants compared with TDCs. Limited differences were observed between ASD and SSD participants in the cuneus regions only. Associations between slow-4 fALFF and higher-level social cognitive scores across the whole sample were observed in the lateral occipitotemporal and temporoparietal junction. Individual variability within the ASD and SSD groups was also significantly higher compared with TDC. Similar patterns of fALFF and individual variability in ASD and SSD suggest some common neurobiological deficits across these related heterogeneous conditions.

Frequency-specific directed connectivity between the hippocampus and parietal cortex during verbal and spatial episodic memory: an intracranial EEG replication

Hippocampus-parietal cortex circuits are thought to play a crucial role in memory and attention, but their neural basis remains poorly understood. We employed intracranial intracranial electroencephalography (iEEG) to investigate the neurophysiological underpinning of these circuits across three memory tasks spanning verbal and spatial domains. We uncovered a consistent pattern of higher causal directed connectivity from the hippocampus to both lateral parietal cortex (supramarginal and angular gyrus) and medial parietal cortex (posterior cingulate cortex) in the delta-theta band during memory encoding and recall. This connectivity was independent of activation or suppression states in the hippocampus or parietal cortex. Crucially, directed connectivity from the supramarginal gyrus to the hippocampus was enhanced in participants with higher memory recall, highlighting its behavioral significance. Our findings align with the attention-to-memory model, which posits that attention directs cognitive resources toward pertinent information during memory formation. The robustness of these results was demonstrated through Bayesian replication analysis of the memory encoding and recall periods across the three tasks. Our study sheds light on the neural basis of casual signaling within hippocampus-parietal circuits, broadening our understanding of their critical roles in human cognition.

Improved fMRI-based pain prediction using Bayesian group-wise functional registration

In recent years, the field of neuroimaging has undergone a paradigm shift, moving away from the traditional brain mapping approach towards the development of integrated, multivariate brain models that can predict categories of mental events. However, large interindividual differences in both brain anatomy and functional localization after standard anatomical alignment remain a major limitation in performing this type of analysis, as it leads to feature misalignment across subjects in subsequent predictive models. This article addresses this problem by developing and validating a new computational technique for reducing misalignment across individuals in functional brain systems by spatially transforming each subject's functional data to a common latent template map. Our proposed Bayesian functional group-wise registration approach allows us to assess differences in brain function across subjects and individual differences in activation topology. We achieve the probabilistic registration with inverse-consistency by utilizing the generalized Bayes framework with a loss function for the symmetric group-wise registration. It models the latent template with a Gaussian process, which helps capture spatial features in the template, producing a more precise estimation. We evaluate the method in simulation studies and apply it to data from an fMRI study of thermal pain, with the goal of using functional brain activity to predict physical pain. We find that the proposed approach allows for improved prediction of reported pain scores over conventional approaches. Received on 2 January 2017. Editorial decision on 8 June 2021.

Development of human visual cortical function: A scoping review of task- and naturalistic-fMRI studies through the interactive specialization and maturational frameworks

Overarching theories such as the interactive specialization and maturational frameworks have been proposed to describe human functional brain development. However, these frameworks have not yet been systematically examined across the fMRI literature. Visual processing is one of the most well-studied fields in neuroimaging, and research in this area has recently expanded to include naturalistic paradigms that facilitate study in younger age ranges, allowing for an in-depth critical appraisal of these frameworks across childhood. To this end, we conducted a scoping review of 94 developmental visual fMRI studies, including both traditional experimental task and naturalistic studies, across multiple sub-domains (early visual processing, category-specific higher order processing, naturalistic visual processing). We found that across domains, many studies reported progressive development, but few studies describe regressive or emergent changes necessary to fit the maturational or interactive specialization frameworks. Our findings suggest a need for the expansion of developmental frameworks and clearer reporting of both progressive and regressive changes, along with well-powered, longitudinal studies.

Task structure tailors the geometry of neural representations in human lateral prefrontal cortex

How do human brains represent tasks of varying structure? The lateral prefrontal cortex (lPFC) flexibly represents task information. However, principles that shape lPFC representational geometry remain unsettled. We use fMRI and pattern analyses to reveal the structure of lPFC representational geometries as humans perform two distinct categorization tasks- one with flat, conjunctive categories and another with hierarchical, context-dependent categories. We show that lPFC encodes task-relevant information with task-tailored geometries of intermediate dimensionality. These geometries preferentially enhance the separability of task-relevant variables while encoding a subset in abstract form. Specifically, in the flat task, a global axis encodes response-relevant categories abstractly, while category-specific local geometries are high-dimensional. In the hierarchy task, a global axis abstractly encodes the higher-level context, while low-dimensional, context-specific local geometries compress irrelevant information and abstractly encode the relevant information. Comparing these task geometries exposes generalizable principles by which lPFC tailors representations to different tasks.

Reconfigurations of cortical manifold structure during reward-based motor learning

Adaptive motor behavior depends on the coordinated activity of multiple neural systems distributed across the brain. While the role of sensorimotor cortex in motor learning has been well established, how higher-order brain systems interact with sensorimotor cortex to guide learning is less well understood. Using functional MRI, we examined human brain activity during a reward-based motor task where subjects learned to shape their hand trajectories through reinforcement feedback. We projected patterns of cortical and striatal functional connectivity onto a low-dimensional manifold space and examined how regions expanded and contracted along the manifold during learning. During early learning, we found that several sensorimotor areas in the dorsal attention network exhibited increased covariance with areas of the salience/ventral attention network and reduced covariance with areas of the default mode network (DMN). During late learning, these effects reversed, with sensorimotor areas now exhibiting increased covariance with DMN areas. However, areas in posteromedial cortex showed the opposite pattern across learning phases, with its connectivity suggesting a role in coordinating activity across different networks over time. Our results establish the neural changes that support reward-based motor learning and identify distinct transitions in the functional coupling of sensorimotor to transmodal cortex when adapting behavior.

Reducing individual differences in task fMRI with OGRE (One-step General Registration and Extraction) preprocessing

Volumetric analysis methods continue to enjoy great popularity in the analysis of task-related functional MRI (fMRI) data. Among these methods, the software package FSL (FMRIB, Oxford, UK) is omnipresent throughout the field. However, it remains unknown what advantages might be gained by integrating FSL with alternative preprocessing tools. Here we developed the One-step General Registration and Extraction (OGRE) pipeline to apply FreeSurfer brain extraction for simultaneous registration and motion correction ("one-step resampling"), for FSL volumetric analysis of fMRI data. We compared three preprocessing approaches (OGRE, FSL, and fMRIPrep) in a dataset wherein adult human volunteers (N=39) performed a precision drawing task during fMRI scanning. The three approaches produced grossly similar results, but OGRE's preprocessing led to lower inter-individual variability across the brain and greater detected activation in primary motor cortex contralateral to movement. This demonstrates that FreeSurfer tools and one-step resampling can improve FSL's volumetric analysis of fMRI data. The OGRE pipeline provides an off-the-shelf method to apply FreeSurfer-based brain extraction and one-step resampling of motion correction and registration for FSL analysis of task fMRI data.

Quality over quantity: powering neuroimaging samples in psychiatry

Neuroimaging has been widely adopted in psychiatric research, with hopes that these non-invasive methods will provide important clues to the underpinnings and prediction of various mental health symptoms and outcomes. However, the translational impact of neuroimaging has not yet reached its promise, despite the plethora of computational methods, tools, and datasets at our disposal. Some have lamented that too many psychiatric neuroimaging studies have been underpowered with respect to sample size. In this review, we encourage this discourse to shift from a focus on sheer increases in sample size to more thoughtful choices surrounding experimental study designs. We propose considerations at multiple decision points throughout the study design, data modeling and analysis process that may help researchers working in psychiatric neuroimaging boost power for their research questions of interest without necessarily increasing sample size. We also provide suggestions for leveraging multiple datasets to inform each other and strengthen our confidence in the generalization of findings to both population-level and clinical samples. Through a greater emphasis on improving the quality of brain-based and clinical measures rather than merely quantity, meaningful and potentially translational clinical associations with neuroimaging measures can be achieved with more modest sample sizes in psychiatry.

Structural connectivity matures along a sensorimotor-association connectional axis in youth

Childhood and adolescence are associated with protracted developmental remodeling of cortico-cortical structural connectivity. However, how heterochronous development in white matter structural connectivity spatially and temporally unfolds across the macroscale human connectome remains unknown. Leveraging non-invasive diffusion MRI data from both cross-sectional (N = 590) and longitudinal (baseline: N = 3,949; two-year follow-up: N = 3,155) developmental datasets, we found that structural connectivity development diverges along a pre-defined sensorimotor-association (S-A) connectional axis from ages 8.1 to 21.9 years. Specifically, we observed a continuum of developmental profiles that spans from an early childhood increase in connectivity strength in sensorimotor-sensorimotor connections to a late adolescent increase in association-association connectional strength. The S-A connectional axis also captured spatial variations in associations between structural connectivity and both higher-order cognition and general psychopathology. Together, our findings reveal a hierarchical axis in the development of structural connectivity across the human connectome.

Functional connectome through the human life span

The lifespan growth of the functional connectome remains unknown. Here, we assemble task-free functional and structural magnetic resonance imaging data from 33,250 individuals aged 32 postmenstrual weeks to 80 years from 132 global sites. We report critical inflection points in the nonlinear growth curves of the global mean and variance of the connectome, peaking in the late fourth and late third decades of life, respectively. After constructing a fine-grained, lifespan-wide suite of system-level brain atlases, we show distinct maturation timelines for functional segregation within different systems. Lifespan growth of regional connectivity is organized along a primary-to-association cortical axis. These connectome-based normative models reveal substantial individual heterogeneities in functional brain networks in patients with autism spectrum disorder, major depressive disorder, and Alzheimer's disease. These findings elucidate the lifespan evolution of the functional connectome and can serve as a normative reference for quantifying individual variation in development, aging, and neuropsychiatric disorders.

Longitudinal associations between language network characteristics in the infant brain and school-age reading abilities are mediated by early-developing phonological skills

Reading acquisition is a prolonged learning process relying on language development starting in utero. Behavioral longitudinal studies reveal prospective associations between infant language abilities and preschool/kindergarten phonological development that relates to subsequent reading performance. While recent pediatric neuroimaging work has begun to characterize the neural network underlying language development in infants, how this neural network scaffolds long-term language and reading acquisition remains unknown. We addressed this question in a 7-year longitudinal study from infancy to school-age. Seventy-six infants completed resting-state fMRI scanning, and underwent standardized language assessments in kindergarten. Of this larger cohort, forty-one were further assessed on their emergent word reading abilities after receiving formal reading instructions. Hierarchical clustering analyses identified a modular infant language network in which functional connectivity (FC) of the inferior frontal module prospectively correlated with kindergarten-age phonological skills and emergent word reading abilities. These correlations were obtained when controlling for infant age at scan, nonverbal IQ and parental education. Furthermore, kindergarten-age phonological skills mediated the relationship between infant FC and school-age reading abilities, implying a critical mid-way milestone for long-term reading development from infancy. Overall, our findings illuminate the neurobiological mechanisms by which infant language capacities could scaffold long-term reading acquisition.

Fast connectivity gradient approximation: maintaining spatially fine-grained connectivity gradients while reducing computational costs

Brain connectome analysis suffers from the high dimensionality of connectivity data, often forcing a reduced representation of the brain at a lower spatial resolution or parcellation. This is particularly true for graph-based representations, which are increasingly used to characterize connectivity gradients, capturing patterns of systematic spatial variation in the functional connectivity structure. However, maintaining a high spatial resolution is crucial for enabling fine-grained topographical analysis and preserving subtle individual differences that might otherwise be lost. Here we introduce a computationally efficient approach to establish spatially fine-grained connectivity gradients. At its core, it leverages a set of landmarks to approximate the underlying connectivity structure at the full spatial resolution without requiring a full-scale vertex-by-vertex connectivity matrix. We show that this approach reduces computational time and memory usage while preserving informative individual features and demonstrate its application in improving brain-behavior predictions. Overall, its efficiency can remove computational barriers and enable the widespread application of connectivity gradients to capture spatial signatures of the connectome. Importantly, maintaining a spatially fine-grained resolution facilitates to characterize the spatial transitions inherent in the core concept of gradients of brain organization.

Individual Brain Charting dataset extension, third release for movie watching and retinotopy data

The Individual Brain Charting (IBC) is a multi-task functional Magnetic Resonance Imaging dataset acquired at high spatial-resolution and dedicated to the cognitive mapping of the human brain. It consists in the deep phenotyping of twelve individuals, covering a broad range of psychological domains suitable for functional-atlasing applications. Here, we present the inclusion of task data from both naturalistic stimuli and trial-based designs, to uncover structures of brain activation. We rely on the Fast Shared Response Model (FastSRM) to provide a data-driven solution for modelling naturalistic stimuli, typically containing many features. We show that data from left-out runs can be reconstructed using FastSRM, enabling the extraction of networks from the visual, auditory and language systems. We also present the topographic organization of the visual system through retinotopy. In total, six new tasks were added to IBC, wherein four trial-based retinotopic tasks contributed with a mapping of the visual field to the cortex. IBC is open access: source plus derivatives imaging data and meta-data are available in public repositories.

Simultaneous cortical, subcortical, and brainstem mapping of sensory activation

Nonpainful tactile sensory stimuli are processed in the cortex, subcortex, and brainstem. Recent functional magnetic resonance imaging studies have highlighted the value of whole-brain, systems-level investigation for examining sensory processing. However, whole-brain functional magnetic resonance imaging studies are uncommon, in part due to challenges with signal to noise when studying the brainstem. Furthermore, differentiation of small sensory brainstem structures such as the cuneate and gracile nuclei necessitates high-resolution imaging. To address this gap in systems-level sensory investigation, we employed a whole-brain, multi-echo functional magnetic resonance imaging acquisition at 3T with multi-echo independent component analysis denoising and brainstem-specific modeling to enable detection of activation across the entire sensory system. In healthy participants, we examined patterns of activity in response to nonpainful brushing of the right hand, left hand, and right foot (n = 10 per location), and found the expected lateralization, with distinct cortical and subcortical responses for upper and lower limb stimulation. At the brainstem level, we differentiated the adjacent cuneate and gracile nuclei, corresponding to hand and foot stimulation respectively. Our findings demonstrate that simultaneous cortical, subcortical, and brainstem mapping at 3T could be a key tool to understand the sensory system in both healthy individuals and clinical cohorts with sensory deficits.

Neuroanatomical correlates of peripersonal space: bridging the gap between perception, action, emotion and social cognition

Peripersonal space (PPS) is a construct referring to the portion of space immediately surrounding our bodies, where most of the interactions between the subject and the environment, including other individuals, take place. Decades of animal and human neuroscience research have revealed that the brain holds a separate representation of this region of space: this distinct spatial representation has evolved to ensure proper relevance to stimuli that are close to the body and prompt an appropriate behavioral response. The neural underpinnings of such construct have been thoroughly investigated by different generations of studies involving anatomical and electrophysiological investigations in animal models, and, recently, neuroimaging experiments in human subjects. Here, we provide a comprehensive anatomical overview of the anatomical circuitry underlying PPS representation in the human brain. Gathering evidence from multiple areas of research, we identified cortical and subcortical regions that are involved in specific aspects of PPS encoding.We show how these regions are part of segregated, yet integrated functional networks within the brain, which are in turn involved in higher-order integration of information. This wide-scale circuitry accounts for the relevance of PPS encoding in multiple brain functions, including not only motor planning and visuospatial attention but also emotional and social cognitive aspects. A complete characterization of these circuits may clarify the derangements of PPS representation observed in different neurological and neuropsychiatric diseases.

Automated individual cortical parcellation via consensus graph representation learning

Cortical parcellation plays a pivotal role in elucidating the brain organization. Despite the growing efforts to develop parcellation algorithms using functional magnetic resonance imaging, achieving a balance between intra-individual specificity and inter-individual consistency proves challenging, making the generation of high-quality, subject-consistent cortical parcellations particularly elusive. To solve this problem, our paper proposes a fully automated individual cortical parcellation method based on consensus graph representation learning. The method integrates spectral embedding with low-rank tensor learning into a unified optimization model, which uses group-common connectivity patterns captured by low-rank tensor learning to optimize subjects' functional networks. This not only ensures consistency in brain representations across different subjects but also enhances the quality of each subject's representation matrix by eliminating spurious connections. More importantly, it achieves an adaptive balance between intra-individual specificity and inter-individual consistency during this process. Experiments conducted on a test-retest dataset from the Human Connectome Project (HCP) demonstrate that our method outperforms existing methods in terms of reproducibility, functional homogeneity, and alignment with task activation. Extensive network-based comparisons on the HCP S900 dataset reveal that the functional network derived from our cortical parcellation method exhibits greater capabilities in gender identification and behavior prediction than other approaches.