Genome-wide association studies (GWASs) have linked numerous genetic loci active in vascular cells to coronary artery disease (CAD), implicating smooth muscle cells (SMCs) and SMC-derived mesenchymal cells as potential mediators. We combined CAD GWAS with single-cell RNA sequencing (scRNA-seq) from human atherosclerotic plaques to identify 20 risk genes with putative action in SMCs, and then performed in vitro perturbation experiments in SMCs driven toward plaque-relevant phenotypes. Although the genes encode diverse proteins, their perturbations converged on shared transcriptional programs regulating contractile machinery, cell-cycle progression, nuclear factor κB (NF-κB), and type I interferon signaling. Integrating GWAS effect-direction with cholesterol- and stretch-responsive gene modules suggest that cholesterol-induced signaling promotes pro-atherogenic SMC states and is differentially modulated by risk versus protective variants. These results delineate polygenic regulation of SMC disease mechanisms and show that GWAS effect-directionality can help prioritize cellular pathways for follow-up functional studies across genes and cell types.

Abstract Background The presence of coronary plaques with high-risk characteristics is strongly associated with adverse cardiac events beyond the identification of coronary stenosis. Testing by coronary computed tomography angiography (CCTA) enables the identification of high-risk plaques (HRP). Referral for CCTA is presently based on pre-test probability estimates including clinical risk factors (CRFs); however, proteomics and/or genetic information could potentially improve patient selection for CCTA and, hence, identification of HRP. We aimed to (1) identify proteomic and genetic features associated with HRP presence and (2) investigate the effect of combining CRFs, proteomics, and genetics to predict HRP presence. Methods Consecutive chest pain patients (n = 1462) undergoing CCTA to diagnose obstructive coronary artery disease (CAD) were included. Coronary plaques were assessed using a semi-automatic plaque analysis tool. Measurements of 368 circulating proteins were obtained with targeted Olink panels, and DNA genotyping was performed in all patients. Imputed genetic variants were used to compute a multi-trait multi-ancestry genome-wide polygenic score (GPSMult). HRP presence was defined as plaques with two or more high-risk characteristics (low attenuation, spotty calcification, positive remodeling, and napkin ring sign). Prediction of HRP presence was performed using the glmnet algorithm with repeated fivefold cross-validation, using CRFs, proteomics, and GPSMult as input features. Results HRPs were detected in 165 (11%) patients, and 15 input features were associated with HRP presence. Prediction of HRP presence based on CRFs yielded a mean area under the receiver operating curve (AUC) ± standard error of 73.2 ± 0.1, versus 69.0 ± 0.1 for proteomics and 60.1 ± 0.1 for GPSMult. Combining CRFs with GPSMult increased prediction accuracy (AUC 74.8 ± 0.1 (P = 0.004)), while the inclusion of proteomics provided no significant improvement to either the CRF (AUC 73.2 ± 0.1, P = 1.00) or the CRF + GPSMult (AUC 74.6 ± 0.1, P = 1.00) models, respectively. Conclusions In patients with suspected CAD, incorporating genetic data with either clinical or proteomic data improves the prediction of high-risk plaque presence. Trial registration https://clinicaltrials.gov/ct2/show/NCT02264717 (September 2014).

Drugs that lower plasma apolipoprotein B (ApoB)-containing lipoproteins are central to treating advanced atherosclerosis and provide partial protection against clinical events. Previous research showed that lowering ApoB-containing lipoproteins stops plaque inflammation, but how these drugs affect the heterogeneous population of plaque cells derived from smooth muscle cells (SMCs) is unknown. SMC-derived cells are the main cellular component of atherosclerotic lesions and the source of structural components that determine the size of plaques and their propensity to rupture and trigger thrombosis, the proximate cause of heart attack and stroke. Using lineage tracing and single-cell techniques to investigate the full SMC-derived cellular compartment in progressing and regressing plaques in mice, here we show that lowering ApoB-containing lipoproteins reduces nuclear factor kappa-light-chain-enhancer of activated B cells signaling in SMC-derived fibromyocytes and chondromyocytes and leads to depletion of these abundant cell types from plaques. These results uncover an important mechanism through which cholesterol-lowering drugs can achieve plaque regression.

Background: Patients with de novo chest pain, referred for evaluation of possible coronary artery disease (CAD), frequently have an absence of CAD resulting in millions of tests not having any clinical impact. The objective of this study was to investigate whether polygenic risk scores and targeted proteomics improve the prediction of absence of CAD in patients with suspected CAD, when added to the PROMISE (Prospective Multicenter Imaging Study for Evaluation of Chest Pain) minimal risk score (PMRS). Methods: Genotyping and targeted plasma proteomics (N=368 proteins) were performed in 1440 patients with symptoms suspected to be caused by CAD undergoing coronary computed tomography angiography. Based on individual genotypes, a polygenic risk score for CAD (PRS CAD ) was calculated. The prediction was performed using combinations of PRS CAD , proteins, and PMRS as features in models using stability selection and machine learning. Results: Prediction of absence of CAD yielded an area under the curve of PRS CAD -model, 0.64±0.03; proteomic-model, 0.58±0.03; and PMRS model, 0.76±0.02. No significant correlation was found between the genetic and proteomic risk scores (Pearson correlation coefficient, −0.04; P =0.13). Optimal predictive ability was achieved by the full model (PRS CAD +protein+PMRS) yielding an area under the curve of 0.80±0.02 for absence of CAD, significantly better than the PMRS model alone ( P <0.001). For reclassification purpose, the full model enabled down-classification of 49% (324 of 661) of the 5% to 15% pretest probability patients and 18% (113 of 611) of >15% pretest probability patients. Conclusions: For patients with chest pain and low-intermediate CAD risk, incorporating targeted proteomics and polygenic risk scores into the risk assessment substantially improved the ability to predict the absence of CAD. Genetics and proteomics seem to add complementary information to the clinical risk factors and improve risk stratification in this large patient group. REGISTRATION: URL: https://www.clinicaltrials.gov ; Unique identifier: NCT02264717

Introduction Current guideline recommend functional imaging for myocardial ischaemia if coronary CT angiography (CTA) has shown coronary artery disease (CAD) of uncertain functional significance. However, diagnostic accuracy of selective myocardial perfusion imaging after coronary CTA is currently unclear. The Danish study of Non-Invasive testing in Coronary Artery Disease 3 trial is designed to evaluate head to head the diagnostic accuracy of myocardial perfusion imaging with positron emission tomography (PET) using the tracers 82Rubidium (82Rb-PET) compared with oxygen-15 labelled water PET (15O-water-PET) in patients with symptoms of obstructive CAD and a coronary CT scan with suspected obstructive CAD.Methods and analysis This prospective, multicentre, cross-sectional study will include approximately 1000 symptomatic patients without previous CAD. Patients are included after referral to coronary CTA. All patients undergo a structured interview and blood is sampled for genetic and proteomic analysis and a coronary CTA. Patients with possible obstructive CAD at coronary CTA are examined with both 82Rb-PET, 15O-water-PET and invasive coronary angiography with three-vessel fractional flow reserve and thermodilution measurements of coronary flow reserve. After enrolment, patients are followed with Seattle Angina Questionnaires and follow-up PET scans in patients with an initially abnormal PET scan and for cardiovascular events in 10 years.Ethics and dissemination Ethical approval was obtained from Danish regional committee on health research ethics. Written informed consent will be provided by all study participants. Results of this study will be disseminated via articles in international peer-reviewed journal.Trial registration number NCT04707859.

Endometriosis is a common condition associated with debilitating pelvic pain and infertility. A genome-wide association study meta-analysis, including 60,674 cases and 701,926 controls of European and East Asian descent, identified 42 genome-wide significant loci comprising 49 distinct association signals. Effect sizes were largest for stage 3/4 disease, driven by ovarian endometriosis. Identified signals explained up to 5.01% of disease variance and regulated expression or methylation of genes in endometrium and blood, many of which were associated with pain perception/maintenance (SRP14/BMF, GDAP1, MLLT10, BSN and NGF). We observed significant genetic correlations between endometriosis and 11 pain conditions, including migraine, back and multisite chronic pain (MCP), as well as inflammatory conditions, including asthma and osteoarthritis. Multitrait genetic analyses identified substantial sharing of variants associated with endometriosis and MCP/migraine. Targeted investigations of genetically regulated mechanisms shared between endometriosis and other pain conditions are needed to aid the development of new treatments and facilitate early symptomatic intervention. Meta-analyses of genome-wide association studies for endometriosis identify 49 distinct association signals. Fine-mapping of causal variants explores functional effects across various tissues. Genetic correlations between endometriosis and other pain conditions are also highlighted.

Objectives Allele counts of sequence variants obtained by whole genome sequencing (WGS) often play a central role in interpreting the results of genetic and genomic research. However, such variant counts are not readily available for individuals in the Danish population. Here, we present a dataset with allele counts for sequence variants (single nucleotide variants (SNVs) and indels) identified from WGS of 8,671 (5,418 females) individuals from the Danish population. The data resource is based on WGS data from three independent research projects aimed at assessing genetic risk factors for cardiovascular, psychiatric, and headache disorders. To enable the sharing of information on sequence variation in Danish individuals, we created summarized statistics on allele counts from anonymized data and made them available through the European Genome-phenome Archive (EGA, https://identifiers.org/ega.dataset:EGAD0 0001009756) and in a dedicated browser, DanMAC5 (available at www.danmac5.dk). The summary level data and the DanMAC5 browser provide insight into the allelic spectrum of sequence variants segregating in the Danish population, which is important in variant interpretation.Data description Three WGS datasets with an average coverage of 30x were processed independently using the same quality control pipeline. Subsequently, we summarized, filtered, and merged allele counts to create a high -quality summary level dataset of sequence variants.

BACKGROUND AND PURPOSE:Understanding migraine in a sex-specific manner is crucial for improving clinical care, diagnosis and therapy for both females and males. Here, data on sex differences are provided in the presentation of migraine in a large European-based population cohort, which is representative of the general population.METHODS:A population-based study of 62,672 Danish blood donors (both present and previous donors), of whom 12,658 had migraine, was performed. All participants completed a 105-item diagnostic migraine questionnaire sent via an electronic mailing system (e-Boks) between May 2020 and August 2020. The questionnaire allowed for correct diagnosis of migraine according to the International Classification of Headache Disorders, third edition.RESULTS:The migraine questionnaire was in-cohort validated and had a positive predictive value of 97% for any migraine, a specificity of 93% and a sensitivity of 93%. There were 9184 females (mean age 45.1 years) and 3434 males (mean age 48.0 years). The 3-month prevalence of migraine without aura was 11% in females and 3.59% in males. The 3-month prevalence of migraine with aura was 1.72% in females and 1.58% in males. In females, the age-related 3-month prevalence of migraine without aura increased markedly during childbearing age. In males, migraine both with and without aura showed less age variation. Females had a higher frequency of migraine attacks (odds ratio [OR] 1.22) but a lower frequency of non-migraine headaches (OR = 0.35). Females also had a greater intensity of pain, more unilateral and pulsatile pain, and exacerbation by physical activity (OR = 1.40-1.49) as well as more associated symptoms (OR = 1.26-1.98). Females carried 79% of the total migraine disease burden, which was almost exclusively driven by migraine without aura (77%), whilst there was no sex difference in the disease burden of migraine with aura.CONCLUSION:Females have more severe disease, resulting in a much higher migraine disease burden than indicated by prevalence alone.

Background: Migraine is the second leading cause of disability worldwide and affects two to three times more females than males. Other differences between sexes are, however, not well-documented. The aim of the present study was to assess differences between sexes in the clinical burden of migraine, i.e. attack frequency, severity, and number of non-headache symptoms, i.e. accompanying symptoms, and to assess if sex-specific genetic risk factors provide early insight into such differences.Methods: A population-based study of 62,672 Danes. All participants completed a 105-item diagnostic migraine questionnaire, of whom 12,658 had migraine. The questionnaire had a positive predictive value of 97% for any migraine. Participants were genotyped using a genome-wide array and imputation was done using a Danish reference genome. Results were replicated in a large clinical material.Findings: The age-related prevalence showed a marked increase during childbearing age for migraine without aura. Females also had a higher headache frequency, longer duration of migraine attacks, greater intensity of pain, more unilateral and pulsatile pain, exacerbation by physical activity, nausea, vomiting, photophobia, phonophobia, osmophobia, and allodynia than males. Females had lower effect of over-the-counter analgesics than males. We found five genome-wide significant loci for females, of which FHL5/UFL1 replicated in females with migraine (OR=1.27, 95%CI [1.12-1.45], P<0.001) and STAT6/LRP1 replicated in females with migraine without aura (OR=1.33, 95%CI [1.08-1.14], P<0.001). Secondary analyses showed a pathway enrichment of nerve growth factor signaling in females with migraine without aura (P<0.001) and glucose metabolism in males with migraine with aura (P<0.001).Interpretation: The clinical burden of migraine in females was considerably more severe than reflected by prevalence only. Genetic findings gave early insight into the sex-specific mechanisms of migraine. The significantly higher severity of migraine in females demands increased focus on sex-specific mechanisms and treatment.Funding Candys Foundation and The Research Funding Pool at Rigshospitalet.Funding Information: The study was financed by a grant from Candys Foundation, the CEHEAD program (J.O.) and The Research Funding Pool at Rigshospitalet (M.A.C.). Declaration of Interests: The authors have no potential conflicts of interest with respect to the research, authorship, and/or publication of this study.Ethics Approval Statement: Written informed consent was obtained from all participants. The DBDS study is an on-going, national study and was approved by the Danish Ethical Standards Committees in the relevant regions of Denmark (DESC) (1-10-72-95-13, SJ-740, 1-90-09-88 and 1-70-04-07) and the Danish Data Protection Agency (DDPA) (P-2019-99). Studies from the Danish Headache Center were approved by the DESC (H-2-2010-122) and the DDPA (01080/GLO-2010-10).

BackgroundThe cause of atrioventricular block (AVB) remains unknown in approximately half of young patients with the diagnosis. Although variants in several genes associated with cardiac conduction diseases have been identified, the contribution of genetic variants in younger patients with AVB is unknown.Methods and ResultsUsing the Danish Pacemaker and Implantable Cardioverter Defibrillator (ICD) Registry, we identified all patients younger than 50 years receiving a pacemaker because of AVB in Denmark in the period from January 1, 1996 to December 31, 2015. From medical records, we identified patients with unknown cause of AVB at time of pacemaker implantation. These patients were invited to a genetic screening using a panel of 102 genes associated with inherited cardiac diseases. We identified 471 living patients with AVB of unknown cause, of whom 226 (48%) accepted participation. Median age at the time of pacemaker implantation was 39 years (interquartile range, 32–45 years), and 123 (54%) were men. We found pathogenic or likely pathogenic variants in genes associated with or possibly associated with AVB in 12 patients (5%). Most variants were found in theLMNAgene (n=5).LMNAvariant carriers all had a family history of either AVB and/or sudden cardiac death.ConclusionsIn young patients with AVB of unknown cause, we found a possible genetic cause in 1 out of 20 participating patients. Variants in theLMNAgene were most common and associated with a family history of AVB and/or sudden cardiac death, suggesting that genetic testing should be a part of the diagnostic workup in these patients to stratify risk and screen family members.

Brain homeostasis depends on the existence of the blood-brain barrier (BBB). Despite decades of research, the factors and signalling pathways for modulating and maintaining BBB integrity are not fully elucidated. Here, we characterise the expression and function of the multifunctional receptor, sortilin, in the cells of the BBB, in vivo and in vitro. We show that sortilin acts as an important regulatory protein of the BBB's tightness. In rats lacking sortilin, the BBB was leaky, which correlated well with relocated distribution of the localisation of zonula occludens-1, VE-cadherin and β-catenin junctional proteins. Furthermore, the absence of sortilin in brain endothelial cells resulted in decreased phosphorylation of Akt signalling protein and increased the level of phospho-ERK1/2. As a putative result of MAPK/ERK pathway activity, the junctions between the brain endothelial cells were disintegrated and the integrity of the BBB became compromised. The identified barrier differences between wild-type and Sort1-/- brain endothelial cells can pave the way for a better understanding of sortilin's role in the healthy and diseased BBB.

Background: Polygenic risk scores (PRSs) are associated with coronary artery disease (CAD), but the clinical potential of using PRSs at the single-patient level for risk stratification has yet to be established. We investigated whether adding a PRS to clinical risk factors (CRFs) improves risk stratification in patients referred to coronary computed tomography angiography on a suspicion of obstructive CAD. Methods: In this prespecified diagnostic substudy of the Dan-NICAD trial (Danish study of Non-Invasive testing in Coronary Artery Disease), we included 1617 consecutive patients with stable chest symptoms and no history of CAD referred for coronary computed tomography angiography. CRFs used for risk stratification were age, sex, symptoms, prior or active smoking, antihypertensive treatment, lipid-lowering treatment, and diabetes. In addition, patients were genotyped, and their PRSs were calculated. All patients underwent coronary computed tomography angiography. Patients with a suspected ≥50% stenosis also underwent invasive coronary angiography with fractional flow reserve. A combined end point of obstructive CAD was defined as a visual invasive coronary angiography stenosis >90%, fractional flow reserve <0.80, or a quantitative coronary analysis stenosis >50% if fractional flow reserve measurements were not feasible. Results: The PRS was associated with obstructive CAD independent of CRFs (adjusted odds ratio, 1.8 [95% CI, 1.5–2.2] per SD). The PRS had an area under the curve of 0.63 (0.59–0.68), which was similar to that for age and sex. Combining the PRS with CRFs led to a CRF+PRS model with area under the curve of 0.75 (0.71–0.79), which was 0.04 more than the CRF model ( P =0.0029). By using pretest probability (pretest probability) cutoffs at 5% and 15%, a net reclassification improvement of 15.8% ( P =3.1×10 −4 ) was obtained, with a down-classification of risk in 24% of patients (211 of 862) in whom the pretest probability was 5% to 15% based on CRFs alone. Conclusions: Adding a PRS improved risk stratification of obstructive CAD beyond CRFs, suggesting a modest clinical potential of using PRSs to guide diagnostic testing in the contemporary clinical setting. Registration: URL: https://www.clinicaltrials.gov ; Unique identifier: NCT02264717.

Genetic variants in the genomic region containing SORT1 (encoding the protein sortilin) are strongly associated with cholesterol levels and the risk of coronary artery disease (CAD). Circulating sortilin has therefore been proposed as a potential biomarker for cardiovascular disease. Multiple studies have reported association between plasma sortilin levels and cardiovascular outcomes. However, the findings are not consistent across studies, and most studies have small sample sizes. The aim of this study was to evaluate sortilin as a biomarker for CAD in a well-characterized cohort with symptoms suggestive of CAD. In total, we enrolled 1,173 patients with suspected stable CAD referred to coronary computed tomography angiography. Sortilin was measured in plasma using two different technologies for quantifying circulating sortilin: a custom-made enzyme-linked immunosorbent assay (ELISA) and OLINK Cardiovascular Panel II. We found a relative poor correlation between the two methods (correlation coefficient = 0.21). In addition, genotyping and whole-genome sequencing were performed on all patients. By whole-genome regression analysis of sortilin levels measured with ELISA and OLINK, two independent cis protein quantitative trait loci (pQTL) on chromosome 1p13.3 were identified, with one of them being a well-established risk locus for CAD. Incorporating rare genetic variants from whole-genome sequence data did not identify any additional pQTLs for plasma sortilin. None of the traditional CAD risk factors, such as sex, age, smoking, and statin use, were associated with plasma sortilin levels. Furthermore, there was no association between circulating sortilin levels and coronary artery calcium score (CACS) or disease severity. Sortilin did not improve discrimination of obstructive CAD, when added to a clinical pretest probability (PTP) model for CAD. Overall, our results indicate that studies using different methodologies for measuring circulating sortilin should be compared with caution. In conclusion, the well-known SORT1 risk locus for CAD is linked to lower sortilin levels in circulation, measured with ELISA; however, the effect sizes are too small for sortilin to be a useful biomarker for CAD in a clinical setting of low- to intermediate-risk chest-pain patients.

ABSTRACT Migraine has a heritability of up to 65%. Genome-wide association studies (GWAS) on migraine have identified 123 risk loci, explaining only 10.6% of migraine heritability. Thus, there is a considerable genetic component not identified with GWAS. Further, the causality of the identified risk loci remains inconclusive. Rare variants contribute to the risk of migraine but GWAS are often underpowered to detect these. Whole genome sequencing is reliable for analyzing rare variants but is not frequently used in large-scale. We assessed if rare variants in the migraine risk loci associated with migraine. We used a large cohort of whole genome sequenced migraine patients (1,040 individuals from 155 families). The findings were replicated in an independent case-control cohort (2,027 migraine patients, 1,650 controls). We found rare variants (minor allele frequency<0.1%) associated with migraine in a Polycomb Response Element in the ASTN2 locus. The association was independent of the GWAS lead risk variant in the locus. The findings place rare variants as risk factors for migraine. We propose a biological mechanism by which epigenetic regulation by Polycomb Response Elements plays a crucial role in migraine etiology.

SORL1 is strongly associated with both sporadic and familial forms of Alzheimer’s disease (AD), but a lack of information about alternatively spliced transcripts currently limits our understanding of the role of SORL1 in AD. Here, we describe a SORL1 transcript ( SORL1-38b ) characterized by inclusion of a novel exon (E38b) that encodes a truncated protein. We identified E38b-containing transcripts in several brain regions, with the highest expression in the cerebellum and showed that SORL1-38b is largely located in neuronal dendrites, which is in contrast to the somatic distribution of transcripts encoding the full-length SORLA protein ( SORL1-fl ). SORL1-38b transcript levels were significantly reduced in AD cerebellum in three independent cohorts of postmortem brains, whereas no changes were observed for SORL1-fl . A trend of lower 38b transcript level in cerebellum was found for individuals carrying the risk variant at rs2282649 (known as SNP24), although not reaching statistical significance. These findings suggest synaptic functions for SORL1-38b in the brain, uncovering novel aspects of SORL1 that can be further explored in AD research.

A major challenge to long read sequencing data is their high error rate of up to 15%. We present Ratatosk, a method to correct long reads with short read data. We demonstrate on 5 human genome trios that Ratatosk reduces the error rate of long reads 6-fold on average with a median error rate as low as 0.22 %. SNP calls in Ratatosk corrected reads are nearly 99 % accurate and indel calls accuracy is increased by up to 37 %. An assembly of Ratatosk corrected reads from an Ashkenazi individual yields a contig N50 of 45 Mbp and less misassemblies than a PacBio HiFi reads assembly.

Hematological malignancies are an aggregate of diverse populations of cells that arise following a complex process of clonal evolution and selection. Recent approaches have facilitated the study of clonal populations and their evolution over time across multiple phenotypic cell populations. In this review, we present current concepts on the role of clonal evolution in leukemic initiation, disease progression, and relapse. We highlight recent advances and unanswered questions about the contribution of the hematopoietic stem cell population to these processes.

Background The development of long read sequencing (LRS) has led to greater access to the human genome. LRS produces long read lengths at the cost of high error rates and has shown to be more useful in calling structural variants than short read sequencing (SRS) data. In this paper we evaluate how to use LRS data from Oxford Nanopore Technologies (ONT) to call small variants in regions in- and outside the reach of SRS. Results Calling single nucleotide polymorphisms (SNPs) with ONT data has comparable accuracy to Illumina when evaluating against the Genome in a Bottle truth set v4.2. In the major histocompatibility complex (MHC) and regions where mapping short reads is difficult, the F-measure of ONT calls exceeds those of short reads by 2-4% when sequence coverage is 20X or greater. We develop recommendations for how to perform small variant calling with LRS data and improve current approaches to the difficult regions by re-genotyping variants to increase the F-measure from 97.24% to 98.78%. Furthermore, we show how LRS can call variants in genomic regions inaccessible to SRS, including medically relevant genes such as STRC and CFC1B . Conclusions Although small variant calling in LRS data is still immature, current methods are clearly useful in difficult and inaccessible regions of the genome, enabling variant calling in medically relevant genes not accessible to SRS. ### Competing Interest Statement Guillaume Holley, Doruk Beyler and Bjarni V. Halldórsson are employees of deCODE genetics/Amgen Inc.

The most recent genome-wide association study of migraine increased the total number of known migraine risk loci to 38. Still, most of the heritability of migraine remains unexplained, and it has been suggested that rare gene dysregulatory variants play an important role in migraine etiology. Addressing the missing heritability of migraine, we aim to fine-map signals from the known migraine risk loci to regulatory mechanisms and associate these to downstream genic targets. We analyzed a large cohort of whole-genome sequenced patients from extended migraine pedigrees (1040 individuals from 155 families). We test for association between rare variants segregating in regulatory regions with migraine. The findings were replicated in an independent case-control cohort (2027 migraineurs, 1650 controls). We report an increased burden of rare variants in one CpG island and three polycomb group response elements near four migraine risk loci. We found that the association is independent of the common risk variants in the loci. The regulatory regions are suggested to affect different genes than those originally tagged by the index SNPs of the migraine loci. Families with familial clustering of migraine have an increased burden of rare variants in regulatory regions near known migraine risk loci, with effects that are independent of the variants in the loci. The possible regulatory targets suggest different genes than those originally tagged by the index SNPs of the migraine loci.