EEG μ -RHYTHM REACTIVITY IN CHILDREN DURING IMITATION OF BIOLOGICAL AND NON-BIOLOGICAL MOTION

The development of brain-computer interfaces based on the use of EEG sensorimotor rhythms reactivity parameters and designed for the rehabilitation of people (including children) with impaired motor functions is currently relevant. The study was aimed to analyse the EEG µ-rhythm in the individual frequency range in children during imitation of biological and non-biological motion. EEG was recorded at frontal, central and parietal cortical regions in 136 normally developing right-handed children aged 4–15, at rest and during the execution and imitation of movements using the computer mouse. When the children moved the computer mouse on their own (F 1, 132 = 31.17; p < 0.001) and executed the concentric moving of the coloured circle (F 1, 132 = 90.34; p < 0.001), the µ-rhythm desynchronization developed in the frontal, central and parietal neocortical regions. The µ-rhythm synchronization was detected during the non-biologocal motion imitation (F 1, 132 = 12.65; p < 0.001), compared to the task on the autonomous movement execution. The µ-rhythm desynchronization was observed during the biologocal motion imitation in relation to autonomous movement execution (F 1, 132 = 9.58; p = 0.002). The described effects had their own features in the groups of children aged 4–6, 7–9, 10–12 and 13–15. The study results demonstrate the desirability of taking into account the µ-rhythm reactivity age-related features and the visual stimuli nature when developing software for the brain-computer interfaces.

. Task queue. А. Concentric colored circle moving on the monitor screen using the computer mouse (Mn.1). B. Imitation of other colored circle's movement (ImNB). C. Children move the computer mouse in circle on their own (Mn.2). D. Children imitate the researcher's movements (ImB) А B C D is important for recognition of movements and associated intentions. Mirror neurons are neurons able to activate in a similar way, both when executing actions, and when watching other individuals executing similar actions [1,2]. It has been suggested that MNS plays an important role in the complex forms of social interaction. [3]. Desynchronization of the EEG sensorimotor rhythm, the µ-rhythm, is considered to be a marker of MNS activation [4]. Since modulations of the α-rhythm in the occipital region can overlap the effects of µ-rhythm desynchronization [5], to determine the individual frequency range and reactivity, the following features are taken into account: unlike the occipital α-rhythm, the µ-rhythm is most pronounced in the frontoparietal regions; the µ-rhythm amplitude decreases when the subject moves, imagines movement or watches the other subjects' movement, but does not change significantly when the subject opens or closes the eyes [6,7].
Sensorimotor rhythm amplitude depression during the biological motion watching is more pronounced than during watching the non-biological motion [8], which is also characteristic of watching social actions, compared to actions outside the social context [9]. Regarding the movement imitation, it is assumed that imitation is associated with the activation of human MNS, and is the result of the comparison of the observed action and the internal motor plan for the execution of action [10].
The study of MNS and µ-rhythm reactivity is of special interest in a view of new methods development for rehabilitation of patients with various motor impairments using braincomputer interfaces [11,12]. In particular, in the treatment of adult patients, the synchronous interfaces are used, based on the analysis of the EEG sensorimotor rhythms reactivity when representing the movement in response to the signal presented [13,14]. Recently, such methods are beginning to be used for rehabilitation of children with cerebral palsy [15]. Symbols or text commands are reported to be used as the signals presented to patients. However, the concept of MNS suggests that stimuli visually representing movements and requiring the simulation of movements could be more effective for triggering reactions in the EEG µ-rhythm range. It should be noted that when working with children it is preferable to use the actions that are in the child's motor repertoire for more effective task execution [16]. The study was aimed to analyse the EEG µ-rhythm under conditions of biological and non-biological motion imitation in children aged 4-15 using the computer pointer device, the mouse. Now, even the preschool children are familiar with the computer mouse operation.

Characteristics of a sample
The study was performed at the Center for Collective Use of Scientific Equipment "Experimental physiology and biophysics" of V.I. Vernadsky Crimean Federal University.
The study included 136 righ-handed children aged 4-15 (69 boys and 67 girls). Inclusion criteria: normal (or corrected to normal) vision and hearing; preferred right hand when operating the computer mouse; sufficient degree of cognitive development (IQ at least 80 points according to the Wechsler scale, variants WPPSI and WISC). Exclusion criteria: taking the CNS affecting drugs; severe chronic somatic diseases. The children were divided into four age groups: 4-6 years (30 people), 7-9 years (46 people), 10-12 years (30 people) and 13-15 years (30 people).

EEG recording
EEG recording was performed using the Neuron-Spectrum-3 EEG System (Neurosoft; Russia). Data were obtained using the WinEEG version 2.8 software (available for free). Independent component analysis was used for the artifacts correction. The 19 monopolar EEG electrodes were used in accordance with the 10-20 system. In our study, the frontal, central and parietal neocortical regions were the area of concern (F3, F4, Fz, С3, С4, Cz, P3, P4, Pz loci). Paired electrodes attached to the the ear lobes were the reference electrode. Cut-off frequencies of the high and low pass filters were 1.5 and 35 Hz, respectively, EEG digitization rate was 250 Hz.
EEG recording was performed while the children performed a queue of sequential tasks, the duration of each task was 30 s. EEG segments were processed using the Fast Fourier Transform with the 4 epoch of analysis and 50% mutual overlapping of epochs.
To imitate the non-biological motion, the following tasks were used: 1) gaze fixation on the video of the computer mouse (baseline); 2) concentric moving of the coloured circle on the monitor screen using the computer mouse (Mn.1) (Fig. 1А); 3) imitation of the other coloured circle's motion (imitation of non-biological motion, ImNB) (Fig. 1B).
When imitating the biological motion, the subject and the researcher were located at the tables next to each other (the researcher on the right), each of tables had a monitor and a computer mouse on it. Using the webcam, the working plane of the researcher's table with the mouse on it was demonstrated on the monitor in front of the subject. The tasks queue was as follows: 1) gaze fixation on the video of the computer mouse (baseline); Fig. 2. EEG µ-rhythm amplitude (А, Ln µV) in children aged 4-6 during imitation of non-biological (А) and biological (B) motion. 1 -baseline, 2 -autonomous movements' execution at arbitrary speed, 3 -motion imitation. Amplitude differences between baseline and autonomous movements' execution: * -p < 0.05; **р ≤ 0.01; *** -р ≤ 0.001; when executing autonomous movements and imitating: ▲ -p < 0.05

А B
2) moving the computer mouse in a circle by children on their own (Mn.2) (Fig. 1C); 3) imitation of the researcher's movements by the children (imitation of biological motion, ImB) (Fig. 1D).
EEG was analysed in the individual µ-rhythm frequency range defined when the subject moved his right hand on his own (С3). The full frequency range of the µ-rhythm (6-13 Hz) was divided into segments of 1 Hz. As an individual frequency range, two adjacent segments were taken with maximum desynchronization in relation to baseline [17]. The µ-rhythm amplitude within the individual frequency range was calculated for each experimental situation. Log transformation was used for normalization of the amplitude values distribution.
Reactivity indices were used for comparison of µ-rhythm parameters under conditions of biological and non-biological movement imitation. These indices were calculated according to the generally accepted scheme [18] using the following formula: k = ln (B/A), where k is the sensorimotor rhythm reactivity index, B is the sensorimotor rhythm amplitude in the major situation, and A is the sensorimotor rhythm amplitude in the initial reference situation (baseline or subjects' moving on their own). Positive reactivity index values corresponded to synchronization of the sensorimotor rhythm, and negative values corresponded to desynchronization.

Statistical analysis
Statistical analysis was performed using the STATISTICA 12.0 software (StatSoft Inc.; USA). To describe the non-normal distributions, median and interquartile range were used, the differences between the groups were evaluated using the Mann-Whitney U-test. For normal data distribution, the mean and standard error of the mean were used. The differences of the amplitude and reactivity indices of the µ-rhythm recorded in different experimental situations were evaluated by the repeated measures ANOVA. The 4×2×9 scheme was used for assessment of the one intersubjective factor (age group, AGE) and two intrasubjective factors (situation, SIT, and locus, LOC) influence. To calculate the statistical significance of the sensorimotor rhythm differences in relation to each of the nine EEG derivations within each age group, the ad-hoc analysis method (F-distribution estimation) was used.

μ-rhythm frequency parameters
The median values of the individual µ-rhythm range lower boundary were 9 Hz (8.5; 10), the extreme values were 6 and 11 Hz. The median values of the individual µ-rhythm range upper boundary were 11 Hz (10.5; 12), and the extreme values were 8 and 13 Hz. The differences between the age groups were not significant.

EEG μ-rhythm amplitude at rest and under condition of motion execution and imitation
The µ-rhythm amplitude differences analysis of variance in the Mn.1 situation in relation to baseline taking into account the age group and EEG locus revealed the sifnificant influence of the SIT (F The µ-rhythm amplitude differences analysis of variance in the Mn.2 situation in relation to baseline taking into account the age group and EEG locus revealed the sifnificant influence of the SIT (F 1, 132 = 31.17; p < 0.001), AGE (F 3, 132 = 6.46; p < 0.001) and LOC (F 8, 1056 = 71.55; p < 0.001), factors, as well as the SIT×LOC (F 8, 1056 = 28.32; p < 0.001) and SIT×AGE (F 3, 132 = 6.35; p < 0.001) interactions. Evaluation of the µ-rhythm amplitude changes in the ImB situation in relation to Mn.2 revealed the significant influence of the SIT (F 1, 132 = 9.58; p = 0.002), AGE (F 3, 132 = 18.63; p < 0.001) and LOC (F 8, 1056 = 54.08; p < 0.001) factors, as well as the SIT×LOC (F 8, 1056 = 3.28; p = 0.001) and SIT×AGE (F 3, 132 = 6.2; p = 0.001) interactions.  In children aged 4-6 executing the concentric moving of the coloured circle, the significant EEG µ-rhythm desynchronization was detected (Mn.1 in relation to baseline) in most studied regions. The µ-rhythm amplitude changes in children imitating the coloured circle movement (ImNB in relation to Mn.1) were not significant (Fig. 2А). When the children of that age moved the computer mouse on their own, the significant EEG µ-rhythm amplitude increase (Mn.2 in relation to baseline) was detected in the the right hemisphere central locus (C4). When the children imitated the researcher's movements, the sensorimotor rhythm desynchronization (ImB in relation to Mn.2) was registered in the в mid-parietal locus (Pz) (Fig. 2B).
In the group of children aged 7-9, the significant depression of µ-rhythm in the Mn.1 situation was observed in most studied regions. In the ImNB situation (in relation to Mn.1) the µ-rhythm amplitude changes were not significant (Fig. 3А).
In the Mn.2 situation the significant sensorimotor rhythm desynchronization was detected in the central (C3 and Cz) and all parietal loci. In the ImB situation (in relation to Mn.2) there were no significant µ-rhythm amplitude changes (Fig. 3B).
In children aged 10-12, in the Mn.1 sitiation the significant µ-rhythm suppression was detected in most studied loci. During the coloured circle movement imitation the significant sensorimotor rhythm synchronization (ImNB in relation to Mn.1) was detected in the mid-frontal locus (Fig. 4А). In the Mn.2 situation the significant decrease in µ-rhythm amplitude was observed in the central (C3 и Cz) and all parietal loci. In the ImB situation (in relation to Mn.2) there were no significant sensorimotor rhythm amplitude changes (Fig. 4B).
In the group of teenagers aged 13-15, in the Mn.1 sitiation the significant µ-rhythm suppression was observed in most studied loci. During the non-biological motion imitation (ImNB in relation to Mn.1) the significant sensorimotor rhythm  synchronization was detected in most studied regions (Fig. 5А). When the subjects moved the computer mouse on their own, significant sensorimotor rhythm desynchronization (Mn.2 in relation to baseline) was registered in all studied regions. In the ImB situation, the additional (compared to previous task) µ-rhythm desynchronization was observed that was significant in all loci (Fig. 5B).

EEG μ-rhythm reactivity comparison under conditions of biological and non-biological motion imitation
To evaluate the differences of µ-rhythm reactivity in the ImNB and ImB situations (compared to execution of movements by children on their own at arbitrary speed), the reactivity indices analysis of variance was performed taking into account the age group and EEG locus. The mean µ-rhythm reactivity index values for children of four age groups are presented in Tables  In the group of pre-school children, the significant µ-rhythm reactivity indices differences in the ImNB and ImB situations were detected in the parietal loci Pz and P4 (p = 0.03). In children aged 7-9, the significant differences were observed in the locus Fz (p = 0.04). In children aged 10-12, no significant µ-rhythm reactivity indices differences were detected. In the group of teenagers aged 13-15, the differences were significant in all studied regions (p ≤ 0.001).

DISCUSSION
According to the study results, the individual sensorimotor rhythm frequency ranges of the 4-15 years old children vary widely, and there are no significant differences in the mean values between different age groups. In the other authors's paper [19] reporting the EEG µ-rhythm reactivity analysis in the selected frequency range in children aged 4-11, the average sencorimotor rhythm band was 9-11 Hz. High sensorimotor rhythm parameters variability among the individuals and no association with the children's age were detected. These indicate the need to determine the children's individual frequency range when studying the sensorimotor rhythm reactivity, as well as when attempting the correction using the µ-rhythm parameters (EEG based neurofeedback training, correction using the brain-computer interface). Analysis of the µ-rhythm amplitude changes demonstrated that in 4-6 years old children arbitrarily moving the coloured circle (Mn.1) the significant EEG µ-rhythm desynchronization in the frontal and central loci of left hemisphere, as well as in the median frontal and all parietal loci (F3, Fz, C3, P3, Pz, P4) could be detected. The results of our study are consistent with the literature data on the sensorimotor rhythm desynchronization during the voluntary movements' execution [20]. When the children moved the computer mouse on their own (Mn.2), no significant µ-rhythm amplitude decrease was observed. It is possible that for children of this age, the task of relatively simple circular movements' execution with a computer mouse was simpler than the task of capturing and moving a colored circle using the computer mouse left button, and it did not require any special motor control. During the task execution, the significant EEG µ-rhythm amplitude increase in the locus C4 was registered, which could be due to inhibition of the ipsilateral hemisphere (in relation to the hand used) [21].
In the groups of children aged 7-9 and 10-12 executing the movements on their own, the µ-rhythm desynchronization was detected in most studied regions. The concentric coloured circle moving (Mn.1) unlike the mouse moving in a circle (Mn.2) was also associated with the sensorimotor rhythm desynchronization in the frontal loci (F3, Fz). It is known that the frontal cortical regions are responsible for planning and preparation of complex movements [22]. It is also assumed that more complex motor actions are accompanied by a more widespread µ-activity desynchronization [23]. Presumably, moving the color circle in the group of 7-12 years old children, as well as in the group of younger children, required considerable effort, which led to the involvement of the cerebral cortex frontal region.
The situations of biological and non-biological motion imitation in children aged 4-6, 7-9 and 10-12 were associated with almost no additional modulation of the µ-rhythm in relation to the arbitrary movements' execution. This may indicate that in children of said age the required for processing In the group of teenagers aged 13-15, the significant sensorimotor rhythm desynchronization during the autonomous movements' execution was detected in all studied regions. The µ-rhythm amplitude decrease above the frontal, central and parietal loci in elder children may be due to development of connections between the neocortical regions involved. During the non-biological motion imitation a smaller drop in the sensorimotor rhythm amplitude was observed than during the autonomous movements execution and biological motion imitation (which is especially pronounced in the frontal and central loci). It can be assumed, that the need to imitate the movements of another object (colored circle) led to the shift of attention to its perception and, as a result, to weakening of one's own movements' motor control. In children of this age, in the ImB situation, the additional (compared to the observed during the Mn.2 task execution) significant µ-rhythm desynchronization in all loci was detected. The more pronouced reaction in the parietal loci is noteworthy. It is known, that parietal cortical regions are involved into the information processing during watching the human's motion (compared to watching the non-biological objects' motion) [24]. The sensorimotor rhythm modulation, revealed by us in the described regions during the biological motion imitation, may be due to involvement of the parietal cortex MNS components responsible for coding of goals underlying the watched movements [25]. The mirror neurons are associated with the cognitive integration of visual, auditory and motor stimuli needed for social interaction in children [26] and adults [27]. Thus, it can be assumed, that the additional µ-rhythm desynchronization during the other man's movements imitation is caused precisely by the social context to which the MNS is sensitive.
Comparison of the µ-rhythm reactivity indices for imitation tasks revealed that the biological motion imitation in elder children was associated with the greater desynchronization, compared with the situation of color circle movements' imitation. As already noted, similar features of the sensorimotor rhythm reactivity during watching the biological and non-biological objects movements were detected in adult volunteers [8].
A sensorimotor rhythm reactivity patterns comparative analysis in children of different ages allows us to come to a number of conclusions. In the group of youngest children (4-6 years), the most pronounced activation of the frontal, central and parietal cortical regions, manifested in the µ-rhythm amplitude decrease (more pronounced in the left hemisphere), is observed during computer mouse operation associated with a non-biological motion (coloured circle) (Mn.1). Random rhythm computer mouse movements' execution (Mn.2) does not lead to the significant decrease in the µ-rhythm amplitude, and the biological motion imitation (the other person's hand movement) is not associated with any additional activation in most loci. Thus, in the described experimental situation, in pre-school children, the cortical center of motor analyzer is especially sensitive to manipulations with biological objects. In elder children (7-9 and 10-12 year), a similar neocortical activation pattern was revealed during execution of movements associated with non-biological objects, computer mouse moving and biological object (researcher's hand) motion imitation. Unlike the previously described groups, the children aged 13-15 demonstrate the significant µ-rhythm desynchronization in the frontal, central and parietal cortical regions of both hemispheres during imitation of the other person's motion.
It stands to reason, that processes of perception and other person's movement imitation in younger children are in their infance, and in teenagers, these processes are rather developed and similar to those in adults. In teenagers, the pronounced µ-rhythm desynchronization in all studied regions during moving on their own or imitating the biological visual images motion may be due to maturation of motor, sensorimotor and associative cortical regions involved in the execution and imitation of movements [28]. The revealed agerelated sensorimotor rhythm reactivity features may be used for improvement of existing rehabilitation techniques based on the EEG-controlled robotic systems for children with cerebral palsy [15].

CONCLUSION
When children aged 4-15 move the computer mouse on their own, the µ-rhythm desynchronization develops in the frontal, central and parietal neocortical regions, which is more pronounced in the left hemisphere. When the children aged 4-6, 7-9 and 10-12 imitate the biological and non-biological motion no significant additional µ-rhythm modulation is revealed, compared to the execution of movements on their own. In children aged 13-15, the highest sensorimotor rhythm desynchronization is observed during the resercher's hand