Genome-wide Association Studies in Ancestrally Diverse Populations: Opportunities, Methods, Pitfalls, and Recommendations.

Genome-wide association studies (GWASs) have focused primarily on populations of European descent, but it is essential that diverse populations become better represented. Increasing diversity among study participants will advance our understanding of genetic architecture in all populations and ensure that genetic research is broadly applicable. To facilitate and promote research in multi-ancestry and admixed cohorts, we outline key methodological considerations and highlight opportunities, challenges, solutions, and areas in need of development. Despite the perception that analyzing genetic data from diverse populations is difficult, it is scientifically and ethically imperative, and there is an expanding analytical toolbox to do it well.

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.

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.

Genetic associations of protein-coding variants in human disease.

Genome-wide association studies (GWAS) have identified thousands of genetic variants linked to the risk of human disease. However, GWAS have so far remained largely underpowered in relation to identifying associations in the rare and low-frequency allelic spectrum and have lacked the resolution to trace causal mechanisms to underlying genes1. Here we combined whole-exome sequencing in 392,814 UK Biobank participants with imputed genotypes from 260,405 FinnGen participants (653,219 total individuals) to conduct association meta-analyses for 744 disease endpoints across the protein-coding allelic frequency spectrum, bridging the gap between common and rare variant studies. We identified 975 associations, with more than one-third being previously unreported. We demonstrate population-level relevance for mutations previously ascribed to causing single-gene disorders, map GWAS associations to likely causal genes, explain disease mechanisms, and systematically relate disease associations to levels of 117 biomarkers and clinical-stage drug targets. Combining sequencing and genotyping in two population biobanks enabled us to benefit from increased power to detect and explain disease associations, validate findings through replication and propose medical actionability for rare genetic variants. Our study provides a compendium of protein-coding variant associations for future insights into disease biology and drug discovery.

Mapping genomic loci implicates genes and synaptic biology in schizophrenia.

Schizophrenia has a heritability of 60-80%1, much of which is attributable to common risk alleles. Here, in a two-stage genome-wide association study of up to 76,755 individuals with schizophrenia and 243,649 control individuals, we report common variant associations at 287 distinct genomic loci. Associations were concentrated in genes that are expressed in excitatory and inhibitory neurons of the central nervous system, but not in other tissues or cell types. Using fine-mapping and functional genomic data, we identify 120 genes (106 protein-coding) that are likely to underpin associations at some of these loci, including 16 genes with credible causal non-synonymous or untranslated region variation. We also implicate fundamental processes related to neuronal function, including synaptic organization, differentiation and transmission. Fine-mapped candidates were enriched for genes associated with rare disruptive coding variants in people with schizophrenia, including the glutamate receptor subunit GRIN2A and transcription factor SP4, and were also enriched for genes implicated by such variants in neurodevelopmental disorders. We identify biological processes relevant to schizophrenia pathophysiology; show convergence of common and rare variant associations in schizophrenia and neurodevelopmental disorders; and provide a resource of prioritized genes and variants to advance mechanistic studies.

Accumulation of the inner nuclear envelope protein Sun1 is pathogenic in progeric and dystrophic laminopathies.

Human LMNA gene mutations result in laminopathies that include Emery-Dreifuss muscular dystrophy (AD-EDMD) and Hutchinson-Gilford progeria, the premature aging syndrome (HGPS). The Lmna null (Lmna(-/-)) and progeroid LmnaΔ9 mutant mice are models for AD-EDMD and HGPS, respectively. Both animals develop severe tissue pathologies with abbreviated life spans. Like HGPS cells, Lmna(-/-) and LmnaΔ9 fibroblasts have typically misshapen nuclei. Unexpectedly, Lmna(-/-) or LmnaΔ9 mice that are also deficient for the inner nuclear membrane protein Sun1 show markedly reduced tissue pathologies and enhanced longevity. Concordantly, reduction of SUN1 overaccumulation in LMNA mutant fibroblasts and in cells derived from HGPS patients corrected nuclear defects and cellular senescence. Collectively, these findings implicate Sun1 protein accumulation as a common pathogenic event in Lmna(-/-), LmnaΔ9, and HGPS disorders.

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.

Identifying drug targets for neurological and psychiatric disease via genetics and the brain transcriptome.

Discovering drugs that efficiently treat brain diseases has been challenging. Genetic variants that modulate the expression of potential drug targets can be utilized to assess the efficacy of therapeutic interventions. We therefore employed Mendelian Randomization (MR) on gene expression measured in brain tissue to identify drug targets involved in neurological and psychiatric diseases. We conducted a two-sample MR using cis-acting brain-derived expression quantitative trait loci (eQTLs) from the Accelerating Medicines Partnership for Alzheimer's Disease consortium (AMP-AD) and the CommonMind Consortium (CMC) meta-analysis study (n = 1,286) as genetic instruments to predict the effects of 7,137 genes on 12 neurological and psychiatric disorders. We conducted Bayesian colocalization analysis on the top MR findings (using P<6x10-7 as evidence threshold, Bonferroni-corrected for 80,557 MR tests) to confirm sharing of the same causal variants between gene expression and trait in each genomic region. We then intersected the colocalized genes with known monogenic disease genes recorded in Online Mendelian Inheritance in Man (OMIM) and with genes annotated as drug targets in the Open Targets platform to identify promising drug targets. 80 eQTLs showed MR evidence of a causal effect, from which we prioritised 47 genes based on colocalization with the trait. We causally linked the expression of 23 genes with schizophrenia and a single gene each with anorexia, bipolar disorder and major depressive disorder within the psychiatric diseases and 9 genes with Alzheimer's disease, 6 genes with Parkinson's disease, 4 genes with multiple sclerosis and two genes with amyotrophic lateral sclerosis within the neurological diseases we tested. From these we identified five genes (ACE, GPNMB, KCNQ5, RERE and SUOX) as attractive drug targets that may warrant follow-up in functional studies and clinical trials, demonstrating the value of this study design for discovering drug targets in neuropsychiatric diseases.

Familial aggregation and shared genetic loading for major psychiatric disorders and type 2 diabetes.

AIMS/HYPOTHESIS: Psychiatric disorders, such as schizophrenia (SCZ), major depressive disorder (MDD) and bipolar disorder (BPD), are highly comorbid with type 2 diabetes. However, the mechanisms underlying such comorbidity are understudied. This study explored the familial aggregation of common psychiatric disorders and type 2 diabetes by testing family history association, and investigated the shared genetic loading between them by testing the polygenic risk score (PRS) association. METHODS: A total of 105,184 participants were recruited from the Taiwan Biobank, and genome-wide genotyping data were available for 95,238 participants. The Psychiatric Genomics Consortium-derived PRS for SCZ, MDD and BPD was calculated. Logistic regression was used to estimate the OR with CIs between a family history of SCZ/MDD/BPD and a family history of type 2 diabetes, and between the PRS and the risk of type 2 diabetes. RESULTS: A family history of type 2 diabetes was associated with a family history of SCZ (OR 1.23, 95% CI 1.08, 1.40), MDD (OR 1.19, 95% CI 1.13, 1.26) and BPD (OR 1.26, 95% CI 1.15, 1.39). Compared with paternal type 2 diabetes, maternal type 2 diabetes was associated with a higher risk of a family history of SCZ. SCZ PRS was negatively associated with type 2 diabetes in women (OR 0.92, 95% CI 0.88, 0.97), but not in men; the effect of SCZ PRS reduced after adjusting for BMI. MDD PRS was positively associated with type 2 diabetes (OR 1.04, 95% CI 1.00, 1.07); the effect of MDD PRS reduced after adjusting for BMI or smoking. BPD PRS was not associated with type 2 diabetes. CONCLUSIONS/INTERPRETATION: The comorbidity of type 2 diabetes with psychiatric disorders may be explained by shared familial factors. The shared polygenic loading between MDD and type 2 diabetes implies not only pleiotropy but also a shared genetic aetiology for the mechanism behind the comorbidity. The negative correlation between polygenic loading for SCZ and type 2 diabetes implies the role of environmental factors.

GWAS Identifies Risk Locus for Erectile Dysfunction and Implicates Hypothalamic Neurobiology and Diabetes in Etiology.

Erectile dysfunction (ED) is a common condition affecting more than 20% of men over 60 years, yet little is known about its genetic architecture. We performed a genome-wide association study of ED in 6,175 case subjects among 223,805 European men and identified one locus at 6q16.3 (lead variant rs57989773, OR 1.20 per C-allele; p = 5.71 × 10-14), located between MCHR2 and SIM1. In silico analysis suggests SIM1 to confer ED risk through hypothalamic dysregulation. Mendelian randomization provides evidence that genetic risk of type 2 diabetes mellitus is a cause of ED (OR 1.11 per 1-log unit higher risk of type 2 diabetes). These findings provide insights into the biological underpinnings and the causes of ED and may help prioritize the development of future therapies for this common disorder.

Causality of abdominal obesity on cognition: a trans-ethnic Mendelian randomization study.

BACKGROUND: Obesity has been associated with cognition in observational studies; however, whether its effect is confounding or a reverse causality remains inconclusive. This study aimed to investigate the causal relationships of overall obesity, measured by body mass index (BMI), and abdominal adiposity, measured by waist-hip ratio adjusted for BMI (WHRadjBMI), and cognition across European and Asian populations using Mendelian randomization (MR) analysis. METHODS: We used publicly available genome-wide association study (GWAS) summary data of European ancestry, including BMI (n = 322,154) and WHRadjBMI (n = 210,088) from the GIANT consortium, and cognition performance (n = 257,828) from the UK Biobank and COGENT consortium. Data for individuals of Asian ancestry were retrieved from Taiwan Biobank to perform GWAS for BMI (n = 65,689), WHRadjBMI (n = 65,683), and Mini-Mental State Examination (MMSE, n = 21,273). MR analysis was carried out using the inverse-variance weighted method for the main results. Further, we examined the overall pleiotropy by MR-Egger intercept, and detected and adjusted for possible outliers using MR PRESSO. RESULTS: No causal effect of BMI on cognition performance (beta [95% CI] = 0.00 [-0.07, 0.07], p value = 0.91) was found for Europeans; however, a 1-SD increase in WHRadjBMI was associated with a 0.07 standardized score decrease in cognition performance (beta [95% CI] = -0.07 [-0.12, -0.02], p value = 0.006). Further, no causal effect of BMI on MMSE (beta [95% CI] = 0.01 [-0.08, 0.10], p = 0.91) was found for Asians; however, a 1-SD increase in WHRadjBMI was associated with a 0.17 standardized score decrease in MMSE (beta [95% CI] = -0.17 [-0.30, -0.03], p = 0.02). In both populations, overall pleiotropy was not detected, and outliers did not affect the robustness of the main findings. CONCLUSIONS: This trans-ethnic MR study reveals that abdominal adiposity, as measured by WHR adjusted for BMI, impairs cognition, whereas weak evidence suggests that BMI impairs cognition.

Correction: Phenome-wide heritability analysis of the UK Biobank.

[This corrects the article DOI: 10.1371/journal.pgen.1006711.].

Large-scale meta-analysis highlights the hypothalamic-pituitary-gonadal axis in the genetic regulation of menstrual cycle length.

The normal menstrual cycle requires a delicate interplay between the hypothalamus, pituitary and ovary. Therefore, its length is an important indicator of female reproductive health. Menstrual cycle length has been shown to be partially controlled by genetic factors, especially in the follicle-stimulating hormone beta-subunit (FSHB) locus. A genome-wide association study meta-analysis of menstrual cycle length in 44 871 women of European ancestry confirmed the previously observed association with the FSHB locus and identified four additional novel signals in, or near, the GNRH1, PGR, NR5A2 and INS-IGF2 genes. These findings not only confirm the role of the hypothalamic-pituitary-gonadal axis in the genetic regulation of menstrual cycle length but also highlight potential novel local regulatory mechanisms, such as those mediated by IGF2.

Prospective study of polygenic risk, protective factors, and incident depression following combat deployment in US Army soldiers.

BACKGROUND: Whereas genetic susceptibility increases the risk for major depressive disorder (MDD), non-genetic protective factors may mitigate this risk. In a large-scale prospective study of US Army soldiers, we examined whether trait resilience and/or unit cohesion could protect against the onset of MDD following combat deployment, even in soldiers at high polygenic risk. METHODS: Data were analyzed from 3079 soldiers of European ancestry assessed before and after their deployment to Afghanistan. Incident MDD was defined as no MDD episode at pre-deployment, followed by a MDD episode following deployment. Polygenic risk scores were constructed from a large-scale genome-wide association study of major depression. We first examined the main effects of the MDD PRS and each protective factor on incident MDD. We then tested the effects of each protective factor on incident MDD across strata of polygenic risk. RESULTS: Polygenic risk showed a dose-response relationship to depression, such that soldiers at high polygenic risk had greatest odds for incident MDD. Both unit cohesion and trait resilience were prospectively associated with reduced risk for incident MDD. Notably, the protective effect of unit cohesion persisted even in soldiers at highest polygenic risk. CONCLUSIONS: Polygenic risk was associated with new-onset MDD in deployed soldiers. However, unit cohesion - an index of perceived support and morale - was protective against incident MDD even among those at highest genetic risk, and may represent a potent target for promoting resilience in vulnerable soldiers. Findings illustrate the value of combining genomic and environmental data in a prospective design to identify robust protective factors for mental health.

The Shared Genetic Basis of Educational Attainment and Cerebral Cortical Morphology.

Individual differences in educational attainment are linked to differences in intelligence, and predict important social, economic, and health outcomes. Previous studies have found common genetic factors that influence educational achievement, cognitive performance and total brain volume (i.e., brain size). Here, in a large sample of participants from the UK Biobank, we investigate the shared genetic basis between educational attainment and fine-grained cerebral cortical morphological features, and associate this genetic variation with a related aspect of cognitive ability. Importantly, we execute novel statistical methods that enable high-dimensional genetic correlation analysis, and compute high-resolution surface maps for the genetic correlations between educational attainment and vertex-wise morphological measurements. We conduct secondary analyses, using the UK Biobank verbal-numerical reasoning score, to confirm that variation in educational attainment that is genetically correlated with cortical morphology is related to differences in cognitive performance. Our analyses relate the genetic overlap between cognitive ability and cortical thickness measurements to bilateral primary motor cortex as well as predominantly left superior temporal cortex and proximal regions. These findings extend our understanding of the neurobiology that connects genetic variation to individual differences in educational attainment and cognitive performance.

Phenome-wide heritability analysis of the UK Biobank.

Heritability estimation provides important information about the relative contribution of genetic and environmental factors to phenotypic variation, and provides an upper bound for the utility of genetic risk prediction models. Recent technological and statistical advances have enabled the estimation of additive heritability attributable to common genetic variants (SNP heritability) across a broad phenotypic spectrum. Here, we present a computationally and memory efficient heritability estimation method that can handle large sample sizes, and report the SNP heritability for 551 complex traits derived from the interim data release (152,736 subjects) of the large-scale, population-based UK Biobank, comprising both quantitative phenotypes and disease codes. We demonstrate that common genetic variation contributes to a broad array of quantitative traits and human diseases in the UK population, and identify phenotypes whose heritability is moderated by age (e.g., a majority of physical measures including height and body mass index), sex (e.g., blood pressure related traits) and socioeconomic status (education). Our study represents the first comprehensive phenome-wide heritability analysis in the UK Biobank, and underscores the importance of considering population characteristics in interpreting heritability.

Association between GLP-1 receptor gene polymorphisms with reward learning, anhedonia and depression diagnosis.

BACKGROUND: Glucagon-like peptide-1 receptors (GLP-1Rs) are widely expressed in the brain. Evidence suggests that they may play a role in reward responses and neuroprotection. However, the association of GLP-1R with anhedonia and depression diagnosis has not been studied. Here, we examined the association of GLP-1R polymorphisms with objective and subjective measures of anhedonia, as well as depression diagnosis. METHODS: Objective [response bias assessed by the probabilistic reward task (PRT)] and subjective [Snaith-Hamilton Pleasure Scale (SHAPS)] measures of anhedonia, clinical variables and DNA samples were collected from 100 controls and 164 patients at McLean Hospital. An independent sample genotyped as part of the Psychiatric Genomics Consortium (PGC) was used to study the effect of putative GLP-1R polymorphisms linked to response bias in PRT on depression diagnosis. RESULTS: The C allele in rs1042044 was significantly associated with increased PRT response bias, when controlling for age, sex, case-control status and PRT discriminability. AA genotype of rs1042044 showed higher anhedonia phenotype based on SHAPS scores. However, analysis of PGC major depressive disorder data showed no association between rs1042044 and depression diagnosis. CONCLUSION: Findings suggest a possible association of rs1042044 with anhedonia but no association with depression diagnosis.

Genome-wide association study of shared liability to anxiety disorders in Army STARRS.

Anxiety disorders (ANX), namely generalized anxiety, panic disorder, and phobias, are common, etiologically complex syndromes that show increasing prevalence and comorbidity throughout adolescence and beyond. Few genome-wide association studies (GWAS) examining ANX risk have been published and almost exclusively in individuals of European ancestry. In this study, we phenotyped participants from the Army Study To Assess Risk and Resilience in Servicemembers (STARRS) to approximate DSM-based ANX diagnoses. We factor-analyzed those to create a single dimensional anxiety score for each subject. GWAS were conducted using that score within each of three ancestral groups (EUR, AFR, LAT) and then meta-analyzed across ancestries (NTotal = 16,510). We sought to (a) replicate prior ANX GWAS findings in ANGST; (b) determine whether results extended to other ancestry groups; and (c) meta-analyze with ANGST for increased power to identify novel susceptibility loci. No reliable genome-wide significant SNP associations were detected in STARRS. However, SNPs within the CAMKMT gene located in region 2p21 associated with shared ANX risk in ANGST were replicated in EUR soldiers but not other ancestry groups. Combining EUR STARRS and ANGST (N = 28,950) yielded a more robust 2p21 association signal (p = 9.08x10-11 ). Gene-based analyses supported three genes within 2p21 and LBX1 on chromosome 10. More powerful ANX genetic studies will be required to identify further loci.

Genome-wide analysis of insomnia disorder.

Insomnia is a worldwide problem with substantial deleterious health effects. Twin studies have shown a heritable basis for various sleep-related traits, including insomnia, but robust genetic risk variants have just recently begun to be identified. We conducted genome-wide association studies (GWAS) of soldiers in the Army Study To Assess Risk and Resilience in Servicemembers (STARRS). GWAS were carried out separately for each ancestral group (EUR, AFR, LAT) using logistic regression for each of the STARRS component studies (including 3,237 cases and 14,414 controls), and then meta-analysis was conducted across studies and ancestral groups. Heritability (SNP-based) for lifetime insomnia disorder was significant (h2g = 0.115, p = 1.78 × 10-4 in EUR). A meta-analysis including three ancestral groups and three study cohorts revealed a genome-wide significant locus on Chr 7 (q11.22) (top SNP rs186736700, OR = 0.607, p = 4.88 × 10-9) and a genome-wide significant gene-based association (p = 7.61 × 10-7) in EUR for RFX3 on Chr 9. Polygenic risk for sleeplessness/insomnia severity in UK Biobank was significantly positively associated with likelihood of insomnia disorder in STARRS. Genetic contributions to insomnia disorder in STARRS were significantly positively correlated with major depressive disorder (rg = 0.44, se = 0.22, p = 0.047) and type 2 diabetes (rg = 0.43, se = 0.20, p = 0.037), and negatively with morningness chronotype (rg = -0.34, se = 0.17, p = 0.039) and subjective well being (rg = -0.59, se = 0.23, p = 0.009) in external datasets. Insomnia associated loci may contribute to the genetic risk underlying a range of health conditions including psychiatric disorders and metabolic disease.