Partitioning variance in cortical morphometry into genetic, environmental, and subject-specific components.

The relative contributions of genetic variation and experience in shaping the morphology of the adolescent brain are not fully understood. Using longitudinal data from 11,665 subjects in the ABCD Study, we fit vertex-wise variance components including family effects, genetic effects, and subject-level effects using a computationally efficient framework. Variance in cortical thickness and surface area is largely attributable to genetic influence, whereas sulcal depth is primarily explained by subject-level effects. Our results identify areas with heterogeneous distributions of heritability estimates that have not been seen in previous work using data from cortical regions. We discuss the biological importance of subject-specific variance and its implications for environmental influences on cortical development and maturation.

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.

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/.

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.

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.

Intelligence Polygenic Score Is More Predictive of Crystallized Measures: Evidence From the Adolescent Brain Cognitive Development (ABCD) Study.

Findings in adults have shown that crystallized measures of intelligence, which are more culturally sensitive than fluid intelligence measures, have greater heritability; however, these results have not been found in children. The present study used data from 8,518 participants between 9 and 11 years old from the Adolescent Brain Cognitive Development (ABCD) Study. We found that polygenic predictors of intelligence test performance (based on genome-wide association meta-analyses of data from 269,867 individuals) and of educational attainment (based on data from 1.1 million individuals) predicted neurocognitive performance. We found that crystallized measures were more strongly associated with both polygenic predictors than were fluid measures. This mirrored heritability differences reported previously in adults and suggests similar associations in children. This may be consistent with a prominent role of gene-environment correlation in cognitive development measured by crystallized intelligence tests. Environmental and experiential mediators may represent malleable targets for improving cognitive outcomes.

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.

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.

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.

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.

TMEM161B modulates radial glial scaffolding in neocortical development.

TMEM161B encodes an evolutionarily conserved widely expressed novel 8-pass transmembrane protein of unknown function in human. Here we identify TMEM161B homozygous hypomorphic missense variants in our recessive polymicrogyria (PMG) cohort. Patients carrying TMEM161B mutations exhibit striking neocortical PMG and intellectual disability. Tmem161b knockout mice fail to develop midline hemispheric cleavage, whereas knock-in of patient mutations and patient-derived brain organoids show defects in apical cell polarity and radial glial scaffolding. We found that TMEM161B modulates actin filopodia, functioning upstream of the Rho-GTPase CDC42. Our data link TMEM161B with human PMG, likely regulating radial glia apical polarity during neocortical development.

Generalization of cortical MOSTest genome-wide associations within and across samples.

Genome-Wide Association studies have typically been limited to univariate analysis in which a single outcome measure is tested against millions of variants. Recent work demonstrates that a Multivariate Omnibus Statistic Test (MOSTest) is well powered to discover genomic effects distributed across multiple phenotypes. Applied to cortical brain MRI morphology measures, MOSTest has resulted in a drastic improvement in power to discover loci when compared to established approaches (min-P). One question that arises is how well these discovered loci replicate in independent data. Here we perform 10 times cross validation within 34,973 individuals from UK Biobank for imaging measures of cortical area, thickness and sulcal depth (>1,000 dimensionality for each). By deploying a replication method that aggregates discovered effects distributed across multiple phenotypes, termed PolyVertex Score (MOSTest-PVS), we demonstrate a higher replication yield and comparable replication rate of discovered loci for MOSTest (# replicated loci: 242-496, replication rate: 96-97%) in independent data when compared with the established min-P approach (# replicated loci: 26-55, replication rate: 91-93%). An out-of-sample replication of discovered loci was conducted with a sample of 4,069 individuals from the Adolescent Brain Cognitive Development® (ABCD) study, who are on average 50 years younger than UK Biobank individuals. We observe a higher replication yield and comparable replication rate of MOSTest-PVS compared to min-P. This finding underscores the importance of using well-powered multivariate techniques for both discovery and replication of high dimensional phenotypes in Genome-Wide Association studies.

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.

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.

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−8). 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.

Generalizing post-stroke prognoses from research data to clinical data.

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. Our results suggest one important confound, in attempting to generalize from research data to clinical data, is the delay between scan acquisition and language assessment. This delay is typically small for research data, where scans and assessments are often acquired contemporaneously. But the most natural, clinical application of these predictions will employ acute prognostic factors to predict much longer-term outcomes. We mitigated this confound by projecting the clinical patients' lesions from the time when their scans were acquired, to the time when their language abilities were assessed; with this projection in place, there was strong evidence that prognoses derived from research data generalized equally well to research and clinical data. These results encourage attention to the confounding role that lesion growth may play in other types of lesion-symptom analysis.