Amplitudes of resting-state functional networks - investigation into their correlates and biophysical properties.

Resting-state fMRI studies have shown that multiple functional networks, which consist of distributed brain regions that share synchronised spontaneous activity, co-exist in the brain. As these resting-state networks (RSNs) have been thought to reflect the brain's intrinsic functional organization, intersubject variability in the networks' spontaneous fluctuations may be associated with individuals' clinical, physiological, cognitive, and genetic traits. Here, we investigated resting-state fMRI data along with extensive clinical, lifestyle, and genetic data collected from 37,842 UK Biobank participants, with the object of elucidating intersubject variability in the fluctuation amplitudes of RSNs. Functional properties of the RSN amplitudes were first examined by analyzing correlations with the well-established between-network functional connectivity. It was found that a network amplitude is highly correlated with the mean strength of the functional connectivity that the network has with the other networks. Intersubject clustering analysis showed the amplitudes are most strongly correlated with age, cardiovascular factors, body composition, blood cell counts, lung function, and sex, with some differences in the correlation strengths between sensory and cognitive RSNs. Genome-wide association studies (GWASs) of RSN amplitudes identified several significant genetic variants reported in previous GWASs for their implications in sleep duration. We provide insight into key factors determining RSN amplitudes and demonstrate that intersubject variability of the amplitudes primarily originates from differences in temporal synchrony between functionally linked brain regions, rather than differences in the magnitude of raw voxelwise BOLD signal changes. This finding additionally revealed intriguing differences between sensory and cognitive RSNs with respect to sex effects on temporal synchrony and provided evidence suggesting that synchronous coactivations of functionally linked brain regions, and magnitudes of BOLD signal changes, may be related to different genetic mechanisms. These results underscore that intersubject variability of the amplitudes in health and disease need to be interpreted largely as a measure of the sum of within-network temporal synchrony and amplitudes of BOLD signals, with a dominant contribution from the former.

Comparison between gradients and parcellations for functional connectivity prediction of behavior.

Resting-state functional connectivity (RSFC) is widely used to predict behavioral measures. To predict behavioral measures, representing RSFC with parcellations and gradients are the two most popular approaches. Here, we compare parcellation and gradient approaches for RSFC-based prediction of a broad range of behavioral measures in the Human Connectome Project (HCP) and Adolescent Brain Cognitive Development (ABCD) datasets. Among the parcellation approaches, we consider group-average "hard" parcellations (Schaefer et al., 2018), individual-specific "hard" parcellations (Kong et al., 2021a), and an individual-specific "soft" parcellation (spatial independent component analysis with dual regression; Beckmann et al., 2009). For gradient approaches, we consider the well-known principal gradients (Margulies et al., 2016) and the local gradient approach that detects local RSFC changes (Laumann et al., 2015). Across two regression algorithms, individual-specific hard-parcellation performs the best in the HCP dataset, while the principal gradients, spatial independent component analysis and group-average "hard" parcellations exhibit similar performance. On the other hand, principal gradients and all parcellation approaches perform similarly in the ABCD dataset. Across both datasets, local gradients perform the worst. Finally, we find that the principal gradient approach requires at least 40 to 60 gradients to perform as well as parcellation approaches. While most principal gradient studies utilize a single gradient, our results suggest that incorporating higher order gradients can provide significant behaviorally relevant information. Future work will consider the inclusion of additional parcellation and gradient approaches for comparison.

A Phenome-Wide Association Study (PheWAS) of Late Onset Alzheimer Disease Genetic Risk in Children of European Ancestry at Middle Childhood: Results from the ABCD Study.

Genetic risk for Late Onset Alzheimer Disease (AD) has been associated with lower cognition and smaller hippocampal volume in healthy young adults. However, whether these and other associations are present during childhood remains unclear. Using data from 5556 genomically-confirmed European ancestry youth who completed the baseline session of the ongoing the Adolescent Brain Cognitive DevelopmentSM Study (ABCD Study®), our phenome-wide association study estimating associations between four indices of genetic risk for late-onset AD (i.e., AD polygenic risk scores (PRS), APOE rs429358 genotype, AD PRS with the APOE region removed (ADPRS-APOE), and an interaction between ADPRS-APOE and APOE genotype) and 1687 psychosocial, behavioral, and neural phenotypes revealed no significant associations after correction for multiple testing (all ps > 0.0002; all pfdr > 0.07). These data suggest that AD genetic risk may not phenotypically manifest during middle-childhood or that effects are smaller than this sample is powered to detect.

Mental health in the UK Biobank: A roadmap to self-report measures and neuroimaging correlates.

The UK Biobank (UKB) is a highly promising dataset for brain biomarker research into population mental health due to its unprecedented sample size and extensive phenotypic, imaging, and biological measurements. In this study, we aimed to provide a shared foundation for UKB neuroimaging research into mental health with a focus on anxiety and depression. We compared UKB self-report measures and revealed important timing effects between scan acquisition and separate online acquisition of some mental health measures. To overcome these timing effects, we introduced and validated the Recent Depressive Symptoms (RDS-4) score which we recommend for state-dependent and longitudinal research in the UKB. We furthermore tested univariate and multivariate associations between brain imaging-derived phenotypes (IDPs) and mental health. Our results showed a significant multivariate relationship between IDPs and mental health, which was replicable. Conversely, effect sizes for individual IDPs were small. Test-retest reliability of IDPs was stronger for measures of brain structure than for measures of brain function. Taken together, these results provide benchmarks and guidelines for future UKB research into brain biomarkers of mental health.

Recent developments in representations of the connectome.

Research into the human connectome (i.e., all connections in the human brain) with the use of resting state functional MRI has rapidly increased in popularity in recent years, especially with the growing availability of large-scale neuroimaging datasets. The goal of this review article is to describe innovations in functional connectome representations that have come about in the past 8 years, since the 2013 NeuroImage special issue on 'Mapping the Connectome'. In the period, research has shifted from group-level brain parcellations towards the characterization of the individualized connectome and of relationships between individual connectomic differences and behavioral/clinical variation. Achieving subject-specific accuracy in parcel boundaries while retaining cross-subject correspondence is challenging, and a variety of different approaches are being developed to meet this challenge, including improved alignment, improved noise reduction, and robust group-to-subject mapping approaches. Beyond the interest in the individualized connectome, new representations of the data are being studied to complement the traditional parcellated connectome representation (i.e., pairwise connections between distinct brain regions), such as methods that capture overlapping and smoothly varying patterns of connectivity ('gradients'). These different connectome representations offer complimentary insights into the inherent functional organization of the brain, but challenges for functional connectome research remain. Interpretability will be improved by future research towards gaining insights into the neural mechanisms underlying connectome observations obtained from functional MRI. Validation studies comparing different connectome representations are also needed to build consensus and confidence to proceed with clinical trials that may produce meaningful clinical translation of connectome insights.

Hierarchical modelling of functional brain networks in population and individuals from big fMRI data.

A major goal of large-scale brain imaging datasets is to provide resources for investigating heterogeneous populations. Characterisation of functional brain networks for individual subjects from these datasets will have an enormous potential for prediction of cognitive or clinical traits. We propose for the first time a technique, Stochastic Probabilistic Functional Modes (sPROFUMO), that is scalable to UK Biobank (UKB) with expected 100,000 participants, and hierarchically estimates functional brain networks in individuals and the population, while allowing for bidirectional flow of information between the two. Using simulations, we show the model's utility, especially in scenarios that involve significant cross-subject variability, or require delineation of fine-grained differences between the networks. Subsequently, by applying the model to resting-state fMRI from 4999 UKB subjects, we mapped resting state networks (RSNs) in single subjects with greater detail than has been possible previously in UKB (>100 RSNs), and demonstrate that these RSNs can predict a range of sensorimotor and higher-level cognitive functions. Furthermore, we demonstrate several advantages of the model over independent component analysis combined with dual-regression (ICA-DR), particularly with respect to the estimation of the spatial configuration of the RSNs and the predictive power for cognitive traits. The proposed model and results can open a new door for future investigations into individualised profiles of brain function from big data.

Modelling subject variability in the spatial and temporal characteristics of functional modes.

Recent work has highlighted the scale and ubiquity of subject variability in observations from functional MRI data (fMRI). Furthermore, it is highly likely that errors in the estimation of either the spatial presentation of, or the coupling between, functional regions can confound cross-subject analyses, making accurate and unbiased representations of functional data essential for interpreting any downstream analyses. Here, we extend the framework of probabilistic functional modes (PFMs) (Harrison et al., 2015) to capture cross-subject variability not only in the mode spatial maps, but also in the functional coupling between modes and in mode amplitudes. A new implementation of the inference now also allows for the analysis of modern, large-scale data sets, and the combined inference and analysis package, PROFUMO, is available from git.fmrib.ox.ac.uk/samh/profumo. A new implementation of the inference now also allows for the analysis of modern, large-scale data sets. Using simulated data, resting-state data from 1000 subjects collected as part of the Human Connectome Project (Van Essen et al., 2013), and an analysis of 14 subjects in a variety of continuous task-states (Kieliba et al., 2019), we demonstrate how PFMs are able to capture, within a single model, a rich description of how the spatio-temporal structure of resting-state fMRI activity varies across subjects. We also compare the new PFM model to the well established independent component analysis with dual regression (ICA-DR) pipeline. This reveals that, under PFM assumptions, much more of the (behaviorally relevant) cross-subject variability in fMRI activity should be attributed to the variability in spatial maps, and that, after accounting for this, functional coupling between modes primarily reflects current cognitive state. This has fundamental implications for the interpretation of cross-sectional studies of functional connectivity that do not capture cross-subject variability to the same extent as PFMs.

Classification of temporal ICA components for separating global noise from fMRI data: Reply to Power.

We respond to a critique of our temporal Independent Components Analysis (ICA) method for separating global noise from global signal in fMRI data that focuses on the signal versus noise classification of several components. While we agree with several of Power's comments, we provide evidence and analysis to rebut his major criticisms and to reassure readers that temporal ICA remains a powerful and promising denoising approach.

Using temporal ICA to selectively remove global noise while preserving global signal in functional MRI data.

Temporal fluctuations in functional Magnetic Resonance Imaging (fMRI) have been profitably used to study brain activity and connectivity for over two decades. Unfortunately, fMRI data also contain structured temporal "noise" from a variety of sources, including subject motion, subject physiology, and the MRI equipment. Recently, methods have been developed to automatically and selectively remove spatially specific structured noise from fMRI data using spatial Independent Components Analysis (ICA) and machine learning classifiers. Spatial ICA is particularly effective at removing spatially specific structured noise from high temporal and spatial resolution fMRI data of the type acquired by the Human Connectome Project and similar studies. However, spatial ICA is mathematically, by design, unable to separate spatially widespread "global" structured noise from fMRI data (e.g., blood flow modulations from subject respiration). No methods currently exist to selectively and completely remove global structured noise while retaining the global signal from neural activity. This has left the field in a quandary-to do or not to do global signal regression-given that both choices have substantial downsides. Here we show that temporal ICA can selectively segregate and remove global structured noise while retaining global neural signal in both task-based and resting state fMRI data. We compare the results before and after temporal ICA cleanup to those from global signal regression and show that temporal ICA cleanup removes the global positive biases caused by global physiological noise without inducing the network-specific negative biases of global signal regression. We believe that temporal ICA cleanup provides a "best of both worlds" solution to the global signal and global noise dilemma and that temporal ICA itself unlocks interesting neurobiological insights from fMRI data.

Investigations into within- and between-subject resting-state amplitude variations.

The amplitudes of spontaneous fluctuations in brain activity may be a significant source of within-subject and between-subject variability, and this variability is likely to be carried through into functional connectivity (FC) estimates (whether directly or indirectly). Therefore, improving our understanding of amplitude fluctuations over the course of a resting state scan and variation in amplitude across individuals is of great relevance to the interpretation of FC findings. We investigate resting state amplitudes in two large-scale studies (HCP and UK Biobank), with the aim of determining between-subject and within-subject variability. Between-subject clustering distinguished between two groups of brain networks whose amplitude variation across subjects were highly correlated with each other, revealing a clear distinction between primary sensory and motor regions ('primary sensory/motor cluster') and cognitive networks. Within subjects, all networks in the primary sensory/motor cluster showed a consistent increase in amplitudes from the start to the end of the scan. In addition to the strong increases in primary sensory/motor amplitude, a large number of changes in FC were found when comparing the two scans acquired on the same day (HCP data). Additive signal change analysis confirmed that all of the observed FC changes could be fully explained by changes in amplitude. Between-subject correlations in UK Biobank data showed a negative correlation between primary sensory/motor amplitude and average sleep duration, suggesting a role of arousal. Our findings additionally reveal complex relationships between amplitude and head motion. These results suggest that network amplitude is a source of significant variability both across subjects, and within subjects on a within-session timescale. Future rfMRI studies may benefit from obtaining arousal-related (self report) measures, and may wish to consider the influence of amplitude changes on measures of (dynamic) functional connectivity.

Hand classification of fMRI ICA noise components.

We present a practical "how-to" guide to help determine whether single-subject fMRI independent components (ICs) characterise structured noise or not. Manual identification of signal and noise after ICA decomposition is required for efficient data denoising: to train supervised algorithms, to check the results of unsupervised ones or to manually clean the data. In this paper we describe the main spatial and temporal features of ICs and provide general guidelines on how to evaluate these. Examples of signal and noise components are provided from a wide range of datasets (3T data, including examples from the UK Biobank and the Human Connectome Project, and 7T data), together with practical guidelines for their identification. Finally, we discuss how the data quality, data type and preprocessing can influence the characteristics of the ICs and present examples of particularly challenging datasets.

Functional Connectivity under Anticipation of Shock: Correlates of Trait Anxious Affect versus Induced Anxiety.

Sustained anxiety about potential future negative events is an important feature of anxiety disorders. In this study, we used a novel anticipation of shock paradigm to investigate individual differences in functional connectivity during prolonged threat of shock. We examined the correlates of between-participant differences in trait anxious affect and induced anxiety, where the latter reflects changes in self-reported anxiety resulting from the shock manipulation. Dissociable effects of trait anxious affect and induced anxiety were observed. Participants with high scores on a latent dimension of anxious affect showed less increase in ventromedial pFC-amygdala connectivity between periods of safety and shock anticipation. Meanwhile, lower levels of induced anxiety were linked to greater augmentation of dorsolateral pFC-anterior insula connectivity during shock anticipation. These findings suggest that ventromedial pFC-amygdala and dorsolateral pFC-insula networks might both contribute to regulation of sustained fear responses, with their recruitment varying independently across participants. The former might reflect an evolutionarily old mechanism for reducing fear or anxiety, whereas the latter might reflect a complementary mechanism by which cognitive control can be implemented to diminish fear responses generated due to anticipation of aversive stimuli or events. These two circuits might provide complementary, alternate targets for exploration in future pharmacological and cognitive intervention studies.

Resting state correlates of subdimensions of anxious affect.

Resting state fMRI may help identify markers of risk for affective disorder. Given the comorbidity of anxiety and depressive disorders and the heterogeneity of these disorders as defined by DSM, an important challenge is to identify alterations in resting state brain connectivity uniquely associated with distinct profiles of negative affect. The current study aimed to address this by identifying differences in brain connectivity specifically linked to cognitive and physiological profiles of anxiety, controlling for depressed affect. We adopted a two-stage multivariate approach. Hierarchical clustering was used to independently identify dimensions of negative affective style and resting state brain networks. Combining the clustering results, we examined individual differences in resting state connectivity uniquely associated with subdimensions of anxious affect, controlling for depressed affect. Physiological and cognitive subdimensions of anxious affect were identified. Physiological anxiety was associated with widespread alterations in insula connectivity, including decreased connectivity between insula subregions and between the insula and other medial frontal and subcortical networks. This is consistent with the insula facilitating communication between medial frontal and subcortical regions to enable control of physiological affective states. Meanwhile, increased connectivity within a frontoparietal-posterior cingulate cortex-precunous network was specifically associated with cognitive anxiety, potentially reflecting increased spontaneous negative cognition (e.g., worry). These findings suggest that physiological and cognitive anxiety comprise subdimensions of anxiety-related affect and reveal associated alterations in brain connectivity.

Study protocol: The Whitehall II imaging sub-study.

BACKGROUND: The Whitehall II (WHII) study of British civil servants provides a unique source of longitudinal data to investigate key factors hypothesized to affect brain health and cognitive ageing. This paper introduces the multi-modal magnetic resonance imaging (MRI) protocol and cognitive assessment designed to investigate brain health in a random sample of 800 members of the WHII study. METHODS/DESIGN: A total of 6035 civil servants participated in the WHII Phase 11 clinical examination in 2012-2013. A random sample of these participants was included in a sub-study comprising an MRI brain scan, a detailed clinical and cognitive assessment, and collection of blood and buccal mucosal samples for the characterisation of immune function and associated measures. Data collection for this sub-study started in 2012 and will be completed by 2016. The participants, for whom social and health records have been collected since 1985, were between 60-85 years of age at the time the MRI study started. Here, we describe the pre-specified clinical and cognitive assessment protocols, the state-of-the-art MRI sequences and latest pipelines for analyses of this sub-study. DISCUSSION: The integration of cutting-edge MRI techniques, clinical and cognitive tests in combination with retrospective data on social, behavioural and biological variables during the preceding 25 years from a well-established longitudinal epidemiological study (WHII cohort) will provide a unique opportunity to examine brain structure and function in relation to age-related diseases and the modifiable and non-modifiable factors affecting resilience against and vulnerability to adverse brain changes.

Opaque Ontology: Neuroimaging Classification of ICD-10 Diagnostic Groups in the UK Biobank.

Background: 1.The use of machine learning to classify diagnostic cases versus controls defined based on diagnostic ontologies such as the ICD-10 from neuroimaging features is now commonplace across a wide range of diagnostic fields. However, transdiagnostic comparisons of such classifications are lacking. Such transdiagnostic comparisons are important to establish the specificity of classification models, set benchmarks, and assess the value of diagnostic ontologies. Results: 2.We investigated case-control classification accuracy in 17 different ICD-10 diagnostic groups from Chapter V (mental and behavioral disorders) and Chapter VI (diseases of the nervous system) using data from the UK Biobank. Classification models were trained using either neuroimaging (structural or functional brain MRI feature sets) or socio-demographic features. Random forest classification models were adopted using rigorous shuffle splits to estimate stability as well as accuracy of case-control classifications. Diagnostic classification accuracies were benchmarked against age classification (oldest versus youngest) from the same feature sets and against additional classifier types (K-nearest neighbors and linear support vector machine). In contrast to age classification accuracy, which was high for all feature sets, few ICD-10 diagnostic groups were classified significantly above chance (namely, demyelinating diseases based on structural neuroimaging features, and depression based on socio-demographic and functional neuroimaging features). Conclusion: 3.These findings highlight challenges with the current disease classification system, leading us to recommend caution with the use of ICD-10 diagnostic groups as target labels in brain-based disease prediction studies.

Multiscale Modes of Functional Brain Connectivity.

Information processing in the brain spans from localised sensorimotor processes to higher-level cognition that integrates across multiple regions. Interactions between and within these subsystems enable multiscale information processing. Despite this multiscale characteristic, functional brain connectivity is often either estimated based on 10-30 distributed modes or parcellations with 100-1000 localised parcels, both missing across-scale functional interactions. We present Multiscale Probabilistic Functional Modes (mPFMs), a new mapping which comprises modes over various scales of granularity, thus enabling direct estimation of functional connectivity within- and across-scales. Crucially, mPFMs emerged from data-driven multilevel Bayesian modelling of large functional MRI (fMRI) populations. We demonstrate that mPFMs capture both distributed brain modes and their co-existing subcomponents. In addition to validating mPFMs using simulations and real data, we show that mPFMs can predict ~900 personalised traits from UK Biobank more accurately than current standard techniques. Therefore, mPFMs can offer a paradigm shift in functional connectivity modelling and yield enhanced fMRI biomarkers for traits and diseases.

Hippocampal volumes in UK Biobank are associated with APOE only in older adults.

INTRODUCTION: The hippocampus atrophies with age and is implicated in neurodegenerative disorders including Alzheimer's disease (AD). We examined the interplay between age and APOE genotype on total hippocampal volume. METHODS: Utilizing neuroimaging data from 37,463 UK Biobank participants, we applied linear regression to quantify the association of age and APOE with hippocampal volume and identified the age when volumes of ε2/ε3, ε3/ε4, and ε4/ε4 carriers significantly deviated from ε3/ε3 using generalized additive modeling. RESULTS: Total hippocampal volume declined with age, with significant differences by APOE genotype emerging after age 60. ε3/ε4 and ε4/ε4 carriers displayed reduced volumes from ages 69 and 61, respectively, while ε2/ε3 showed delayed decline starting at age 76. DISCUSSION: The association of APOE and hippocampal volume is age-dependent, with differences in volumes of ε4/ε4 carriers detected as early as age 61. This work underscores the importance of APOE genotype in determining when to begin screening for AD.

Associations between prenatal caffeine exposure and child development: Longitudinal results from the Adolescent Brain Cognitive Development (ABCD) Study.

Objective: Though caffeine use during pregnancy is common, its longitudinal associations with child behavioral and physical health outcomes remain poorly understood. Here, we estimated associations between prenatal caffeine exposure, body mass index (BMI), and behavior as children enter adolescence. Method: Longitudinal data and caregiver-reported prenatal caffeine exposure were obtained from the ongoing Adolescent Brain and Cognitive Development (ABCD) SM Study, which recruited 11,875 children aged 9-11 years at baseline from 21 sites across the United States starting June 1, 2016. Prenatal caffeine exposure was analyzed as a 4-level categorical variable, and further group contrasts were used to characterize "any exposure" and "daily exposure" groups. Outcomes included psychopathology characteristics in children, sleep problems, and BMI. Potentially confounding covariates included familial (e.g., income, familial psychopathology), pregnancy (e.g., prenatal substance exposure), and child (e.g., caffeine use) variables. Results: Among 10,873 children (5,686 boys [52.3%]; mean [SD] age, 9.9 [0.6] years) with nonmissing prenatal caffeine exposure data, 6,560 (60%) were exposed to caffeine prenatally. Relative to no exposure, daily caffeine exposure was associated with higher child BMI (ß=0.08; FDR-corrected p=0.02), but was not associated with child behavior. Those exposed to two or more cups of caffeine daily (n=1,028) had greater sleep problems than those with lower/no exposure (ß>0.92; FDR-corrected p<0.04). Conclusion: Daily prenatal caffeine exposure is associated with heightened childhood BMI, and when used multiple times a day greater sleep problems even after accounting for potential confounds. Whether this relationship is a consequence of prenatal caffeine exposure or its correlated factors remains unknown.

Associations between COVID-19 and putative markers of neuroinflammation: A diffusion basis spectrum imaging study.

COVID-19 remains a significant international public health concern. Yet, the mechanisms through which symptomatology emerges remain poorly understood. While SARS-CoV-2 infection may induce prolonged inflammation within the central nervous system, the evidence primarily stems from limited small-scale case investigations. To address this gap, our study capitalized on longitudinal UK Biobank neuroimaging data acquired prior to and following COVID-19 testing (N = 416 including n = 224 COVID-19 cases; Mage = 58.6). Putative neuroinflammation was assessed in gray matter structures and white matter tracts using non-invasive Diffusion Basis Spectrum Imaging (DBSI), which estimates inflammation-related cellularity (DBSI-restricted fraction; DBSI-RF) and vasogenic edema (DBSI-hindered fraction; DBSI-HF). We hypothesized that COVID-19 case status would be associated with increases in DBSI markers after accounting for potential confound (age, sex, race, body mass index, smoking frequency, and data acquisition interval) and multiple testing. COVID-19 case status was not significantly associated with DBSI-RF (|ß|'s < 0.28, pFDR >0.05), but with greater DBSI-HF in left pre- and post-central gyri and right middle frontal gyrus (ß's > 0.3, all pFDR = 0.03). Intriguingly, the brain areas exhibiting increased putative vasogenic edema had previously been linked to COVID-19-related functional and structural alterations, whereas brain regions displaying subtle differences in cellularity between COVID-19 cases and controls included regions within or functionally connected to the olfactory network, which has been implicated in COVID-19 psychopathology. Nevertheless, our study might not have captured acute and transitory neuroinflammatory effects linked to SARS-CoV-2 infection, possibly due to symptom resolution before the imaging scan. Future research is warranted to explore the potential time- and symptom-dependent neuroinflammatory relationship with COVID-19.

Investigating the Relationship Between Smoking Behavior and Global Brain Volume.

Background: Previous studies have shown that brain volume is negatively associated with cigarette smoking, but there is an ongoing debate about whether smoking causes lowered brain volume or a lower brain volume is a risk factor for smoking. We address this debate through multiple methods that evaluate directionality: Bradford Hill's criteria, which are commonly used to understand a causal relationship in epidemiological studies, and mediation analysis. Methods: In 32,094 participants of European descent from the UK Biobank dataset, we examined the relationship between a history of daily smoking and brain volumes, as well as an association of genetic risk score to ever smoking with brain volume. Results: A history of daily smoking was strongly associated with decreased brain volume, and a history of heavier smoking was associated with a greater decrease in brain volume. The strongest association was between total gray matter volume and a history of daily smoking (effect size = -2964 mm3, p = 2.04 × 10-16), and there was a dose-response relationship with more pack years smoked associated with a greater decrease in brain volume. A polygenic risk score for smoking initiation was strongly associated with a history of daily smoking (effect size = 0.05, p = 4.20 × 10-84), but only modestly associated with total gray matter volume (effect size = -424 mm3, p = .01). Mediation analysis indicated that a history of daily smoking mediated the relationship between the smoking initiation polygenic risk score and total gray matter volume. Conclusions: A history of daily smoking is strongly associated with a decreased total brain volume.