Principles and methods for transferring polygenic risk scores across global populations.

Polygenic risk scores (PRSs) summarize the genetic predisposition of a complex human trait or disease and may become a valuable tool for advancing precision medicine. However, PRSs that are developed in populations of predominantly European genetic ancestries can increase health disparities due to poor predictive performance in individuals of diverse and complex genetic ancestries. We describe genetic and modifiable risk factors that limit the transferability of PRSs across populations and review the strengths and weaknesses of existing PRS construction methods for diverse ancestries. Developing PRSs that benefit global populations in research and clinical settings provides an opportunity for innovation and is essential for health equity.

Mapping genomic loci implicates genes and synaptic biology in schizophrenia.

Schizophrenia has a heritability of 60-80%1, much of which is attributable to common risk alleles. Here, in a two-stage genome-wide association study of up to 76,755 individuals with schizophrenia and 243,649 control individuals, we report common variant associations at 287 distinct genomic loci. Associations were concentrated in genes that are expressed in excitatory and inhibitory neurons of the central nervous system, but not in other tissues or cell types. Using fine-mapping and functional genomic data, we identify 120 genes (106 protein-coding) that are likely to underpin associations at some of these loci, including 16 genes with credible causal non-synonymous or untranslated region variation. We also implicate fundamental processes related to neuronal function, including synaptic organization, differentiation and transmission. Fine-mapped candidates were enriched for genes associated with rare disruptive coding variants in people with schizophrenia, including the glutamate receptor subunit GRIN2A and transcription factor SP4, and were also enriched for genes implicated by such variants in neurodevelopmental disorders. We identify biological processes relevant to schizophrenia pathophysiology; show convergence of common and rare variant associations in schizophrenia and neurodevelopmental disorders; and provide a resource of prioritized genes and variants to advance mechanistic studies.

Population heterogeneity in clinical cohorts affects the predictive accuracy of brain imaging.

Brain imaging research enjoys increasing adoption of supervised machine learning for single-participant disease classification. Yet, the success of these algorithms likely depends on population diversity, including demographic differences and other factors that may be outside of primary scientific interest. Here, we capitalize on propensity scores as a composite confound index to quantify diversity due to major sources of population variation. We delineate the impact of population heterogeneity on the predictive accuracy and pattern stability in 2 separate clinical cohorts: the Autism Brain Imaging Data Exchange (ABIDE, n = 297) and the Healthy Brain Network (HBN, n = 551). Across various analysis scenarios, our results uncover the extent to which cross-validated prediction performances are interlocked with diversity. The instability of extracted brain patterns attributable to diversity is located preferentially in regions part of the default mode network. Collectively, our findings highlight the limitations of prevailing deconfounding practices in mitigating the full consequences of population diversity.

Heritability of individualized cortical network topography.

Human cortex is patterned by a complex and interdigitated web of large-scale functional networks. Recent methodological breakthroughs reveal variation in the size, shape, and spatial topography of cortical networks across individuals. While spatial network organization emerges across development, is stable over time, and is predictive of behavior, it is not yet clear to what extent genetic factors underlie interindividual differences in network topography. Here, leveraging a nonlinear multidimensional estimation of heritability, we provide evidence that individual variability in the size and topographic organization of cortical networks are under genetic control. Using twin and family data from the Human Connectome Project (n = 1,023), we find increased variability and reduced heritability in the size of heteromodal association networks (h2 : M = 0.34, SD = 0.070), relative to unimodal sensory/motor cortex (h2 : M = 0.40, SD = 0.097). We then demonstrate that the spatial layout of cortical networks is influenced by genetics, using our multidimensional estimation of heritability (h2-multi; M = 0.14, SD = 0.015). However, topographic heritability did not differ between heteromodal and unimodal networks. Genetic factors had a regionally variable influence on brain organization, such that the heritability of network topography was greatest in prefrontal, precuneus, and posterior parietal cortex. Taken together, these data are consistent with relaxed genetic control of association cortices relative to primary sensory/motor regions and have implications for understanding population-level variability in brain functioning, guiding both individualized prediction and the interpretation of analyses that integrate genetics and neuroimaging.

Impact of selection bias on polygenic risk score estimates in healthcare settings.

BACKGROUND: Hospital-based biobanks are being increasingly considered as a resource for translating polygenic risk scores (PRS) into clinical practice. However, since these biobanks originate from patient populations, there is a possibility of bias in polygenic risk estimation due to overrepresentation of patients with higher frequency of healthcare interactions. METHODS: PRS for schizophrenia, bipolar disorder, and depression were calculated using summary statistics from the largest available genomic studies for a sample of 24 153 European ancestry participants in the Mass General Brigham (MGB) Biobank. To correct for selection bias, we fitted logistic regression models with inverse probability (IP) weights, which were estimated using 1839 sociodemographic, clinical, and healthcare utilization features extracted from electronic health records of 1 546 440 non-Hispanic White patients eligible to participate in the Biobank study at their first visit to the MGB-affiliated hospitals. RESULTS: Case prevalence of bipolar disorder among participants in the top decile of bipolar disorder PRS was 10.0% (95% CI 8.8-11.2%) in the unweighted analysis but only 6.2% (5.0-7.5%) when selection bias was accounted for using IP weights. Similarly, case prevalence of depression among those in the top decile of depression PRS was reduced from 33.5% (31.7-35.4%) to 28.9% (25.8-31.9%) after IP weighting. CONCLUSIONS: Non-random selection of participants into volunteer biobanks may induce clinically relevant selection bias that could impact implementation of PRS in research and clinical settings. As efforts to integrate PRS in medical practice expand, recognition and mitigation of these biases should be considered and may need to be optimized in a context-specific manner.

Age differences in the functional architecture of the human brain.

The intrinsic functional organization of the brain changes into older adulthood. Age differences are observed at multiple spatial scales, from global reductions in modularity and segregation of distributed brain systems, to network-specific patterns of dedifferentiation. Whether dedifferentiation reflects an inevitable, global shift in brain function with age, circumscribed, experience-dependent changes, or both, is uncertain. We employed a multimethod strategy to interrogate dedifferentiation at multiple spatial scales. Multi-echo (ME) resting-state fMRI was collected in younger (n = 181) and older (n = 120) healthy adults. Cortical parcellation sensitive to individual variation was implemented for precision functional mapping of each participant while preserving group-level parcel and network labels. ME-fMRI processing and gradient mapping identified global and macroscale network differences. Multivariate functional connectivity methods tested for microscale, edge-level differences. Older adults had lower BOLD signal dimensionality, consistent with global network dedifferentiation. Gradients were largely age-invariant. Edge-level analyses revealed discrete, network-specific dedifferentiation patterns in older adults. Visual and somatosensory regions were more integrated within the functional connectome; default and frontoparietal control network regions showed greater connectivity; and the dorsal attention network was more integrated with heteromodal regions. These findings highlight the importance of multiscale, multimethod approaches to characterize the architecture of functional brain aging.

Effects of copy number variations on brain structure and risk for psychiatric illness: Large-scale studies from the ENIGMA working groups on CNVs.

The Enhancing NeuroImaging Genetics through Meta-Analysis copy number variant (ENIGMA-CNV) and 22q11.2 Deletion Syndrome Working Groups (22q-ENIGMA WGs) were created to gain insight into the involvement of genetic factors in human brain development and related cognitive, psychiatric and behavioral manifestations. To that end, the ENIGMA-CNV WG has collated CNV and magnetic resonance imaging (MRI) data from ~49,000 individuals across 38 global research sites, yielding one of the largest studies to date on the effects of CNVs on brain structures in the general population. The 22q-ENIGMA WG includes 12 international research centers that assessed over 533 individuals with a confirmed 22q11.2 deletion syndrome, 40 with 22q11.2 duplications, and 333 typically developing controls, creating the largest-ever 22q11.2 CNV neuroimaging data set. In this review, we outline the ENIGMA infrastructure and procedures for multi-site analysis of CNVs and MRI data. So far, ENIGMA has identified effects of the 22q11.2, 16p11.2 distal, 15q11.2, and 1q21.1 distal CNVs on subcortical and cortical brain structures. Each CNV is associated with differences in cognitive, neurodevelopmental and neuropsychiatric traits, with characteristic patterns of brain structural abnormalities. Evidence of gene-dosage effects on distinct brain regions also emerged, providing further insight into genotype-phenotype relationships. Taken together, these results offer a more comprehensive picture of molecular mechanisms involved in typical and atypical brain development. This "genotype-first" approach also contributes to our understanding of the etiopathogenesis of brain disorders. Finally, we outline future directions to better understand effects of CNVs on brain structure and behavior.

Integrative analysis of genomic and exposomic influences on youth mental health.

BACKGROUND: Understanding complex influences on mental health problems in young people is needed to inform early prevention strategies. Both genetic and environmental factors are known to influence youth mental health, but a more comprehensive picture of their interplay, including wide-ranging environmental exposures - that is, the exposome - is needed. We perform an integrative analysis of genomic and exposomic data in relation to internalizing and externalizing symptoms in a cohort of 4,314 unrelated youth from the Adolescent Brain and Cognitive Development (ABCD) Study. METHODS: Using novel GREML-based approaches, we model the variance in internalizing and externalizing symptoms explained by additive and interactive influences from the genome (G) and modeled exposome (E) consisting of up to 133 variables at the family, peer, school, neighborhood, life event, and broader environmental levels, including genome-by-exposome (G × E) and exposome-by-exposome (E × E) effects. RESULTS: A best-fitting integrative model with G, E, and G × E components explained 35% and 63% of variance in youth internalizing and externalizing symptoms, respectively. Youth in the top quintile of model-predicted risk accounted for the majority of individuals with clinically elevated symptoms at follow-up (60% for internalizing; 72% for externalizing). Of note, different domains of environmental exposures were most impactful for internalizing (life events) and externalizing (contextual including family, school, and peer-level factors) symptoms. In addition, variance explained by G × E contributions was substantially larger for externalizing (33%) than internalizing (13%) symptoms. CONCLUSIONS: Advanced statistical genetic methods in a longitudinal cohort of youth can be leveraged to address fundamental questions about the role of 'nature and nurture' in developmental psychopathology.

Individual-Specific Areal-Level Parcellations Improve Functional Connectivity Prediction of Behavior.

Resting-state functional magnetic resonance imaging (rs-fMRI) allows estimation of individual-specific cortical parcellations. We have previously developed a multi-session hierarchical Bayesian model (MS-HBM) for estimating high-quality individual-specific network-level parcellations. Here, we extend the model to estimate individual-specific areal-level parcellations. While network-level parcellations comprise spatially distributed networks spanning the cortex, the consensus is that areal-level parcels should be spatially localized, that is, should not span multiple lobes. There is disagreement about whether areal-level parcels should be strictly contiguous or comprise multiple noncontiguous components; therefore, we considered three areal-level MS-HBM variants spanning these range of possibilities. Individual-specific MS-HBM parcellations estimated using 10 min of data generalized better than other approaches using 150 min of data to out-of-sample rs-fMRI and task-fMRI from the same individuals. Resting-state functional connectivity derived from MS-HBM parcellations also achieved the best behavioral prediction performance. Among the three MS-HBM variants, the strictly contiguous MS-HBM exhibited the best resting-state homogeneity and most uniform within-parcel task activation. In terms of behavioral prediction, the gradient-infused MS-HBM was numerically the best, but differences among MS-HBM variants were not statistically significant. Overall, these results suggest that areal-level MS-HBMs can capture behaviorally meaningful individual-specific parcellation features beyond group-level parcellations. Multi-resolution trained models and parcellations are publicly available (https://github.com/ThomasYeoLab/CBIG/tree/master/stable_projects/brain_parcellation/Kong2022_ArealMSHBM).

Correction: Phenome-wide heritability analysis of the UK Biobank.

[This corrects the article DOI: 10.1371/journal.pgen.1006711.].

Polygenic risk for major depression is associated with lifetime suicide attempt in US soldiers independent of personal and parental history of major depression.

Suicide is a major public health problem. The contribution of common genetic variants for major depressive disorder (MDD) independent of personal and parental history of MDD has not been established. Polygenic risk score (using PRS-CS) for MDD was calculated for US Army soldiers of European ancestry. Associations between polygenic risk for MDD and lifetime suicide attempt (SA) were tested in models that also included parental or personal history of MDD. Models were adjusted for age, sex, tranche (where applicable), and 10 principal components reflecting ancestry. In the first cohort, 417 (6.3%) of 6,573 soldiers reported a lifetime history of SA. In a multivariable model that included personal [OR = 3.83, 95% CI:3.09-4.75] and parental history of MDD [OR = 1.43, 95% CI:1.13-1.82 for one parent and OR = 1.64, 95% CI:1.20-2.26 for both parents), MDD PRS was significantly associated with SA (OR = 1.22 [95% CI:1.10-1.36]). In the second cohort, 204 (4.2%) of 4,900 soldiers reported a lifetime history of SA. In a multivariable model that included personal [OR = 3.82, 95% CI:2.77-5.26] and parental history of MDD [OR = 1.42, 95% CI:0.996-2.03 for one parent and OR = 2.21, 95% CI:1.33-3.69 for both parents) MDD PRS continued to be associated (at p = .0601) with SA (OR = 1.15 [95% CI:0.994-1.33]). A soldier's PRS for MDD conveys information about likelihood of a lifetime SA beyond that conveyed by two predictors readily obtainable by interview: personal or parental history of MDD. Results remain to be extended to prospective prediction of incident SA. These findings portend a role for PRS in risk stratification for suicide attempts.

The Shared Genetic Basis of Educational Attainment and Cerebral Cortical Morphology.

Individual differences in educational attainment are linked to differences in intelligence, and predict important social, economic, and health outcomes. Previous studies have found common genetic factors that influence educational achievement, cognitive performance and total brain volume (i.e., brain size). Here, in a large sample of participants from the UK Biobank, we investigate the shared genetic basis between educational attainment and fine-grained cerebral cortical morphological features, and associate this genetic variation with a related aspect of cognitive ability. Importantly, we execute novel statistical methods that enable high-dimensional genetic correlation analysis, and compute high-resolution surface maps for the genetic correlations between educational attainment and vertex-wise morphological measurements. We conduct secondary analyses, using the UK Biobank verbal-numerical reasoning score, to confirm that variation in educational attainment that is genetically correlated with cortical morphology is related to differences in cognitive performance. Our analyses relate the genetic overlap between cognitive ability and cortical thickness measurements to bilateral primary motor cortex as well as predominantly left superior temporal cortex and proximal regions. These findings extend our understanding of the neurobiology that connects genetic variation to individual differences in educational attainment and cognitive performance.

A Polygenic Score for Higher Educational Attainment is Associated with Larger Brains.

People who score higher on intelligence tests tend to have larger brains. Twin studies suggest the same genetic factors influence both brain size and intelligence. This has led to the hypothesis that genetics influence intelligence partly by contributing to the development of larger brains. We tested this hypothesis using four large imaging genetics studies (combined N = 7965) with polygenic scores derived from a genome-wide association study (GWAS) of educational attainment, a correlate of intelligence. We conducted meta-analysis to test associations among participants' genetics, total brain volume (i.e., brain size), and cognitive test performance. Consistent with previous findings, participants with higher polygenic scores achieved higher scores on cognitive tests, as did participants with larger brains. Participants with higher polygenic scores also had larger brains. We found some evidence that brain size partly mediated associations between participants' education polygenic scores and their cognitive test performance. Effect sizes were larger in the population-based samples than in the convenience-based samples. Recruitment and retention of population-representative samples should be a priority for neuroscience research. Findings suggest promise for studies integrating GWAS discoveries with brain imaging to understand neurobiology linking genetics with cognitive performance.

Heritability analysis with repeat measurements and its application to resting-state functional connectivity.

Heritability, defined as the proportion of phenotypic variation attributable to genetic variation, provides important information about the genetic basis of a trait. Existing heritability analysis methods do not discriminate between stable effects (e.g., due to the subject's unique environment) and transient effects, such as measurement error. This can lead to misleading assessments, particularly when comparing the heritability of traits that exhibit different levels of reliability. Here, we present a linear mixed effects model to conduct heritability analyses that explicitly accounts for intrasubject fluctuations (e.g., due to measurement noise or biological transients) using repeat measurements. We apply the proposed strategy to the analysis of resting-state fMRI measurements-a prototypic data modality that exhibits variable levels of test-retest reliability across space. Our results reveal that the stable components of functional connectivity within and across well-established large-scale brain networks can be considerably heritable. Furthermore, we demonstrate that dissociating intra- and intersubject variation can reveal genetic influence on a phenotype that is not fully captured by conventional heritability analyses.

Phenome-wide heritability analysis of the UK Biobank.

Heritability estimation provides important information about the relative contribution of genetic and environmental factors to phenotypic variation, and provides an upper bound for the utility of genetic risk prediction models. Recent technological and statistical advances have enabled the estimation of additive heritability attributable to common genetic variants (SNP heritability) across a broad phenotypic spectrum. Here, we present a computationally and memory efficient heritability estimation method that can handle large sample sizes, and report the SNP heritability for 551 complex traits derived from the interim data release (152,736 subjects) of the large-scale, population-based UK Biobank, comprising both quantitative phenotypes and disease codes. We demonstrate that common genetic variation contributes to a broad array of quantitative traits and human diseases in the UK population, and identify phenotypes whose heritability is moderated by age (e.g., a majority of physical measures including height and body mass index), sex (e.g., blood pressure related traits) and socioeconomic status (education). Our study represents the first comprehensive phenome-wide heritability analysis in the UK Biobank, and underscores the importance of considering population characteristics in interpreting heritability.

Global signal regression strengthens association between resting-state functional connectivity and behavior.

Global signal regression (GSR) is one of the most debated preprocessing strategies for resting-state functional MRI. GSR effectively removes global artifacts driven by motion and respiration, but also discards globally distributed neural information and introduces negative correlations between certain brain regions. The vast majority of previous studies have focused on the effectiveness of GSR in removing imaging artifacts, as well as its potential biases. Given the growing interest in functional connectivity fingerprinting, here we considered the utilitarian question of whether GSR strengthens or weakens associations between resting-state functional connectivity (RSFC) and multiple behavioral measures across cognition, personality and emotion. By applying the variance component model to the Brain Genomics Superstruct Project (GSP), we found that behavioral variance explained by whole-brain RSFC increased by an average of 47% across 23 behavioral measures after GSR. In the Human Connectome Project (HCP), we found that behavioral variance explained by whole-brain RSFC increased by an average of 40% across 58 behavioral measures, when GSR was applied after ICA-FIX de-noising. To ensure generalizability, we repeated our analyses using kernel regression. GSR improved behavioral prediction accuracies by an average of 64% and 12% in the GSP and HCP datasets respectively. Importantly, the results were consistent across methods. A behavioral measure with greater RSFC-explained variance (using the variance component model) also exhibited greater prediction accuracy (using kernel regression). A behavioral measure with greater improvement in behavioral variance explained after GSR (using the variance component model) also enjoyed greater improvement in prediction accuracy after GSR (using kernel regression). Furthermore, GSR appeared to benefit task performance measures more than self-reported measures. Since GSR was more effective at removing motion-related and respiratory-related artifacts, GSR-related increases in variance explained and prediction accuracies were unlikely the result of motion-related or respiratory-related artifacts. However, it is worth emphasizing that the current study focused on whole-brain RSFC, so it remains unclear whether GSR improves RSFC-behavioral associations for specific connections or networks. Overall, our results suggest that at least in the case for young healthy adults, GSR strengthens the associations between RSFC and most (although not all) behavioral measures. Code for the variance component model and ridge regression can be found here: https://github.com/ThomasYeoLab/CBIG/tree/master/stable_projects/preprocessing/Li2019_GSR.

Accelerated estimation and permutation inference for ACE modeling.

There are a wealth of tools for fitting linear models at each location in the brain in neuroimaging analysis, and a wealth of genetic tools for estimating heritability for a small number of phenotypes. But there remains a need for computationally efficient neuroimaging genetic tools that can conduct analyses at the brain-wide scale. Here we present a simple method for heritability estimation on twins that replaces a variance component model-which requires iterative optimisation-with a (noniterative) linear regression model, by transforming data to squared twin-pair differences. We demonstrate that the method has comparable bias, mean squared error, false positive risk, and power to best practice maximum-likelihood-based methods, while requiring a small fraction of the computation time. Combined with permutation, we call this approach "Accelerated Permutation Inference for the ACE Model (APACE)" where ACE refers to the additive genetic (A) effects, and common (C), and unique (E) environmental influences on the trait. We show how the use of spatial statistics like cluster size can dramatically improve power, and illustrate the method on a heritability analysis of an fMRI working memory dataset.

Transcriptional profiles of supragranular-enriched genes associate with corticocortical network architecture in the human brain.

The human brain is patterned with disproportionately large, distributed cerebral networks that connect multiple association zones in the frontal, temporal, and parietal lobes. The expansion of the cortical surface, along with the emergence of long-range connectivity networks, may be reflected in changes to the underlying molecular architecture. Using the Allen Institute's human brain transcriptional atlas, we demonstrate that genes particularly enriched in supragranular layers of the human cerebral cortex relative to mouse distinguish major cortical classes. The topography of transcriptional expression reflects large-scale brain network organization consistent with estimates from functional connectivity MRI and anatomical tracing in nonhuman primates. Microarray expression data for genes preferentially expressed in human upper layers (II/III), but enriched only in lower layers (V/VI) of mouse, were cross-correlated to identify molecular profiles across the cerebral cortex of postmortem human brains (n = 6). Unimodal sensory and motor zones have similar molecular profiles, despite being distributed across the cortical mantle. Sensory/motor profiles were anticorrelated with paralimbic and certain distributed association network profiles. Tests of alternative gene sets did not consistently distinguish sensory and motor regions from paralimbic and association regions: (i) genes enriched in supragranular layers in both humans and mice, (ii) genes cortically enriched in humans relative to nonhuman primates, (iii) genes related to connectivity in rodents, (iv) genes associated with human and mouse connectivity, and (v) 1,454 gene sets curated from known gene ontologies. Molecular innovations of upper cortical layers may be an important component in the evolution of long-range corticocortical projections.

Morphometricity as a measure of the neuroanatomical signature of a trait.

Complex physiological and behavioral traits, including neurological and psychiatric disorders, often associate with distributed anatomical variation. This paper introduces a global metric, called morphometricity, as a measure of the anatomical signature of different traits. Morphometricity is defined as the proportion of phenotypic variation that can be explained by macroscopic brain morphology. We estimate morphometricity via a linear mixed-effects model that uses an anatomical similarity matrix computed based on measurements derived from structural brain MRI scans. We examined over 3,800 unique MRI scans from nine large-scale studies to estimate the morphometricity of a range of phenotypes, including clinical diagnoses such as Alzheimer's disease, and nonclinical traits such as measures of cognition. Our results demonstrate that morphometricity can provide novel insights about the neuroanatomical correlates of a diverse set of traits, revealing associations that might not be detectable through traditional statistical techniques.

Effects of Angiotensin Receptor Neprilysin Inhibitors on Inducibility of Ventricular Arrhythmias in Rats with Ischemic Cardiomyopathy.

The aims of the present study were to investigate the effects of angiotensin receptor neprilysin inhibitors (ARNi) on the susceptibility of ventricular arrhythmias (VAs) in rats with myocardial infarction (MI) and to explore the related mechanisms.A total of 32 adult male Sprague-Dawley rats were divided into 3 groups: a control group, MI group, and MI+ARNi group. MI was generated by ligation of the left anterior descending coronary artery. ARNi was given at 68 mg/kg/day for 4 weeks after MI surgery. At 4 weeks after MI, electrical programmed stimulation (EPS) was performed in all groups for the evaluation of VAs, and echocardiography was used to evaluate cardiac function. Indicators of sympathetic neural remodeling and cardiac remodeling were detected to further explore the related mechanisms.Four weeks after MI, rats in the ARNi group exhibited low susceptibility of VAs in comparison with that in the MI group, which was coincident with the attenuation of sympathetic nerve remodeling, amelioration of cardiac fibrosis, and regulation of Cx43 expression.ARNi is effective in reducing VAs in rats with ischemic cardiomyopathy, which is associated with attenuating sympathetic nerve remodeling and myocardial fibrosis.