Brain-phenotype models fail for individuals who defy sample stereotypes.

Individual differences in brain functional organization track a range of traits, symptoms and behaviours1-12. So far, work modelling linear brain-phenotype relationships has assumed that a single such relationship generalizes across all individuals, but models do not work equally well in all participants13,14. A better understanding of in whom models fail and why is crucial to revealing robust, useful and unbiased brain-phenotype relationships. To this end, here we related brain activity to phenotype using predictive models-trained and tested on independent data to ensure generalizability15-and examined model failure. We applied this data-driven approach to a range of neurocognitive measures in a new, clinically and demographically heterogeneous dataset, with the results replicated in two independent, publicly available datasets16,17. Across all three datasets, we find that models reflect not unitary cognitive constructs, but rather neurocognitive scores intertwined with sociodemographic and clinical covariates; that is, models reflect stereotypical profiles, and fail when applied to individuals who defy them. Model failure is reliable, phenotype specific and generalizable across datasets. Together, these results highlight the pitfalls of a one-size-fits-all modelling approach and the effect of biased phenotypic measures18-20 on the interpretation and utility of resulting brain-phenotype models. We present a framework to address these issues so that such models may reveal the neural circuits that underlie specific phenotypes and ultimately identify individualized neural targets for clinical intervention.

Mapping Early Brain-Body Interactions: Associations of Fetal Heart Rate Variation with Newborn Brainstem, Hypothalamic, and Dorsal Anterior Cingulate Cortex Functional Connectivity.

The autonomic nervous system (ANS) regulates the body's physiology, including cardiovascular function. As the ANS develops during the second to third trimester, fetal heart rate variability (HRV) increases while fetal heart rate (HR) decreases. In this way, fetal HR and HRV provide an index of fetal ANS development and future neurobehavioral regulation. Fetal HR and HRV have been associated with child language ability and psychomotor development behavior in toddlerhood. However, their associations with postbirth autonomic brain systems, such as the brainstem, hypothalamus, and dorsal anterior cingulate cortex (dACC), have yet to be investigated even though brain pathways involved in autonomic regulation are well established in older individuals. We assessed whether fetal HR and HRV were associated with the brainstem, hypothalamic, and dACC functional connectivity in newborns. Data were obtained from 60 pregnant individuals (ages 14-42) at 24-27 and 34-37 weeks of gestation using a fetal actocardiograph to generate fetal HR and HRV. During natural sleep, their infants (38 males and 22 females) underwent a fMRI scan between 40 and 46 weeks of postmenstrual age. Our findings relate fetal heart indices to brainstem, hypothalamic, and dACC connectivity and reveal connections with widespread brain regions that may support behavioral and emotional regulation. We demonstrated the basic physiologic association between fetal HR indices and lower- and higher-order brain regions involved in regulatory processes. This work provides the foundation for future behavioral or physiological regulation research in fetuses and infants.

Improving power in functional magnetic resonance imaging by moving beyond cluster-level inference.

Inference in neuroimaging typically occurs at the level of focal brain areas or circuits. Yet, increasingly, well-powered studies paint a much richer picture of broad-scale effects distributed throughout the brain, suggesting that many focal reports may only reflect the tip of the iceberg of underlying effects. How focal versus broad-scale perspectives influence the inferences we make has not yet been comprehensively evaluated using real data. Here, we compare sensitivity and specificity across procedures representing multiple levels of inference using an empirical benchmarking procedure that resamples task-based connectomes from the Human Connectome Project dataset (∼1,000 subjects, 7 tasks, 3 resampling group sizes, 7 inferential procedures). Only broad-scale (network and whole brain) procedures obtained the traditional 80% statistical power level to detect an average effect, reflecting >20% more statistical power than focal (edge and cluster) procedures. Power also increased substantially for false discovery rate- compared with familywise error rate-controlling procedures. The downsides are fairly limited; the loss in specificity for broad-scale and FDR procedures was relatively modest compared to the gains in power. Furthermore, the broad-scale methods we introduce are simple, fast, and easy to use, providing a straightforward starting point for researchers. This also points to the promise of more sophisticated broad-scale methods for not only functional connectivity but also related fields, including task-based activation. Altogether, this work demonstrates that shifting the scale of inference and choosing FDR control are both immediately attainable and can help remedy the issues with statistical power plaguing typical studies in the field.

Sex differences in default mode network connectivity in healthy aging adults.

Women show an increased lifetime risk of Alzheimer's disease (AD) compared with men. Characteristic brain connectivity changes, particularly within the default mode network (DMN), have been associated with both symptomatic and preclinical AD, but the impact of sex on DMN function throughout aging is poorly understood. We investigated sex differences in DMN connectivity over the lifespan in 595 cognitively healthy participants from the Human Connectome Project-Aging cohort. We used the intrinsic connectivity distribution (a robust voxel-based metric of functional connectivity) and a seed connectivity approach to determine sex differences within the DMN and between the DMN and whole brain. Compared with men, women demonstrated higher connectivity with age in posterior DMN nodes and lower connectivity in the medial prefrontal cortex. Differences were most prominent in the decades surrounding menopause. Seed-based analysis revealed higher connectivity in women from the posterior cingulate to angular gyrus, which correlated with neuropsychological measures of declarative memory, and hippocampus. Taken together, we show significant sex differences in DMN subnetworks over the lifespan, including patterns in aging women that resemble changes previously seen in preclinical AD. These findings highlight the importance of considering sex in neuroimaging studies of aging and neurodegeneration.

A generalizable connectome-based marker of in-scan sustained attention in neurodiverse youth.

Difficulty with attention is an important symptom in many conditions in psychiatry, including neurodiverse conditions such as autism. There is a need to better understand the neurobiological correlates of attention and leverage these findings in healthcare settings. Nevertheless, it remains unclear if it is possible to build dimensional predictive models of attentional state in a sample that includes participants with neurodiverse conditions. Here, we use 5 datasets to identify and validate functional connectome-based markers of attention. In dataset 1, we use connectome-based predictive modeling and observe successful prediction of performance on an in-scan sustained attention task in a sample of youth, including participants with a neurodiverse condition. The predictions are not driven by confounds, such as head motion. In dataset 2, we find that the attention network model defined in dataset 1 generalizes to predict in-scan attention in a separate sample of neurotypical participants performing the same attention task. In datasets 3-5, we use connectome-based identification and longitudinal scans to probe the stability of the attention network across months to years in individual participants. Our results help elucidate the brain correlates of attentional state in youth and support the further development of predictive dimensional models of other clinically relevant phenotypes.

Simultaneous cortex-wide fluorescence Ca2+ imaging and whole-brain fMRI.

Achieving a comprehensive understanding of brain function requires multiple imaging modalities with complementary strengths. We present an approach for concurrent widefield optical and functional magnetic resonance imaging. By merging these modalities, we can simultaneously acquire whole-brain blood-oxygen-level-dependent (BOLD) and whole-cortex calcium-sensitive fluorescent measures of brain activity. In a transgenic murine model, we show that calcium predicts the BOLD signal, using a model that optimizes a gamma-variant transfer function. We find consistent predictions across the cortex, which are best at low frequency (0.009-0.08 Hz). Furthermore, we show that the relationship between modality connectivity strengths varies by region. Our approach links cell-type-specific optical measurements of activity to the most widely used method for assessing human brain function.

Predicting depressed and elevated mood symptomatology in bipolar disorder using brain functional connectomes.

BACKGROUND: The study is aimed to identify brain functional connectomes predictive of depressed and elevated mood symptomatology in individuals with bipolar disorder (BD) using the machine learning approach Connectome-based Predictive Modeling (CPM). METHODS: Functional magnetic resonance imaging data were obtained from 81 adults with BD while they performed an emotion processing task. CPM with 5000 permutations of leave-one-out cross-validation was applied to identify functional connectomes predictive of depressed and elevated mood symptom scores on the Hamilton Depression and Young Mania rating scales. The predictive ability of the identified connectomes was tested in an independent sample of 43 adults with BD. RESULTS: CPM predicted the severity of depressed [concordance between actual and predicted values (r = 0.23, pperm (permutation test) = 0.031) and elevated (r = 0.27, pperm = 0.01) mood. Functional connectivity of left dorsolateral prefrontal cortex and supplementary motor area nodes, with inter- and intra-hemispheric connections to other anterior and posterior cortical, limbic, motor, and cerebellar regions, predicted depressed mood severity. Connectivity of left fusiform and right visual association area nodes with inter- and intra-hemispheric connections to the motor, insular, limbic, and posterior cortices predicted elevated mood severity. These networks were predictive of mood symptomatology in the independent sample (r ⩾ 0.45, p = 0.002). CONCLUSIONS: This study identified distributed functional connectomes predictive of depressed and elevated mood severity in BD. Connectomes subserving emotional, cognitive, and psychomotor control predicted depressed mood severity, while those subserving emotional and social perceptual functions predicted elevated mood severity. Identification of these connectome networks may help inform the development of targeted treatments for mood symptoms.

Large-scale functional brain networks of maladaptive childhood aggression identified by connectome-based predictive modeling.

Disruptions in frontoparietal networks supporting emotion regulation have been long implicated in maladaptive childhood aggression. However, the association of connectivity between large-scale functional networks with aggressive behavior has not been tested. The present study examined whether the functional organization of the connectome predicts severity of aggression in children. This cross-sectional study included a transdiagnostic sample of 100 children with aggressive behavior (27 females) and 29 healthy controls without aggression or psychiatric disorders (13 females). Severity of aggression was indexed by the total score on the parent-rated Reactive-Proactive Aggression Questionnaire. During fMRI, participants completed a face emotion perception task of fearful and calm faces. Connectome-based predictive modeling with internal cross-validation was conducted to identify brain networks that predicted aggression severity. The replication and generalizability of the aggression predictive model was then tested in an independent sample of children from the Adolescent Brain Cognitive Development (ABCD) study. Connectivity predictive of aggression was identified within and between networks implicated in cognitive control (medial-frontal, frontoparietal), social functioning (default mode, salience), and emotion processing (subcortical, sensorimotor) (r = 0.31, RMSE = 9.05, p = 0.005). Out-of-sample replication (p < 0.002) and generalization (p = 0.007) of findings predicting aggression from the functional connectome was demonstrated in an independent sample of children from the ABCD study (n = 1791; n = 1701). Individual differences in large-scale functional networks contribute to variability in maladaptive aggression in children with psychiatric disorders. Linking these individual differences in the connectome to variation in behavioral phenotypes will advance identification of neural biomarkers of maladaptive childhood aggression to inform targeted treatments.

Predicting the future of neuroimaging predictive models in mental health.

Predictive modeling using neuroimaging data has the potential to improve our understanding of the neurobiology underlying psychiatric disorders and putatively information interventions. Accordingly, there is a plethora of literature reviewing published studies, the mathematics underlying machine learning, and the best practices for using these approaches. As our knowledge of mental health and machine learning continue to evolve, we instead aim to look forward and "predict" topics that we believe will be important in current and future studies. Some of the most discussed topics in machine learning, such as bias and fairness, the handling of dirty data, and interpretable models, may be less familiar to the broader community using neuroimaging-based predictive modeling in psychiatry. In a similar vein, transdiagnostic research and targeting brain-based features for psychiatric intervention are modern topics in psychiatry that predictive models are well-suited to tackle. In this work, we target an audience who is a researcher familiar with the fundamental procedures of machine learning and who wishes to increase their knowledge of ongoing topics in the field. We aim to accelerate the utility and applications of neuroimaging-based predictive models for psychiatric research by highlighting and considering these topics. Furthermore, though not a focus, these ideas generalize to neuroimaging-based predictive modeling in other clinical neurosciences and predictive modeling with different data types (e.g., digital health data).

Atypical Intrinsic Hemispheric Interaction Associated with Autism Spectrum Disorder Is Present within the First Year of Life.

Autism spectrum disorder (ASD) is characterized by atypical connectivity lateralization of functional networks. However, previous studies have not directly investigated if differences in specialization between ASD and typically developing (TD) peers are present in infancy, leaving the timing of onset of these differences relatively unknown. We studied the hemispheric asymmetries of connectivity in children with ASD and infants later meeting the diagnostic criteria for ASD. Analyses were performed in 733 children with ASD and TD peers and in 71 infants at high risk (HR) or normal risk (NR) for ASD, with data collected at 1 month and 9 months of age. Comparing children with ASD (n = 301) to TDs (n = 432), four regions demonstrated group differences in connectivity: posterior cingulate cortex (PCC), posterior superior temporal gyrus, extrastriate cortex, and anterior prefrontal cortex. At 1 month, none of these regions exhibited group differences between ASD (n = 10), HR-nonASD (n = 15), or NR (n = 18) infants. However, by 9 months, the PCC and extrastriate exhibited atypical connectivity in ASD (n = 11) and HR-nonASD infants (n = 24) compared to NR infants (n = 22). Connectivity did not correlate with symptoms in either sample. Our results demonstrate that differences in network asymmetries associated with ASD risk are observable prior to the age of a reliable clinical diagnosis.

Functional connectivity for the language network in the developing brain: 30 weeks of gestation to 30 months of age.

Although the neural scaffolding for language is putatively present before birth, the maturation of functional connections among the key nodes of the language network, Broca's and Wernicke's areas, is less known. We leveraged longitudinal and cross-sectional data from three sites collected through six studies to track the development of functional circuits between Broca's and Wernicke's areas from 30 weeks of gestation through 30 months of age in 127 unique participants. Using resting-state fMRI data, functional connectivity was calculated as the correlation between fMRI time courses from pairs of regions, defined as Broca's and Wernicke's in both hemispheres. The primary analysis evaluated 23 individuals longitudinally imaged from 30 weeks postmenstrual age (fetal) through the first postnatal month (neonatal). A secondary analysis in 127 individuals extended these curves into older infants and toddlers. These data demonstrated significant growth of interhemispheric connections including left Broca's and its homolog and left Wernicke's and its homolog from 30 weeks of gestation through the first postnatal month. In contrast, intrahemispheric connections did not show significant increases across this period. These data represent an important baseline for language systems in the developing brain against which to compare those neurobehavioral disorders with the potential fetal onset of disease.

Functional connectivity predicts changes in attention observed across minutes, days, and months.

The ability to sustain attention differs across people and changes within a single person over time. Although recent work has demonstrated that patterns of functional brain connectivity predict individual differences in sustained attention, whether these same patterns capture fluctuations in attention within individuals remains unclear. Here, across five independent studies, we demonstrate that the sustained attention connectome-based predictive model (CPM), a validated model of sustained attention function, generalizes to predict attentional state from data collected across minutes, days, weeks, and months. Furthermore, the sustained attention CPM is sensitive to within-subject state changes induced by propofol as well as sevoflurane, such that individuals show functional connectivity signatures of stronger attentional states when awake than when under deep sedation and light anesthesia. Together, these results demonstrate that fluctuations in attentional state reflect variability in the same functional connectivity patterns that predict individual differences in sustained attention.

Associations among Household and Neighborhood Socioeconomic Disadvantages, Resting-state Frontoamygdala Connectivity, and Internalizing Symptoms in Youth.

Exposure to socioeconomic disadvantages (SED) can have negative impacts on mental health, yet SED are a multifaceted construct and the precise processes by which SED confer deleterious effects are less clear. Using a large and diverse sample of preadolescents (ages 9-10 years at baseline, n = 4038, 49% female) from the Adolescent Brain Cognitive Development Study, we examined associations among SED at both household (i.e., income-needs and material hardship) and neighborhood (i.e., area deprivation and neighborhood unsafety) levels, frontoamygdala resting-state functional connectivity, and internalizing symptoms at baseline and 1-year follow-up. SED were positively associated with internalizing symptoms at baseline and indirectly predicted symptoms 1 year later through elevated symptoms at baseline. At the household level, youth in households characterized by higher disadvantage (i.e., lower income-to-needs ratio) exhibited more strongly negative frontoamygdala coupling, particularly between the bilateral amygdala and medial OFC (mOFC) regions within the frontoparietal network. Although more strongly positive amygdala-mOFC coupling was associated with higher levels of internalizing symptoms at baseline and 1-year follow-up, it did not mediate the association between income-to-needs ratio and internalizing symptoms. However, at the neighborhood level, amygdala-mOFC functional coupling moderated the effect of neighborhood deprivation on internalizing symptoms. Specifically, higher neighborhood deprivation was associated with higher internalizing symptoms for youth with more strongly positive connectivity, but not for youth with more strongly negative connectivity, suggesting a potential buffering effect. Findings highlight the importance of capturing multilevel socioecological contexts in which youth develop to identify youth who are most likely to benefit from early interventions.

Leveraging edge-centric networks complements existing network-level inference for functional connectomes.

The human connectome is modular with distinct brain regions clustering together to form large-scale communities, or networks. This concept has recently been leveraged in novel inferencing procedures by averaging the edge-level statistics within networks to induce more powerful inferencing at the network level. However, these networks are constructed based on the similarity between pairs of nodes. Emerging work has described novel edge-centric networks, which instead use the similarity between pairs of edges to construct networks. In this work, we use these edge-centric networks in a network-level inferencing procedure and compare this novel method to traditional inferential procedures and the network-level procedure using node-centric networks. We use data from the Human Connectome Project, the Healthy Brain Network, and the Philadelphia Neurodevelopmental Cohort and use a resampling technique with various sample sizes (n=40, 80, 120) to probe the power and specificity of each method. Across datasets and sample sizes, using the edge-centric networks outperforms using node-centric networks for inference as well as edge-level FDR correction and NBS. Additionally, the edge-centric networks were found to be more consistent in clustering effect sizes of similar values as compared to node-centric networks, although node-centric networks often had a lower average within-network effect size variability. Together, these findings suggest that using edge-centric networks for network-level inference can procure relatively powerful results while remaining similarly accurate to the underlying edge-level effects across networks, complementing previous inferential methods.

A cognitive state transformation model for task-general and task-specific subsystems of the brain connectome.

The human brain flexibly controls different cognitive behaviors, such as memory and attention, to satisfy contextual demands. Much progress has been made to reveal task-induced modulations in the whole-brain functional connectome, but we still lack a way to model context-dependent changes. Here, we present a novel connectome-to-connectome (C2C) transformation framework that enables us to model the brain's functional reorganization from one connectome state to another in response to specific task goals. Using functional magnetic resonance imaging data from the Human Connectome Project, we demonstrate that the C2C model accurately generates an individual's task-related connectomes from their task-free (resting-state) connectome with a high degree of specificity across seven different cognitive states. Moreover, the C2C model amplifies behaviorally relevant individual differences in the task-free connectome, thereby improving behavioral predictions with increased power, achieving similar performance with just a third of the subjects needed when relying on resting-state data alone. Finally, the C2C model reveals how the brain reorganizes between cognitive states. Our observations support the existence of reliable state-specific subsystems in the brain and demonstrate that we can quantitatively model how the connectome reconfigures to different cognitive states, enabling more accurate predictions of behavior with fewer subjects.

Large-scale differences in functional organization of left- and right-handed individuals using whole-brain, data-driven analysis of connectivity.

Handedness influences differences in lateralization of language areas as well as dominance of motor and somatosensory cortices. However, differences in whole-brain functional connectivity (i.e., functional connectomes) due to handedness have been relatively understudied beyond pre-specified networks of interest. Here, we compared functional connectomes of left- and right-handed individuals at the whole brain level. We explored differences in functional connectivity of previously established regions of interest, and showed differences between primarily left- and primarily right-handed individuals in the motor, somatosensory, and language areas using functional connectivity. We then proceeded to investigate these differences in the whole brain and found that the functional connectivity of left- and right-handed individuals are not specific to networks of interest, but extend across every region of the brain. In particular, we found that connections between and within the cerebellum show distinct patterns of connectivity. To put these effects into context, we show that the effect sizes associated with handedness differences account for a similar amount of individual differences in the connectome as sex differences. Together these results shed light on regions of the brain beyond those traditionally explored that contribute to differences in the functional organization of left- and right-handed individuals and underscore that handedness effects are neurobiologically meaningful in addition to being statistically significant.

Arousal impacts distributed hubs modulating the integration of brain functional connectivity.

Even when subjects are at rest, it is thought that brain activity is organized into distinct brain states during which reproducible patterns are observable. Yet, it is unclear how to define or distinguish different brain states. A potential source of brain state variation is arousal, which may play a role in modulating functional interactions between brain regions. Here, we use simultaneous resting state functional magnetic resonance imaging (fMRI) and pupillometry to study the impact of arousal levels indexed by pupil area on the integration of large-scale brain networks. We employ a novel sparse dictionary learning-based method to identify hub regions participating in between-network integration stratified by arousal, by measuring k-hubness, the number (k) of functionally overlapping networks in each brain region. We show evidence of a brain-wide decrease in between-network integration and inter-subject variability at low relative to high arousal, with differences emerging across regions of the frontoparietal, default mode, motor, limbic, and cerebellum networks. State-dependent changes in k-hubness relate to the actual patterns of network integration within these hubs, suggesting a brain state transition from high to low arousal characterized by global synchronization and reduced network overlaps. We demonstrate that arousal is not limited to specific brain areas known to be directly associated with arousal regulation, but instead has a brain-wide impact that involves high-level between-network communications. Lastly, we show a systematic change in pairwise fMRI signal correlation structures in the arousal state-stratified data, and demonstrate that the choice of global signal regression could result in different conclusions in conventional graph theoretical analysis and in the analysis of k-hubness when studying arousal modulations. Together, our results suggest the presence of global and local effects of pupil-linked arousal modulations on resting state brain functional connectivity.

(Un)common space in infant neuroimaging studies: A systematic review of infant templates.

In neuroimaging, spatial normalization is an important step that maps an individual's brain onto a template brain permitting downstream statistical analyses. Yet, in infant neuroimaging, there remain several technical challenges that have prevented the establishment of a standardized template for spatial normalization. Thus, many different approaches are used in the literature. To quantify the popularity and variability of these approaches in infant neuroimaging studies, we performed a systematic review of infant magnetic resonance imaging (MRI) studies from 2000 to 2020. Here, we present results from 834 studies meeting inclusion criteria. Studies were classified into (a) processing data in single subject space, (b) using an off the shelf, or "off the shelf," template, (c) creating a study specific template, or (d) using a hybrid of these methods. We found that across the studies in the systematic review, single subject space was the most used (no common space). This was the most used common space for diffusion-weighted imaging and structural MRI studies while functional MRI studies preferred off the shelf atlases. We found a pattern such that more recently published studies are more commonly using off the shelf atlases. When considering special populations, preterm studies most used single subject space while, when no special populations were being analyzed, an off the shelf template was most common. The most used off the shelf templates were the UNC Infant Atlases (24%). Using a systematic review of infant neuroimaging studies, we highlight a lack of an established "standard" template brain in these studies.

Age-related topographic map of magnetic resonance diffusion metrics in neonatal brains.

Accelerated maturation of brain parenchyma close to term-equivalent age leads to rapid changes in diffusion-weighted imaging (DWI) and diffusion tensor imaging (DTI) metrics of neonatal brains, which can complicate the evaluation and interpretation of these scans. In this study, we characterized the topography of age-related evolution of diffusion metrics in neonatal brains. We included 565 neonates who had MRI between 0 and 3 months of age, with no structural or signal abnormality-including 162 who had DTI scans. We analyzed the age-related changes of apparent diffusion coefficient (ADC) values throughout brain and DTI metrics (fractional anisotropy [FA] and mean diffusivity [MD]) along white matter (WM) tracts. Rate of change in ADC, FA, and MD values across 5 mm cubic voxels was calculated. There was significant reduction of ADC and MD values and increase of FA with increasing gestational age (GA) throughout neonates' brain, with the highest temporal rates in subcortical WM, corticospinal tract, cerebellar WM, and vermis. GA at birth had significant effect on ADC values in convexity cortex and corpus callosum as well as FA/MD values in corpus callosum, after correcting for GA at scan. We developed online interactive atlases depicting age-specific normative values of ADC (ages 34-46 weeks), and FA/MD (35-41 weeks). Our results show a rapid decrease in diffusivity metrics of cerebral/cerebellar WM and vermis in the first few weeks of neonatal age, likely attributable to myelination. In addition, prematurity and low GA at birth may result in lasting delay in corpus callosum myelination and cerebral cortex cellularity.

Associations between grip strength, brain structure, and mental health in > 40,000 participants from the UK Biobank.

BACKGROUND: Grip strength is a widely used and well-validated measure of overall health that is increasingly understood to index risk for psychiatric illness and neurodegeneration in older adults. However, existing work has not examined how grip strength relates to a comprehensive set of mental health outcomes, which can detect early signs of cognitive decline. Furthermore, whether brain structure mediates associations between grip strength and cognition remains unknown. METHODS: Based on cross-sectional and longitudinal data from over 40,000 participants in the UK Biobank, this study investigated the behavioral and neural correlates of handgrip strength using a linear mixed effect model and mediation analysis. RESULTS: In cross-sectional analysis, we found that greater grip strength was associated with better cognitive functioning, higher life satisfaction, greater subjective well-being, and reduced depression and anxiety symptoms while controlling for numerous demographic, anthropometric, and socioeconomic confounders. Further, grip strength of females showed stronger associations with most behavioral outcomes than males. In longitudinal analysis, baseline grip strength was related to cognitive performance at ~9 years follow-up, while the reverse effect was much weaker. Further, baseline neuroticism, health, and financial satisfaction were longitudinally associated with subsequent grip strength. The results revealed widespread associations between stronger grip strength and increased grey matter volume, especially in subcortical regions and temporal cortices. Moreover, grey matter volume of these regions also correlated with better mental health and considerably mediated their relationship with grip strength. CONCLUSIONS: Overall, using the largest population-scale neuroimaging dataset currently available, our findings provide the most well-powered characterization of interplay between grip strength, mental health, and brain structure, which may facilitate the discovery of possible interventions to mitigate cognitive decline during aging.