Integrative common and rare variant analyses provide insights into the genetic architecture of liver cirrhosis.

We report a multi-ancestry genome-wide association study on liver cirrhosis and its associated endophenotypes, alanine aminotransferase (ALT) and γ-glutamyl transferase. Using data from 12 cohorts, including 18,265 cases with cirrhosis, 1,782,047 controls, up to 1 million individuals with liver function tests and a validation cohort of 21,689 cases and 617,729 controls, we identify and validate 14 risk associations for cirrhosis. Many variants are located near genes involved in hepatic lipid metabolism. One of these, PNPLA3 p.Ile148Met, interacts with alcohol intake, obesity and diabetes on the risk of cirrhosis and hepatocellular carcinoma (HCC). We develop a polygenic risk score that associates with the progression from cirrhosis to HCC. By focusing on prioritized genes from common variant analyses, we find that rare coding variants in GPAM associate with lower ALT, supporting GPAM as a potential target for therapeutic inhibition. In conclusion, this study provides insights into the genetic underpinnings of cirrhosis.

Genome-wide analyses identify 21 infertility loci and over 400 reproductive hormone loci across the allele frequency spectrum.

Genome-wide association studies (GWASs) may help inform treatments for infertility, whose causes remain unknown in many cases. Here we present GWAS meta-analyses across six cohorts for male and female infertility in up to 41,200 cases and 687,005 controls. We identified 21 genetic risk loci for infertility (P≤5E-08), of which 12 have not been reported for any reproductive condition. We found positive genetic correlations between endometriosis and all-cause female infertility (rg=0.585, P=8.98E-14), and between polycystic ovary syndrome and anovulatory infertility (rg=0.403, P=2.16E-03). The evolutionary persistence of female infertility-risk alleles in EBAG9 may be explained by recent directional selection. We additionally identified up to 269 genetic loci associated with follicle-stimulating hormone (FSH), luteinising hormone, oestradiol, and testosterone through sex-specific GWAS meta-analyses (N=6,095-246,862). While hormone-associated variants near FSHB and ARL14EP colocalised with signals for anovulatory infertility, we found no rg between female infertility and reproductive hormones (P>0.05). Exome sequencing analyses in the UK Biobank (N=197,340) revealed that women carrying testosterone-lowering rare variants in GPC2 were at higher risk of infertility (OR=2.63, P=1.25E-03). Taken together, our results suggest that while individual genes associated with hormone regulation may be relevant for fertility, there is limited genetic evidence for correlation between reproductive hormones and infertility at the population level. We provide the first comprehensive view of the genetic architecture of infertility across multiple diagnostic criteria in men and women, and characterise its relationship to other health conditions.

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.

Meta-analysis of ACE inhibitor-induced angioedema identifies novel risk locus.

BACKGROUND: Angioedema is a rare but potentially life-threatening adverse drug reaction in patients receiving angiotensin-converting enzyme inhibitors (ACEis). Research suggests that susceptibility to ACEi-induced angioedema (ACEi-AE) involves both genetic and nongenetic risk factors. Genome- and exome-wide studies of ACEi-AE have identified the first genetic risk loci. However, understanding of the underlying pathophysiology remains limited. OBJECTIVE: We sought to identify further genetic factors of ACEi-AE to eventually gain a deeper understanding of its pathophysiology. METHODS: By combining data from 8 cohorts, a genome-wide association study meta-analysis was performed in more than 1000 European patients with ACEi-AE. Secondary bioinformatic analyses were conducted to fine-map associated loci, identify relevant genes and pathways, and assess the genetic overlap between ACEi-AE and other traits. Finally, an exploratory cross-ancestry analysis was performed to assess shared genetic factors in European and African-American patients with ACEi-AE. RESULTS: Three genome-wide significant risk loci were identified. One of these, located on chromosome 20q11.22, has not been implicated previously in ACEi-AE. Integrative secondary analyses highlighted previously reported genes (BDKRB2 [bradykinin receptor B2] and F5 [coagulation factor 5]) as well as biologically plausible novel candidate genes (PROCR [protein C receptor] and EDEM2 [endoplasmic reticulum degradation enhancing alpha-mannosidase like protein 2]). Lead variants at the risk loci were found with similar effect sizes and directions in an African-American cohort. CONCLUSIONS: The present results contributed to a deeper understanding of the pathophysiology of ACEi-AE by (1) providing further evidence for the involvement of bradykinin signaling and coagulation pathways and (2) suggesting, for the first time, the involvement of the fibrinolysis pathway in this adverse drug reaction. An exploratory cross-ancestry comparison implicated the relevance of the associated risk loci across diverse ancestries.

Blood donation and migraine relief: A national population cohort study in Denmark.

INTRODUCTION: Migraine is a prevalent neurological headache disorder. Due to challenges associated with finding effective treatment, many individuals with migraine feel compelled to explore alternative treatment strategies, such as blood donation, hypothesized to provide migraine relief. METHODS: Through logistic, Poisson, and Cox regression methods, we examined the links between migraine and blood donation activities in two population cohorts: Danish blood donors in the Scandinavian Donations and Transfusions Database (SCANDAT-DK, N >1 million) and the Danish Blood Donor Study (N ~ 100,000). RESULTS: SCANDAT-DK analyses showed no link between migraine and the propensity to become a blood donor among males (odds ratio [OR]Males = 0.95 [95% Confidence Interval: 0.86-1.04], and a reduced propensity among females ORFemales = 0.88 [0.83-0.93]). The incidence of migraine was not reduced upon blood donation (standardized incidence ratio [SIR]Males = 0.94 [0.83-1.06]; SIRFemales = 1.04 [0.99-1.10]). Donors with migraine demonstrated longer intervals between donations (hazard ratio [HR]Males = 0.87 [0.85-0.91], HRFemales = 0.80 [0.78-0.82]), and an increased risk of donor lapse (ORMales = 1.23 [1.14-1.32]; ORFemales = 1.28 [1.22-1.33]). Results were corroborated in DBDS using self-reported migraine. Genetic predisposition to migraine associated with longer intervals in females (HRFemales = 0.98 [0.97-0.99]), but not in males. DISCUSSION: Our findings do not support the hypothesis that blood donation serves as a viable treatment strategy among migraine patients. Future prospective investigations may help to elucidate the underlying biological mechanisms by which blood donation may influence migraine pathology.

Variants at the Interleukin 1 Gene Locus and Pericarditis.

Importance: Recurrent pericarditis is a treatment challenge and often a debilitating condition. Drugs inhibiting interleukin 1 cytokines are a promising new treatment option, but their use is based on scarce biological evidence and clinical trials of modest sizes, and the contributions of innate and adaptive immune processes to the pathophysiology are incompletely understood. Objective: To use human genomics, transcriptomics, and proteomics to shed light on the pathogenesis of pericarditis. Design, Setting, and Participants: This was a meta-analysis of genome-wide association studies of pericarditis from 5 countries. Associations were examined between the pericarditis-associated variants and pericarditis subtypes (including recurrent pericarditis) and secondary phenotypes. To explore mechanisms, associations with messenger RNA expression (cis-eQTL), plasma protein levels (pQTL), and CpG methylation of DNA (ASM-QTL) were assessed. Data from Iceland (deCODE genetics, 1983-2020), Denmark (Copenhagen Hospital Biobank/Danish Blood Donor Study, 1977-2022), the UK (UK Biobank, 1953-2021), the US (Intermountain, 1996-2022), and Finland (FinnGen, 1970-2022) were included. Data were analyzed from September 2022 to August 2023. Exposure: Genotype. Main Outcomes and Measures: Pericarditis. Results: In this genome-wide association study of 4894 individuals with pericarditis (mean [SD] age at diagnosis, 51.4 [17.9] years, 2734 [67.6%] male, excluding the FinnGen cohort), associations were identified with 2 independent common intergenic variants at the interleukin 1 locus on chromosome 2q14. The lead variant was rs12992780 (T) (effect allele frequency [EAF], 31%-40%; odds ratio [OR], 0.83; 95% CI, 0.79-0.87; P = 6.67 × 10-16), downstream of IL1B and the secondary variant rs7575402 (A or T) (EAF, 45%-55%; adjusted OR, 0.89; 95% CI, 0.85-0.93; adjusted P = 9.6 × 10-8). The lead variant rs12992780 had a smaller odds ratio for recurrent pericarditis (0.76) than the acute form (0.86) (P for heterogeneity = .03) and rs7575402 was associated with CpG methylation overlapping binding sites of 4 transcription factors known to regulate interleukin 1 production: PU.1 (encoded by SPI1), STAT1, STAT3, and CCAAT/enhancer-binding protein ß (encoded by CEBPB). Conclusions and Relevance: This study found an association between pericarditis and 2 independent sequence variants at the interleukin 1 gene locus. This finding has the potential to contribute to development of more targeted and personalized therapy of pericarditis with interleukin 1-blocking drugs.

Obesity variants in the GIPR gene do not associate with risk of fracture or bone mineral density.

OBJECTIVE: The objective of this study was to investigate the risk of fracture and bone mineral density (BMD) of sequence variants in GIPR that reduce the activity of the GIPR receptor and have been associated with reduced body mass index (BMI). METHODS: We analysed the association of three missense variants in GIPR, a common variant, rs1800437 (p.Glu354Gln), and two rare variants, rs139215588 (p.Arg190Gln) and rs143430880 (p.Glu288Gly), as well as a burden of predicted loss of function (LoF) variants with risk of fracture and with BMD in a large meta-analysis of up to 1.2 million participants. We analysed associations with fractures at different skeletal sites in the general population; any fractures, hip fractures, vertebral fractures and forearm fractures, and specifically non-vertebral and osteoporotic fractures in postmenopausal women. We also evaluated associations with BMD at the lumbar spine, femoral neck, and total body measured with dual-energy X-ray absorptiometry (DXA), and with BMD estimated from heel ultrasound (eBMD). RESULTS: None of the three missense variants in GIPR associated significantly with increased risk of fractures or with lower BMD. Burden of LoF variants in GIPR were not associated with fractures or with BMD measured with clinically validated DXA, but associated with eBMD. CONCLUSION: Missense variants in GIPR, or burden of LoF variants in the gene, do not associate with risk of fractures or with lower BMD.

Genetic prediction of 33 blood group phenotypes using an existing genotype dataset.

BACKGROUND: Accurate blood type data are essential for blood bank management, but due to costs, few of 43 blood group systems are routinely determined in Danish blood banks. However, a more comprehensive dataset of blood types is useful in scenarios such as rare blood type allocation. We aimed to investigate the viability and accuracy of predicting blood types by leveraging an existing dataset of imputed genotypes for two cohorts of approximately 90,000 each (Danish Blood Donor Study and Copenhagen Biobank) and present a more comprehensive overview of blood types for our Danish donor cohort. STUDY DESIGN AND METHODS: Blood types were predicted from genome array data using known variant determinants. Prediction accuracy was confirmed by comparing with preexisting serological blood types. The Vel blood group was used to test the viability of using genetic prediction to narrow down the list of candidate donors with rare blood types. RESULTS: Predicted phenotypes showed a high balanced accuracy >99.5% in most cases: A, B, C/c, Coa /Cob , Doa /Dob , E/e, Jka /Jkb , Kna /Knb , Kpa /Kpb , M/N, S/s, Sda , Se, and Yta /Ytb , while some performed slightly worse: Fya /Fyb , K/k, Lua /Lub , and Vel ~99%-98% and CW and P1 ~96%. Genetic prediction identified 70 potential Vel negatives in our cohort, 64 of whom were confirmed correct using polymerase chain reaction (negative predictive value: 91.5%). DISCUSSION: High genetic prediction accuracy in most blood groups demonstrated the viability of generating blood types using preexisting genotype data at no cost and successfully narrowed the pool of potential individuals with the rare Vel-negative phenotype from 180,000 to 70.

An atlas of genetic determinants of forearm fracture.

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.

Rare variants with large effects provide functional insights into the pathology of migraine subtypes, with and without aura.

Migraine is a complex neurovascular disease with a range of severity and symptoms, yet mostly studied as one phenotype in genome-wide association studies (GWAS). Here we combine large GWAS datasets from six European populations to study the main migraine subtypes, migraine with aura (MA) and migraine without aura (MO). We identified four new MA-associated variants (in PRRT2, PALMD, ABO and LRRK2) and classified 13 MO-associated variants. Rare variants with large effects highlight three genes. A rare frameshift variant in brain-expressed PRRT2 confers large risk of MA and epilepsy, but not MO. A burden test of rare loss-of-function variants in SCN11A, encoding a neuron-expressed sodium channel with a key role in pain sensation, shows strong protection against migraine. Finally, a rare variant with cis-regulatory effects on KCNK5 confers large protection against migraine and brain aneurysms. Our findings offer new insights with therapeutic potential into the complex biology of migraine and its subtypes.

Complex effects of sequence variants on lipid levels and coronary artery disease.

Many sequence variants have additive effects on blood lipid levels and, through that, on the risk of coronary artery disease (CAD). We show that variants also have non-additive effects and interact to affect lipid levels as well as affecting variance and correlations. Variance and correlation effects are often signatures of epistasis or gene-environmental interactions. These complex effects can translate into CAD risk. For example, Trp154Ter in FUT2 protects against CAD among subjects with the A1 blood group, whereas it associates with greater risk of CAD in others. His48Arg in ADH1B interacts with alcohol consumption to affect lipid levels and CAD. The effect of variants in TM6SF2 on blood lipids is greatest among those who never eat oily fish but absent from those who often do. This work demonstrates that variants that affect variance of quantitative traits can allow for the discovery of epistasis and interactions of variants with the environment.

First cases of SARS-CoV-2 BA.2.86 in Denmark, 2023.

We describe 10 cases of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant BA.2.86 detected in Denmark, including molecular characteristics and results from wastewater surveillance that indicate that the variant is circulating in the country at a low level. This new variant with many spike gene mutations was classified as a variant under monitoring by the World Health Organization on 17 August 2023. Further global monitoring of COVID-19, BA.2.86 and other SARS-CoV-2 variants is highly warranted.

Pregnancy-Associated Bleeding and Genetics: Five Sequence Variants in the Myometrium and Progesterone Signaling Pathway are associated with postpartum hemorrhage.

Bleeding in early pregnancy and postpartum hemorrhage (PPH) bear substantial risks, with the former closely associated with pregnancy loss and the latter being the foremost cause of maternal death, underscoring the severity of these complications in maternal-fetal health. Here, we investigated the genetic variation underlying aspects of pregnancy-associated bleeding and identified five loci associated with PPH through a meta-analysis of 21,512 cases and 259,500 controls. Functional annotation analysis indicated candidate genes, HAND2, TBX3, and RAP2C/FRMD7, at three loci and showed that at each locus, associated variants were located within binding sites for progesterone receptors (PGR). Furthermore, there were strong genetic correlations with birth weight, gestational duration, and uterine fibroids. Early bleeding during pregnancy (28,898 cases and 302,894 controls) yielded no genome-wide association signals, but showed strong genetic correlation with a variety of human traits, indicative of polygenic and pleiotropic effects. Our results suggest that postpartum bleeding is related to myometrium dysregulation, whereas early bleeding is a complex trait related to underlying health and possibly socioeconomic status.

Medium term moderate to low-level air pollution exposure is associated with higher C-reactive protein among healthy Danish blood donors.

Air pollution is a significant contributor to the global burden of disease with a plethora of associated health effects such as pulmonary and systemic inflammation. C-reactive protein (CRP) is associated with a wide range of diseases and is associated with several exposures. Studies on the effect of air pollution exposure on CRP levels in low to moderate pollution settings have shown inconsistent results. In this cross-sectional study high sensitivity CRP measurements on 18,463 Danish blood donors were linked to modelled air pollution data for NOx, NO2, O3, CO, SO2, NH3, mineral dust, black carbon, organic carbon, sea salt, secondary inorganic aerosols and its components, primary PM2.5, secondary organic aerosols, total PM2.5, and total PM10 at their residential address over the previous month. Associations were analysed using ordered logistic regression with CRP quartile as individuals outcome and air pollution exposure as scaled deciles. Analyses were adjusted for health related and socioeconomic covariates using health questionnaires and Danish register data. Exposure to different air pollution components was generally associated with higher CRP (odds ratio estimates ranging from 1.11 to 1.67), while exposure to a few air pollution components was associated with lower CRP. For example, exposure to NO2 increased the odds of high CRP 1.32-fold (95%CI 1.16-1.49), while exposure to NH3 decreased the odds of high CRP 0.81-fold (95%CI 0.73-0.89). This large study among healthy individuals found air pollution exposure to be associated with increased levels of CRP even in a setting with low to moderate air pollution levels.

Population-Based Characterization of Menstrual Migraine and Proposed Diagnostic Criteria.

Importance: There is a need for better recognition and more extensive research into menstrual migraine (MM) in the general population, and a revision of the diagnostic criteria for MM is warranted to move the field forward. Increased understanding of MM is crucial for improving clinical care, diagnosis, and therapy for MM. Objectives: To assess the clinical characteristics of MM, including severity and treatment response, and to propose new diagnostic criteria for pure MM and menstrually related migraine. Design, Setting, and Participants: This is a case-control study of Danish individuals with migraine. All individuals completed a 105-item validated diagnostic migraine questionnaire, sent via the Danish electronic mailing system (e-Boks) between May and August 2020, allowing diagnosis of pure MM and menstrually related migraine by the International Classification of Headache Disorders, Third Edition (ICHD-3). Data analysis was performed from September 2021 to November 2022. Exposure: Diagnosis of migraine. Main Outcomes and Measures: Clinical characteristics of women with MM and women with nonmenstrual migraine (non-MM) were compared using the ICHD-3 diagnostic criteria. A simulation of the risk of randomly misclassifying MM was based on number of migraine attacks during 3 menstrual cycles (3 × 28 days), and simulation analyses were performed using 100 000 permutations of random migraine attacks in migraine patients. Results: A total of 12 618 individuals, including 9184 women, with migraine participated in the study. Among the women with migraine, the prevalence of MM was 16.6% (1532 women), and the prevalence of non-MM was 45.9% (4216 women). The mean (SD) age was 38.7 (8.7) years for women with MM and 37.0 (9.2) years for women with non-MM. Of the 1532 women with MM, 410 (26.8%) fulfilled ICHD-3 diagnostic criteria for pure MM, 1037 (67.7%) fulfilled ICHD-3 diagnostic criteria for menstrually related migraine, and 152 (9.9%) fulfilled proposed diagnostic criteria for rare pure MM. MM was associated with a higher frequency of migraine-accompanying symptoms (odds ratio [OR], 1.98; 95% CI, 1.71-2.29), more frequent (OR, 7.21; 95% CI, 5.77-9.03) and more severe (OR, 1.17; 95% CI, 1.13-1.21) migraine attacks, lower frequency of nonmigraine headache (OR, 0.31; 95% CI, 0.18-0.49), an overall greater response to treatment with triptans (OR, 1.66; 95% CI, 1.24-2.24), better improvement of migraine attacks during late pregnancy (OR, 5.10; 95% CI, 2.17-14.00), and faster reappearance of migraine attacks post partum (OR, 3.19; 95% CI, 2.40-4.25). Hormonal contraceptive-related MM was associated with a higher prevalence of migraine without aura than migraine related to spontaneous menstruation (OR, 1.82; 95% CI, 1.62-2.06). Otherwise, no differences between hormonal and spontaneous MM were observed. The risk of random diagnostic misclassification of ICHD-3 menstrually related migraine in women with high frequency episodic migraine was 43%. This risk was reduced to 3% when applying the proposed criteria for menstrually related migraine. Conclusions and Relevance: In this case-control study, MM in the general population had clinical characteristics that were quantitively different from those of non-MM. Detailed descriptive data and suggested improved diagnostic criteria for pure MM and menstrually related migraine were provided.

Developmental language disorder - a comprehensive study of more than 46,000 individuals.

Developmental language disorder (DLD) is characterized by enduring low language abilities with a significant functional impact, in the absence of biomedical conditions in which language impairment is part of a complex of impairments. There is a lack of awareness of DLD even among healthcare professionals. Here we estimated the prevalence of DLD and its links to reading and learning difficulties and physical and mental health in the Danish Blood Donor Study (N = 46,547), where DLD-related information is based on questionnaires (self-report). We compared the questionnaire-derived DLD status with the relevant language-related diagnoses from hospital registers. We also investigated the genetic architecture of DLD in a subset of the cohort (N = 18,380). DLD was significantly associated with reading and learning difficulties and poorer mental and physical health. DLD prevalence was 3.36%-3.70% based on questionnaires, compared with 0.04% in hospital registers. Our genetic analyses identified one genome-wide significant locus, but not a significant heritability estimate. Our study shows that DLD has health-related implications that may last into adulthood, and that DLD may be undiagnosed in general healthcare. Furthermore, DLD is likely more genetically heterogeneous than narrower developmental language phenotypes. Our results emphasize the need to raise awareness of DLD and consider criteria for molecular studies of DLD to reduce case heterogeneity.

DNA-methylation and immunological response in medication overuse headache.

OBJECTIVE: To investigate whether medication-overuse headache patients have differential DNA-methylation pattern. METHODS: We collected blood samples from 120 medication-overuse headache-patients, 57 controls (29 episodic migraine patients and 28 healthy controls) in a hypothesis-generating cross-sectional case-control pilot study; 100 of the medication-overuse headache-patients were followed for six months and samples were collected at two and six months for the longitudinal methylation analyses. Blood cell proportions of leucocytes (neutrophils, NK-cells, monocytes, CD8+ and CD4+ T-cells, and B-cells) and the neutrophile-lymphocyte ratio were estimated using methylation data as a measure for immunological analysis and a cell type-specific epigenome wide association study was conducted between medication-overuse headache-patients and controls, and longitudinally for reduction in headache days/month among medication-overuse headache-patients. RESULTS: We found a higher neutrophile-lymphocyte ratio in medication-overuse headache-patients compared to controls, indicating a higher immunological response in medication-overuse headache-patients (false discovery rate (adjusted p-value)<0.001). Reduction in headache days/month (9.8; 95% CI 8.1-11.5) was associated with lower neutrophile-lymphocyte ratio (false discovery rate adjusted p-value = 0.041).Three genes (CORIN, CCKBR and CLDN9) were hypermethylated in specific cell types in medication-overuse headache-patients compared to controls. No methylation differences were associated with reduction in headache days in medication-overuse headache-patients after six months. CONCLUSION: This pilot study was consistent with higher immunological response in medication-overuse headache-patients which decreased with a reduction in headache days in longitudinal analysis. medication-overuse headache-patients exhibited differential methylation in innate immune cells but did not exhibit longitudinal differences with alterations in headache days. Our study creates hypotheses for further biomarker searches.ClinicalTrials.gov Identifier: NCT02993289.