Search

As there are effective treatments to reduce hip fractures, identification of patients at high risk of hip fracture is important to inform efficient intervention strategies. To obtain a new tool for hip fracture prediction, we developed a protein-based risk score in the Cardiovascular Health Study using an aptamer-based proteomic platform. The proteomic risk score predicted incident hip fractures and improved hip fracture discrimination in two Trøndelag Health Study validation cohorts using the same aptamer-based platform. When transferred to an antibody-based proteomic platform in a UK Biobank validation cohort, the proteomic risk score was strongly associated with hip fractures (hazard ratio per s.d. increase, 1.64; 95% confidence interval 1.53-1.77). The proteomic risk score, but not available polygenic risk scores for fractures or bone mineral density, improved the C-index beyond the fracture risk assessment tool (FRAX), which integrates information from clinical risk factors (C-index, FRAX 0.735 versus FRAX + proteomic risk score 0.776). The developed proteomic risk score constitutes a new tool for stratifying patients according to hip fracture risk; however, its improvement in hip fracture discrimination is modest and its clinical utility beyond FRAX with information on femoral neck bone mineral density remains to be determined., (© 2024. The Author(s).)

Hip fractures are associated with significant disability, high cost, and mortality. However, the exact biological mechanisms underlying susceptibility to hip fractures remain incompletely understood. In an exploratory search of the underlying biology as reflected through the circulating proteome, we performed a comprehensive Circulating Proteome Association Study (CPAS) meta-analysis for incident hip fractures. Analyses included 6430 subjects from two prospective cohort studies (Cardiovascular Health Study and Trøndelag Health Study) with circulating proteomics data (aptamer-based 5 K SomaScan version 4.0 assay; 4979 aptamers). Associations between circulating protein levels and incident hip fractures were estimated for each cohort using age and sex-adjusted Cox regression models. Participants experienced 643 incident hip fractures. Compared with the individual studies, inverse-variance weighted meta-analyses yielded more statistically significant associations, identifying 23 aptamers associated with incident hip fractures (conservative Bonferroni correction 0.05/4979, P < 1.0 × 10-5). The aptamers most strongly associated with hip fracture risk corresponded to two proteins of the growth hormone/insulin growth factor system (GHR and IGFBP2), as well as GDF15 and EGFR. High levels of several inflammation-related proteins (CD14, CXCL12, MMP12, ITIH3) were also associated with increased hip fracture risk. Ingenuity pathway analysis identified reduced LXR/RXR activation and increased acute phase response signaling to be overrepresented among those proteins associated with increased hip fracture risk. These analyses identified several circulating proteins and pathways consistently associated with incident hip fractures. These findings underscore the usefulness of the meta-analytic approach for comprehensive CPAS in a similar manner as has previously been observed for large-scale human genetic studies. Future studies should investigate the underlying biology of these potential novel drug targets., (© The Author(s) 2024. Published by Oxford University Press on behalf of the American Society for Bone and Mineral Research.)

Background: Data on the proportion of venous thromboembolism (VTE) risk attributed to prothrombotic genotypes in men and women are limited., Objectives: We aimed to estimate the population attributable fraction (PAF) of VTE for recognized, common prothrombotic genotypes in men and women using a population-based case cohort., Methods: Cases with incident VTE ( n  = 1493) and a randomly sampled subcohort ( n  = 13,069) were derived from the Tromsø study (1994-2012) and the Trøndelag Health Study (1995-2008) cohorts. DNA samples were genotyped for 17 single-nucleotide polymorphisms (SNPs) previously associated with VTE. PAFs with 95% bias-corrected CIs (based on 10,000 bootstrap samples) were estimated for SNPs significantly associated with VTE, and a 6-SNP cumulative model was constructed for both sexes., Results: In women, the individual PAFs for SNPs included in the cumulative model were 16.9% for ABO (rs8176719), 17.6% for F11 (rs2036914), 15.1% for F11 (rs2289252), 8.7% for FVL (rs6025), 6.0% for FGG (rs2066865), and 0.2% for F2 (rs1799963). The cumulative PAF for this 6-SNP model was 37.8%. In men, the individual PAFs for SNPs included in the cumulative model were 21.3% for ABO , 12.2% for F11 (rs2036914), 10.4% for F11 (rs2289252), 7.5% for FVL , 7.8% for FGG , and 1.1% for F2 . This resulted in a cumulative PAF in men of 51.9%., Conclusion: Our findings in a Norwegian population suggest that 52% and 38% of the VTEs can be attributed to known prothrombotic genotypes in men and women, respectively., (© 2024 The Author(s).)

Osteoporotic fracture is among the most common and costly of diseases. While reasonably heritable, its genetic determinants have remained elusive. Forearm fractures are the most common clinically recognized osteoporotic fractures with a relatively high heritability. To establish an atlas of the genetic determinants of forearm fractures, we performed genome-wide association analyses including 100,026 forearm fracture cases. We identified 43 loci, including 26 new fracture loci. Although most fracture loci associated with bone mineral density, we also identified loci that primarily regulate bone quality parameters. Functional studies of one such locus, at TAC4, revealed that Tac4 -/- mice have reduced mechanical bone strength. The strongest forearm fracture signal, at WNT16, displayed remarkable bone-site-specificity with no association with hip fractures. Tall stature and low body mass index were identified as new causal risk factors for fractures. The insights from this atlas may improve fracture prediction and enable therapeutic development to prevent fractures., (© 2023. The Author(s).)

Objective: The objective of this study was to aggregate data for the first genomewide association study meta-analysis of cluster headache, to identify genetic risk variants, and gain biological insights., Methods: A total of 4,777 cases (3,348 men and 1,429 women) with clinically diagnosed cluster headache were recruited from 10 European and 1 East Asian cohorts. We first performed an inverse-variance genomewide association meta-analysis of 4,043 cases and 21,729 controls of European ancestry. In a secondary trans-ancestry meta-analysis, we included 734 cases and 9,846 controls of East Asian ancestry. Candidate causal genes were prioritized by 5 complementary methods: expression quantitative trait loci, transcriptome-wide association, fine-mapping of causal gene sets, genetically driven DNA methylation, and effects on protein structure. Gene set and tissue enrichment analyses, genetic correlation, genetic risk score analysis, and Mendelian randomization were part of the downstream analyses., Results: The estimated single nucleotide polymorphism (SNP)-based heritability of cluster headache was 14.5%. We identified 9 independent signals in 7 genomewide significant loci in the primary meta-analysis, and one additional locus in the trans-ethnic meta-analysis. Five of the loci were previously known. The 20 genes prioritized as potentially causal for cluster headache showed enrichment to artery and brain tissue. Cluster headache was genetically correlated with cigarette smoking, risk-taking behavior, attention deficit hyperactivity disorder (ADHD), depression, and musculoskeletal pain. Mendelian randomization analysis indicated a causal effect of cigarette smoking intensity on cluster headache. Three of the identified loci were shared with migraine., Interpretation: This first genomewide association study meta-analysis gives clues to the biological basis of cluster headache and indicates that smoking is a causal risk factor. ANN NEUROL 2023;94:713-726., (© 2023 The Authors. Annals of Neurology published by Wiley Periodicals LLC on behalf of American Neurological Association.)

Stroke is the second leading cause of death and disability worldwide. Stroke prevalence varies by sex and ancestry, possibly due to genetic heterogeneity between subgroups. We performed a genome-wide meta-analysis of 16 biobanks across multiple ancestries to study the genetics of ischemic stroke (60,176 cases, 1,310,725 controls) as part of the Global Biobank Meta-analysis Initiative (GBMI) and further combined the results with previously published MegaStroke. Five novel loci for ischemic stroke ( LAMC1 , CALCRL , PLSCR1 , CDKN1A , and SWAP70 ) were identified after replication in four additional datasets. One previously reported locus showed significant ancestry heterogeneity ( ABO ), and one showed significant sex heterogeneity ( ALDH2 ). The ALDH2 association was male specific (males p = 1.67e-24, females p = 0.126) and was additionally observed only in the East Asian ancestry (male) samples. These findings emphasize the need for more diverse datasets with large sample sizes to further understand the genetic predisposition of stroke in different ancestry and sex groups., Competing Interests: N.R.S. is an advisor for Abbott, Philips, and Shockwave and has received honoraria for speaking from Zoll, Cordis. C.J.W.’s spouse works at Regeneron Pharmaceuticals., (© 2023 The Authors.)

Appendicular lean mass (ALM) associates with mobility and bone mineral density (BMD). While associations between gut microbiota composition and ALM have been reported, previous studies rely on relatively small sample sizes. Here, we determine the associations between prevalent gut microbes and ALM in large discovery and replication cohorts with information on relevant confounders within the population-based Norwegian HUNT cohort (n = 5196, including women and men). We show that the presence of three bacterial species - Coprococcus comes, Dorea longicatena, and Eubacterium ventriosum - are reproducibly associated with higher ALM. When combined into an anabolic species count, participants with all three anabolic species have 0.80 kg higher ALM than those without any. In an exploratory analysis, the anabolic species count is positively associated with femoral neck and total hip BMD. We conclude that the anabolic species count may be used as a marker of ALM and BMD. The therapeutic potential of these anabolic species to prevent sarcopenia and osteoporosis needs to be determined., (© 2023. The Author(s).)

The timing of parturition is crucial for neonatal survival and infant health. Yet, its genetic basis remains largely unresolved. We present a maternal genome-wide meta-analysis of gestational duration (n = 195,555), identifying 22 associated loci (24 independent variants) and an enrichment in genes differentially expressed during labor. A meta-analysis of preterm delivery (18,797 cases, 260,246 controls) revealed six associated loci and large genetic similarities with gestational duration. Analysis of the parental transmitted and nontransmitted alleles (n = 136,833) shows that 15 of the gestational duration genetic variants act through the maternal genome, whereas 7 act both through the maternal and fetal genomes and 2 act only via the fetal genome. Finally, the maternal effects on gestational duration show signs of antagonistic pleiotropy with the fetal effects on birth weight: maternal alleles that increase gestational duration have negative fetal effects on birth weight. The present study provides insights into the genetic effects on the timing of parturition and the complex maternal-fetal relationship between gestational duration and birth weight., (© 2023. The Author(s).)

Tobacco and alcohol use are heritable behaviours associated with 15% and 5.3% of worldwide deaths, respectively, due largely to broad increased risk for disease and injury 1-4 . These substances are used across the globe, yet genome-wide association studies have focused largely on individuals of European ancestries 5 . Here we leveraged global genetic diversity across 3.4 million individuals from four major clines of global ancestry (approximately 21% non-European) to power the discovery and fine-mapping of genomic loci associated with tobacco and alcohol use, to inform function of these loci via ancestry-aware transcriptome-wide association studies, and to evaluate the genetic architecture and predictive power of polygenic risk within and across populations. We found that increases in sample size and genetic diversity improved locus identification and fine-mapping resolution, and that a large majority of the 3,823 associated variants (from 2,143 loci) showed consistent effect sizes across ancestry dimensions. However, polygenic risk scores developed in one ancestry performed poorly in others, highlighting the continued need to increase sample sizes of diverse ancestries to realize any potential benefit of polygenic prediction., (© 2022. The Author(s).)

Hip fracture is the clinically most important fracture, but the genetic architecture of hip fracture is unclear. Here, we perform a large-scale hip fracture genome-wide association study meta-analysis and Mendelian randomization study using five cohorts from European biobanks. The results show that five genetic signals associate with hip fractures. Among these, one signal associates with falls, but not with bone mineral density (BMD), while four signals are in loci known to be involved in bone biology. Mendelian randomization analyses demonstrate a strong causal effect of decreased femoral neck BMD and moderate causal effects of Alzheimer's disease and having ever smoked regularly on risk of hip fractures. The substantial causal effect of decreased femoral neck BMD on hip fractures in both young and old subjects and in both men and women supports the use of change in femoral neck BMD as a surrogate outcome for hip fractures in clinical trials., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 The Author(s). Published by Elsevier Inc. All rights reserved.)

The Trøndelag Health Study (HUNT) is a population-based cohort of ∼229,000 individuals recruited in four waves beginning in 1984 in Trøndelag County, Norway. Approximately 88,000 of these individuals have available genetic data from array genotyping. HUNT participants were recruited during four community-based recruitment waves and provided information on health-related behaviors, self-reported diagnoses, family history of disease, and underwent physical examinations. Linkage via the Norwegian personal identification number integrates digitized health care information from doctor visits and national health registries including death, cancer and prescription registries. Genome-wide association studies of HUNT participants have provided insights into the mechanism of cardiovascular, metabolic, osteoporotic, and liver-related diseases, among others. Unique features of this cohort that facilitate research include nearly 40 years of longitudinal follow-up in a motivated and well-educated population, family data, comprehensive phenotyping, and broad availability of DNA, RNA, urine, fecal, plasma, and serum samples., Competing Interests: G.R.A. works for Regeneron Pharmaceuticals. C.J.W.’s spouse works for Regeneron Pharmaceuticals., (© 2022 The Author(s).)

Background:  The proportion of venous thromboembolism (VTE) events that can be attributed to established prothrombotic genotypes has been scarcely investigated in the general population. We aimed to estimate the proportion of VTEs in the population that could be attributed to established prothrombotic genotypes using a population-based case-cohort., Methods:  Cases with incident VTE ( n  = 1,493) and a randomly sampled subcohort ( n  = 13,069) were derived from the Tromsø Study (1994-2012) and the Nord-Trøndelag Health (HUNT) study (1995-2008). DNA samples were genotyped for 17 single-nucleotide polymorphisms (SNPs) associated with VTE. Hazard ratios with 95% confidence intervals (CIs) were estimated in Cox regression models. Population-attributable fractions (PAFs) with 95% bias-corrected CIs (based on 10,000 bootstrap samples) were estimated using a cumulative model where SNPs significantly associated with VTE were added one by one in ranked order of the individual PAFs., Results:  Six SNPs were significantly associated with VTE (rs1799963 [Prothrombin], rs2066865 [FGG], rs6025 [FV Leiden], rs2289252 [F11], rs2036914 [F11], and rs8176719 [ABO]). The cumulative PAF for the six-SNP model was 45.3% (95% CI: 19.7-71.6) for total VTE and 61.7% (95% CI: 19.6-89.3) for unprovoked VTE. The PAF for prothrombotic genotypes was higher for deep vein thrombosis (DVT; 52.9%) than for PE (33.8%), and higher for those aged <70 years (66.1%) than for those aged ≥70 years (24.9%)., Conclusion:  Our findings suggest that 45 to 62% of all VTE events in the population can be attributed to known prothrombotic genotypes. The PAF of established prothrombotic genotypes was higher in DVT than in PE, and higher in the young than in the elderly., Competing Interests: None declared., (Thieme. All rights reserved.)

Polygenic risk scores (PRS) measure genetic disease susceptibility by combining risk effects across the genome. For coronary artery disease (CAD), type 2 diabetes (T2D), and breast and prostate cancer, we performed cross-ancestry evaluation of genome-wide PRSs in six biobanks in Europe, the United States, and Asia. We studied transferability of these highly polygenic, genome-wide PRSs across global ancestries, within European populations with different health-care systems, and local population substructures in a population isolate. All four PRSs had similar accuracy across European and Asian populations, with poorer transferability in the smaller group of individuals of African ancestry. The PRSs had highly similar effect sizes in different populations of European ancestry, and in early- and late-settlement regions with different recent population bottlenecks in Finland. Comparing genome-wide PRSs to PRSs containing a smaller number of variants, the highly polygenic, genome-wide PRSs generally displayed higher effect sizes and better transferability across global ancestries. Our findings indicate that in the populations investigated, the current genome-wide polygenic scores for common diseases have potential for clinical utility within different health-care settings for individuals of European ancestry, but that the utility in individuals of African ancestry is currently much lower., Competing Interests: A.P. is a member of the Pfizer Genetics Scientific Advisory Panel. J.W.S is an unpaid member of the Bipolar/Depression Research Community Advisory Panel of 23andMe, a member of the Leon Levy Foundation Neuroscience Advisory Board, and received an honorarium for an internal seminar at Biogen, Inc. He is principal investigator of a collaborative study of the genetics of depression and bipolar disorder sponsored by 23andMe for which 23andMe provides analysis time as in-kind support but no payments. B.M.N. is a member of the Deep Genomics Scientific Advisory Board and serves as a consultant for the Camp4 Therapeutics Corporation, Takeda Pharmaceutical, and Biogen. The remaining authors declare no conflict of interests., (© 2022 The Author(s).)

Background: The impact of the combination of obesity and multiple prothrombotic genotypes on venous thromboembolism (VTE) risk remains unclear., Objective: To investigate the joint effect of obesity and a genetic risk score (GRS) composed of established prothrombotic single nucleotide polymorphisms (SNPs) on VTE risk using a population-based case-cohort., Methods: Cases with incident VTE ( n  = 1,470) and a subcohort ( n  = 12,826) were derived from the Tromsø Study (1994-2012) and the Trøndelag Health Study (HUNT) (1995-2008). Participants were genotyped for ABO (rs8176719), F5 (rs6025), F2 (rs1799963), FGG (rs2066865), and F11 (rs2036914) SNPs. Age- and sex-adjusted hazard ratios (HRs) were estimated according to body mass index (BMI) categories and number of risk alleles for individual SNPs and the GRS (0-1, 2, 3, ≥4 alleles)., Results: The combination of obesity (BMI ≥ 30kg/m 2 ) and risk alleles, either as individual SNPs or as a GRS, had an additive effect on VTE risk (i.e., no biological interaction). Obese subjects who were carriers of ≥4 risk alleles had a 2.85-fold (95% confidence interval [CI]: 2.05-3.96) increased risk of overall VTE compared with those with BMI <25 kg/m 2 and 0 to 1 risk allele. However, in subgroups, the combination of obesity and ≥4 risk alleles was more pronounced for deep vein thrombosis (DVT) (HR: 3.20; 95% CI: 2.09-4.90) and unprovoked VTE (HR: 3.82; 95% CI: 2.25-6.47), suggesting a supra-additive effect., Conclusion: Our findings indicate that the combination of obesity and GRS has an additive effect on the risk of overall VTE. However, it may have a supra-additive effect on the risk of DVT and unprovoked VTE., Competing Interests: None declared., (Thieme. All rights reserved.)

Background: About one third of patients with chronic polyneuropathy have no obvious underlying etiology and are classified as having idiopathic polyneuropathy. The lack of knowledge about pathomechanisms and predisposing factors limits the development of effective prevention and treatment for these patients. We report the first genome-wide association study (GWAS) of idiopathic polyneuropathy. Methods: Cases with idiopathic polyneuropathy and healthy controls were identified by linkage to hospital records. We performed genome-wide association studies using genetic data from two large population-based health studies, the Trøndelag Health Study (HUNT, 1,147 cases and 62,204 controls) and UK Biobank (UKB, 946 cases and 383,052 controls). In a two-stage analysis design, we first treated HUNT as a discovery cohort and UK Biobank as a replication cohort. Secondly, we combined the two studies in a meta-analysis. Downstream analyses included genetic correlation to other traits and diseases. We specifically examined previously reported risk loci, and genes known to cause hereditary polyneuropathy. Results: No replicable risk loci were identified in the discovery-replication stage, in line with the limited predicted power of this approach. When combined in a meta-analysis, two independent loci reached statistical significance (rs7294354 in B4GALNT3, P -value 4.51 × 10 -8 ) and (rs147738081 near NR5A2, P -value 4.75 × 10 -8 ). Idiopathic polyneuropathy genetically correlated with several anthropometric measures, most pronounced for height, and with several pain-related traits. Conclusions: In this first GWAS of idiopathic polyneuropathy we identify two risk-loci that indicate possible pathogenetic mechanisms. Future collaborative efforts are needed to replicate and expand on these findings., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Winsvold, Kitsos, Thomas, Skogholt, Gabrielsen, Zwart and Nilsen.)

Background: The causal nature of the observed associations between serum lipids and apolipoproteins and kidney function are unclear., Methods: Using two-sample and multivariable Mendelian randomization (MR), we examined the causal effects of serum lipids and apolipoproteins on kidney function, indicated by the glomerular-filtration rate estimated using creatinine (eGFRcrea) or cystatin C (eGFRcys) and the urinary albumin-to-creatinine ratio (UACR). We obtained lipid- and apolipoprotein-associated genetic variants from the Global Lipids Genetics Consortium (n = 331 368) and UK Biobank (n = 441 016), respectively, and kidney-function markers from the Trøndelag Health Study (HUNT; n = 69 736) and UK Biobank (n = 464 207). The reverse causal direction was examined using variants associated with kidney-function markers selected from recent genome-wide association studies., Results: There were no strong associations between genetically predicted lipid and apolipoprotein levels with kidney-function markers. Some, but inconsistent, evidence suggested a weak association of higher genetically predicted atherogenic lipid levels [indicated by low-density lipoprotein cholesterol (LDL-C), triglycerides and apolipoprotein B] with increased eGFR and UACR. For high-density lipoprotein cholesterol (HDL-C), results differed between eGFRcrea and eGFRcys, but neither analysis suggested substantial effects. We found no clear evidence of a reverse causal effect of eGFR on lipid or apolipoprotein traits, but higher UACR was associated with higher LDL-C, triglyceride and apolipoprotein B levels., Conclusion: Our MR estimates suggest that serum lipid and apolipoprotein levels do not cause substantial changes in kidney function. A possible weak effect of higher atherogenic lipids on increased eGFR and UACR warrants further investigation. Processes leading to higher UACR may lead to more atherogenic lipid levels., (© The Author(s) 2021. Published by Oxford University Press on behalf of the International Epidemiological Association.)

The rs58542926C >T (E167K) variant of the transmembrane 6 superfamily member 2 gene ( TM6SF2 ) is associated with increased risks for nonalcoholic fatty liver disease (NAFLD) and type 2 diabetes (T2D). Nevertheless, the role of the TM6SF2 rs58542926 variant in glucose metabolism is poorly understood. We performed a sex-stratified analysis of the association between the rs58542926C >T variant and T2D in multiple cohorts. The E167K variant was significantly associated with T2D, especially in males. Using an E167K knockin (KI) mouse model, we found that male but not the female KI mice exhibited impaired glucose tolerance. As an ER membrane protein, TM6SF2 was found to interact with inositol-requiring enzyme 1 α (IRE1α), a primary ER stress sensor. The male Tm6sf2 KI mice exhibited impaired IRE1α signaling in the liver. In conclusion, the E167K variant of TM6SF2 is associated with glucose intolerance primarily in males, both in humans and mice., Competing Interests: As of June 2020, A.M. is an employee of Genentech and a holder of Roche stock. The spouse of C.J.W. works at Regeneron Pharmaceuticals. J.B.N. works at Regeneron Pharmaceuticals., (© 2021 The Authors.)

Background: Studies have reported that the combination of some prothrombotic genotypes and overt cancer yields a synergistic effect on VTE risk. Whether individual prothrombotic genotypes or number of risk alleles in a genetic risk score (GRS) affect VTE risk in occult cancer have not been addressed. The aim of this study was to investigate the joint effect of five prothrombotic genotypes and occult cancer on VTE risk., Methods: Cases with incident VTE (n = 1566) and a subcohort (n = 14,537) were sampled from the Scandinavian Thrombosis and Cancer Cohort (1993-2012). Five single nucleotide polymorphisms previously reported in a GRS were genotyped: ABO (rs8176719), F5 (rs6025), F2 (rs1799963), FGG (rs2066865) and F11 (rs2036914). Hazard ratios (HRs) for VTE by individual SNPs and GRS were estimated according to non-cancer and occult cancer (one year preceding a cancer diagnosis) exposure., Results: Occult cancer occurred in 1817 subjects, and of these, 93 experienced a VTE. The VTE risk was 4-fold higher (HR 4.05, 95% CI 3.28-5.00) in subjects with occult cancer compared with those without cancer. Among subjects with occult cancer, those with VTE had a higher proportion of prothrombotic and advanced cancers than those without VTE. The VTE risk increased according to individual prothrombotic genotypes and GRS in cancer-free subjects, while no such effect was observed in subjects with occult cancer (HR for ≥4 versus ≤1 risk alleles in GRS: 1.14, 95% CI 0.61-2.11)., Conclusions: Five well-established prothrombotic genotypes, individually or combined, were not associated with increased risk of VTE in individuals with occult cancer., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

Late-onset Alzheimer's disease is a prevalent age-related polygenic disease that accounts for 50-70% of dementia cases. Currently, only a fraction of the genetic variants underlying Alzheimer's disease have been identified. Here we show that increased sample sizes allowed identification of seven previously unidentified genetic loci contributing to Alzheimer's disease. This study highlights microglia, immune cells and protein catabolism as relevant to late-onset Alzheimer's disease, while identifying and prioritizing previously unidentified genes of potential interest. We anticipate that these results can be included in larger meta-analyses of Alzheimer's disease to identify further genetic variants that contribute to Alzheimer's pathology., (© 2021. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Objective: Haemorrhoidal disease (HEM) affects a large and silently suffering fraction of the population but its aetiology, including suspected genetic predisposition, is poorly understood. We report the first genome-wide association study (GWAS) meta-analysis to identify genetic risk factors for HEM to date., Design: We conducted a GWAS meta-analysis of 218 920 patients with HEM and 725 213 controls of European ancestry. Using GWAS summary statistics, we performed multiple genetic correlation analyses between HEM and other traits as well as calculated HEM polygenic risk scores (PRS) and evaluated their translational potential in independent datasets. Using functional annotation of GWAS results, we identified HEM candidate genes, which differential expression and coexpression in HEM tissues were evaluated employing RNA-seq analyses. The localisation of expressed proteins at selected loci was investigated by immunohistochemistry., Results: We demonstrate modest heritability and genetic correlation of HEM with several other diseases from the GI, neuroaffective and cardiovascular domains. HEM PRS validated in 180 435 individuals from independent datasets allowed the identification of those at risk and correlated with younger age of onset and recurrent surgery. We identified 102 independent HEM risk loci harbouring genes whose expression is enriched in blood vessels and GI tissues, and in pathways associated with smooth muscles, epithelial and endothelial development and morphogenesis. Network transcriptomic analyses highlighted HEM gene coexpression modules that are relevant to the development and integrity of the musculoskeletal and epidermal systems, and the organisation of the extracellular matrix., Conclusion: HEM has a genetic component that predisposes to smooth muscle, epithelial and connective tissue dysfunction., Competing Interests: Competing interests: Vladimir Vacic and Olga V. Sazonova are/were employed by and hold stock or stock options in 23andMe, Inc. All other authors have nothing to declare., (© Author(s) (or their employer(s)) 2021. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.)

Pharmaceutical drugs targeting dyslipidemia and cardiovascular disease (CVD) may increase the risk of fatty liver disease and other metabolic disorders. To identify potential novel CVD drug targets without these adverse effects, we perform genome-wide analyses of participants in the HUNT Study in Norway (n = 69,479) to search for protein-altering variants with beneficial impact on quantitative blood traits related to cardiovascular disease, but without detrimental impact on liver function. We identify 76 (11 previously unreported) presumed causal protein-altering variants associated with one or more CVD- or liver-related blood traits. Nine of the variants are predicted to result in loss-of-function of the protein. This includes ZNF529:p.K405X, which is associated with decreased low-density-lipoprotein (LDL) cholesterol (P = 1.3 × 10 -8 ) without being associated with liver enzymes or non-fasting blood glucose. Silencing of ZNF529 in human hepatoma cells results in upregulation of LDL receptor and increased LDL uptake in the cells. This suggests that inhibition of ZNF529 or its gene product should be prioritized as a novel candidate drug target for treating dyslipidemia and associated CVD.

Background: The role of combined prothrombotic genotypes in cancer-related venous thromboembolism (VTE) is scarcely studied. We aimed to investigate the impact of a 5-single nucleotide polymorphism (SNP) score on the risk of VTE in patients with and without cancer using a population-based case-cohort., Methods: Cases with a first VTE (n = 1493) and a subcohort (n = 13 072) were derived from the Tromsø Study (1994-2012) and the Nord-Trøndelag Health Study (1995-2008). Five SNPs previously reported as a risk score were genotyped: ABO (rs8176719), F5 (rs6025), F2 (rs1799963), FGG (rs2066865), and F11 (rs2036914). Hazard ratios (HRs) for VTE were estimated according to cancer status and the number of risk alleles in the 5-SNP score (0-1, 2-3, and ≥4 alleles)., Results: During a median follow-up of 12.3 years, 1496 individuals were diagnosed with cancer, of whom 232 experienced VTE. The VTE risk increased with the number of risk alleles in the 5-SNP score among subjects without and with cancer. In cancer-free subjects, the HR was 2.17 (95% confidence interval [CI] 1.79-2.62) for ≥4 versus 0-1 risk alleles. In cancer patients, the corresponding HR was 1.93 (95% CI 1.28-2.91). The combination of cancer and ≥4 risk alleles yielded a 17-fold (HR 17.1, 95% CI 12.5-23.4) higher risk of VTE compared with cancer-free subjects with 0-1 risk alleles., Conclusion: The risk of VTE increases with the number of prothrombotic risk alleles in subjects with and without cancer, and the combination of prothrombotic risk alleles and cancer leads to a highly elevated risk of VTE., (© 2020 The Authors. Journal of Thrombosis and Haemostasis published by Wiley Periodicals LLC on behalf of International Society on Thrombosis and Haemostasis.)

Background: Abdominal aortic aneurysm (AAA) is an important cause of cardiovascular mortality; however, its genetic determinants remain incompletely defined. In total, 10 previously identified risk loci explain a small fraction of AAA heritability., Methods: We performed a genome-wide association study in the Million Veteran Program testing ≈18 million DNA sequence variants with AAA (7642 cases and 172 172 controls) in veterans of European ancestry with independent replication in up to 4972 cases and 99 858 controls. We then used mendelian randomization to examine the causal effects of blood pressure on AAA. We examined the association of AAA risk variants with aneurysms in the lower extremity, cerebral, and iliac arterial beds, and derived a genome-wide polygenic risk score (PRS) to identify a subset of the population at greater risk for disease., Results: Through a genome-wide association study, we identified 14 novel loci, bringing the total number of known significant AAA loci to 24. In our mendelian randomization analysis, we demonstrate that a genetic increase of 10 mm Hg in diastolic blood pressure (odds ratio, 1.43 [95% CI, 1.24-1.66]; P =1.6×10 -6 ), as opposed to systolic blood pressure (odds ratio, 1.06 [95% CI, 0.97-1.15]; P =0.2), likely has a causal relationship with AAA development. We observed that 19 of 24 AAA risk variants associate with aneurysms in at least 1 other vascular territory. A 29-variant PRS was strongly associated with AAA (odds ratio PRS , 1.26 [95% CI, 1.18-1.36]; P PRS =2.7×10 -11 per SD increase in PRS), independent of family history and smoking risk factors (odds ratio PRS+family history+smoking , 1.24 [95% CI, 1.14-1.35]; P PRS =1.27×10 -6 ). Using this PRS, we identified a subset of the population with AAA prevalence greater than that observed in screening trials informing current guidelines., Conclusions: We identify novel AAA genetic associations with therapeutic implications and identify a subset of the population at significantly increased genetic risk of AAA independent of family history. Our data suggest that extending current screening guidelines to include testing to identify those with high polygenic AAA risk, once the cost of genotyping becomes comparable with that of screening ultrasound, would significantly increase the yield of current screening at reasonable cost.

A major challenge in genetic association studies is that most associated variants fall in the non-coding part of the human genome. We searched for variants associated with bone mineral density (BMD) after enriching the discovery cohort for loss-of-function (LoF) mutations by sequencing a subset of the Nord-Trøndelag Health Study, followed by imputation in the remaining sample (N = 19,705), and identified ten known BMD loci. However, one previously unreported variant, LoF mutation in MEPE, p.(Lys70IlefsTer26, minor allele frequency [MAF] = 0.8%), was associated with decreased ultradistal forearm BMD (P-value = 2.1 × 10 -18 ), and increased osteoporosis (P-value = 4.2 × 10 -5 ) and fracture risk (P-value = 1.6 × 10 -5 ). The MEPE LoF association with BMD and fractures was further evaluated in 279,435 UK (MAF = 0.05%, heel bone estimated BMD P-value = 1.2 × 10 -16 , any fracture P-value = 0.05) and 375,984 Icelandic samples (MAF = 0.03%, arm BMD P-value = 0.12, forearm fracture P-value = 0.005). Screening for the MEPE LoF mutations before adulthood could potentially prevent osteoporosis and fractures due to the lifelong effect on BMD observed in the study. A key implication for precision medicine is that high-impact functional variants missing from the publicly available cosmopolitan panels could be clinically more relevant than polygenic risk scores.

Attention-deficit/hyperactivity disorder (ADHD) is a common neurodevelopmental disorder with a complex genetic background, hampering identification of underlying genetic risk factors. We hypothesized that combining linkage analysis and whole-exome sequencing (WES) in multi-generation pedigrees with multiple affected individuals can point toward novel ADHD genes. Three families with multiple ADHD-affected members (N total  = 70) and apparent dominant inheritance pattern were included in this study. Genotyping was performed in 37 family members, and WES was additionally carried out in 10 of those. Linkage analysis was performed using multi-point analysis in Superlink Online SNP 1.1. From prioritized linkage regions with a LOD score ≥ 2, a total of 24 genes harboring rare variants were selected. Those genes were taken forward and were jointly analyzed in gene-set analyses of exome-chip data using the MAGMA software in an independent sample of patients with persistent ADHD and healthy controls (N = 9365). The gene-set including all 24 genes together, and particularly the gene-set from one of the three families (12 genes), were significantly associated with persistent ADHD in this sample. Among the latter, gene-wide analysis for the AAED1 gene reached significance. A rare variant (rs151326868) within AAED1 segregated with ADHD in one of the families. The analytic strategy followed here is an effective approach for identifying novel ADHD risk genes. Additionally, this study suggests that both rare and more frequent variants in multiple genes act together in contributing to ADHD risk, even in individual multi-case families.

Thyroid stimulating hormone (TSH) is critical for normal development and metabolism. To better understand the genetic contribution to TSH levels, we conduct a GWAS meta-analysis at 22.4 million genetic markers in up to 119,715 individuals and identify 74 genome-wide significant loci for TSH, of which 28 are previously unreported. Functional experiments show that the thyroglobulin protein-altering variants P118L and G67S impact thyroglobulin secretion. Phenome-wide association analysis in the UK Biobank demonstrates the pleiotropic effects of TSH-associated variants and a polygenic score for higher TSH levels is associated with a reduced risk of thyroid cancer in the UK Biobank and three other independent studies. Two-sample Mendelian randomization using TSH index variants as instrumental variables suggests a protective effect of higher TSH levels (indicating lower thyroid function) on risk of thyroid cancer and goiter. Our findings highlight the pleiotropic effects of TSH-associated variants on thyroid function and growth of malignant and benign thyroid tumors.

Venous thromboembolism (VTE) is a frequent complication in patients with cancer. Homozygous carriers of the fibrinogen gamma gene (FGG) rs2066865 have a moderately increased risk of VTE, but the effect of the FGG variant in cancer is unknown. We aimed to investigate the effect of the FGG variant and active cancer on the risk of VTE. Cases with incident VTE (n=640) and a randomly selected age-weighted sub-cohort (n=3,734) were derived from a population-based cohort (the Tromsø study). Cox-regression was used to estimate hazard ratios (HR) with 95% confidence intervals (CI) for VTE according to categories of cancer and FGG In those without cancer, homozygosity at the FGG variant was associated with a 70% (HR 1.7, 95% CI: 1.2-2.3) increased risk of VTE compared to non-carriers. Cancer patients homozygous for the FGG variant had a two-fold (HR 2.0, 95% CI: 1.1-3.6) higher risk of VTE than cancer patients without the variant. Moreover, the six-months cumulative incidence of VTE among cancer patients was 6.4% (95% CI: 3.5-11.6) in homozygous carriers of FGG and 3.1% (95% CI: 2.3-4.7) in those without risk alleles. A synergistic effect was observed between rs2066865 and active cancer on the risk of VTE (synergy index: 1.81, 95% CI: 1.02-3.21, attributable proportion: 0.43, 95% CI: 0.11-0.74). In conclusion, homozygosity at the FGG variant and active cancer yielded a synergistic effect on the risk of VTE., (Copyright© 2020 Ferrata Storti Foundation.)

Objective: To develop and externally validate prediction models for incident hand osteoarthritis (OA) in a large population-based cohort of middle aged and older men and women., Design: We included 17,153 men and 18,682 women from a population-based cohort, aged 35-70 years at baseline (1995-1997). Incident hand OA were obtained from diagnostic codes in the Norwegian National Patient Register (1995-2018). We studied whether a range of self-reported and clinically measured predictors could predict hand OA, using the Area Under the receiver-operating Curve (AUC) from logistic regression. External validation of an existing prediction model for male hand OA was tested on discrimination in a sample of men. Bootstrapping was used to avoid overfitting., Results: The model for men showed modest discriminatory ability (AUC = 0.67, 95% CI 0.62-0.71). Adding a genetic risk score did not improve prediction. Similar discrimination was observed in the model for women (AUC = 0.62, 95% CI 0.59-0.64). Prediction was not improved by adding a genetic risk score or hormonal and reproductive factors. Applying external validation, similar results were observed among men in HUNT (The Nord-Trøndelag Health Study) as in the developmental sample (AUC = 0.62, 95% CI 0.57-0.65)., Conclusion: We developed prediction models for incident hand OA in men and women. For women, the model included body mass index (BMI), heavy physical work, high physical activity and perceived poor health. The model showed moderate discrimination. For men, we have shown that a prediction model including BMI, education and information on sleep can predict incident hand OA in several populations with moderate discriminative ability., (Copyright © 2020 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Risk factors that contribute to inter-individual differences in the age-of-onset of allergic diseases are poorly understood. The aim of this study was to identify genetic risk variants associated with the age at which symptoms of allergic disease first develop, considering information from asthma, hay fever and eczema. Self-reported age-of-onset information was available for 117,130 genotyped individuals of European ancestry from the UK Biobank study. For each individual, we identified the earliest age at which asthma, hay fever and/or eczema was first diagnosed and performed a genome-wide association study (GWAS) of this combined age-of-onset phenotype. We identified 50 variants with a significant independent association (P<3x10-8) with age-of-onset. Forty-five variants had comparable effects on the onset of the three individual diseases and 38 were also associated with allergic disease case-control status in an independent study (n = 222,484). We observed a strong negative genetic correlation between age-of-onset and case-control status of allergic disease (rg = -0.63, P = 4.5x10-61), indicating that cases with early disease onset have a greater burden of allergy risk alleles than those with late disease onset. Subsequently, a multivariate GWAS of age-of-onset and case-control status identified a further 26 associations that were missed by the univariate analyses of age-of-onset or case-control status only. Collectively, of the 76 variants identified, 18 represent novel associations for allergic disease. We identified 81 likely target genes of the 76 associated variants based on information from expression quantitative trait loci (eQTL) and non-synonymous variants, of which we highlight ADAM15, FOSL2, TRIM8, BMPR2, CD200R1, PRKCQ, NOD2, SMAD4, ABCA7 and UBE2L3. Our results support the notion that early and late onset allergic disease have partly distinct genetic architectures, potentially explaining known differences in pathophysiology between individuals., Competing Interests: C. Tian and D. A. Hinds both report support from 23andMe during the conduct of the study. C. Almqvist received grant 2017-00641 from Swedish Research Council and Swedish Initiative for Research on Microdata in the Social And Medical Sciences (SIMSAM) framework grant (340-2013-5867) for this work. G. H. Koppelman's institution received grants from the Lung Foundation of The Netherlands, the Ubbo Emmius Foundation, TEVA (The Netherlands), GlaxoSmithKline, Vertex, and the Tetri Foundation for other works. L. Paternoster received grant MR/J012165/1 from the UK Medical Research Council for this work and personal fees from Merck. The rest of the authors declare that they have no relevant conflicts of interest.

With very large sample sizes, biobanks provide an exciting opportunity to identify genetic components of complex traits. To analyze rare variants, region-based multiple-variant aggregate tests are commonly used to increase power for association tests. However, because of the substantial computational cost, existing region-based tests cannot analyze hundreds of thousands of samples while accounting for confounders such as population stratification and sample relatedness. Here we propose a scalable generalized mixed-model region-based association test, SAIGE-GENE, that is applicable to exome-wide and genome-wide region-based analysis for hundreds of thousands of samples and can account for unbalanced case-control ratios for binary traits. Through extensive simulation studies and analysis of the HUNT study with 69,716 Norwegian samples and the UK Biobank data with 408,910 White British samples, we show that SAIGE-GENE can efficiently analyze large-sample data (N > 400,000) with type I error rates well controlled.

Background: Variation in mitochondrial DNA (mtDNA) has been indicated in migraine pathogenesis, but genetic studies to date have focused on candidate variants, with sparse findings. We aimed to perform the first mitochondrial genome-wide association study of migraine, examining both single variants and mitochondrial haplogroups., Methods: In total, 71,860 participants from the population-based Nord-Trøndelag Health Study were genotyped. We excluded samples not passing quality control for nuclear genotypes, in addition to samples with low call rate and closely maternally related. We analysed 775 mitochondrial DNA variants in 4021 migraine cases and 14,288 headache-free controls, using logistic regression. In addition, we analysed 3831 cases and 13,584 controls who could be reliably assigned to a mitochondrial haplogroup. Lastly, we attempted to replicate previously reported mitochondrial DNA candidate variants., Results: Neither of the mitochondrial variants or haplogroups were associated with migraine. In addition, none of the previously reported mtDNA candidate variants replicated in our data., Conclusions: Our findings do not support a major role of mitochondrial genetic variation in migraine pathophysiology, but a larger sample is needed to detect rare variants and future studies should also examine heteroplasmic variation, epigenetic changes and copy-number variation.

Venous thromboembolism (VTE) is a significant contributor to morbidity and mortality. To advance our understanding of the biology contributing to VTE, we conducted a genome-wide association study (GWAS) of VTE and a transcriptome-wide association study (TWAS) based on imputed gene expression from whole blood and liver. We meta-analyzed GWAS data from 18 studies for 30 234 VTE cases and 172 122 controls and assessed the association between 12 923 718 genetic variants and VTE. We generated variant prediction scores of gene expression from whole blood and liver tissue and assessed them for association with VTE. Mendelian randomization analyses were conducted for traits genetically associated with novel VTE loci. We identified 34 independent genetic signals for VTE risk from GWAS meta-analysis, of which 14 are newly reported associations. This included 11 newly associated genetic loci (C1orf198, PLEK, OSMR-AS1, NUGGC/SCARA5, GRK5, MPHOSPH9, ARID4A, PLCG2, SMG6, EIF5A, and STX10) of which 6 replicated, and 3 new independent signals in 3 known genes. Further, TWAS identified 5 additional genetic loci with imputed gene expression levels differing between cases and controls in whole blood (SH2B3, SPSB1, RP11-747H7.3, RP4-737E23.2) and in liver (ERAP1). At some GWAS loci, we found suggestive evidence that the VTE association signal for novel and previously known regions colocalized with expression quantitative trait locus signals. Mendelian randomization analyses suggested that blood traits may contribute to the underlying risk of VTE. To conclude, we identified 16 novel susceptibility loci for VTE; for some loci, the association signals are likely mediated through gene expression of nearby genes.

Background: Smoking is an important cause of mortality and recent studies have suggested that even low-intensity smoking might be associated with increased mortality. Still, smoking is associated with lower socio-economic status as well as other potential risk factors, and disease onset might motivate smoking cessation, thus residual confounding and reverse causality might bias results. We aimed to assess the evidence of a causal relationship between smoking intensity and cause-specific as well as all-cause-mortality using Mendelian randomization analyses., Methods: We included 56 019 participants from the Norwegian HUNT2 Study and 337 103 participants from UK Biobank, linked to national registry data on causes of death. We estimated associations of self-reported smoking as well as the genetic variant rs1051730 as an instrument for smoking intensity with all-cause and cause-specific mortality. We subsequently meta-analysed the results from the two cohorts., Results: Each effect allele of the rs1051730 was associated with a 9% increased hazard of all-cause mortality [95% confidence interval (CI) 6-11] among ever smokers. Effect alleles were also associated with death by neoplasms [hazard ratio (HR) 1.11, 95% CI 1.06-1.15], circulatory diseases (HR 1.06, 95% CI 1.01-1.11) and respiratory diseases (HR 1.15, 95% CI 1.05-1.26) among ever smokers. The association was stronger among ever than never smokers for all-cause mortality (p < 0.001), neoplasms (p = 0.001) and respiratory diseases (p = 0.038)., Conclusions: Our results indicate a causal effect of smoking intensity on all-cause mortality and death by neoplasms and respiratory diseases. There was weaker evidence of a causal effect of smoking intensity on death by circulatory diseases., (© The Author(s) 2019; all rights reserved. Published by Oxford University Press on behalf of the International Epidemiological Association.)

Background: Family history of myocardial infarction (FHMI) is known to increase the risk of venous thromboembolism (VTE)., Objectives: To investigate the effect of prothrombotic genotypes on the association between FHMI and VTE in a case-cohort recruited from a general population., Methods: Cases with a first VTE (n = 1493) and a subcohort (n = 13 072) were sampled from the Tromsø study (1994-2012) and the Nord-Trøndelag health (HUNT) study (1995-2008). The DNA samples were genotyped for rs8176719 (ABO), rs6025 (F5), rs1799963 (F2), rs2066865 (FGG), and rs2036914 (F11). Participants with missing information on risk alleles (n = 175), FHMI (n = 2769), and BMI (n = 52) were excluded. Cox regression models were used to estimate hazard ratios (HRs) with 95% confidence intervals (CI) for VTE. To explore the role of prothrombotic genotypes for the association between FHMI and VTE, we (a) included the genotypes in the multivariable-adjusted models and (b) assessed the joint effects between FHMI and genotypes on VTE risk., Results: The FHMI was associated with a 1.3-fold increased risk of VTE (HR 1.32, 95% CI 1.16-1.50) and 1.5-fold increased risk of unprovoked VTE (HR 1.47, 95% CI 1.22-1.78). The risk of VTE by FHMI did not alter after adjustment for the five genotypes. The combination of FHMI and the different prothrombotic genotypes did not result in an excess VTE risk (i.e. no biological interaction)., Conclusions: Our findings suggest that the risk of VTE by FHMI is not explained by rs8176719 (ABO), rs6025 (F5), rs1799963 (F2), rs2066865 (FGG), and rs2036914 (F11). The combination of FHMI with prothrombotic genotypes had an additive effect on VTE risk., (© 2019 International Society on Thrombosis and Haemostasis.)

Chronic kidney disease (CKD) is a growing health burden currently affecting 10-15% of adults worldwide. Estimated glomerular filtration rate (eGFR) as a marker of kidney function is commonly used to diagnose CKD. We analyze eGFR data from the Nord-Trøndelag Health Study and Michigan Genomics Initiative and perform a GWAS meta-analysis with public summary statistics, more than doubling the sample size of previous meta-analyses. We identify 147 loci (53 novel) associated with eGFR, including genes involved in transcriptional regulation, kidney development, cellular signaling, metabolism, and solute transport. Additionally, sex-stratified analysis identifies one locus with more significant effects in women than men. Using genetic risk scores constructed from these eGFR meta-analysis results, we show that associated variants are generally predictive of CKD with only modest improvements in detection compared with other known clinical risk factors. Collectively, these results yield additional insight into the genetic factors underlying kidney function and progression to CKD.

Tobacco and alcohol use are leading causes of mortality that influence risk for many complex diseases and disorders 1 . They are heritable 2,3 and etiologically related 4,5 behaviors that have been resistant to gene discovery efforts 6-11 . In sample sizes up to 1.2 million individuals, we discovered 566 genetic variants in 406 loci associated with multiple stages of tobacco use (initiation, cessation, and heaviness) as well as alcohol use, with 150 loci evidencing pleiotropic association. Smoking phenotypes were positively genetically correlated with many health conditions, whereas alcohol use was negatively correlated with these conditions, such that increased genetic risk for alcohol use is associated with lower disease risk. We report evidence for the involvement of many systems in tobacco and alcohol use, including genes involved in nicotinic, dopaminergic, and glutamatergic neurotransmission. The results provide a solid starting point to evaluate the effects of these loci in model organisms and more precise substance use measures.

To identify genetic variation underlying atrial fibrillation, the most common cardiac arrhythmia, we performed a genome-wide association study of >1,000,000 people, including 60,620 atrial fibrillation cases and 970,216 controls. We identified 142 independent risk variants at 111 loci and prioritized 151 functional candidate genes likely to be involved in atrial fibrillation. Many of the identified risk variants fall near genes where more deleterious mutations have been reported to cause serious heart defects in humans (GATA4, MYH6, NKX2-5, PITX2, TBX5) 1 , or near genes important for striated muscle function and integrity (for example, CFL2, MYH7, PKP2, RBM20, SGCG, SSPN). Pathway and functional enrichment analyses also suggested that many of the putative atrial fibrillation genes act via cardiac structural remodeling, potentially in the form of an 'atrial cardiomyopathy' 2 , either during fetal heart development or as a response to stress in the adult heart.

In genome-wide association studies (GWAS) for thousands of phenotypes in large biobanks, most binary traits have substantially fewer cases than controls. Both of the widely used approaches, the linear mixed model and the recently proposed logistic mixed model, perform poorly; they produce large type I error rates when used to analyze unbalanced case-control phenotypes. Here we propose a scalable and accurate generalized mixed model association test that uses the saddlepoint approximation to calibrate the distribution of score test statistics. This method, SAIGE (Scalable and Accurate Implementation of GEneralized mixed model), provides accurate P values even when case-control ratios are extremely unbalanced. SAIGE uses state-of-art optimization strategies to reduce computational costs; hence, it is applicable to GWAS for thousands of phenotypes by large biobanks. Through the analysis of UK Biobank data of 408,961 samples from white British participants with European ancestry for > 1,400 binary phenotypes, we show that SAIGE can efficiently analyze large sample data, controlling for unbalanced case-control ratios and sample relatedness.