Abstract
Biological rhythms are natural, endogenous cycles with period lengths ranging from less than 24 h (ultradian rhythms) to more than 24 h (infradian rhythms). The impact of the circadian rhythm (approximately 24 h) and ultradian rhythms on spectral characteristics of electroencephalographic (EEG) signals has been investigated for more than half a century. Yet, only little is known on how biological rhythms influence the properties of EEG-derived evolving functional brain networks. Here, we derive such networks from multiday, multichannel EEG recordings and use different centrality concepts to assess the time-varying importance hierarchy of the networks’ vertices and edges as well as the various aspects of their structural integration in the network. We observe strong circadian and ultradian influences that highlight distinct subnetworks in the evolving functional brain networks. Our findings indicate the existence of a vital and fundamental subnetwork that is rather generally involved in ongoing brain activities during wakefulness and sleep.
1 Introduction
While describing natural complex dynamical systems is a notoriously difficult endeavor, the network approach (; ; Newman, 2018) has been repeatedly shown to provide novel and important insights into such systems in various research areas ranging from neurosciences (; ) via genomics (Tyler et al., 2009) and proteomics (Uetz et al., 2000) to ecology (; Olesen et al., 2011; ; ), climatology (; Zhou et al., 2015), and sociology (Onnela et al., 2007; Palla et al., 2007). This broad applicability is not least explained by the large manifold of network metrics, describing global aspects to local aspects in network terms, which in principle can be directly related to the properties of the described system. Identifying key network constituents is highly relevant when it comes to improving the understanding and control of networks, as it allows us to gain insights about the importance hierarchy of its constituents with respect to the network structure and dynamics. The characterization of a constituent’s role in the network structure and dynamics can be achieved through different concepts and a growing number of metrics such as centralities. Most of these concepts focus on the description of vertices or groups of such [e.g., hubs (Newman, 2003), hub regions in the brain (Stanley et al., 2013; ), and k-core decompositions ()], while only a few metrics assess the centrality of an edge. We recently proposed a novel strength-based edge-centrality concept () and introduced modifications to vertex closeness and vertex eigenvector centrality concepts yielding corresponding edge centrality concepts (). We demonstrated that these edge centralities—together with edge betweenness centrality (; )—provide additional information about the network constituents for various topologies. These four centrality concepts, while different in their definition, can be considered complementary in the description of a constituent’s structural integration in the network.
Many early studies assumed networks to be static; however, the recent paradigm shift toward time-dependent (or evolving) networks (; ) allows one to describe many systems more accurately. This particularly holds true for biological networks, such as the brain, for which time dependencies on different scales have been observed. Both exogenous and endogenous biological rhythms are expected to assert influences on the level of the network description (; ; Mitsis et al., 2020; ; ) and therefore on network metrics such as centralities (; ; ; ).
For more than 50 years, it has been known that the circadian rhythm and ultradian rhythms impact electroencephalographic (EEG) signals [see for a recent overview]. Many former studies, however, were based on EEG recordings that either assessed the dynamics of few brain regions only or/and covered timescales ranging only from few seconds to hours. Here, we extend the recent studies and observations (Spoormaker et al., 2011; ; Park et al., 2012; ; ) and investigate how biological rhythms, particularly the circadian rhythm (with a period length of approximately 24 h), influence the importance hierarchies of the constituents of evolving functional brain networks. Therefore, we focus on both the networks’ vertices that are associated with the sampled brain regions and networks’ edges that represent time-evolving interactions between brain regions.
2 Materials and methods
2.1 Data
We analyzed electroencephalographic signals obtained from eight subjects (three females, age 19–81 years) with (five subjects) and without disorders (three subjects) of the central nervous system (CNS). All subjects were under stable CNS medication (if taking any). The EEG data were recorded continuously over 4 to 8 days from 19 electrodes placed according to the 10–20 EEG system () (Cz served as a physical reference) with a sampling rate of 256 Hz, using a 16-bit analog-to-digital converter (Micromed, S.p.A., Mogliano Veneto, Italy). Data were band-pass filtered offline (bandwidth: 1–45 Hz; fourth-order Butterworth characteristic), and a notch filter (third order) was used to suppress contributions at the line frequency (50 Hz). Data used in this study were visually inspected to remove segments containing strong artifacts (e.g., subject movements or amplifier saturation).
2.2 Deriving evolving functional brain networks
Time-dependent, fully connected, weighted functional brain networks were constructed through a time-resolved synchronization analysis of an EEG recording (Mormann et al., 2000; Osterhage et al., 2007; ; ) to track the changes in the importance hierarchies of network constituents (; Rings et al., 2019; ) possibly related to biological rhythms (). To perform this, network vertices were associated with brain regions whose dynamics were sampled by electrodes and network edges were associated with time-varying estimates of the strength of interactions between the dynamics of the pairs of those brain regions, regardless of their anatomical connections. As an estimate of the strength of the interaction, we employed mean phase coherence (R) (Mormann et al., 2000), which assesses the degree of synchronization between two phase time series (R = 1 indicates fully phase-synchronized brain regions, and R = 0 indicates no phase synchronization). A non-overlapping sliding window with a duration of 20 s (5,120 data points) was used to calculate R in a time-resolved manner. The chosen duration of a window can be considered a compromise between the required statistical accuracy for the calculation of R and the approximate stationarity within the window’s length (). For each window, we derived the instantaneous phase time series via the Hilbert transform of the EEG time series. An important property of this analytic signal approach (particularly in the case of two or more superimposed oscillatory components) is that the instantaneous frequency relates to the predominant frequency in the Fourier spectrum (). The predominant frequency may be subjected to fluctuations in the EEG time series. Thus, the instantaneous frequency can vary rhythmically around the predominant frequency, which results in spurious estimates of the instantaneous phase. By taking the temporal average, these effects can be reduced. From an electrophysiological point of view, we consider it more reasonable to look adaptively (e.g., via the Hilbert transform) at interactions between predominant rhythms in EEG rather than to look at interactions in some a priori fixed frequency bands (e.g., via wavelet transform) for which there is no power in the time series (Osterhage et al., 2007; ). Following these steps of analysis for each window, we end up with a temporal sequence of snapshot functional brain networks, each of which consists of V vertices and E edges and can be described by a weight matrix , where refers to the edge weight (strength of the interaction) between vertices i and j. The number of actual windows per subject depended on the respective recording duration, thus yielding approximately 9,500–26,000 windows.
2.3 Estimating the importance of network constituents
In order to further investigate the temporal sequence of snapshot functional brain networks, different approaches may be adopted. Estimating distance or (dis-)similarity between two networks might be one such approach, although finding suitable distance metrics still remains a challenge (Mheich et al., 2020). Another approach consists of the so-far insufficiently studied concept of multilayer networks (). Yet, due to some fundamental limitations, a meaningful interpretation of multilayer brain networks remains to be explored (). Therefore, here, we pursue the investigation of the time series of networks’ characteristics (; ) and utilize the following centrality metrics to estimate importance of each network’s vertices and edges.
Strength centrality of a vertex i is the sum of the weights of all edges connected to this vertex:The higher the vertex’s degree/strength, the more central it is considered to be. A related metric for edges is nearest-neighbor centrality, which considers an edge to be more central when its weight is larger and the strengths of the vertices that are connected by that edge are more similar and higher. Nearest-neighbor centrality of an edge k between vertices a and b is defined as follows ():where , (a, b) ∈ {1, …, V}, and denotes the weight of edge k connecting vertices a and b. Similar to strength centrality of a vertex, nearest-neighbor centrality of an edge is only influenced by its adjacent constituents. Hence, vertices (brain regions) and edges (interactions between pairs of brain regions) that have high values are largely interconnected with the neighboring vertices and edges. When structurally viewed, these neighboring constituents would be located in either the same areas of the brain or an adjacent area; from a functional perspective, they might be associated with any area of the brain. Hence, constituents with high values are highly interconnected within the functional network, although not necessarily allowing a structural interpretation.
Eigenvector centrality considers the influence of a vertex/edge on the network as a whole, where a vertex/edge is considered central if the vertex/edge adjacent to it is also central, and it is defined asIn case of vertices, (i, l) ∈ {1, …, V} and M denotes the weight matrix , where denotes the weight of an edge between vertices i and l. We define with . In the case of edges, (i, l) ∈ {1, …, E} and M denotes the weight matrix , whose entries are assigned to the average weight of edges i and l, if these edges are connected to the same vertex, and 0 otherwise. As mentioned previously, we define with . Equation 3 is applied iteratively until eigenvector centrality values remain stable. Hence, eigenvector centrality can be considered a strength-based centrality concept, which, due to its recursive definition, relates a vertex/edge to all the other vertices/edges in the network. Similar to strength/nearest-neighbor centrality, constituents with high values are gradiently stronger connected to closer constituents than to those that are far off. Again, distance-related descriptions, such as “close” or “far-off,” relate to the functional network and, hence, do not necessarily allow a structural interpretation, meaning constituents with large values are highly connected to the functional network in a rather general sense. This high inter-connectedness refers to many and/or possibly strong interactions with constituents either associated with the same brain area or possibly with any other area.
Closeness centrality considers the distance between a vertex/edge and all the other vertices/edges in the network. A vertex/edge with high closeness centrality is considered central as information obtained from this vertex/edge can reach all the other constituents in the network via short paths and so the vertex/edge can exert a more direct influence on the network. Closeness centrality of vertex k is defined as ()where and dnm represents the length of the shortest path between vertices n and m, which is calculated as the sum of the inverse of all edge weights on the path. Closeness centrality of edge k between vertices a and b can be defined as (where and (a, b, i) ∈ {1, …, V}. Hence, closeness centrality can be considered a path-based centrality concept, which is therefore influenced by the network as a whole. High closeness centrality points toward a constituent, which is associated with any brain area that is functionally “close” to any other constituent associated with any brain area. Hence, certain parts (with high ) of certain brain areas interact strongly with many other parts in the same area, while also interacting with many other brain areas.
Betweenness centrality is a measure of how frequently the shortest path traverses a given vertex/edge. A vertex/edge with high betweenness centrality is considered central because it acts as a bridge between other brain regions. Vertex/edge betweenness centrality (of vertex/edge i) can be defined as (; ; )where , , and ; qnm(i) represents the number of shortest paths between vertices n and m running through vertex/edge i, and Gnm represents the total number of shortest paths between vertices n and m. The length of a path is calculated as the sum of the inverse of all edge weights on that path. The normalization factor is given as F = (V−1) (V−2) in the case of vertices and F = V(V−1) in the case of edges. Hence, betweenness centrality can be considered a path-based centrality concept, which is therefore influenced by the network as a whole. Constituents with high betweenness centrality are likely to be part of bottleneck-like structures spanning between brain areas, both in a structural and functional sense.
In order to facilitate a qualitative comparision of the results obtained with the different centrality concepts, we utilized centrality value-based importance ranking of constituents (). A vertex/edge with the largest centrality value gets assigned to rank 1. The rank increases in increments of 1 for the vertex/edge with the second largest centrality value, third largest centrality value, etc., yielding an increasing rank as the centrality values decrease. This ranking can be further normalized, yielding a relative ranking with the highest relative rank being 1 (most important) and the lowest relative rank being 0 (least important). Hence, we can deduce an importance hierarchy for the vertices and edges of each snapshot functional brain network.
2.4 Characterizing the influence of biological rhythms on the importance of network constituents
The aforementioned steps of analysis provide us with a temporal sequence of vertex/edge importance hierarchies of an evolving functional brain network and enable the investigation of how biological rhythms impact this hierarchy. To this end, we proceed as described in and estimate the power spectral densities [Lomb–Scargle periodogram (Press et al., 1989)] of the respective temporal sequences. Eventually, we quantified the influence of the circadian rhythm on each such sequence as the portion of power for period lengths in the range of 20–28 h relative to the total power in the range of 1–36 h. We refer to this ratio as P24 in the following sections.
3 Results
We observe contributions of rhythms with period lengths of approximately 24 h (and to a lesser extent from rhythms of approximately 12 h or shorter) in all temporal sequences of importance (centrality values and ranks) of the respective network constituents from each subject. However, and contrary to the expectation, we found that these circadian and ultradian contributions are more pronounced in the sequences of some vertices and edges, i.e., some brain regions, as well as the interactions between their dynamics (see Figure 1).
FIGURE 1
Interestingly, we also observe that the strength of circadian contributions differs for temporal sequences derived with different centrality metrics. This is to be expected, at least to some extent, since the metrics highlight different structural aspects of a network, such as the path-structure or strength distribution. Nevertheless, these observations suggest that the circadian rhythm affects these different structural aspects. Figure 2 demonstrates that this rather unspecific relation regarding path/strength-based centrality metrics for vertices and edges can be observed in the data from all the investigated subjects. Moreover, it becomes quite apparent from this figure that by combining the results yielded by different centrality concepts, almost all network constituents are impacted by the circadian rhythm in all subjects.
FIGURE 2
Yet, it can also be observed that there is no trivial relation between the influence of the circadian rhythm (estimated with P24) and a constituent’s importance (cf.Figure 3). Neither for the most nor for the least important network constituents do we observe a generally specific influence of the circadian rhythm (as well as for ultradian rhythms (data not shown)).
FIGURE 3
Furthermore, we find that different centrality metrics identify different constituents as the most important (on average over the whole observation time, in line with previous studies) (see, e.g.,
FIGURE 4

Relation between constituents’ average relative rank over the total recording time and their fraction of the recording time, for which the constituents are deemed the most important. Marker colors encode the different subjects, and marker sizes encode the relative power corresponding to the 24-h peak related to the circadian rhythm.
In order to improve on the findings achieved so far, we investigate whether there exists a day/night pattern in the temporal evolution of the importance of vertices and edges (cf.Figures 1A, C). Interestingly, we observe that the largest differences in the importance of network constituents between night- and daytimes are related to very distinct brain areas along with their interactions (see Figure 5). These vertices and edges not only exhibit the largest change in centrality values when functional brain networks transit from night- to daytimes but are also further identified as the most important constituent on average and for the largest fraction of the total recording time. While betweenness centrality highlights bilateral frontotemporal vertices and edges, closeness centrality highlights the predominantly left temporoparietal vertices and edges. Both strength-based centralities (eigenvector centrality and nearest-neighbor centrality) predominantly highlight the left temporoparietal and left parietooccipital vertices and edges. Independent of the employed centrality metric, temporoparietal network constituents (vertices T7, T8, and P7 as well as their associated edges) are identified as the most important during daytime (here, 12:00 to 16:00 h). In contrast, during nighttime (here, 24:00 to 4:00 h), the importance shifts to parietal network constituents (vertices P3, P4, and PZ as well as their associated edges). Apart from this night–daytime-related spatial shift of the importance of a few network components, our findings also point toward a key arrangement of connected vertices and edges, relating to a subnetwork that comprises vertices T7, P7, P3 PZ, and P4 together with their associated edges. This subnetwork, which is slightly more dominantly located on the left brain hemisphere, is possibly involved in ongoing activities during wakefulness and sleep.
FIGURE 5

(A) Fraction of the recording time during which a network constituent is the most important (color-coded) and the average relative rank of the constituent (size-coded). Importance assessed with betweenness centrality , closeness centrality , eigenvector centrality , and strength/nearest-neighbor centrality . Middle column: Absolute night–day differences in the fraction of the recording time during which the respective constituent is the most important (color-coded, green/purple indicating a higher fraction of the recording time during night/day) and the absolute value of the absolute difference in the constituents’ relative rankings (size-coded). (B) Absolute night–day difference in constituents’ relative rankings (color-coded, red/blue indicating a higher relative ranking during night/day) and the absolute value of the absolute difference in the fraction of the recording time during which a network constituent is the most important (size-coded). The respective data are averaged over all subjects and their respective night/day periods.
4 Discussion
We investigated how the circadian rhythm impacts the time-dependent importance hierarchy of the vertices and edges of an evolving functional brain network. We employed different path- and strength-based centrality metrics for vertices and edges to comprehensively characterize the importance hierarchy of these network constituents. At the single-constituent level, we observed that their time-dependent changes in importance are clearly modulated by the circadian rhythm although to a varying degree. Irrespective of their interindividual variabilities, we observed pronounced differences in the constituents’ importance hierarchy when contrasting data obtained during daytimes and nighttimes. This may point to a local, circadian rhythm-driven modulation of the dynamics of various brain regions alongside their interactions. These brain regions form vital and fundamental subnetworks within the evolving functional brain networks.
The subnetwork highlighted with betweenness centrality comprises temporofrontal brain regions from both hemispheres. It is rather unexpected that the subnetworks, as highlighted with closeness, eigenvector, and strength/nearest-neighbor centrality metrics (but not with betweenness centrality), are largely overlapping, despite the fact that the different centrality metrics assess different structural aspects of network constituents. This subnetwork is predominantly restricted to the temporoparietal brain regions, with a left-hemispheric dominance during the nighttime. However, whether the subnetworks observed here are related to the resting-state network needs further investigation (Raichle, 2015).
Studies revealed that the hippocampus, deep inside the temporal lobe, and the visual cortex are simultaneously involved in the reactivation of coherent memory traces during sleep, which points toward a contribution to the memory consolidation process (Prabhakaran et al., 2000;
During the daytime, the vertices and edges associated with the parietal lobes (PZ, P3, and P4) are deemed more important and for a larger fraction of time compared to the nighttime. These areas consolidate spatial and visual information and integrate perceptions with other sensory inputs, resulting in the recognition of the trajectories of moving objects. These areas also mediate proprioception (perception of the position of the body in space) and are involved in skills such as arithmetic, writing, left–right orientation, and finger perception (see Rizzolatti et al. (1997) for an overview). Since these functions may also be involved during dream phases (which account for approximately 25% of the sleep period), these vertices and edges are nonetheless important during nighttimes although for shorter periods of time.
Overall, we observe that circadian (and ultradian) biological rhythms strongly influence the importance hierarchy, as assessed with different centrality concepts, of the constituents in time-dependent functional brain networks. These observations highlight, for each employed centrality concept, distinct subnetworks in evolving functional brain networks. The structural composition of these networks, however, largely coincides, which points toward the existence of a vital and fundamental subnetwork that is rather generally involved in ongoing brain activities.
Statements
Data availability statement
The data analyzed in this study are subject to the following licenses/restrictions: the datasets presented in this article are not readily available because they contain information that could compromise the privacy of research participants. Requests to access the datasets should be directed to the corresponding author.
Ethics statement
The studies involving human participants were reviewed and approved by the Ethics Committee of the University of Bonn. The patients/participants provided their written informed consent to participate in this study.
Author contributions
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Funding
This work was supported by the Verein zur Foerderung der Epilepsieforschung e.V. (Bonn).
Acknowledgments
The authors thank Christoph Helmstaedter and Thorsten Rings for fruitful discussions.
Conflict of interest
The authors RvW and KL declared that they were editorial board members of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
The remaining author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AlbouyG.KingB. R.MaquetP.DoyonJ. (2013). Hippocampus and striatum: dynamics and interaction during acquisition and sleep-related motor sequence memory consolidation. Hippocampus23, 985–1004. 10.1002/hipo.22183
2
BavelasA. (1950). Communication patterns in task-oriented groups. J. Acoust. Soc. Am.22, 725–730. 10.1121/1.1906679
3
BoashashB. (1992). Time frequency signal analysis: Methods and applications. Melbourne: Longman Cheshire.
4
BoccalettiS.BianconiG.CriadoR.Del GenioC. I.Gómez-GardenesJ.RomanceM.et al (2014). The structure and dynamics of multilayer networks. Phys. Rep.544, 1–122. 10.1016/j.physrep.2014.07.001
5
BoccalettiS.LatoraV.MorenoY.ChavezM.HwangD.-U. (2006). Complex networks: structure and dynamics. Phys. Rep.424, 175–308. 10.1016/j.physrep.2005.10.009
6
BrandesU. (2001). A faster algorithm for betweenness centrality. J. Math. Sociol.25, 163–177. 10.1080/0022250X.2001.9990249
7
BröhlT.LehnertzK. (2022). A straightforward edge centrality concept derived from generalizing degree and strength. Sci. Rep.12, 4407. 10.1038/s41598-022-08254-5
8
BröhlT.LehnertzK. (2019). Centrality-based identification of important edges in complex networks. Chaos29, 033115. 10.1063/1.5081098
9
BuldúJ. M.PapoD. (2018). Can multilayer brain networks be a real step forward?Comment on “network science of biological systems at different scales: A review” by M. Gosak et al. Phys. Life Rev.24, 153–155. 10.1016/j.plrev.2017.12.007
10
BullmoreE.SpornsO. (2009). Complex brain networks: graph theoretical analysis of structural and functional systems. Nat. Rev. Neurosci.10, 186–198. 10.1038/nrn2575
11
ChuC. J.KramerM. A.PathmanathanJ.BianchiM. T.WestoverM. B.WizonL.et al (2012). Emergence of stable functional networks in long-term human electroencephalography. J. Neurosci.32, 2703–2713. 10.1523/JNEUROSCI.5669-11.2012
12
ChungM. K. (2019). Brain network analysis. Cambridge, UK: Cambridge University Press.
13
De DomenicoM. (2017). Multilayer modeling and analysis of human brain networks. GIGASCIENCE6, 1–8. 10.1093/gigascience/gix004
14
DelmasE.BessonM.BriceM.-H.BurkleL. A.Dalla RivaG. V.FortinM.-J.et al (2019). Analysing ecological networks of species interactions. Biol. Rev.94, 16–36. 10.1111/brv.12433
15
DongesJ. F.ZouY.MarwanN.KurthsJ. (2009). Complex networks in climate dynamics. Eur. Phys. J.-Spec. Top.174, 157–179. 10.1140/epjst/e2009-01098-2
16
FarahaniF. V.FafrowiczM.KarwowskiW.BohaterewiczB.SobczakA. M.CeglarekA.et al (2021). Identifying diurnal variability of brain connectivity patterns using graph theory. Brain Sci.11, 111. 10.3390/brainsci11010111
17
FreemanL. C. (1977). A set of measures of centrality based on betweenness. Sociometry40, 35–41. 10.2307/3033543
18
FreiM. G.ZaveriH. P.ArthursS.BergeyG. K.JounyC.LehnertzK.et al (2010). Controversies in epilepsy: debates held during the fourth international workshop on seizure prediction. Epilepsy Behav.19, 4–16. 10.1016/j.yebeh.2010.06.009
19
FruengelR.BröhlT.RingsT.LehnertzK. (2020). Reconfiguration of human evolving large-scale epileptic brain networks prior to seizures: an evaluation with node centralities. Sci. Rep.10, 21921. 10.1038/s41598-020-78899-7
20
GeierC.LehnertzK.BialonskiS. (2015). Time-dependent degree-degree correlations in epileptic brain networks: from assortative to dissortative mixing. Front. Hum. Neurosci.9, 462. 10.3389/fnhum.2015.00462
21
GeierC.LehnertzK. (2017). Long-term variability of importance of brain regions in evolving epileptic brain networks. Chaos27, 043112. 10.1063/1.4979796
22
GirvanM.NewmanM. E. J. (2002). Community structure in social and biological networks. Proc. Natl. Acad. Sci. U.S.A.99, 7821–7826. 10.1073/pnas.122653799
23
GoodfellowM.AndrzejakR.MasollerC.LehnertzK. (2022). What models and tools can contribute to a better understanding of brain activity?Front. Netw. Physiol.2, 907995. 10.3389/fnetp.2022.907995
24
HalekotteL.FeudelU. (2020). Minimal fatal shocks in multistable complex networks. Sci. Rep.10, 11783. 10.1038/s41598-020-68805-6
25
HeglandS. J.NielsenA.LázaroA.BjerknesA.-L.TotlandØ. (2009). How does climate warming affect plant-pollinator interactions?Ecol. Lett.12, 184–195. 10.1111/j.1461-0248.2008.01269.x
26
HolmeP.SaramäkiJ. (2012). Temporal networks. Phys. Rep.519, 97–125. 10.1016/j.physrep.2012.03.001
27
KiveläM.ArenasA.BarthelemyM.GleesonJ. P.MorenoY.PorterM. A. (2014). Multilayer networks. J. Complex Netw.2, 203–271. 10.1093/comnet/cnu016
28
KlemG.LüdersH.JasperH.ElgerC. (1999). The ten-twenty electrode system of the international federation. The international federation of clinical neurophysiology. Electroencephalogr. Clin. Neurophysiol. Suppl.52, 3–6.
29
KongY.-X.ShiG.-Y.WuR.-J.ZhangY.-C. (2019). k-core: theories and applications. Phys. Rep.832, 1–32. 10.1016/j.physrep.2019.10.004
30
KuhnertM.-T.ElgerC. E.LehnertzK. (2010). Long-term variability of global statistical properties of epileptic brain networks. Chaos20, 043126. 10.1063/1.3504998
31
KuhnertM.-T.GeierC.ElgerC. E.LehnertzK. (2012). Identifying important nodes in weighted functional brain networks: A comparison of different centrality approaches. Chaos22, 023142. 10.1063/1.4729185
32
KurthJ. G.RingsT.LehnertzK. (2021). Testing jump-diffusion in epileptic brain dynamics: impact of daily rhythms. Entropy23, 309. 10.3390/e23030309
33
LehnertzK.AnsmannG.BialonskiS.DicktenH.GeierC.PorzS. (2014). Evolving networks in the human epileptic brain. Phys. D.267, 7–15. 10.1016/j.physd.2013.06.009
34
LehnertzK.GeierC.RingsT.StahnK. (2017). Capturing time-varying brain dynamics. EPJ Nonlin. Biomed. Phys.5, 2. 10.1051/epjnbp/2017001
35
LehnertzK.RingsT.BröhlT. (2021). Time in brain: how biological rhythms impact on EEG signals and on EEG-derived brain networks. Front. Netw. Physiol.1, 755016. 10.3389/fnetp.2021.755016
36
LiaoH.MarianiM. S.MedoM.ZhangY.-C.ZhouM.-Y. (2017). Ranking in evolving complex networks. Phys. Rep.689, 1–54. 10.1016/j.physrep.2017.05.001
37
LiuK. K.BartschR. P.LinA.MantegnaR. N.IvanovP. C. (2015). Plasticity of brain wave network interactions and evolution across physiologic states. Front. Neural Circ.9, 62. 10.3389/fncir.2015.00062
38
MheichA.WendlingF.HassanM. (2020). Brain network similarity: methods and applications. Netw. Neurosci.4, 507–527. 10.1162/netn_a_00133
39
MitsisG. D.AnastasiadouM. N.ChristodoulakisM.PapathanasiouE. S.PapacostasS. S.HadjipapasA. (2020). Functional brain networks of patients with epilepsy exhibit pronounced multiscale periodicities, which correlate with seizure onset. Hum. Brain Mapp.41, 2059–2076. 10.1002/hbm.24930
40
MormannF.LehnertzK.DavidP.ElgerC. E. (2000). Mean phase coherence as a measure for phase synchronization and its application to the EEG of epilepsy patients. Phys. D.144, 358–369. 10.1016/S0167-2789(00)00087-7
41
NewmanM. E. J. (2003). The structure and function of complex networks. SIAM Rev.45, 167–256. 10.1002/3527602755.ch2
42
NewmanM. (2018). Networks. Oxford University Press.
43
OlesenJ. M.BascompteJ.DupontY. L.ElberlingH.RasmussenC.JordanoP. (2011). Missing and forbidden links in mutualistic networks. Proc. Roy. Soc. B Biol. Sci.278, 725–732. 10.1098/rspb.2010.1371
44
OnnelaJ. P.SaramäkiJ.HyvönenJ.SzáboG.LazerD.KaskiK.et al (2007). Structure and tie strengths in mobile communication networks. Proc. Natl. Acad. Sci. U.S.A.104, 7332–7336. 10.1073/pnas.0610245104
45
OsterhageH.MormannF.StaniekM.LehnertzK. (2007). Measuring synchronization in the epileptic brain: A comparison of different approaches. Int. J. Bifurc. Chaos Appl. Sci. Eng.17, 3539–3544. 10.1142/s0218127407019330
46
PallaG.BarabásiA.-L.VicsekT. (2007). Quantifying social group evolution. Nature446, 664–667. 10.1038/nature05670
47
ParkB.KimJ. I.LeeD.JeongS.-O.LeeJ. D.ParkH.-J. (2012). Are brain networks stable during a 24-hour period?Neuroimage59, 456–466. 10.1016/j.neuroimage.2011.07.049
48
PrabhakaranV.NarayananK.ZhaoZ.GabrieliJ. D. E. (2000). Integration of diverse information in working memory within the frontal lobe. Nat. Neurosci.3, 85–90. 10.1038/71156
49
PressW. H.FlanneryB.TeukolskyS.VetterlingW. (1989). Numerical recipes in Pascal: The art of scientific computing. Cambridge, UK: Cambridge Univ. Press.
50
RaichleM. E. (2015). The brain’s default mode network. Annu. Rev. Neurosci.38, 433–447. 10.1146/annurev-neuro-071013-014030
51
RingsT.von WredeR.LehnertzK. (2019). Precursors of seizures due to specific spatial-temporal modifications of evolving large-scale epileptic brain networks. Sci. Rep.9, 10623. 10.1038/s41598-019-47092-w
52
RizzolattiG.FogassiL.GalleseV. (1997). Parietal cortex: from sight to action. Curr. Opin. Neurobiol.7, 562–567. 10.1016/s0959-4388(97)80037-2
53
SpoormakerV. I.CzischM.MaquetP.JänckeL. (2011). Large-scale functional brain networks in human non-rapid eye movement sleep: insights from combined electroencephalographic/functional magnetic resonance imaging studies. Phil. Trans. Roy. Soc. A Math. Phys. Eng. Sci.369, 3708–3729. 10.1098/rsta.2011.0078
54
StanleyM. L.MoussaM. N.PaoliniB.LydayR. G.BurdetteJ. H.LaurientiP. J. (2013). Defining nodes in complex brain networks. Front. Comput. Neurosci.7, 169. 10.3389/fncom.2013.00169
55
TylerA. L.AsselbergsF. W.WilliamsS. M.MooreJ. H. (2009). Shadows of complexity: what biological networks reveal about epistasis and pleiotropy. Bioessays31, 220–227. 10.1002/bies.200800022
56
UetzP.GiotL.CagneyG.MansfieldT. A.JudsonR. S.KnightJ. R.et al (2000). A comprehensive analysis of protein–protein interactions in Saccharomyces cerevisiae. Nature403, 623–627. 10.1038/35001009
57
ZhouD.GozolchianiA.AshkenazyY.HavlinS. (2015). Teleconnection paths via climate network direct link detection. Phys. Rev. Lett.115, 268501. 10.1103/PhysRevLett.115.268501
Summary
Keywords
functional brain network, vertex centrality, edge centrality, circadian rhythm, electroencephalographic signals
Citation
Bröhl T, von Wrede R and Lehnertz K (2023) Impact of biological rhythms on the importance hierarchy of constituents in time-dependent functional brain networks. Front. Netw. Physiol. 3:1237004. doi: 10.3389/fnetp.2023.1237004
Received
08 June 2023
Accepted
09 August 2023
Published
29 August 2023
Volume
3 - 2023
Edited by
Plamen Ch. Ivanov, Boston University, United States
Reviewed by
Milan Paluš, Institute of Computer Science, Czechia
Marko Gosak, University of Maribor, Slovenia
Updates

Check for updates
Copyright
© 2023 Bröhl, von Wrede and Lehnertz.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Timo Bröhl, timo.bröhl@uni-bonn.de
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.