Longitudinal changes in functional connectivity and pain-induced brain activations in patients with migraine: a functional MRI study pre- and post- treatment with Erenumab
The Journal of Headache and Pain volume 23, Article number: 159 (2022)
Migraine involves central and peripheral nervous system mechanisms. Erenumab, an anti-calcitonin gene-related peptide (CGRP) receptor monoclonal antibody with little central nervous system penetrance, is effective for migraine prevention. The objective of this study was to determine if response to erenumab is associated with alterations in brain functional connectivity and pain-induced brain activations.
Adults with 6–25 migraine days per month during a 4-week headache diary run-in phase underwent pre-treatment brain functional MRI (fMRI) that included resting-state functional connectivity and BOLD measurements in response to moderately painful heat stimulation to the forearm. This was followed by two treatments with 140 mg erenumab, at baseline and 4 weeks later. Post-treatment fMRI was performed 2 weeks and 8 weeks following the first erenumab treatment. A longitudinal Sandwich estimator analysis was used to identify pre- to post-treatment changes in resting-state functional connectivity and brain activations in response to thermal pain. fMRI findings were compared between erenumab treatment-responders vs. erenumab non-responders.
Pre- and post-treatment longitudinal imaging data were available from 32 participants. Average age was 40.3 (+/− 13) years and 29 were female. Pre-treatment average migraine day frequency was 13.8 (+/− 4.7) / 28 days and average headache day frequency was 15.8 (+/− 4.4) / 28 days. Eighteen of 32 (56%) were erenumab responders. Compared to erenumab non-responders, erenumab responders had post-treatment differences in 1) network functional connectivity amongst pain-processing regions, including higher global efficiency, clustering coefficient, node degree, regional efficiency, and modularity, 2) region-to-region functional connectivity between several regions including temporal pole, supramarginal gyrus, and hypothalamus, and 3) pain-induced activations in the middle cingulate, posterior cingulate, and periaqueductal gray matter.
Reductions in migraine day frequency accompanying erenumab treatment are associated with changes in resting state functional connectivity and central processing of extracranial painful stimuli that differ from erenumab non-responders.
Prior publications have demonstrated abnormal brain structure, functional connectivity, and stimulus-induced brain responses amongst those with migraine [1,2,3,4,5,6]. Other studies have demonstrated cycling activity and functional connectivity of brain regions such as the hypothalamus, pons, and trigeminal nucleus, that correlate with the occurrence of migraine attacks [3, 4, 7]. Only a few studies have investigated changes in resting-state functional connectivity or pain-induced brain activations associated with response to migraine preventive treatment [8,9,10,11,12]. These publications demonstrate that migraine preventive treatments, even those that are unlikely to directly access the brain, may be associated with changes in brain functional connectivity and pain-induced activations.
Erenumab (erenumab-aooe in the U.S.), a monoclonal antibody (mAb) that targets the calcitonin gene-related peptide (CGRP) receptor, is effective for the prevention of episodic and chronic migraine [13,14,15]. The anti-CGRP pathway monoclonal antibodies are very large molecules that have little penetrance through the blood-brain barrier and therefore are not thought to exert their effects directly in the central nervous system. None-the-less, they are effective migraine preventive medications, possibly through their action at locations outside of the blood-brain barrier such as the trigeminal ganglia, trigeminal nerves, or dura mater [16,17,18]. A previously published functional magnetic resonance imaging (fMRI) study demonstrated that erenumab treatment is associated with changes in central processing of trigeminal pain and brain resting state functional connectivity at 2 weeks after starting treatment . Another study investigated the impact of treatment with galcanezumab, another anti-CGRP pathway mAb, on brain activation in response to trigeminal nociception measured 2–3 weeks after starting galcanezumab . The objective of our study was to determine if participants with migraine who respond to erenumab, an anti-CGRP pathway mAb targeting the receptor, have post-treatment changes in functional connectivity and pain-induced brain activations that differ from erenumab non-responders. To address this objective, we interrogated 1) resting-state functional connectivity amongst brain regions involved with pain processing, including global and regional graph theory network measures and region-to-region connectivity, and 2) brain responses to extracranial painful thermal stimuli. Furthermore, our study investigated the timing of these fMRI changes, with measurements available from 2 weeks and 8 weeks following the first erenumab treatment.
Informed consent and study registration
This study was approved by the Mayo Clinic Institutional Review Board. Each research participant reviewed and signed a consent form after participating in the informed consent process. The study was registered at clinicaltrials.gov (NCT03773562).
Adults between the ages of 18–65 years who had episodic migraine or chronic migraine and self-reported a history of 6–25 migraine days per month during the prior 3 months were eligible for participation. Migraine diagnoses were made according to the International Classification of Headache Disorders 3rd edition criteria . Individuals were excluded if they were older than 50 years of age at migraine onset, had migraine onset within the prior 12 months, had cluster headache or hemiplegic migraine, continuous headache (i.e. no headache free periods during the 1 month prior to screening), used opioids or butalbital on 6 or more days per month, had no therapeutic response to adequate trials of migraine preventive medications from four or more different medication classes, were currently taking more than two migraine preventive medications, had received botulinum toxin treatment within the prior 4 months or nerve blocks used for treatment of headache within 2 months, had a history of myocardial infarction, stroke, transient ischemic attack, unstable angina, coronary artery bypass surgery, or other revascularization procedures within the prior 12 months, had contraindications to MRI, were pregnant or lactating, were not willing to use a reliable form of contraception, or had received a CGRP pathway mAb within the prior 4 months. After completing the four-week run-in diary phase, those who had between 6 and 25 migraine days and were at least 80% compliant with providing headache diary data were eligible to continue in the study.
For those participants who qualified for continued study participation after the run-in phase, there were a total of six research visits during a 16-week period. (Table 1) Questionnaires were completed at each visit and a headache diary was maintained for the entire 16 weeks. Brain MRIs and quantitative sensory testing (QST) were performed three times. Each participant received two treatments with 140 mg of erenumab. Adverse events (AEs) were collected through 12 weeks following the first erenumab treatment.
Questionnaires and Structured Interviews: Structured interviews and questionnaires were used to collect data on participant demographics, migraine history and characteristics, medications, medical history, migraine-related disability (Migraine Disability Assessment (MIDAS)), depression (Beck Depression Inventory (BDI)), and cutaneous allodynia (Allodynia Symptom Checklist – 12 (ASC-12)) [20,21,22]. Adverse events were collected at each study visit.
Headache Diary: The headache diary was an e-diary that was developed within REDCap . Each day, participants were prompted to provide information about the presence of headache and associated symptoms. The headache diary was used to determine the number of headache days and migraine days experienced by each participant. A migraine day was defined as any calendar day during which a person experienced a qualified migraine headache (onset, continuation, or recurrence of the migraine headache). A qualified migraine headache was defined as a migraine with or without aura, lasting for at least 30 minutes, and meeting at least one of the following criteria (a and/or b): a) at least 2 of the following pain features: unilateral, throbbing, moderate to severe, exacerbated with exercise/physical activity; b) at least 1 of the following associated symptoms: nausea and/or vomiting, photophobia and phonophobia. If the participant took a migraine-specific as-needed medication (e.g., triptan or ergotamine) on a calendar day, then it was counted as a migraine day regardless of the duration and pain features/associated symptoms.
Quantitative Sensory Testing (QST): A Medoc Pathway platform with a 30 mm × 30 mm thermode was used for QST to determine the temperature required to cause moderately intense heat pain. Moderately intense heat pain was defined by a participant rating their pain intensity between 4 and 7, ratings consistent with moderate intensity pain, using an 11-point scale in which zero indicates no pain and 10 indicates the most severe pain. Initially, heat pain threshold was determined for each participant using the method of limits. The thermode was placed on the left forearm and fastened with a Velcro strap. The thermode increased in temperature from 32 degrees Celsius at a rate of one degree Celsius per second until the participant pushed a button indicating that the stimulus changed from a heat sensation to a painful sensation. The average of three trials was considered the heat pain threshold. Participants were then stimulated for 7.5 seconds with a temperature equal to their heat pain threshold plus 1 degree Celsius. Individuals then rated the pain intensity. If pain intensity was rated below 4, the temperature was increased in 0.5 degree Celsius increments until rated between 4 and 7. Conversely, if the pain intensity was rated greater than 7, the stimulation temperature was decreased in 0.5 degree Celsius increments until the participant rated the pain between 4 and 7/10 in intensity. QST methods were consistent with those previously reported .
Brain MRI: Brain MRIs were completed on a Siemens Skyra 3 T scanner. Sequences included a high-resolution 3D T1-weighted MPRAGE sequence (TE 3.03 ms, TR 2400 ms, voxel size 1x1x1.25 mm, field of view 256x256x160 mm, flip angle 8 degrees), two five-minute runs of resting-state blood-oxygen-level-dependent (BOLD) collection with the participant relaxing with eyes closed (TE 27 ms, TR 2500 ms, voxel size 4x4x4 mm, field of view 256x256x152 mm, flip angle 90 degrees), and three runs of event-related BOLD collection (TE 27 ms, TR 2500 ms, voxel size 4x4x4 mm, field of view 256x256x152 mm, flip angle 90 degrees) during which participants received painful thermal stimuli. During the event-related paradigm, the thermode was attached to the left forearm with a Velcro strap. An auditory cue was followed by moderately intense heat pain, using the temperature determined to cause moderately intense pain via testing of the participant just before the MRI. Heat pain stimuli were delivered nine times for each patient, with each stimulus lasting for 7.5 seconds and the interstimulus interval varying from 44 to 46 seconds. Auditory cue followed by no thermal stimulus (i.e., no change in the thermode temperature) was randomly inserted into the event-related paradigm to account for brain activations associated with pain anticipation or alerting responses. The event-related paradigm methods were consistent with those previously reported .
Erenumab Treatment: Participants received treatment with 140 mg of erenumab via subcutaneous injection into the upper arm. Injections were completed in the office during two research visits separated by 4 weeks.
Pain-induced activation MRI data processing
All preprocessing and General Linear Model (GLM) estimation of whole-brain activation patterns were performed using SPM12 software (Wellcome Department of Cognitive Neurology, Institute of Neurology, London, UK) interfaced with MATLAB version 11.0 (Mathworks, Natick, MA, USA). Functional images were realigned to the mean volume in the series, motion-corrected, realigned to each individual’s structural images, and smoothed using an 8 mm Full Width Half Maximum (FWHM) Gaussian kernel. Data were checked for excessive motion and all included scans showed < 3 mm movement in any direction. First-level single subject analysis was performed within SPM12 using a GLM approach with regression of 6 motion parameters. The immediate BOLD response to pain stimuli used an event-related design with the measured duration of hemodynamic response equal to 12.5 s (i.e. five TRs), starting with the MRI frame during which the painful heat stimulus began and ending 12.5 s later. The delayed response to pain was also investigated, with measurements starting 12.5 s after onset of painful heat stimulus and ending 12.5 s later . Brain regions activated in response to painful stimuli were identified by generating contrast maps representing activations associated with auditory cue followed by painful stimulation vs. activations associated with auditory cue followed by no stimulation. Each individual’s contrast map was normalized to standard stereotaxic space using the Montreal Neurological Institute template (MNI 152).
Statistical analysis was performed using the Sandwich estimator (SWE) for neuroimaging longitudinal data toolbox version 1.2.8 SWE (SWE, Guillaume 2014 & 2015) interfaced with MATLAB . All contrast maps were added to the SWE longitudinal model with default settings for small sample adjustments and estimated degrees of freedom. Responder vs non-responder group-visit effects were examined using F-statistic to include positive and negative relationships. For the pain vs no stimulation contrast, a false discovery rate with p < 0.05 with a cluster forming threshold of 10 voxels was set for the main effect and uncorrected p < 0.005 with a cluster forming threshold of 50 voxels for the group visit effect . The delayed response used an uncorrected p < 0.001 with a cluster forming threshold of 10 voxels for the main effect and uncorrected p < 0.005 with a cluster forming threshold of 25 voxels for the group visit effect.
Resting functional connectivity MRI data processing
Resting-state functional MRI data were pre-processed using standard procedures in SPM12 including the following steps: slice-time correction, motion correction, re-alignment to first image of the set, re-alignment to an average MNI template, and smoothing with a 6 mm FWHM Gaussian kernel. Further post-processing in AFNI 3dTproject included band-pass filtering (0.01–0.1 HZ) after removal of nuisance signals from framewise displacement, cerebrospinal fluid signal, and linear drift .
A region of interest approach (ROI) was used to interrogate functional connectivity. Thirty-one bilateral ROIs and one midline ROI were selected based on previous findings demonstrating that these regions participate in pain processing, migraine, and/or multisensory integration [3, 7, 10, 29,30,31,32,33,34,35,36]. [Table 2] Each ROI was an 8 mm diameter sphere created using MarsBaR toolbox in SPM12. The Pearson correlation coefficient was calculated between the signals of the ROI pairs where motion did not exceed 2 mm. The first five frames were excluded to allow the signal to reach steady state. Following the correlation analysis, Fisher r-z transforms were calculated for each ROI-to-ROI pair and third visit mean functional connectivity was compared between erenumab-responders and non-responders using two tailed t-tests (uncorrected p-value less than 0.005 considered significant). An undirected weighted adjacency matrix was then created for each participant with ROI-ROI z-scores as the edge weights. The diagonal was set to zero to exclude self-connections. Nodes with z-scores less than 0.1 were set to zero to construct undirectional weighted functional connectomes. Global graph theory measurements were performed for global efficiency . Local measurements were performed for efficiency, clustering coefficient, degrees centrality (betweenness and degree), and modularity [38,39,40]. Two sample t-tests (2 tailed) were performed comparing erenumab-responders to non-responders at the 3rd visit with significance defined as an uncorrected p-value less than 0.05.
Erenumab responder definition
Treatment response was defined as at least a 50% reduction in the frequency of migraine days during weeks 5–8 compared to the 4-week pre-treatment run-in phase as recorded in the headache diary. For months with eight or fewer days of missing diary data, data were imputed via simple proportional imputation. Larger amounts of missing diary data were imputed using diary data from the current 4-week period and the previous 4-week period combined.
Comparing Erenumab responders to Erenumab non-responders
Participant demographics, headache characteristics, ASC-12 scores, BDI scores, MIDAS scores, headache and migraine presence at the time of research testing, and pain stimulus temperatures were compared using two-sided t-tests or Fisher exact test, as appropriate, with p < 0.05 being considered significant. Changes in headache frequency, migraine frequency, and ASC-12 scores were compared between erenumab responders and non-responders using two-sided t-tests with p < 0.05 being considered significant.
The flow of participants through the study is illustrated in Fig. 1. Amongst the 50 participants who entered the run-in phase, 10 were excluded because they did not experience 6–25 migraine days during the run-in phase and 4 were lost to follow-up prior to having a baseline MRI and first erenumab treatment. Amongst the 36 remaining participants, one was withdrawn due to abnormal brain MRI findings and three were lost to follow-up. Thirty-two participants had pre- and post-erenumab fMRI and were thus included in this analysis, 29 of whom at least completed the pre-treatment and 8-week post-treatment MRI. MRI data from all three timepoints (pre-treatment, 2-week post-treatment, and 8-week post-treatment) were available from 25 participants. On average, data from 27.5 diary days was provided during the four-week run-in period and 24.6 days during the 5–8 week period.
The participants (n = 32) average age was 40.3 (+/− 13) and 29 were female. 31 were white, 1 was Asian, and 4 were Hispanic. Pre-treatment, they had migraine for an average of 21.9 (+/− 12.1) years, had an average of 13.8 (+/− 4.7) migraine days per month and 15.8 (+/− 4.4) headache days per month (which includes the migraine days) (Table 3). Pre-treatment, 11 had episodic migraine and 21 had chronic migraine. All participants were using migraine abortive medications. Twenty-four were using prescription abortive medications and eight were using over-the-counter abortive medications only. Nineteen of thirty-two participants had medication overuse at baseline. Migraine prescription preventive medications, in addition to erenumab, were being used by ten participants, and included propranolol (n = 4 participants), topiramate (n = 2), gabapentin (n = 2), venlafaxine (n = 2), zonisamide (n = 1), and amitriptyline (n = 1) (two patients were using two migraine preventive medications in addition to erenumab).
Amongst all pre-treatment research visits for testing, headache was present during 15/32 (46.9%) visits and migraine was present on 5/32 (15.6%) visits. Amongst participants who went on to be erenumab-responders, headache was present during 9/18 (50%) pre-treatment visits compared to those who did not go on to be erenumab responders in whom headache was present during 6/14 (42.9%) pre-treatment visits (p = 0.73). Amongst participants who went on to be erenumab-responders, migraine was present during 4/18 (22.2%) pre-treatment visits and during 1/14 (7.1%) pre-treatment visits for those who went on to be erenumab non-responders (p = 0.35) .
Amongst all post-treatment research visits for testing, headache was present during 12/59 visits and migraine was present during none of the visits. Amongst erenumab responders, headache was present during 5/59 (8.5%) of post-treatment visits compared to 7/59 (11.9%) post-treatment visits among erenumab non-responders (p = 0.76).
Eighteen (56%) participants were erenumab responders and 14 (44%) were non-responders. Changes in headache day frequency, migraine day frequency, and ASC-12 scores are shown in Table 4. There were no differences in post-treatment thermode temperatures required to cause moderate intensity pain. Post-treatment amongst all participants the average temperature was 45.3 °C (+/− 2.2 °C); 45.6 °C (+/− 2.3 °C) in erenumab responders vs. 45.0 °C (+/− 2.2 °C) in erenumab non-responders (p = 0.32).
AEs were reported by 25 participants who reported a total of 36 AEs. There were no serious AEs. Eighteen were considered mild, 13 moderate, and 5 severe. Fourteen AEs were considered unrelated to study medication (dizziness × 2, and one each of flank pain, loss of balance, weight gain, twitching of arms, ear infection, urinary tract infection, agitation, communication issues, insomnia, COVID-19 infection, allergic reaction, West Nile Virus infection), 5 were unlikely to be related (itching × 2, and one each of dizziness, skin sensitivity, tinnitus), 1 was possibly related (worsening headache), 15 were probably related (constipation × 14, abdominal pain × 1), and 1 was definitely related (injection-site erythema). Two participants withdrew from the study due to adverse events, including one patient who had constipation and abdominal pain and one patient who reported worsening headaches.
Pain-induced brain activations
The main effect of painful heat stimulation to the left forearm at the baseline visit is demonstrated in Fig. 2. This main effect demonstrates activation of ‘pain matrix’ regions, i.e., those brain regions that have consistently shown to be activated by painful stimulation and provides evidence that the pain stimuli resulted in the expected brain activations.
When comparing erenumab responders to non-responders, at post-treatment visits there were significant differences in brain activations occurring immediately after onset of painful stimuli in the left middle cingulate, left posterior cingulate, and right putamen (Fig. 3).
The left middle cingulate and left posterior cingulate showed significant differences in pain-induced activations between erenumab-responders and non-responders at visit 3 (p < 0.001) whereas the putamen showed a transient difference at visit 2 (p = 0.002). In all three regions the baseline activation is similar between groups.
Analyses of delayed BOLD activations in response to painful stimulation demonstrated differences between erenumab-responders vs. non-responders for the periaqueductal gray and left frontal supplemental area (Fig. 4). The periaqueductal gray matter response to pain was greater in erenumab-responders vs. non-responders at the 3rd visit (p = 0.003). The frontal supplemental region had a greater response in the 1st visit in non-responders compared to responders (p = 0.005).
Resting state functional connectivity
Pre-treatment, the global efficiency amongst all a priori selected ROIs did not differ between erenumab responders (24.6% ± 1%) vs. non-responders (24.3% ± 0.7%, p = 0.81). However, there was a difference in global efficiency at 8 weeks after the first erenumab treatment in responders (26.2% ± 1.1%) vs. non-responders (23.3% ± 0.7%, p = 0.04). The longitudinal measures of global efficiency are shown in Fig. 5.
The local network measures for each ROI that differed between erenumab-responders vs. non-responders are shown in Table 5. Differences for several network measures were identified, including regional efficiency, betweenness, clustering coefficient, node degree, and modularity.
The pairwise differences in ROI-to-ROI correlations between erenumab-responders vs. non-responders at 8 weeks following the first erenumab treatment are shown in Table 6. These functional connections included regions such as the supramarginal gyrus, inferior lateral parietal, hypothalamus, temporal pole, middle temporal, middle occipital, and middle frontal lobes, dorsolateral and ventromedial prefrontal cortices, and the pulvinar.
The main finding of this study is that effective migraine treatment with erenumab is associated with changes in brain activations in response to extracranial thermal pain and with changes in resting-state functional connectivity. fMRI differences between erenumab-responders vs non-responders were seen within 8 weeks of initiating erenumab treatment, with some differences starting to emerge at 2 weeks. Brain regions impacted across the different analyses in this study included the hypothalamus, inferior lateral parietal, temporal pole, supramarginal gyrus, amygdala, and periaqueductal gray, amongst others. As discussed further below, these regions are known to have important roles in migraine pain processing including their likely roles in migraine attack generation, pain perception, pain modulation, and multisensory integration.
Our pain-induced activation analyses identified differences between erenumab responders vs. non-responders in activations of regions in the cingulate cortex, putamen, frontal lobe, and periaqueductal gray matter. Brain activations occurring immediately (i.e., within 12.5 seconds) after painful stimuli and delayed activations occurring between 12.5 s to 25 seconds after painful stimuli were investigated. Prior studies have demonstrated a biphasic hemodynamic brain response to painful heat stimuli, consisting of a stimulus-locked response and a second peak delayed by about 12.5 seconds . It has been suggested that the initial BOLD response might be related to pain processing via myelinated A-delta fibers, while the delayed BOLD response could be associated with processing via the slower C fibers [25, 41]. It has also been theorized that the initial BOLD response might represent a fast, non-conscious processing of pain that helps to quickly determine the threat level, while the delayed response might represent a more conscious processing of the painful stimulus [25, 41, 42]. Our study identified differences in immediate and delayed BOLD responses to pain when comparing erenumab-responders vs. non-responders.
For immediate BOLD responses, there were differences in the middle and posterior cingulate at 8 weeks following the first erenumab treatment and in the putamen only at 2 weeks following the first erenumab treatment. Changes in BOLD responses in the middle cingulate and posterior cingulate were already seen at 2 weeks after the first erenumab treatment, but the differences between erenumab-responders and non-responders were not significant until the eight-week MRI. The middle cingulate participates in rapid behavioral adaptive responses to the threat associated with pain, as well as cognitive and affective components of pain processing [43, 44]. A prior fMRI study of migraine demonstrated a strong positive correlation between pain-induced activation in the middle cingulate and headache frequency (r = 0.627, p = .001), providing relatively strong evidence for its role in migraine . In our study, it is possible that the increase in pain-induced activation of the middle cingulate amongst erenumab-responders is due to the pain stimulus as being perceived as more novel and thus as a greater threat compared to the erenumab non-responders who are experiencing more frequent pain, with the painful stimulus thus losing its novelty. The posterior cingulate, a key region of the default mode network, participates in self-referential processing of stimuli, including externally generated pain . Numerous pain and migraine studies have demonstrated atypical functional connectivity of default mode network regions [46,47,48,49]. A prior migraine study demonstrated that migraine improvement (i.e. reduced time with headache) was associated with cortical thickness changes in the left posterior cingulate . The putamen participates in sensory-discriminative aspects of pain processing including the determination of pain sensitivity . Prior migraine studies have demonstrated atypical functional connectivity and structure of the putamen [52, 53]. In our study, differences in putamen activation were detected at two-weeks following the first erenumab treatment, but not at 8 weeks, a pattern that is difficult to interpret.
Delayed BOLD response differences between erenumab-responders and non-responders were identified in the periaqueductal gray at 8 weeks following initiation of erenumab treatment and in the frontal supplementary area prior to starting erenumab (and thus not related to treatment). The periaqueductal gray is a key region of the pain modulating pathway, predominantly involved in pain inhibition. Numerous studies have identified atypical periaqueductal gray structure and function associated with migraine, the extent to which might correlate with migraine disease severity [54,55,56,57,58]. Less effective pain modulation by the periaqueductal gray could be a mechanism by which those with migraine experience more severe pain, allodynia, and increased frequency of migraine attacks [54,55,56,57]. Increased pain-induced activation of the periaqueductal gray amongst erenumab-responders might reflect a stronger pain-inhibitory response.
Overall, our pain activation studies suggest that effective erenumab treatment is associated with changes in cognitive, affective, sensory-discriminative, and modulating aspects of pain processing.
Resting-state functional connectivity
Our resting-state functional connectivity analyses included investigation of graph theory network parameter differences and ROI-to-ROI static functional connectivity differences between erenumab-responders vs. non-responders. Graph theory provides a method for quantitatively describing the topological organization of brain networks . Graph theory measures in our analyses included global and regional efficiency, betweenness, clustering coefficient, node degree, and modularity [39, 40]. Efficiency is the inverse of the minimum path length between an ROI and all other ROIs in a network . Global and local efficiency reflect a network’s ability to transmit information at the global and local levels . Betweenness of a ROI is the number of shortest paths between any two ROIs that run through the ROI . Clustering coefficient describes the local connectedness of a network and is calculated by determining the number of connections that a ROI has with its immediate neighbors divided by all of its possible connections . Node degree is a measure of the number of connections to an ROI . Modularity measures how connected ROIs are to members of their own group, identifying subnetworks within a larger network.
Erenumab responders had an increase in global efficiency from pre-treatment to eight-weeks post-treatment and greater global efficiency at eight-weeks compared to erenumab non-responders. An increase in global efficiency suggests that erenumab-response was associated with an improvement in the ability of the studied regions, those that participate in various aspects of pain processing and modulation, to functionally communicate on a global scale. At eight-weeks post-treatment, compared to erenumab non-responders, responders had ROIs with higher efficiency, cluster coefficient, node degree, and modularity, findings suggesting a greater ‘small world’ quality of the network. Regions most strongly highlighted by these differences included those in the hypothalamus, amygdala, and inferior parietal lobe. The hypothalamus likely plays an important role in the generation of migraine attacks [59, 60]. The amygdala contributes to affective and attentional responses to pain and pain modulation . The inferior parietal region, like the posterior cingulate discussed above, is a key region of the default mode network . Other regions included the temporal pole, supramarginal gyrus, caudate, cuneus, trigeminal nucleus, lingual gyrus, and middle cingulate.
ROI-to-ROI differences in functional connectivity at 8 weeks between erenumab-responders vs. non-responders included regions that were also identified as having differences in graph theory network measures, including the supramarginal gyrus, inferior lateral parietal, hypothalamus, and temporal pole. Additionally, several other regions were involved in these functional connections including regions in the middle temporal, middle occipital, and middle frontal lobes, dorsolateral and ventromedial prefrontal cortices, and the pulvinar.
Study results in context of prior studies
Prior studies have interrogated the impact of migraine treatment on brain functional connectivity and brain activations. Krebs and colleagues demonstrated that treatment with sphenopalatine ganglion blocks was associated with changes in functional connectivity amongst salience network and executive network regions . Russo and colleagues investigated pain-induced brain activations before and after 60 days of treatment with external trigeminal neurostimulation . Neurostimulation treatment was associated with a reduced BOLD response in the anterior cingulate cortex. Acupuncture treatment has been associated with an increase in periaqueductal gray functional connectivity with anterior cingulate cortex among individuals who have migraine without aura . A pre-treatment and 2–3 week post-treatment study of galcanezumab, a CGRP ligand mAB, for migraine demonstrated that galcanezumab decreased hypothalamic activation in response to nociceptive trigeminal stimulation to a greater extent in galcanezumab-responders vs. non-responders . There were also responder-specific decreases in BOLD activation in the inferior parietal lobule, insula and parahippocampal gyrus. Spinal trigeminal nucleus functional connectivity changes from the pre-treatment to post-treatment scans were interrogated for all treated patients and demonstrated weakened connectivity with hypothalamus and superior temporal gyrus and stronger connectivity with the cerebellum, middle temporal gyrus, and insula at the post-treatment timepoint. The study most closely related to ours investigated the impact of erenumab on brain activations in response to nociceptive trigeminal stimulation . In that study, 27 individuals with migraine underwent fMRI prior to and 2 weeks after treatment with 70 mg of erenumab. During the fMRI paradigm, intranasal ammonia was used as a painful stimulus. 63% of participants were considered erenumab treatment responders, which was defined as at least a 30% reduction in headache days during the first month following treatment. Amongst all patients there were post-treatment decreases in pain-induced activations in the thalamus, lingual gyrus, middle temporal gyrus, operculum, and cerebellum. Compared to non-responders, erenumab-responders had a significant reduction of activation in the hypothalamus, insula, and cerebellum. Analysis of hypothalamic functional connectivity amongst all treated patients demonstrated a reduction in connectivity strength with the temporal lobe, hippocampus, parahippocampus, fusiform gyrus, cerebellum, red nucleus, and spinal trigeminal nucleus, and an increase in connectivity strength with the anterior insula. This study and ours complement one another, both demonstrating changes in pain-induced activation and resting functional connectivity associated with erenumab treatment and response. Our study adds to the literature since it determined treatment response during weeks 5–8 after starting treatment, studied 140 mg of erenumab rather than 70 mg, investigated early (i.e. 2 weeks) and later (i.e. 8 weeks) fMRI changes after initiating treatment, utilized two treatments with erenumab rather than one, used an extra-trigeminal painful heat stimulus for the event-related paradigm, and interrogated graph theory network measures of functional connectivity.
Relationship between fMRI changes and Erenumab treatment
How treatment with erenumab is associated with changes in pain-induced brain activations and functional connectivity is a matter of debate. It is perhaps unlikely that the small amount of erenumab that might cross the blood-brain barrier could exert a meaningful central effect and have a direct impact on brain processing of painful stimuli. Alternatively, erenumab might alter brain pain processing indirectly, via its impact on peripheral structures such as the trigeminal ganglia, trigeminal nerves, or dura mater [16,17,18]. Finally, it is possible that the changes in pain processing demonstrated in this study are attributable to the reduction in migraine days associated with erenumab response, but not specifically attributable to the mechanisms by which erenumab exerts therapeutic effects.
Study considerations and limitations
Considerations and limitations of our study include: 1) Although the number of MRIs completed and included in this study is relatively large (n = 86), the number of erenumab-responders (n = 18) and non-responders (n = 14) is relatively small. Larger sample sizes might allow for more stringent statistical corrections for multiplicity. 2) Our study is not able to determine if the pre- to post-erenumab changes in pain-induced brain activations and functional connectivity are directly attributable to erenumab or if there would be similar findings associated with reductions in migraine frequency regardless of the specific reason for such a reduction. Optimally, future studies would include migraine treatments that work via different mechanisms and individuals who have longitudinal reductions in migraine frequency in the absence of treatment. 3) Inclusion of a healthy control group would allow for better interpretation of the pre-to-post treatment changes, whether the changes are consistent with a “normalization” of brain function or adaptive changes, for example. 4) Neuroimaging research studies use different statistical corrections and cluster forming thresholds for determining the significance of results. More stringent methods increase the likelihood for type II error, while less stringent methods increase the likelihood for type I error. The approaches taken in this study should be considered when interpreting the results. Like with all neuroimaging studies, replication of results would further strengthen the assessment of their validity. 5) The sample size of research participants included in this study prevented us from performing additional subgroup analyses, such as differences that might be present based on headache frequency (e.g. episodic migraine vs. chronic migraine), participant sex, use of concurrent migraine preventive medications, and frequency of using migraine as-needed medications. 6) Prior to starting erenumab, a larger proportion of participants who became erenumab responders were concurrently taking other migraine preventive medications. The use of concurrent migraine preventive medications could have an impact on pre-treatment fMRI comparisons between responder and non-responder groups but is unlikely to have impacted differences that were first seen at 2 weeks or 8 weeks following initiation of erenumab. 7) We did not limit the use of medications within 48 hours of the MRI and QST. Thus, some participants had used abortive medications within that time-period. No participants were using opiates, which might directly impact pain sensation during QST and thermal stimulation. 8) QST and fMRI results could be impacted by the participants headache and migraine state, meaning that findings might differ according to whether headache or migraine is present during the test. However, in our study and as presented in the Results, the frequency of testing during headache or migraine was similar between erenumab-responders vs. non-responders and there were few post-treatment QSTs and MRIs collected during headache and none during migraine. Thus, there was not ample justification or sample sizes for analyzing QST and MRI data according to headache and migraine status. 9) Future fMRI studies could use other timepoints for determining treatment response, such as the 9–12-week period after starting erenumab. It is possible that longer durations of response to erenumab could be associated with more substantial changes in fMRI measurements.
Compared to erenumab non-responders, response to erenumab for migraine prevention is associated with post-treatment differences in pain-induced brain activations and resting state functional connectivity. Whether direct or indirect, results suggest that erenumab has effects on brain function, likely impacting central nervous system migraine mechanisms.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Allodynia Symptom Checklist 12
Beck Depression Inventory
calcitonin gene-related peptide
dorsolateral prefrontal cortex
functional magnetic resonance imaging
full width half maximum
general linear model
Migraine Disability Assessment
Montreal Neurological Institute
quantitative sensory testing
region of interest
ventromedial prefrontal cortex
Chong CD, Schwedt TJ, Hougaard A (2019) Brain functional connectivity in headache disorders: a narrative review of MRI investigations. J Cereb Blood Flow Metab 39:650–669
Schwedt TJ, Krauss MJ, Frey K, Gereau RWt. (2011) Episodic and chronic migraineurs are hypersensitive to thermal stimuli between migraine attacks. Cephalalgia 31:6–12
Stankewitz A, Aderjan D, Eippert F, May A (2011) Trigeminal nociceptive transmission in migraineurs predicts migraine attacks. J Neurosci 31:1937–1943
Schulte LH, May A (2016) The migraine generator revisited: continuous scanning of the migraine cycle over 30 days and three spontaneous attacks. Brain 139:1987–1993
Moulton EA, Burstein R, Tully S, Hargreaves R, Becerra L, Borsook D (2008) Interictal dysfunction of a brainstem descending modulatory center in migraine patients. PLoS One 3:e3799
Marciszewski KK, Meylakh N, Di Pietro F et al (2018) Changes in brainstem pain modulation circuitry function over the migraine cycle. J Neurosci 38:10479–10488
Schulte LH, Menz MM, Haaker J, May A (2020) The migraineur's brain networks: continuous resting state fMRI over 30 days. Cephalalgia 40:1614–1621
Krebs K, Rorden C, Androulakis XM (2018) Resting state functional connectivity after sphenopalatine ganglion blocks in chronic migraine with medication overuse headache: a pilot longitudinal fMRI study. Headache 58:732–743
Russo A, Tessitore A, Esposito F et al (2017) Functional changes of the Perigenual part of the anterior cingulate cortex after external trigeminal Neurostimulation in migraine patients. Front Neurol 8:282
Ziegeler C, Mehnert J, Asmussen K, May A (2020) Central effects of erenumab in migraine patients: an event-related functional imaging study. Neurology 95:e2794–ee802
Hebestreit JM, May A (2017) Topiramate modulates trigeminal pain processing in thalamo-cortical networks in humans after single dose administration. PLoS One 12:e0184406
Basedau H, Sturm LM, Mehnert J, Peng KP, Schellong M, May A (2022) Migraine monoclonal antibodies against CGRP change brain activity depending on ligand or receptor target - an fMRI study. Elife 11
Sun H, Dodick DW, Silberstein S et al (2016) Safety and efficacy of AMG 334 for prevention of episodic migraine: a randomised, double-blind, placebo-controlled, phase 2 trial. Lancet Neurol 15:382–390
Goadsby PJ, Reuter U, Hallstrom Y et al (2017) A controlled trial of Erenumab for episodic migraine. N Engl J Med 377:2123–2132
Tepper S, Ashina M, Reuter U et al (2017) Safety and efficacy of erenumab for preventive treatment of chronic migraine: a randomised, double-blind, placebo-controlled phase 2 trial. Lancet Neurol 16:425–434
Noseda R, Schain AJ, Melo-Carrillo A et al (2020) Fluorescently-labeled fremanezumab is distributed to sensory and autonomic ganglia and the dura but not to the brain of rats with uncompromised blood brain barrier. Cephalalgia 40:229–240
Edvinsson JCA, Warfvinge K, Krause DN et al (2019) C-fibers may modulate adjacent Adelta-fibers through axon-axon CGRP signaling at nodes of Ranvier in the trigeminal system. J Headache Pain 20:105
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:338–350
Headache classification Committee of the International Headache Society (IHS) (2018) The international classification of headache disorders, 3rd edition. Cephalalgia 38:1–211
Lipton RB, Bigal ME, Ashina S et al (2008) Cutaneous allodynia in the migraine population. Ann Neurol 63:148–158
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:S20–S28
Beck AT, Steer RA, Brown GK (1996) Manual for Beck depression inventory II. (BDI-II). Pscyhology Corp, San Antonio, TX
Harris PA, Taylor R, Thielke R, Payne J, Gonzalez N, Conde JG (2009) Research electronic data capture (REDCap)--a metadata-driven methodology and workflow process for providing translational research informatics support. J Biomed Inform 42:377–381
Schwedt TJ, Chong CD, Chiang CC, Baxter L, Schlaggar BL, Dodick DW (2014) Enhanced pain-induced activity of pain-processing regions in a case-control study of episodic migraine. Cephalalgia 34:947–958
Upadhyay J, Pendse G, Anderson J et al (2010) Improved characterization of BOLD responses for evoked sensory stimuli. Neuroimage 49:2275–2286
Guillaume B, Hua X, Thompson PM, Waldorp L, Nichols TE (2014) Alzheimer's disease neuroimaging I. fast and accurate modelling of longitudinal and repeated measures neuroimaging data. Neuroimage 94:287–302
Forman SD, Cohen JD, Fitzgerald M, Eddy WF, Mintun MA, Noll DC (1995) Improved assessment of significant activation in functional magnetic resonance imaging (fMRI): use of a cluster-size threshold. Magn Reson Med 33:636–647
Cox RW (1996) AFNI: software for analysis and visualization of functional magnetic resonance neuroimages. Comput Biomed Res 29:162–173
Chong CD, Gaw N, Fu Y, Li J, Wu T, Schwedt TJ (2017) Migraine classification using magnetic resonance imaging resting-state functional connectivity data. Cephalalgia 37:828–844
Duerden EG, Albanese MC (2013) Localization of pain-related brain activation: a meta-analysis of neuroimaging data. Hum Brain Mapp 34:109–149
Maniyar FH, Sprenger T, Monteith T, Schankin C, Goadsby PJ (2014) Brain activations in the premonitory phase of nitroglycerin-triggered migraine attacks. Brain 137:232–241
Mickleborough MJ, Ekstrand C, Gould L et al (2016) Attentional network differences between Migraineurs and non-migraine controls: fMRI evidence. Brain Topogr 29:419–428
Moulton EA, Becerra L, Maleki N et al (2011) Painful heat reveals hyperexcitability of the temporal pole in interictal and ictal migraine states. Cereb Cortex 21:435–448
Peyron R, Laurent B, Garcia-Larrea L (2000) Functional imaging of brain responses to pain. A review and meta-analysis (2000). Neurophysiol Clin 30:263–288
Wager TD, Atlas LY, Lindquist MA, Roy M, Woo CW, Kross E (2013) An fMRI-based neurologic signature of physical pain. N Engl J Med 368:1388–1397
Amin FM, Hougaard A, Magon S et al (2018) Altered thalamic connectivity during spontaneous attacks of migraine without aura: a resting-state fMRI study. Cephalalgia 38:1237–1244
Rubinov M, Sporns O (2010) Complex network measures of brain connectivity: uses and interpretations. Neuroimage 52:1059–1069
He Y, Evans A (2010) Graph theoretical modeling of brain connectivity. Curr Opin Neurol 23:341–350
Lv H, Wang Z, Tong E et al (2018) Resting-state functional MRI: everything that nonexperts have always wanted to know. AJNR Am J Neuroradiol 39:1390–1399
Wang J, Zuo X, He Y (2010) Graph-based network analysis of resting-state functional MRI. Front Syst Neurosci 4:16
Chen JI, Ha B, Bushnell MC, Pike B, Duncan GH (2002) Differentiating noxious- and innocuous-related activation of human somatosensory cortices using temporal analysis of fMRI. J Neurophysiol 88:464–474
Becerra L, Breiter HC, Wise R, Gonzalez RG, Borsook D (2001) Reward circuitry activation by noxious thermal stimuli. Neuron 32:927–946
Frot M, Mauguiere F, Magnin M, Garcia-Larrea L (2008) Parallel processing of nociceptive A-delta inputs in SII and midcingulate cortex in humans. J Neurosci 28:944–952
Kulkarni B, Bentley DE, Elliott R et al (2005) Attention to pain localization and unpleasantness discriminates the functions of the medial and lateral pain systems. Eur J Neurosci 21:3133–3142
Alshelh Z, Marciszewski KK, Akhter R et al (2018) Disruption of default mode network dynamics in acute and chronic pain states. Neuroimage Clin 17:222–231
Baliki MN, Geha PY, Apkarian AV, Chialvo DR (2008) Beyond feeling: chronic pain hurts the brain, disrupting the default-mode network dynamics. J Neurosci 28:1398–1403
Tessitore A, Russo A, Giordano A et al (2013) Disrupted default mode network connectivity in migraine without aura. J Headache Pain 14:89
Coppola G, Di Renzo A, Tinelli E et al (2018) Resting state connectivity between default mode network and insula encodes acute migraine headache. Cephalalgia 38:846–854
Zhang J, Su J, Wang M et al (2016) Increased default mode network connectivity and increased regional homogeneity in migraineurs without aura. J Headache Pain 17:98
Amaral VCG, Tukamoto G, Kubo T, Luiz RR, Gasparetto E, Vincent MB (2018) Migraine improvement correlates with posterior cingulate cortical thickness reduction. Arq Neuropsiquiatr 76:150–157
Starr CJ, Sawaki L, Wittenberg GF et al (2011) The contribution of the putamen to sensory aspects of pain: insights from structural connectivity and brain lesions. Brain 134:1987–2004
Wei HL, Chen J, Chen YC et al (2020) Impaired effective functional connectivity of the sensorimotor network in interictal episodic migraineurs without aura. J Headache Pain 21:111
Chen XY, Chen ZY, Dong Z, Liu MQ, Yu SY (2020) Regional volume changes of the brain in migraine chronification. Neural Regen Res 15:1701–1708
Dominguez C, Lopez A, Ramos-Cabrer P et al (2019) Iron deposition in periaqueductal gray matter as a potential biomarker for chronic migraine. Neurology 92:e1076–e1e85
Solstrand Dahlberg L, Linnman CN, Lee D, Burstein R, Becerra L, Borsook D (2018) Responsivity of periaqueductal gray connectivity is related to headache frequency in episodic migraine. Front Neurol 9:61
Schwedt TJ, Larson-Prior L, Coalson RS et al (2014) Allodynia and descending pain modulation in migraine: a resting state functional connectivity analysis. Pain Med 15:154–165
Mainero C, Boshyan J, Hadjikhani N (2011) Altered functional magnetic resonance imaging resting-state connectivity in periaqueductal gray networks in migraine. Ann Neurol 70:838–845
Li Z, Zhou J, Lan L et al (2020) Concurrent brain structural and functional alterations in patients with migraine without aura: an fMRI study. J Headache Pain 21:141
Stankewitz A, Keidel L, Rehm M et al (2021) Migraine attacks as a result of hypothalamic loss of control. Neuroimage Clin 32:102784
Schulte LH, Mehnert J, May A (2020) Longitudinal neuroimaging over 30 days: temporal characteristics of migraine. Ann Neurol 87:646–651
Veinante P, Yalcin I, Barrot M (2013) The amygdala between sensation and affect: a role in pain. J Mol Psychiatry 1:9
Alves PN, Foulon C, Karolis V et al (2019) An improved neuroanatomical model of the default-mode network reconciles previous neuroimaging and neuropathological findings. Commun Biol 2:370
Li Z, Liu M, Lan L et al (2016) Altered periaqueductal gray resting state functional connectivity in migraine and the modulation effect of treatment. Sci Rep 6:20298
This study was funded by Amgen, Inc. as an investigator sponsored study, ISS 20187183.
Ethics approval and consent to participate
This study was approved by the Mayo Clinic Institutional Review Board. Each research participant reviewed and signed a consent form after participating in the informed consent process. The study was registered at clinicaltrials.gov (NCT03773562).
All authors received funding for their work on this research from an Investigator Sponsored Study grant from Amgen, Inc. Within the prior 24 months: TJS has received personal compensation for consulting with Abbvie, Allergan, Axsome, Biodelivery Science, Biohaven, Click Therapeutics, Collegium, Eli Lilly, Ipsen, Linpharma, Lundbeck, Novartis, Satsuma, and Tonix Pharma. His employer received funds for his consulting work with Amgen. He holds stock options in Aural Analytics and Nocira. He receives royalties from UpToDate. SN has nothing to report. GD has received an investigator-initiated research grant from Amgen, Inc. JL has nothing to report. TW has nothing to report. CDC has nothing to report.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Schwedt, T.J., Nikolova, S., Dumkrieger, G. et al. Longitudinal changes in functional connectivity and pain-induced brain activations in patients with migraine: a functional MRI study pre- and post- treatment with Erenumab. J Headache Pain 23, 159 (2022). https://doi.org/10.1186/s10194-022-01526-5