Citation: Bayazit Karaman, R. Murat Demirer, Coskun Bayrak, M. Mert Su. Modeling the Antipodal Connectivity Structure of Neural Communities[J]. AIMS Neuroscience, 2016, 3(2): 163-180. doi: 10.3934/Neuroscience.2016.2.163
[1] | Shinichi Tamura, Yoshi Nishitani, Chie Hosokawa . Feasibility of Multiplex Communication in a 2D Mesh Asynchronous Neural Network with Fluctuations. AIMS Neuroscience, 2016, 3(4): 385-397. doi: 10.3934/Neuroscience.2016.4.385 |
[2] | Shun Sakuma, Yuko Mizuno-Matsumoto, Yoshi Nishitani, Shinichi Tamura . Simulation of Spike Wave Propagation and Two-to-one Communication with Dynamic Time Warping. AIMS Neuroscience, 2016, 3(4): 474-486. doi: 10.3934/Neuroscience.2016.4.474 |
[3] | Yoshi Nishitani, Chie Hosokawa, Yuko Mizuno-Matsumoto, Tomomitsu Miyoshi, Shinichi Tamura . Effect of correlating adjacent neurons for identifying communications: Feasibility experiment in a cultured neuronal network. AIMS Neuroscience, 2018, 5(1): 18-31. doi: 10.3934/Neuroscience.2018.1.18 |
[4] | Shun Sakuma, Yuko Mizuno-Matsumoto, Yoshi Nishitani, Shinichi Tamura . Learning Times Required to Identify the Stimulated Position and Shortening of Propagation Path by Hebb’s Rule in Neural Network. AIMS Neuroscience, 2017, 4(4): 238-253. doi: 10.3934/Neuroscience.2017.4.238 |
[5] | Tien-Wen Lee, Gerald Tramontano . Connectivity changes following transcranial alternating current stimulation at 5-Hz: an EEG study. AIMS Neuroscience, 2024, 11(4): 439-448. doi: 10.3934/Neuroscience.2024026 |
[6] | Zeynep Kaya, Mohammadreza Soltanipour, Alessandro Treves . Non-hexagonal neural dynamics in vowel space. AIMS Neuroscience, 2020, 7(3): 275-298. doi: 10.3934/Neuroscience.2020015 |
[7] | Chacchu Bhattarai, Phanindra Prasad Poudel, Arnab Ghosh, Sneha Guruprasad Kalthur . The RET gene encodes RET protein, which triggers intracellular signaling pathways for enteric neurogenesis, and RET mutation results in Hirschsprung's disease. AIMS Neuroscience, 2022, 9(1): 128-149. doi: 10.3934/Neuroscience.2022008 |
[8] | Dirk Roosterman, Graeme S. Cottrell . Astrocytes and neurons communicate via a monocarboxylic acid shuttle. AIMS Neuroscience, 2020, 7(2): 94-106. doi: 10.3934/Neuroscience.2020007 |
[9] | Yoshi Nishitani, Chie Hosokawa, Yuko Mizuno-Matsumoto, Tomomitsu Miyoshi, Shinichi Tamura . Classification of Spike Wave Propagations in a Cultured Neuronal Network: Investigating a Brain Communication Mechanism. AIMS Neuroscience, 2017, 4(1): 1-13. doi: 10.3934/Neuroscience.2017.1.1 |
[10] | Ernest Greene . New encoding concepts for shape recognition are needed. AIMS Neuroscience, 2018, 5(3): 162-178. doi: 10.3934/Neuroscience.2018.3.162 |
A neural network, which has the property of a modular/colonial structure and controls the information behavior of a human, is governed with a set of organizational principles. The structure in this modular/colonial environment is hierarchical and connections are established through the mutually exclusive nodes. We hypothesize that the decision mechanism behind a cognitive process takes place at the highest hierarchical and spatial level. The hierarchy is present between both functional and structural connections such that the neural signal paths emerges in white matter and concludes at gray matter. Such understanding is analogous to the wiring principles of the cerebrum that white matter accommodates the myelinated axons, which carry impulses to the cerebral cortex. (This generalization is also supported by the widely-accepted convention that the white matter accommodates the myelinated axons, which carry impulses to the cerebral cortex.) In this context, we are trying to unify the cognitive decision-making mechanism rather than trying to distinguish it between certain regions of the brain. The complexity of the cognitive decision mechanism can be observed in patients after a brain surgery. Although removing some portions of the brain may stop the progression of diseases, the brain loses some of its previous capabilities that were not related to that region, after the surgery. This suggests that there is a higher-level hierarchy in decision-making. Note: The decision-making mentioned here should not be confused with the decisions related to conscious choices. We assume that storing or sending information between neural populations is an internal cognitive decision-making process that starts at the lowest hierarchical structural level and concludes at the highest level in the hierarchy before any related task is completed. The rationale behind this approach lies at the very foundation of how the brain operates. Brain is constantly active whether it is in resting or active state because of the continuous blood-flow even in the deep sleep stage known as N3, which can be observed with EEG. Therefore, the brain may operate in a much more complicated and unified fashion rather than depending on single regional interactions for cognitive decision-making.
Existing methods in neuroscience have made important advances in the last two decades with the introduction of functional Magnetic Resonance Imaging (fMRI), diffusion-MRI, and arterial spin labeling. Inductive reasoning constitutes the foundational background for these technologies that the final conclusion depends on examining the relationship between smaller fragments of a larger implication. The complexity of the relationship between structural and functional connectivity amongst the neural communities makes inductive reasoning extremely challenging in providing a solid understanding of the cognitive functionality. Although the methods have been proven to be effective in neural analysis, one fundamental disadvantage is in the defining of the size of the fragments. This is a tedious task because of the immense number (~100 billion) of neurons in the mammalian brain. Moreover, the brain operates non-linearly, inferring that the size of the communities is dynamic and the contributions of these communities do not have identical effects on the final conclusion. We believe the size, boundary, hierarchy and assignments/responsibilities of neural communities are dynamic, in opposition to the general concept of predetermined size of modules. Therefore, in spite of the general convention of identifying the closely coupled neurons as “neural modules”, we prefer a new term; “neural communities” because the word “community” not only symbolizes the dynamic relationships and associations but also emphasizes the hierarchical foundation of it. Therefore, a new approach based on deductive reasoning is developed for neural analysis. Such an approach will represent the activity of the neural populations more realistically because of the following reasons:
1. Hemispherical shape of the brain and tendency of the electric current to flow on the surface of a curved conductor
2. Proven modular connectivity between neural populations
3. Overall decision-making mechanism is at the highest spatial level due to Reason 1.
A schematic diagram showing the differences between inductive and deductive reasoning is provided in Fig. 1.
For the reasons explained above, EEG emerges as an excellent candidate for utilizing the deductive approach because of the data collection method. However, the existing EEG data analysis methods are not reliable for this task and an enhanced and innovative data analysis models associated with a proper tool is needed in order to:
1. Contemplate the voltage emerging from neural electrical activity via the limited cross-sectional area of the electrode in a single dimension, and
2. Collect EEG data from separate electrodes (also called channels or sensors) for simultaneous correlation in the third-dimensional complex domain.
Understanding physiological process in the brain requires formal system models. Therefore, the connectivity is a core component of such models. According to a recent study [5], “70% of all information within the cortical regions in the brain passed through only 20% of these region’s neurons”. In other words, the 20% of this ration are referred to as “Hub-neurons”. Recent studies conducted at Indiana University [5] Hub-neurons not only aid in understanding the information processing behavior of the cortex but also reveals some feedback on the impact of neurogenerative diseases that effects the network. Hub-neurons are known to be information rich and have a direct or indirect link in the process of forming the “antipodal connections”. Although the existence of Hub-neurons, which carry the majority of information between cortical regions in the brain, have been reported in [6,7,8,9], the first known mathematical model showing their existence was introduced by [8,9]. This study, therefore, proposes the introduction of a potential cognitive biomarker called Antipodal Hub Neurons (AHNs) along with the associated hypothesis called antipodal connectivity of neural communities, in an effort to provide a better perspective in cognitive dynamics (Fig. 2).
The complexity of the brain is complemented with the activities resulting from neural communication and coordination. Extracting details of the process is inadequate if a simple analysis is conducted. Therefore, a comprehensive analysis with a solid theoretical foundation should be conducted at different levels to capture the details. Hence, a theoretical model, based on the potentials of Gauss' Law and Borsuk-Ulam theorem, were adapted to project single dimensional EEG data to third-dimension by complex analysis thus, the detailed communication and coordination of neural communities at the highest spatial and hierarchical level can be examined.
Gauss’ Law states that "the total electric flux (ϕ) out of a closed surface is equal to the charge (Q) enclosed divided by the electrical permittivity (ϵ)", i.e. ϕ=Q/ϵ [10,11] (Fig. 3). The electrical charges inside a solid closed surface (for example; skull) will move freely until the excess charges gather on the surface of the conductor. The charged particles will continue their motion until the electric field inside the conductor will be zero, which is known as electrostatic equilibrium.
Generally speaking, Gauss’ Law [12] states that the positive electrical charge density inside a conductor sphere creates an electrical field pointing outward. A similar assumption can be made for the brain, which resembles a hemisphere, since a certain amount of electrical current is flowing through the surface in a given time due to the continuous activity of the neural communities. However, the electrostatic equilibrium can only be held for a very short period of time as a result of this perpetual activity. The cerebral cortex and the cerebrospinal fluid are electrically conductive [13]. Therefore, Gauss’ Law is an excellent representation of how the electrical charges gather on the outer surface of the conductive margin of the brain, creating an electric field (Fig. 3).
The antipodal connectivity hypothesis is based on the antipodal-points concept of the Borsuk-Ulam theorem[14]. The theorem describes the existence of a pair of points, which are 180 degree apart on a spherical surface, having the identical shared value at a given time. Accordingly, we hypothesize that this principle is also true between the neural communities. Hence, there exists a pair of neural communities connected through a set of hub-neurons that are apart at the highest spatial level and share identical information simultaneously. This is credible evidence indicating that a paired connection based on shared value exists between those communities. These pairs are defined as the AHNs and the associated sets of communications are defined as “antipodal connectivity of neural communities”. A schematic representation of this hypothesis is provided in Fig. 4.
Borsuk-Ulam Theorem is the one most commonly used in complex analysis and algebraic topology to extract details [15,16]. The formal definition of the theorem, which was introduced in 1933, states, "Every continuous mapping of n-dimensional sphere, Sn into n-dimensional Euclidean space, Rn, identifies a pair of antipodes". More specifically, for every continuous mapping f:Sn→Rn there exists a point x∈Sn such that f(x)=f(−x).
The behavior of the neural populations is intricate to understand with raw EEG data. Therefore, discovering the antipodal points to understand the dynamically cohesive nature of coupled neurons needs a preprocessing step, which is crucial for neural interactivity analysis and mapping [17,18]. For this step, the data needs to be projected into a higher dimension through the use of the Borsuk-Ulam theorem.
The identification of antipodal structures is a significant part of the dynamically cohesive nature of coupled neurons. For instance, if there is a change in the number of antipodal pairs for different activity stages of the brain, it suggests that there exists a pattern of cohesive cooperation between the communities. This is calculated by evaluating the existence of antipodal points between the different EEG channels at a given time (t). Furthermore, the existence of antipodal points is evaluated in the proceeding discreet-time-step (t+1) until all the selected window is processed (t∈{1,2,3,...,W}) where W is the window length. The analytic signals from multichannel EEG data are compared at every time-step (t=1,2,3,...,W) by using a mathematical model called PNT, which is based on Borsuk-Ulam theorem. The selection of W is based on the given stimulus (i.e. "resting" and "active").
Expanding on the Borsuk-Ulam theorem, the PNT model provides a mathematical formulation for analyzing the antipodal connectivity at the highest hierarchical level. PNT is designed to extract the structural patterns in the cognitive decision mechanism using EEG test data recorded with different stimuli. PNT is a three-step processing algorithm. Order of the steps is as follows; 1- Transformation, 2- Normalization, and 3- Projection. The phase information is obtained by using transformation, which is required for assessment of the synchronization between sensors. Additionally, Hilbert transformation yields to an analytic signal (R2), which has only positive frequencies. The analytic signal is free of negative frequency content which is physically meaningless. Normalization is used to confine the transformed data into a unit disk which enables uniform processing for all data sets. Finally, the normalized EEG signal is mapped on the unit sphere with the help of Projection. The resulting mathematical formulation for PNT process is given below.
In Eq. 1, EEG data, which includes the amplitudes, is converted into analytical signal using Hilbert transformation (H) to extract complex attributes channel by channel, f∈R→f∈C(D), i.e. [19].
f(x,t)H⟶f(x+iy,t) | (1) |
Driving from the formal definition, the function G, a closed unit disk on the upper hemisphere, can be defined by the function f in the complex domain as [15]:
f(x+iy,t)=G(x,y,√1−x2−y2,t) | (2) |
The general formula of the PNT method is given in Eq. (3).
GP{e}t{→N{e}W{→U{e}W{→X{e}W}}} | (3) |
where U:R→C, W is the window size, and t is the discrete time-step, which is between 1 and W. i.e., 1≤t≤W∈t={1,2,3,...,W}. e denotes the electrode number out of total number of electrodes. i.e., e={k,l|Fp1,Fp2,Cz,...,Oz,O2}.
The data collected with 14-channel EEG system to test the capabilities of PNT method is processed in a segmented fashion, where each data (signal) segment is composed of data points (p) (in our case we used 1024 data points, due to the stimulus used in the experimentation). Thus, X is a matrix with rows (r) (in our case we used 14 rows) and p columns, as shown in Eq. (4).
→X{e}W={x{e}1,x{e}2,...,x{e}W} | (4) |
where →X{e}W, raw EEG data, is transformed vector by vector using Hilbert Transform. Therefore, the output signal, U, is also a matrix, which has the same dimension as in Eq. (5).
→U{e}W={x{e}1+y{e}1,x{e}2+y{e}2,...,x{e}W+y{e}W} | (5) |
where x{e}t+y{e}t was defined as u{e}t , therefore the general formula of Transformed data as given in Eq. (6).
→U{e}W={u{e}1,u{e}2,...,u{e}W} | (6) |
According to Eq. (2), one fundamental constraints of projection onto the upper hemisphere is that the real and imaginary data must be smaller than 1. Therefore, the normalization is performed to prepare the data for projection. In other words, performing normalization for each vector produces an output signal, which is a unit circle in complex plane, as shown in Eq. (7).
→N{e}W={uη{e}1,uη{e}2,...,uη{e}W} | (7) |
Calculating the third dimension using the real and imaginary parts of data is to monitor the closely coupled neural activities. Therefore, function G:R2→S2 is defined. The main objective is to discover two electrode pairs having equal PNT values in three-dimensions at time, t. Hence,
Gt(xη{e}t,yη{e}t,√1−xη2t−yη2t) | (8) |
where
x{e}t : The real part of the data collected from sensor e at time t.
xη{e}t : The real part of the data collected from sensor e at time t “after normalization”.
y{e}t : The imaginary part of the data from sensor e at time t.
yη{e}t : The imaginary part of the data from sensor e at time t “after normalization”.
The GPt function is the simplified version of Eq. (8). The function requires two input values, which are real and imaginary parts of the specific data to calculate the third dimension value, as shown in Eq. (9.
GPt=√1−xη2t−yη2t | (9) |
Once the complex attributes of EEG data is extracted using Hilbert transform, they are projected onto the upper hemisphere using Eq. (9). This process is repeated for every channel at each time-step. Thus, the 3-D values are calculated in spatial-temporal domain from 1-D EEG amplitudes. When the degree of similarity between 3D values for channel-pairs are closely identical (≤0.00001) at a time-step, t, two remote communities at the highest spatial level (right below the electrodes) are assumed to have bridged a connection at the highest hierarchical level. Such neural communities are defined as AHNs. Hence, the following formal definition holds.
Definition-1: For all electrodes, there exist two different sensors (k and l), which have approximately the same 3-D values at time t (Eq. (10)).
∀e∃k,l∋GPt(uη{e|k}t)≅GPt(uη{e|l}t) | (10) |
With respect to the formal foundation established for PNT model, it is safe to make the following observation.
Observation-1: In the spatial-temporal complex domain, there exists a pair of sensors whose projected 3-D values are equal with 5-digit precision in terms of real and imaginary parts.
Since the approach and the mathematical model provided in this study have not been tested before, evaluation of the measurable outcomes is purely based on this experiment. We initially hypothesized that the amount of Antipodal Hub Neurons (AHN) varies according to the brain's current stage, which is “active” or “resting”. In order to test this hypothesis, we conducted this comprehensive experiment on normal subjects with different age groups and sexes. In our study we also hypothesized that the number and meaning of the AHNs can change with the frequency bands of the EEG data. The best known and used frequency bands of the EEG data is between 1-60 Hz, dividing the EEG data into five main bands, i.e., delta (1-4 Hz), theta (4-8 Hz), alpha (8-13 Hz), beta (13-30 Hz) and gamma (30-60 Hz) [20]. The waves, delta through beta, are associated with problem solving, memory functions, relaxation, and creativity [21].
For the experimental setup, a EEG system with 14 channels is used to record the data. It has the same properties of an international 10-20 system with a channel configuration as follows; AF3, F7, F3, FC5, T7, P7, O1, O2, P8, T8, FC6, F4, F8, and AF4. The EEGLAB which is an open-source and the most widely used toolbox in the world to analyze the EEG data [22] was used to filter the data .
The data collection environment is designed in such a way to avoid any noise and artifacts. Therefore, without causing any distraction a list containing 26 letters were laid out randomly on a cardboard display that was placed in front of the test subjects at a distance of approximately five feet (see Fig. 6).
After obtaining the IRB certification, 18 individuals of different age groups and sexes (11 males and 7 females) were recruited to participate in the data recording. Subjects were given small condiments 45 minutes prior to the experiment. Thus, the subjects are assumed to accomplish the tasks under similar physiological states. The age range for the male and female participants are 25-56 (average 28.1) and 21-45 (average 27.9), respectively.
As shown in Fig. 6 subjects were asked to follow a certain procedure to complete experiments. The steps were explained to the subjects before the experimentation in order to unify the testing procedure.
1. Subjects start the procedure with their eyes closed while listening to background music.
2. Data recording starts after a short while ensuring that the subjects were only concentrating on the music.
3. After 10 seconds, a trigger signal goes off notifying the subjects to open their eyes and start searching for the letters of a predetermined 5-letter meaningful word in order in the least possible time. Up to this point, subjects had not seen the letters on the cardboard hence; the process would require a certain level of challenge.
4. Upon finding the word on the cardboard subjects snap their fingers as a confirmation of accomplishing the task.
Calm classical music is used as background noise to help subjects relax and to foster an idle state of the brain as much as possible [23,24]. The initial 10-second period was determined as the control data where brain activity was assumed to be minimum/default. In this context, a quiet environment was not desired since the subjects might indulge in other thought processes. Three unique 5-letter words were chosen for the test, which consisted dissimilar letters for each word.
In total, 54 experiments were conducted. The participants were asked to find letters in the correct order under 1 minute. Although, there was no limitation for the duration of experiments, finding each of the 5-letter words lasted 8 seconds on average. During Experiment-3, participants 1, 3, 5, and 18 completed the task in an amount of time relatively shorter than required (≥8secs) therefore extracting meaningful information from that data set was not possible. Hence, they were excluded from the results due to the unreliable information they might carry. The collected data were analyzed using a FIR filter implemented in EEGLab [25] for noise reduction and band selection for the related brain activity.
Using PNT method, a total of 50 experiments were analyzed to show the existence of AHN network structure in the brain network. The significant increase in the number of AHNs in active stage for Experiment-2 with Subject-5, and Experiments 2 and 3 with Subject-6 were detected, these findings were defined as outliers and excluded from the result. Fig. 7 shows the histogram of PNT results for 47 experiments where significant increase in the number of AHNs were detected in “active” brain stage relative to the “resting” stage. The majority of AHNs are encountered in frontal and parietal cortex. The total number of AHNs for all the test subjects during the resting stage and active stage are 1564 and 2957, respectively. In other words, the average number of AHNs per experiment during the resting stage and active stage are 33.277 and 62.915, respectively. These results confirm that the AHNs increased approximately 2 times on average during the active stage. Additionally, the standard deviation of AHN counts for each experiment during the resting stage is 23.01, while it is 40.786 at active stage.
Considering the results presented in Fig. 8 and in Table 1, we encounter a similar conclusion in that all of the test subjects (mean) showed an increase in activity and number of AHNs while their brain was engaged with some tasks. Additionally, a randomly generated (Gaussian distribution) data set was processed with PNT method. For consistency, each data set was produced to contain an equal number of data points as the real EEG waveforms collected during the experiments. The AHN counts for the data set were significantly less than that of the EEG data collected from the subjects which verified that the occurrence of antipodal connectivity is not random but depends on the level of activity. This confirms the reliability of the proposed formal methodology. provided in this study.
Due to the non-consistency pattern, i.e., there were no consistent number of AHNs detected for the test groups (see Table 1) used in the experimentation, a statistical analysis was conducted to compare the significance of average numbers of AHNs encountered during the resting stage with active stage and resting stage with randomly generated data per subject.
Number of Antipodal Hub Neurons | |||
Subject Number | Randomly Generated | Resting Stage | Active Stage |
1 | 8 | 19 | 36 |
2 | 17 | 55 | 77 |
3 | 14 | 33 | 48 |
4 | 11 | 43 | 65 |
5 | 18 | 52 | 84 |
6 | 3 | 96 | 280 |
7 | 6 | 49 | 67 |
8 | 2 | 21 | 36 |
9 | 7 | 13 | 69 |
10 | 3 | 25 | 58 |
11 | 19 | 24 | 86 |
12 | 11 | 68 | 83 |
13 | 19 | 23 | 35 |
14 | 6 | 26 | 43 |
15 | 8 | 23 | 32 |
16 | 3 | 22 | 54 |
17 | 9 | 25 | 51 |
18 | 4 | 31 | 72 |
Since differences between data sets(randomly generated, active, and resting stages) individually were not normally distributed, the non-parametric statistical method (Mann-Whitney U test) was used. Firstly, the Mann-Whitney U test was applied to the resting and active stages data sets to reject the null hypothesis, since it states that two data sets are not significantly different at p≤0.01(probability).The default critical value of U is 81 when p≤0.01 and n=19. With respect to the comparative analysis performed on different data (resting-active) sets, the U-value calculated to be 54.5 which is smaller than the critical value of U at p≤0.01, hence the null hypothesis was rejected. Similarly, the random generated and resting stage data sets revealed U-value of 6 which is smaller than the critical value of U at p≤0.01. Nonetheless, the final outcome of the analysis indicated a significance differences between randomly generated data and resting stage data and the resting stage data with active stage.
It is worth noting that, an EEG system with more channels may be more efficient to capture more AHNs during the active stage, since the brain is operating in a non-linear fashion as mentioned “…the non-stationarity of the EEG signals can be observed during the change in alertness and wakefulness …” [26]. In other words, electrical activity of the brain varies with the level of processing.
Out of hundreds of AHN pairs detected by the PNT algorithm, an AHN pair captured by the Sensors, F4 and O2 is illustrated in Fig. 9. PNT detected this AHN pair at t = 1.66 seconds. Therefore, in order to test the validity of the PNT method, the amplitudes of the concurrent EEG values are compared to detect any similarities for the same sensors as a control measure.
The compared data sets show that the concurrent raw EEG amplitudes are dissimilar while the PNT values are identical with a precision of ten thousandths. Table 2 shows the comparison between the values of three AHNs detected with PNT (t=1.66, 12.65, 14.48 seconds) and the amplitudes of the raw EEG data for the same sensor pair (F4 and O2). These results confirm that raw EEG data is not sufficient to discover AHNs without using the PNT. Additionally, the increase in AHN counts in the active stage compared to that of the resting stage strongly suggests that there is a pattern. Although, such a pattern is not fully understood yet, the increase in the correlated pairs detected by PNT presents strong evidence for its existence.
1-31 Hz EEG Amplitudes | 3-D Values from PNT | |||
time (sec) | O2 | F4 | O2 | F4 |
1.66 | -8.029 | 5.654 | 0.976607143616050 | 0.976608882839799 |
12.65 | -2.001 | -4.622 | 0.984382629815184 | 0.984375119249569 |
14.48 | -0.3727 | 0.2089 | 0.999775830582554 | 0.999781126238373 |
As another control measure, the raw EEG values from all the channels are cross-correlated at each time-step to identify the potential similarities. As an example, raw EEG values are displayed for the concurrent AHN pair (Fig. 10). No correlation was detected between the raw EEG channels, which verifies that the PNT is capable of detecting the connections established between the remote neural communities, using complex analysis. These deductions were satisfactory to conclude that connections are established between the remote communities, which AHNs not only exist but also have significant importance.
There are a few other deductions worth mentioning. The processed data which was filtered between 1-31 Hz, is the most suitable band for the PNT algorithm, showed the difference in engagement of the neural communities in resting and active brain stages. The difference in AHN counts are not consistent with that of the 1-31 Hz band when the 1-60 Hz band is examined (i.e. gamma wave band is included). Therefore, our results support the idea that not all neuroscientists are convinced of the existence of gamma waves, which are thought to be the artifact of the electromyographic activity [27,28,29].
It was observed and previously mentioned that processing a shorter EEG signal might bias the results. Since EEG data is non-stationary [26], selection of the Hilbert window plays a crucial role in calculating the AHN counts. Selecting a shorter Hilbert window results in drastically fewer AHNs because it fails to accurately detect the changes in the frequency content of the EEG data for a specific stimulus. Since the method intended to compare the AHN counts during the resting and active stages, the resting period is characterized as one stimulus while the active period is characterized as another. Therefore, the Hilbert window size was selected according to the length of these stimuli for the experimentation. For example, doubling the size of the Hilbert window increased the AHN counts by three to four times in both stages. Therefore, the mathematical model provided in this study to validate the existence of the antipodal connectivity structure is satisfactory because not only does raw EEG data show no correlation but also antipodal connectivity structure conforms to the nonlinearity principle [26] of the brain's cognitive processing.
Neural activities are evolving in an organized fashion. In recent studies [6,8], confirmed this evolution in focusing on the cortical functional connectivity, known as default network [30,31,32,33,34,35,36] and the interactions in the spatial-temporal domain. While investigating the regional connectedness by modeling the interactions, we realized a unique set of dynamic coupling reflecting the existence of antipodal structures among regions.
Borsuk-Ulam’s antipodal symmetry theory was adapted to establish the mathematical foundation for the PNT model. The PNT was tested on the EEG data collected in 50 experiments. A test procedure was established to distinguish the AHN counts during resting and active stages of the brain. According to these results, AHNs increase when brain is engaged with tasks such as problem solving and/or focused attention. There are no certain number of AHNs for each task due to two immediate factors;
1. the limited number of EEG data sets,
2. unique processing characteristics of each brain.
There are possibly additional factors that play a key role in AHN counts. However; this will require more testing and more detailed investigation. The precision of the system is also a key factor in determining the occurrence rate of antipodal connectivity, which was kept at hundred thousandths (10-5).
The occurrence rate of antipodal connectivity is not random but depends on the stimuli and the selection of Hilbert transform window size. This is due to the non-stationary behavior of the electrical activity, which also causes antipodal connectivity to be a nonlinear process. Nonlinearity and lack of randomness of antipodal connectivity is substantially correlated with the previously proven theories. Our results indicated that the AHN counts are more reliable between the 1-31 Hz band. This frequency band corresponds to the frequency band of the EEG signals, delta through beta, associated with the different stages, such as mediation through problem solving, respectively, but the upper limit of this band is usually at 38-40 Hz (beta waves). This suggests that the activity bandwidth may actually be shorter.
The authors thank Dr. Eric Kaufmann in the Department of Mathematics and Statistics at UALR for his assistance in mathematical concepts. For this study, Human Subject Use authorization was obtained for the data recording.
[1] |
Chen ZJ, He Y, Rosa-Neto P, et al. (2008) Revealing modular architecture of human brain structural networks by using cortical thickness from MRI. Cerebral cortex 18: 2374-2381. doi: 10.1093/cercor/bhn003
![]() |
[2] |
Ferri R, Rundo F, Bruni O, et al. (2008) The functional connectivity of di erent EEG bands moves towards small-world network organization during sleep. Clinical Neurophysiology 119: 2026-2036. doi: 10.1016/j.clinph.2008.04.294
![]() |
[3] |
Meunier D, Achard S, Morcom A, et al. (2009) Age-related changes in modular organization of human brain functional networks. Neuroimage 44: 715-723. doi: 10.1016/j.neuroimage.2008.09.062
![]() |
[4] | Niedermeyer E, da Silva FL (2005) Electroencephalography: basic principles, clinical applications, and related fields: Lippincott Williams & Wilkins. |
[5] |
Nigam S, Shimono M, Ito S, et al. (2016) Rich-Club Organization in E ective Connectivity among Cortical Neurons. The Journal of Neuroscience 36: 670-684. doi: 10.1523/JNEUROSCI.2177-15.2016
![]() |
[6] |
Hagmann P, Cammoun L, Gigandet X, et al. (2008) Mapping the structural core of human cerebral cortex. PLoS Biol 6: e159. doi: 10.1371/journal.pbio.0060159
![]() |
[7] |
Rubinov M, Sporns O (2010) Complex network measures of brain connectivity: uses and interpretations. Neuroimage 52: 1059-1069. doi: 10.1016/j.neuroimage.2009.10.003
![]() |
[8] | Sporns O, Tononi G, Kotter R (2005) The human connectome: a structural description of the human. Brain PLos Comput Biol(2005) 1. |
[9] | Leergaard TB (2012) Mapping the connectome: multi-level analysis of brain connectivity: Frontiers E-books. |
[10] | Balanis CA (2012) Advanced engineering electromagnetics: Wiley Online Library. |
[11] | Feynman RP, Leighton RB, Sands M (2013) The Feynman Lectures on Physics, Desktop Edition Volume I: Basic Books. |
[12] | Serway R, Jewett J (2013) Physics for scientists and engineers with modern physics: Cengage learning. |
[13] |
Baumann SB, Wozny DR, Kelly SK, et al. (1997) The electrical conductivity of human cerebrospinal fluid at body temperature. Biomedical Engineering, IEEE Transactions on 44: 220-223. doi: 10.1109/10.554770
![]() |
[14] | Prescott T (2002) Extensions of the Borsuk-Ulam Theorem. Harvey Mudd College, http://citeseerx.ist.psu.edu. |
[15] | Browder A (2006) Complex Numbers and the Ham Sandwich Theorem. The American Mathematical Monthly: 935-937. |
[16] | Matousek J (2008) Using the Borsuk-Ulam theorem: lectures on topological methods in combinatorics and geometry: Springer Science & Business Media. |
[17] |
Aler R, Galvn IM, Valls JM (2012) Applying evolution strategies to preprocessing EEG signals for braincomputer interfaces. Information Sciences 215: 53-66. doi: 10.1016/j.ins.2012.05.012
![]() |
[18] | Yang S, Deravi F. Wavelet-based EEG preprocessing for biometric applications; 2013. IEEE. pp. 43-46. |
[19] | Hahn SL (1996) Hilbert transforms in signal processing: Artech House on Demand. |
[20] |
Yuvaraj R, Murugappan M, Ibrahim NM, et al. (2014) On the analysis of EEG power, frequency and asymmetry in Parkinsons disease during emotion processing. Behav brain Funct 10: 12. doi: 10.1186/1744-9081-10-12
![]() |
[21] | Sanei S (2013) Adaptive processing of brain signals: John Wiley & Sons. |
[22] | Iversen JR, Makeig S (2014) MEG/EEG Data Analysis Using EEGLAB. Magnetoencephalography: Springer. pp. 199-212. |
[23] |
Nater UM, Abbruzzese E, Krebs M, et al. (2006) Sex di erences in emotional and psychophysiological responses to musical stimuli. International journal of psychophysiology 62: 300-308. doi: 10.1016/j.ijpsycho.2006.05.011
![]() |
[24] |
KABUTO M, KAGEYAMA T, NITTA H (1993) EEG Power Spectrum Changes due to Listening to Pleasant Musics and Their Relation to Relaxation E ects. Nippon Eiseigaku Zasshi 48: 807- 818. doi: 10.1265/jjh.48.807
![]() |
[25] |
Delorme A, Makeig S (2004) EEGLAB: an open source toolbox for analysis of single-trial EEG dynamics including independent component analysis. Journal of neuroscience methods 134: 9-21. doi: 10.1016/j.jneumeth.2003.10.009
![]() |
[26] | Sanei S, Chambers JA (2013) EEG signal processing: John Wiley & Sons. |
[27] |
Whitham EM, Lewis T, Pope KJ, et al. (2008) Thinking activates EMG in scalp electrical recordings. Clinical Neurophysiology 119: 1166-1175. doi: 10.1016/j.clinph.2008.01.024
![]() |
[28] |
Whitham EM, Pope KJ, Fitzgibbon SP, et al. (2007) Scalp electrical recording during paralysis: quantitative evidence that EEG frequencies above 20Hz are contaminated by EMG. Clinical Neurophysiology 118: 1877-1888. doi: 10.1016/j.clinph.2007.04.027
![]() |
[29] |
Yuval-Greenberg S, Tomer O, Keren AS, et al. (2008) Transient induced gamma-band response in EEG as a manifestation of miniature saccades. Neuron 58: 429-441. doi: 10.1016/j.neuron.2008.03.027
![]() |
[30] |
Fox MD, Raichle ME (2007) Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging. Nature Reviews Neuroscience 8: 700-711. doi: 10.1038/nrn2201
![]() |
[31] | Shulman GL, Fiez JA, Corbetta M, et al. (1997) Common blood flow changes across visual tasks: II. Decreases in cerebral cortex. Journal of cognitive neuroscience 9: 648-663. |
[32] |
Raichle ME, MacLeod AM, Snyder AZ, et al. (2001) A default mode of brain function. Proceedings of the National Academy of Sciences 98: 676-682. doi: 10.1073/pnas.98.2.676
![]() |
[33] |
Greicius MD, Krasnow B, Reiss AL, et al. (2003) Functional connectivity in the resting brain: a network analysis of the default mode hypothesis. Proceedings of the National Academy of Sciences 100: 253-258. doi: 10.1073/pnas.0135058100
![]() |
[34] |
Fox MD, Snyder AZ, Zacks JM, et al. (2006) Coherent spontaneous activity accounts for trial-totrial variability in human evoked brain responses. Nature neuroscience 9: 23-25. doi: 10.1038/nn1616
![]() |
[35] |
Fox MD, Snyder AZ, Vincent JL, et al. (2005) The human brain is intrinsically organized into dynamic, anticorrelated functional networks. Proceedings of the National Academy of Sciences of the United States of America 102: 9673-9678. doi: 10.1073/pnas.0504136102
![]() |
[36] |
Fox MD, Corbetta M, Snyder AZ, et al. (2006) Spontaneous neuronal activity distinguishes human dorsal and ventral attention systems. Proceedings of the National Academy of Sciences 103: 10046-10051. doi: 10.1073/pnas.0604187103
![]() |
1. | Arturo Tozzi, The multidimensional brain, 2019, 31, 15710645, 86, 10.1016/j.plrev.2018.12.004 |
Number of Antipodal Hub Neurons | |||
Subject Number | Randomly Generated | Resting Stage | Active Stage |
1 | 8 | 19 | 36 |
2 | 17 | 55 | 77 |
3 | 14 | 33 | 48 |
4 | 11 | 43 | 65 |
5 | 18 | 52 | 84 |
6 | 3 | 96 | 280 |
7 | 6 | 49 | 67 |
8 | 2 | 21 | 36 |
9 | 7 | 13 | 69 |
10 | 3 | 25 | 58 |
11 | 19 | 24 | 86 |
12 | 11 | 68 | 83 |
13 | 19 | 23 | 35 |
14 | 6 | 26 | 43 |
15 | 8 | 23 | 32 |
16 | 3 | 22 | 54 |
17 | 9 | 25 | 51 |
18 | 4 | 31 | 72 |
1-31 Hz EEG Amplitudes | 3-D Values from PNT | |||
time (sec) | O2 | F4 | O2 | F4 |
1.66 | -8.029 | 5.654 | 0.976607143616050 | 0.976608882839799 |
12.65 | -2.001 | -4.622 | 0.984382629815184 | 0.984375119249569 |
14.48 | -0.3727 | 0.2089 | 0.999775830582554 | 0.999781126238373 |
Number of Antipodal Hub Neurons | |||
Subject Number | Randomly Generated | Resting Stage | Active Stage |
1 | 8 | 19 | 36 |
2 | 17 | 55 | 77 |
3 | 14 | 33 | 48 |
4 | 11 | 43 | 65 |
5 | 18 | 52 | 84 |
6 | 3 | 96 | 280 |
7 | 6 | 49 | 67 |
8 | 2 | 21 | 36 |
9 | 7 | 13 | 69 |
10 | 3 | 25 | 58 |
11 | 19 | 24 | 86 |
12 | 11 | 68 | 83 |
13 | 19 | 23 | 35 |
14 | 6 | 26 | 43 |
15 | 8 | 23 | 32 |
16 | 3 | 22 | 54 |
17 | 9 | 25 | 51 |
18 | 4 | 31 | 72 |
1-31 Hz EEG Amplitudes | 3-D Values from PNT | |||
time (sec) | O2 | F4 | O2 | F4 |
1.66 | -8.029 | 5.654 | 0.976607143616050 | 0.976608882839799 |
12.65 | -2.001 | -4.622 | 0.984382629815184 | 0.984375119249569 |
14.48 | -0.3727 | 0.2089 | 0.999775830582554 | 0.999781126238373 |