Predicting response to tVNS in patients with migraine using functional MRI: A voxels-based machine learning analysis

Background Migraine is a common disorder, affecting many patients. However, for one thing, lacking objective biomarkers, misdiagnosis, and missed diagnosis happen occasionally. For another, though transcutaneous vagus nerve stimulation (tVNS) could alleviate migraine symptoms, the individual difference of tVNS efficacy in migraineurs hamper the clinical application of tVNS. Therefore, it is necessary to identify biomarkers to discriminate migraineurs as well as select patients suitable for tVNS treatment. Methods A total of 70 patients diagnosed with migraine without aura (MWoA) and 70 matched healthy controls were recruited to complete fMRI scanning. In study 1, the fractional amplitude of low-frequency fluctuation (fALFF) of each voxel was calculated, and the differences between healthy controls and MWoA were compared. Meaningful voxels were extracted as features for discriminating model construction by a support vector machine. The performance of the discriminating model was assessed by accuracy, sensitivity, and specificity. In addition, a mask of these significant brain regions was generated for further analysis. Then, in study 2, 33 of the 70 patients with MWoA in study 1 receiving real tVNS were included to construct the predicting model in the generated mask. Discriminative features of the discriminating model in study 1 were used to predict the reduction of attack frequency after a 4-week tVNS treatment by support vector regression. A correlation coefficient between predicted value and actual value of the reduction of migraine attack frequency was conducted in 33 patients to assess the performance of predicting model after tVNS treatment. We vislized the distribution of the predictive voxels as well as investigated the association between fALFF change (post-per treatment) of predict weight brain regions and clinical outcomes (frequency of migraine attack) in the real group. Results A biomarker containing 3,650 features was identified with an accuracy of 79.3%, sensitivity of 78.6%, and specificity of 80.0% (p < 0.002). The discriminative features were found in the trigeminal cervical complex/rostral ventromedial medulla (TCC/RVM), thalamus, medial prefrontal cortex (mPFC), and temporal gyrus. Then, 70 of 3,650 discriminative features were identified to predict the reduction of attack frequency after tVNS treatment with a correlation coefficient of 0.36 (p = 0.03). The 70 predictive features were involved in TCC/RVM, mPFC, temporal gyrus, middle cingulate cortex (MCC), and insula. The reduction of migraine attack frequency had a positive correlation with right TCC/RVM (r = 0.433, p = 0.021), left MCC (r = 0.451, p = 0.016), and bilateral mPFC (r = 0.416, p = 0.028), and negative with left insula (r = −0.473, p = 0.011) and right superior temporal gyrus/middle temporal gyrus (r = −0.684, p < 0.001), respectively. Conclusions By machine learning, the study proposed two potential biomarkers that could discriminate patients with MWoA and predict the efficacy of tVNS in reducing migraine attack frequency. The pivotal features were mainly located in the TCC/RVM, thalamus, mPFC, and temporal gyrus.