Finemap-MiXeR: A variational Bayesian approach for genetic finemapping.

Genome-wide association studies (GWAS) implicate broad genomic loci containing clusters of highly correlated genetic variants. Finemapping techniques can select and prioritize variants within each GWAS locus which are more likely to have a functional influence on the trait. Here, we present a novel method, Finemap-MiXeR, for finemapping causal variants from GWAS summary statistics, controlling for correlation among variants due to linkage disequilibrium. Our method is based on a variational Bayesian approach and direct optimization of the Evidence Lower Bound (ELBO) of the likelihood function derived from the MiXeR model. After obtaining the analytical expression for ELBO's gradient, we apply Adaptive Moment Estimation (ADAM) algorithm for optimization, allowing us to obtain the posterior causal probability of each variant. Using these posterior causal probabilities, we validated Finemap-MiXeR across a wide range of scenarios using both synthetic data, and real data on height from the UK Biobank. Comparison of Finemap-MiXeR with two existing methods, FINEMAP and SuSiE RSS, demonstrated similar or improved accuracy. Furthermore, our method is computationally efficient in several aspects. For example, unlike many other methods in the literature, its computational complexity does not increase with the number of true causal variants in a locus and it does not require any matrix inversion operation. The mathematical framework of Finemap-MiXeR is flexible and may also be applied to other problems including cross-trait and cross-ancestry finemapping.

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.

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

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.

Genomic insights into the shared and distinct genetic architecture of cognitive function and schizophrenia.

Cognitive impairment is a major determinant of functional outcomes in schizophrenia, however, understanding of the biological mechanisms underpinning cognitive dysfunction in the disorder remains incomplete. Here, we apply Genomic Structural Equation Modelling to identify latent cognitive factors capturing genetic liabilities to 12 cognitive traits measured in the UK Biobank. We identified three broad factors that underly the genetic correlations between the cognitive tests. We explore the overlap between latent cognitive factors, schizophrenia, and schizophrenia symptom dimensions using a complementary set of statistical approaches, applied to data from the latest schizophrenia genome-wide association study (Ncase = 53,386, Ncontrol = 77,258) and the Thematically Organised Psychosis study (Ncase = 306, Ncontrol = 1060). Global genetic correlations showed a significant moderate negative genetic correlation between each cognitive factor and schizophrenia. Local genetic correlations implicated unique genomic regions underlying the overlap between schizophrenia and each cognitive factor. We found substantial polygenic overlap between each cognitive factor and schizophrenia and biological annotation of the shared loci implicated gene-sets related to neurodevelopment and neuronal function. Lastly, we show that the common genetic determinants of the latent cognitive factors are not predictive of schizophrenia symptoms in the Norwegian Thematically Organized Psychosis cohort. Overall, these findings inform our understanding of cognitive function in schizophrenia by demonstrating important differences in the shared genetic architecture of schizophrenia and cognitive abilities.

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.

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.

COSGAP: COntainerized Statistical Genetics Analysis Pipelines.

Summary: The collection and analysis of sensitive data in large-scale consortia for statistical genetics is hampered by multiple challenges, due to their non-shareable nature. Time-consuming issues in installing software frequently arise due to different operating systems, software dependencies, and limited internet access. For federated analysis across sites, it can be challenging to resolve different problems, including format requirements, data wrangling, setting up analysis on high-performance computing (HPC) facilities, etc. Easier, more standardized, automated protocols and pipelines can be solutions to overcome these issues. We have developed one such solution for statistical genetic data analysis using software container technologies. This solution, named COSGAP: "COntainerized Statistical Genetics Analysis Pipelines," consists of already established software tools placed into Singularity containers, alongside corresponding code and instructions on how to perform statistical genetic analyses, such as genome-wide association studies, polygenic scoring, LD score regression, Gaussian Mixture Models, and gene-set analysis. Using provided helper scripts written in Python, users can obtain auto-generated scripts to conduct the desired analysis either on HPC facilities or on a personal computer. COSGAP is actively being applied by users from different countries and projects to conduct genetic data analyses without spending much effort on software installation, converting data formats, and other technical requirements. Availability and implementation: COSGAP is freely available on GitHub (https://github.com/comorment/containers) under the GPLv3 license.

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.

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.

Distinct biological signature and modifiable risk factors underlie the comorbidity between major depressive disorder and cardiovascular disease.

Major depressive disorder (MDD) and cardiovascular disease (CVD) are often comorbid, resulting in excess morbidity and mortality. Here we show that CVDs share most of their genetic risk factors with MDD. Multivariate genome-wide association analysis of shared genetic liability between MDD and atherosclerotic CVD revealed seven loci and distinct patterns of tissue and brain cell-type enrichments, suggesting the involvement of the thalamus. Part of the genetic overlap was explained by shared inflammatory, metabolic and psychosocial or lifestyle risk factors. Our data indicated causal effects of genetic liability to MDD on CVD risk, but not from most CVDs to MDD, and showed that the causal effects were partly explained by metabolic and psychosocial or lifestyle factors. The distinct signature of MDD-atherosclerotic CVD comorbidity suggests an immunometabolic subtype of MDD that is more strongly associated with CVD than overall MDD. In summary, we identified biological mechanisms underlying MDD-CVD comorbidity and modifiable risk factors for prevention of CVD in individuals with MDD.

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

BACKGROUND: Carriers of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants exhibit regional and global brain differences compared with noncarriers. However, interpreting regional differences is challenging if a global difference drives the regional brain differences. Intraindividual variability measures can be used to test for regional differences beyond global differences in brain structure. METHODS: Magnetic resonance imaging data were used to obtain regional brain values for 1q21.1 distal deletion (n = 30) and duplication (n = 27) and 15q11.2 BP1-BP2 deletion (n = 170) and duplication (n = 243) carriers and matched noncarriers (n = 2350). Regional intra-deviation scores, i.e., the standardized difference between an individual's regional difference and global difference, were used to test for regional differences that diverge from the global difference. RESULTS: For the 1q21.1 distal deletion carriers, cortical surface area for regions in the medial visual cortex, posterior cingulate, and temporal pole differed less and regions in the prefrontal and superior temporal cortex differed more than the global difference in cortical surface area. For the 15q11.2 BP1-BP2 deletion carriers, cortical thickness in regions in the medial visual cortex, auditory cortex, and temporal pole differed less and the prefrontal and somatosensory cortex differed more than the global difference in cortical thickness. CONCLUSIONS: We find evidence for regional effects beyond differences in global brain measures in 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants. The results provide new insight into brain profiling of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants, with the potential to increase understanding of the mechanisms involved in altered neurodevelopment.

Identification of Novel Loci and Cross-Disorder Pleiotropy Through Multi-Ancestry Genome-Wide Analysis of Alcohol Use Disorder in Over One Million Individuals.

Alcohol use disorder (AUD) is highly heritable and burdensome worldwide. Genome-wide association studies (GWASs) can provide new evidence regarding the aetiology of AUD. We report a multi-ancestry GWASs across diverse ancestries focusing on a narrow AUD phenotype, using novel statistical tools in a total sample of 1,041,450 individuals [102,079 cases; European, 75,583; African, 20,689 (mostly African-American); Hispanic American, 3,449; East Asian, 2,254; South Asian, 104; descent]. Cross-ancestry functional analyses were performed with European and African samples. Thirty-seven genome-wide significant loci were identified, of which seven were novel for AUD and six for other alcohol phenotypes. Loci were mapped to genes enriched for brain regions relevant for AUD (striatum, hypothalamus, and prefrontal cortex) and potential drug targets (GABAergic, dopaminergic and serotonergic neurons). African-specific analysis yielded a unique pattern of immune-related gene sets. Polygenic overlap and positive genetic correlations showed extensive shared genetic architecture between AUD and both mental and general medical phenotypes, suggesting they are not only complications of alcohol use but also share genetic liability with AUD. Leveraging a cross-ancestry approach allowed identification of novel genetic loci for AUD and underscores the value of multi-ancestry genetic studies. These findings advance our understanding of AUD risk and clinically-relevant comorbidities.