Dietary and Lifestyle Factors of Brain Iron Accumulation and Parkinson's Disease Risk

Purpose

Iron is an essential nutrient which can only be absorbed through an individual's diet. Excess iron accumulates in organs throughout the body including the brain. Iron dysregulation in the brain is commonly associated with neurodegenerative diseases like Alzheimer's disease and Parkinson's Disease (PD). Our previous research has shown that a pattern of iron accumulation in motor regions of the brain related to a genetic iron-storage disorder called hemochromatosis is associated with an increased risk of PD. To understand how diet and lifestyle factors relate to this brain endophenotype and risk of PD we analyzed the relationship between these measures, estimates of nutrient intake, and diet and lifestyle preference using data from UK Biobank.

Methods

Using distinct imaging and non-imaging samples (20,477 to 28,388 and 132,023 to 150,603 participants, respectively), we performed linear and logistic regression analyses using estimated dietary nutrient intake and food preferences to predict a) brain iron accumulation score (derived from T2-Weighted Magnetic Resonance Imaging) and b) PD risk. In addition, we performed a factor analysis of diet and lifestyle preferences to investigate if latent lifestyle factors explained significant associations. Finally, we performed an instrumental variable regression of our results related to iron accumulation and PD risk to identify if there were common dietary and lifestyle factors that were jointly associated with differences in brain iron accumulation and PD risk.

Results

We found multiple highly significant associations with measures of brain iron accumulation and preferences for alcohol (factor 7: t=4.02, pFDR=0.0003), exercise (factor 11: t=-4.31, pFDR=0.0001), and high-sugar foods (factor 2: t=-3.73, pFDR=0.0007). Preference for alcohol (factor 7: t=-5.83, pFDR<1×10-8), exercise (factor 11: t=-7.66, pFDR<1×10-13), and high sugar foods (factor 2: t=6.03, pFDR<1×10-8) were also associated with PD risk. Instrumental variable regression of individual preferences revealed a significant relationship in which dietary preferences associated with higher brain iron levels also appeared to be linked to a lower risk for PD (p=0.004). A similar relationship was observed for estimates of nutrient intake (p=0.0006). Voxel-wise analysis of i) high-sugar and ii) alcohol factors confirmed T2-weighted signal differences consistent with iron accumulation patterns in motor regions of the brain including the cerebellum and basal ganglia.

Conclusion

Dietary and lifestyle factors and preferences, especially those related to carbohydrates, alcohol, and exercise, are related to detectable differences in brain iron accumulation and alterations in risk of PD, suggesting a potential avenue for lifestyle interventions that could influence risk.

Polygenic profiles define aspects of clinical heterogeneity in attention deficit hyperactivity disorder

Attention deficit hyperactivity disorder (ADHD) is a complex disorder that manifests variability in long-term outcomes and clinical presentations. The genetic contributions to such heterogeneity are not well understood. Here we show several genetic links to clinical heterogeneity in ADHD in a case-only study of 14,084 diagnosed individuals. First, we identify one genome-wide significant locus by comparing cases with ADHD and autism spectrum disorder (ASD) to cases with ADHD but not ASD. Second, we show that cases with ASD and ADHD, substance use disorder and ADHD, or first diagnosed with ADHD in adulthood have unique polygenic score (PGS) profiles that distinguish them from complementary case subgroups and controls. Finally, a PGS for an ASD diagnosis in ADHD cases predicted cognitive performance in an independent developmental cohort. Our approach uncovered evidence of genetic heterogeneity in ADHD, helping us to understand its etiology and providing a model for studies of other disorders.

FEMA: Fast and efficient mixed-effects algorithm for large sample whole-brain imaging data

The linear mixed-effects model (LME) is a versatile approach to account for dependence among observations. Many large-scale neuroimaging datasets with complex designs have increased the need for LME; however LME has seldom been used in whole-brain imaging analyses due to its heavy computational requirements. In this paper, we introduce a fast and efficient mixed-effects algorithm (FEMA) that makes whole-brain vertex-wise, voxel-wise, and connectome-wide LME analyses in large samples possible. We validate FEMA with extensive simulations, showing that the estimates of the fixed effects are equivalent to standard maximum likelihood estimates but obtained with orders of magnitude improvement in computational speed. We demonstrate the applicability of FEMA by studying the cross-sectional and longitudinal effects of age on region-of-interest level and vertex-wise cortical thickness, as well as connectome-wide functional connectivity values derived from resting state functional MRI, using longitudinal imaging data from the Adolescent Brain Cognitive DevelopmentSM Study release 4.0. Our analyses reveal distinct spatial patterns for the annualized changes in vertex-wise cortical thickness and connectome-wide connectivity values in early adolescence, highlighting a critical time of brain maturation. The simulations and application to real data show that FEMA enables advanced investigation of the relationships between large numbers of neuroimaging metrics and variables of interest while considering complex study designs, including repeated measures and family structures, in a fast and efficient manner. The source code for FEMA is available via: https://github.com/cmig-research-group/cmig_tools/.

Beyond the Global Brain Differences: Intraindividual Variability Differences in 1q21.1 Distal and 15q11.2 BP1-BP2 Deletion Carriers

BACKGROUND

Carriers of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants exhibit regional and global brain differences compared with noncarriers. However, interpreting regional differences is challenging if a global difference drives the regional brain differences. Intraindividual variability measures can be used to test for regional differences beyond global differences in brain structure.

METHODS

Magnetic resonance imaging data were used to obtain regional brain values for 1q21.1 distal deletion (n = 30) and duplication (n = 27) and 15q11.2 BP1-BP2 deletion (n = 170) and duplication (n = 243) carriers and matched noncarriers (n = 2350). Regional intra-deviation scores, i.e., the standardized difference between an individual's regional difference and global difference, were used to test for regional differences that diverge from the global difference.

RESULTS

For the 1q21.1 distal deletion carriers, cortical surface area for regions in the medial visual cortex, posterior cingulate, and temporal pole differed less and regions in the prefrontal and superior temporal cortex differed more than the global difference in cortical surface area. For the 15q11.2 BP1-BP2 deletion carriers, cortical thickness in regions in the medial visual cortex, auditory cortex, and temporal pole differed less and the prefrontal and somatosensory cortex differed more than the global difference in cortical thickness.

CONCLUSIONS

We find evidence for regional effects beyond differences in global brain measures in 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants. The results provide new insight into brain profiling of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants, with the potential to increase understanding of the mechanisms involved in altered neurodevelopment.

Improved multimodal prediction of progression from MCI to Alzheimer's disease combining genetics with quantitative brain MRI and cognitive measures

INTRODUCTION

There is a pressing need for non-invasive, cost-effective tools for early detection of Alzheimer's disease (AD).

METHODS

Using data from the Alzheimer's Disease Neuroimaging Initiative (ADNI), Cox proportional models were conducted to develop a multimodal hazard score (MHS) combining age, a polygenic hazard score (PHS), brain atrophy, and memory to predict conversion from mild cognitive impairment (MCI) to dementia. Power calculations estimated required clinical trial sample sizes after hypothetical enrichment using the MHS. Cox regression determined predicted age of onset for AD pathology from the PHS.

RESULTS

The MHS predicted conversion from MCI to dementia (hazard ratio for 80th versus 20th percentile: 27.03). Models suggest that application of the MHS could reduce clinical trial sample sizes by 67%. The PHS alone predicted age of onset of amyloid and tau.

DISCUSSION

The MHS may improve early detection of AD for use in memory clinics or for clinical trial enrichment.

HIGHLIGHTS

A multimodal hazard score (MHS) combined age, genetics, brain atrophy, and memory. The MHS predicted time to conversion from mild cognitive impairment to dementia. MHS reduced hypothetical Alzheimer's disease (AD) clinical trial sample sizes by 67%. A polygenic hazard score predicted age of onset of AD neuropathology.

Estimating the Total Variance Explained by Whole-Brain Imaging for Zero-inflated Outcomes

Zero-inflated outcomes are very common in behavioral data, particularly for responses to psychological questionnaires. Modeling these challenging distributions is further exacerbated by the absence of established statistical models capable of characterizing total signals attributed to whole-brain imaging features, making the accurate assessment of brain-behavior relationships particularly formidable. Given this critical need, 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 robustly estimates the fraction of variance explained (FVE) and the proportion of non-null effects from large-scale imaging data. In simulations, ZIV outperformed other linear prediction algorithms. Applying ZIV to data from one of the largest neuroimaging studies, the Adolescent Brain Cognitive Development SM (ABCD) Study, we found that whole-brain imaging features have a larger FVE for externalizing compared to internalizing behavior. We also demonstrate that the ZIV estimator, especially applied to focal sub-scales, can localize key neurocircuitry associated with human behavior.

Heritability Estimation of Cognitive Phenotypes in the ABCD Study® Using Mixed Models

Twin and family studies have historically aimed to partition phenotypic variance into components corresponding to additive genetic effects (A), common environment (C), and unique environment (E). Here we present the ACE Model and several extensions in the Adolescent Brain Cognitive Development℠ Study (ABCD Study®), employed using the new Fast Efficient Mixed Effects Analysis (FEMA) package. In the twin sub-sample (n = 924; 462 twin pairs), heritability estimates were similar to those reported by prior studies for height (twin heritability = 0.86) and cognition (twin heritability between 0.00 and 0.61), respectively. Incorporating SNP-derived genetic relatedness and using the full ABCD Study® sample (n = 9,742) led to narrower confidence intervals for all parameter estimates. By leveraging the sparse clustering method used by FEMA to handle genetic relatedness only for participants within families, we were able to take advantage of the diverse distribution of genetic relatedness within the ABCD Study® sample.

Genotype Data and Derived Genetic Instruments of Adolescent Brain Cognitive Development Study® for Better Understanding of Human Brain Development

The data release of Adolescent Brain Cognitive Development® (ABCD) Study represents an extensive resource for investigating factors relating to child development and mental wellbeing. The genotype data of ABCD has been used extensively in the context of genetic analysis, including genome-wide association studies and polygenic score predictions. However, there are unique opportunities provided by ABCD genetic data that have not yet been fully tapped. The diverse genomic variability, the enriched relatedness among ABCD subsets, and the longitudinal design of the ABCD challenge researchers to perform novel analyses to gain deeper insight into human brain development. Genetic instruments derived from the ABCD genetic data, such as genetic principal components, can help to better control confounds beyond the context of genetic analyses. To facilitate the use genomic information in the ABCD for inference, we here detail the processing procedures, quality controls, general characteristics, and the corresponding resources in the ABCD genotype data of release 4.0.

Comparing Pruning and Thresholding with Continuous Shrinkage Polygenic Score Methods in a Large Sample of Ancestrally Diverse Adolescents from the ABCD Study®

Using individuals' genetic data researchers can generate Polygenic Scores (PS) that are able to predict risk for diseases, variability in different behaviors as well as anthropomorphic measures. This is achieved by leveraging models learned from previously published large Genome-Wide Association Studies (GWASs) associating locations in the genome with a phenotype of interest. Previous GWASs have predominantly been performed in European ancestry individuals. This is of concern as PS generated in samples with a different ancestry to the original training GWAS have been shown to have lower performance and limited portability, and many efforts are now underway to collect genetic databases on individuals of diverse ancestries. In this study, we compare multiple methods of generating PS, including pruning and thresholding and Bayesian continuous shrinkage models, to determine which of them is best able to overcome these limitations. To do this we use the ABCD Study, a longitudinal cohort with deep phenotyping on individuals of diverse ancestry. We generate PS for anthropometric and psychiatric phenotypes using previously published GWAS summary statistics and examine their performance in three subsamples of ABCD: African ancestry individuals (n = 811), European ancestry Individuals (n = 6703), and admixed ancestry individuals (n = 3664). We find that the single ancestry continuous shrinkage method, PRScs (CS), and the multi ancestry meta method, PRScsx Meta (CSx Meta), show the best performance across ancestries and phenotypes.

Multi-ancestry phenome-wide association of complement component 4 variation with psychiatric and brain phenotypes in youth

BACKGROUND

Increased expression of the complement component 4A (C4A) gene is associated with a greater lifetime risk of schizophrenia. In the brain, C4A is involved in synaptic pruning; yet, it remains unclear the extent to which upregulation of C4A alters brain development or is associated with the risk for psychotic symptoms in childhood. Here, we perform a multi-ancestry phenome-wide association study in 7789 children aged 9-12 years to examine the relationship between genetically regulated expression (GREx) of C4A, childhood brain structure, cognition, and psychiatric symptoms.

RESULTS

While C4A GREx is not related to childhood psychotic experiences, cognition, or global measures of brain structure, it is associated with a localized reduction in regional surface area (SA) of the entorhinal cortex. Furthermore, we show that reduced entorhinal cortex SA at 9-10 years predicts a greater number and severity of psychosis-like events at 1-year and 2-year follow-up time points. We also demonstrate that the effects of C4A on the entorhinal cortex are independent of genome-wide polygenic risk for schizophrenia.

CONCLUSIONS

Our results suggest neurodevelopmental effects of C4A on childhood medial temporal lobe structure, which may serve as a biomarker for schizophrenia risk prior to symptom onset.

Unique prediction of developmental psychopathology from genetic and familial risk

BACKGROUND

Early detection is critical for easing the rising burden of psychiatric disorders. However, the specificity of psychopathological measurements and genetic predictors is unclear among youth.

METHODS

We measured associations between genetic risk for psychopathology (polygenic risk scores (PRS) and family history (FH) measures) and a wide range of behavioral measures in a large sample (n = 5,204) of early adolescent participants (9-11 years) from the Adolescent Brain and Cognitive Development StudySM . Associations were measured both with and without accounting for shared variance across measures of genetic risk.

RESULTS

When controlling for genetic risk for other psychiatric disorders, polygenic risk for problematic opioid use (POU) is uniquely associated with lower behavioral inhibition. Attention deficit hyperactivity disorder (ADHD), depression (DEP), and attempted suicide (SUIC) PRS shared many significant associations with externalizing, internalizing, and psychosis-related behaviors. However, when accounting for all measures of genetic and familial risk, these PRS also showed clear, unique patterns of association. Polygenic risk for ASD, BIP, and SCZ, and attempted suicide uniquely predicted variability in cognitive performance. FH accounted for unique variability in behavior above and beyond PRS and vice versa, with FH measures explaining a greater proportion of unique variability compared to the PRS.

CONCLUSION

Our results indicate that, among youth, many behaviors show shared genetic influences; however, there is also specificity in the profile of emerging psychopathologies for individuals with high genetic risk for particular disorders. This may be useful for quantifying early, differential risk for psychopathology in development.

Association of Genetic Variant Linked to Hemochromatosis With Brain Magnetic Resonance Imaging Measures of Iron and Movement Disorders

Importance

Hereditary hemochromatosis (HH) is an autosomal recessive genetic disorder that leads to iron overload. Conflicting results from previous research has led some to believe the brain is spared the toxic effects of iron in HH.

Objective

To test the association of the strongest genetic risk variant for HH on brainwide measures sensitive to iron deposition and the rates of movement disorders in a substantially larger sample than previous studies of its kind.

Design, Setting, and Participants

This cross-sectional retrospective study included participants from the UK Biobank, a population-based sample. Genotype, health record, and neuroimaging data were collected from January 2006 to May 2021. Data analysis was conducted from January 2021 to April 2022. Disorders tested included movement disorders (International Statistical Classification of Diseases and Related Health Problems, Tenth Revision [ICD-10], codes G20-G26), abnormalities of gait and mobility (ICD-10 codes R26), and other disorders of the nervous system (ICD-10 codes G90-G99).

Exposures

Homozygosity for p.C282Y, the largest known genetic risk factor for HH.

Main Outcomes and Measures

T2-weighted and T2* signal intensity from brain magnetic resonance imaging scans, measures sensitive to iron deposition, and clinical diagnosis of neurological disorders.

Results

The total cohort consisted of 488 288 individuals (264 719 female; ages 49-87 years, largely northern European ancestry), 2889 of whom were p.C282Y homozygotes. The neuroimaging analysis consisted of 836 individuals: 165 p.C282Y homozygotes (99 female) and 671 matched controls (399 female). A total of 206 individuals were excluded from analysis due to withdrawal of consent. Neuroimaging analysis showed that p.C282Y homozygosity was associated with decreased T2-weighted and T2* signal intensity in subcortical motor structures (basal ganglia, thalamus, red nucleus, and cerebellum; Cohen d >1) consistent with substantial iron deposition. Across the whole UK Biobank (2889 p.C282Y homozygotes, 485 399 controls), we found a significantly increased prevalence for movement disorders in male homozygotes (OR, 1.80; 95% CI, 1.28-2.55; P = .001) but not female individuals (OR, 1.09; 95% CI, 0.70-1.73; P = .69). Among the 31 p.C282Y male homozygotes with a movement disorder, only 10 had a concurrent HH diagnosis.

Conclusions and Relevance

These findings indicate increased iron deposition in subcortical motor circuits in p.C282Y homozygotes and confirm an increased association with movement disorders in male homozygotes. Early treatment in HH effectively prevents the negative consequences of iron overload in the liver and heart. Our work suggests that screening for p.C282Y homozygosity in high-risk individuals also has the potential to reduce brain iron accumulation and to reduce the risk of movement disorders among male individuals who are homozygous for this mutation.

Multivariate genome-wide association study on tissue-sensitive diffusion metrics highlights pathways that shape the human brain

The molecular determinants of tissue composition of the human brain remain largely unknown. Recent genome-wide association studies (GWAS) on this topic have had limited success due to methodological constraints. Here, we apply advanced whole-brain analyses on multi-shell diffusion imaging data and multivariate GWAS to two large scale imaging genetic datasets (UK Biobank and the Adolescent Brain Cognitive Development study) to identify and validate genetic association signals. We discover 503 unique genetic loci that have impact on multiple regions of human brain. Among them, more than 79% are validated in either of two large-scale independent imaging datasets. Key molecular pathways involved in axonal growth, astrocyte-mediated neuroinflammation, and synaptogenesis during development are found to significantly impact the measured variations in tissue-specific imaging features. Our results shed new light on the biological determinants of brain tissue composition and their potential overlap with the genetic basis of neuropsychiatric disorders.

The genetic architecture of human cortical folding

The folding of the human cerebral cortex is a highly genetically regulated process that allows for a much larger surface area to fit into the cranial vault and optimizes functional organization. Sulcal depth is a robust yet understudied measure of localized folding, previously associated with multiple neurodevelopmental disorders. Here, we report the first genome-wide association study of sulcal depth. Through the multivariate omnibus statistical test (MOSTest) applied to vertex-wise measures from 33,748 U.K. Biobank participants (mean age, 64.3 years; 52.0% female), we identified 856 genome-wide significant loci (P < 5 × 10⁻⁸). Comparisons with cortical thickness and surface area indicated that sulcal depth has higher locus yield, heritability, and effective sample size. There was a large amount of genetic overlap between these traits, with gene-based analyses indicating strong associations with neurodevelopmental processes. Our findings demonstrate sulcal depth is a promising neuroimaging phenotype that may enhance our understanding of cortical morphology.

Vertex-wise multivariate genome-wide association study identifies 780 unique genetic loci associated with cortical morphology

Brain morphology has been shown to be highly heritable, yet only a small portion of the heritability is explained by the genetic variants discovered so far. Here we extended the Multivariate Omnibus Statistical Test (MOSTest) and applied it to genome-wide association studies (GWAS) of vertex-wise structural magnetic resonance imaging (MRI) cortical measures from N=35,657 participants in the UK Biobank. We identified 695 loci for cortical surface area and 539 for cortical thickness, in total 780 unique genetic loci associated with cortical morphology robustly replicated in 8,060 children of mixed ethnicity from the Adolescent Brain Cognitive Development (ABCD) Study®. This reflects more than 8-fold increase in genetic discovery at no cost to generalizability compared to the commonly used univariate GWAS methods applied to region of interest (ROI) data. Functional follow up including gene-based analyses implicated 10% of all protein-coding genes and pointed towards pathways involved in neurogenesis and cell differentiation. Power analysis indicated that applying the MOSTest to vertex-wise structural MRI data triples the effective sample size compared to conventional univariate GWAS approaches. The large boost in power obtained with the vertex-wise MOSTest together with pronounced replication rates and highlighted biologically meaningful pathways underscores the advantage of multivariate approaches in the context of highly distributed polygenic architecture of the human brain.

Distinct Regionalization Patterns of Cortical Morphology are Associated with Cognitive Performance Across Different Domains

Cognitive performance in children is predictive of academic and social outcomes; therefore, understanding neurobiological mechanisms underlying individual differences in cognition during development may be important for improving quality of life. The belief that a single, psychological construct underlies many cognitive processes is pervasive throughout society. However, it is unclear if there is a consistent neural substrate underlying many cognitive processes. Here, we show that a distributed configuration of cortical surface area and apparent thickness, when controlling for global imaging measures, is differentially associated with cognitive performance on different types of tasks in a large sample (N = 10 145) of 9-11-year-old children from the Adolescent Brain and Cognitive DevelopmentSM (ABCD) study. The minimal overlap in these regionalization patterns of association has implications for competing theories about developing intellectual functions. Surprisingly, not controlling for sociodemographic factors increased the similarity between these regionalization patterns. This highlights the importance of understanding the shared variance between sociodemographic factors, cognition and brain structure, particularly with a population-based sample such as ABCD.

Generalizing post-stroke prognoses from research data to clinical data

•

Lesion-symptom models for aphasia must generalize, if they are to be useful.

•

We show that these models can generalize across countries, native languages, neuroimaging technology and from chronic to acute scans.

•

The result depends on modelling lesion growth over years post-stroke.

•

Lesion growth may be a more significant confound, in previous lesion-symptom analyses, than previously thought.

Around a third of stroke survivors suffer from acquired language disorders (aphasia), but current medicine cannot predict whether or when they might recover. Prognostic research in this area increasingly draws on datasets associating structural brain imaging data with outcome scores for ever-larger samples of stroke patients. The aim is to learn brain-behaviour trends from these data, and generalize those trends to predict outcomes for new patients. The practical significance of this work depends on the expected breadth of that generalization. Here, we show that these models can generalize across countries and native languages (from British patients tested in English to Chilean patients tested in Spanish), across neuroimaging technology (from MRI to CT), and from scans collected months or years after stroke for research purposes, to scans collected days or weeks after stroke for clinical purposes.

Differences in mitochondrial DNA inheritance and function align with body conformation in genetically lean and fat sheep

Body weight and adiposity are determined by the balance between energy intake, energy expenditure, and nutrient deposition. We have identified differences in appetite-regulating peptides in sheep selectively bred to be either lean or fat, wherein gene expression for orexin and melanin-concentrating hormone are elevated in the lean group. Despite this, the underlying mechanisms leading to differences in body composition in the lean and fat lines remains unknown. We measured postprandial temperature in adipose tissue and muscle to ascertain whether a difference in thermogenesis is associated with the difference in body composition in genetically lean (n = 8) and fat (n = 12) ewes. Body weight was higher (P < 0.01) but percent fat mass was lower (P < 0.001) in the lean group. The percent lean mass was similar in lean and fat groups. Animals received intracerebroventricular cannulae and temperature probes implanted into the retroperitoneal fat and the hind-limb skeletal muscle (vastus lateralis). Animals were meal fed (1100-1600 h) to entrain postprandial thermogenesis. Food intake was similar between lean and fat animals. Postprandial thermogenesis was greater (P < 0.05) in the retroperitoneal adipose tissue of lean animals but not in skeletal muscle. Intracerebroventricular infusion of leptin reduced (P< 0.05) food intake by an equal extent in both groups. Postprandial expression of UCP1 mRNA was greater (P < 0.05) in retroperitoneal fat of lean animals, with similar UCP3 expression in skeletal muscle. Mitochondrial genome sequencing indicated haplotypic clustering in lean and fat animals within both the encoding and nonencoding regions. This demonstrates that differences in body composition may be underpinned by differences in thermogenesis, specifically within adipose tissue. Furthermore, thermogenic differences may be associated with specific mitochondrial DNA haplotypes, suggesting a strong genetic component inherited through the maternal lineage.