An integrated framework for local genetic correlation analysis

The Genetic Architecture of the Human Corpus Callosum and its Subregions

The corpus callosum (CC) is the largest set of white matter fibers connecting the two hemispheres of the brain. In humans, it is essential for coordinating sensorimotor responses, performing associative/executive functions, and representing information in multiple dimensions. Understanding which genetic variants underpin corpus callosum morphometry, and their shared influence on cortical structure and susceptibility to neuropsychiatric disorders, can provide molecular insights into the CC's role in mediating cortical development and its contribution to neuropsychiatric disease. To characterize the morphometry of the midsagittal corpus callosum, we developed a publicly available artificial intelligence based tool to extract, parcellate, and calculate its total and regional area and thickness. Using the UK Biobank (UKB) and the Adolescent Brain Cognitive Development study (ABCD), we extracted measures of midsagittal corpus callosum morphometry and performed a genome-wide association study (GWAS) meta-analysis of European participants (combined N = 46,685). We then examined evidence for generalization to the non-European participants of the UKB and ABCD cohorts (combined N = 7,040). Post-GWAS analyses implicate prenatal intracellular organization and cell growth patterns, and high heritability in regions of open chromatin, suggesting transcriptional activity regulation in early development. Results suggest programmed cell death mediated by the immune system drives the thinning of the posterior body and isthmus. Global and local genetic overlap, along with causal genetic liability, between the corpus callosum, cerebral cortex, and neuropsychiatric disorders such as attention-deficit/hyperactivity and bipolar disorders were identified. These results provide insight into variability of corpus callosum development, its genetic influence on the cerebral cortex, and biological mechanisms related to neuropsychiatric dysfunction.

Shared genetic etiology of vessel diseases: A genome-wide multi-traits association analysis

BACKGROUND

The comorbidity among vascular diseases has been widely reported, however, the contribution of shared genetic components remains ambiguous.

METHODS

Based on genome-wide association study summary statistics, we employed statistical genetics methodologies to explore the shared genetic basis of eight vascular diseases: coronary artery disease, abdominal aortic aneurysm, ischemic stroke, peripheral artery disease, thoracic aortic aneurysm, phlebitis, varicose veins, and venous thromboembolism. We assessed global and local genetic correlations among these disorders by linkage disequilibrium score regression, high-definition likelihood, and local analysis of variant association. Cross-trait analyses conducted with CPASSOC identified pleiotropic variants and loci. Further, biological pathways at the multi-omics level were explored using multimarker analysis of genomic annotation, transcriptome-wide and proteome-wide association studies. Causal associations among the vascular diseases were evaluated by mendelian randomization and latent causal variable to assess vertical pleiotropic effects.

RESULTS

We found significant global genetic associations in 18 pairs of vascular diseases. Additionally, we discovered 317 unique genomic regions where at least one pair of traits demonstrated significant correlation. Multi-trait association analysis identified 19,361 significant potential pleiotropic variants in 274 independent pleiotropic loci. Multi-trait colocalization analysis revealed 56 colocalized loci in specific disease sets. Gene-based analysis identified 700 potential pleiotropic genes, which were subsequently validated at both transcriptome and protein levels. Gene-set enrichment analysis supports the role of biological pathways such as vessel wall structure, coagulation and lipid transport in vascular disease. Additionally, 7 pairs of vascular diseases have a causal relationship.

CONCLUSIONS

Our study indicates a shared genetic basis and the presence of common risk genes among vascular diseases. These findings offer novel insights into potential mechanisms underlying the association between vascular diseases, as well as provide guidance for interventions and treatments of multi-vascular conditions.

Using Genetics to Investigate Relationships between Phenotypes: Application to Endometrial Cancer

Genome-wide association studies (GWAS) have accelerated the exploration of genotype-phenotype associations, facilitating the discovery of replicable genetic markers associated with specific traits or complex diseases. This narrative review explores the statistical methodologies developed using GWAS data to investigate relationships between various phenotypes, focusing on endometrial cancer, the most prevalent gynecological malignancy in developed nations. Advancements in analytical techniques such as genetic correlation, colocalization, cross-trait locus identification, and causal inference analyses have enabled deeper exploration of associations between different phenotypes, enhancing statistical power to uncover novel genetic risk regions. These analyses have unveiled shared genetic associations between endometrial cancer and many phenotypes, enabling identification of novel endometrial cancer risk loci and furthering our understanding of risk factors and biological processes underlying this disease. The current status of research in endometrial cancer is robust; however, this review demonstrates that further opportunities exist in statistical genetics that hold promise for advancing the understanding of endometrial cancer and other complex diseases.

Trait selection strategy in multi-trait GWAS: Boosting SNP discoverability

Since the first genome-wide association studies (GWASs), thousands of variant-trait associations have been discovered. However, comprehensively mapping the genetic determinant of complex traits through univariate testing can require prohibitive sample sizes. Multi-trait GWAS can circumvent this issue and improve statistical power by leveraging the joint genetic architecture of human phenotypes. Although many methodological hurdles of multi-trait testing have been solved, the strategy to select traits has been overlooked. In this study, we conducted multi-trait GWAS on approximately 20,000 combinations of 72 traits using an omnibus test as implemented in the Joint Analysis of Summary Statistics. We assessed which genetic features of the sets of traits analyzed were associated with an increased detection of variants compared with univariate screening. Several features of the set of traits, including the heritability, the number of traits, and the genetic correlation, drive the multi-trait test gain. Using these features jointly in predictive models captures a large fraction of the power gain of the multi-trait test (Pearson's r between the observed and predicted gain equals 0.43, p < 1.6 × 10-60). Applying an alternative multi-trait approach (Multi-Trait Analysis of GWAS), we identified similar features of interest, but with an overall 70% lower number of new associations. Finally, selecting sets based on our data-driven models systematically outperformed the common strategy of selecting clinically similar traits. This work provides a unique picture of the determinant of multi-trait GWAS statistical power and outlines practical strategies for multi-trait testing.

A parametric bootstrap approach for computing confidence intervals for genetic correlations with application to genetically determined protein-protein networks

Genetic correlation refers to the correlation between genetic determinants of a pair of traits. When using individual-level data, it is typically estimated based on a bivariate model specification where the correlation between the two variables is identifiable and can be estimated from a covariance model that incorporates the genetic relationship between individuals, e.g., using a pre-specified kinship matrix. Inference relying on asymptotic normality of the genetic correlation parameter estimates may be inaccurate when the sample size is low, when the genetic correlation is close to the boundary of the parameter space, and when the heritability of at least one of the traits is low. We address this problem by developing a parametric bootstrap procedure to construct confidence intervals for genetic correlation estimates. The procedure simulates paired traits under a range of heritability and genetic correlation parameters, and it uses the population structure encapsulated by the kinship matrix. Heritabilities and genetic correlations are estimated using the close-form, method of moment, Haseman-Elston regression estimators. The proposed parametric bootstrap procedure is especially useful when genetic correlations are computed on pairs of thousands of traits measured on the same exact set of individuals. We demonstrate the parametric bootstrap approach on a proteomics dataset from the Jackson Heart Study.

Shared Genetic Architecture Among Gastrointestinal Diseases, Schizophrenia, and Brain Subcortical Volumes

BACKGROUND AND HYPOTHESIS

The gut-brain axis plays important roles in both gastrointestinal diseases (GI diseases) and schizophrenia (SCZ). Moreover, both GI diseases and SCZ exhibit notable abnormalities in brain subcortical volumes. However, the genetic mechanisms underlying the comorbidity of these diseases and the shared alterations in brain subcortical volumes remain unclear.

STUDY DESIGN

Using the genome-wide association studies data of SCZ, 14 brain subcortical volumes, and 8 GI diseases, the global polygenic overlap and local genetic correlations were identified, as well as the shared genetic variants among those phenotypes. Furthermore, we conducted multi-trait colocalization analyses to bolster our findings. Functional annotations, cell-type enrichment, and protein-protein interaction (PPI) analyses were carried out to reveal the critical etiology and pathology mechanisms.

STUDY RESULTS

The global polygenic overlap and local genetic correlations informed the close relationships between SCZ and both GI diseases and brain subcortical volumes. Moreover, 84 unique lead-shared variants were identified. The associated genes were linked to vital biological processes within the immune system. Additionally, significant correlations were observed with key immune cells and the PPI analysis identified several histone-associated hub genes. These findings highlighted the pivotal roles played by the immune system for both SCZ and GI diseases, along with the shared alterations in brain subcortical volumes.

CONCLUSIONS

These findings revealed the shared genetic architecture contributing to SCZ and GI diseases, as well as their shared alterations in brain subcortical volumes. These insights have substantial implications for the concurrent development of intervention and therapy targets for these diseases.

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.

Investigating the shared genetic architecture between attention-deficit/hyperactivity disorder and risk taking behavior: A large-scale genomewide cross-trait analysis

BACKGROUND

This study aims to explore the genetic architecture shared between Attention-Deficit/Hyperactivity Disorder (ADHD) and risk behavior.

METHODS

Based on the latest large-scale Genome-wide association studies (GWAS), we firstly employed Linkage disequilibrium score regression (LDSC) and Local Analysis of Variant Association (LAVA) to investigate the genetic correlation between risk behavior and ADHD. Then, we conducted cross-trait analysis to identified the Pleiotropic loci. Finally, bidirectional Mendelian randomization analysis (MR) was applied to examine the causal relationship.

RESULTS

We found a significant positive genetic correlation between ADHD and risk-taking behavior (rg = 0.351, p = 6.50E-37). The cross-trait meta-analysis identified 27 significant SNPs shared between ADHD and risk behavior. The most significant locus, located near the CADM2 gene on chromosome 3, had been identified associated with this two trait (pADHD = 3.07E-05 and prisk-taking behavior = 2.47E-30). The same situation can also be observed near the FOXP2 gene on chromosome 7 (rs8180817, pmeta = 5.72E-21). We found CCDC171 gene and other genes played a significant role in ADHD and risk behavior in mRNA level. Bidirectional MR analysis found a causal relationship between them.

LIMITATION

The majority of our data sources were of European origin, which may limit the generalizability of our findings to other ethnic populations.

CONCLUSION

This article reveals in depth the shared genetic structure between ADHD and risk-taking behavior, finding a significant positive genetic correlation between ADHD and risk-taking behavior. Providing insights for the future treatment and management of these two traits.

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.

Shared genetic aetiology of respiratory diseases: a genome-wide multitraits association analysis

OBJECTIVE

This study aims to explore the common genetic basis between respiratory diseases and to identify shared molecular and biological mechanisms.

METHODS

This genome-wide pleiotropic association study uses multiple statistical methods to systematically analyse the shared genetic basis between five respiratory diseases (asthma, chronic obstructive pulmonary disease, idiopathic pulmonary fibrosis, lung cancer and snoring) using the largest publicly available genome wide association studies summary statistics. The missions of this study are to evaluate global and local genetic correlations, to identify pleiotropic loci, to elucidate biological pathways at the multiomics level and to explore causal relationships between respiratory diseases. Data were collected from 27 November 2022 to 30 March 2023 and analysed from 14 April 2023 to 13 July 2023.

MAIN OUTCOMES AND MEASURES

The primary outcomes are shared genetic loci, pleiotropic genes, biological pathways and estimates of genetic correlations and causal effects.

RESULTS

Significant genetic correlations were found for 10 paired traits in 5 respiratory diseases. Cross-Phenotype Association identified 12 400 significant potential pleiotropic single-nucleotide polymorphism at 156 independent pleiotropic loci. In addition, multitrait colocalisation analysis identified 15 colocalised loci and a subset of colocalised traits. Gene-based analyses identified 432 potential pleiotropic genes and were further validated at the transcriptome and protein levels. Both pathway enrichment and single-cell enrichment analyses supported the role of the immune system in respiratory diseases. Additionally, five pairs of respiratory diseases have a causal relationship.

CONCLUSIONS AND RELEVANCE

This study reveals the common genetic basis and pleiotropic genes among respiratory diseases. It provides strong evidence for further therapeutic strategies and risk prediction for the phenomenon of respiratory disease comorbidity.

Shared Genetic Architectures between Coronary Artery Disease and Type 2 Diabetes Mellitus in East Asian and European Populations

Coronary artery disease (CAD) is a common comorbidity of type 2 diabetes mellitus (T2DM). However, the pathophysiology connecting these two phenotypes remains to be further understood. Combined analysis in multi-ethnic populations can help contribute to deepening our understanding of biological mechanisms caused by shared genetic loci. We applied genetic correlation analysis and then performed conditional and joint association analyses in Chinese, Japanese, and European populations to identify the genetic variants jointly associated with CAD and T2DM. Next, the associations between genes and the two traits were also explored. Finally, fine-mapping and functional enrichment analysis were employed to identify the potential causal variants and pathways. Genetic correlation results indicated significant genetic overlap between CAD and T2DM in the three populations. Over 10,000 shared signals were identified, and 587 were shared by East Asian and European populations. Fifty-six novel shared genes were found to have significant effects on both CAD and T2DM. Most loci were fine-mapped to plausible causal variant sets. Several similarities and differences of the involved genes in GO terms and KEGG pathways were revealed across East Asian and European populations. These findings highlight the importance of immunoregulation, neuroregulation, heart development, and the regulation of glucose metabolism in shared etiological mechanisms between CAD and T2DM.

Brain-wide pleiotropy investigation of alcohol drinking and tobacco smoking behaviors

To investigate the pleiotropic mechanisms linking brain structure and function to alcohol drinking and tobacco smoking, we integrated genome-wide data generated by the GWAS and Sequencing Consortium of Alcohol and Nicotine use (GSCAN; up to 805,431 participants) with information related to 3,935 brain imaging-derived phenotypes (IDPs) available from UK Biobank (N=33,224). We observed global genetic correlation of smoking behaviors with white matter hyperintensities, the morphology of the superior longitudinal fasciculus, and the mean thickness of pole-occipital. With respect to the latter brain IDP, we identified a local genetic correlation with age at which the individual began smoking regularly (hg38 chr2:35,895,678-36,640,246: rho=1, p=1.01×10 -5 ). This region has been previously associated with smoking initiation, educational attainment, chronotype, and cortical thickness. Our genetically informed causal inference analysis using both latent causal variable approach and Mendelian randomization linked the activity of prefrontal and premotor cortex and that of superior and inferior precentral sulci, and cingulate sulci to the number of alcoholic drinks per week (genetic causality proportion, gcp=0.38, p=8.9×10 -4 , rho=-0.18±0.07; inverse variance weighting, IVW beta=-0.04, 95%CI=-0.07 - -0.01). This relationship could be related to the role of these brain regions in the modulation of reward-seeking motivation and the processing of social cues. Overall, our brain-wide investigation highlighted that different pleiotropic mechanisms likely contribute to the relationship of brain structure and function with alcohol drinking and tobacco smoking, suggesting decision-making activities and chemosensory processing as modulators of propensity towards alcohol and tobacco consumption.

Genome-wide cross-disease analyses highlight causality and shared biological pathways of type 2 diabetes with gastrointestinal disorders

Studies suggest links between diabetes and gastrointestinal (GI) traits; however, their underlying biological mechanisms remain unclear. Here, we comprehensively assess the genetic relationship between type 2 diabetes (T2D) and GI disorders. Our study demonstrates a significant positive global genetic correlation of T2D with peptic ulcer disease (PUD), irritable bowel syndrome (IBS), gastritis-duodenitis, gastroesophageal reflux disease (GERD), and diverticular disease, but not inflammatory bowel disease (IBD). We identify several positive local genetic correlations (negative for T2D - IBD) contributing to T2D's relationship with GI disorders. Univariable and multivariable Mendelian randomisation analyses suggest causal effects of T2D on PUD and gastritis-duodenitis and bidirectionally with GERD. Gene-based analyses reveal a gene-level genetic overlap between T2D and GI disorders and identify several shared genes reaching genome-wide significance. Pathway-based study implicates leptin (T2D - IBD), thyroid, interferon, and notch signalling (T2D - IBS), abnormal circulating calcium (T2D - PUD), cardiovascular, viral, proinflammatory and (auto)immune-mediated mechanisms in T2D and GI disorders. These findings support a risk-increasing genetic overlap between T2D and GI disorders (except IBD), implicate shared biological pathways with putative causality for certain T2D - GI pairs, and identify targets for further investigation.

Systems biology dissection of PTSD and MDD across brain regions, cell types, and blood

The molecular pathology of stress-related disorders remains elusive. Our brain multiregion, multiomic study of posttraumatic stress disorder (PTSD) and major depressive disorder (MDD) included the central nucleus of the amygdala, hippocampal dentate gyrus, and medial prefrontal cortex (mPFC). Genes and exons within the mPFC carried most disease signals replicated across two independent cohorts. Pathways pointed to immune function, neuronal and synaptic regulation, and stress hormones. Multiomic factor and gene network analyses provided the underlying genomic structure. Single nucleus RNA sequencing in dorsolateral PFC revealed dysregulated (stress-related) signals in neuronal and non-neuronal cell types. Analyses of brain-blood intersections in >50,000 UK Biobank participants were conducted along with fine-mapping of the results of PTSD and MDD genome-wide association studies to distinguish risk from disease processes. Our data suggest shared and distinct molecular pathology in both disorders and propose potential therapeutic targets and biomarkers.

Improving multi-population genomic prediction accuracy using multi-trait GBLUP models which incorporate global or local genetic correlation information

In the application of genomic prediction, a situation often faced is that there are multiple populations in which genomic prediction (GP) need to be conducted. A common way to handle the multi-population GP is simply to combine the multiple populations into a single population. However, since these populations may be subject to different environments, there may exist genotype-environment interactions which may affect the accuracy of genomic prediction. In this study, we demonstrated that multi-trait genomic best linear unbiased prediction (MTGBLUP) can be used for multi-population genomic prediction, whereby the performances of a trait in different populations are regarded as different traits, and thus multi-population prediction is regarded as multi-trait prediction by employing the between-population genetic correlation. Using real datasets, we proved that MTGBLUP outperformed the conventional multi-population model that simply combines different populations together. We further proposed that MTGBLUP can be improved by partitioning the global between-population genetic correlation into local genetic correlations (LGC). We suggested two LGC models, LGC-model-1 and LGC-model-2, which partition the genome into regions with and without significant LGC (LGC-model-1) or regions with and without strong LGC (LGC-model-2). In analysis of real datasets, we demonstrated that the LGC models could increase universally the prediction accuracy and the relative improvement over MTGBLUP reached up to 163.86% (25.64% on average).

Genetic evidence for the causal effects of C-reactive protein on self-reported habitual sleep duration

Inflammatory responses to acute stimuli are proposed to regulate sleep, but the relationship between chronic inflammation and habitual sleep duration is elusive. Here, we study this relation using genetically predicted level of chronic inflammation, indexed by CRP and IL6 signaling, and self-reported sleep duration. By Mendelian randomization analysis, we show that elevated CRP level within <10 mg/L has a homeostatic effect that facilitates maintaining 7-8 h sleep duration per day - making short-sleepers sleep longer (p = 2.42 × 10-2) and long-sleepers sleep shorter (1.87 × 10-7); but it is not associated with the overall sleep duration (p = 0.17). This homeostatic effect replicated in an independent CRP dataset. We observed causal effects of the soluble interleukin 6 receptor and gp130 on overall sleep duration (p = 1.62 × 10-8, p = 2.61 × 10-58, respectively), but these effects disappeared when CRP effects were accounted for in the model. Using polygenic score analysis, we found that the homeostatic effect of CRP on sleep duration stems primarily from the genetic variants within the CRP gene region: when genetic variants outside of this region were used to predict CRP levels, the opposite direction of effect was observed. In conclusion, we show that elevated CRP level may causally facilitate maintaining an optimal sleep duration that is beneficial to health, thus updating our current knowledge of immune regulation on sleep.

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.

Unraveling COVID-19 relationship with anxiety disorders and symptoms using genome-wide data

BACKGROUND

There is still a limited understanding of the dynamics contributing to the comorbidity of COVID-19 and anxiety outcomes.

METHODS

To dissect the pleiotropic mechanisms contributing to COVID-19/anxiety comorbidity, we used genome-wide data from UK Biobank (up to 420,531 participants), FinnGen Project (up to 329,077 participants), Million Veteran Program (175,163 participants), and COVID-19 Host Genetics Initiative (up to 122,616 cases and 2,475,240 controls). Specifically, we assessed global and local genetic correlation and genetically inferred effects linking COVID-19 outcomes (infection, hospitalization, and severe respiratory symptoms) to anxiety disorders and symptoms.

RESULTS

We observed a strong genetic correlation of anxiety disorder with COVID-19 positive status (rg = 0.35, p = 2×10-4) and COVID-19 hospitalization (rg = 0.31, p = 7.2×10-4). Among anxiety symptoms, "Tense, sore, or aching muscles during worst period of anxiety" was genetically correlated with COVID-19 positive status (rg = 0.33, p = 0.001), while "Frequent trouble falling or staying asleep during worst period of anxiety" was genetically correlated with COVID-19 hospitalization (rg = 0.24, p = 0.004). Through a latent causal variable analysis, we observed that COVID-19 outcomes have statistically significant genetic causality proportion (gcp) on anxiety symptoms (e.g., COVID-19 positive status→"Recent easy annoyance or irritability" │gcp│ = 0.18, p = 6.72×10-17). Conversely, anxiety disorders appear to have a possible causal effect on COVID-19 (│gcp│ = 0.38, p = 3.17×10-9). Additionally, we also identified multiple loci with evidence of local genetic correlation between anxiety and COVID-19. These appear to be related to genetic effects shared with lung function, brain morphology, alcohol and tobacco use, and hematologic parameters.

CONCLUSIONS

This study provided insights into the pleiotropic mechanisms linking COVID-19 and anxiety outcomes, suggesting differences between dynamics related to anxiety disorders and those related to anxiety symptoms.

Insights into the genetic architecture of cerebellar lobules derived from the UK Biobank

In this work we endeavor to further understand the genetic architecture of the cerebellum by examining the genetic underpinnings of the different cerebellar lob(ul)es, identifying their genetic relation to cortical and subcortical regions, as well as to psychiatric disorders, as well as traces of their evolutionary trajectories. We confirm the moderate heritability of cerebellar volumes, and reveal genetic clustering and variability across their different substructures, which warranted a detailed analysis using this higher structural resolution. We replicated known genetic correlations with several subcortical volumes, and report new cortico-cerebellar genetic correlations, including negative genetic correlations between anterior cerebellar lobules and cingulate, and positive ones between lateral Crus I and lobule VI with cortical measures in the fusiform region. Heritability partitioning for evolutionary annotations highlighted that the vermis of Crus II has depleted heritability in genomic regions of "archaic introgression deserts", but no enrichment/depletion of heritability in any other cerebellar regions. Taken together, these findings reveal novel insights into the genetic underpinnings of the different cerebellar lobules.

Common risk alleles for schizophrenia within the major histocompatibility complex predict white matter microstructure

Recent research has highlighted the role of complement genes in shaping the microstructure of the brain during early development, and in contributing to common allele risk for Schizophrenia. We hypothesised that common risk variants for schizophrenia within complement genes will associate with structural changes in white matter microstructure within tracts innervating the frontal lobe. Results showed that risk alleles within the complement gene set, but also intergenic alleles, significantly predict axonal density in white matter tracts connecting frontal cortex with parietal, temporal and occipital cortices. Specifically, risk alleles within the Major Histocompatibility Complex region in chromosome 6 appeared to drive these associations. No significant associations were found for the orientation dispersion index. These results suggest that changes in axonal packing - but not in axonal coherence - determined by common risk alleles within the MHC genomic region - including variants related to the Complement system - appear as a potential neurobiological mechanism for schizophrenia.

Characterizing the polygenic overlap and shared loci between rheumatoid arthritis and cardiovascular diseases

BACKGROUND

Despite substantial research revealing that patients with rheumatoid arthritis (RA) have excessive morbidity and mortality of cardiovascular disease (CVD), the mechanism underlying this association has not been fully known. This study aims to systematically investigate the phenotypic and genetic correlation between RA and CVD.

METHODS

Based on UK Biobank, we conducted two cohort studies to evaluate the phenotypic relationships between RA and CVD, including atrial fibrillation (AF), coronary artery disease (CAD), heart failure (HF), and stroke. Next, we used linkage disequilibrium score regression, Local Analysis of [co]Variant Association, and bivariate causal mixture model (MiXeR) methods to examine the genetic correlation and polygenic overlap between RA and CVD, using genome-wide association summary statistics. Furthermore, we explored specific shared genetic loci by conjunctional false discovery rate analysis and association analysis based on subsets.

RESULTS

Compared with the general population, RA patients showed a higher incidence of CVD (hazard ratio [HR] = 1.21, 95% confidence interval [CI]: 1.15-1.28). We observed positive genetic correlations of RA with AF and stroke, and a mixture of negative and positive local genetic correlations underlying the global genetic correlation for CAD and HF, with 13 ~ 33% of shared genetic variants for these trait pairs. We further identified 23 pleiotropic loci associated with RA and at least one CVD, including one novel locus (rs7098414, TSPAN14, 10q23.1). Genes mapped to these shared loci were enriched in immune and inflammatory-related pathways, and modifiable risk factors, such as high diastolic blood pressure.

CONCLUSIONS

This study revealed the shared genetic architecture of RA and CVD, which may facilitate drug target identification and improved clinical management.

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.

Genetic influences and causal pathways shared between cannabis use disorder and other substance use traits

Cannabis use disorder (CanUD) has increased with the legalization of the use of cannabis. Around 20% of individuals using cannabis develop CanUD, and the number of users has grown with increasing ease of access. CanUD and other substance use disorders (SUDs) are associated phenotypically and genetically. We leveraged new CanUD genomics data to undertake genetically-informed analyses with unprecedented power, to investigate the genetic architecture and causal relationships between CanUD and lifetime cannabis use with risk for developing SUDs and substance use traits. Analyses included calculating local and global genetic correlations, genomic structural equation modeling (genomicSEM), and Mendelian Randomization (MR). Results from the genetic correlation and genomicSEM analyses demonstrated that CanUD and cannabis use differ in their relationships with SUDs and substance use traits. We found significant causal effects of CanUD influencing all the analyzed traits: opioid use disorder (OUD) (Inverse variant weighted, IVW β = 0.925 ± 0.082), problematic alcohol use (PAU) (IVW β = 0.443 ± 0.030), drinks per week (DPW) (IVW β = 0.182 ± 0.025), Fagerström Test for Nicotine Dependence (FTND) (IVW β = 0.183 ± 0.052), cigarettes per day (IVW β = 0.150 ± 0.045), current versus former smokers (IVW β = 0.178 ± 0.052), and smoking initiation (IVW β = 0.405 ± 0.042). We also found evidence of bidirectionality showing that OUD, PAU, smoking initiation, smoking cessation, and DPW all increase risk of developing CanUD. For cannabis use, bidirectional relationships were inferred with PAU, smoking initiation, and DPW; cannabis use was also associated with a higher risk of developing OUD (IVW β = 0.785 ± 0.266). GenomicSEM confirmed that CanUD and cannabis use load onto different genetic factors. We conclude that CanUD and cannabis use can increase the risk of developing other SUDs. This has substantial public health implications; the move towards legalization of cannabis use may be expected to increase other kinds of problematic substance use. These harmful outcomes are in addition to the medical harms associated directly with CanUD.

Genetic overlap between Alzheimer's disease and immune-mediated diseases: an atlas of shared genetic determinants and biological convergence

The occurrence of immune disease comorbidities in Alzheimer's disease (AD) has been observed in both epidemiological and molecular studies, suggesting a neuroinflammatory basis in AD. However, their shared genetic components have not been systematically studied. Here, we composed an atlas of the shared genetic associations between 11 immune-mediated diseases and AD by analyzing genome-wide association studies (GWAS) summary statistics. Our results unveiled a significant genetic overlap between AD and 11 individual immune-mediated diseases despite negligible genetic correlations, suggesting a complex shared genetic architecture distributed across the genome. The shared loci between AD and immune-mediated diseases implicated several genes, including GRAMD1B, FUT2, ADAMTS4, HBEGF, WNT3, TSPAN14, DHODH, ABCB9, and TNIP1, all of which are protein-coding genes and thus potential drug targets. Top biological pathways enriched with these identified shared genes were related to the immune system and cell adhesion. In addition, in silico single-cell analyses showed enrichment of immune and brain cells, including neurons and microglia. In summary, our results suggest a genetic relationship between AD and the 11 immune-mediated diseases, pinpointing the existence of a shared however non-causal genetic basis. These identified protein-coding genes have the potential to serve as a novel path to therapeutic interventions for both AD and immune-mediated diseases and their comorbidities.

Analyzing bivariate cross-trait genetic architecture in GWAS summary statistics with the BIGA cloud computing platform

As large-scale biobanks provide increasing access to deep phenotyping and genomic data, genome-wide association studies (GWAS) are rapidly uncovering the genetic architecture behind various complex traits and diseases. GWAS publications typically make their summary-level data (GWAS summary statistics) publicly available, enabling further exploration of genetic overlaps between phenotypes gathered from different studies and cohorts. However, systematically analyzing high-dimensional GWAS summary statistics for thousands of phenotypes can be both logistically challenging and computationally demanding. In this paper, we introduce BIGA (https://bigagwas.org/), a website that aims to offer unified data analysis pipelines and processed data resources for cross-trait genetic architecture analyses using GWAS summary statistics. We have developed a framework to implement statistical genetics tools on a cloud computing platform, combined with extensive curated GWAS data resources. Through BIGA, users can upload data, submit jobs, and share results, providing the research community with a convenient tool for consolidating GWAS data and generating new insights.

Genetic correlations and causal relationships between cardio-metabolic traits and sepsis

Cardio-metabolic traits have been reported to be associated with the development of sepsis. It is, however, unclear whether these co-morbidities reflect causal associations, shared genetic heritability, or are confounded by environmental factors. We performed three analyses to explore the relationships between cardio-metabolic traits and sepsis. Mendelian randomization (MR) study to evaluate the causal effects of multiple cardio-metabolic traits on sepsis. Global genetic correlation analysis to explore the correlations between cardio-metabolic traits and sepsis. Local genetic correlation (GC) analysis to explore shared genetic heritability between cardio-metabolic traits and sepsis. Some loci were further examined for related genes responsible for the causal relationships. Genetic associations were obtained from the UK Biobank data or published large-scale genome-wide association studies with sample sizes between 200,000 to 750,000. In MR, we found causality between BMI and sepsis (OR: 1.53 [1.4-1.67]; p < 0.001). Body mass index (BMI), which is confirmed by sensitivity analyses and multivariable MR adjusting for confounding factors. Global GC analysis showed a significant correlation between BMI and sepsis (rg = 0.55, p < 0.001). More cardio-metabolic traits were identified to be correlated to the sepsis onset such as CRP (rg = 0.37, p = 0.035), type 2 diabetes (rg = 0.33, p < 0.001), HDL (rg = - 0.41, p < 0.001), and coronary artery disease (rg = 0.43, p < 0.001). Local GC revealed some shared genetic loci responsible for the causality. The top locus 1126 was located at chromosome 7 and comprised genes HIBADH, JAZF1, and CREB5. The present study provides evidence for an independent causal effect of BMI on sepsis. Further detailed analysis of the shared genetic heritability between cardio-metabolic traits and sepsis provides the opportunity to improve the preventive strategies for sepsis.

Local patterns of genetic sharing challenge the boundaries between neuropsychiatric and insulin resistance-related conditions

The co-occurrence of insulin resistance (IR)-related metabolic conditions with neuropsychiatric disorders is a complex public health challenge. Evidence of the genetic links between these phenotypes is emerging, but little is currently known about the genomic regions and biological functions that are involved. To address this, we performed Local Analysis of [co]Variant Association (LAVA) using large-scale (N=9,725-933,970) genome-wide association studies (GWASs) results for three IR-related conditions (type 2 diabetes mellitus, obesity, and metabolic syndrome) and nine neuropsychiatric disorders. Subsequently, positional and expression quantitative trait locus (eQTL)-based gene mapping and downstream functional genomic analyses were performed on the significant loci. Patterns of negative and positive local genetic correlations (|rg|=0.21-1, pFDR<0.05) were identified at 109 unique genomic regions across all phenotype pairs. Local correlations emerged even in the absence of global genetic correlations between IR-related conditions and Alzheimer's disease, bipolar disorder, and Tourette's syndrome. Genes mapped to the correlated regions showed enrichment in biological pathways integral to immune-inflammatory function, vesicle trafficking, insulin signalling, oxygen transport, and lipid metabolism. Colocalisation analyses further prioritised 10 genetically correlated regions for likely harbouring shared causal variants, displaying high deleterious or regulatory potential. These variants were found within or in close proximity to genes, such as SLC39A8 and HLA-DRB1, that can be targeted by supplements and already known drugs, including omega-3/6 fatty acids, immunomodulatory, antihypertensive, and cholesterol-lowering drugs. Overall, our findings underscore the complex genetic landscape of IR-neuropsychiatric multimorbidity, advocating for an integrated disease model and offering novel insights for research and treatment strategies in this domain.

Using encrypted genotypes and phenotypes for collaborative genomic analyses to maintain data confidentiality

To adhere to and capitalize on the benefits of the FAIR (findable, accessible, interoperable, and reusable) principles in agricultural genome-to-phenome studies, it is crucial to address privacy and intellectual property issues that prevent sharing and reuse of data in research and industry. Direct sharing of genotype and phenotype data is often prohibited due to intellectual property and privacy concerns. Thus, there is a pressing need for encryption methods that obscure confidential aspects of the data, without affecting the outcomes of certain statistical analyses. A homomorphic encryption method for genotypes and phenotypes (HEGP) has been proposed for single-marker regression in genome-wide association studies (GWAS) using linear mixed models with Gaussian errors. This methodology permits frequentist likelihood-based parameter estimation and inference. In this paper, we extend HEGP to broader applications in genome-to-phenome analyses. We show that HEGP is suited to commonly used linear mixed models for genetic analyses of quantitative traits including genomic best linear unbiased prediction (GBLUP) and ridge-regression best linear unbiased prediction (RR-BLUP), as well as Bayesian variable selection methods (e.g. those in Bayesian Alphabet), for genetic parameter estimation, genomic prediction, and GWAS. By advancing the capabilities of HEGP, we offer researchers and industry professionals a secure and efficient approach for collaborative genomic analyses while preserving data confidentiality.

The genome-wide association study of serum IgE levels demonstrated a shared genetic background in allergic diseases

Immunoglobulin E (IgE) synthessis is highly related to a variety of atopic diseases, and several genome-wide association studies (GWASs) have demonstrated the association between genes and IgE level. In this study, we conducted the largest genome-wide association study of IgE involving a Taiwanese Han population. Eight independent variants exhibited genome-wide significance. Among them, an intronic SNP of CD28, rs1181388, and an intergenic SNP, rs1002957030, on 11q23.2 were identified as novel signals for IgE. Seven of the loci were replicated successfully in a meta-analysis using data on Japanese population. Among all the human leukocyte antigen (HLA) regions, HLA-DQA1*03:02 - HLA-DQB1*03:03 was the most significant haplotype (OR = 1.25, SE = 0.02, FDR = 1.6 × 10-14), corresponding to HLA-DQA1 Asp160 and HLA-DQB1 Leu87 amino acid residues. The genetic correlation showed significance between IgE and allergic diseases including asthma, atopic dermatitis, and pollinosis. IgE PRS was significantly correlated with total IgE levels. Furthermore, the top decile IgE polygenic risk score (PRS) group had the highest risk of asthma for the Taiwan Biobank and Biobank Japan cohorts. IgE PRS may be used to aid in predicting the occurrence of allergic reactions before symptoms occur and biomarkers are detectable. Our study provided a more comprehensive understanding of the impact of genomic variants, including complex HLA alleles, on serum IgE levels.

Shared genetic architectures of educational attainment in East Asian and European populations

Educational attainment (EduYears), a heritable trait often used as a proxy for cognitive ability, is associated with various health and social outcomes. Previous genome-wide association studies (GWASs) on EduYears have been focused on samples of European (EUR) genetic ancestries. Here we present the first large-scale GWAS of EduYears in people of East Asian (EAS) ancestry (n = 176,400) and conduct a cross-ancestry meta-analysis with EduYears GWAS in people of EUR ancestry (n = 766,345). EduYears showed a high genetic correlation and power-adjusted transferability ratio between EAS and EUR. We also found similar functional enrichment, gene expression enrichment and cross-trait genetic correlations between two populations. Cross-ancestry fine-mapping identified refined credible sets with a higher posterior inclusion probability than single population fine-mapping. Polygenic prediction analysis in four independent EAS and EUR cohorts demonstrated transferability between populations. Our study supports the need for further research on diverse ancestries to increase our understanding of the genetic basis of educational attainment.

Genetic correlations of alcohol consumption and alcohol use disorder with sex hormone levels in females and males

Background

Alcohol consumption behaviors and alcohol use disorder risk and presentation differ by sex, and these complex traits are associated with blood concentrations of the steroid sex hormones, testosterone and estradiol, and their regulatory binding proteins, sex hormone binding globulin (SHBG) and albumin. Genetic variation is associated with alcohol consumption and alcohol use disorder, as well as levels of steroid sex hormones and their binding proteins.

Methods

To assess the contribution of genetic factors to previously described phenotypic associations between alcohol-use traits and sex-hormone levels, we estimated genetic correlations (rg) using summary statistics from prior published, large sample size genome-wide association studies (GWAS) of alcohol consumption, alcohol dependence, testosterone, estradiol, SHBG, and albumin.

Results

For alcohol consumption, we observed positive genetic correlation (i.e. genetic effects in the same direction) with total testosterone in males (rg = 0.084, p = 0.007) and trends toward positive genetic correlation with bioavailable testosterone (rg = 0.060, p = 0.084) and SHBG in males (rg = 0.056, p = 0.086) and with albumin in a sex-combined cohort (rg = 0.082, p = 0.015); however in females, we observed positive genetic correlation with SHBG (rg = 0.089, p = 0.004) and a trend toward negative genetic correlation (i.e. genetic effects in opposite directions) with bioavailable testosterone (rg = -0.064, p = 0.032). For alcohol dependence, we observed a trend toward negative genetic correlation with total testosterone in females (rg = -0.106, p = 0.024) and positive genetic correlation with BMI-adjusted SHBG in males (rg = 0.119, p = 0.017). Several of these genetic correlations differed between females and males and were not in the same direction as the corresponding phenotypic associations.

Conclusions

Findings suggest that shared genetic effects may contribute to positive associations of alcohol consumption with albumin in both sexes, as well as positive associations between alcohol consumption and bioavailable testosterone and between alcohol dependence and SHBG in males. However, relative contributions of heritable and environmental factors to associations between alcohol-use traits and sex-hormone levels may differ by sex, with genetic factors contributing more in males and environmental factors contributing more in females.

Gene Discovery and Biological Insights into Anxiety Disorders from a Multi-Ancestry Genome-wide Association Study of >1.2 Million Participants

We leveraged information from more than 1.2 million participants to investigate the genetics of anxiety disorders across five continental ancestral groups. Ancestry-specific and cross-ancestry genome-wide association studies identified 51 anxiety-associated loci, 39 of which are novel. Additionally, polygenic risk scores derived from individuals of European descent were associated with anxiety in African, Admixed-American, and East Asian groups. The heritability of anxiety was enriched for genes expressed in the limbic system, the cerebral cortex, the cerebellum, the metencephalon, the entorhinal cortex, and the brain stem. Transcriptome- and proteome-wide analyses highlighted 115 genes associated with anxiety through brain-specific and cross-tissue regulation. We also observed global and local genetic correlations with depression, schizophrenia, and bipolar disorder and putative causal relationships with several physical health conditions. Overall, this study expands the knowledge regarding the genetic risk and pathogenesis of anxiety disorders, highlighting the importance of investigating diverse populations and integrating multi-omics information.

Heritability and genetic correlations for sleep apnea, insomnia, and hypersomnia in a large clinical biobank

RATIONALE

Comorbid insomnia and sleep apnea is reported to have worse outcomes than either condition alone. The local genetic correlations of these disorders are unknown.

OBJECTIVES

To identify local genomic regions with heritability for clinically diagnosed sleep apnea and insomnia, and to identify local genetic correlations between these disorders and/or hypersomnia.

METHODS

Fifty thousand two hundred seventeen patients of European ancestry were examined. Global and local heritability and genetic correlations for independent regions were calculated, adjusting for obesity and other covariates.

RESULTS

Sleep apnea and insomnia were significantly globally heritable and had 118 and 168 genetic regions with local heritability p-values <.05, respectively. One region had a significant genetic correlation for sleep apnea and hypersomnia (p-value = 9.85 × 10-4).

CONCLUSIONS

Clinically diagnosed sleep apnea and insomnia have minimal shared genetic architecture, supporting genetically distinct comorbid insomnia and sleep apnea components. However, additional correlated regions may be identified with additional sample size and methodological improvements.

Investigating the shared genetic architecture between COVID-19 and obesity: a large-scale genome wide cross-trait analysis

Observational studies have reported high comorbidity between obesity and severe COVID-19. The aim of this study is to explore whether genetic factors are involved in the co-occurrence of the two traits. Based on the available genome-wide association studies (GWAS) summary statistics, we explored the genetic correlation and performed cross-trait meta-analysis (CPASSOC) and colocalization analysis (COLOC) to detect pleiotropic single nucleotide polymorphisms (SNPs). At the genetic level, we obtained genes detected by Functional mapping and annotation (FUMA) and the Multi-marker Analysis of GenoMic Annotation (MAGMA). Potential functional genes were further investigated by summary-data-based Mendelian randomization (SMR). Finally, the casualty was identiied using the latent causal variable model (LCV). A significant positive genetic correlation was revealed between obesity and COVID-19. We found 331 shared genetic SNPs by CPASSOC and 13 shared risk loci by COLOC. At the genetic level, We obtained 3546 pleiotropic genes, among which 107 genes were found to be significantly expressed by SMR. Lastly, we observed these genes were mainly enriched in immune pathways and signaling transduction. These indings could provide new insights into the etiology of comorbidity and have implications for future therapeutic trial.

Distinct genomic signatures and modifiable risk factors underly 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. Using genomic data, this study elucidates biological mechanisms, key risk factors, and causal pathways underlying their comorbidity. We show that CVDs share a large proportion of their genetic risk factors with MDD. Multivariate genome-wide association analysis of the shared genetic liability between MDD and atherosclerotic CVD (ASCVD) revealed seven novel loci and distinct patterns of tissue and brain cell-type enrichments, suggesting a role for the thalamus. Part of the genetic overlap was explained by shared inflammatory, metabolic, and psychosocial/lifestyle risk factors. Finally, we found support for causal effects of genetic liability to MDD on CVD risk, but not from most CVDs to MDD, and demonstrated that the causal effects were partly explained by metabolic and psychosocial/lifestyle factors. The distinct signature of MDD-ASCVD comorbidity aligns with the idea of an immunometabolic sub-type of MDD more strongly associated with CVD than overall MDD. In summary, we identify plausible biological mechanisms underlying MDD-CVD comorbidity, as well as key modifiable risk factors for prevention of CVD in individuals with MDD.

A genome-wide investigation into the underlying genetic architecture of personality traits and overlap with psychopathology

Personality is influenced by both genetic and environmental factors and is associated with other psychiatric traits such as anxiety and depression. The "Big Five" personality traits, which include neuroticism, extraversion, agreeableness, conscientiousness, and openness, are a widely accepted and influential framework for understanding and describing human personality. Of the big five personality traits, neuroticism has most often been the focus of genetic studies and is linked to various mental illnesses including depression, anxiety, and schizophrenia. Our knowledge of the genetic architecture of the other four personality traits is more limited. Utilizing the Million Veteran Program (MVP) cohort we conducted a genome-wide association study (GWAS) in individuals of European and African ancestry. Adding other published data, we performed GWAS meta-analysis for each of the five personality traits with sample sizes ranging from 237,390 to 682,688. We identified 158, 14, 3, 2, and 7 independent genome-wide significant (GWS) loci associated with neuroticism, extraversion, agreeableness, conscientiousness, and openness, respectively. These findings represent 55 novel loci for neuroticism, as well as the first GWS loci discovered for extraversion and agreeableness. Gene-based association testing revealed 254 genes showing significant association with at least one of the five personality traits. Transcriptome-wide and proteome-wide analysis identified altered expression of genes and proteins such as CRHR1, SLC12A5, MAPT, and STX4. Pathway enrichment and drug perturbation analyses identified complex biology underlying human personality traits. We also studied the inter-relationship of personality traits with 1,437 other traits in a phenome-wide genetic correlation analysis, identifying new associations. Mendelian randomization showed positive bidirectional effects between neuroticism and depression and anxiety while a negative bidirectional effect was observed for agreeableness and these psychiatric traits. This study improves our comprehensive understanding of the genetic architecture underlying personality traits and their relationship to other complex human traits.

Genomic regions and biological pathways associated with sex-limited reproductive traits in bovine species

Many animal species exhibit sex-limited traits, where certain phenotypes are exclusively expressed in one sex. Yet, the genomic regions that contribute to these sex-limited traits in males and females remain a subject of debate. Reproductive traits are ideal phenotypes to study sexual differences since they are mostly expressed in a sex-limited way. Therefore, this study aims to use local correlation analyses to identify genomic regions and biological pathways significantly associated with male and female sex-limited traits in two distinct cattle breeds (Brahman [BB] and Tropical Composite [TC]). We used the Correlation Scan method to perform local correlation analysis on 42 trait pairs consisting of six female and seven male reproductive traits recorded on ~1,000 animals for each sex in each breed. To pinpoint a specific region associated with these sex-limited reproductive traits, we investigated the genomic region(s) consistently identified as significant across the 42 trait pairs in each breed. The genes found in the identified regions were subjected to Quantitative Trait Loci (QTL) colocalization, QTL enrichment analyses, and functional analyses to gain biological insight into sexual differences. We found that the genomic regions associated with the sex-limited reproductive phenotypes are widely distributed across all the chromosomes. However, no single region across the genome was associated with all the 42 reproductive trait pairs in the two breeds. Nevertheless, we found a region on the X-chromosome to be most significant for 80% to 90% (BB: 33 and TC: 38) of the total 42 trait pairs. A considerable number of the genes in this region were regulatory genes. By considering only genomic regions that were significant for at least 50% of the 42 trait pairs, we observed more regions spread across the autosomes and the X-chromosome. All genomic regions identified were highly enriched for trait-specific QTL linked to sex-limited traits (percentage of normal sperm, metabolic weight, average daily gain, carcass weight, age at puberty, etc.). The gene list created from these identified regions was enriched for biological pathways that contribute to the observed differences between sexes. Our results demonstrate that genomic regions associated with male and female sex-limited reproductive traits are distributed across the genome. Yet, chromosome X seems to exert a relatively larger effect on the phenotypic variation observed between the sexes.

Genetic Overlap Between Global Cortical Brain Structure, C-Reactive Protein, and White Blood Cell Counts

BACKGROUND

For many brain disorders, a subset of patients jointly exhibit alterations in cortical brain structure and elevated levels of circulating immune markers. This may be driven in part by shared genetic architecture. Therefore, we investigated the phenotypic and genetic associations linking global cortical surface area and thickness with blood immune markers (i.e., white blood cell counts and plasma C-reactive protein levels).

METHODS

Linear regression was used to assess phenotypic associations in 30,823 UK Biobank participants. Genome-wide and local genetic correlations were assessed using linkage disequilibrium score regression and local analysis of covariance annotation. The number of shared trait-influencing genetic variants was estimated using MiXeR. Shared genetic architecture was assessed using a conjunctional false discovery rate framework, and mapped genes were included in gene-set enrichment analyses.

RESULTS

Cortical structure and blood immune markers exhibited predominantly inverse phenotypic associations. There were modest genome-wide genetic correlations, the strongest of which were for C-reactive protein levels (rg_surface_area = -0.13, false discovery rate-corrected p = 4.17 × 10-3; rg_thickness = -0.13, false discovery rate-corrected p = 4.00 × 10-2). Meanwhile, local genetic correlations showed a mosaic of positive and negative associations. White blood cells shared on average 46.24% and 38.64% of trait-influencing genetic variants with surface area and thickness, respectively. Additionally, surface area shared 55 unique loci with the blood immune markers while thickness shared 15. Overall, monocyte count exhibited the largest genetic overlap with cortical brain structure. A series of gene enrichment analyses implicated neuronal-, astrocytic-, and schizophrenia-associated genes.

CONCLUSIONS

The findings indicate shared genetic underpinnings for cortical brain structure and blood immune markers, with implications for neurodevelopment and understanding the etiology of brain-related disorders.

Evidence for a shared genetic contribution to loneliness and borderline personality disorder

Loneliness, influenced by genetic and environmental factors such as childhood maltreatment, is one aspect of interpersonal dysfunction in Borderline Personality Disorder (BPD). Numerous studies link loneliness and BPD and twin studies indicate a genetic contribution to this association. The aim of our study was to investigate whether genetic predisposition for loneliness and BPD risk overlap and whether genetic risk for loneliness contributes to higher loneliness reported by BPD patients, using genome-wide genotype data. We assessed the genetic correlation of genome-wide association studies (GWAS) of loneliness and BPD using linkage disequilibrium score regression and tested whether a polygenic score for loneliness (loneliness-PGS) was associated with case-control status in two independent genotyped samples of BPD patients and healthy controls (HC; Witt2017-sample: 998 BPD, 1545 HC; KFO-sample: 187 BPD, 261 HC). In the KFO-sample, we examined associations of loneliness-PGS with reported loneliness, and whether the loneliness-PGS influenced the association between childhood maltreatment and loneliness. We found a genetic correlation between the GWAS of loneliness and BPD in the Witt2017-sample (rg = 0.23, p = 0.015), a positive association of loneliness-PGS with BPD case-control status (Witt2017-sample: NkR² = 2.3%, p = 2.7*10-12; KFO-sample: NkR² = 6.6%, p = 4.4*10-6), and a positive association between loneliness-PGS and loneliness across patient and control groups in the KFO-sample (β = 0.186, p = 0.002). The loneliness-PGS did not moderate the association between childhood maltreatment and loneliness in BPD. Our study is the first to use genome-wide genotype data to show that the genetic factors underlying variation in loneliness in the general population and the risk for BPD overlap. The loneliness-PGS was associated with reported loneliness. Further research is needed to investigate which genetic mechanisms and pathways are involved in this association and whether a genetic predisposition for loneliness contributes to BPD risk.

Inferring disease architecture and predictive ability with LDpred2-auto

LDpred2 is a widely used Bayesian method for building polygenic scores (PGSs). LDpred2-auto can infer the two parameters from the LDpred model, the SNP heritability h2 and polygenicity p, so that it does not require an additional validation dataset to choose best-performing parameters. The main aim of this paper is to properly validate the use of LDpred2-auto for inferring multiple genetic parameters. Here, we present a new version of LDpred2-auto that adds an optional third parameter α to its model, for modeling negative selection. We then validate the inference of these three parameters (or two, when using the previous model). We also show that LDpred2-auto provides per-variant probabilities of being causal that are well calibrated and can therefore be used for fine-mapping purposes. We also introduce a formula to infer the out-of-sample predictive performance r2 of the resulting PGS directly from the Gibbs sampler of LDpred2-auto. Finally, we extend the set of HapMap3 variants recommended to use with LDpred2 with 37% more variants to improve the coverage of this set, and we show that this new set of variants captures 12% more heritability and provides 6% more predictive performance, on average, in UK Biobank analyses.

Genetic relationship between the immune system and autism

Autism spectrum disorder (ASD) is a common and complex neurodevelopmental condition. The pathophysiology of ASD is poorly defined; however, it includes a strong genetic component and there is increasing evidence to support a role of immune dysregulation. Nonetheless, it is unclear which immune phenotypes link to ASD through genetics. Hence, we investigated the genetic correlation between ASD and diverse classes of immune conditions and markers; and if these immune-related genetic factors link to specific autistic-like traits in the population. We estimated global and local genetic correlations between ASD (n = 55,420) and 11 immune phenotypes (n = 14,256-755,406) using genome-wide association study summary statistics. Subsequently, polygenic scores (PGS) for these immune phenotypes were calculated in a population-based sample (n = 2487) and associated to five autistic-like traits (i.e., attention to detail, childhood behaviour, imagination, rigidity, social skills), and a total autistic-like traits score. Sex-stratified PGS analyses were also performed. At the genome-wide level, ASD was positively correlated with allergic diseases (ALG), and negatively correlated with lymphocyte count, rheumatoid arthritis (RA), and systemic lupus erythematosus (SLE) (FDR-p = 0.01-0.02). At the local genetic level, ASD was correlated with RA, C-reactive protein, and granulocytes and lymphocyte counts (p = 5.8 × 10-6-0.002). In the general population sample, increased genetic liability for SLE, RA, ALG, and lymphocyte levels, captured by PGS, was associated with the total autistic score and with rigidity and childhood behaviour (FDR-p = 0.03). In conclusion, we demonstrated a genetic relationship between ASD and immunity that depends on the type of immune phenotype considered; some increase likelihood whereas others may potentially help build resilience. Also, this relationship may be restricted to specific genetic loci and link to specific autistic dimensions (e.g., rigidity).

The genetic relationships between brain structure and schizophrenia

Genetic risks for schizophrenia are theoretically mediated by genetic effects on brain structure but it has been unclear which genes are associated with both schizophrenia and cortical phenotypes. We accessed genome-wide association studies (GWAS) of schizophrenia (N = 69,369 cases; 236,642 controls), and of three magnetic resonance imaging (MRI) metrics (surface area, cortical thickness, neurite density index) measured at 180 cortical areas (N = 36,843, UK Biobank). Using Hi-C-coupled MAGMA, 61 genes were significantly associated with both schizophrenia and one or more MRI metrics. Whole genome analysis with partial least squares demonstrated significant genetic covariation between schizophrenia and area or thickness of most cortical regions. Genetic similarity between cortical areas was strongly coupled to their phenotypic covariance, and genetic covariation between schizophrenia and brain phenotypes was strongest in the hubs of structural covariance networks. Pleiotropically associated genes were enriched for neurodevelopmental processes and positionally concentrated in chromosomes 3p21, 17q21 and 11p11. Mendelian randomization analysis indicated that genetically determined variation in a posterior cingulate cortical area could be causal for schizophrenia. Parallel analyses of GWAS on bipolar disorder, Alzheimer's disease and height showed that pleiotropic association with MRI metrics was stronger for schizophrenia compared to other disorders.

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, and efforts to understand the biological basis of cognitive dysfunction in the disorder are ongoing. Previous studies have suggested genetic overlap between global cognitive ability and schizophrenia, but further work is needed to delineate the shared genetic architecture. Here, we apply genomic structural equation modelling to identify latent cognitive factors capturing genetic liabilities to 12 cognitive traits measured in the UK Biobank (UKB). 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). We identified three broad factors (visuo-spatial, verbal analytic and decision/reaction time) that underly the genetic correlations between the UKB cognitive tests. Global genetic correlations showed a significant but 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 evidence of substantial polygenic overlap between each cognitive factor and schizophrenia but show that most loci shared between the latent cognitive factors and schizophrenia have unique patterns of association with the cognitive factors. Biological annotation of the shared loci implicated gene-sets related to neurodevelopment and neuronal function. Lastly, we find that the common genetic determinants of the latent cognitive factors are not predictive of schizophrenia symptom dimensions. 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.

Factorizing polygenic epistasis improves prediction and uncovers biological pathways in complex traits

Epistasis is central in many domains of biology, but it has not yet been proven useful for understanding the etiology of complex traits. This is partly because complex-trait epistasis involves polygenic interactions that are poorly captured in current models. To address this gap, we developed a model called Epistasis Factor Analysis (EFA). EFA assumes that polygenic epistasis can be factorized into interactions between a few epistasis factors (EFs), which represent latent polygenic components of the observed complex trait. The statistical goals of EFA are to improve polygenic prediction and to increase power to detect epistasis, while the biological goal is to unravel genetic effects into more-homogeneous units. We mathematically characterize EFA and use simulations to show that EFA outperforms current epistasis models when its assumptions approximately hold. Applied to predicting yeast growth rates, EFA outperforms the additive model for several traits with large epistasis heritability and uniformly outperforms the standard epistasis model. We replicate these prediction improvements in a second dataset. We then apply EFA to four previously characterized traits in the UK Biobank and find statistically significant epistasis in all four, including two that are robust to scale transformation. Moreover, we find that the inferred EFs partly recover pre-defined biological pathways for two of the traits. Our results demonstrate that more realistic models can identify biologically and statistically meaningful epistasis in complex traits, indicating that epistasis has potential for precision medicine and characterizing the biology underlying GWAS results.

The shared genetic architecture and evolution of human language and musical rhythm

Rhythm and language-related traits are phenotypically correlated, but their genetic overlap is largely unknown. Here, we leveraged two large-scale genome-wide association studies performed to shed light on the shared genetics of rhythm (N=606,825) and dyslexia (N=1,138,870). Our results reveal an intricate shared genetic and neurobiological architecture, and lay groundwork for resolving longstanding debates about the potential co-evolution of human language and musical traits.

Shifting the trajectory of therapeutic development for neurological and psychiatric disorders

Clinical trials for central nervous system disorders often enroll patients with unrecognized heterogeneous diseases, leading to costly trials that have high failure rates. Here, we discuss the potential of emerging technologies and datasets to elucidate disease mechanisms and identify biomarkers to improve patient stratification and monitoring of disease progression in clinical trials for neuropsychiatric disorders. Greater efforts must be centered on rigorously standardizing data collection and sharing of methods, datasets, and analytical tools across sectors. To address health care disparities in clinical trials, diversity of genetic ancestries and environmental exposures of research participants and associated biological samples must be prioritized.

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.

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.

Trait selection strategy in multi-trait GWAS: Boosting SNPs discoverability

Since the first Genome-Wide Association Studies (GWAS), thousands of variant-trait associations have been discovered. However, the sample size required to detect additional variants using standard univariate association screening is increasingly prohibitive. Multi-trait GWAS offers a relevant alternative: it can improve statistical power and lead to new insights about gene function and the joint genetic architecture of human phenotypes. Although many methodological hurdles of multi-trait testing have been discussed, the strategy to select trait, among overwhelming possibilities, has been overlooked. In this study, we conducted extensive multi-trait tests using JASS (Joint Analysis of Summary Statistics) and assessed which genetic features of the analysed sets were associated with an increased detection of variants as compared to univariate screening. Our analyses identified multiple factors associated with the gain in the association detection in multi-trait tests. Together, these factors of the analysed sets are predictive of the gain of the multi-trait test (Pearson's ρ equal to 0.43 between the observed and predicted gain, P < 1.6 × 10-60). Applying an alternative multi-trait approach (MTAG, multi-trait analysis of GWAS), we found that in most scenarios but particularly those with larger numbers of traits, JASS outperformed MTAG. Finally, we benchmark several strategies to select set of traits including the prevalent strategy of selecting clinically similar traits, which systematically underperformed selecting clinically heterogenous traits or selecting sets that issued from our data-driven models. This work provides a unique picture of the determinant of multi-trait GWAS statistical power and outline practical strategies for multi-trait testing.

A parametric bootstrap approach for computing confidence intervals for genetic correlations with application to genetically-determined protein-protein networks

Genetic correlation refers to the correlation between genetic determinants of a pair of traits. When using individual-level data, it is typically estimated based on a bivariate model specification where the correlation between the two variables is identifiable and can be estimated from a covariance model that incorporates the genetic relationship between individuals, e.g., using a pre-specified kinship matrix. Inference relying on asymptotic normality of the genetic correlation parameter estimates may be inaccurate when the sample size is low, when the genetic correlation is close to the boundary of the parameter space, and when the heritability of at least one of the traits is low. We address this problem by developing a parametric bootstrap procedure to construct confidence intervals for genetic correlation estimates. The procedure simulates paired traits under a range of heritability and genetic correlation parameters, and it uses the population structure encapsulated by the kinship matrix. Heritabilities and genetic correlations are estimated using the close-form, method of moment, Haseman-Elston regression estimators. The proposed parametric bootstrap procedure is especially useful when genetic correlations are computed on pairs of thousands of traits measured on the same exact set of individuals. We demonstrate the parametric bootstrap approach on a proteomics dataset from the Jackson Heart Study.