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.

Blood-based DNA methylation and exposure risk scores predict PTSD with high accuracy in military and civilian cohorts.

Background: Incorporating genomic data into risk prediction has become an increasingly useful approach for rapid identification of individuals most at risk for complex disorders such as PTSD. Our goal was to develop and validate Methylation Risk Scores (MRS) using machine learning to distinguish individuals who have PTSD from those who do not. Methods: Elastic Net was used to develop three risk score models using a discovery dataset (n = 1226; 314 cases, 912 controls) comprised of 5 diverse cohorts with available blood-derived DNA methylation (DNAm) measured on the Illumina Epic BeadChip. The first risk score, exposure and methylation risk score (eMRS) used cumulative and childhood trauma exposure and DNAm variables; the second, methylation-only risk score (MoRS) was based solely on DNAm data; the third, methylation-only risk scores with adjusted exposure variables (MoRSAE) utilized DNAm data adjusted for the two exposure variables. The potential of these risk scores to predict future PTSD based on pre-deployment data was also assessed. External validation of risk scores was conducted in four independent cohorts. Results: The eMRS model showed the highest accuracy (92%), precision (91%), recall (87%), and f1-score (89%) in classifying PTSD using 3730 features. While still highly accurate, the MoRS (accuracy = 89%) using 3728 features and MoRSAE (accuracy = 84%) using 4150 features showed a decline in classification power. eMRS significantly predicted PTSD in one of the four independent cohorts, the BEAR cohort (beta = 0.6839, p-0.003), but not in the remaining three cohorts. Pre-deployment risk scores from all models (eMRS, beta = 1.92; MoRS, beta = 1.99 and MoRSAE, beta = 1.77) displayed a significant (p < 0.001) predictive power for post-deployment PTSD. Conclusion: Results, especially those from the eMRS, reinforce earlier findings that methylation and trauma are interconnected and can be leveraged to increase the correct classification of those with vs. without PTSD. Moreover, our models can potentially be a valuable tool in predicting the future risk of developing PTSD. As more data become available, including additional molecular, environmental, and psychosocial factors in these scores may enhance their accuracy in predicting the condition and, relatedly, improve their performance in independent cohorts.

Whole-exome sequencing in UK Biobank reveals rare genetic architecture for depression.

Nearly two hundred common-variant depression risk loci have been identified by genome-wide association studies (GWAS). However, the impact of rare coding variants on depression remains poorly understood. Here, we present whole-exome sequencing analyses of depression with seven different definitions based on survey, questionnaire, and electronic health records in 320,356 UK Biobank participants. We showed that the burden of rare damaging coding variants in loss-of-function intolerant genes is significantly associated with risk of depression with various definitions. We compared the rare and common genetic architecture across depression definitions by genetic correlation and showed different genetic relationships between definitions across common and rare variants. In addition, we demonstrated that the effects of rare damaging coding variant burden and polygenic risk score on depression risk are additive. The gene set burden analyses revealed overlapping rare genetic variant components with developmental disorder, autism, and schizophrenia. Our study provides insights into the contribution of rare coding variants, separately and in conjunction with common variants, on depression with various definitions and their genetic relationships with neurodevelopmental disorders.

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.

Large-Scale Mendelian Randomization Study Reveals Circulating Blood-based Proteomic Biomarkers for Psychopathology and Cognitive Task Performance.

Background: Research on peripheral (e.g., blood-based) biomarkers for psychiatric illness has typically been low-throughput in terms of both the number of subjects and the range of assays performed. Moreover, traditional case-control studies examining blood-based biomarkers are subject to potential confounds of treatment and other exposures common to patients with psychiatric illnesses. Our research addresses these challenges by leveraging large-scale, high-throughput proteomics data and Mendelian Randomization (MR) to examine the causal impact of circulating proteins on psychiatric phenotypes and cognitive task performance. Methods: We utilized plasma proteomics data from the UK Biobank (3,072 proteins assayed in 34,557 European-ancestry individuals) and deCODE Genetics (4,719 proteins measured across 35,559 Icelandic individuals). Significant proteomic quantitative trait loci (both cis-pQTLs and trans-pQTLs) served as MR instruments, with the most recent GWAS for schizophrenia, bipolar disorder, major depressive disorder, and cognitive task performance (all excluding overlapping UK Biobank participants) as phenotypic outcomes. Results: MR revealed 109 Bonferroni-corrected causal associations (44 novel) involving 88 proteins across the four phenotypes. Several immune-related proteins, including interleukins and complement factors, stood out as pleiotropic across multiple outcome phenotypes. Drug target enrichment analysis identified several novel potential pharmacologic repurposing opportunities, including anti-inflammatory agents for schizophrenia and bipolar disorder and duloxetine for cognitive performance. Conclusions: Identification of causal effects for these circulating proteins suggests potential biomarkers for these conditions and offers insights for developing innovative therapeutic strategies. The findings also indicate substantial evidence for the pleiotropic effects of many proteins across different phenotypes, shedding light on the shared etiology among psychiatric conditions and cognitive ability.

Analysis across Taiwan Biobank, Biobank Japan, and UK Biobank identifies hundreds of novel loci for 36 quantitative traits.

Genome-wide association studies (GWASs) have identified tens of thousands of genetic loci associated with human complex traits. However, the majority of GWASs were conducted in individuals of European ancestries. Failure to capture global genetic diversity has limited genomic discovery and has impeded equitable delivery of genomic knowledge to diverse populations. Here we report findings from 102,900 individuals across 36 human quantitative traits in the Taiwan Biobank (TWB), a major biobank effort that broadens the population diversity of genetic studies in East Asia. We identified 968 novel genetic loci, pinpointed novel causal variants through statistical fine-mapping, compared the genetic architecture across TWB, Biobank Japan, and UK Biobank, and evaluated the utility of cross-phenotype, cross-population polygenic risk scores in disease risk prediction. These results demonstrated the potential to advance discovery through diversifying GWAS populations and provided insights into the common genetic basis of human complex traits in East Asia.

Plasma proteomic associations with genetics and health in the UK Biobank.

The Pharma Proteomics Project is a precompetitive biopharmaceutical consortium characterizing the plasma proteomic profiles of 54,219 UK Biobank participants. Here we provide a detailed summary of this initiative, including technical and biological validations, insights into proteomic disease signatures, and prediction modelling for various demographic and health indicators. We present comprehensive protein quantitative trait locus (pQTL) mapping of 2,923 proteins that identifies 14,287 primary genetic associations, of which 81% are previously undescribed, alongside ancestry-specific pQTL mapping in non-European individuals. The study provides an updated characterization of the genetic architecture of the plasma proteome, contextualized with projected pQTL discovery rates as sample sizes and proteomic assay coverages increase over time. We offer extensive insights into trans pQTLs across multiple biological domains, highlight genetic influences on ligand-receptor interactions and pathway perturbations across a diverse collection of cytokines and complement networks, and illustrate long-range epistatic effects of ABO blood group and FUT2 secretor status on proteins with gastrointestinal tissue-enriched expression. We demonstrate the utility of these data for drug discovery by extending the genetic proxied effects of protein targets, such as PCSK9, on additional endpoints, and disentangle specific genes and proteins perturbed at loci associated with COVID-19 susceptibility. This public-private partnership provides the scientific community with an open-access proteomics resource of considerable breadth and depth to help to elucidate the biological mechanisms underlying proteo-genomic discoveries and accelerate the development of biomarkers, predictive models and therapeutics1.

Discovery of 95 PTSD loci provides insight into genetic architecture and neurobiology of trauma and stress-related disorders.

Posttraumatic 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 novel). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (e.g., GRIA1, GRM8, CACNA1E ), developmental, axon guidance, and transcription factors (e.g., FOXP2, EFNA5, DCC ), synaptic structure and function genes (e.g., PCLO, NCAM1, PDE4B ), and endocrine or immune regulators (e.g., ESR1, TRAF3, 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.

Causal influence of sleeping phenotypes on the risk of coronary artery disease and sudden cardiac arrest: A Mendelian randomization analysis.

OBJECTIVES: To assess the causal influence of sleep and circadian traits on coronary artery disease and sudden cardiac arrest with adjustment for obesity through a two-sample Mendelian randomization study. METHODS: We used summary statistics of 5 sleep and circadian traits for genome-wide association studies, including chronotype, sleep duration, long sleep (≥9 h a day), short sleep (<7 h a day), and insomnia (sample size range: 237,622-651,295). Coronary artery disease genome-wide association studies with 60,801 cases and 123,504 controls, sudden cardiac arrest genome-wide association studies with 3939 cases and 25,989 controls, and obesity genome-wide association studies with 806,834 individuals were also used. Multivariable Mendelian randomization was performed to estimate the causality. RESULTS: After adjusting for obesity, genetically predicted short sleep (odds ratio = 1.87 and p = .02), and genetically predicted insomnia (odds ratio = 1.17 and p = .001) were causally associated with increased odds of coronary artery disease. Genetically predicted long sleep (odds ratio = 0.06 and p = .02) and genetically predicted longer sleep duration (odds ratio = 0.36 for per-hour increase in sleep duration and p = .0006) were causally associated with decreased odds of sudden cardiac arrest. CONCLUSIONS: The findings of this Mendelian randomization study indicate that insomnia and short sleep contribute to the development of coronary artery disease, whereas a longer sleep duration protects from sudden cardiac arrest, independent of the influence of obesity. The mechanisms underlying these associations warrant further investigation.

Rare genetic variants impact muscle strength.

Muscle strength is highly heritable and predictive for multiple adverse health outcomes including mortality. Here, we present a rare protein-coding variant association study in 340,319 individuals for hand grip strength, a proxy measure of muscle strength. We show that the exome-wide burden of rare protein-truncating and damaging missense variants is associated with a reduction in hand grip strength. We identify six significant hand grip strength genes, KDM5B, OBSCN, GIGYF1, TTN, RB1CC1, and EIF3J. In the example of the titin (TTN) locus we demonstrate a convergence of rare with common variant association signals and uncover genetic relationships between reduced hand grip strength and disease. Finally, we identify shared mechanisms between brain and muscle function and uncover additive effects between rare and common genetic variation on muscle strength.

The impact of rare protein coding genetic variation on adult cognitive function.

Compelling evidence suggests that human cognitive function is strongly influenced by genetics. Here, we conduct a large-scale exome study to examine whether rare protein-coding variants impact cognitive function in the adult population (n = 485,930). We identify eight genes (ADGRB2, KDM5B, GIGYF1, ANKRD12, SLC8A1, RC3H2, CACNA1A and BCAS3) that are associated with adult cognitive function through rare coding variants with large effects. Rare genetic architecture for cognitive function partially overlaps with that of neurodevelopmental disorders. In the case of KDM5B we show how the genetic dosage of one of these genes may determine the variability of cognitive, behavioral and molecular traits in mice and humans. We further provide evidence that rare and common variants overlap in association signals and contribute additively to cognitive function. Our study introduces the relevance of rare coding variants for cognitive function and unveils high-impact monogenic contributions to how cognitive function is distributed in the normal adult population.

Brain expression quantitative trait locus and network analyses reveal downstream effects and putative drivers for brain-related diseases.

Identification of therapeutic targets from genome-wide association studies (GWAS) requires insights into downstream functional consequences. We harmonized 8,613 RNA-sequencing samples from 14 brain datasets to create the MetaBrain resource and performed cis- and trans-expression quantitative trait locus (eQTL) meta-analyses in multiple brain region- and ancestry-specific datasets (n ≤ 2,759). Many of the 16,169 cortex cis-eQTLs were tissue-dependent when compared with blood cis-eQTLs. We inferred brain cell types for 3,549 cis-eQTLs by interaction analysis. We prioritized 186 cis-eQTLs for 31 brain-related traits using Mendelian randomization and co-localization including 40 cis-eQTLs with an inferred cell type, such as a neuron-specific cis-eQTL (CYP24A1) for multiple sclerosis. We further describe 737 trans-eQTLs for 526 unique variants and 108 unique genes. We used brain-specific gene-co-regulation networks to link GWAS loci and prioritize additional genes for five central nervous system diseases. This study represents a valuable resource for post-GWAS research on central nervous system diseases.

FinnGen provides genetic insights from a well-phenotyped isolated population.

Population isolates such as those in Finland benefit genetic research because deleterious alleles are often concentrated on a small number of low-frequency variants (0.1% ≤ minor allele frequency < 5%). These variants survived the founding bottleneck rather than being distributed over a large number of ultrarare variants. Although this effect is well established in Mendelian genetics, its value in common disease genetics is less explored1,2. FinnGen aims to study the genome and national health register data of 500,000 Finnish individuals. Given the relatively high median age of participants (63 years) and the substantial fraction of hospital-based recruitment, FinnGen is enriched for disease end points. Here we analyse data from 224,737 participants from FinnGen and study 15 diseases that have previously been investigated in large genome-wide association studies (GWASs). We also include meta-analyses of biobank data from Estonia and the United Kingdom. We identified 30 new associations, primarily low-frequency variants, enriched in the Finnish population. A GWAS of 1,932 diseases also identified 2,733 genome-wide significant associations (893 phenome-wide significant (PWS), P < 2.6 × 10-11) at 2,496 (771 PWS) independent loci with 807 (247 PWS) end points. Among these, fine-mapping implicated 148 (73 PWS) coding variants associated with 83 (42 PWS) end points. Moreover, 91 (47 PWS) had an allele frequency of <5% in non-Finnish European individuals, of which 62 (32 PWS) were enriched by more than twofold in Finland. These findings demonstrate the power of bottlenecked populations to find entry points into the biology of common diseases through low-frequency, high impact variants.

GWAS quality score for evaluating associated regions in GWAS analyses.

MOTIVATION: The number of significantly associated regions reported in genome-wide association studies (GWAS) for polygenic traits typically increases with sample size. A traditional tool for quality control and identification of significant regions has been a visual inspection of how significant and correlated genetic variants cluster within a region. However, while inspecting hundreds of regions, this subjective method can misattribute significance to some loci or neglect others that are significant. RESULTS: The GWAS quality score (GQS) identifies suspicious regions and prevents erroneous interpretations with an objective, quantitative and automated method. The GQS assesses all measured single nucleotide polymorphisms (SNPs) that are linked by inheritance to each other [linkage disequilibrium (LD)] and compares the significance of trait association of each SNP to its LD value for the reported index SNP. A GQS value of 1.0 ascribes a high level of confidence to the entire region and its underlying gene(s), while GQS values <1.0 indicate the need to closely inspect the outliers. We applied the GQS to published and non-published genome-wide summary statistics and report suspicious regions requiring secondary inspection while supporting the majority of reported regions from large-scale published meta-analyses. AVAILABILITY AND IMPLEMENTATION: The GQS code/scripts can be cloned from GitHub (https://github.com/Xswapnil/GQS/). The analyst can use whole-genome summary statistics to estimate GQS for each defined region. We also provide an online tool (http://35.227.18.38/) that gives access to the GQS. The quantitative measure of quality attributes by GQS and its visualization is an objective method that enhances the confidence of each genomic hit. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Fine-mapping across diverse ancestries drives the discovery of putative causal variants underlying human complex traits and diseases.

Genome-wide association studies (GWAS) of human complex traits or diseases often implicate genetic loci that span hundreds or thousands of genetic variants, many of which have similar statistical significance. While statistical fine-mapping in individuals of European ancestries has made important discoveries, cross-population fine-mapping has the potential to improve power and resolution by capitalizing on the genomic diversity across ancestries. Here we present SuSiEx, an accurate and computationally efficient method for cross-population fine-mapping, which builds on the single-population fine-mapping framework, Sum of Single Effects (SuSiE). SuSiEx integrates data from an arbitrary number of ancestries, explicitly models population-specific allele frequencies and LD patterns, accounts for multiple causal variants in a genomic region, and can be applied to GWAS summary statistics. We comprehensively evaluated SuSiEx using simulations, a range of quantitative traits measured in both UK Biobank and Taiwan Biobank, and schizophrenia GWAS across East Asian and European ancestries. In all evaluations, SuSiEx fine-mapped more association signals, produced smaller credible sets and higher posterior inclusion probability (PIP) for putative causal variants, and captured population-specific causal variants.

Associations of polygenic risks, depression, and obesity-related traits in Taiwan Biobank.

BACKGROUND: The comorbidity of obesity and major depressive disorder (MDD) may be attributable to a bidirectional relationship and shared genetic influence. We aimed to examine the polygenic associations between obesity and MDD and to characterize their corresponding impacts on the obesity mechanism. METHODS: Genome-wide genotyping was available in 106,604 unrelated individuals from Taiwan Biobank. Polygenic risk score (PRS) for body mass index (BMI) and MDD was derived to evaluate their effects on obesity-related traits. Stratified analyses were performed for the modified effect of depression on the polygenic associations. RESULTS: The MDD PRS was positively associated with waistline (beta in per SD increase in PRS = 0.12), hipline (beta = 0.08), waist-hip ratio (WHR) (beta = 0.05), body fat rate (beta = 0.08), BMI (beta = 0.05), overweight (OR = 1.02 for BMI ≥ 25), and obesity (OR = 1.05 for BMI ≥ 30). For the synergism between depression and BMI PRS, the presence of active depression symptoms defined by the PHQ-4 (p for interaction < 0.05 for waistline, WHR, and BMI) was more salient than lifetime MDD. LIMITATIONS: Limitations include recall bias for MDD due to a retrospective self-reporting questionnaire, a low response rate of the PHQ-4 for evaluating active psychological symptoms, and limited generalizability to non-Taiwanese ancestries. CONCLUSIONS: The shared genetic etiology of obesity and depression was demonstrated. The amplified effect of BMI polygenic effect on obesity for individuals with active depressive symptoms was also characterized. The study may be helpful for designing public health interventions to reduce the disease burden caused by obesity and depression.

The accuracy of polygenic score models for anthropometric traits and Type II Diabetes in the Native Hawaiian Population.

Polygenic scores (PGS) are promising in stratifying individuals based on the genetic susceptibility to complex diseases or traits. However, the accuracy of PGS models, typically trained in European- or East Asian-ancestry populations, tend to perform poorly in other ethnic minority populations, and their accuracies have not been evaluated for Native Hawaiians. Using body mass index, height, and type-2 diabetes as examples of highly polygenic traits, we evaluated the prediction accuracies of PGS models in a large Native Hawaiian sample from the Multiethnic Cohort with up to 5,300 individuals. We evaluated both publicly available PGS models or genome-wide PGS models trained in this study using the largest available GWAS. We found evidence of lowered prediction accuracies for the PGS models in some cases, particularly for height. We also found that using the Native Hawaiian samples as an optimization cohort during training did not consistently improve PGS performance. Moreover, even the best performing PGS models among Native Hawaiians would have lowered prediction accuracy among the subset of individuals most enriched with Polynesian ancestry. Our findings indicate that factors such as admixture histories, sample size and diversity in GWAS can influence PGS performance for complex traits among Native Hawaiian samples. This study provides an initial survey of PGS performance among Native Hawaiians and exposes the current gaps and challenges associated with improving polygenic prediction models for underrepresented minority populations.

Collective genomic segments with differential pleiotropic patterns between cognitive dimensions and psychopathology.

Cognitive deficits are known to be related to most forms of psychopathology. Here, we perform local genetic correlation analysis as a means of identifying independent segments of the genome that show biologically interpretable pleiotropic associations between cognitive dimensions and psychopathology. We identify collective segments of the genome, which we call "meta-loci", showing differential pleiotropic patterns for psychopathology relative to either cognitive task performance (CTP) or performance on a non-cognitive factor (NCF) derived from educational attainment. We observe that neurodevelopmental gene sets expressed during the prenatal-early childhood period predominate in CTP-relevant meta-loci, while post-natal gene sets are more involved in NCF-relevant meta-loci. Further, we demonstrate that neurodevelopmental gene sets are dissociable across CTP meta-loci with respect to their spatial distribution across the brain. Additionally, we find that GABA-ergic, cholinergic, and glutamatergic genes drive pleiotropic relationships within dissociable meta-loci.

Evaluation of the causal relationship between smoking and schizophrenia in East Asia.

Cigarette smoking has been suggested to be associated with the risk of schizophrenia in observational studies. A significant causal effect of smoking on schizophrenia has been reported in European populations using the Mendelian randomization approach; however, no evidence of causality was found in participants from East Asia. Using Taiwan Biobank (TWBB), we conducted genome-wide association studies (GWAS) to identify susceptibility loci for smoking behaviors, including smoking initiation (N = 79,989) and the onset age (N = 15,582). We then meta-analyzed GWAS from TWBB and Biobank Japan (BBJ) with the total sample size of 245,425 for smoking initiation and 46,000 for onset age of smoking. The GWAS for schizophrenia was taken from the East Asia Psychiatric Genomics Consortium, which included 22,778 cases and 35,362 controls. We performed a two-sample Mendelian randomization to estimate the causality of smoking behaviors on schizophrenia in East Asia. In TWBB, we identified one locus that met genome-wide significance for onset age. In a meta-analysis of TWBB and BBJ, we identified two loci for smoking initiation. In Mendelian randomization, genetically predicted smoking initiation (odds ratio (OR) = 4.00, 95% confidence interval (CI) = 0.89-18.01, P = 0.071) and onset age (OR for a per-year increase = 0.96, 95% CI = 0.91-1.01, P = 0.098) were not significantly associated with schizophrenia; the direction of effect was consistent with European Ancestry samples, which had higher statistical power. These findings provide tentative evidence consistent with a causal role of smoking on the development of schizophrenia in East Asian populations.