Modeling genotype-protein interaction and correlation for Alzheimer's disease: a multi-omics imaging genetics study.

Integrating and analyzing multiple omics data sets, including genomics, proteomics and radiomics, can significantly advance researchers' comprehensive understanding of Alzheimer's disease (AD). However, current methodologies primarily focus on the main effects of genetic variation and protein, overlooking non-additive effects such as genotype-protein interaction (GPI) and correlation patterns in brain imaging genetics studies. Importantly, these non-additive effects could contribute to intermediate imaging phenotypes, finally leading to disease occurrence. In general, the interaction between genetic variations and proteins, and their correlations are two distinct biological effects, and thus disentangling the two effects for heritable imaging phenotypes is of great interest and need. Unfortunately, this issue has been largely unexploited. In this paper, to fill this gap, we propose $\textbf{M}$ulti-$\textbf{T}$ask $\textbf{G}$enotype-$\textbf{P}$rotein $\textbf{I}$nteraction and $\textbf{C}$orrelation disentangling method ($\textbf{MT-GPIC}$) to identify GPI and extract correlation patterns between them. To ensure stability and interpretability, we use novel and off-the-shelf penalties to identify meaningful genetic risk factors, as well as exploit the interconnectedness of different brain regions. Additionally, since computing GPI poses a high computational burden, we develop a fast optimization strategy for solving MT-GPIC, which is guaranteed to converge. Experimental results on the Alzheimer's Disease Neuroimaging Initiative data set show that MT-GPIC achieves higher correlation coefficients and classification accuracy than state-of-the-art methods. Moreover, our approach could effectively identify interpretable phenotype-related GPI and correlation patterns in high-dimensional omics data sets. These findings not only enhance the diagnostic accuracy but also contribute valuable insights into the underlying pathogenic mechanisms of AD.

Genomic loci influence patterns of structural covariance in the human brain.

Normal and pathologic neurobiological processes influence brain morphology in coordinated ways that give rise to patterns of structural covariance (PSC) across brain regions and individuals during brain aging and diseases. The genetic underpinnings of these patterns remain largely unknown. We apply a stochastic multivariate factorization method to a diverse population of 50,699 individuals (12 studies and 130 sites) and derive data-driven, multi-scale PSCs of regional brain size. PSCs were significantly correlated with 915 genomic loci in the discovery set, 617 of which are newly identified, and 72% were independently replicated. Key pathways influencing PSCs involve reelin signaling, apoptosis, neurogenesis, and appendage development, while pathways of breast cancer indicate potential interplays between brain metastasis and PSCs associated with neurodegeneration and dementia. Using support vector machines, multi-scale PSCs effectively derive imaging signatures of several brain diseases. Our results elucidate genetic and biological underpinnings that influence structural covariance patterns in the human brain.

Bayesian mixed model inference for genetic association under related samples with brain network phenotype.

Genetic association studies for brain connectivity phenotypes have gained prominence due to advances in noninvasive imaging techniques and quantitative genetics. Brain connectivity traits, characterized by network configurations and unique biological structures, present distinct challenges compared to other quantitative phenotypes. Furthermore, the presence of sample relatedness in the most imaging genetics studies limits the feasibility of adopting existing network-response modeling. In this article, we fill this gap by proposing a Bayesian network-response mixed-effect model that considers a network-variate phenotype and incorporates population structures including pedigrees and unknown sample relatedness. To accommodate the inherent topological architecture associated with the genetic contributions to the phenotype, we model the effect components via a set of effect network configurations and impose an inter-network sparsity and intra-network shrinkage to dissect the phenotypic network configurations affected by the risk genetic variant. A Markov chain Monte Carlo (MCMC) algorithm is further developed to facilitate uncertainty quantification. We evaluate the performance of our model through extensive simulations. By further applying the method to study, the genetic bases for brain structural connectivity using data from the Human Connectome Project with excessive family structures, we obtain plausible and interpretable results. Beyond brain connectivity genetic studies, our proposed model also provides a general linear mixed-effect regression framework for network-variate outcomes.

Developmental dyslexia susceptibility genes DNAAF4, DCDC2, and NRSN1 are associated with brain function in fluently reading adolescents and young adults.

Reading skills and developmental dyslexia, characterized by difficulties in developing reading skills, have been associated with brain anomalies within the language network. Genetic factors contribute to developmental dyslexia risk, but the mechanisms by which these genes influence reading skills remain unclear. In this preregistered study (https://osf.io/7sehx), we explored if developmental dyslexia susceptibility genes DNAAF4, DCDC2, NRSN1, and KIAA0319 are associated with brain function in fluently reading adolescents and young adults. Functional MRI and task performance data were collected during tasks involving written and spoken sentence processing, and DNA sequence variants of developmental dyslexia susceptibility genes previously associated with brain structure anomalies were genotyped. The results revealed that variation in DNAAF4, DCDC2, and NRSN1 is associated with brain activity in key language regions: the left inferior frontal gyrus, middle temporal gyrus, and intraparietal sulcus. Furthermore, NRSN1 was associated with task performance, but KIAA0319 did not yield any significant associations. Our findings suggest that individuals with a genetic predisposition to developmental dyslexia may partly employ compensatory neural and behavioral mechanisms to maintain typical task performance. Our study highlights the relevance of these developmental dyslexia susceptibility genes in language-related brain function, even in individuals without developmental dyslexia, providing valuable insights into the genetic factors influencing language processing.

Feature aggregation graph convolutional network based on imaging genetic data for diagnosis and pathogeny identification of Alzheimer's disease.

The roles of brain regions activities and gene expressions in the development of Alzheimer's disease (AD) remain unclear. Existing imaging genetic studies usually has the problem of inefficiency and inadequate fusion of data. This study proposes a novel deep learning method to efficiently capture the development pattern of AD. First, we model the interaction between brain regions and genes as node-to-node feature aggregation in a brain region-gene network. Second, we propose a feature aggregation graph convolutional network (FAGCN) to transmit and update the node feature. Compared with the trivial graph convolutional procedure, we replace the input from the adjacency matrix with a weight matrix based on correlation analysis and consider common neighbor similarity to discover broader associations of nodes. Finally, we use a full-gradient saliency graph mechanism to score and extract the pathogenetic brain regions and risk genes. According to the results, FAGCN achieved the best performance among both traditional and cutting-edge methods and extracted AD-related brain regions and genes, providing theoretical and methodological support for the research of related diseases.

IHGC-GAN: influence hypergraph convolutional generative adversarial network for risk prediction of late mild cognitive impairment based on imaging genetic data.

Predicting disease progression in the initial stage to implement early intervention and treatment can effectively prevent the further deterioration of the condition. Traditional methods for medical data analysis usually fail to perform well because of their incapability for mining the correlation pattern of pathogenies. Therefore, many calculation methods have been excavated from the field of deep learning. In this study, we propose a novel method of influence hypergraph convolutional generative adversarial network (IHGC-GAN) for disease risk prediction. First, a hypergraph is constructed with genes and brain regions as nodes. Then, an influence transmission model is built to portray the associations between nodes and the transmission rule of disease information. Third, an IHGC-GAN method is constructed based on this model. This method innovatively combines the graph convolutional network (GCN) and GAN. The GCN is used as the generator in GAN to spread and update the lesion information of nodes in the brain region-gene hypergraph. Finally, the prediction accuracy of the method is improved by the mutual competition and repeated iteration between generator and discriminator. This method can not only capture the evolutionary pattern from early mild cognitive impairment (EMCI) to late MCI (LMCI) but also extract the pathogenic factors and predict the deterioration risk from EMCI to LMCI. The results on the two datasets indicate that the IHGC-GAN method has better prediction performance than the advanced methods in a variety of indicators.

A novel generation adversarial network framework with characteristics aggregation and diffusion for brain disease classification and feature selection.

Imaging genetics provides unique insights into the pathological studies of complex brain diseases by integrating the characteristics of multi-level medical data. However, most current imaging genetics research performs incomplete data fusion. Also, there is a lack of effective deep learning methods to analyze neuroimaging and genetic data jointly. Therefore, this paper first constructs the brain region-gene networks to intuitively represent the association pattern of pathogenetic factors. Second, a novel feature information aggregation model is constructed to accurately describe the information aggregation process among brain region nodes and gene nodes. Finally, a deep learning method called feature information aggregation and diffusion generative adversarial network (FIAD-GAN) is proposed to efficiently classify samples and select features. We focus on improving the generator with the proposed convolution and deconvolution operations, with which the interpretability of the deep learning framework has been dramatically improved. The experimental results indicate that FIAD-GAN can not only achieve superior results in various disease classification tasks but also extract brain regions and genes closely related to AD. This work provides a novel method for intelligent clinical decisions. The relevant biomedical discoveries provide a reliable reference and technical basis for the clinical diagnosis, treatment and pathological analysis of disease.

A multivariate to multivariate approach for voxel-wise genome-wide association analysis.

The joint analysis of imaging-genetics data facilitates the systematic investigation of genetic effects on brain structures and functions with spatial specificity. We focus on voxel-wise genome-wide association analysis, which may involve trillions of single nucleotide polymorphism (SNP)-voxel pairs. We attempt to identify underlying organized association patterns of SNP-voxel pairs and understand the polygenic and pleiotropic networks on brain imaging traits. We propose a bi-clique graph structure (ie, a set of SNPs highly correlated with a cluster of voxels) for the systematic association pattern. Next, we develop computational strategies to detect latent SNP-voxel bi-cliques and an inference model for statistical testing. We further provide theoretical results to guarantee the accuracy of our computational algorithms and statistical inference. We validate our method by extensive simulation studies, and then apply it to the whole genome genetic and voxel-level white matter integrity data collected from 1052 participants of the human connectome project. The results demonstrate multiple genetic loci influencing white matter integrity measures on splenium and genu of the corpus callosum.

Preference matrix guided sparse canonical correlation analysis for mining brain imaging genetic associations in Alzheimer's disease.

Investigating the relationship between genetic variation and phenotypic traits is a key issue in quantitative genetics. Specifically for Alzheimer's disease, the association between genetic markers and quantitative traits remains vague while, once identified, will provide valuable guidance for the study and development of genetics-based treatment approaches. Currently, to analyze the association of two modalities, sparse canonical correlation analysis (SCCA) is commonly used to compute one sparse linear combination of the variable features for each modality, giving a pair of linear combination vectors in total that maximizes the cross-correlation between the analyzed modalities. One drawback of the plain SCCA model is that the existing findings and knowledge cannot be integrated into the model as priors to help extract interesting correlations as well as identify biologically meaningful genetic and phenotypic markers. To bridge this gap, we introduce preference matrix guided SCCA (PM-SCCA) that not only takes priors encoded as a preference matrix but also maintains computational simplicity. A simulation study and a real-data experiment are conducted to investigate the effectiveness of the model. Both experiments demonstrate that the proposed PM-SCCA model can capture not only genotype-phenotype correlation but also relevant features effectively.

Multimodal mapping of regional brain vulnerability to focal cortical dysplasia.

Focal cortical dysplasia (FCD) type II is a highly epileptogenic developmental malformation and a common cause of surgically treated drug-resistant epilepsy. While clinical observations suggest frequent occurrence in the frontal lobe, mechanisms for such propensity remain unexplored. Here, we hypothesized that cortex-wide spatial associations of FCD distribution with cortical cytoarchitecture, gene expression and organizational axes may offer complementary insights into processes that predispose given cortical regions to harbour FCD. We mapped the cortex-wide MRI distribution of FCDs in 337 patients collected from 13 sites worldwide. We then determined its associations with (i) cytoarchitectural features using histological atlases by Von Economo and Koskinas and BigBrain; (ii) whole-brain gene expression and spatiotemporal dynamics from prenatal to adulthood stages using the Allen Human Brain Atlas and PsychENCODE BrainSpan; and (iii) macroscale developmental axes of cortical organization. FCD lesions were preferentially located in the prefrontal and fronto-limbic cortices typified by low neuron density, large soma and thick grey matter. Transcriptomic associations with FCD distribution uncovered a prenatal component related to neuroglial proliferation and differentiation, likely accounting for the dysplastic makeup, and a postnatal component related to synaptogenesis and circuit organization, possibly contributing to circuit-level hyperexcitability. FCD distribution showed a strong association with the anterior region of the antero-posterior axis derived from heritability analysis of interregional structural covariance of cortical thickness, but not with structural and functional hierarchical axes. Reliability of all results was confirmed through resampling techniques. Multimodal associations with cytoarchitecture, gene expression and axes of cortical organization indicate that prenatal neurogenesis and postnatal synaptogenesis may be key points of developmental vulnerability of the frontal lobe to FCD. Concordant with a causal role of atypical neuroglial proliferation and growth, our results indicate that FCD-vulnerable cortices display properties indicative of earlier termination of neurogenesis and initiation of cell growth. They also suggest a potential contribution of aberrant postnatal synaptogenesis and circuit development to FCD epileptogenicity.

Identifying the joint signature of brain atrophy and gene variant scores in Alzheimer's Disease.

The joint modeling of genetic data and brain imaging information allows for determining the pathophysiological pathways of neurodegenerative diseases such as Alzheimer's disease (AD). This task has typically been approached using mass-univariate methods that rely on a complete set of Single Nucleotide Polymorphisms (SNPs) to assess their association with selected image-derived phenotypes (IDPs). However, such methods are prone to multiple comparisons bias and, most importantly, fail to account for potential cross-feature interactions, resulting in insufficient detection of significant associations. Ways to overcome these limitations while reducing the number of traits aim at conveying genetic information at the gene level and capturing the integrated genetic effects of a set of genetic variants, rather than looking at each SNP individually. Their associations with brain IDPs are still largely unexplored in the current literature, though they can uncover new potential genetic determinants for brain modulations in the AD continuum. In this work, we explored an explainable multivariate model to analyze the genetic basis of the grey matter modulations, relying on the AD Neuroimaging Initiative (ADNI) phase 3 dataset. Cortical thicknesses and subcortical volumes derived from T1-weighted Magnetic Resonance were considered to describe the imaging phenotypes. At the same time the genetic counterpart was represented by gene variant scores extracted by the Sequence Kernel Association Test (SKAT) filtering model. Moreover, transcriptomic analysis was carried on to assess the expression of the resulting genes in the main brain structures as a form of validation. Results highlighted meaningful genotype-phenotype interactionsas defined by three latent components showing a significant difference in the projection scores between patients and controls. Among the significant associations, the model highlighted EPHX1 and BCAS1 gene variant scores involved in neurodegenerative and myelination processes, hence relevant for AD. In particular, the first was associated with decreased subcortical volumes and the second with decreasedtemporal lobe thickness. Noteworthy, BCAS1 is particularly expressed in the dentate gyrus. Overall, the proposed approach allowed capturing genotype-phenotype interactions in a restricted study cohort that was confirmed by transcriptomic analysis, offering insights into the underlying mechanisms of neurodegeneration in AD in line with previous findings and suggesting new potential disease biomarkers.

A brain-wide genome-wide association study of candidate quantitative trait loci associated with structural and functional phenotypes of pain sensitivity.

Individual pain sensitivity is modulated by the brain's structural and functional features, but its heritability remains unclear. This paper conducted a brain-wide genome-wide association study (GWAS) to explore the genetic bases of neuroimage phenotypes of pain sensitivity. In total, 432 normal participants were divided into high and low pain sensitivity groups according to the laser quantitative test threshold. Then, the brain's gray matter density (GMD) features correlated with pain sensitivity were identified. Next, GWAS was performed on each GMD phenotype using quality-controlled genotypes. Based on the heatmap and hierarchical clustering results, the right insula was identified for further refined analysis in terms of subregions GMD and resting-state functional connectivity (rs-FC) phenotypes. The results indicate that the right insula GMD in the high sensitivity group is significantly lower than that in the low sensitivity group. Also, the TT/TC group at locus rs187974 has lower right insula GMD than the CC group. Further, loci at gene CYP2D6 may lead to a variation of rs-FC between the right insula and left putamen. In conclusion, our study suggests that the right insula and multiple candidate loci may be importantly involved in pain sensitivity modulation, which may guide the future development of precision pain therapeutics.

Association of BDNF risk variant and dorsolateral cortical thickness with long-term treatment response to valproate in type I bipolar disorder: An exploratory study.

Valproate is among the most prescribed drugs for bipolar disorder; however, 87% of patients do not report full long-term treatment response (LTTR) to this medication. One of valproate's suggested mechanisms of action involves the brain-derived neurotrophic factor (BDNF), expressed in the brain areas regulating emotions, such as the prefrontal cortex. Nonetheless, data about the role of BDNF in LTTR and its implications in the structure of the dorsolateral prefrontal cortex (dlPFC) is scarce. We explore the association of BDNF variants and dorsolateral cortical thickness (CT) with LTTR to valproate in bipolar disorder type I (BDI). Twenty-eight BDI patients were genotyped for BDNF polymorphisms rs1519480, rs6265, and rs7124442, and T1-weighted 3D brain scans were acquired. LTTR to valproate was evaluated with Alda's scale. A logistic regression analysis was conducted to evaluate LTTR according to BDNF genotypes and CT. We evaluated CT differences by genotypes with analysis of covariance. LTTR was associated with BDNF rs1519480 and right dlPFC thickness. Insufficient responders with the CC genotype had thicker right dlPFC than TC and TT genotypes. Full responders reported thicker right dlPFC in TC and TT genotypes. In conclusion, different patterns of CT related to BDNF genotypes were identified, suggesting a potential biomarker of LTTR to valproate in our population.

Neuroimaging feature extraction using a neural network classifier for imaging genetics.

BACKGROUND: Dealing with the high dimension of both neuroimaging data and genetic data is a difficult problem in the association of genetic data to neuroimaging. In this article, we tackle the latter problem with an eye toward developing solutions that are relevant for disease prediction. Supported by a vast literature on the predictive power of neural networks, our proposed solution uses neural networks to extract from neuroimaging data features that are relevant for predicting Alzheimer's Disease (AD) for subsequent relation to genetics. The neuroimaging-genetic pipeline we propose is comprised of image processing, neuroimaging feature extraction and genetic association steps. We present a neural network classifier for extracting neuroimaging features that are related with the disease. The proposed method is data-driven and requires no expert advice or a priori selection of regions of interest. We further propose a multivariate regression with priors specified in the Bayesian framework that allows for group sparsity at multiple levels including SNPs and genes. RESULTS: We find the features extracted with our proposed method are better predictors of AD than features used previously in the literature suggesting that single nucleotide polymorphisms (SNPs) related to the features extracted by our proposed method are also more relevant for AD. Our neuroimaging-genetic pipeline lead to the identification of some overlapping and more importantly some different SNPs when compared to those identified with previously used features. CONCLUSIONS: The pipeline we propose combines machine learning and statistical methods to benefit from the strong predictive performance of blackbox models to extract relevant features while preserving the interpretation provided by Bayesian models for genetic association. Finally, we argue in favour of using automatic feature extraction, such as the method we propose, in addition to ROI or voxelwise analysis to find potentially novel disease-relevant SNPs that may not be detected when using ROIs or voxels alone.

Brain-wide genome-wide colocalization study for integrating genetics, transcriptomics and brain morphometry in Alzheimer's disease.

Alzheimer's disease (AD) is one of the most common neurodegenerative diseases. However, the AD mechanism has not yet been fully elucidated to date, hindering the development of effective therapies. In our work, we perform a brain imaging genomics study to link genetics, single-cell gene expression data, tissue-specific gene expression data, brain imaging-derived volumetric endophenotypes, and disease diagnosis to discover potential underlying neurobiological pathways for AD. To do so, we perform brain-wide genome-wide colocalization analyses to integrate multidimensional imaging genomic biobank data. Specifically, we use (1) the individual-level imputed genotyping data and magnetic resonance imaging (MRI) data from the UK Biobank, (2) the summary statistics of the genome-wide association study (GWAS) from multiple European ancestry cohorts, and (3) the tissue-specific cis-expression quantitative trait loci (cis-eQTL) summary statistics from the GTEx project. We apply a Bayes factor colocalization framework and mediation analysis to these multi-modal imaging genomic data. As a result, we derive the brain regional level GWAS summary statistics for 145 brain regions with 482,831 single nucleotide polymorphisms (SNPs) followed by posthoc functional annotations. Our analysis yields the discovery of a potential AD causal pathway from a systems biology perspective: the SNP chr10:124165615:G>A (rs6585827) mutation upregulates the expression of BTBD16 gene in oligodendrocytes, a specialized glial cells, in the brain cortex, leading to a reduced risk of volumetric loss in the entorhinal cortex, resulting in the protective effect on AD. We substantiate our findings with multiple evidence from existing imaging, genetic and genomic studies in AD literature. Our study connects genetics, molecular and cellular signatures, regional brain morphologic endophenotypes, and AD diagnosis, providing new insights into the mechanistic understanding of the disease. Our findings can provide valuable guidance for subsequent therapeutic target identification and drug discovery in AD.

Bayesian sparse heritability analysis with high-dimensional neuroimaging phenotypes.

Heritability analysis plays a central role in quantitative genetics to describe genetic contribution to human complex traits and prioritize downstream analyses under large-scale phenotypes. Existing works largely focus on modeling single phenotype and currently available multivariate phenotypic methods often suffer from scaling and interpretation. In this article, motivated by understanding how genetic underpinning impacts human brain variation, we develop an integrative Bayesian heritability analysis to jointly estimate heritabilities for high-dimensional neuroimaging traits. To induce sparsity and incorporate brain anatomical configuration, we impose hierarchical selection among both regional and local measurements based on brain structural network and voxel dependence. We also use a nonparametric Dirichlet process mixture model to realize grouping among single nucleotide polymorphism-associated phenotypic variations, providing biological plausibility. Through extensive simulations, we show the proposed method outperforms existing ones in heritability estimation and heritable traits selection under various scenarios. We finally apply the method to two large-scale imaging genetics datasets: the Alzheimer's Disease Neuroimaging Initiative and United Kingdom Biobank and show biologically meaningful results.

Impact of autism genetic risk on brain connectivity: a mechanism for the female protective effect.

The biological mechanisms underlying the greater prevalence of autism spectrum disorder in males than females remain poorly understood. One hypothesis posits that this female protective effect arises from genetic load for autism spectrum disorder differentially impacting male and female brains. To test this hypothesis, we investigated the impact of cumulative genetic risk for autism spectrum disorder on functional brain connectivity in a balanced sample of boys and girls with autism spectrum disorder and typically developing boys and girls (127 youth, ages 8-17). Brain connectivity analyses focused on the salience network, a core intrinsic functional connectivity network which has previously been implicated in autism spectrum disorder. The effects of polygenic risk on salience network functional connectivity were significantly modulated by participant sex, with genetic load for autism spectrum disorder influencing functional connectivity in boys with and without autism spectrum disorder but not girls. These findings support the hypothesis that autism spectrum disorder risk genes interact with sex differential processes, thereby contributing to the male bias in autism prevalence and proposing an underlying neurobiological mechanism for the female protective effect.

Identifying genetic risk variants associated with brain volumetric phenotypes via K-sample Ball Divergence method.

Regional human brain volumes including total area, average thickness, and total volume are heritable and associated with neurological disorders. However, the genetic architecture of brain structure and function is still largely unknown and worthy of exploring. The Pediatric Imaging, Neurocognition, and Genetics (PING) data set provides an excellent resource with genome-wide genetic data and related neuroimaging data. In this study, we perform genome-wide association studies (GWAS) of 315 brain volumetric phenotypes from the PING data set including 1036 samples with 539,865 single-nucleotide polymorphisms (SNPs). We introduce a nonparametric test based on K-sample Ball Divergence (KBD) to identify genetic risk variants that influence regional brain volumes. We carry out simulations to demonstrate that KBD is a powerful test for identifying significant SNPs associated with multivariate phenotypes although controlling the type I error rate. We successfully identify nine SNPs below a significance level of 5 × 10-5 for the PING data. Among the nine identified genetic variants, two SNPs rs486179 and rs562110 are located in the ADRA1A gene that is a well-known risk factor of mental illness, such as schizophrenia and attention deficit hyperactivity disorder. Our study suggests that the nonparametric test KBD is an effective method for identifying genetic variants associated with complex diseases in large-scale GWAS of multiple phenotypes.

Associations between brain imaging and polygenic scores of mental health and educational attainment in children aged 9-11.

Psychiatric disorders are highly heritable and polygenic, and many have their peak onset in late childhood and adolescence, a period of tremendous changes. Although the neurodevelopmental antecedents of mental illness are widely acknowledged, research in youth population cohorts is still scarce, preventing our progress towards the early characterization of these disorders. We included 7,124 children (9-11 years old) from the Adolescent Brain and Cognitive Development Study to map the associations of structural and diffusion brain imaging with common genetic variants and polygenic scores for psychiatric disorders and educational attainment. We used principal component analysis to derive imaging components, and calculated their heritability. We then assessed the relationship of imaging components with genetic and clinical psychiatric risk with univariate models and Canonical correlation analysis (CCA). Most imaging components had moderate heritability. Univariate models showed limited evidence and small associations of polygenic scores with brain structure at this age. CCA revealed two significant modes of covariation. The first mode linked higher polygenic scores for educational attainment with less externalizing problems and larger surface area. The second mode related higher polygenic scores for schizophrenia, bipolar disorder, and autism spectrum disorder to higher global cortical thickness, smaller white matter volumes of the fornix and cingulum, larger medial occipital surface area and smaller surface area of lateral and medial temporal regions. While cross-validation suggested limited generalizability, our results highlight the potential of multivariate models to better understand the transdiagnostic and distributed relationships between mental health and brain structure in late childhood.

The genetic architecture of language functional connectivity.

Language is a unique trait of the human species, of which the genetic architecture remains largely unknown. Through language disorders studies, many candidate genes were identified. However, such complex and multifactorial trait is unlikely to be driven by only few genes and case-control studies, suffering from a lack of power, struggle to uncover significant variants. In parallel, neuroimaging has significantly contributed to the understanding of structural and functional aspects of language in the human brain and the recent availability of large scale cohorts like UK Biobank have made possible to study language via image-derived endophenotypes in the general population. Because of its strong relationship with task-based fMRI (tbfMRI) activations and its easiness of acquisition, resting-state functional MRI (rsfMRI) have been more popularised, making it a good surrogate of functional neuronal processes. Taking advantage of such a synergistic system by aggregating effects across spatially distributed traits, we performed a multivariate genome-wide association study (mvGWAS) between genetic variations and resting-state functional connectivity (FC) of classical brain language areas in the inferior frontal (pars opercularis, triangularis and orbitalis), temporal and inferior parietal lobes (angular and supramarginal gyri), in 32,186 participants from UK Biobank. Twenty genomic loci were found associated with language FCs, out of which three were replicated in an independent replication sample. A locus in 3p11.1, regulating EPHA3 gene expression, is found associated with FCs of the semantic component of the language network, while a locus in 15q14, regulating THBS1 gene expression is found associated with FCs of the perceptual-motor language processing, bringing novel insights into the neurobiology of language.