Genetic polymorphisms in IL-2, IL-10 and FOXP3 are associated with autoimmune neutropaenia in early childhood and autoantibody specificity in a Danish cohort.

Autoimmune neutropaenia (AIN) in early childhood is characterized by chronic neutropaenia and positivity for human neutrophil antibodies (HNA), resulting in the excessive destruction of neutrophils. The association between regulatory T cells (Tregs) and AIN has been described, and in this study, we investigated three Treg-associated genes, IL-2, IL-10 and FOXP3. The frequencies of three single nucleotide polymorphisms (SNPs) in IL-2 -330T>G (rs2069762), +114G>T (rs2069763) and IVS3-116 A>G (rs2069772), four SNPs in IL-10 -3575T>A (rs1800890), -1082G>A (rs1800896), -819 C>T (rs1800871) and -592 C>A (rs1800872) and three SNPs in FOXP3 -3499 A>G (rs3761547), -3279 C>A (rs3761548) and -924 A>G (rs2232365) were compared between 166 Danish AIN patients and 358 healthy controls. Disease association was observed for IL-2 IVS3-116 GG (p = 0.0081, OR = 0.35 [0.15-0.80]), IL-10 -3575 TT (p = 0.0078, OR = 1.71 [1.16-2.54]) and IL-10 -1082 AA (p = 0.014, OR = 1.76 [1.14-2.72]) in all patients and FOXP3 -924 (p = 0.0005, A OR = 0.41 [0.25-0.68] and G OR = 2.42 [1.46-4.01]) in male patients. None of the associations were linked to antibody specificity. Disease-associated haplotypes were observed in IL-2 and FOXP3. IL-2 -330T/+114 T/IVS3-116A was associated with anti-FcγRIIIb-positive patients (p = 0.012, OR = 2.07 [1.18-3.62]). FOXP3 -3499A/-3279C/-924A was associated with anti-HNA-1a-positive male patients (p = 0.016, OR = 0.41 [0.20-0.83]), and ACG was associated with female patients, both in the combined group (p = 0.006, OR = NA) and the anti-FcγRIIIb-positive group (p = 0.002, OR = NA). We conclude that our findings reveal a correlation between SNP in Treg-associated genes and AIN, indicating that AIN could be driven by dysfunction of immune homeostatic-evolving Tregs.

CD21 deficiency in 2 siblings and frequency of the associated mutation in the Danish population.

The clinical presentation of CD21 deficiency in 2 siblings caused by a novel mutation in the CD21 gene is reported, and the frequency of this mutation in the Danish population is explored. Successful treatment with IgG replacement in both patients with CD21 deficiency is also reported.

High-resolution HLA genotyping identifies risk alleles in both class I and II for primary autoimmune neutropenia in early childhood in a Danish cohort.

HLA studies in patients with autoimmune neutropenia (AIN) have shown very consistent results for the association with HLA class II alleles at low resolution. This study aimed to examine the association of both HLA class I and class II at high resolution to clarify the contribution of risk alleles to the disease. A total of 107 AIN patients were genotyped for six loci of HLA class I (HLA-A, -B and -C) and class II (HLA-DRB1, -DQB1, and -DPB1) genes by a high-resolution (3-field, 6-digit) analysis and compared with HLA typing of 1000 healthy controls. Compared with the controls, the allele frequencies were significantly higher in AIN patients for A*02:17:01G, C*01:02:01G, DRB1*10:01:01G, DRB1*14:01:01G, DRB1*16:01:01G, DQB1*05:02:01G, and DQB1*05:03:01G but lower significant for C*03:04:01G, DRB1*04:01:01G, DRB1*13:02:01G, DQB1*03:02:01G, and DQB1*06:04:01G. Frequently associated two-locus haplotypes were found to be DRB1*10:01:01G-DQB1*05:01:01G and DRB1*16:01:01G-DQB1*05:02:01G, while the S2 (Q- or D-KRAA) shared epitope (SE) was associated with lower risk. A unique association with HLA alleles was observed between patients with specific anti-HNA-1a antibodies and broad-reacting anti-FcγRIIIb. Anti-HNA-1a antibody-positive patients were associated with C*01:02:01G, DRB1*01:01:01G, DRB1*16:01:01G, DQB1*05:01:01G, DQB1*05:02:01G, DQB1*06:04:01G, and DPB1*10:01:01G; the two-locus haplotypes DRB1*01:01:01G-DQB1*05:01:01G and DRB1*16:01:01G-DQB1*05:02:01G; and the S3P (Q- or R-RRAA) SE. Anti-FcγRIIIb antibody-positive patients were associated with the alleles A*02:17:01G, DRB1*10:01:01G, and DQB1*05:02:01G; the haplotypes DRB1*10:01:01G-DQB1*05:01:01G and DRB1*11:01:02G-DQB1*05:02:01G; and the S3D (DRRAA) SE. The different associations regarding FcγRIIIb antibody specificities could indicate disease heterogeneity.

Lifestyle and demographic associations with 47 inflammatory and vascular stress biomarkers in 9876 blood donors.

BACKGROUND: The emerging use of biomarkers in research and tailored care introduces a need for information about the association between biomarkers and basic demographics and lifestyle factors revealing expectable concentrations in healthy individuals while considering general demographic differences. METHODS: A selection of 47 biomarkers, including markers of inflammation and vascular stress, were measured in plasma samples from 9876 Danish Blood Donor Study participants. Using regression models, we examined the association between biomarkers and sex, age, Body Mass Index (BMI), and smoking. RESULTS: Here we show that concentrations of inflammation and vascular stress biomarkers generally increase with higher age, BMI, and smoking. Sex-specific effects are observed for multiple biomarkers. CONCLUSION: This study provides comprehensive information on concentrations of 47 plasma biomarkers in healthy individuals. The study emphasizes that knowledge about biomarker concentrations in healthy individuals is critical for improved understanding of disease pathology and for tailored care and decision support tools.

Blood-based biomarkers are circulating molecules that can help to indicate health or disease. Biomarker levels may vary depending on demographic and lifestyle factors such as age, sex, smoking status, and body mass index. Here, we examine the effects of these demographic and lifestyle factors on levels of biomarkers related to activation of the immune system and cardiovascular stress. Measurements of 47 different proteins were performed on blood samples from nearly 10,000 healthy Danish blood donors. Measurement data were linked with questionnaire data to assess effects of lifestyle. We found that immune activation and vascular stress generally increased with age, BMI, and smoking. As these measurements are from healthy blood donors they can serve as a reference for expectable effects and inflammation levels in healthy individuals. Knowledge about the healthy state is important for understanding disease progression and optimizing care.

A genome-wide association study of social trust in 33,882 Danish blood donors.

Social trust is a heritable trait that has been linked with physical health and longevity. In this study, we performed genome-wide association studies of self-reported social trust in n = 33,882 Danish blood donors. We observed genome-wide and local evidence of genetic similarity with other brain-related phenotypes and estimated the single nucleotide polymorphism-based heritability of trust to be 6% (95% confidence interval = (2.1, 9.9)). In our discovery cohort (n = 25,819), we identified one significantly associated locus (lead variant: rs12776883) in an intronic enhancer region of PLPP4, a gene highly expressed in brain, kidneys, and testes. However, we could not replicate the signal in an independent set of donors who were phenotyped a year later (n = 8063). In the subsequent meta-analysis, we found a second significantly associated variant (rs71543507) in an intergenic enhancer region. Overall, our work confirms that social trust is heritable, and provides an initial look into the genetic factors that influence it.

World human neutrophil antigens investigation survey.

BACKGROUND AND OBJECTIVES: Human neutrophil antigens (HNAs) are categorized into five systems: HNA-1 to HNA-5. Given the importance of neutrophils in immunity, we sought to create awareness of the role of HNA diagnostic services in managing immune neutropenia and transfusion-related acute lung injury. To provide health communities all around the world with access to these services, we conducted a survey to create a directory of these HNA diagnostic services. MATERIALS AND METHODS: An Excel table-based survey was created to capture information on the laboratory's location and was emailed to 55 individuals with known or possible HNA investigation activity. The collected data were then summarized and analysed. RESULTS: Of contacted laboratories, the surveys were returned from 23 (38.2%) laboratories; 17 have already established HNA diagnostic (of them 12 were regular participants of the International Granulocyte Immunobiology Workshop [ISBT-IGIW]), 4 laboratories were in the process of establishing their HNA investigation and the remaining 2 responder laboratories, did not conduct HNA investigations. In established laboratories, investigation for autoimmune neutropenia (infancies and adults) was the most frequently requested, and antibodies against HNA-1a and HNA-1b were the most commonly detected. CONCLUSION: The directory of survey respondents provides a resource for health professionals wanting to access HNA diagnostic services. The present study offers a comprehensive picture of HNA diagnostics (typing and serology), identifying weak points and areas for improvement for the first time. Identifying more laboratories involved in HNA diagnostics with limited access to international societies in the field will globally improve HNA diagnostics.

Deficit of homozygosity among 1.52 million individuals and genetic causes of recessive lethality.

Genotypes causing pregnancy loss and perinatal mortality are depleted among living individuals and are therefore difficult to find. To explore genetic causes of recessive lethality, we searched for sequence variants with deficit of homozygosity among 1.52 million individuals from six European populations. In this study, we identified 25 genes harboring protein-altering sequence variants with a strong deficit of homozygosity (10% or less of predicted homozygotes). Sequence variants in 12 of the genes cause Mendelian disease under a recessive mode of inheritance, two under a dominant mode, but variants in the remaining 11 have not been reported to cause disease. Sequence variants with a strong deficit of homozygosity are over-represented among genes essential for growth of human cell lines and genes orthologous to mouse genes known to affect viability. The function of these genes gives insight into the genetics of intrauterine lethality. We also identified 1077 genes with homozygous predicted loss-of-function genotypes not previously described, bringing the total set of genes completely knocked out in humans to 4785.

Symptoms of attention deficit hyperactivity disorder are associated with Hidradenitis suppurativa in Danish blood donors.

Hidradenitis suppurativa (HS) is a chronic inflammatory skin disease associated with psychiatric comorbidity. Attention deficit hyperactivity disorder (ADHD) is a mental disorder associated with systemic and skin inflammation such as psoriasis and atopic dermatitis. Whether HS symptoms are associated with ADHD symptoms remains unexplored. Thus, the aim of this study was to explore the possible association between HS and ADHD. Participants in the Danish Blood Donor Study (DBDS) were included in this cross-sectional study during 2015-2017. The participants provided questionnaire data on screening items of HS, ADHD symptoms (ASRS-score), and depressive symptoms, smoking and body mass index (BMI). A logistic regression with HS symptoms as a binary outcome predicted by ADHD adjusted for age, sex, smoking, BMI, and depression was conducted to investigate the association between HS and ADHD. A total of 52,909 Danish blood donors were included in the study. Of these were 1004/52,909 (1.9%) considered participants with HS. Of the participants with HS, 74/996 (7.4%) screened positive of ADHD symptoms, while only 1786/51,129 (3.5%) of the participants without HS screened positive of ADHD. Adjusted for confounders, ADHD was positively associated with HS, odds ratio 1.85 (95% confidence interval: 1.43-2.37). Psychiatric comorbidity of HS is not limited to depression and anxiety. This study shows a positive association between HS and ADHD. Further research on the biological mechanisms behind this association is warranted.

Genetic variants associated with syncope implicate neural and autonomic processes.

AIMS: Syncope is a common and clinically challenging condition. In this study, the genetics of syncope were investigated to seek knowledge about its pathophysiology and prognostic implications. METHODS AND RESULTS: This genome-wide association meta-analysis included 56 071 syncope cases and 890 790 controls from deCODE genetics (Iceland), UK Biobank (United Kingdom), and Copenhagen Hospital Biobank Cardiovascular Study/Danish Blood Donor Study (Denmark), with a follow-up assessment of variants in 22 412 cases and 286 003 controls from Intermountain (Utah, USA) and FinnGen (Finland). The study yielded 18 independent syncope variants, 17 of which were novel. One of the variants, p.Ser140Thr in PTPRN2, affected syncope only when maternally inherited. Another variant associated with a vasovagal reaction during blood donation and five others with heart rate and/or blood pressure regulation, with variable directions of effects. None of the 18 associations could be attributed to cardiovascular or other disorders. Annotation with regard to regulatory elements indicated that the syncope variants were preferentially located in neural-specific regulatory regions. Mendelian randomization analysis supported a causal effect of coronary artery disease on syncope. A polygenic score (PGS) for syncope captured genetic correlation with cardiovascular disorders, diabetes, depression, and shortened lifespan. However, a score based solely on the 18 syncope variants performed similarly to the PGS in detecting syncope risk but did not associate with other disorders. CONCLUSION: The results demonstrate that syncope has a distinct genetic architecture that implicates neural regulatory processes and a complex relationship with heart rate and blood pressure regulation. A shared genetic background with poor cardiovascular health was observed, supporting the importance of a thorough assessment of individuals presenting with syncope.

Genetic variations in low-to-medium-affinity Fcγ receptors and autoimmune neutropenia in early childhood in a Danish cohort.

Autoimmune neutropenia (AIN) in early childhood is caused by autoantibodies directed against antigens on the neutrophil membrane and is a frequent cause of neutropenia in children. Association of AIN with Fcγ receptor (FCGR) 3B variants is well described. In this study, we investigate genetic variations in the FCGR locus and copy number variation of FCGR3B. A total of 130 antibody-positive AIN patients, 64 with specific anti-HNA-1a antibodies and 66 with broad-reacting anti-FcγRIIIb antibodies, were genotyped with a multiplex ligation probe assay and compared with healthy controls. Positive findings were confirmed with real-time q-PCR. We determined copy numbers of the FCGR2 and FCGR3 genes and the following SNPs: FCGR2A Q62W (rs201218628), FCGR2A H166R (rs1801274), FCGR2B I232T (rs1050501), FCGR3A V176F (rs396991), haplotypes for FCGR2B/C promoters (rs3219018/rs780467580), FCGR2C STOP/ORF and HNA-1 genotypes in FCGR3B (rs447536, rs448740, rs52820103, rs428888 and rs2290834). Generally, associations were antibody specific, with all associations being representative of the anti-HNA-1a-positive group, while the only association found in the anti-FcγRIIIb group was with the HNA-1 genotype. An increased risk of AIN was observed for patients with one copy of FCGR3B; the HNA genotypes HNA-1a, HNA-1aa or HNA-1aac; the FCGR2A 166H and FCGR2B 232I variations; and no copies of FCGR2B 2B.4. A decreased risk was observed for HNA genotype HNA-1bb; FCGR2A 166R; FCGR2B 232T; and one copy of FCGR2B promoter 2B.4. We conclude that in our Danish cohort, there was a strong association between variation in the FCGR locus and AIN. The findings of different genetic associations between autoantibody groups could indicate the presence of two different disease entities and disease heterogeneity.

Genome-wide meta-analysis identifies 93 risk loci and enables risk prediction equivalent to monogenic forms of venous thromboembolism.

We report a genome-wide association study of venous thromboembolism (VTE) incorporating 81,190 cases and 1,419,671 controls sampled from six cohorts. We identify 93 risk loci, of which 62 are previously unreported. Many of the identified risk loci are at genes encoding proteins with functions converging on the coagulation cascade or platelet function. A VTE polygenic risk score (PRS) enabled effective identification of both high- and low-risk individuals. Individuals within the top 0.1% of PRS distribution had a VTE risk similar to homozygous or compound heterozygous carriers of the variants G20210A (c.*97 G > A) in F2 and p.R534Q in F5. We also document that F2 and F5 mutation carriers in the bottom 10% of the PRS distribution had a risk similar to that of the general population. We further show that PRS improved individual risk prediction beyond that of genetic and clinical risk factors. We investigated the extent to which venous and arterial thrombosis share clinical risk factors using Mendelian randomization, finding that some risk factors for arterial thrombosis were directionally concordant with VTE risk (for example, body mass index and smoking) whereas others were discordant (for example, systolic blood pressure and triglyceride levels).

Development and validation of circulating protein signatures as diagnostic biomarkers for biliary tract cancer.

Background & Aims: Biliary tract cancer (BTC) is associated with a dismal prognosis, partly because it is typically diagnosed late, highlighting the need for diagnostic biomarkers. The purpose of this project was to identify and validate multiprotein signatures that could differentiate patients with BTC from non-cancer controls. Methods: In this study, we included treatment-naïve patients with BTC, healthy controls, and patients with benign conditions including benign biliary tract disease. Participants were divided into three non-overlapping cohorts: a case-control-based discovery cohort (BTC = 186, controls = 249); a case-control-based validation cohort (validation cohort 1: BTC = 113, controls = 241); and a cohort study-based validation cohort including participants (BTC = 8, controls = 132) referred for diagnostic work-up for suspected cancer (validation cohort 2). Immuno-Oncology (I-O)-related proteins were measured in serum and plasma using a proximity extension assay (Olink Proteomics). Lasso and Ridge regressions were used to generate protein signatures of I-O-related proteins and carbohydrate antigen 19-9 (CA19-9) in the discovery cohort. Results: Sixteen protein signatures, including 2 to 82 proteins, were generated. All signatures included CA19-9 and chemokine C-C motif ligand 20. Signatures discriminated between patients with BTC vs. controls, with AUCs ranging from 0.95 to 0.99 in the discovery cohort and 0.94 to 0.97 in validation cohort 1. In validation cohort 2, AUCs ranged from 0.84 to 0.94. Nine signatures achieved a specificity of 82% to 84% while keeping a sensitivity of 100% in validation cohort 2. All signatures performed better than CA19-9, and signatures including >15 proteins showed the best performance. Conclusion: The study demonstrated that it is possible to generate protein signatures that can successfully differentiate patients with BTC from non-cancer controls. Impact and implications: We attempted to find blood sample-based protein profiles that could differentiate patients with biliary tract cancer from those without cancer. Several profiles were found and tested in different groups of patients. The profiles were successful at identifying most patients with biliary tract cancer, pointing towards the utility of multiprotein signatures in this context.

Genome-wide association meta-analysis of knee and hip osteoarthritis uncovers genetic differences between patients treated with joint replacement and patients without joint replacement.

OBJECTIVES: Osteoarthritis is a common and severe, multifactorial disease with a well-established genetic component. However, little is known about how genetics affect disease progression, and thereby the need for joint placement. Therefore, we aimed to investigate whether the genetic associations of knee and hip osteoarthritis differ between patients treated with joint replacement and patients without joint replacement. METHODS: We included knee and hip osteoarthritis cases along with healthy controls, altogether counting >700 000 individuals. The cases were divided into two groups based on joint replacement status (surgical vs non-surgical) and included in four genome-wide association meta-analyses: surgical knee osteoarthritis (N = 22 525), non-surgical knee osteoarthritis (N = 38 626), surgical hip osteoarthritis (N = 20 221) and non-surgical hip osteoarthritis (N = 17 847). In addition, we tested for genetic correlation between the osteoarthritis groups and the pain phenotypes intervertebral disc disorder, dorsalgia, fibromyalgia, migraine and joint pain. RESULTS: We identified 52 sequence variants associated with knee osteoarthritis (surgical: 17, non-surgical: 3) or hip osteoarthritis (surgical: 34, non-surgical: 1). For the surgical phenotypes, we identified 10 novel variants, including genes involved in autophagy (rs2447606 in ATG7) and mechanotransduction (rs202127176 in PIEZO1). One variant, rs13107325 in SLC39A8, associated more strongly with non-surgical knee osteoarthritis than surgical knee osteoarthritis. For all other variants, significance and effect sizes were higher for the surgical phenotypes. In contrast, genetic correlations with pain phenotypes tended to be stronger in the non-surgical groups. CONCLUSIONS: Our results indicate differences in genetic associations between knee and hip osteoarthritis depending on joint replacement status.

The Social and Psychological Impact of Acne Treatment: A Cross-Sectional Study of Blood Donors.

BACKGROUND: Acne in adolescence and adulthood is believed to have a long-term impact on socioeconomic status (SES) and health-related quality-of-life (HRQoL) in adults. OBJECTIVE: To estimate the cross-sectional prevalence of medically treated (MedTreAc) and untreated acne (UnTreAc) and to characterize its long-term impact in adults. METHODS: A nationwide cross-sectional study on 17 428 blood donors aged 18-35 was performed. Associations among acne and HRQoL, depressive symptoms, total income, and SES were investigated via linear/logistic/multinomial logistic regression analyses adjusted for relevant covariables. HRQoL was measured by the Short Form-12, and depressive symptoms by the Major Depression Inventory. The data were self-reported. RESULTS: Of the participants, 3591 (20.6%) and 1354 (7.8%) identified as the MedTreAc and UnTreAc phenotype, respectively. Neither phenotype was associated with a long-term impact on total income, but the MedTreAc group was associated with being an apprentice/student (OR = 1.26; 95% CI: 1.12, 1.42; P = 1.3×10-4) or high skill-level employee (OR = 1.22, 95% CI: 1.07; 1.39, P = .0023), while self-employment was more common for those with UnTreAc (OR = 1.53; 95% CI: 1.12, 2.06, P = .0061). Additionally, the UnTreAc group was associated with a lower mental HRQoL (SF-12 mental component summary score -1.05, 95% CI: -1.56, -0.54; P = 1.4×10-9) and increased odds ratio of depressive symptoms (OR = 1.44; 95% CI: 1.00, 2.02, P = .046). CONCLUSION: In this population of blood donors, the cumulative prevalence of MedTreAc and UnTreAc were 20.6% and 7.8%, respectively. Untreated acne had a long-term impact on psychosocial well-being in adulthood. It was associated with lower mental HRQoL and higher occurrence of depressive symptoms. Acne was not associated with a lower salary or SES.

ABO, secretor, and Lewis carbohydrate histo-blood groups are associated with autoimmune neutropenia of early childhood in Danish patients.

BACKGROUND: Autoimmune neutropenia of early childhood (AIN) is caused by autoantibodies directed against antigens on the neutrophil membrane. The ABO, secretor, and Lewis histo-blood group systems control the expression of carbohydrate antigens and have previously been linked to autoimmune diseases. We aimed to investigate the association between genotypes and the risk of AIN in Danish patients. STUDY DESIGN AND METHODS: One hundred fifty-four antibody-positive AIN patients were included. Controls (n = 400) were healthy unrelated Danish blood donors. Molecular determination of ABO, secretor (FUT2), and Lewis (FUT3) genotypes were determined using real-time polymerase chain reaction (qPCR) or Sanger sequencing to infer the prevalence of Lewis antigens (Lea and Leb ) and secretor (SeSe or Sese) or nonsecretor (sese) phenotypes. RESULTS: Blood type O was more common in controls (46.8%) than in AIN patients (36.4%) (OR = 0.65; p = 0.028). Secretors of H Leb antigens were less frequent among AIN patients (25.2%) than controls (35.0%) (OR = 0.62; p = 0.037). DISCUSSION: ABO blood group antigens and the secretion of these antigens are associated with a diagnosis of AIN. The mechanism underlying the association between autoimmunity and interaction among ABO, secretor, and Lewis genotypes has not yet been elucidated, but several studies indicate a connection to the gut microbiota.

[Myocarditis and development of dilated cardiomyopathy].

Inflammation is increasingly recognised as a causal factor in the development and progression of cardiovascular disease. With the introduction of immune checkpoint inhibitors in oncology and the ongoing COVID-19 pandemic the role of the immune system in myocardial inflammation (myocarditis) and subsequent inflammatory cardiomyopathy has once again regained attention. In this review, we want to bring myocardial inflammation to the clinician's attention and provide up-to-date knowledge on its diagnostic workup, prognostication, and current management recommendations.

Rare SLC13A1 variants associate with intervertebral disc disorder highlighting role of sulfate in disc pathology.

Back pain is a common and debilitating disorder with largely unknown underlying biology. Here we report a genome-wide association study of back pain using diagnoses assigned in clinical practice; dorsalgia (119,100 cases, 909,847 controls) and intervertebral disc disorder (IDD) (58,854 cases, 922,958 controls). We identify 41 variants at 33 loci. The most significant association (ORIDD = 0.92, P = 1.6 × 10-39; ORdorsalgia = 0.92, P = 7.2 × 10-15) is with a 3'UTR variant (rs1871452-T) in CHST3, encoding a sulfotransferase enzyme expressed in intervertebral discs. The largest effects on IDD are conferred by rare (MAF = 0.07 - 0.32%) loss-of-function (LoF) variants in SLC13A1, encoding a sodium-sulfate co-transporter (LoF burden OR = 1.44, P = 3.1 × 10-11); variants that also associate with reduced serum sulfate. Genes implicated by this study are involved in cartilage and bone biology, as well as neurological and inflammatory processes.

Genome-wide association study of febrile seizures implicates fever response and neuronal excitability genes.

Febrile seizures represent the most common type of pathological brain activity in young children and are influenced by genetic, environmental and developmental factors. In a minority of cases, febrile seizures precede later development of epilepsy. We conducted a genome-wide association study of febrile seizures in 7635 cases and 83 966 controls identifying and replicating seven new loci, all with P < 5 × 10-10. Variants at two loci were functionally related to altered expression of the fever response genes PTGER3 and IL10, and four other loci harboured genes (BSN, ERC2, GABRG2, HERC1) influencing neuronal excitability by regulating neurotransmitter release and binding, vesicular transport or membrane trafficking at the synapse. Four previously reported loci (SCN1A, SCN2A, ANO3 and 12q21.33) were all confirmed. Collectively, the seven novel and four previously reported loci explained 2.8% of the variance in liability to febrile seizures, and the single nucleotide polymorphism heritability based on all common autosomal single nucleotide polymorphisms was 10.8%. GABRG2, SCN1A and SCN2A are well-established epilepsy genes and, overall, we found positive genetic correlations with epilepsies (rg = 0.39, P = 1.68 × 10-4). Further, we found that higher polygenic risk scores for febrile seizures were associated with epilepsy and with history of hospital admission for febrile seizures. Finally, we found that polygenic risk of febrile seizures was lower in febrile seizure patients with neuropsychiatric disease compared to febrile seizure patients in a general population sample. In conclusion, this largest genetic investigation of febrile seizures to date implicates central fever response genes as well as genes affecting neuronal excitability, including several known epilepsy genes. Further functional and genetic studies based on these findings will provide important insights into the complex pathophysiological processes of seizures with and without fever.