Improved functional mapping of complex trait heritability with GSA-MiXeR implicates biologically specific gene sets.

While genome-wide association studies are increasingly successful in discovering genomic loci associated with complex human traits and disorders, the biological interpretation of these findings remains challenging. Here we developed the GSA-MiXeR analytical tool for gene set analysis (GSA), which fits a model for the heritability of individual genes, accounting for linkage disequilibrium across variants and allowing the quantification of partitioned heritability and fold enrichment for small gene sets. We validated the method using extensive simulations and sensitivity analyses. When applied to a diverse selection of complex traits and disorders, including schizophrenia, GSA-MiXeR prioritizes gene sets with greater biological specificity compared to standard GSA approaches, implicating voltage-gated calcium channel function and dopaminergic signaling for schizophrenia. Such biologically relevant gene sets, often with fewer than ten genes, are more likely to provide insights into the pathobiology of complex diseases and highlight potential drug targets.

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.

Dissecting the Shared Genetic Architecture of Common Epilepsies With Cortical Brain Morphology.

Background and Objectives: Epilepsies are associated with differences in cortical thickness (TH) and surface area (SA). However, the mechanisms underlying these relationships remain elusive. We investigated the extent to which these phenotypes share genetic influences. Methods: We analyzed genome-wide association study data on common epilepsies (n = 69,995) and TH and SA (n = 32,877) using Gaussian mixture modeling MiXeR and conjunctional false discovery rate (conjFDR) analysis to quantify their shared genetic architecture and identify overlapping loci. We biologically interrogated the loci using a variety of resources and validated in independent samples. Results: The epilepsies (2.4 k-2.9 k variants) were more polygenic than both SA (1.8 k variants) and TH (1.3 k variants). Despite absent genome-wide genetic correlations, there was a substantial genetic overlap between SA and genetic generalized epilepsy (GGE) (1.1 k), all epilepsies (1.1 k), and juvenile myoclonic epilepsy (JME) (0.7 k), as well as between TH and GGE (0.8 k), all epilepsies (0.7 k), and JME (0.8 k), estimated with MiXeR. Furthermore, conjFDR analysis identified 15 GGE loci jointly associated with SA and 15 with TH, 3 loci shared between SA and childhood absence epilepsy, and 6 loci overlapping between SA and JME. 23 loci were novel for epilepsies and 11 for cortical morphology. We observed a high degree of sign concordance in the independent samples. Discussion: Our findings show extensive genetic overlap between generalized epilepsies and cortical morphology, indicating a complex genetic relationship with mixed-effect directions. The results suggest that shared genetic influences may contribute to cortical abnormalities in epilepsies.

Unraveling the shared genetics of common epilepsies and general cognitive ability.

Objective: Cognitive impairment is prevalent among individuals with epilepsy, and it is possible that genetic factors can underlie this relationship. Here, we investigated the potential shared genetic basis of common epilepsies and general cognitive ability (COG). Methods: We applied linkage disequilibrium score (LDSC) regression, MiXeR and conjunctional false discovery rate (conjFDR) to analyze different aspects of genetic overlap between COG and epilepsies. We used the largest available genome-wide association study data on COG (n = 269,867) and common epilepsies (n = 27,559 cases, 42,436 controls), including the broad phenotypes 'all epilepsy', focal epilepsies and genetic generalized epilepsies (GGE), and as well as specific subtypes. We functionally annotated the identified loci using a variety of biological resources and validated the results in independent samples. Results: Using MiXeR, COG (11.2k variants) was estimated to be almost four times more polygenic than 'all epilepsy', GGE, juvenile myoclonic epilepsy (JME), and childhood absence epilepsy (CAE) (2.5k - 2.9k variants). The other epilepsy phenotypes were insufficiently powered for analysis. We show extensive genetic overlap between COG and epilepsies with significant negative genetic correlations (-0.23 to -0.04). COG was estimated to share 2.9k variants with both GGE and 'all epilepsy', and 2.3k variants with both JME and CAE. Using conjFDR, we identified 66 distinct loci shared between COG and epilepsies, including novel associations for GGE (27), 'all epilepsy' (5), JME (5) and CAE (5). The implicated genes were significantly expressed in multiple brain regions. The results were validated in independent samples (COG: p = 1.0 × 10-14; 'all epilepsy': p = 5.6 × 10-3). Significance: Our study demonstrates a substantial genetic basis shared between epilepsies and COG and identifies novel overlapping genomic loci. Enhancing our understanding of the relationship between epilepsies and COG may lead to the development of novel comorbidity-targeted epilepsy treatments.

Genome-wide association analyses identify 95 risk loci and provide insights into the neurobiology of post-traumatic stress disorder.

Post-traumatic stress disorder (PTSD) genetics are characterized by lower discoverability than most other psychiatric disorders. The contribution to biological understanding from previous genetic studies has thus been limited. We performed a multi-ancestry meta-analysis of genome-wide association studies across 1,222,882 individuals of European ancestry (137,136 cases) and 58,051 admixed individuals with African and Native American ancestry (13,624 cases). We identified 95 genome-wide significant loci (80 new). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (for example, GRIA1, GRM8 and CACNA1E), developmental, axon guidance and transcription factors (for example, FOXP2, EFNA5 and DCC), synaptic structure and function genes (for example, PCLO, NCAM1 and PDE4B) and endocrine or immune regulators (for example, ESR1, TRAF3 and TANK). Additional top genes influence stress, immune, fear and threat-related processes, previously hypothesized to underlie PTSD neurobiology. These findings strengthen our understanding of neurobiological systems relevant to PTSD pathophysiology, while also opening new areas for investigation.

Associations of plasma neurofilament light chain with cognition and neuroimaging measures in community-dwelling early old age men.

BACKGROUND: Plasma neurofilament light chain (NfL) is a promising biomarker of neurodegeneration with potential clinical utility in monitoring the progression of neurodegenerative diseases. However, the cross-sectional associations of plasma NfL with measures of cognition and brain have been inconsistent in community-dwelling populations. METHODS: We examined these associations in a large community-dwelling sample of early old age men (N = 969, mean age = 67.57 years, range = 61-73 years), who are either cognitively unimpaired (CU) or with mild cognitive impairment (MCI). Specifically, we investigated five cognitive domains (executive function, episodic memory, verbal fluency, processing speed, visual-spatial ability), as well as neuroimaging measures of gray and white matter. RESULTS: After adjusting for age, health status, and young adult general cognitive ability, plasma NfL level was only significantly associated with processing speed and white matter hyperintensity (WMH) volume, but not with other cognitive or neuroimaging measures. The association with processing speed was driven by individuals with MCI, as it was not detected in CU individuals. CONCLUSIONS: These results suggest that in early old age men without dementia, plasma NfL does not appear to be sensitive to cross-sectional individual differences in most domains of cognition or neuroimaging measures of gray and white matter. The revealed plasma NfL associations were limited to WMH for all participants and processing speed only within the MCI cohort. Importantly, considering cognitive status in community-based samples will better inform the interpretation of the relationships of plasma NfL with cognition and brain and may help resolve mixed findings in the literature.

Dissecting the genetic overlap between three complex phenotypes with trivariate MiXeR.

Comorbidities are an increasing global health challenge. Accumulating evidence suggests overlapping genetic architectures underlying comorbid complex human traits and disorders. The bivariate causal mixture model (MiXeR) can quantify the polygenic overlap between complex phenotypes beyond global genetic correlation. Still, the pattern of genetic overlap between three distinct phenotypes, which is important to better characterize multimorbidities, has previously not been possible to quantify. Here, we present and validate the trivariate MiXeR tool, which disentangles the pattern of genetic overlap between three phenotypes using summary statistics from genome-wide association studies (GWAS). Our simulations show that the trivariate MiXeR can reliably reconstruct different patterns of genetic overlap. We further demonstrate how the tool can be used to estimate the proportions of genetic overlap between three phenotypes using real GWAS data, providing examples of complex patterns of genetic overlap between diverse human traits and diseases that could not be deduced from bivariate analyses. This contributes to a better understanding of the etiology of complex phenotypes and the nature of their relationship, which may aid in dissecting comorbidity patterns and their biological underpinnings.

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

Quantitative MRI biomarker for classification of clinically significant prostate cancer: calibration for reproducibility across echo times.

Background: Restriction Spectrum Imaging restriction score (RSIrs) is a quantitative biomarker for detecting clinically significant prostate cancer (csPCa). However, the quantitative value of the RSIrs is affected by imaging parameters such as echo time (TE). Purpose: The purpose of the present study is to develop a calibration method to account for differences in echo times and facilitate use of RSIrs as a quantitative biomarker for the detection of csPCa. Methods: This study included 197 consecutive patients who underwent MRI and biopsy examination; 97 were diagnosed with csPCa (grade group ≥ 2). RSI data were acquired three times during the same session: twice at minimum TE∼75ms and once at TE=90ms (TEmin 1 , TEmin 2 , and TE90, respectively). A proposed calibration method, trained on patients without csPCa, estimated a linear scaling factor (f) for each of the four diffusion compartments (C) of the RSI signal model. A linear regression model was determined to match C-maps of TE90 to the reference C-maps of TEmin 1 within the interval ranging from 95 th to 99 th percentile of signal intensity within the prostate. RSIrs comparisons were made at 98 th percentile within each patient's prostate. We compared RSIrs from calibrated TE90 (RSIrs TE90corr ) and uncorrected TE90 (RSIrs TE90 ) to RSIrs from reference TEmin 1 (RSIrs TEmin1 ) and repeated TEmin 2 (RSIrs TEmin2 ). Calibration performance was evaluated with sensitivity, specificity, area under the ROC curve, positive predicted value, negative predicted value, and F1-score. Results: Scaling factors for C 1 , C 2 , C 3 , and C 4 were estimated as 1.70, 1.38, 1.03, and 1.19, respectively. In non-csPCa cases, the 98 th percentile of RSIrs TEmin2 and RSIrs TEmin1 differed by 0.27±0.86SI (mean±standard deviation), whereas RSIrs TE90 differed from RSIrs TEmin1 by 1.81±1.20SI. After calibration, this bias was reduced to -0.41±1.20SI, representing a 78% reduction in absolute error. For patients with csPCa, the difference was 0.54±1.98SI between RSIrs TEmin2 and RSIrs TEmin1 and 2.28±2.06SI between RSIrs TE90 and RSIrs TEmin1 . After calibration, the mean difference decreased to -0.86SI, a 38% reduction in absolute error. At the Youden index for patient-level classification of csPCa (8.94SI), RSIrs TEmin1 has a sensitivity of 66% and a specificity of 72%. Prior to calibration, RSIrs TE90 at the same threshold tended to over-diagnose benign cases (sensitivity 44%, specificity 88%). Post-calibration, RSIrs TE90corr performs more similarly to the reference (sensitivity 71%, specificity 62%). Conclusion: The proposed linear calibration method produces similar quantitative biomarker values for acquisitions with different TE, reducing TE-induced error by 78% and 38% for non-csPCa and csPCa, respectively.

Genome-wide analyses reveal shared genetic architecture and novel risk loci between opioid use disorder and general cognitive ability.

BACKGROUND: Opioid use disorder (OUD), a serious health burden worldwide, is associated with lower cognitive function. Recent studies have demonstrated a negative genetic correlation between OUD and general cognitive ability (COG), indicating a shared genetic basis. However, the specific genetic variants involved, and the underlying molecular mechanisms remain poorly understood. Here, we aimed to quantify and identify the genetic basis underlying OUD and COG. METHODS: We quantified the extent of genetic overlap between OUD and COG using a bivariate causal mixture model (MiXeR) and identified specific genetic loci applying conditional/conjunctional FDR. Finally, we investigated biological function and expression of implicated genes using available resources. RESULTS: We estimated that ~94% of OUD variants (4.8k out of 5.1k variants) also influence COG. We identified three novel OUD risk loci and one locus shared between OUD and COG. Loci identified implicated biological substrates in the basal ganglia. CONCLUSION: We provide new insights into the complex genetic risk architecture of OUD and its genetic relationship with COG.

Uncovering population contributions to the extracellular potential in the mouse visual system using Laminar Population Analysis.

The local field potential (LFP), the low-frequency part of the extracellular potential, reflects transmembrane currents in the vicinity of the recording electrode. Thought mainly to stem from currents caused by synaptic input, it provides information about neural activity complementary to that of spikes, the output of neurons. However, the many neural sources contributing to the LFP, and likewise the derived current source density (CSD), can often make it challenging to interpret. Efforts to improve its interpretability have included the application of statistical decomposition tools like principal component analysis (PCA) and independent component analysis (ICA) to disentangle the contributions from different neural sources. However, their underlying assumptions of, respectively, orthogonality and statistical independence are not always valid for the various processes or pathways generating LFP. Here, we expand upon and validate a decomposition algorithm named Laminar Population Analysis (LPA), which is based on physiological rather than statistical assumptions. LPA utilizes the multiunit activity (MUA) and LFP jointly to uncover the contributions of different populations to the LFP. To perform the validation of LPA, we used data simulated with the large-scale, biophysically detailed model of mouse V1 developed by the Allen Institute. We find that LPA can identify laminar positions within V1 and the temporal profiles of laminar population firing rates from the MUA. We also find that LPA can estimate the salient current sinks and sources generated by feedforward input from the lateral geniculate nucleus (LGN), recurrent activity in V1, and feedback input from the lateromedial (LM) area of visual cortex. LPA identifies and distinguishes these contributions with a greater accuracy than the alternative statistical decomposition methods, PCA and ICA. Lastly, we also demonstrate the application of LPA on experimentally recorded MUA and LFP from 24 animals in the publicly available Visual Coding dataset. Our results suggest that LPA can be used both as a method to estimate positions of laminar populations and to uncover salient features in LFP/CSD contributions from different populations.

How Real-World Data Can Facilitate the Development of Precision Medicine Treatment in Psychiatry.

Precision medicine has the ambition to improve treatment response and clinical outcomes through patient stratification and holds great potential for the treatment of mental disorders. However, several important factors are needed to transform current practice into a precision psychiatry framework. Most important are 1) the generation of accessible large real-world training and test data including genomic data integrated from multiple sources, 2) the development and validation of advanced analytical tools for stratification and prediction, and 3) the development of clinically useful management platforms for patient monitoring that can be integrated into health care systems in real-life settings. This narrative review summarizes strategies for obtaining the key elements-well-powered samples from large biobanks integrated with electronic health records and health registry data using novel artificial intelligence algorithms-to predict outcomes in severe mental disorders and translate these models into clinical management and treatment approaches. Key elements are massive mental health data and novel artificial intelligence algorithms. For the clinical translation of these strategies, we discuss a precision medicine platform for improved management of mental disorders. We use cases to illustrate how precision medicine interventions could be brought into psychiatry to improve the clinical outcomes of mental disorders.

Childhood Disadvantage Moderates Late Midlife Default Mode Network Cortical Microstructure and Visual Memory Association.

BACKGROUND: Childhood disadvantage is a prominent risk factor for cognitive and brain aging. Childhood disadvantage is associated with poorer episodic memory in late midlife and functional and structural brain abnormalities in the default mode network (DMN). Although age-related changes in DMN are associated with episodic memory declines in older adults, it remains unclear if childhood disadvantage has an enduring impact on this later-life brain-cognition relationship earlier in the aging process. Here, within the DMN, we examined whether its cortical microstructural integrity-an early marker of structural vulnerability that increases the risk for future cognitive decline and neurodegeneration-is associated with episodic memory in adults at ages 56-66, and whether childhood disadvantage moderates this association. METHODS: Cortical mean diffusivity (MD) obtained from diffusion magnetic resonance imaging was used to measure microstructural integrity in 350 community-dwelling men. We examined both visual and verbal episodic memory in relation to DMN MD and divided participants into disadvantaged and nondisadvantaged groups based on parental education and occupation. RESULTS: Higher DMN MD was associated with poorer visual memory but not verbal memory (ß = -0.11, p = .040 vs ß = -0.04, p = .535). This association was moderated by childhood disadvantage and was significant only in the disadvantaged group (ß = -0.26, p = .002 vs ß = -0.00, p = .957). CONCLUSIONS: Lower DMN cortical microstructural integrity may reflect visual memory vulnerability in cognitively normal adults earlier in the aging process. Individuals who experienced childhood disadvantage manifested greater vulnerability to cortical microstructure-related visual memory dysfunction than their nondisadvantaged counterparts who exhibited resilience in the face of low cortical microstructural integrity.

Identification of Novel Genomic Loci for Anxiety and Extensive Genetic Overlap with Psychiatric Disorders.

Background: Anxiety disorders are prevalent and anxiety symptoms co-occur with many psychiatric disorders. We aimed to identify genomic risk loci associated with anxiety, characterize its genetic architecture, and genetic overlap with psychiatric disorders. Methods: We used the GWAS of anxiety symptoms, schizophrenia, bipolar disorder, major depression, and attention deficit hyperactivity disorder (ADHD). We employed MiXeR and LAVA to characterize the genetic architecture and genetic overlap between the phenotypes. Conditional and conjunctional false discovery rate analyses were performed to boost the identification of genomic loci associated with anxiety and those shared with psychiatric disorders. Gene annotation and gene set analyses were conducted using OpenTargets and FUMA, respectively. Results: Anxiety was polygenic with 12.9k estimated genetic risk variants and overlapped extensively with psychiatric disorders (4.1-11.4k variants). MiXeR and LAVA revealed predominantly positive genetic correlations between anxiety and psychiatric disorders. We identified 114 novel loci for anxiety by conditioning on the psychiatric disorders. We also identified loci shared between anxiety and major depression (n = 47), bipolar disorder (n = 33), schizophrenia (n = 71), and ADHD (n = 20). Genes annotated to anxiety loci exhibit enrichment for a broader range of biological pathways and differential tissue expression in more diverse tissues than those annotated to the shared loci. Conclusions: Anxiety is a highly polygenic phenotype with extensive genetic overlap with psychiatric disorders. These genetic overlaps enabled the identification of novel loci for anxiety. The shared genetic architecture may underlie the extensive cross-disorder comorbidity of anxiety, and the identified genetic loci implicate molecular pathways that may lead to potential drug targets.

Restriction spectrum imaging with elastic image registration for automated evaluation of response to neoadjuvant therapy in breast cancer.

Purpose: Dynamic contrast-enhanced MRI (DCE) and apparent diffusion coefficient (ADC) are currently used to evaluate treatment response of breast cancer. The purpose of the current study was to evaluate the three-component Restriction Spectrum Imaging model (RSI3C), a recent diffusion-weighted MRI (DWI)-based tumor classification method, combined with elastic image registration, to automatically monitor breast tumor size throughout neoadjuvant therapy. Experimental design: Breast cancer patients (n=27) underwent multi-parametric 3T MRI at four time points during treatment. Elastically-registered DWI images were used to generate an automatic RSI3C response classifier, assessed against manual DCE tumor size measurements and mean ADC values. Predictions of therapy response during treatment and residual tumor post-treatment were assessed using non-pathological complete response (non-pCR) as an endpoint. Results: Ten patients experienced pCR. Prediction of non-pCR using ROC AUC (95% CI) for change in measured tumor size from pre-treatment time point to early-treatment time point was 0.65 (0.38-0.92) for the RSI3C classifier, 0.64 (0.36-0.91) for DCE, and 0.45 (0.16-0.75) for change in mean ADC. Sensitivity for detection of residual disease post-treatment was 0.71 (0.44-0.90) for the RSI3C classifier, compared to 0.88 (0.64-0.99) for DCE and 0.76 (0.50-0.93) for ADC. Specificity was 0.90 (0.56-1.00) for the RSI3C classifier, 0.70 (0.35-0.93) for DCE, and 0.50 (0.19-0.81) for ADC. Conclusion: The automatic RSI3C classifier with elastic image registration suggested prediction of response to treatment after only three weeks, and showed performance comparable to DCE for assessment of residual tumor post-therapy. RSI3C may guide clinical decision-making and enable tailored treatment regimens and cost-efficient evaluation of neoadjuvant therapy of breast cancer.

Cross-trait genome-wide association analysis of C-reactive protein level and psychiatric disorders.

C-reactive protein (CRP) tends to be elevated in individuals with psychiatric disorders. Recent findings have suggested a protective effect of the genetic liability to elevated CRP on schizophrenia risk and a causative effect on depression despite weak genetic correlations, while causal relationships with bipolar disorder were inconclusive. We investigated the shared genetic underpinnings of psychiatric disorders and variation in CRP levels. Genome-wide association studies for CRP (n = 575,531), bipolar disorder (n = 413,466), depression (n = 480,359), and schizophrenia (n = 130,644) were used in causal mixture models to compare CRP with psychiatric disorders based on polygenicity, discoverability, and genome-wide genetic overlap. The conjunctional false discovery rate method was used to identify specific shared genetic loci. Shared variants were mapped to putative causal genes, which were tested for overrepresentation among gene ontology gene-sets. CRP was six to ten times less polygenic (n = 1400 vs 8600-14,500 variants) and had a discoverability one to two orders of magnitude higher than psychiatric disorders. Most CRP-associated variants were overlapping with psychiatric disorders. We identified 401 genetic loci jointly associated with CRP and psychiatric disorders with mixed effect directions. Gene-set enrichment analyses identified predominantly CNS-related gene sets for CRP and each of depression and schizophrenia, and basic cellular processes for CRP and bipolar disorder. In conclusion, CRP has a markedly different genetic architecture to psychiatric disorders, but the majority of CRP associated variants are also implicated in psychiatric disorders. Shared genetic loci implicated CNS-related processes to a greater extent than immune processes, which may have implications for how we conceptualise causal relationships between CRP and psychiatric disorders.

Genome-wide analysis of anorexia nervosa and major psychiatric disorders and related traits reveals genetic overlap and identifies novel risk loci for anorexia nervosa.

Anorexia nervosa (AN) is a heritable eating disorder (50-60%) with an array of commonly comorbid psychiatric disorders and related traits. Although significant genetic correlations between AN and psychiatric disorders and related traits have been reported, their shared genetic architecture is largely understudied. We investigated the shared genetic architecture of AN and schizophrenia (SCZ), bipolar disorder (BIP), major depression (MD), mood instability (Mood), neuroticism (NEUR), and intelligence (INT). We applied the conditional false discovery rate (FDR) method to identify novel risk loci for AN, and conjunctional FDR to identify loci shared between AN and related phenotypes, to summarize statistics from relevant genome-wide association studies (GWAS). Individual GWAS samples varied from 72,517 to 420,879 participants. Using conditional FDR we identified 58 novel AN loci. Furthermore, we identified 38 unique loci shared between AN and major psychiatric disorders (SCZ, BIP, and MD) and 45 between AN and psychological traits (Mood, NEUR, and INT). In line with genetic correlations, the majority of shared loci showed concordant effect directions. Functional analyses revealed that the shared loci are involved in 65 unique pathways, several of which overlapped across analyses, including the "signal by MST1" pathway involved in Hippo signaling. In conclusion, we demonstrated genetic overlap between AN and major psychiatric disorders and related traits, and identified novel risk loci for AN by leveraging this overlap. Our results indicate that some shared characteristics between AN and related disorders and traits may have genetic underpinnings.

Psychiatric disorders and brain white matter exhibit genetic overlap implicating developmental and neural cell biology.

Improved understanding of the shared genetic architecture between psychiatric disorders and brain white matter may provide mechanistic insights for observed phenotypic associations. Our objective is to characterize the shared genetic architecture of bipolar disorder (BD), major depression (MD), and schizophrenia (SZ) with white matter fractional anisotropy (FA) and identify shared genetic loci to uncover biological underpinnings. We used genome-wide association study (GWAS) summary statistics for BD (n = 413,466), MD (n = 420,359), SZ (n = 320,404), and white matter FA (n = 33,292) to uncover the genetic architecture (i.e., polygenicity and discoverability) of each phenotype and their genetic overlap (i.e., genetic correlations, overlapping trait-influencing variants, and shared loci). This revealed that BD, MD, and SZ are at least 7-times more polygenic and less genetically discoverable than average FA. Even in the presence of weak genetic correlations (range = -0.05 to -0.09), average FA shared an estimated 42.5%, 43.0%, and 90.7% of trait-influencing variants as well as 12, 4, and 28 shared loci with BD, MD, and SZ, respectively. Shared variants were mapped to genes and tested for enrichment among gene-sets which implicated neurodevelopmental expression, neural cell types, myelin, and cell adhesion molecules. For BD and SZ, case vs control tract-level differences in FA associated with genetic correlations between those same tracts and the respective disorder (rBD = 0.83, p = 4.99e-7 and rSZ = 0.65, p = 5.79e-4). Genetic overlap at the tract-level was consistent with average FA results. Overall, these findings suggest a genetic basis for the involvement of brain white matter aberrations in the pathophysiology of psychiatric disorders.

Bidirectional genetic overlap between autism spectrum disorder and cognitive traits.

Autism spectrum disorder (ASD) is a highly heritable condition with a large variation in cognitive function. Here we investigated the shared genetic architecture between cognitive traits (intelligence (INT) and educational attainment (EDU)), and risk loci jointly associated with ASD and the cognitive traits. We analyzed data from genome-wide association studies (GWAS) of INT (n = 269,867), EDU (n = 766,345) and ASD (cases n = 18,381, controls n = 27,969). We used the bivariate causal mixture model (MiXeR) to estimate the total number of shared genetic variants, local analysis of co-variant annotation (LAVA) to estimate local genetic correlations, conditional false discovery rate (cond/conjFDR) to identify specific overlapping loci. The MiXeR analyses showed that 12.7k genetic variants are associated with ASD, of which 12.0k variants are shared with EDU, and 11.1k are shared with INT with both positive and negative relationships within overlapping variants. The majority (59-68%) of estimated shared loci have concordant effect directions, with a positive, albeit modest, genetic correlation between ASD and EDU (rg = 0.21, p = 2e-13) and INT (rg = 0.22, p = 4e-12). We discovered 43 loci jointly associated with ASD and cognitive traits (conjFDR<0.05), of which 27 were novel for ASD. Functional analysis revealed significant differential expression of candidate genes in the cerebellum and frontal cortex. To conclude, we quantified the genetic architecture shared between ASD and cognitive traits, demonstrated mixed effect directions, and identified the associated genetic loci and molecular pathways. The findings suggest that common genetic risk factors for ASD can underlie both better and worse cognitive functioning across the ASD spectrum, with different underlying biology.