The power features of Masseter muscle activity in tension-type and migraine without aura headache during open-close clench cycles

Introduction

The generation of semi-automatic movements like chewing is in the central pattern generator (CPG). CPG is a group of neurons that are organized in the spinal cord to induce repeated movement. CPG is modulated by the supra-spinal descending pathway (Vogt, Pfeifer & Banzer, 2003). The proprioception afferents and pain projections also have an important role to modulate the CPG (Westberg, 1999). The projection of CPG and peripheral sensory afferents such as proprioception and pain are integrated to gamma motor neuron to control sensitivity of motor neuron (Capra, Hisley & Masri, 2007; Ellaway, Taylor & Durbaba, 2015). The change of motor unit activity is reflected and it is visible in some features of EMG (Farina & Holobar, 2016).

Some studies showed the cervico-cephalic reflexes are sensitized in headache due to central pain facilitation (Reshkova, Bogdanova & Milanov, 2015). Electromyography recording of muscles in headache patients especially in tension type headache, (TTH) showed increase of intensity of rest and duration of masseter muscle activity during chewing-like motion (Hatef et al., 2012; Jensen, Fuglsang-Frederiksen & Olesen, 1994; Shimada, Baad-Hansen & Svensson, 2015; Sohn, Choi & Jun, 2013). Some models had been proposed based on the results of many studies about interaction of pain and motor activity (Sessle et al., 2008). Modified pain adaptation model explains that the central motor command to control pattern of agonist and antagonist muscle activity during a rhythmic movement changes to limit the range of mobility and variation in movement (Minami et al., 2013; Peck, Murray & Gerzina, 2008).

It seems that several types of headaches influence motor planning differently. This means that whereas the regions of pain of TTH and migraine are very close to each other, the changes of intensity and duration of masseter muscle activity and responses to motor reflex were significantly different between them (Hatef et al., 2012; Jensen, Fuglsang-Frederiksen & Olesen, 1994; Proietti Cecchini et al., 2003). Previous studies had assessed the static or dynamic activity of muscles in headache and reported the mean of changes not the trend during muscular activity. This means that they did not evaluate the pattern of activity during a period of time. To evaluate the effect of pain and proprioception disturbance on the output of CPG in the EMG, the aim of this study was to investigate the effect of TTH and migraine without aura (MOA) headache on the pattern of the energy distribution of the masseter muscle activity during repeated open-close-clench cyclesinteracted with TMJ click.

Methods

Noise removing

Noise removing from EMG signals must be done by considering the signal properties. EMG is a signal in which various noises and multiple sources of noises appear (Kale & Dudul, 2009; Katariya & Devasahayam, 2006). Noise removing has been done by two steps as follows:

Removing unwanted frequencies

EMG has a frequency range between 10 and 500 Hz, so a high-pass filter was used to remove the frequencies below the threshold by applying a Butterworth filter which creates the pass band as flat as possible (Thede, 2005). Based on Nyquist theory, using the sampling rate of 1024 Hz in data collection, there is no need to use a low-pass filter.

Suppressing sparse areas of EMG signal

EMG of masseter in sequential jaw pressing has a semi-regular pattern of sequential activity and silence (Fig. 1). The silence area (small ellipses) where noises exist without a real activity must be removed. One of the methods that can be used in this situation is Overlapping Group Sparse Denoising (OGSD) that is based on convex optimization and non- convex regularization. The method acts using the combination of convex and non-convex optimization and has the benefits of both approaches simultaneously (Chen & Selesnick, 2014).

In order to use OGSD for noise removing, the signal must be in the form of Group—Sparse in which the values of the signals tend to create separate clusters, without predefined clear boundaries and also their start and finish points are not regular (Chen & Selesnick, 2014). As shown in Fig. 1, EMG signal of masseter muscle in jaw pressing has such properties. In this study, OGSD was used to remove sparse area noises.

Statistical analysis

To compare the baseline characteristics such as the age between three groups, year and the frequency of headache between TTH and MOA groups, the Mann–Whitney test was used. Chi-square test was done to compare the distribution of the class of occlusion and existence of click in TMJ. The feature extraction was done for all the signals, and then to compare the three groups, six parameters of energy were computed on (Table 1). Mean, Variance, Skewness, Kurtosis, and first and second half energy over all signal energy were calculated for the three groups of signals on both left and right sides.

Table 1:

Six parameters of energy.

Parameter	Description	Interpretation
Mean (1st moment)	Estimating the value around which central clustering occurs	High value indicates high energy in signal overall.
Variance (2nd moment)	The measurement of the spread between the numbers in a data set. The variance measures how far each number in the set is from the mean.	High value indicates high difference of energy level.
Skewness (3rd moment)	Characterizing the degree of asymmetry of a distribution around its mean	Positive skewness indicates that the data are positively skewed or skewed right, meaning that the right tail of the distribution is longer than the left (decreasing energy by time) and vice versa.
Kurtosis (4th moment)	Measuring the relative peakedness or flatness of a distribution	High kurtosis value shows the concentration of activity in the middle of the signal to the tails.
First half energy over all signal energy	The sum of signal activity in the first half of signal over the whole signal	The high value of the parameters shows a high portion of signal energy occurred in the first part of the signal.
Second half energy over all signal energy	The sum of signal activity in the second half of signal over the whole signal	The high value of the parameters shows a high portion of signal energy occurred in the second part of the signal.

DOI: 10.7717/peerj.3556/table-1

It was useful to characterize the set by a few numbers related to its moments below general relation for n’th moments (Press et al., 2015). ${\displaystyle {\mathrm{\mu }}_n=\mathrm{E}\left[{\left(X-\mathrm{E}\left[X\right]\right)}^n\right].}$

Clicked TMJ and side were considered as factors in addition to the three groups. Then, the cases were divided into two categories of clicked or no-clicked TMJ. Another categorization was based on the masseter of left or right side.

Subsequently, a three-way ANOVA analysis (Control/Migraine/Tension Headache, left/right, and TMJ with/without click) was conducted upon the above parameters since they had a normal distribution. The flow chart of all the steps is shown in Fig. 2. A P-value less than 0.05 was considered as significant.

Results

Table 2 gives a comparative statistics of the demographic and clinical characteristics for the participants of the three groups. The prevalence of occlusion type classes and TMJ with click (click or crepitus) category were not significantly different among the three groups. The patients of TTH group were significantly younger than the two other groups. The mean of the frequency of headache was around nine headaches per month in both headache groups but the years of headache in MOA group were twice as more as TTH group.

Tables 3 and 4 represent the descriptive statistics of EMG variables (mean, lower and upper band of 95% confidence interval of the mean) of the three groups (Control, MOA and TTH groups), two sides (right and left) with or without clicked TMJ.

The data were distributed normally and the equal variances were assumed. The results of the three-way ANOVA and the following Bonferroni collection test evaluating the interaction effect of groups, side of EMG recording and the existence of clicked TMJ as main factors were presented in Table 5. When the time of cycles was divided into two parts, the energy of muscle activity can be compared between the main factors (group, TMJ click and side) in the first and last part. The power of masseter activity of TTH group had less activity in the first part to all the cycles and more power of activity in the second part to all than the two other groups independently to side and TMJ click (P-value <  0.04). In the left side, the masseter muscle activity of MOA groupshowed a higher mean power of activity than the two other groups in the clicked TMJ (P-value <  0.002) and less mean power than that in the no-clicked TMJ (P-value <  0.02) (Fig. 3). Variation from mean power significantly decreased in the MOA group than the control group in both sides in the no-clicked TMJ. Skewness which showed the pattern of power during cycles was affected by main factors. The results showed that the Skewness of the left masseter muscle in MOA group was less than the two other groups in the clicked TMJ (P-value <  0.002) and it was more only in the control group in the no-clicked TMJ (P-value <  0.009) (Fig. 4). Kurtosis of masseter activity did not show a significant difference between the groups and sides and the existence of TMJ click (P-value >  0.2).

Discussion and Conclusion

The results showed that the effect of the type of headache on the distribution of the energy of EMG activity was different in masseter muscle during the short time of chewing-like movement. The comparison of the MOA group to other groups showed that the greater the mean of the power with delayed height of energy of the left masseter muscle during cycles in the clicked TMJ and inverted pattern in the no-clicked TMJ. On one side, the variation of energy was significantly less in MOA group than the two other groups in the no-clicked TMJ. On the other side, the relative first or second time parts to all time showed that the TTH group had less energy in the first part and more energy in the second part in comparison with the other groups. The TMJ click had a critical effect on the change of muscle activity in the MOA group. Previous studies revealed that TMJ click changes the masseter muscle timing during chewing-like movement significantly (Hatef et al., 2007). The finding showed that the differences were more on the left side. Some other potential conclusions can be presented. It might be related to right-hand participants or the MOA may also have a lateralization effect. The left side might be affected by MOA greater than the right side. It was not clearwhy the TMJ click inverted the distribution of power of masseter muscle activity in the MOA subjects during chewing-like movement. However, the difference between TTH and other groups was not related to the side of recording and the existence of TMJ click.

The experimental studies showed that the pain afferents and CPG inputs are projected and integrated to interneurons. A premotor neuron of caudaltrigeminal nucleus modulates the gama motor neuron and alpha motor neuron activity (Capra, Hisley & Masri, 2007; Westberg, 1999). Nevertheless, the pain descending pathway from anterior cingulate cortex (ACC)- Periaqueductal gray matter (PAG)- Rostral ventromedial medulla (RVM)- dorsal horn (DH) axis and the projection of dorso-caudal of ACC to motor regions reveal extensive evidence to show the significant role of pain to modulate motor activity (Rainville, 2002). Previous studies demonstrated that headaches affect the timing and intensity of masseter muscle activity (Burnett et al., 2000; Hatef et al., 2012; Jensen, Fuglsang-Frederiksen & Olesen, 1994). Some of them compared the effect of two tension-type and migraine headache on EMG. They found that the level of EMG of temporalis muscle at rest or frontalis muscle at mental work in Tension-type were higher than migraine headache (Sandrini et al., 1994). The timing pattern of the masseter muscle activity during chewing showed more changes in TTH than MOA patients (Hatef et al., 2012). The evaluation of reflexes such as terigemocervical or blink reflex showed that chronic TTH and migraine had more influence on them than acute TTH (Nardone & Tezzon, 2003a; Nardone & Tezzon, 2003b; Proietti Cecchini et al., 2003). When the headache became chronic, the pain matrix included supra-spinal regions especially anterior cingulate cortex was sensitized to pain stimulus and showed more effect in motor planning (Millan, 2002).

In the present study, the years of headache in the MOA group were twice as much as the TTH group. Also, some changes in MOA were seen which did not occur in the TTH group. However, muscle pain afferents, type III fibers, are involved more in TTH than in migraines (De Tommaso & Fernandez-de Las-Penas, 2016) and type III fiber afferents in muscle had a more significant effect on the sensitivity of the motor neuron and its muscular reflex (Schomburg et al., 2012). Some of the findings are defensible because the involvement of type III fibers is more in TTH than in migraines. TTH showed that the level of the energy of masseter activity was lower in earlier cycles and it became higher in later cycles in comparison to the two other groups. The findings revealed the new and different aspects of the modulation of TTH and migraine headaches on the activity of masseter muscle during repeated jaw motion like chewing. The change of energy rising in MOA was more complex. There was no significant difference in the dominance of headache side in the MOA patients by asking about headache occurrence lately, although there was a significant difference in the pattern of energy changing during the cycles in left side with or without TMJ click. Furthermore, there is an open question about the lateralization effect of migraine on the masseter muscle behavior which did not occur in TTH.

The sensory inputs such as proprioception from muscles and most closely joints implicated in controlling the movement interacted with CPG (Prochazka & Ellaway, 2012). The disturbance of the joint kinematics could change the proprioception inputs and CPG pattern to activate the muscles (Martin et al., 2012). Click in TMJ, as a sign of jaw dysfunction, can affect jaw CPG and masseter muscle activity. The current study confirmed this sequence which showed the existence of TMJ click interacted with the effect of headache on the mean and skewness of masseter muscle power during the chewing-like movement.

Supplemental Information