Frequency-specific genetic influence on inferior parietal lobule activation commonly observed during action observation and execution

Brain activity relating to recognition of action varies among subjects. These differences have been hypothesised to originate from genetic and environmental factors although the extent of their effect remains unclear. Effects of these factors on brain activity during action recognition were evaluated by comparing magnetoencephalography (MEG) signals in twins. MEG signals of 20 pairs of elderly monozygotic twins and 11 pairs of elderly dizygotic twins were recorded while they observed finger movements and copied them. Beamformer and group statistical analyses were performed to evaluate spatiotemporal differences in cortical activities. Significant event-related desynchronisation (ERD) of the β band (13–25 Hz) at the left inferior parietal lobule (IPL) was observed for both action observation and execution. Moreover, β-band ERD at the left IPL during action observation was significantly better correlated among monozygotic twins compared to unrelated pairs (Z-test, p = 0.027). β-band ERD heritability at the left IPL was 67% in an ACE model. These results demonstrate that β-band ERD at the IPL, which is commonly observed during action recognition and execution, is affected by genetic rather than environmental factors. The effect of genetic factors on the cortical activity of action recognition may depend on anatomical location and frequency characteristics.

genes. In particular, elderly twins are good subjects to investigate the genetic and environmental effects on brain activity, because long-term differences in living environments experienced after twins aged and started living separately may have caused differences in brain function. If monozygotic twins are more similar in a trait than dizygotic twins, we can speculate that the trait is influenced by genetic rather than environmental factors.
Previous twin studies have reported similarities in brain activity between monozygotic and dizygotic twins using electroencephalography (EEG) 11,12 , magnetoencephalography (MEG) 13,14 and functional magnetic resonance imaging (fMRI) 15,16 . A comparison of EEGs in twins demonstrated that much of the variance in EEG power during the resting state is explained by genetic factors 17 . Similarly, MEG studies involving twins demonstrated that genetic factors have a large effect on MEG responses during somatosensory 18 and motor tasks 19 .
Using a source reconstruction technique for MEG signals, task-related cortical activities can be evaluated for each anatomical area. In particular, we have shown that changes in the signal power of neural oscillations, i.e. event-related synchronisation (ERS) and event-related desynchronisation (ERD) 20-22 , represent the cortical activity relevant to the task. Moreover, we have recently shown that genetic factors largely explain the alteration of frequency-specific cortical activity during word recognition 23 . Both for action recognition and action execution, β-band ERD was mainly enhanced in the sensorimotor and parietal areas [24][25][26] . Therefore, we hypothesised that the variations of β-band ERD during action observation, which are commonly observed during action execution, are due to genetic factors. In the present study, we quantified the effect of genetic and environmental factors on cortical activity related to action observation by comparing MEG responses of elderly monozygotic and dizygotic twins during the observation of hand movements.

Results
Action observation and task execution. Twenty monozygotic twins and 11 dizygotic twins were enrolled in this study after obtaining written informed consent in accordance with the protocols approved by the ethics committee of the Osaka University Graduate School of Medicine. MEG signals were recorded for all subjects lying in the recumbent position on a bed in a magnetically shielded room. During the scan, the subjects were instructed to fix their gaze on the centre of the screen in front of them. After being presented with an image of a resting left hand for 1,200 ms, a video clip of the movement of the index or middle finger being lifted was presented for 300 ms. Thereafter, a black dot was presented on the image between the index and middle fingers as a signal for the subject to reproduce the movement. When this dot appeared, the subjects had to lift their own finger immediately to imitate the movement presented in the preceding video clip (Fig. 1).

Neural oscillations during action observation and execution. Significant ERS/ERD was observed
during the observation and execution of hand movements. Figure 2 shows the spatiotemporal distribution of significant ERS/ERD based on the statistical analysis during action observation and execution. The peaks of the ERS/ ERD are distributed widely, including the middle temporal gyrus (MTG), middle occipital gyrus (MOG), inferior temporal gyrus (IFG), inferior parietal lobule (IPL) and superior parietal lobule (SPL) (Fig. 2a, Tables 1,2). Moreover, the peak varied among frequency bands and with time after presentation of the stimulus. The red points in Fig. 2a show the spatiotemporal peaks of the significant ERS/ERD for each frequency band (Tables 1 and  2). Similarly, significant ERS/ERD was observed during action execution. Interestingly, the spatiotemporal peak of β-band ERD was also observed at the left IPL during the execution of an action. Notably, at the spatiotemporal Figure 1. Schematic illustration of the task and analysis. After presenting a static picture of the hand for 1,200 ms, a video clip showing the movement of an index finger or middle finger being lifted was presented for 300 ms. Next, the same static picture was presented for 1,000 ms. Finally, a black dot was presented on the static picture for 500 ms. When the black dot appeared, the subjects were tasked with mirroring the finger movement from the video clip with the right hand. The start of the video clip was defined as 0 ms. The period from −500 ms to −200 ms was referred to as the resting period. MEG signals from 0 ms to 900 ms were analysed in three time windows of 300 ms each for action observation, and signals from 1,300 ms to 2,100 ms were analysed in three time windows of 300 ms each for action execution.
peak in the left IPL, ERD was dominantly observed for the β-band (Fig. 2b). Therefore, commonly for the action observation and execution, the left IPL showed a significant ERD of the β-band.

Correlation of ERS/ERD between twins.
At the peak location of ERD/ERS among all subjects, the powers of each subject were compared among subject pairs in monozygotic and dizygotic twins and in unrelated pairs. The correlation coefficients of the powers between the pairs were evaluated for each of the three groups (intra-class correlations). Table 3 summarizes the intra-class correlations for each frequency band at each peak during three time windows of action observation. At 300-600 ms, the ERD of the β-band was significantly correlated among the monozygotic twin pairs (p = 0.036, Bonferroni-corrected, Fig. S1). The intra-class correlation of β-band ERD was significantly higher for monozygotic twins than for the unrelated pairs (Z-test, p = 0.027, Fig. 3). Also, for ERDs of the α-band, the intra-class correlations were significant at the right SPL for 300-600 ms (p = 0.00068, Bonferroni-corrected) and at the left IPL for 600-900 ms (p = 0.00048, Bonferroni-corrected) among monozygotic twins (Supplementary Table S1). These intra-class correlations among monozygotic twins were significantly higher than those of the unrelated pairs and dizygotic pairs (Z-test, p < 0.001). Notably, for other frequency bands of θ, low γ and high γ, the ERS/ERDs were not correlated among monozygotic twin pairs. The intra-class correlations were significantly high for the ERD of α and β bands among the monozygotic twins.
Similarly, Fig. 4 shows the intraclass correlations of three groups during action execution. The intra-class correlation was significantly high for ERDs of the α-band at the left postcentral gyrus (postCG) (p = 0.00090, Bonferroni-corrected) and for ERD of the β-band at the left IPL (p = 0.022, Bonferroni-corrected) among monozygotic twins. The intra-class correlations among monozygotic twins were significantly higher than those of the unrelated pairs and dizygotic pairs (Z-test, p < 0.001).
In addition to intra-class correlation, we examined absolute differences of ERS/ERD power between the twin pairs during action observation and execution (Supplementary Tables S2 and S3). During action observation, the absolute difference in β-band ERD between monozygotic twin pairs was significantly smaller than that of unrelated pairs (p = 0.0091, Bonferroni-corrected). Other differences, in the θ-, α-, low γand high γ-bands, between monozygotic pairs tend to be smaller than those between dizygotic and unrelated pairs, but the difference is not significant. Also, during action execution, absolute differences of the αand β-band ERDs between monozygotic twin pairs were significantly smaller than those of unrelated pairs (α-band ERDs: p = 0.000069, β-band ERDs: p = 0.0033).

Discussion
In this study, we identified significant ERSs and ERDs in each frequency band at various locations during action observation and execution. Among them, the βand low γ-band ERDs at the IPL during action observation had higher similarity among monozygotic twins than among dizygotic twins and unrelated pairs. Similarly, the α-band ERD at the postcentral gyrus and the β-band ERD at the IPL during action execution had high similarity among monozygotic twins. These results demonstrate that genetic factors influence brain activity during action observation and execution depending on the frequency and spatial differences of the activities. We observed significant ERSs/ERDs of the θand α-bands at the pre-and postcentral gyrus during action execution. These ERSs/ERDs are the cortical oscillatory changes involved in movement. In particular, the heritability of the powers was high for α-band ERD at the left postcentral gyrus (77%) and for β-band ERD at the left IPL (65%). In a previous twin study, we found more similarities in the latency and current source location of movement-evoked field 1 in monozygotic twins compared with unrelated pairs 19 . In the present study, we compared spatial and frequency patterns of MEG signals among monozygotic twins, dizygotic twins and unrelated pairs, to evaluate more precisely the genetic and environmental influences on movement-related cortical activities. By combining these results, we have evidence strongly supporting that movement-related cortical activities are affected by genetic influences.
It is noteworthy that the ERD of the left IPL is commonly observed both during action observation and execution in the β-band. It is known that the left IPL is part of the mirror neuron system 3,27 . This system contributes to the imitation of other people's actions 5,28 as well as to the perception of the intentions of others 4 . Thus, the task of imitating finger movements in this study may have activated the mirror neuron system in the left IPL. In addition, the heritability of β-band ERD powers at the left IPL was about 65% during both action observation and execution. This result suggests that the power of β-band ERD, which represents the activities of the mirror neuron system, is more influenced by genetic factors than by environmental factors.
Previous studies demonstrated that EEG powers are significantly correlated among monozygotic twins regardless of the cortical area 17,29 . However, the similarity in EEG signals may be due to the similarity of brain construction between twins (e.g. skull and cerebrospinal fluid) because EEG signals are influenced by differences in conductivity as determined by brain constructions. However, the MEG signal is theoretically not influenced by differences in conductivity. Moreover, the correlation was significant for the β-power band at the left IPL but not for other ERS/ERD. However, the other frequency band powers showed significant ERS/ERD in other cortical areas. Therefore, our results demonstrate that cortical activation during the cognitive task is influenced by genetic factors and depends on the cortical areas and frequency bands. The genetic influence is especially significant for the β power at the IPL during action observation.
Genes acting as genetic factors for the cortical activities have so far not been identified. A previous study has demonstrated that catechol-O-methyltransferase (COMT) Val/Met polymorphism was associated with variations in ERD/ERS responses using MEG during a working memory task 30 . In that study, Val homozygotes, as compared to Met homozygotes, increased ERD in the β-band in the left frontal-parietotemporal area. It is reasonable to hypothesise that such genetic factors influence the β-band ERD associated with the mirror neuron system.
Since we did not evaluate motor performance in this study, the relationship between brain oscillations and the performance of actual movement is unknown. In the future, to unveil the genetic factors of motor performance correlated with the cortical activities, detailed evaluation of the motor performance would be necessary.
Recently, we showed that the power of the low-γ ERD in the left frontal area during a language task was more homogeneous between monozygotic twins than between dizygotic twins 23 . However, the intra-class correlation of the β power at the left IPL during action observation in the present study (intra-class correlation of monozygotic twins = 0.60, intra-class correlation of dizygotic twins = 0.21) was larger than that of the previous study (intra-class correlation of monozygotic twins = 0.42, intra-class correlation of dizygotic twins = 0.24). These results indicate that cognitive function, such as action observation, is more innate than language function and is affected by genetics. The heritability of cortical activity varies depending on frequency, anatomical location and task. MEG studies involving twins can reveal the hereditary relationship between brain functions and the corresponding cortical activity. Our results demonstrate that genetic factors strongly affect the cortical activities for higher cognitive function.
It should be noted that the sample size of this study is smaller than in other twin studies. However, we have previously demonstrated that the cortical activities measured by MEG were highly correlated among twins during simple movement 19 and language tasks 23 . In this study, the ERD during the observation of action was dominant for β bands (Fig. 2b). These results strongly support that the sample size of this study is sufficient to evaluate the particular MEG responses for the present task. Although the signal-to-noise ratio (SNR) might be improved with a larger sample, and this would allow showing the hereditability of other frequency bands, it might be difficult to evaluate the power changes of frequency bands higher than the β and low γ bands due to the low SNR in MEG signals. Our analysis of the hereditability of the β ERD was done with a sufficient SNR and sample size to detect genetic influences.

Subjects. Subjects participated voluntarily and registered with the Center for Twin Research at Osaka
University 31 . We assessed 20 monozygotic twins (six males and 14 females; mean age, 58.1 ± 11.0 years) and 11 dizygotic twins (five males and six females; mean age, 61.4 ± 17.3 years). In all dizygotic twin pairs, both individuals of each twin pair were of the same sex. Zygosity was confirmed using the 15 loci of short tandem repeat markers derived from the blood 32 . Twin pairs that agreed completely for these short tandem repeat markers were designated as monozygotic, whereas all other pairs were designated as dizygotic.
SCientifiC RepoRts | 7: 17660 | DOI:10.1038/s41598-017-17662-x All subjects demonstrated normal comprehension and intellectual capacity. All subjects were right-hand dominant per the Edinburgh Handedness Inventory 33 . Written informed consent was obtained from all subjects after an explanation of the purpose and possible consequences of the study. All experimental procedures were performed in accordance with the protocols approved by the ethics committee of the Osaka University Graduate School of Medicine (No. 10190). The study was conducted in accordance with the Ethical Guidelines for Medical and Health Research Involving Human Subjects issued by the Ministries of Education, Culture, Sports, Science and Technology (MEXT) and Health, Labour and Welfare (MHLW) in Japan.

Measurements. A 160-channel whole-head neuromagnetometer device was used (MEG vision NEO; Ricoh
Company, Ltd., Yokohama, Japan). Visual stimuli were displayed on a projection screen 325 mm from the subject's eyes using a visual presentation system (Presentation; Neurobehavioral Systems, Berkeley, CA, USA) and a liquid-crystal projector (LVP-HC6800; Mitsubishi Electric Corporation, Tokyo, Japan).
The MEG data were recorded using an online band-pass filter between 0.01 and 200 Hz and sampled at 1,000 Hz. We also recorded electrooculograms and electrocardiograms (ECG) to remove the artefacts offline and surface electromyograms (EMG) of the right extensor indicis and extensor digitorum muscles to monitor finger movements. A 3D-T1 volume brain MRI was performed for all subjects using a 3.0-T magnetic resonance scanner.
A 3D digitiser (FastSCAN Corba, Polhemus, Colchester, VT, USA) and five head localisation coils attached to the subject's head were used to align the MEG data with the individual brain MRI results. The data of the 3D digitiser were superimposed on the individual's MRI data.
Tasks. After the presentation of the resting left hand image for 1,200 ms, the subjects observed the video clip of the index finger or middle finger being lifted without moving their own fingers for 300 ms. The resting hand image was then presented again for 1,000 ms; thereafter, a black dot was presented on the image between index and middle finger. When this dot appeared, the subjects lifted their own finger immediately to emulate the movements presented in the preceding video clip (Fig. 1). Using EMG activity, we assessed whether the subject moved the correct finger. This sequence was repeated 70 times. The order in which to move either finger was randomized.
Data analyses. MEG data were analysed using standard MEG software (MEG Laboratory; Ricoh Company, Ltd., Yokohama, Japan). We used an offline band-pass filter between 0.1 and 100 Hz and a band-elimination filter between 58 and 62 Hz (AC components). Eye blink artefacts on the data were separated and rejected by signal-space projection, a method for separating external disturbances, used in the Brainstorm software (version 3.2, http://neuroimage.usc.edu/brainstorm). Using signal-space projection, we detected obvious eye blinks from channels of the bilateral frontal base and rejected their artefacts from the data. We defined the onset time of the finger being lifted as 0 ms. The time window of a trial was defined as −500 to 2,100 ms. Among the 70 trials for each subject, unsuitable trials that included excessive muscle activity were not used in the analyses.
Beamformer. To calculate ERS/ERD, the data of each subject were analysed using a narrow-band adaptive beamformer, which is a spatial filtering method to improve the spatial resolution of neuromagnetic sources 34 . We used the beamformer software (coohea; Ricoh Company, Ltd., Yokohama, Japan). Estimates of differences in the source power between baseline and the period of interest for the selected time windows and frequency bands were computed as pseudo-t-values 20 . The distribution of the pseudo-t-values was then superimposed onto the individual anatomical MRI results coregistered to the MEG data. Beamformer created a volume that covered the whole brain of each subject with a voxel size of 5 × 5 × 5 mm 3 . The baseline was from −500 to −200 ms. The periods of interest during action observation were three 300-ms windows totalling 0 to 900 ms (0-300, 300-600 and 600-900 ms). The periods of interest during action execution were three 300-ms windows totalling 1,300 to 2,100 ms (1,300-1,600, 1,600-1,900 and 1,900-2,100 ms). The MEG data were divided into θ-(3-8 Hz), α-(8-13 Hz), β-(13-25 Hz), low γ-(25-50 Hz) and high γ-(50-100 Hz) bands by fast Fourier transformation (FFT).
Group statistical analysis. The cortical areas commonly activated among subjects were identified by group statistical maps. We used the statistical parametric mapping software package SPM8 (http://www.fil.ion.ucl.ac.uk/ spm/) to examine the statistical significance of ERS/ERD across subjects. The coregistered MRI data were normalised to the Montreal Neurological Institute (MNI) template with a voxel size of 2 × 2 × 2 mm 3 . A one-sample t-test at each voxel level was performed using a t-value incorporating variance smoothing with a 10-mm Gaussian kernel width. A p-value of < 0.001 (family-wise error rate) in each voxel was considered a statistically significant difference. Significant voxels were superimposed onto the MNI template brain using SPM8. We analysed only the surface voxels of the brain.

Virtual sensors and time-frequency analysis.
After identifying the spatiotemporal peaks of significant ERS/ERD during action observation using group statistical analysis, the cortical currents at the peak location in each frequency band were estimated by a virtual sensor toolbox (BRAM; Ricoh Company, Ltd., Yokohama, Japan). Next, a time-frequency analysis was applied to the estimated cortical currents.
The virtual sensor projects the signals exclusively from the targeted voxel using the adaptive beam-forming technique 35 . The beamformer is constructed to project signals exclusively from targeted voxels while eliminating residual noise to suppress signals from other parts of the brain. Thus, the virtual sensors provide data regarding neural activity at the target voxels with a considerably higher SNR than that of raw MEG data. The virtual sensor data were estimated by the individual MEG data divided into θ (3-8 Hz), α (8-13 Hz), β (13-25 Hz), low γ (25-50 Hz) and high γ (50-100 Hz) bands by FFT.
First, we determined the peak coordinates of the largest t-value for each frequency band using group statistical analysis. Next, the selected coordinate was used as the target location for the virtual sensor. The virtual sensor's SCientifiC RepoRts | 7: 17660 | DOI:10.1038/s41598-017-17662-x coordinates were extracted from individual MRI results based on the MNI coordinates and the warping parameters were calculated by SPM8 using an MRI-T1 template and individual MRI-T1 images. A tomographic reconstruction of the data was created by generating a single-sphere head model based on the head shape obtained from the structural MRI results for each subject.
A time-frequency analysis was applied to these signals to determine the oscillatory changes at the location of the virtual sensor in each subject, which was performed using FFT in the EEGLAB software (version 12.0.2.0, https://sccn.ucsd.edu/eeglab/) with MATLAB [version 8.2.0.701 (R2013b); The MathWorks Inc.]). Event-related spectral perturbations (ERSPs) were calculated while the window of 0.511 s overlapped 400 points (every 0.4 s) with the 1000 Hz sampling. Values at each time-frequency point were averaged over the trials. The baseline period extended from −500 to −200 ms. The frequency range of 1-100 Hz and the data windows were zero-padded with a ratio of 4. During the observation of action, the average power of ERSPs was calculated for each of the five frequency bands for each estimated cortical current from the 0-to 900-ms poststimulus time window with three 300-ms intervals. In the same way, during action execution, the average power of ERSPs was calculated from the 1,300-to 2,100-ms poststimulus time window with three 300-ms intervals.

Calculation of correlation coefficients and genetic analysis.
To compare the similarities in the power of ERS/ERD within twin pairs from both groups (monozygotic and dizygotic twins), we calculated the intra-class correlation between ERS/ERD power for the members of each twin pair and unrelated pairs. To apply the ACE model, the averages and variances of ERS/ERD power in twin1 (x axis) and twin2 (y axis) should not be significantly different 36 . We, therefore, randomly selected twin1 and twin2 so that there was no significant difference in the average and variance value of ERS/ERD power in twin1 and twin2 pairs. The group of 27 unrelated pairs was constructed by randomly selecting 54 subjects from all subjects of the same sex with age differences within 5 years. The unrelated pairs were generated 100 times to evaluate the intra-class correlation. Next, the intra-class correlations were averaged to estimate the intra-class correlation of the sex-and age-matched unrelated pairs. In addition, we calculated absolute differences of ERS/ERD power between the twin pairs. We compared these absolute differences among monozygotic or dizygotic twins to the differences of all unrelated pairs from each group (363 pairs in monozygotic twins and 124 pairs in dizygotic twins of same sex).
The genetic influences on individual differences were estimated by an ACE model to the powers of ERS/ERD in which the observed variance of the phenotype of interest is decomposed into genetic and environmental components. These components include variance caused by additive genetics effects (A), common environmental effects (C) and individual environmental effects including measurement error (E). The relative contribution of these latent factors is expressed by the path coefficients a, c and e, which are estimated by maximum likelihood. To estimate the values of a, c and e, the ERS/ERD data of the monozygotic and dizygotic twin pairs were analysed using the OpenMX program (http://openmx.psyc.virginia.edu/) and R statistical software.