- Open Access
A study of differential microRNA expression profile in migraine: the microMIG exploratory study
The Journal of Headache and Pain volume 24, Article number: 11 (2023)
Several studies have described potential microRNA (miRNA) biomarkers associated with migraine, but studies are scarcely reproducible primarily due to the heterogeneous variability of participants. Increasing evidence shows that disease-related intrinsic factors together with lifestyle (environmental factors), influence epigenetic mechanisms and in turn, diseases. Hence, the main objective of this exploratory study was to find differentially expressed miRNAs (DE miRNA) in peripheral blood mononuclear cells (PBMC) of patients with migraine compared to healthy controls in a well-controlled homogeneous cohort of non-menopausal women.
Patients diagnosed with migraine according to the International Classification of Headache Disorders (ICHD-3) and healthy controls without familial history of headache disorders were recruited. All participants completed a very thorough questionnaire and structured-interview in order to control for environmental factors. RNA was extracted from PBMC and a microarray system (GeneChip miRNA 4.1 Array chip, Affymetrix) was used to determine the miRNA profiles between study groups. Principal components analysis and hierarchical clustering analysis were performed to study samples distribution and random forest (RF) algorithms were computed for the classification task. To evaluate the stability of the results and the prediction error rate, a bootstrap (.632 + rule) was run through all the procedure. Finally, a functional enrichment analysis of selected targets was computed through protein–protein interaction networks.
After RF classification, three DE miRNA distinguished study groups in a very homogeneous female cohort, controlled by factors such as demographics (age and BMI), life-habits (physical activity, caffeine and alcohol consumptions), comorbidities and clinical features associated to the disease: miR-342-3p, miR-532-3p and miR-758-5p. Sixty-eight target genes were predicted which were linked mainly to enriched ion channels and signaling pathways, neurotransmitter and hormone homeostasis, infectious diseases and circadian entrainment.
A 3-miRNA (miR-342-3p, miR-532-3p and miR-758-5p) novel signature has been found differentially expressed between controls and patients with migraine. Enrichment analysis showed that these pathways are closely associated with known migraine pathophysiology, which could lead to the first reliable epigenetic biomarker set. Further studies should be performed to validate these findings in a larger and more heterogeneous sample.
Migraine is a genetically-driven, chronic brain disease that ranks second on the most disabling neurological disorders of the Global Burden of Disease Study [1, 2]. If not adequately diagnosed and treated, migraine can increase in frequency and transform, through a chronification process, into chronic migraine (CM), which is defined by having 15 or more headache days per month , where patients’ burden increases exponentially . Multiple factors, including both genetic and environmental mechanisms seem to play a role in migraine and its chronification [5,6,7] process, although the susceptibility of headache attack recurrence mainly depends on the brain-sensorial adaptability of both internal and external stimuli [8, 9]. Moreover, migraine stills represent a clinical challenge because, up to date, disease diagnosis, evolution and treatment response cannot be predicted due to the lack of reliable biomarkers.
From large-scale genome-wide association studies (GWAS) performed in the last decades, we have learnt that migraine has a polygenic heritability pattern and, up to now, we have identified more than 120 risk loci associated to migraine , mainly codifying for neuronal and vascular functions. However, the effect size of each individual reported gene is relatively modest and most of these genes have regulatory effects on gene expression rather than protein coding. Environmental factors may also modulate the expression of DNA. Therefore, the study of epigenetics which are the environmental influenced genetic regulatory mechanisms  could provide further understanding on migraine pathophysiology as well as, give insight to the dynamic aspect of the disease specifically contributing to its chronification .
Among the epigenetic regulatory mechanisms, microRNAs (miRNA) are a subtype of noncoding RNA of ∼22 nucleotides long that act as key regulators of genetic expression by inhibiting transcription or promoting degradation of selected messenger RNAs (mRNAs). A single miRNA sequence could regulate multiple mRNA, exerting a pleiotropic modulation of cellular processes. miRNA have been found intracellularly and in all body fluids, including plasma, saliva and urine . Accessing the central nervous system (CNS) is too invasive and there is an increasing evidence that a relationship between peripheral blood mononuclear cells (PBMC) and brain-based epigenetic activity exists . Thus, PBMC gene expression may reflect both central and peripheral neurological mechanisms in migraine, induced, for example in migraine by the activation and neuroinflammation of the trigeminocervical pain pathway . In fact, previous studies in headache disorders have shown that transcriptional changes of PBMC can reflect alterations in the sensitization of the peripheral and/or central pain pathway [16, 17].
Several pilot studies have described miRNAs as potential peripheral biomarkers of migraine disease [18,19,20,21,22,23,24,25,26]. Nevertheless, studies are scarcely comparable and reproducible not only because of the lack of standardized protocols, biological matrix or sample size but mainly because of the heterogeneous population and clinical variability included. Increasing evidence shows that disease-related intrinsic factors altogether with lifestyle factors may influence epigenetic mechanisms and identifying them is crucial to find unbiased reliable factors associated with a specific disease [27,28,29].
Hence, the goal of this exploratory study is to find differentially expressed miRNAs (DE miRNA) in PBMC of patients with migraine compared to healthy controls in a well-controlled homogeneous non-menopausal female cohort.
Participants and protocol approval
The study was approved by our Institutional Review Board (PR(IR)224/2016). All participants provided written informed consent before participation. We recruited patients treated for migraine in the Headache Outpatient Clinic. Patients had a diagnosis of migraine with or without aura, confirmed by a neurologist according to the International Classification of Headache Disorders, 3rd edition (ICHD-3) . Healthy controls (HC) were recruited among hospital staff and non-related acquaintances of patients who denied past or first-degree familial history of any recurrent primary or secondary headache disorder. To diminish confounding factors, the inclusion criteria were: female sex, 18–45 years old (non-menopausal) and Mediterranean ethnicity for all participants. Moreover, patients could not have received prophylactic treatment for migraine prevention or if they received it, it had been done for a brief period of time more than one year or more before the study and without following national guidelines of effective dosage and duration. Due to the strict inclusion criteria, participants were included from August 2016 to September 2017.
As miRNAs can be modulated by different lifestyle habits and disease characteristics, recruited participants completed a very thorough and complete questionnaire and structured-interview to collect demographic information including: age, sex, ethnicity, educational level, marital status, gross income; clinical information including: habits of caffeine, tobacco and drugs consumption, physical activity level using the International Physical Activity Questionnaire (IPAQ-SF) ; life-time comorbidities, medications use, menstrual status and migraine information including disease duration, acute medication use, monthly headache frequency and family history of headache. We registered headache pain, accompanying symptoms and menstrual cycle the day of the blood sample extraction. We only analyzed samples from migraine patients with at least 48 h of a headache-free period prior day of the blood extraction. In addition, all patients had to fill in the following scales: Allodynia Symptom Checklist (ASC-12) , Headache Impact Test-6 (HIT-6) , Migraine Disability Assessment Test (MIDAS) , Beck Depression Inventory 2nd Edition (BDI-II) , State-Trait Anxiety Inventory (STAI)  and Perceived stress scale (PSS) . Patients were asked to complete headache diary 30 days before extraction to confirm migraine diagnosis (screening period). We also excluded patients who fulfilled acute medication overuse headache criteria during the screening period.
Blood extraction, PBMC isolation and RNA extraction
For blood extraction patients were asked to fast for 6 h, minimize their exercise/activity and not take analgesic medication 24 h prior to blood extraction. Peripheral blood was collected by antebrachial vein puncture using three, 8 mL, BD Vacutainer® CPT™ Cell Preparation Tubes with Sodium Citrate (BD, Franklin Lakes, NJ, USA) and processed within 2 h following the manufacturer’s instruction to isolate PBMC pellets that were immediately stored at -80ºC till miRNA extraction. To minimize variability due to manipulation, the miRNA extraction was made by the same laboratory technician in the same facilities. For miRNA extraction, we used the Qiagen miRNeasy Mini Kit (Qiagen, Valencia, CA, USA), following the manufacturer’s instruction. RNA concentration was measured using Nanodrop 2000 (NanoDrop Technologies, Wilmington, DE). RNA integrity number (RIN) was determined the kit Agilent RNA 6000 Nano by the Bioanalyzer Agilent 2100 (Agilent Technologies, Santa Clara, CA, USA). We excluded from analysis samples with RIN under 6.
Microarrays service was outsourced and carried out in at High Technology Unit at Vall d’Hebron Research Institute (HTU-VHIR), Barcelona (Spain). For this study, we used the Affymetrix GeneTitan microarray platform and the Genechip miRNA array plate 4.1. Starting material was 550 ng of total RNA of each sample. MiRNA in the sample was labelled using Flash Tag Biotin HSR RNA Labelling kit following the manufacturer’s instructions.
All statistical analysis were conducted in R v4.2.0  and figures were produced using the package ggplot2 v.3.3.6 . A statistical power calculation was not conducted prior to the study because the sample size was based on the available data for this exploratory analysis. However, to increase the likelihood of finding differential miRNA expression profiles for this study, we selected participants from a well-controlled homogeneous non-menopausal female cohort. P-values < 0.05 were considered as statistically significant and are reported for a two-tailed test.
Nominal variables (aura and presence of allodynia) were reported as frequencies (percentages) while median and interquartile range (IQR) were reported for quantitative variables: age, BMI, IPAQ-SF, caffeine and alcohol consumptions, STAI, BDI-II, PSS, disease evolution time, ASC-12, monthly headache frequency (MHD), monthly acute medication intake (MAMI), MIDAS and HIT-6. Normality assumption of quantitative variables was checked through visual methods (Q-Q plots) and normality tests (Shapiro–Wilk test). Statistical significance between study groups (HC vs. patients with migraine) was assessed by unpaired t-test for BMI, caffeine consumption, STAI and PSS was used and unpaired Wilcoxon rank-sum test was used for the other quantitative variables that did not follow any normality assumption.
The statistical analysis of the Genechip miRNA array plate 4.1 was performed using libraries developed for the microarray analysis in the Bioconductor Project . At the moment of the analysis, there were no specific packages for miRNA 4.1 arrays available, thus we programmed our own package. Quality control (QA) was assessed with the package arrayQualityMetrics before the normalization process . Additionally, principal components analysis (PCA) and hierarchical clustering analysis (HCA) were performed to study samples distribution and detect technical problems, including batch effect. Each array (CEL file) was preprocessed and background corrected, normalized and summarized using RMA (Robust Multiarray Average) through Bioconductor packages [41, 42].
The identification of differentially expressed miRNA (DE miRNAs) was based on adjusting linear model with empirical Bayes moderated t-statistic, conducted by the limma package  in the comparison between controls and patients with migraine. Models were adjusted by covariates using surrogate variables estimation . DE miRNAs that had more than 0.5-fold increase/decrease relative to controls were chosen for further analysis. The minimum of 0.5-fold expression difference cut-off has been chosen based on the several miRNA profiling studies using microarrays; confirming that subtle changes in miRNA expression, such as a 0.5-fold difference, have a significant biological impact [45,46,47,48]. Results were further corrected for multi-testing using the Benjamini–Hochberg procedure for control the false discovery rate (FDR) . Features that met both criteria (adjusted t-test and ± 0.5 fold change) were combined in a single matrix of DE miRNAs that underwent classification analysis.
We used Random Forests (RF) models for the classification task, a machine learning method widely used in the microarray analysis [50,51,52]. RF algorithm uses an ensemble of classification trees, internally and randomly constructed using a bootstrap sample of the data. Then, the algorithm merges all constructed tree altogether (forest) in order to get an accurate and stable prediction [53, 54]. Once we obtained the forest, the construction of the classifier is performed by a feature reduction algorithm, where the less important features are successively eliminated and out of the bag (OOB) error is continuously analyzed. To evaluate the stability of the results and the prediction error rate, a bootstrap (0.632 + rule) was run through all the procedure. All classification task and bootstrap were performed through the varSelRF R package [55, 56].
Target genes of selected DE miRNAs were predicted by using the miRDB database through the multimiR package  and functional and pathway enrichment analysis was performed using the pathfindR package . PathfindR identifies gene sets that form active subnetworks in a protein–protein interaction network (PIN) in order to identify distinct active subnetworks and then perform enrichment analyses on these subnetworks. For assigning a significance value for each predicted target gene (pathfindR input), all Homo sapiens miRNA-target gene predictions were obtained: 1,954,039 miRNA-gene predictions from all H. sapiens miRNAs listed in the miRDB (2,656 miRNAs). Then, the significance for each predicted gene was computed as the probability of observing a score greater than or equal to the score of this target gene over all H. sapiens miRNA-target gene scores . For genes that were targeted by more than one miRNA, the lowest significance was kept.
Clinical characteristics of the patients
Participant recruitment is represented in Fig. 1. Since August 2016, a total of 76 (49 patients with migraine and 27 HC) participants were recruited. For this exploratory study, RNA was extracted from 69 samples that fulfilled inclusion/exclusion criteria. Among these, 37 showed a RIN < 6 and were discarded. Therefore, 32 participants (12 HC, 20 patients) were considered for the analysis with a median age of 33.2 [28.0, 39.8] years old. No statistically significant differences were found neither in demographical variables nor life-habits measurements between groups. In terms of comorbidities, we only found that patients with migraine presented a statistically significantly greater trait-anxiety score (HC: 31.0 [27.0, 38.2] vs. migraine: 39.5 [37.2, 55.8]; p = 0.007). In patients, median MHD was 10.5 [8.0, 24.2] days/month (d/mo) and median MAMI was 7.5 [4.0, 9.0] d/mo. Main demographic and clinical characteristics are shown in Table 1.
No statistically significant differences were found in terms of RNA quality (RIN, median [IQR]) between study groups: 6.9 [6.5, 7.8] HC vs. 6.7 [6.3, 7.2] patients; p = 0.420.
MiRNA microarray expression profiling
During microarrays QA, we distinguished two samples (from patients) with higher values of intensity in all their probe sets. Similarly, normalized unscaled standard errors values and relative log expression values for the same two samples diverged from the rest, indicating aberrant expression. Hence, we discarded these 2 samples from further analysis (Fig. 1). Control and migraine samples were perfectly segregated along PCA and HCA (Fig. 2).
Adjusting by STAI-trait values, we found 191 DE miRNA between study groups: 59 over-expressed and 132 under-expressed transcripts in patients (Fig. 3A). Then, RF analysis combined with a feature selecting algorithm was used to build the classifier and to identify the best targets that distinguished between controls and migraine. The OOB for the initial model was 0.067 (91.7% sensitivity, 94.4% specificity) and a 3-miRNA signature was found as best variables without dropping the OBB error substantially: miR-342-3p, miR-532-3p and miR-758-5p. Then, as RF algorithms are deterministic, performing it on different training samples is the only way to generate diversity in the selection. To this end, we performed bootstrapping (+ 0.632) in which the random forest constructed for a certain number of variables was subsampled and compared. The prediction error rate among the bootstrap samples was 0.142. The same 3-miRNA signature were consistently selected (stability) in all the sub-samples, with miR-532-3p being selected in the 60% of cases. Figure 3B shows the heatmap from the selected DE miRNA between groups where the expression from all of them were under-expressed in patients: miR-342-3p (log2FC: -0.551), miR-532-3p (log2FC: -1.109) and miR-758-5p (log2FC: -0.508).
Targets prediction and functional enrichment analysis
Predicted target genes associated with the 3-miRNA signature were retrieved from the miRBD database through the multimiR package. Each target significance was computed according to the prediction score retrieved from the miRDB (see Methods section). Then, a total of 68 target genes were identified: 42 genes from miR-342-3p, 26 genes from miR-532-3p and 9 genes from miR-758-5p. While four common target genes were found between miR-342-3p and miR-532-3p: ATXN7 (ataxin 7), CDKL5 (cyclin dependent kinase like 5), PAPPA (Pappalysin 1) and RRM2 (Ribonucleotide Reductase Regulatory Subunit M2); two common targets were found between miR-342-3p, miR-758-5p: DHRSX (Dehydrogenase/Reductase X-Linked), NUFIP2 (Nuclear FMR1 Interacting Protein 2); and 3 common targets between miR-532-3p and miR-758-5p: DSG3 (Desmoglein 3), PTP4A1 (Protein Tyrosine Phosphatase 4A1), TET3 (Tet Methylcytosine Dioxygenase 3). No common target gene was found for all three miRNAs.
Functional enrichment analysis of miRNA-target genes predicted from the miRDB was performed through active-subnetwork-oriented analysis via pathfindR and a total of 24 enriched KEGG pathways were identified (Table 2). Top pathway was calcium signaling pathway and a total of 9 different target genes were involve in all significantly enriched pathways: CACNA1C, HTR2C, EP300, PTGER4, ID4, RICTOR and EEA1 from miR-342-3p; PRKCA and DTX4 from miR-532-3p and ID2 from miR-758-5p.
In this exploratory study, we found 191 differentially expressed miRNAs (59 over-expressed and 132 under-expressed transcripts) between controls and patients with migraine in a very homogeneous female cohort, controlled by factors as demographics, life-habits and comorbidities (Table 1, Fig. 3A). After transcripts classification task, we finally selected three DE miRNAs that distinguished study groups: miR-342-3p, miR-532-3p and miR-758-5p. All three DE miRNAs were under-expressed in patients compared to controls (Fig. 3B). Sixty-eight target genes were predicted from this 3-miRNA signature and in the functional analysis, we found that they mainly enriched ion channels and signaling pathways (calcium, notch, TRP, TGF-beta, GnRH, HIF-1, thyroid, oxytocin, mTOR and Wnt), neurotransmitter and hormone homeostasis (aldosterone, serotonin, growth hormone, dopamine and adrenaline), infectious diseases (Hepatitis, Influenza and tuberculosis) and circadian entrainment (Table 2).
MiRNAs may have the potential to accurately diagnose disorders at an early stage. For instance, it has been found that levels of miR-182 in the peripheral blood of patients with gliomas increase with disease grade, which had allowed to researchers to patent a diagnostic kit for early diagnosis of human glioma (patent number: CN101792793A) . This is beneficial for patients in order to receive earlier and more individual targeted treatment therapies, a milestone in migraine since we already know that if it is not adequately diagnosed and treated, migraine frequency increases until disease chronification . Moreover, there is increasing evidence that the later we start treating migraine, less effective is the treatment and more difficult it becomes to reverse the disease [5, 62]. In regards to CNS diseases, specific DE miRNAs have also been identified and validated in order to diagnose multiple sclerosis (US9187785B2) , Alzheimer’s disease (EP2733220A1) [64, 65] or major depressive disorder (EP2529222A4) . MiRNAs have also been given considerable attention as a therapeutic targets/agent. In recent years, miRNA-based therapeutic approaches have been developed in CNS diseases, especially pre-clinical studies in neurodegenerative disorders , although clinical trials have not started yet . Hence, these patents and previous studies in other diseases highlight that investigating miRNA profiles in migraine will lead us not only to a better understanding of their pathogenesis, but also to help develop diagnostic kits, identify specific cluster of patients, explore disease onset, monitor migraine attacks, personalize therapeutic treatment and revert disease chronification.
Several authors have made an effort in previous studies to identify potential miRNAs as peripheral biomarkers of migraine . Numerous miRNAs have been described associated with pain-free patients with migraine in comparison to healthy controls: miR-382-5p (over-expressed in serum and PBMC) [21, 23]; miR-34a-5p (over-expressed in serum, saliva and PBMC) [21, 24]; miR-27b (over-expressed in plasma) ; miR-181a, let-7b, and miR-22 (under-expressed in PBMC) ; miR-155, miR-126, and let-7 g (over-expressed in plasma) ; miR-30 (over-expressed in PBMC) ; miR-5189-3p, miR-3613-5p, miR-99a-3p, miR-542-3p (over-expressed in PBMC) and miR-96-5p (under-expressed) . Results are heterogenous and study-dependent but not only because there are differences in lab procedures, biological matrix or sample size but mainly because of the large degree of heterogeneity between study groups. Published studies only controlled by demographic factors (age and sex) but did not consider neither participants selection nor analysis with an adjustment by covariates other factors such as comorbidities and lifestyle habits. Epigenetics have been established as a molecular link between genetic function and the environment. Hence, epigenetic marks constitute a malleable language by which internal and external factors can influence the genome and, for that reason, it is crucial to be able to separate correlative from causal changes in order to understand and quantify the effect of the epigenetic function in diseases . Nowadays, we have more evidence that lifestyle factors such as diet , smoking , alcohol [73, 74], stress , physical activity [29, 76], treatments [77,78,79] and comorbidities [80, 81] have a great impact on the dynamisms of epigenetic marks. Controlling them will allow us to identify transcripts associated with a disease that are not unmasked for other factors.
Functional analysis showed that targeted genes associated to the 3-miRNA signature enriched pathways have been previously described and associated with migraine pathophysiology. One of the most enriched targets is CACNA1C which encodes for the alpha-1 subunit of a voltage-dependent calcium channel. Calcium channels mediate the influx of calcium ions into the cell upon membrane polarization and they are encoded by different subunits of the CACNA1 gene (CACNA1A-S) [82, 83]. From migraine pathophysiology, we already know that the trigeminovascular system (TVS) can be activated by cortical spreading depression (CSD). Mutations in CACNA1A can disrupt ionic homeostasis, increasing the susceptibility to CSD and inducing migraine attacks . Moreover, a CACNA1A variant has also been associated with migraine with aura in the last migraine association study . So, a down-regulation of CACNA1C gene could be also associated with these ionic disturbances in neuronal synapses; indicating that miR-342-3p has a reliable effect on migraine pathophysiology.
Other enriched pathways codifying for vascular functions which may have and important role in migraine etiopathogenesis were the notch signaling pathway, the inflammatory mediator regulation of TRP channels and the HIF-1 signaling pathway. Notch receptors (Notch1-4) are a family of type-1 transmembrane proteins which are essential for embryonic development. Mutations in the Notch pathway genes, like NOTCH2 and JAG1, are associated with vascular abnormalities . In migraine, there is evidence that distinct polymorphisms of NOTCH4 gene significantly modify clinical characteristics of migraine phenotype . The Transient Receptor Potential (TRP) channels are a superfamily of integral membrane proteins well-described in pain conditions that functions as cation channels: TRPA, TRPC, TRPM, TRPP, TRPL, and TRPV [87,88,89]. These channels contribute to a several physiological processes ranging from thermosensation and pain to the regulation of Ca2 + levels in the endoplasmic reticulum. TRP channels are also expressed in the meningeal nociceptors and it has been proved that their activation, promotes the release of the neuropeptide calcitonin gene related peptide (CGRP) . The role of CGRP in migraine has been exhaustively described in both preclinical and clinical studies during the last decades  and the CGRP antagonism has shown to be efficacious for the treatment of migraine [92,93,94]. However, the mechanism by which CGRP is released during migraine is unclear, and for that reason, TRP channels remain a focus of interest for their potential contribution to attacks . Moreover, from our experience in clinical trials and real-world studies there are patients with migraine who do not respond to treatments, even target-driven one such as anti-CGRP therapies , suggesting that CGRP is not the only mechanism driving migraine attacks [97, 98]. In regards to hypoxia-inducible factor 1 (HIF-1), it is a transcription factor that functions as a master regulator of oxygen homeostasis. Since HIF-1α is implicated in neuroprotection and inflammation inhibition, a recent study has demonstrated that roxadustat administration, an anti-anemia medication, ameliorated migraine-like behaviors and inhibited central pain sensitization in nitroglycerin-injected mice, which was mainly mediated by HIF-1α/NF-κB/inflammation pathway .
The relationship between hormones and neurotransmitters with migraine is not entirely clear, but we do know it exists and plays a fundamental role in its pathophysiology. On one hand, thyroid hormones, along with PTX3 and, insulin resistance share common mechanisms in the development of endothelial dysfunction , a vascular pathophysiological genetically driven mechanism of migraine . Moreover, differentially pattern levels have been found between controls and patients with migraine of the hypothalamic-tuberoinfundibular system (prolactin and growth hormone) . On the other hand, the gonadotropin releasing hormone (GnRH) controls the release of sexual hormones although estrogens can negatively regulate the GnRH from the hypothalamus at the time of the preovulatory surge, suggesting that estrogen could be related to the increase in migraine prevalence in relation with puberty and ovulation in women [103, 104]. Other hypothalamic hormone linked to migraine is the Oxytocin (OT). OT has been used to relieve migraine attacks through intranasal or intravenous administration  since OT receptors included in the trigeminal ganglia and neurons can inhibit the release of CGRP .
Altogether, the 3-miRNA signature found could be a reliable biomarker for migraine diagnosis not only because they are differentially expressed between female patients and controls, but they are also involved in enriched biological processes which have been previously described in migraine pathophysiology. So, the under-expression of these transcripts may cause a dysregulated expression of genes associated with the vascular dysfunction and neurohormonal imbalances characterized by patients with migraine.
This study has limitations. First, RNA-seq is the most commonly used NGS technique to explore differences at transcript level to date. RNA-seq is able to detect novel transcripts with higher specificity and sensitivity. However, at the moment of designing and planning the current study, arrays where the more standardized technology. In addition, the RNA-Seq procedure has a few disadvantages compared to arrays: a lack of optimized and standardized protocols for analysis [107, 108]; the size of RNA-Seq files is considerably larger; and RNA-seq requires highly intensive and expensive computation infrastructure and analytics, as well as, longer analysis times [108, 109]. Nevertheless, arrays technologies are still a reliable technology to study gene expression with results comparable to RNA-seq . Second, we did not evaluate the hypothesis that the miRNA profile could be related to the phenotype of circulating PBMC phenotype, rather than a differential miRNA expression. Several changes in PBMC profile related to migraine have been reported but, as different methodologies were used, the overall results are inconclusive. For instance, a meta-analysis of studies from 1966 to 1999 addressing immunological changes related to migraine did not find clear alterations in circulating leukocytes or other immunological parameters . More recently, a relative lymphocytosis has been reported in chronic migraine associated with medication overuse headache, that was considered probably related to a chronic inflammatory state and stress response rather than specifically to migraine [112, 113]. Finally, the sample was small, and highly selected and homogeneous, only considering non-menopausal women in order to reduce sample size and confounding factors. Moreover, the group of individuals diagnosed with migraine was not totally homogeneous in terms of disease pathophysiology (e.g., either episodic or chronic migraine and either the presence or absence of aura). Hence, this 3-miRNA signature should be validated in a larger sample size and more heterogeneous cohort. However, data gathered in this exploratory study will be used to calculate the required sample size for the subsequent discovery and validation (RT-qPCR) studies. Due to lack of power in small sample size experiments, true effects are often missed and many of the detected effects will probably not be further validated. Therefore, the objective of this study was to to perform an exploratory study using the same in- and exclusion criteria for a miRNA high-throughput screen .
A 3-miRNA signature (miR-342-3p, miR-532-3p and miR-758-5p) has been found differentially expressed between controls and patients with migraine in a very homogeneous cohort of non-menopausal women, controlled by factors such as demographics, life-habits and comorbidities. Enrichment analysis show that these three miRNAs are closely associated with migraine pathophysiology, maybe becoming the first reliable biomarker of disease. Further studies will be performed to validate them in a larger and more heterogeneous sample.
Allodynia Symptom Checklist
Beck Depression Inventory 2nd Edition
Central Nervous System
False Discovery Rate
Genome-Wide Association Studies
Hierarchical Clustering Analysis
Headache Impact Test-6
International Classification of Headache Disorders (3rd edition)
International Physical Activity Questionnaire (Short Form)
Monthly Acute Medication Intake
Monthly headache frequency
Migraine Disability Assessment Test
Out Of the Bag
Peripheral Blood Mononuclear Cells
Principal Components Analysis
Protein-protein Interaction Network
Perceived Stress Scale
RNA Integrity Number
Robust Multiarray Average
State-Trait Anxiety Inventory
Vall d’Hebron Research Institute
Stovner LJ, Nichols E, Steiner TJ, Abd-Allah F, Abdelalim A, Al-Raddadi RM et al (2018) Global, regional, and national burden of migraine and tension-type headache, 1990–2016: a systematic analysis for the Global Burden of Disease Study 2016. Lancet Neurol 17(11):954–976
Stovner LJ, Hagen K, Linde M, Steiner TJ (2022) The global prevalence of headache: an update, with analysis of the influences of methodological factors on prevalence estimates. J Headache Pain 23(1):1–17
Headache Classification Committee of the International Headache Society (2018) The International Classification of Headache Disorders, 3rd edition. Cephalalgia 38(1):1–211
Lanteri-Minet M (2014) Economic burden and costs of chronic migraine. Curr Pain Headache Rep 18(1):1–6
May A, Schulte LH (2016) Chronic migraine: risk factors, mechanisms and treatment. Nat Rev Neurol 12(8):455–464
Goadsby PJ, Charbit AR, Andreou AP, Akerman S, Holland PR (2009) Neurobiology of migraine. Neuroscience 161(2):327–341
Torres-Ferrús M, Ursitti F, Alpuente A, Brunello F, Chiappino D, de Vries T et al (2020) From transformation to chronification of migraine: pathophysiological and clinical aspects. J Headache Pain 21(1):1–12
Bernstein C, Burstein R (2012) Sensitization of the trigeminovascular pathway: perspective and implications to migraine pathophysiology. J Clin Neurol 8(2):89–99
Schwedt TJ (2013) Multisensory integration in migraine. Curr Opin Neurol 26(3):248
Hautakangas H, Winsvold BS, Ruotsalainen SE, Bjornsdottir G, Harder AVE, Kogelman LJA, et al (2022) Genome-wide analysis of 102,084 migraine cases identifies 123 risk loci and subtype-specific risk alleles. Nat Genet 54(2):152–60
Tammen SA, Friso S, Choi SW (2013) Epigenetics: The link between nature and nurture. Mol Aspects Med 34(4):753–64. https://doi.org/10.1016/j.mam.2012.07.018
Gallardo VJ, Vila‐Pueyo M, Pozo‐Rosich P (2023) The impact of epigenetic mechanisms in migraine : Current knowledge and future directions. Cephalalgia 43(2):3331024221145916
Ha M, Kim VN (2014) Regulation of microRNA biogenesis. Nat Rev Mol Cell Biol 15(8):509–524
Kos MZ, Puppala S, Cruz D, Neary JL, Kumar A, Dalan E et al (2022) Blood-Based miRNA Biomarkers as Correlates of Brain-Based miRNA Expression. Front Mol Neurosci 15(March):1–13
Buzzi MG, Moskowitz MA (2005) The pathophysiology of migraine: Year 2005. J Headache Pain 6(3):105–111
Gerring Z, Rodriguez-Acevedo AJ, Powell JE, Griffiths LR, Montgomery GW, Nyholt DR (2016) Blood gene expression studies in migraine: potential and caveats. Cephalalgia 36(7):669–678
Aczél T, Kun J, Szőke É, Rauch T, Junttila S, Gyenesei A et al (2018) Transcriptional alterations in the trigeminal ganglia, nucleus and peripheral blood mononuclear cells in a rat orofacial pain model. Front Mol Neurosci 11(June):1–14
Cheng CY, Chen SP, Liao YC, Fuh JL, Wang YF, Wang SJ (2018) Elevated circulating endothelial-specific microRNAs in migraine patients: a pilot study. Cephalalgia 38(9):1585–1591
Zhai Y, Zhu YY (2018) MiR-30a relieves migraine by degrading CALCA. Eur Rev Med Pharmacol Sci 22(7):2022–2028
Chen S, Chang Y, Chou C, Juan C, Lee H, Chen L et al (2021) Circulating microRNAs associated with reversible cerebral vasoconstriction syndrome. Ann Neurol 89(3):459–473
Greco R, de Icco R, Demartini C, Zanaboni AM, Tumelero E, Sances G et al (2020) Plasma levels of CGRP and expression of specific microRNAs in blood cells of episodic and chronic migraine subjects: towards the identification of a panel of peripheral biomarkers of migraine? J Headache Pain 21(1):1–12
Chen YH, Wang H (2021) The association between migraine and depression based on miRNA biomarkers and cohort studies. Curr Med Chem 28(27):5648–5656
Andersen HH, Duroux M, Gazerani P (2016) Serum microRNA signatures in migraineurs during attacks and in pain-free periods. Mol Neurobiol 53(3):1494–1500
Gallelli L, Cione E, Peltrone F, Siviglia S, Verano A, Chirchiglia D et al (2019) Hsa-miR-34a-5p and hsa-miR-375 as biomarkers for monitoring the effects of drug treatment for migraine pain in children and adolescents: a pilot study. J Clin Med 8(7):928
Tafuri E, Santovito D, de Nardis V, Marcantonio P, Paganelli C, Affaitati G et al (2015) MicroRNA profiling in migraine without aura: pilot study. Ann Med 47(6):468–473
Aczél T, Benczik B, Ágg B, Körtési T, Urbán P, Bauer W et al (2022) Disease-and headache-specific microRNA signatures and their predicted mRNA targets in peripheral blood mononuclear cells in migraineurs: role of inflammatory signalling and oxidative stress. J Headache Pain 23(1):1–20
Pérez RF, Santamarina P, Fernández AF, Fraga MF (2019) Epigenetics and lifestyle: the impact of stress, diet, and social habits on tissue homeostasis. Epigenetics Regen 11:461-89
Manokaran S, Binoy A, Bhat D, Babu S, Bhat JG, AH MR (2021) Nutrigenomics in lifestyle disorders: a review. ECS Transactions 107(1):9249
Alegría-Torres JA, Baccarelli A, Bollati V (2011) Epigenetics and lifestyle. Epigenomics 3(3):267–277
Lee PH, Macfarlane DJ, Lam TH, Stewart SM (2011) Validity of the international physical activity questionnaire short form (IPAQ-SF): A systematic review. Int J Behav Nutr Phys Act 8(1):1–11
Melhado EM, Thiers Rister HL, Galego DR, de Oliveira AB, Buttarello IA, Belucio IS et al (2020) Allodynia in menstrually related migraine: score assessment by Allodynia Symptom Checklist (ASC-12). Headache 60(1):162–70
Kosinski M, Bayliss MS, Bjorner JB, Ware JEJ, Garber WH, Batenhorst A et al (2003) A six-item short-form survey for measuring headache impact: the HIT-6. Qual Life Res 12(8):963–974
Stewart WF, Lipton RB, Dowson AJ, Sawyer J (2001) Development and testing of the Migraine Disability Assessment (MIDAS) Questionnaire to assess headache-related disability. Neurology 56(6 Suppl 1):S20–S28
Beck A, Steer R, Brown G (1996) Manual for the Beck Depression Inventory-II. In: Gellman MD, Turner JR (eds) San Antonio Psychol Corp. Springer New York, New York, NY. https://doi.org/10.1007/978-1-4419-1005-9_441
Oei TPS, Evans L, Crook GM (1990) Utility and validity of the STAI with anxiety disorder patients. Br J Clin Psychol 29(4):429–432
Sanz-Carrillo C, Garcıa-Campayo J, Rubio A, Santed MA, Montoro M (2002) Validation of the Spanish version of the Perceived Stress Questionnaire. J Psychosom Res 52(3):167–172
Team RC (2013) A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria (http://www.R-project.org/)
Wickham H (2016) ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag, New York. ISBN 978-3-319-24277-4. https://ggplot2.tidyverse.org
Gentleman RC, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2010;11:202. 2004;5(10):R80.
Kauffmann A, Gentleman R, Huber W (2009) arrayQualityMetrics - A bioconductor package for quality assessment of microarray data. Bioinformatics 25(3):415–416
Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP (2003) Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res 31(4):e15–e15
Bolstad BM, Irizarry RA, Astrand M, Speed TP (2003) A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics 19(2):185–193
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W et al (2015) limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res 43(7):e47
Leek JT, Johnson WE, Parker HS, Fertig EJ, Jaffe AE, Storey JD et al (2019) sva: Surrogate variable analysis. R package version 3:882–883
Lei B, Zhou J, Xuan X, Tian Z, Zhang M, Gao W et al (2019) Circular RNA expression profiles of peripheral blood mononuclear cells in hepatocellular carcinoma patients by sequence analysis. Cancer Med 8(4):1423–1433
Wang M, Liang L, Li L, Han K, Li Q, Peng Y et al (2017) Increased miR-424-5p expression in peripheral blood mononuclear cells from patients with pemphigus. Mol Med Rep 15(6):3479–3484
Smith-Vikos T, Liu Z, Parsons C, Gorospe M, Ferrucci L, Gill TM et al (2016) A serum miRNA profile of human longevity: findings from the Baltimore Longitudinal Study of Aging (BLSA). Aging (Albany NY) 8(11):2971
Freedman JE, Ercan B, Morin KM, Liu CT, Tamer L, Ayaz L et al (2012) The distribution of circulating microRNA and their relation to coronary disease. F1000Res 1:50
Benjamini Y, Yekutieli D (2001) The control of the false discovery rate in multiple testing under dependency. Ann Statist 29(4):1165–88
Pérez LO, González-José R, García PP (2016) Prediction of non-genotoxic carcinogenicity based on genetic profiles of short term exposure assays. Toxicol Res 32(4):289–300
Gunther EC, Stone DJ, Gerwien RW, Bento P, Heyes MP (2003) Prediction of clinical drug efficacy by classification of drug-induced genomic expression profiles in vitro. Proc Natl Acad Sci 100(16):9608–9613
Man MZ, Dyson G, Johnson K, Liao B (2004) Evaluating methods for classifying expression data. J Biopharm Stat 14(4):1065–1084
Breiman L, Friedman JH, Olshen RA, Stone CJ (2017) Classification and regression trees. Routledge. Chapman & Hall/CRC, Boca Raton. https://doi.org/10.1201/9781315139470
Breiman L (1996) Bagging predictors. Mach Learn 24(2):123–140
Díaz-Uriarte R, Alvarez de Andrés S (2006) Gene selection and classification of microarray data using random forest. BMC Bioinform 7(1):1–13
Diaz-Uriarte R (2007) GeneSrF and varSelRF: a web-based tool and R package for gene selection and classification using random forest. BMC Bioinform 8(1):1–7
Ru Y, Kechris KJ, Tabakoff B, Hoffman P, Radcliffe RA, Bowler R et al (2014) The multiMiR R package and database: Integration of microRNA-target interactions along with their disease and drug associations. Nucleic Acids Res 42(17):e133
Ulgen E, Ozisik O, Sezerman OU (2019) pathfindR: an R package for comprehensive identification of enriched pathways in omics data through active subnetworks. Front Genet 10:858
Dogan B, Gumusoglu E, Ulgen E, Sezerman OU, Gunel T (2022) Integrated bioinformatics analysis of validated and circulating miRNAs in ovarian cancer. Genomics Inform 20(2):e20
Mondal I, Kulshreshtha R (2021) Potential of microRNA based diagnostics and therapeutics in glioma: a patent review. Expert Opin Ther Pat 31(1):91–106. https://doi.org/10.1080/13543776.2021.1837775
Manack AN, Buse DC, Lipton RB (2011) Chronic migraine: Epidemiology and disease burden. Curr Pain Headache Rep 15(1):70–78
Alpuente A, Gallardo VJ, Torres-Ferrús M, Álvarez-Sabin J, Pozo-Rosich P (2020) Short and mid-term predictors of response to OnabotulinumtoxinA: real-life experience observational study. Headache 60(4):677–685
Keller A, Leidinger P, Lange J, Borries A, Schroers H, Scheffler M et al (2009) Multiple sclerosis: microRNA expression profiles accurately differentiate patients with relapsing-remitting disease from healthy controls. PLoS ONE 4(10):e7440
Geekiyanage H, Jicha GA, Nelson PT, Chan C (2012) Blood serum miRNA: non-invasive biomarkers for Alzheimer’s disease. Exp Neurol 235(2):491–496
Schipper HM, Maes OC, Chertkow HM, Wang E (2007) MicroRNA expression in Alzheimer blood mononuclear cells. Gene Regul Syst Bio. 1:GRSB-S361
Schroeter ML, Abdul-Khaliq H, Krebs M, Diefenbacher A, Blasig IE (2008) Serum markers support disease-specific glial pathology in major depression. J Affect Disord 111(2–3):271–280
Paul S, Vázquez LAB, Uribe SP, Reyes-Pérez PR, Sharma A (2020) Current status of microrna-based therapeutic approaches in neurodegenerative disorders. Cells 9(7):1–26
Chakraborty C, Sharma AR, Sharma G, Lee SS (2021) Therapeutic advances of miRNAs: a preclinical and clinical update. J Adv Res 28:127–138
Zhai Y, Zhu YY (2018) MiR-30a relieves migraine by degrading CALCA. Eur Rev Med Pharmacol Sci 22(7):2022–2028
Feil R, Fraga MF (2012) Epigenetics and the environment: emerging patterns and implications. Nat Rev Genet 13(2):97–109
Waterland RA, Jirtle RL (2003) Transposable elements: targets for early nutritional effects on epigenetic gene regulation. Mol Cell Biol 23(15):5293–5300
Lee KW, Pausova Z (2013) Cigarette smoking and DNA methylation. Front Genet 4:132 (PubMed; 2013)
Ponomarev I, Wang S, Zhang L, Harris RA, Mayfield RD (2012) Gene coexpression networks in human brain identify epigenetic modifications in alcohol dependence. J Neurosci 32(5):1884–1897
Kamat PK, Mallonee CJ, George AK, Tyagi SC, Tyagi N (2016) Homocysteine, alcoholism, and its potential epigenetic mechanism. Alcohol Clin Exp Res 40(12):2474–2481
Franklin TB, Russig H, Weiss IC, Gräff J, Linder N, Michalon A et al (2010) Epigenetic transmission of the impact of early stress across generations. Biol Psychiatry 68(5):408–415
Wallace RG, Twomey LC, Custaud MA, Turner JD, Moyna N, Cummins PM et al (2018) The role of epigenetics in cardiovascular health and ageing: A focus on physical activity and nutrition. Mech Ageing Dev 174:76–85
Menke A, Binder EB (2022) Epigenetic alterations in depression and antidepressant treatment. Dialogues Clin Neurosci 16(3):395–404
Terlizzi R, Bacalini MG, Pirazzini C, Giannini G, Pierangeli G, Garagnani P et al (2018) Epigenetic DNA methylation changes in episodic and chronic migraine. Neurol Sci 39(1):67–68
Bocchio-Chiavetto L, Maffioletti E, Bettinsoli P, Giovannini C, Bignotti S, Tardito D et al (2013) Blood microRNA changes in depressed patients during antidepressant treatment. Eur Neuropsychopharmacol 23(7):602–611
Blažeković A, Borovečki F (2021) Psychiatric Comorbidities in Parkinson’s Disease Seen through the Prism of Genomics and Epigenetics. Psychiatr Danub 33(suppl 4):456–462
Citraro R, Leo A, de Caro C, Nesci V, Gallo Cantafio ME, Amodio N et al (2020) Effects of histone deacetylase inhibitors on the development of epilepsy and psychiatric comorbidity in WAG/Rij rats. Mol Neurobiol 57(1):408–421
Cain SM, Snutch TP (2011) Voltage-gated calcium channels and disease. BioFactors 37(3):197–205
Pietrobon D (2013) Calcium channels and migraine. Biochimica et Biophysica Acta (BBA)-Biomembranes 1828(7):1655–65
Anne D, Elisabeth TL, Marie-Germaine B (2002) The genetics of Migraine. Lancet Neurol 1:285–293
Turnpenny PD, Ellard S (2012) Alagille syndrome: pathogenesis, diagnosis and management. Eur J Hum Genet 20(3):251–257
Rubino E, Fenoglio P, Gallone S, Govone F, Vacca A, de Martino P et al (2013) Genetic variants in the NOTCH4 gene influence the clinical features of migraine. J Headache Pain 14(1):1–5
Vriens J, Appendino G, Nilius B (2009) Pharmacology of vanilloid transient receptor potential cation channels. Mol Pharmacol 75(6):1262–1279
Julius D (2013) TRP channels and pain. Annu Rev Cell Dev Biol 29:355–384
Jardín I, López JJ, Diez R, Sánchez-Collado J, Cantonero C, Albarrán L et al (2017) TRPs in pain sensation. Front Physiol 8:392
Russell FA, King R, Smillie SJ, Kodji X, Brain SD (2014) Calcitonin gene-related peptide: physiology and pathophysiology. Physiol Rev 94(4):1099–1142
Goadsby PJ, Edvinsson L (1993) The trigeminovascular system and migraine: studies characterizing cerebrovascular and neuropeptide changes seen in humans and cats. Ann Neurol 33(1):48–56
Tso AR, Goadsby PJ (2017) Anti-CGRP monoclonal antibodies: the next era of migraine prevention? Curr Treat Options Neurol 19(8):1–11
Edvinsson L, Haanes KA, Warfvinge K, Krause DN (2018) CGRP as the target of new migraine therapies—successful translation from bench to clinic. Nat Rev Neurol 14(6):338–350
Charles A, Pozo-Rosich P (2019) Targeting calcitonin gene-related peptide: a new era in migraine therapy. Lancet 394(10210):1765–1774
Benemei S, Dussor G (2019) TRP channels and migraine: recent developments and new therapeutic opportunities. Pharmaceuticals 12(2):54
Torres-Ferrús M, Gallardo VJ, Alpuente A, Caronna E, Gine-Cipres E, Pozo-Rosich P (2021) The impact of anti-CGRP monoclonal antibodies in resistant migraine patients: a real-world evidence observational study. J Neurol 268(10):3789–98
Alpuente A, Gallardo VJ, Asskour L, Caronna E, Torres-Ferrus M, Pozo-Rosich P (2022) Salivary CGRP can monitor the different migraine phases: CGRP (in) dependent attacks. Cephalalgia 42(3):186–96
Alpuente A, Gallardo VJ, Asskour L, Caronna E, Torres-Ferrus M, Pozo-Rosich P (2022) Salivary CGRP and Erenumab Treatment Response: Towards Precision Medicine in Migraine. Ann Neurol 92(5):846–59
Yang D gang, Gao Y yao, Yin Z qun, Wang X rui, Meng X she, Zou T feng, et al (2022) Roxadustat alleviates nitroglycerin-induced migraine in mice by regulating HIF-1α/NF-κB/inflammation pathway. Acta Pharmacol Sin (2):308–20
Mirouliaei M, Fallah R, Bashardoost N, Partovee M, Ordooei M (2012) Efficacy of levothyroxine in migraine headaches in children with subclinical hypothyroidism. Iran J Child Neurol 6(4):23
Paolucci M, Altamura C, Vernieri F (2021) The role of endothelial dysfunction in the pathophysiology and cerebrovascular effects of migraine: a narrative review. J Clin Neurol 17(2):164
Peres MFP, del Rio MS, Seabra ML, Tufik S, Abucham J, Cipolla-Neto J et al (2001) Hypothalamic involvement in chronic migraine. J Neurol Neurosurg Psychiatry. 71(6):747–51
Borsook D, Erpelding N, Lebel A, Linnman C, Veggeberg R, Grant PE et al (2014) Sex and the migraine brain. Neurobiol Dis 68:200–214
MacGregor EA, Frith A, Ellis J, Aspinall L, Hackshaw A (2006) Incidence of migraine relative to menstrual cycle phases of rising and falling estrogen. Neurology 67(12):2154–2158
Phillips WJ, Ostrovsky O, Galli RL, Dickey S (2006) Relief of acute migraine headache with intravenous oxytocin: report of two cases. J Pain Palliat Care Pharmacother 20(3):25–28
Tzabazis A, Kori S, Mechanic J, Miller J, Pascual C, Manering N et al (2017) Oxytocin and migraine headache. Headache 57:64–75
Chandramohan R, Wu PY, Phan JH, Wang MD (2013) Benchmarking RNA-Seq quantification tools. Annu Int Conf IEEE Eng Med Biol Soc 2013:647–50. https://doi.org/10.1109/EMBC.2013.6609583
Hayer KE, Pizarro A, Lahens NF, Hogenesch JB, Grant GR (2015) Benchmark analysis of algorithms for determining and quantifying full-length mRNA splice forms from RNA-seq data. Bioinformatics 31(24):3938–3945
Robinson MD, Oshlack A (2010) A scaling normalization method for differential expression analysis of RNA-seq data. Genome Biol 11(3):1–9
Rao MS, van Vleet TR, Ciurlionis R, Buck WR, Mittelstadt SW, Blomme EAG et al (2019) Comparison of RNA-Seq and microarray gene expression platforms for the toxicogenomic evaluation of liver from short-term rat toxicity studies. Front Genet 9:636
Kemper RHA, Meijler WJ, Korf J, ter Horst GJ (2001) Migraine and function of the immune system: a meta-analysis of clinical literature published between 1966 and 1999. Cephalalgia 21(5):549–557
Forcelini CM, Dantas DCM, Luz C, Santin R, Stein AT, Barros HMT et al (2011) Analysis of leukocytes in medication-overuse headache, chronic migraine, and episodic migraine. Headache 51(8):1228–1238
Grazzi L, Corsini E, Ciusani E, Usai S, Vasco C, Bussone G (2014) Evaluation of immune parameters in chronic migraine with medication overuse. Neurol Sci 35(SUPPL. 1):171–173
Kok MGM, de Ronde MWJ, Moerland PD, Ruijter JM, Creemers EE, Pinto-Sietsma SJ (2018) Small sample sizes in high-throughput miRNA screens: a common pitfall for the identification of miRNA biomarkers. Biomol Detect Quantif 15:1–5
This study was solely funded by the Migraine Research Foundation.
Ethics approval and consent to participate
The study was approved by our Institutional Review Board (PR(IR)224/2016) and all participants consented to participated.
VJ.G., JB.G.G. and A.L. have nothing to declare. M.T.F. has received honoraria as a consultant or speaker for Allergan-Abbvie, Almirall, Chiesi, Eli Lilly, Novartis and Teva. A.A. has received honoraria from Allergan-Abbvie, Novartis, Chiesi. E.C. reports no diclosures. P.P.R. has received honoraria as a consultant and speaker for: AbbVie, Amgen, Biohaven, Chiesi, Eli Lilly, Lundbeck, Medscape, Novartis, Pfizer and Teva. Her research group has received research grants from AbbVie, Novartis and Teva; and, has received funding for clinical trials from Alder, AbbVie, Amgen, Biohaven, Electrocore, Eli Lilly, Lundbeck, Novartis and Teva. She is the Honorary Secretary of the International Headache Society. She is in the editorial board of Revista de Neurologia, associate editor for Cephalalgia, Headache, Neurologia, Frontiers of Neurology. She is a member of the Clinical Trials Guidelines Committee and Scientific Committee of the International Headache Society. She has edited the Guidelines for the Diagnosis and Treatment of Headache of the Spanish Neurological Society. She is the founder of www.midolordecabeza.org. Pozo-Rosich P does not own stocks from any pharmaceutical company.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Gallardo, V.J., Gómez-Galván, J.B., Asskour, L. et al. A study of differential microRNA expression profile in migraine: the microMIG exploratory study. J Headache Pain 24, 11 (2023). https://doi.org/10.1186/s10194-023-01542-z
- Epigenetic mechanisms
- Personalized medicine