Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE
PreviousNext
Research Articles, Neurobiology of Disease

Three Dimensions of Association Link Migraine Symptoms and Functional Connectivity

Samuel R. Krimmel, Danielle D. DeSouza, Michael L. Keaser, Bharati M. Sanjanwala, Robert P. Cowan, Martin A. Lindquist, Jennifer A. Haythornthwaite and David A. Seminowicz
Journal of Neuroscience 3 August 2022, 42 (31) 6156-6166; DOI: https://doi.org/10.1523/JNEUROSCI.1796-21.2022
Samuel R. Krimmel
1Department of Neural and Pain Sciences, School of Dentistry, and Center to Advance Chronic Pain Research, University of Maryland, Baltimore, Maryland 21201
2Program in Neuroscience, University of Maryland School of Medicine, Baltimore, Maryland 21201
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Danielle D. DeSouza
3Department of Neurology and Neurological Sciences, Headache and Facial Pain Program, Stanford University, California 94305
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Michael L. Keaser
1Department of Neural and Pain Sciences, School of Dentistry, and Center to Advance Chronic Pain Research, University of Maryland, Baltimore, Maryland 21201
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Bharati M. Sanjanwala
3Department of Neurology and Neurological Sciences, Headache and Facial Pain Program, Stanford University, California 94305
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Robert P. Cowan
3Department of Neurology and Neurological Sciences, Headache and Facial Pain Program, Stanford University, California 94305
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Martin A. Lindquist
4Department Biostatistics, Johns Hopkins Bloomberg School of Public Health, Baltimore, Maryland 21205
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jennifer A. Haythornthwaite
5Department of Psychiatry and Behavioral Sciences, Johns Hopkins University School of Medicine, Baltimore, Maryland 21224
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
David A. Seminowicz
1Department of Neural and Pain Sciences, School of Dentistry, and Center to Advance Chronic Pain Research, University of Maryland, Baltimore, Maryland 21201
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for David A. Seminowicz
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Migraine is a heterogeneous disorder with variable symptoms and responsiveness to therapy. Because of previous analytic shortcomings, variance in migraine symptoms has been inconsistently related to brain function. In the current analysis, we used data from two sites (n = 143, male and female humans), and performed canonical correlation analysis, relating resting-state functional connectivity (RSFC) with a broad range of migraine symptoms, ranging from headache characteristics to sleep abnormalities. This identified three dimensions of covariance between symptoms and RSFC. The first dimension related to headache intensity, headache frequency, pain catastrophizing, affect, sleep disturbances, and somatic abnormalities, and was associated with frontoparietal and dorsal attention network connectivity, both of which are major cognitive networks. Additionally, RSFC scores from this dimension, both the baseline value and the change from baseline to postintervention, were associated with responsiveness to mind-body therapy. The second dimension was related to an inverse association between pain and anxiety, and to default mode network connectivity. The final dimension was related to pain catastrophizing, and salience, sensorimotor, and default mode network connectivity. In addition to performing canonical correlation analysis, we evaluated the current clustering of migraine patients into episodic and chronic subtypes, and found no evidence to support this clustering. However, when using RSFC scores from the three significant dimensions, we identified a novel clustering of migraine patients into four biotypes with unique functional connectivity patterns. These findings provide new insight into individual variability in migraine, and could serve as the foundation for novel therapies that take advantage of migraine heterogeneity.

SIGNIFICANCE STATEMENT Using a large multisite dataset of migraine patients, we identified three dimensions of multivariate association between symptoms and functional connectivity. This analysis revealed neural networks that relate to all measured symptoms, but also to specific symptom ensembles, such as patient propensity to catastrophize painful events. Using these three dimensions, we found four biotypes of migraine informed by clinical and neural variation together. Such findings pave the way for precision medicine therapy for migraine.

  • chronic pain
  • fMRI
  • individual variability
  • mindfulness-based stress reduction
  • networks
  • pain catastrophizing

Introduction

Migraine is a common, poorly managed, and disabling disorder, affecting >10% of the global population (Rasmussen, 1995; Lipton et al., 2001, 2007; Leonardi and Raggi, 2013). Even with improved available treatments (Yuan et al., 2019; Seminowicz et al., 2020), continued research is warranted. Although diagnosis is based on migraine attacks, migraine is a very complex disorder, often featuring sleep, affective, cognitive, and general health abnormalities (Breslau et al., 1994, 2000; Karthik et al., 2012; Lin et al., 2016). Importantly, the presentation of these symptoms is heterogeneous across subjects, which could have important implications for treatment. One approach the field has taken to address clinical variability is to subtype migraine patients. A common subtyping uses a hard cutoff in the frequency of headaches, dividing patients into groups based on higher (chronic migraine) and lower (episodic migraine) frequency (International Headache Society, 2018). Such categorization has important consequences, as treatments can be approved for specific subtypes of migraine and not others (Escher et al., 2017). The current approach is limited for several reasons. First, it is unclear how well episodic versus chronic subtyping captures clinical variability, especially for symptoms other than headache frequency. Second, grouping of patients is based on a limited range of symptoms, ignoring the complex presentation of multiple symptoms in migraine. Third, common migraine subtyping is done independently of biology. Improved analysis of clinical variability, especially when done in conjunction with biological variability, could pave the way for precision medicine migraine therapy, where patients are treated based on unique clinical and biological variables.

The brain can be segregated into functional networks (Yeo et al., 2011), which have been linked to migraine symptoms (Schwedt et al., 2015). A common method for identifying brain networks is resting-state functional connectivity (RSFC), a measure of functional coupling between brain regions (McIntosh, 2000), typically estimated through correlating timeseries from functional units in fMRI acquired during rest (i.e., in the absence of a task or stimulus). Heterogeneity in neural networks could be a useful avenue to explain variability in migraine symptoms, and links between brain function and clinical phenotype would provide potential targets for intervention. Importantly, a better understanding of brain-symptom relationships is critical for developing biologically defined subtypes of illnesses. While studies have examined the relationship between brain function and migraine symptoms, they typically associate RSFC with individual migraine symptoms but do not account for interdependence of symptoms or biological variables (including RSFC).

Canonical correlation analysis (CCA) (Hotelling, 1935) is a multivariate statistical approach that linearly transforms variables from two sets of data to maximize the correlation between them. CCA has become increasingly popular in neuroimaging (H. T. Wang et al., 2020) and has identified dimensions of association between brain connectivity and demographic information (Smith et al., 2015), psychopathology (Drysdale et al., 2017; Xia et al., 2018; Mihalik et al., 2019), and mind-wandering (H. T. Wang et al., 2018). Here we sought to identify dimensions of covariance between diverse migraine symptoms and RSFC. In a large sample from two research centers, we revealed three modes of association between migraine symptoms and RSFC. Unlike standard mass univariate analyses, these three dimensions were each informed by many clinical and RSFC features simultaneously. These dimensions related to symptoms globally and also to more specific combinations of symptoms, and reveal potential targets for brain stimulation or other interventions. Additionally, we found that one of these dimensions was associated with improvements from enhanced mindfulness-based stress reduction (MBSR+), an emerging nonpharmacologic therapy for chronic pain. Finally, we were able to show that the current clustering of migraine into episodic or chronic subtypes based on headache frequency alone is not supported, and propose a new data-driven biotyping of migraine patients into four clusters that is based on RSFC and symptoms together. Our findings offer a new path forward for migraine research where biological and clinical variability is leveraged to provide patient-tailored therapy.

Materials and Methods

Overview

We took clinical measures and RSFC from two study centers; 143 subjects had complete RSFC and clinical data, and were used as input data for the CCA; 166 subjects had complete clinical data and were used to assess the empirical support of episodic versus chronic subtypes of migraine. Multisite harmonization was used to correct for effects of site for RSFC, while still maintaining clinical information. We first applied principal components analysis (PCA) to both clinical and RSFC data to address multicollinearity while reducing dimensionality. Regularized CCA was then performed on components explaining 80% of variance and permutation tests were used to assess statistical significance. Canonical correlation is simply a correlation between canonical variates, where canonical variates are created from two domains of data (clinical and RSFC) so that the canonical correlation is maximized. Stated differently, CCA allows for dimensions of association to be found between clinical scores and RSFC, where these associations are formed through unique combinations of clinical and RSFC variables. Each canonical correlation consists of a clinical and RSFC canonical variate that are created to be associated with one another. For significant canonical correlations, RSFC canonical variates were then tested for association with motion (a major confound for RSFC studies) and clinical improvement from mind-body therapy. Finally, these scores underwent k-means clustering to identify biotypes in our sample (for study overview, see Fig. 1).

Participants

Data from two migraine neuroimaging studies were analyzed, referred to here as the chronic-Stanford and episodic-UMB datasets. In the chronic-Stanford dataset, participants were recruited to participate in research investigating biomarkers of chronic daily headache using clinical, behavioral, and MRI data at Stanford University. All subjects were over 18 years of age and met the International Classification of Headache Disorders ICHD-3 diagnostic criteria for chronic migraine as determined by a physician for a minimum of 3 months (International Headache Society, 2018). Patients were excluded if they had any MRI contraindications or a history of severe neurologic/psychiatric disorders. Subjects provided written informed consent in accordance with Stanford University guidelines. MRI data were collected using a General Electric 3T using an 8-channel phased array head coil. MRI scans used in this analysis were a structural T1-weighted (T1w) scan (TR 5.9 ms, TE minimum, flip angle 15°, voxel size 0.9 × 0.9 × 1 mm), and an eyes closed resting-state functional scan (8 min, gradient echo spiral-pulse, TR 2000 ms, TE 30 ms, slice thickness 4 mm, FOV 220 mm, flip angle 80°, voxel size 3.4 × 3.4 × 4.5 mm). For analyses of clinical data only, subjects were excluded because of missing data, leaving a total of 46 subjects for analysis (average age of 39, SD = 13, 38 female). For MRI analyses, subjects were excluded because of poor MRI quality, excessive motion during resting state (average framewise displacement >0.4 mm), and missing clinical data, leaving a total of 41 subjects for analysis (average age of 38, SD = 13, 34 women).

The second dataset, episodic-UMB, consisted of data from a previously described clinical trial assessing enhanced MBSR+ treatment of migraine (Seminowicz et al., 2020). Briefly, participants were between the ages of 18-65 and met the ICHD-3 criteria for episodic migraine with or without aura for more than 1 year (International Headache Society, 2018). Participants were excluded if they had a history of mindfulness training, reported severe or unstable psychiatric symptoms, and/or used opioids. All patients provided written informed consent in accordance with University of Maryland Baltimore standards. All baseline (pretreatment) MRI scans used a Siemens Tim-Trio 3T MRI scanner with a 32-channel head coil. A scanner upgrade to a Siemens 3T Prisma Fit MRI with a 64-channel head coil occurred in the final year of the clinical trial which affected 1 cohort of 8. MRI scans used in this analysis were a structural T1w scan (TR 2300 ms, TE 2.98 ms, slice thickness 1 mm, FOV 256 mm, flip angle 9°, and voxel size 1 × 1 × 1 mm) and an eyes open resting-state scan (10 min, EPI, TR 2000 ms, TE 28 ms, slice thickness 4 mm, FOV 220 mm, flip angle 77°, and voxel size 3.4 × 3.4 × 4 mm). For analyses of clinical data only, subjects were excluded because of missing data, leaving a total of 120 subjects for analysis (average age of 38, SD = 12, 106 female). For analyses of MRI data, subjects were excluded because of poor MRI quality, excessive motion during resting state (average framewise displacement >0.4 mm), and missing clinical data, for a total of 102 subjects (average age of 38, SD = 12, 91 women); 83 of these subjects were also randomized to a clinical trial receiving either MBSR+ (n = 42) or an active control condition (stress management for headache; n = 41). These randomized participants received additional scans 10 and 20 weeks after the baseline scan (Seminowicz et al., 2020). Two MBSR+ subjects were excluded when analyzing changes in RSFC canonical variate scores leaving 40 subjects.

MRI preprocessing

Preprocessing was performed using FMRIPREP version stable (Esteban et al., 2020), a Nipype (Gorgolewski et al., 2017) based tool. Each T1w volume was corrected for intensity nonuniformity using N4BiasFieldCorrection version 2.1.0 (Tustison et al., 2010) and skull-stripped using antsBrainExtraction.sh version 2.1.0 (using the OASIS template). Spatial normalization to the ICBM 152 Nonlinear Asymmetrical template version 2009c (Fonov et al., 2009) was performed through nonlinear registration with the antsRegistration tool of ANTs version 2.1.0 (Avants et al., 2008), using brain-extracted versions of both T1w volume and template. Brain tissue segmentation of CSF, white matter, and gray matter was performed on the brain-extracted T1w using fast (FSL version 5.0.9) (Zhang et al., 2001).

Functional data were slice time-corrected using 3dTshift from AFNI version 16.2.07 (Cox, 1996) and motion-corrected using mcflirt (FSL version 5.0.9) (Jenkinson et al., 2002). This was followed by coregistration to the corresponding T1w using boundary-based registration (Greve and Fischl, 2009) with 6 degrees of freedom, using flirt (FSL). Motion-correcting transformations, BOLD-to-T1w transformation and T1w-to-template (MNI) warp, were concatenated and applied in a single step using antsApplyTransforms (ANTs version 2.1.0) using Lanczos interpolation.

Framewise displacement (Power et al., 2014) was calculated for each functional run using the implementation of Nipype. ICA-based Automatic Removal Of Motion Artifacts (AROMA) was used to generate aggressive noise regressors as well as to create a variant of data that is nonaggressively denoised (Pruim et al., 2015). Many internal operations of FMRIPREP use Nilearn (Abraham et al., 2014), principally within the BOLD-processing workflow. For more details of the pipeline, see https://fmriprep.readthedocs.io/en/stable/workflows.html.

RSFC processing

Following nonaggressive implementation of AROMA (see MRI preprocessing), resting-state data underwent several additional denoising steps. Resting-state denoising and seed-based analysis were conducted in the CONN functional connectivity toolbox version 19f (Whitfield-Gabrieli and Nieto-Castanon, 2012) (http://www.nitrc.org/projects/conn). Given sustained controversy concerning global signal, we elected not to remove global signal in our analyses (Murphy and Fox, 2017). We performed aCompCor (Behzadi et al., 2007; Muschelli et al., 2014) to remove noise captured in white matter and CSF. We calculated the first five eigenvectors from four erosions of subject-specific white matter masks and two erosions of subject-specific CSF masks, as these levels of erosion no longer contain global signal (Power et al., 2017). Variance in resting-state data related to these 10 anatomic timeseries was removed via linear regression at a voxel-wise level. Given that AROMA has been shown to remove motion artifact from RSFC well relative to other common methods (Ciric et al., 2017), no additional steps were taken to directly remove motion artifact. To avoid noise being reintroduced while removing frequency specific noise, voxel-wise and regressor data were simultaneously bandpass filtered (0.008-0.15 Hz), and all data underwent linear detrending.

Following resting-state functional denoising, timeseries were extracted from 32 ROIs based on an existing independent component analysis of the Human Connectome Project available from the CONN functional connectivity toolbox. These ROIs were organized into default mode, sensorimotor, visual, salience/cingulo-opercular, dorsal attention, frontoparietal, language, and cerebellar networks. To estimate functional connectivity, the 32 ROI timeseries from the entire scanning session were correlated, and the resulting Pearson correlation coefficients were converted to z scores. We excluded functional connections with visual and language network seeds, as we did not believe that these would be associated with the clinical symptoms in our analysis. We do note that some studies have found a relationship between neural abnormalities in migraine and the visual system that are more than effects of aura (Burke et al., 2020). Language network regions included parts of the posterior superior temporal gyrus and the inferior frontal gyrus (consistent with Wernicke's and Broca's area). Visual network regions included nodes in the occipital cortex. While the occipital cortex has been implicated in migraine, it is often associated with symptoms of aura (Hadjikhani et al., 2001), a symptom domain that was not measured in our two datasets. This left a total of 24 ROIs with 276 unique functional connections.

Study site has been shown to influence fMRI data and should be adjusted for when possible (Van Horn and Toga, 2009; Noble et al., 2017; Yu et al., 2018). To adjust for effects of site while controlling for clinical variables, age, and sex, we used Combining Batches of microarray data (ComBat) (Johnson et al., 2007). ComBat was developed to deal with batch effects in high dimensional data and uses an empirical Bayes method to remove batch effects. ComBat was applied following removal of irrelevant ROIs to create site-corrected RSFC, which was used in all analyses, except for when determining changes in canonical variate scores with therapy (see Canonical variate association with clinical outcomes of treatment).

Clinical features

Participants completed multiple self-report questionnaires to assess clinical characteristics. The two headache measures acquired were the frequency of headaches and the average pain associated with headaches. We also examined the exaggerated mental set applied to pain or anticipation of pain using the Pain Catastrophizing Scale (Sullivan et al., 1995). The Pittsburg Sleep Quality Index was used to measure sleep quality as well as sleep disturbances (Buysse et al., 1989). We used the Patient Health Questionnaire 15 to quantify somatization and somatic symptom severity (Kroenke et al., 2002). Finally, we examined affective health using the Generalized Anxiety Disorder-7 measure of anxiety (Spitzer et al., 2006) and the Patient Health Questionnaire 9 measure of depression (Kroenke et al., 2001). The total score for each questionnaire (where applicable) was used for all analyses.

PCA

We used PCA to address multicollinearity and to reduce the dimensionality of our data. PCA creates a set of orthogonal components that can be linearly combined to explain all of the variance in the original dataset. We used singular value decomposition and selected the first four clinical components (see Fig. 2A), explaining >80% of the clinical variance, and the first 35 RSFC components, explaining >80% of the RSFC variance (see Fig. 2B). These components were used as input data for the CCA (see Regularized CCA).

Regularized CCA

CCA is a method for identifying multivariate relationships between sets of variables (Hotelling, 1935). It is used to study the relationship between variables from two domains of interest (e.g., clinical symptoms and biological data) measured on the same set of individuals. In CCA, linear combinations of the input variables from both domains are created (called canonical variates) so that the correlation between the combined variables (canonical correlation, Rc) is maximized. Similar to PCA, it is possible to continue to seek additional pairs of canonical variates that are orthogonal to the preceding canonical variate(s). For each additional canonical variate pair, the corresponding canonical correlation Rc decreases. The total number of possible pairs of canonical variates is given by the smaller of the number of variables in each domain, which in our analysis equaled 4. CCA assumes minimal multicollinearity and is prone to overfitting, especially when the ratio of subjects to input variables is small. We chose CCA over other multivariate approaches, such as multivariate regression, because we had no response-predictor relationship assumed. Several articles provide an excellent introduction with more detail for the interested reader (Tabachnick et al., 2007; Dinga et al., 2019; H. T. Wang et al., 2020). Previous attempts to associate clinical features with RSFC using CCA (Drysdale et al., 2017) may have been severely overfit (Dinga et al., 2019). In our analysis, the starting subject to variable ratio was 0.5, which is likely to overfit, meaning that our results would be unlikely to generalize. Using PCA (described above), we reduced the subject to variable ratio to 3.4. While improved, this ratio was still worse than the generally recommended ratio of 10:1 (Tabachnick et al., 2007; Pituch and Stevens, 2015). To mitigate overfitting, we used an L2 norm regularized CCA (Gonzalez et al., 2008), and optimized regularization parameters (λ) for clinical and RSFC inputs based on held-out testing data from 10-fold cross-validation stratified for study site. The CCA model was built using 9 folds, and the resultant canonical weights were used to create canonical variates in the held-out fold. We repeated this process 10 times; and for each repeat; a new set of 10 folds was created, making in total 100 training sets with unique canonical weights and 100 corresponding testing sets to evaluate performance. For our final model, we used the regularization parameters that yielded the best average Rc value in held-out testing data for the first canonical correlation.

Reconstructing model weights

CCA was performed following PCA of clinical and RSFC data. However, the interpretation of the rotation matrix is challenging, and we therefore transformed the CCA weights back into their original space using the CCA model weights and the rotation matrix. Following this transformation, we were able to use the CCA model weights in the original clinical and RSFC spaces. Applying these weights to clinical and RSFC data yielded clinical and RSFC canonical variates, respectively. We did not use these reconstructed weights directly for interpretation, and instead used canonical loadings for the purposes of interpretation (see Statistics; Fig. 4). We did, however, use these weights to extract canonical variates over time.

Canonical variate association with clinical outcomes of treatment

We examined whether the canonical variates obtained from CCA were associated with outcomes to MBSR+ in the episodic-UMB dataset. To do this, we examined subjects randomized to the enhanced MBSR condition and correlated baseline canonical variates for RSFC with the reduction in headache frequency after 20 weeks, the primary outcome of the trial (Seminowicz et al., 2020). We also sought to determine whether changes in canonical variates were associated with reductions in headache frequency. To do this, we calculated the canonical variates using model weights from the CCA model described above, and calculated difference scores before and after therapy (20 week gap in between scans). In this analysis, we did not use ComBat-corrected RSFC data when determining changes in canonical variate scores, since the data were from a single site and we felt the comparison would be more appropriate when done for data processed in an identical manner. We would like to emphasize though that the CCA model weights used for canonical variate extraction came from ComBat-corrected RSFC data exclusively. The relationship between the change in canonical variate scores and reductions in headache frequency was assessed by means of correlation analysis over the corresponding time period.

Determining support for episodic versus chronic migraine clustering on migraine symptoms

Chronic versus episodic migraine can be thought of as a clustering solution to migraine and is defined based on a 15 headache day per month frequency cutoff. We sought to determine how well this clustering solution was supported by a more holistic catalog of clinical data. We therefore determined the silhouette index for episodic versus chronic distinction using seven clinical variables in the combined sample (n = 166). Should the episodic/chronic clustering have a silhouette index no better than chance (i.e., likely coming from a multivariate normal distribution), then it would argue against the current field norm of episodic and chronic migraine subtypes (see Statistics).

Data-driven clustering of RSFC canonical variates

We sought to identify biotypes in migraine based on the significant RSFC canonical variates from CCA using k-means clustering. We used Euclidean distance as our distance metric for all clustering. We compared 26 clustering metrics using the NBclust package (Charrad et al., 2014) with cluster number (k) ranging from 2 to 10. For each value of k, a score for all 26 clustering metrics was computed, and cluster performance was ranked per metric (e.g., the value of k with the largest silhouette index received a rank of 1). Per metric, a value of k “won” if it received the first rank and therefore had the best performance for that specific clustering metric. We chose a final value for k using a majority rule, where the cluster number that had the most total victories (came in first for a clustering metric) was deemed the best overall. We paid special attention to the silhouette index, as this metric has been shown to accurately identify clusters of data (Petrovic, 2006). This metric compares the relative similarity within a data point's cluster (cohesion) to outside clusters (separation). Larger positive values (maximum of 1) suggest a better clustering solution than values approaching zero.

Statistics

To determine statistical significance of our CCA, we performed a permutation test by constructing 5000 bootstraps with replacement (Efron, 2000) using the final model λ parameters (see Three modes of association were found between RSFC and migraine symptoms). For each bootstrap, we calculated Rc values for four canonical correlations and then compared this null distribution against true Rc values using a one-sided t test for each of the four canonical correlations in a process previously described (Dinga et al., 2019). To interpret CCA results, we calculated canonical loading by correlating canonical variates with their respective model inputs (e.g., RSFC values). We used a false discovery rate multiple comparisons correction on the resultant p values (Benjamini et al., 2001). Determining significance for clustering analysis can be challenging (Liu et al., 2008; Dinga et al., 2019). For our clustering, we tested the true silhouette index scores for the final clustering solution relative to 5000 multivariate Gaussian distributions based on the mean and covariance of the input data using a one-sample t test, in a process similar to ones previously described (Liu et al., 2008; Dinga et al., 2019). This approach allows us to test against the null model of a multivariate normal distribution which does not support clustering. We used logistic regression to predict episodic migraine diagnosis using significant RSFC canonical variates. Based on the predicted versus true diagnosis, a receiver operating characteristic curve was made using the pROC package (Robin et al., 2011) with an area under the curve of 0.86, indicating RSFC canonical variates displayed sensitivity and specificity in separating chronic and episodic migraine. A Spearman rank correlation was used for all correlation analyses unless stated otherwise. For all analyses, significance was determined using a 0.05 false discovery rate-corrected threshold.

Results

Three modes of association were found between RSFC and migraine symptoms

We used CCA to identify dimensions of covariance between a diverse array of clinical symptoms and RSFC from 24 nodes organized into six networks (Fig. 1). We optimized λ values for regularization using performance in held-out testing data to maximize canonical correlation (Rc) for the first pair of canonical variates. The best performance in testing data for the first canonical correlation was an average Rc = 0.6 with λ values of 0.3 for clinical data and 0.9 for RSFC, and we used these regularization parameters for our final model. In the final model, the first canonical correlation was Rc = 0.71, N = 143, p = 0.0002 (Fig. 2C,D), the second was Rc = 0.49, N = 143, p = 0.0023 (Fig. 2E,F), the third was Rc = 0.43, N = 143, p = 0.006 (Fig. 2G,H), and the fourth was Rc = 0.35, N = 143, p = 0.098. The average Rc values for held-out testing data are as follows: first canonical correlation Rc = 0.6, SD = 0.19, second Rc = 0.14, SD = 0.25, third Rc 0.19, SD = 0.23, and fourth Rc = 0.03, SD = 0.23. We focused interpretation and further analysis on the first three significant canonical correlations. Each canonical correlation is between a pair of normally distributed clinical and RSFC canonical variates. Therefore, from the three significant canonical correlations, each subject has three RSFC and three clinical canonical variates.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Study overview. For CCA, 143 migraine subjects were collected from two sites. Seven clinical measures were acquired along with RSFC estimates from 24 regions spanning six functional networks. After site harmonization of RSFC and PCA, CCA was used to associate clinical and biological data. CCA is used to study the relationship between variables from two domains of interest measured in one sample. In CCA, linear combinations of the input variables from both domains are created (called canonical variates) so that the correlation between the combined variables (canonical correlation, Rc) is maximized. RSFC canonical variates were associated with clinical outcomes to mind-body therapy and underwent clustering analysis to reveal biotypes of migraine. Rc, Canonical correlation; CVRSFC, canonical variates for RSFC; CVclinical, canonical variates for clinical data.

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Using CCA to reveal three dimensions of association between clinical features and RSFC. Cumulative variance explained from clinical (A) and RSFC (B) principal components. Red lines indicate 80% variance explained. Canonical correlation (Rc) between RSFC and clinical canonical variates for first (C), second (E), and third (G) canonical correlations. Permutation test for significance of Rc value (shown in red) relative to 5000 permutations for first (D), second (F), and third (H) canonical correlations. CV, Canonical variate.

Association of RSFC canonical variates and scanner motion

Given that motion has profound effects on estimates of RSFC (Power et al., 2012), we ensured that the canonical variates for RSFC were unrelated to motion. We did not find a significant correlation between average framewise displacement during the resting-state scan with the first (rs = 0.12, N = 143, p = 0.14; Fig. 3A), second (rs = 0.08, N = 143, p = 0.3; Fig. 3B), nor third (rs = −0.04, N = 143, p = 0.6; Fig. 3C) RSFC canonical variates. We therefore concluded that there was no evidence for a relationship between RSFC canonical variates and subject movement.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

No association between canonical variates and subject motion. Scatter plot for first (A), second (B), and third (C) RSFC canonical variate relationship with average framewise displacement (FD). CV, Canonical variate.

Interpretation of CCA model

We correlated canonical variates with their respective RSFC/clinical symptoms (i.e., canonical loadings) to interpret the CCA model. While each of the three resultant dimensions consisted of clinical and RSFC canonical variates, we chose to name each dimension with a clinically derived name for ease of discussion. The first canonical variate for symptoms correlated in the same direction with all seven clinical features (Fig. 4A). Hence, we interpret the first clinical canonical correlation as the “global symptom dimension.” Patients with positive clinical canonical variate values for the first canonical correlation had globally better symptoms (less severe symptoms overall). The associated first RSFC canonical variate primarily reflected frontoparietal network (FPN) and dorsal attention network (DAN) connectivity (Fig. 4B). Patients with positive RSFC canonical variate values had less FPN and DAN connectivity than patients with negative RSFC canonical variate scores.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Interpretation of canonical variates. Canonical variates were correlated with input data to facilitate interpretation. A, For clinical data, correlations shown for the three identified dimensions. Correlation between RSFC and global dimension (B), inverse pain/anxiety dimension (C), and pain catastrophizing dimension (D) scores organized over 24 ROIs in six networks. The global dimension was most associated with FPN and DAN connectivity with other networks. The inverse pain/anxiety dimension was most related to DMN and sensorimotor connectivity. The pain catastrophizing dimension was most related to salience, sensorimotor, and DMN connectivity. The same color coding for networks and correlation color range is used in B–D. All off-white correlation values are significant for all four panels. Mean Pain, Average headache pain; NHA, headache frequency; PHQ-15, Patient Healthy Questionnaire 15 item; PSQI, Pittsburgh Sleep Quality Index; PCS, Pain Catastrophizing Scale; GAD-7, Generalized Anxiety Disorder; PHQ-9, Patient Health Questionnaire 9 item; r, Pearson correlation coefficient.

The second clinical canonical variate reflected an inverse relationship between average headache pain and anxiety primarily, referred to here as the “inverse pain/anxiety dimension” (Fig. 4A). Patients with positive clinical canonical variate values had above-average headache pain and below-average anxiety. The pattern of functional connectivity was quite sparse for this but was primarily related to anticorrelation between the default mode network (DMN) and other networks (Fig. 4C). Patients with positive RSFC canonical variate values had increased DMN connectivity relative to negative scoring patients. The third symptom canonical variate primarily related to pain catastrophizing (Fig. 4A) and related to salience, sensorimotor, and default mode connectivity with one another (Fig. 4D). Patients with positive clinical values had below-average pain catastrophizing. Patients with positive RSFC canonical variate scores had above-average sensorimotor connectivity with DMN and salience network. We refer to this here as a “pain catastrophizing dimension.”

RSFC canonical variates are associated with headache reduction in MBSR+

We sought to determine whether RSFC canonical variates were associated with improvements from mindfulness therapy for migraine. We correlated baseline canonical variates from functional connectivity with the reduction in frequency of headaches and found a significant association for the global dimension (rs = −0.44, N = 42, p = 0.0038) (Fig. 5). RSFC canonical variates from neither the inverse anxiety/pain nor the pain catastrophizing dimensions were significantly associated with clinical improvement. Given that baseline global dimension RSFC canonical variate associated with clinical improvements, we next examined whether changes in this dimension associated with headache frequency reduction following MBSR+ in the episodic-UMB cohort and found a significant association (rs= −0.40, N = 40, p = 0.014).

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

Association between baseline (pre-intervention) global symptom RSFC canonical variate and change in headache frequency after therapy (positive numbers indicate a reduction in headache days). Subjects with lower scores and more severe symptoms overall had better responses to MBSR+ than subjects with less severe symptoms. CV, Canonical variate.

Testing the support for current migraine clustering

We tested whether the current clustering of migraine patients into episodic (<15 headaches/mo) and chronic (≥15 headaches/mo) subtypes differs from a multivariate Gaussian distribution. Chronic and episodic migraine patients were collected at unique sites, creating a potential confound (see Discussion). Using the seven clinical features used for CCA, we determined the silhouette index in 166 patients using the chronic/episodic distinction, yielding a value of 0.23. We then again compared this with 5000 multivariate Gaussian distribution based on the mean and covariance of the seven clinical measures. This analysis indicated that the chronic versus episodic clustering was not significantly different from a multivariate Gaussian distribution (N = 166, p = 0.66) (Fig. 6A). Therefore, there is insufficient evidence to reject the null hypothesis that the episodic versus chronic clustering solution comes from a multivariate Gaussian distribution. This means that using a holistic catalog of symptoms, we were unable to find support for subtyping migraine patients into episodic and chronic migraine categories.

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

Lack of support for episodic/chronic migraine clustering and identification of a 4 biotype clustering of migraine. A, Lack of significance for clustering migraine into chronic and episodic subtypes based on silhouette index over 5000 multivariate Gaussian distributions of clinical data. B, Silhouette index for various cluster numbers (k) from three RSFC canonical variates indicating that k = 4 had the best performance. C, Significance of clustering into four biotypes based on silhouette index over 5000 multivariate Gaussian distributions of RSFC canonical variate scores. D, Patterns of RSFC canonical variates for three dimensions over the four identified biotypes. Color coding represents the dimension. Canonical variates are distributed with a mean of zero; therefore, negative values indicate below-average scores and positive values indicate above-average scores. CV, Canonical variate; k, number of clusters.

Identifying biotypes in migraine

In order to identify subtypes of migraine based on biological data, called biotypes, we performed k-means clustering based on the canonical variates for RSFC for the three dimensions. The best performing solution was k = 4, which had the best performance in 27% of metrics, although there were other solutions with similar performance (k = 2 was best in 23%, no other solutions were consistently as strong; for silhouette index scores over various values of k, see Fig. 6B). To determine significance of a four cluster solution, we estimated 5000 multivariate Gaussian distributions based on the mean and covariance of the three RSFC canonical variates. We calculated the silhouette index, a measure of clustering quality, in each of these 5000 distributions and compared the true four cluster silhouette index of 0.28 against them (Fig. 6C). The true silhouette index of 0.28 was significant relative to the multivariate Gaussian distributions (N = 143, p = 0.049) (Fig. 6C). These clusters had unique patterns of RSFC canonical variate scores (Fig. 6D) and clinical data (Table 1). Cluster 1 had negative and therefore below-average RSFC canonical variate values for the global symptom dimension, indicating patients in this group had worse than average symptoms overall. Cluster 2 had below-average RSFC canonical variate scores for the pain catastrophizing dimension, indicating patients in this group had above-average catastrophizing. Cluster 4 had below-average pain/anxiety inverse dimension scores, indicating patients in this group had below-average mean pain and above-average anxiety. Cluster 3 had above-average scores in all three dimensions, indicating patients in this group had the opposite pattern of symptoms as all other biotypes.

View this table:
  • View inline
  • View popup
Table 1.

Clinical profiles for migraine biotypesa

Discussion

Migraine is a heterogeneous disorder with variable response to therapy and expression of symptoms. This variability has been ignored in many studies that treat migraineurs as a single group in between-subjects analyses contrasting patients and healthy controls. When variability has been factored in, subjects are often subtyped into categories not derived from biological data that have also not been rigorously tested. When studies do perform within-subjects analyses to associate symptoms with biological data, they often examine only one variable of interest, ignoring the complex relationship between migraine symptoms. In the current study, we addressed these existing limitations by using a large sample size, including a wealth of clinical data from well-phenotyped subjects, and by relating these clinical profiles to RSFC in a single multivariate statistical step, as opposed to a standard mass univariate analysis. This approach identified three dimensions of covariation between clinical and RSFC data. The first was a global-symptom dimension that was associated with all symptoms and was primarily related to FPN and DAN functional connectivity with other networks. The FPN includes the posterior parietal cortex and the lateral prefrontal cortex and is consistently associated with cognitive control (Bush and Shin, 2006; Scolari et al., 2015) as well as acute and chronic pain (Seminowicz et al., 2011; Seminowicz and Moayedi, 2017; Y. Wang et al., 2017). The DAN contains the intraparietal sulcus and the frontal eye fields and is consistently involved in attention (Vossel et al., 2014). Our findings suggest that cognitive networks relate to symptoms generally, and that they may function as targets to simultaneously treat all symptoms at once. Consistent with this, the normalization of FPN functional connectivity has been associated with clinical improvement in migraine (Li et al., 2015). The second dimension reflected inverse pain and anxiety scores and was related to default mode (DMN) and sensorimotor network connectivity. The DMN contains the angular gyrus, posterior cingulate cortex, and the mPFC; and the sensorimotor network contains the precentral and postcentral gyri. Increasing or decreasing DMN-sensorimotor connectivity might therefore be a fruitful avenue to modify anxiety and pain, although the inverse relationship between these two clinical variables in the second canonical correlation suggests that such targeting would be side effect prone. Finally, the third dimension was a pain catastrophizing dimension and was primarily associated with functional connectivity of salience network (anterior insula and midcingulate cortex) to DMN and sensorimotor networks, and to DMN connectivity with other networks. Previous literature has shown a relationship between the salience network and pain catastrophizing (Hubbard et al., 2014), and our results suggest that targeting this network might prove to be an especially useful tool for targeting pain catastrophizing.

We were able to show that the RSFC scores for the global symptom dimension were associated with improvements from the nonpharmacologic MBSR+ treatment of migraine, which has recently been established efficacy for migraine (Seminowicz et al., 2020). Subjects with lower scores and thus more severe symptoms showed a larger reduction in headache frequency following MBSR+ than subjects who had comparatively less severe symptoms. In other words, patients with more severe symptoms overall experienced more pronounced improvements from MBSR+. Additionally, we observed that changes in global dimension RSFC scores toward the positive pole (less symptom severe) were associated with improvements from MBSR+ as well. This might indicate that reducing FPN connectivity to salience, sensorimotor, and DANs could enhance the effects of MBSR+, given that the first dimension was negatively associated with this FPN functional connectivity.

The migraine RSFC literature uses small sample sizes and inconsistent analyses, making it difficult to synthesize (Skorobogatykh et al., 2019). However, there do appear to be important points of convergence between our findings and previous work, which further supports the validity of our findings. Abnormal connectivity in migraine is often observed in the FPN, salience, sensorimotor, DMN, and DAN (Russo et al., 2012; Xue et al., 2012; Jin et al., 2013; Hubbard et al., 2014; Ke et al., 2020). We therefore anticipated that these networks would be associated with illness characteristics. Two clinical variables included in our analysis that have been the subject of multiple studies in migraine are headache frequency and pain catastrophizing. Frequency was most related to our first identified global dimension, and associations between headache frequency and functional connectivity of salience and FPN nodes (Mainero et al., 2011; Maleki et al., 2012; Hubbard et al., 2014). Pain catastrophizing is often linked to attentional processes (Quartana et al., 2009), and Pain Catastrophizing Scale scores are related to anterior insula (key node of the salience network) and sensorimotor connectivity in migraine patients (Hubbard et al., 2014).

Our study has major implications to migraine research and clinical trial design. Migraine is currently split into episodic and chronic classifications based on a 15 headache per month cutoff (International Headache Society, 2018). There have been criticisms of the concept of chronic migraine in the literature (Medrea and Christi, 2018), and in particular over the exact cutoff point between episodic and chronic patients (Torres-Ferrús et al., 2017). We tested the quality of the current two cluster solution for migraine patients in our data based on seven clinical features, finding no evidence against the null hypothesis of a multivariate normal distribution (which does not have clusters). This is especially interesting because the episodic and chronic migraine patient cohorts in our study were collected at two separate sites, which if anything should bias the results toward a chronic versus episodic clustering solution. Additionally, we identified four biotypes of migraine that could be used for predicting therapeutic outcomes in future studies. Patients falling into a given cluster had unique brain-symptom relationships arguing for targeted therapy (Fig. 6D). For instance, patients in Cluster 1 had the most extreme global symptom dimension score, which was associated with FPN and DAN connectivity, and perhaps these patients would benefit the most from interventions targeting these networks. On the other hand, Cluster 2 patients who had more extreme pain catastrophizing might benefit the most from interventions that reduce catastrophizing and also target sensorimotor and salience networks. Together, these clustering results argue against the current episodic/chronic clustering that dominates migraine research and clinical trial design, and suggest an alternative based on clinical and biological association. Future work could build on this biotyping by using multimodal data while also examining the potential benefits of different clustering approaches.

This study is not without limitations. A major danger when conducting CCA is overfitting, and previously exciting results in depression may have fallen victim to this issue (Dinga et al., 2019). In the current study, we took several steps to mitigate overfitting. First, we reduced the number of features in the model through parcellations (instead of voxels) for RSFC, feature selection, and PCA. Together, this allowed us to greatly increase the ratio of subjects to model features, attenuating overfitting. We took the additional step of performing regularized CCA and tuned λ parameters based on the held-out strength of the first canonical correlation. This resulted in very minimal overfitting for the first canonical correlation. However, overfitting was larger for the second and third canonical correlations, indicating that some caution is warranted when interpreting and using them. Future work could expand on our results by testing this model in new data and benefit from data sharing to create larger sample sizes that are less prone to overfitting. Additionally, statistical inference for clustering is challenging (Liu et al., 2008; Dinga et al., 2019). In the current study, we tested against a null distribution of a multivariate Gaussian distribution, which when true, would not merit clustering. However, even with rejection of this null hypothesis, it is possible that clustering is still unwarranted because not all multivariate non-Gaussian data are necessarily organized into clusters. Future work could examine the utility of this clustering solution by using the four biotypes to predict therapeutic outcomes from migraine treatment, and by examining whether these biotypes exist in additional chronic pain disorders. An additional limitation comes from potential site effects because of unique populations, MR systems, and scanning protocol decisions (open or closed eyes during resting state). While we reduced effects of study site by excluding certain ROIs and performing multisite harmonization, site effects are potentially still present in RSFC data.

In conclusion, we identify an association between migraine symptoms (including headache severity, quality of life, affective measures, and coping) and whole-brain RSFC, yielding three dimensions of association between these two domains. These results may facilitate the development of personalized medicine, which is limited by treating illnesses as homogeneous groupings. Moreover, the biological association with our clinical data identifies potential targets for therapy and research. Finally, our data argue against the current clustering of migraine patients into chronic and episodic classifications and instead offer an alternative that is grounded in clinical presentation and biology simultaneously.

Footnotes

  • This work was supported by National Center for Complementary and Integrative Health/National Institutes of Health R01 AT007176 to D.A.S.; the Wings Foundation; the Wintercreek Foundation; and the Higgins Family Trust to B.M.S. and R.P.C.

  • D.A.S. is an advisor to Empower Therapeutics. The remaining authors declare no competing financial interests.

  • Correspondence should be addressed to David A. Seminowicz at dseminow{at}uwo.ca

SfN exclusive license.

References

  1. ↵
    1. Abraham A,
    2. Pedregosa F,
    3. Eickenberg M,
    4. Gervais P,
    5. Mueller A,
    6. Kossaifi J,
    7. Gramfort A,
    8. Thirion B,
    9. Varoquaux G
    (2014) Machine learning for neuroimaging with scikit-learn. Front Neuroinform 8:14. doi:10.3389/fninf.2014.00014 pmid:24600388
    OpenUrlCrossRefPubMed
  2. ↵
    1. Avants BB,
    2. Epstein CL,
    3. Grossman M,
    4. Gee JC
    (2008) Symmetric diffeomorphic image registration with cross-correlation: evaluating automated labeling of elderly and neurodegenerative brain. Med Image Anal 12:26–41. doi:10.1016/j.media.2007.06.004 pmid:17659998
    OpenUrlCrossRefPubMed
  3. ↵
    1. Behzadi Y,
    2. Restom K,
    3. Liau J,
    4. Liu TT
    (2007) A component based noise correction method (CompCor) for BOLD and perfusion based fMRI. Neuroimage 37:90–101. pmid:17560126
    OpenUrlCrossRefPubMed
  4. ↵
    1. Benjamini Y,
    2. Drai D,
    3. Elmer G,
    4. Kafkafi N,
    5. Golani I
    (2001) Controlling the false discovery rate in behavior genetics research. Behav Brain Res 125:279–284. pmid:11682119
    OpenUrlCrossRefPubMed
  5. ↵
    1. Breslau N,
    2. Davis GC,
    3. Schultz LR,
    4. Paterson EL
    (1994) Migraine and major depression: a longitudinal study. Headache 34:387–393. doi:10.1111/j.1526-4610.1994.hed3407387.x
    OpenUrlCrossRefPubMed
  6. ↵
    1. Breslau N,
    2. Schultz LR,
    3. Stewart WF,
    4. Lipton RB,
    5. Lucia VC,
    6. Welch KM
    (2000) Headache and major depression: is the association specific to migraine? Neurology 54:308–308. pmid:10668688
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Burke MJ,
    2. Joutsa J,
    3. Cohen AL,
    4. Soussand L,
    5. Cooke D,
    6. Burstein R,
    7. Fox MD
    (2020) Mapping migraine to a common brain network. Brain 143:541–553. doi:10.1093/brain/awz405
    OpenUrlCrossRef
  8. ↵
    1. Bush G,
    2. Shin LM
    (2006) The Multi-Source Interference Task: an fMRI task that reliably activates the cingulo-frontal-parietal cognitive/attention network. Nat Protoc 1:308–313. doi:10.1038/nprot.2006.48 pmid:17406250
    OpenUrlCrossRefPubMed
  9. ↵
    1. Buysse DJ,
    2. Reynolds CF,
    3. Monk TH,
    4. Berman SR,
    5. Kupfer DJ
    (1989) The Pittsburgh Sleep Quality Index: a new instrument for psychiatric practice and research. Psychiatry Res 28:193–213. doi:10.1016/0165-1781(89)90047-4 pmid:2748771
    OpenUrlCrossRefPubMed
  10. ↵
    1. Charrad M,
    2. Ghazzali N,
    3. Boiteau V,
    4. Niknafs A
    (2014) NbClust: An R package for determining the relevant number of clusters in a data set. J Stat Softw 61:1–36.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Ciric R,
    2. Wolf DH,
    3. Power JD,
    4. Roalf DR,
    5. Baum GL,
    6. Ruparel K,
    7. Shinohara RT,
    8. Elliott MA,
    9. Eickhoff SB,
    10. Davatzikos C,
    11. Gur RC,
    12. Gur RE,
    13. Bassett DS,
    14. Satterthwaite TD
    (2017) Benchmarking of participant-level confound regression strategies for the control of motion artifact in studies of functional connectivity. Neuroimage 154:174–187. pmid:28302591
    OpenUrlCrossRefPubMed
  12. ↵
    1. Cox RW
    (1996) AFNI: software for analysis and visualization of functional magnetic resonance neuroimages. Comput Biomed Res 29:162–173. doi:10.1006/cbmr.1996.0014 pmid:8812068
    OpenUrlCrossRefPubMed
  13. ↵
    1. Dinga R,
    2. Schmaal L,
    3. Penninx BW,
    4. van Tol MJ,
    5. Veltman DJ,
    6. van Velzen L,
    7. Mennes M,
    8. van der Wee NJ,
    9. Marquand AF
    (2019) Evaluating the evidence for biotypes of depression: methodological replication and extension of Drysdale et al. (2017). Neuroimage Clin 22:101796. doi:10.1016/j.nicl.2019.101796 pmid:30935858
    OpenUrlCrossRefPubMed
  14. ↵
    1. Drysdale AT, et al
    . (2017) Resting-state connectivity biomarkers define neurophysiological subtypes of depression. Nat Med 23:28–38. doi:10.1038/nm.4246 pmid:27918562
    OpenUrlCrossRefPubMed
  15. ↵
    1. Efron B
    (2000) The bootstrap and modern statistics. J Am Stat Assoc 95:1293–1296. doi:10.1080/01621459.2000.10474333
    OpenUrlCrossRef
  16. ↵
    1. Escher CM,
    2. Paracka L,
    3. Dressler D,
    4. Kollewe K
    (2017) Botulinum toxin in the management of chronic migraine: clinical evidence and experience. Ther Adv Neurol Disord 10:127–135. pmid:28382110
    OpenUrlPubMed
  17. ↵
    1. Esteban O,
    2. Markiewicz CJ,
    3. Goncalves M,
    4. DuPre E,
    5. Kent JD,
    6. Salo T,
    7. Ciric R,
    8. Pinsard B,
    9. Blair RW,
    10. Poldrack RA,
    11. Gorgolewski KJ
    (2020) fMRIPrep: a robust preprocessing pipeline for functional MRI. Zenodo. Available at https://doi.org/10.5281/zenodo.4055773.
  18. ↵
    1. Fonov V,
    2. Evans A,
    3. McKinstry R,
    4. Almli C,
    5. Collins D
    (2009) Unbiased nonlinear average age-appropriate brain templates from birth to adulthood. Neuroimage 47:S102. doi:10.1016/S1053-8119(09)70884-5
    OpenUrlCrossRef
  19. ↵
    1. Gonzalez I,
    2. Déjean S,
    3. Martin P,
    4. Baccini A
    (2008) CCA: an R package to extend canonical correlation analysis. J Stat Soft 23:1–14. doi:10.18637/jss.v023.i12
    OpenUrlCrossRef
  20. ↵
    1. Gorgolewski KJ, et al
    . (2017) Nipype: a flexible, lightweight and extensible neuroimaging data processing framework in Python. 0.13.1. Zenodo. Available at https://doi.org/10.5281/zenodo.581704.
  21. ↵
    1. Greve DN,
    2. Fischl B
    (2009) Accurate and robust brain image alignment using boundary-based registration. Neuroimage 48:63–72. doi:10.1016/j.neuroimage.2009.06.060 pmid:19573611
    OpenUrlCrossRefPubMed
  22. ↵
    1. Hadjikhani N,
    2. Rio MS,
    3. Wu O,
    4. Schwartz D,
    5. Bakker D,
    6. Fischl B,
    7. Kwong KK,
    8. Cutrer FM,
    9. Rosen BR,
    10. Tootell RB,
    11. Sorensen AG,
    12. Moskowitz MA
    (2001) Mechanisms of migraine aura revealed by functional MRI in human visual cortex. Proc Natl Acad Sci USA 98:4687–4692. doi:10.1073/pnas.071582498 pmid:11287655
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Hotelling H
    (1935) Canonical correlation analysis (CCA). J Educ Psychol
  24. ↵
    1. Hubbard CS,
    2. Khan SA,
    3. Keaser ML,
    4. Mathur VA,
    5. Goyal M,
    6. Seminowicz DA
    (2014) Altered brain structure and function correlate with disease severity and pain catastrophizing in migraine patients. eNeuro 1:ENEURO.0006-14.2014. doi:10.1523/ENEURO.0006-14.2014
    OpenUrlCrossRefPubMed
  25. ↵
    International Headache Society (2018) Headache Classification Committee of the International Headache Society (IHS): The International Classification of Headache Disorders, 3rd edition. Cephalalgia 38:1–211.
    OpenUrlCrossRefPubMed
  26. ↵
    1. Jenkinson M,
    2. Bannister P,
    3. Brady M,
    4. Smith S
    (2002) Improved optimization for the robust and accurate linear registration and motion correction of brain images. Neuroimage 17:825–841. doi:10.1006/nimg.2002.1132 pmid:12377157
    OpenUrlCrossRefPubMed
  27. ↵
    1. Jin C,
    2. Yuan K,
    3. Zhao L,
    4. Zhao L,
    5. Yu D,
    6. von Deneen KM,
    7. Zhang M,
    8. Qin W,
    9. Sun W,
    10. Tian J
    (2013) Structural and functional abnormalities in migraine patients without aura. NMR Biomed 26:58–64. pmid:22674568
    OpenUrlCrossRefPubMed
  28. ↵
    1. Johnson WE,
    2. Li C,
    3. Rabinovic A
    (2007) Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics 8:118–127. doi:10.1093/biostatistics/kxj037 pmid:16632515
    OpenUrlCrossRefPubMed
  29. ↵
    1. Karthik N,
    2. Kulkarni GB,
    3. Taly AB,
    4. Rao S,
    5. Sinha S
    (2012) Sleep disturbances in 'migraine without aura': a questionnaire based study. J Neurol Sci 321:73–76. pmid:22877507
    OpenUrlCrossRefPubMed
  30. ↵
    1. Ke J,
    2. Yu Y,
    3. Zhang X,
    4. Su Y,
    5. Wang X,
    6. Hu S,
    7. Dai H,
    8. Hu C,
    9. Zhao H,
    10. Dai L
    (2020) Functional alterations in the posterior insula and cerebellum in migraine without aura: a resting-state MRI study. Front Behav Neurosci 14:567588. doi:10.3389/fnbeh.2020.567588 pmid:33132860
    OpenUrlCrossRefPubMed
  31. ↵
    1. Kroenke K,
    2. Spitzer RL,
    3. Williams JB
    (2001) The PHQ-9: validity of a brief depression severity measure. J Gen Intern Med 16:606–613. pmid:11556941
    OpenUrlCrossRefPubMed
  32. ↵
    1. Kroenke K,
    2. Spitzer RL,
    3. Williams JB
    (2002) The PHQ-15: validity of a new measure for evaluating the severity of somatic symptoms. Psychosom Med 64:258–266. pmid:11914441
    OpenUrlAbstract/FREE Full Text
  33. ↵
    1. Leonardi M,
    2. Raggi A
    (2013) Burden of migraine: international perspectives. Neurol Sci 34 Suppl 1:S117–S118. doi:10.1007/s10072-013-1387-8 pmid:23695058
    OpenUrlCrossRefPubMed
  34. ↵
    1. Li K,
    2. Zhang Y,
    3. Ning Y,
    4. Zhang H,
    5. Liu H,
    6. Fu C,
    7. Ren Y,
    8. Zou Y
    (2015) The effects of acupuncture treatment on the right frontoparietal network in migraine without aura patients. J Headache Pain 16:33. doi:10.1186/s10194-015-0518-4
    OpenUrlCrossRef
  35. ↵
    1. Lin YK,
    2. Lin GY,
    3. Lee JT,
    4. Lee MS,
    5. Tsai CK,
    6. Hsu YW,
    7. Lin YZ,
    8. Tsai YC,
    9. Yang FC
    (2016) Associations between sleep quality and migraine frequency: a cross-sectional case-control study. J Headache Pain 17:97. doi:10.1097/MD.0000000000003554 pmid:27770406
    OpenUrlCrossRefPubMed
  36. ↵
    1. Lipton RB,
    2. Stewart WF,
    3. Diamond S,
    4. Diamond ML,
    5. Reed M
    (2001) Prevalence and burden of migraine in the United States: data from the American Migraine Study II. Headache 41:646–657. doi:10.1046/j.1526-4610.2001.041007646.x pmid:11554952
    OpenUrlCrossRefPubMed
  37. ↵
    1. Lipton RB,
    2. Bigal ME,
    3. Diamond M,
    4. Freitag F,
    5. Reed ML,
    6. Stewart WF
    , AMPP Advisory Group (2007) Migraine prevalence, disease burden, and the need for preventive therapy. Neurology 68:343–349. pmid:17261680
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Liu Y,
    2. Hayes DN,
    3. Nobel A,
    4. Marron JS
    (2008) Statistical significance of clustering for high-dimension, low–sample size data. J Am Stat Assoc 103:1281–1293. doi:10.1198/016214508000000454
    OpenUrlCrossRef
  39. ↵
    1. Mainero C,
    2. Boshyan J,
    3. Hadjikhani N
    (2011) Altered functional magnetic resonance imaging resting-state connectivity in periaqueductal gray networks in migraine. Ann Neurol 70:838–845. pmid:22162064
    OpenUrlCrossRefPubMed
  40. ↵
    1. Maleki N,
    2. Becerra L,
    3. Brawn J,
    4. Bigal M,
    5. Burstein R,
    6. Borsook D
    (2012) Concurrent functional and structural cortical alterations in migraine. Cephalalgia 32:607–620. doi:10.1177/0333102412445622
    OpenUrlCrossRefPubMed
  41. ↵
    1. McIntosh AR
    (2000) Towards a network theory of cognition. Neural Netw 13:861–870. doi:10.1016/S0893-6080(00)00059-9
    OpenUrlCrossRefPubMed
  42. ↵
    1. Medrea I,
    2. Christi S
    (2018) Chronic migraine: evolution of the concept and clinical implications. Headache 58:1495–1500. doi:10.1111/head.13380 pmid:30178876
    OpenUrlCrossRefPubMed
  43. ↵
    1. Mihalik A, et al
    . (2019) Brain-behaviour modes of covariation in healthy and clinically depressed young people. Sci Rep 9:1–11. doi:10.1038/s41598-019-47277-3
    OpenUrlCrossRefPubMed
  44. ↵
    1. Murphy K,
    2. Fox MD
    (2017) Towards a consensus regarding global signal regression for resting state functional connectivity MRI. Neuroimage 154:169–173. pmid:27888059
    OpenUrlCrossRefPubMed
  45. ↵
    1. Muschelli J,
    2. Nebel MB,
    3. Caffo BS,
    4. Barber AD,
    5. Pekar JJ,
    6. Mostofsky SH
    (2014) Reduction of motion-related artifacts in resting state fMRI using aCompCor. Neuroimage 96:22–35. pmid:24657780
    OpenUrlCrossRefPubMed
  46. ↵
    1. Noble S, et al
    . (2017) Multisite reliability of MR-based functional connectivity. Neuroimage 146:959–970. doi:10.1016/j.neuroimage.2016.10.020 pmid:27746386
    OpenUrlCrossRefPubMed
  47. ↵
    1. Petrovic S
    (2006) A comparison between the silhouette index and the Davies-Bouldin Index in labelling IDS clusters. Proceedings of the 11th Nordic Workshop of Secure IT Systems, pp 53–64.
  48. ↵
    1. Pituch KA,
    2. Stevens JP
    (2015) Applied multivariate statistics for the social sciences: analyses with SAS and IBM's SPSS. Milton Park, UK: Routledge.
  49. ↵
    1. Power JD,
    2. Barnes KA,
    3. Snyder AZ,
    4. Schlaggar BL,
    5. Petersen SE
    (2012) Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. Neuroimage 59:2142–2154. doi:10.1016/j.neuroimage.2011.10.018 pmid:22019881
    OpenUrlCrossRefPubMed
  50. ↵
    1. Power JD,
    2. Mitra A,
    3. Laumann TO,
    4. Snyder AZ,
    5. Schlaggar BL,
    6. Petersen SE
    (2014) Methods to detect, characterize, and remove motion artifact in resting state fMRI. Neuroimage 84:320–341. doi:10.1016/j.neuroimage.2013.08.048 pmid:23994314
    OpenUrlCrossRefPubMed
  51. ↵
    1. Power JD,
    2. Plitt M,
    3. Laumann TO,
    4. Martin A
    (2017) Sources and implications of whole-brain fMRI signals in humans. Neuroimage 146:609–625. pmid:27751941
    OpenUrlCrossRefPubMed
  52. ↵
    1. Pruim RH,
    2. Mennes M,
    3. van Rooij D,
    4. Llera A,
    5. Buitelaar JK,
    6. Beckmann CF
    (2015) ICA-AROMA: a robust ICA-based strategy for removing motion artifacts from fMRI data. Neuroimage 112:267–277. doi:10.1016/j.neuroimage.2015.02.064 pmid:25770991
    OpenUrlCrossRefPubMed
  53. ↵
    1. Quartana PJ,
    2. Campbell CM,
    3. Edwards RR
    (2009) Pain catastrophizing: a critical review. Expert Rev Neurother 9:745–758. pmid:19402782
    OpenUrlCrossRefPubMed
  54. ↵
    1. Rasmussen BK
    (1995) Epidemiology of migraine. Biomed Pharmacother 49:452–455. doi:10.1016/0753-3322(96)82689-8
    OpenUrlCrossRefPubMed
  55. ↵
    1. Robin X,
    2. Turck N,
    3. Hainard A,
    4. Tiberti N,
    5. Lisacek F,
    6. Sanchez JC,
    7. Müller M
    (2011) PROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics 12:1–8. doi:10.1186/1471-2105-12-77
    OpenUrlCrossRefPubMed
  56. ↵
    1. Russo A,
    2. Tessitore A,
    3. Giordano A,
    4. Corbo D,
    5. Marcuccio L,
    6. De Stefano M,
    7. Salemi F,
    8. Conforti R,
    9. Esposito F,
    10. Tedeschi G
    (2012) Executive resting-state network connectivity in migraine without aura. Cephalalgia 32:1041–1048. pmid:22908362
    OpenUrlCrossRefPubMed
  57. ↵
    1. Schwedt TJ,
    2. Chiang CC,
    3. Chong CD,
    4. Dodick DW
    (2015) Functional MRI of migraine. Lancet Neurol 14:81–91. pmid:25496899
    OpenUrlCrossRefPubMed
  58. ↵
    1. Scolari M,
    2. Seidl-Rathkopf KN,
    3. Kastner S
    (2015) Functions of the human frontoparietal attention network: evidence from neuroimaging. Curr Opin Behav Sci 1:32–39. pmid:27398396
    OpenUrlCrossRefPubMed
  59. ↵
    1. Seminowicz DA,
    2. Moayedi M
    (2017) The dorsolateral prefrontal cortex in acute and chronic pain. J Pain 18:1027–1035. pmid:28400293
    OpenUrlCrossRefPubMed
  60. ↵
    1. Seminowicz DA,
    2. Wideman TH,
    3. Naso L,
    4. Hatami-Khoroushahi Z,
    5. Fallatah S,
    6. Ware MA,
    7. Jarzem P,
    8. Bushnell MC,
    9. Shir Y,
    10. Ouellet JA,
    11. Stone LS
    (2011) Effective treatment of chronic low back pain in humans reverses abnormal brain anatomy and function. J Neurosci 31:7540–7550. doi:10.1523/JNEUROSCI.5280-10.2011
    OpenUrlAbstract/FREE Full Text
  61. ↵
    1. Seminowicz DA,
    2. Burrowes SA,
    3. Kearson A,
    4. Zhang J,
    5. Krimmel SR,
    6. Samawi L,
    7. Furman AJ,
    8. Keaser ML,
    9. Gould NF,
    10. Magyari T,
    11. White L,
    12. Goloubeva O,
    13. Goyal M,
    14. Peterlin BL,
    15. Haythornthwaite JA
    (2020) Enhanced mindfulness-based stress reduction in episodic migraine: a randomized clinical trial with magnetic resonance imaging outcomes. Pain 161:1837–1846. doi:10.1097/j.pain.0000000000001860 pmid:32701843
    OpenUrlCrossRefPubMed
  62. ↵
    1. Skorobogatykh K,
    2. van Hoogstraten WS,
    3. Degan D,
    4. Prischepa A,
    5. Savitskaya A,
    6. Ileen BM,
    7. Bentivegna E,
    8. Skiba I,
    9. D'Acunto L,
    10. Ferri L,
    11. Sacco S,
    12. Hansen JM,
    13. Amin FM
    , European Headache Federation School of Advanced Studies (EHF-SAS) (2019) Functional connectivity studies in migraine: what have we learned? J Headache Pain 20:108. pmid:31747874
    OpenUrlPubMed
  63. ↵
    1. Smith SM,
    2. Nichols TE,
    3. Vidaurre D,
    4. Winkler AM,
    5. Behrens TE,
    6. Glasser MF,
    7. Ugurbil K,
    8. Barch DM,
    9. Van Essen DC,
    10. Miller KL
    (2015) A positive-negative mode of population covariation links brain connectivity, demographics and behavior. Nat Neurosci 18:1565–1567. pmid:26414616
    OpenUrlCrossRefPubMed
  64. ↵
    1. Spitzer RL,
    2. Kroenke K,
    3. Williams JB,
    4. Löwe B
    (2006) A brief measure for assessing generalized anxiety disorder: the GAD-7. Arch Intern Med 166:1092–1097. pmid:16717171
    OpenUrlCrossRefPubMed
  65. ↵
    1. Sullivan MJ,
    2. Bishop SR,
    3. Pivik J
    (1995) The pain catastrophizing scale: development and validation. Psychol Assess 7:524–532. doi:10.1037/1040-3590.7.4.524
    OpenUrlCrossRef
  66. ↵
    1. Tabachnick BG,
    2. Fidell LS,
    3. Ullman JB
    (2007) Using multivariate statistics. Boston: Pearson.
  67. ↵
    1. Torres-Ferrús M,
    2. Quintana M,
    3. Fernandez-Morales J,
    4. Alvarez-Sabin J,
    5. Pozo-Rosich P
    (2017) When does chronic migraine strike? A clinical comparison of migraine according to the headache days suffered per month. Cephalalgia 37:104–113. doi:10.1177/0333102416636055 pmid:26961321
    OpenUrlCrossRefPubMed
  68. ↵
    1. Tustison NJ,
    2. Avants BB,
    3. Cook PA,
    4. Zheng Y,
    5. Egan A,
    6. Yushkevich PA,
    7. Gee JC
    (2010) N4ITK: improved N3 bias correction. IEEE Trans Med Imaging 29:1310–1320. doi:10.1109/TMI.2010.2046908 pmid:20378467
    OpenUrlCrossRefPubMed
  69. ↵
    1. Van Horn JD,
    2. Toga AW
    (2009) Multi-site neuroimaging trials. Curr Opin Neurol 22:370–378. doi:10.1097/WCO.0b013e32832d92de pmid:19506479
    OpenUrlCrossRefPubMed
  70. ↵
    1. Vossel S,
    2. Geng JJ,
    3. Fink GR
    (2014) Dorsal and ventral attention systems: distinct neural circuits but collaborative roles. Neuroscientist 20:150–159. doi:10.1177/1073858413494269 pmid:23835449
    OpenUrlCrossRefPubMed
  71. ↵
    1. Wang HT,
    2. Poerio G,
    3. Murphy C,
    4. Bzdok D,
    5. Jefferies E,
    6. Smallwood J
    (2018) Dimensions of experience: exploring the heterogeneity of the wandering mind. Psychol Sci 29:56–71. pmid:29131720
    OpenUrlCrossRefPubMed
  72. ↵
    1. Wang HT,
    2. Smallwood J,
    3. Mourao-Miranda J,
    4. Xia CH,
    5. Satterthwaite TD,
    6. Bassett DS,
    7. Bzdok D
    (2020) Finding the needle in a high-dimensional haystack: canonical correlation analysis for neuroscientists. Neuroimage 216:116745. doi:10.1016/j.neuroimage.2020.116745 pmid:32278095
    OpenUrlCrossRefPubMed
  73. ↵
    1. Wang Y,
    2. Cao D,
    3. Remeniuk B,
    4. Krimmel S,
    5. Seminowicz DA,
    6. Zhang M
    (2017) Altered brain structure and function associated with sensory and affective components of classic trigeminal neuralgia. Pain 158:1561–1570. pmid:28520647
    OpenUrlPubMed
  74. ↵
    1. Whitfield-Gabrieli S,
    2. Nieto-Castanon A
    (2012) Conn: a functional connectivity toolbox for correlated and anticorrelated brain networks. Brain Connect 2:125–141. doi:10.1089/brain.2012.0073 pmid:22642651
    OpenUrlCrossRefPubMed
  75. ↵
    1. Xia CH, et al
    . (2018) Linked dimensions of psychopathology and connectivity in functional brain networks. Nat Commun 9:1–14. doi:10.1038/s41467-018-05317-y
    OpenUrlCrossRefPubMed
  76. ↵
    1. Xue T,
    2. Yuan K,
    3. Zhao L,
    4. Yu D,
    5. Zhao L,
    6. Dong T,
    7. Cheng P,
    8. von Deneen KM,
    9. Qin W,
    10. Tian J
    (2012) Intrinsic brain network abnormalities in migraines without aura revealed in resting-state fMRI. PLoS One 7:e52927. pmid:23285228
    OpenUrlCrossRefPubMed
  77. ↵
    1. Yeo BT,
    2. Krienen FM,
    3. Sepulcre J,
    4. Sabuncu MR,
    5. Lashkari D,
    6. Hollinshead M,
    7. Roffman JL,
    8. Smoller JW,
    9. Zöllei L,
    10. Polimeni JR
    (2011) The organization of the human cerebral cortex estimated by intrinsic functional connectivity. J Neurophysiol 106:1125–1165.
    OpenUrlCrossRefPubMed
  78. ↵
    1. Yu M,
    2. Linn KA,
    3. Cook PA,
    4. Phillips ML,
    5. McInnis M,
    6. Fava M,
    7. Trivedi MH,
    8. Weissman MM,
    9. Shinohara RT,
    10. Sheline YI
    (2018) Statistical harmonization corrects site effects in functional connectivity measurements from multi-site fMRI data. Hum Brain Mapp 39:4213–4227. doi:10.1002/hbm.24241 pmid:29962049
    OpenUrlCrossRefPubMed
  79. ↵
    1. Yuan H,
    2. Spare NM,
    3. Silberstein SD
    (2019) Targeting CGRP for the prevention of migraine and cluster headache: a narrative review. Headache 59:20–32. doi:10.1111/head.13583
    OpenUrlCrossRefPubMed
  80. ↵
    1. Zhang Y,
    2. Brady M,
    3. Smith S
    (2001) Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm. IEEE Trans Med Imaging 20:45–57. doi:10.1109/42.906424 pmid:11293691
    OpenUrlCrossRefPubMed
Back to top

In this issue

The Journal of Neuroscience: 42 (31)
Journal of Neuroscience
Vol. 42, Issue 31
3 Aug 2022
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Masthead (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Three Dimensions of Association Link Migraine Symptoms and Functional Connectivity
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Three Dimensions of Association Link Migraine Symptoms and Functional Connectivity
Samuel R. Krimmel, Danielle D. DeSouza, Michael L. Keaser, Bharati M. Sanjanwala, Robert P. Cowan, Martin A. Lindquist, Jennifer A. Haythornthwaite, David A. Seminowicz
Journal of Neuroscience 3 August 2022, 42 (31) 6156-6166; DOI: 10.1523/JNEUROSCI.1796-21.2022

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
Three Dimensions of Association Link Migraine Symptoms and Functional Connectivity
Samuel R. Krimmel, Danielle D. DeSouza, Michael L. Keaser, Bharati M. Sanjanwala, Robert P. Cowan, Martin A. Lindquist, Jennifer A. Haythornthwaite, David A. Seminowicz
Journal of Neuroscience 3 August 2022, 42 (31) 6156-6166; DOI: 10.1523/JNEUROSCI.1796-21.2022
del.icio.us logo Digg logo Reddit logo Twitter logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • chronic pain
  • fMRI
  • individual variability
  • mindfulness-based stress reduction
  • networks
  • pain catastrophizing

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Research Articles

  • Cortically-evoked movement in humans reflects history of prior executions, not plan for upcoming movement
  • Neuronally-derived soluble Abeta evokes cell-wide astrocytic calcium dysregulation in absence of amyloid plaques in vivo
  • Effect of aging and a dual orexin receptor antagonist on sleep architecture and NREM oscillations including a REM Behavior Disorder phenotype in the PS19 mouse model of tauopathy
Show more Research Articles

Neurobiology of Disease

  • Neuronally-derived soluble Abeta evokes cell-wide astrocytic calcium dysregulation in absence of amyloid plaques in vivo
  • Effect of aging and a dual orexin receptor antagonist on sleep architecture and NREM oscillations including a REM Behavior Disorder phenotype in the PS19 mouse model of tauopathy
  • NMDA Receptors at Primary Afferent–Excitatory Neuron Synapses Differentially Sustain Chemotherapy- and Nerve Trauma-Induced Chronic Pain
Show more Neurobiology of Disease
  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
(JNeurosci logo)
(SfN logo)

Copyright © 2023 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.