- 1Inner Mongolia Key Laboratory of Pattern Recognition and Intelligent Image Processing, School of Information Engineering, Inner Mongolia University of Science and Technology, Baotou, China
- 2School of Life Science and Technology, Xidian University, Xi'an, China
- 3Engineering Research Center of Molecular and Neuro Imaging, Ministry of Education, Xi'an, China
- 4Guangxi Key Laboratory of Multi-source Information Mining and Security, Guangxi Normal University, Guilin, China
Previous electroencephalogram (EEG) studies revealed reduced spectral power during the resting state in smokers. However, few studies investigated the changes of global brain networks during the resting state in young smokers by EEG. In the present study, we used minimum spanning tree (MST) to assess the differences of global network efficiency between young smoker (n = 20) and nonsmokers (n = 20). Compared with healthy nonsmokers, young smokers showed decreased leaf fraction, kappa value, increased diameter and eccentricity value in alpha band (r = 0.574, p = 0.008), which suggested the global network efficiency was decreased in young smokers. We also found positive correlation between leaf fraction and onset time of smoking in smokers. These results provided more scientific evidence of the abnormal neural oscillations of young smokers and improved our understanding of smoking addiction.
Introduction
Annual smoking-attributable deaths was about 6 million in the world including one million in China (World Health Organization, 2013). “China Report on the Health Hazards of Smoking” in 2012 indicated China had more than 350 million smokers, of which 14 millions were young smokers. The prevalence of smoking during adolescence was higher than other age stages. People during the period from adolescence to young adulthood are undergoing a series of continuous and important physiological development. Therefore, young people are relative sensitive to outside influences, especially the effects of nicotine exposure (Gulley and Juraska, 2013; Lydon et al., 2014).
During the resting state, the brain may be conceived as a complex functional network constituted by correlated spontaneous fluctuations occur within brain regions (Seeley et al., 2009). Electroencephalography (EEG) had been employed to measure ongoing oscillation of neural electrical activity in the absence of an overt task or stimulation (Delbeuck et al., 2003; Stam et al., 2007). By using spectral analysis of neural oscillations in EEG, previous resting state studies had found decreased theta and delta power, increased beta power and inconsistent results in alpha band.in smokers (Rass et al., 2015). Graph theoretical network analyses (GTA) of EEG signal during the resting state was used to investigate the integrity and efficiency of global brain networks in brain diseases, by assessing the topology characteristics of global brain networks (Delbeuck et al., 2003; Lee et al., 2006; Leistedt et al., 2009; Fraga González et al., 2016). The findings improved our understanding of neural mechanisms of brain diseases from whole brain network. However, few study used EEG to investigate the changes of global brain networks during the resting state in young smokers.
Spectral analysis provided us a way to detect the changes of spectral power in different neural oscillations, which were related with the cortical arousal or activation (Rass et al., 2015). Compared with traditional EEG methods, the phase lag index (PLI) assessed the consistency of phase lags between time series, which may avoid the effect of spurious interactions caused by common reference effects or volume conduction (Stam et al., 2007; Yuan et al., 2016). Moreover, the Minimum Spanning Tree (MST) analysis can be used to derive the connectivity matrix of functional network topology, which contains the strongest connections to avoid the problem of arbitrary thresholding of connectivity values for network reconstruction and preserve important information about network organization (Boersma et al., 2013; Tewarie et al., 2015). Therefore, in the current study, we constructed global network of EEG in different spectrum bands to investigate the changes of functional network topology during the resting state in young smokers. We used PLI to construct the matrix of functional connectivity and employed MST to reconstruct whole-brain functional connectivity networks of young smokers in different frequency bands (Brody et al., 2007; van Dellen et al., 2014, 2015; Tewarie et al., 2015; Fraga González et al., 2016). Main four EEG bands (delta, theta, alpha, and beta) with global brain recording sites were analyzed to examine the differences of global brain network topology between young smokers and nonsmokers. It is hoped that our study could improve the understanding of neural mechanisms in smoking.
Materials and Methods
Ethics Statement
The current study was consistent with the Declaration of Helsinki and was approved by the Medical Ethical Committee of the First Affiliated Hospital of Baotou Medical College, Inner Mongolia University of Science and Technology. All subjects wrote informed consent when they had understood the purpose of our study.
Participants
All participants were recruited from Inner Mongolia University of Science and Technology by advertisement. The exclusion criteria were any neurological illness, psychopathological diseases or epilepsy, drug abuse/dependence (excluding nicotine) based on DSM-V criteria and any current or past medications that may affect cognitive function, claustrophobia or not right-handed. Moreover, smokers were enrolled with regular smoked for more than 2 years and met the requirements of the DSM-IV criteria to nicotine dependence. The nicotine dependence was assessed by the Fagerström Test for Nicotine Dependence (FTND) (Fagerstrom and Schneider, 1989). Nonsmokers smoked less than 3 cigarettes in life. At last, twenty young male smokers (age: 20.95 ± 1.2 years) and 20 matched nonsmokers (age: 20.29 ± 1.1years) participated in our study. The smokers were asked to refrain from smoking about 60 min between the last cigarette and the acquisition of EEG data (average duration of abstinence before acquisition: 50.1 ± 5.8 min) (Yuan et al., 2016). The demographic characteristics of participants were presented in Table 1.
EEG Recording and Signal Processing
Prior to EEG recording, the expiratory carbon monoxide (CO) levels of participants were measured by using the Smokerlyzer system (Bedfont Scientific Ltd., Rochester, UK). The EEG data acquisitions were recorded using digital BrainAmp MR plus amplifiers (Brain Products GmbH. Munich. Germany) with Ag/AgCl electrodes at the 64 scalp sites according to the International 10/20 System. The horizontal electro-oculogram (HEOG) and vertical electro-oculogram (VEOG) recording were using two active Ag/AgCl electrodes placed at outer canthus of right eye and above the left eye. Electrical signal were digitized with a sample rate of 512 Hz and all electrodes impedance was kept below 5KΩ. Offline data was analyzed by Brain Vision Analyzer 2 (Brain Products Inc. Munich, Germany). The reference electrodes were the bilateral mastoids (TP9 and TP10). EEG data were filtered by using 0.5–70 Hz (IIR filter 24 dB/double frequency) band pass filter. In addition, independent component analysis (ICA) is used for the automatic correction function in Brain Visual Analyzer 2. The 2 min EEG data was segmented into 8s epochs (4,096 sample points per epoch) which were visually inspected for eye blinks or other artifacts. For each subject, 10 artifact-free epochs were chosen and exported to ASCII files.
After an initial stage of processing for 2 min EEG data, frequency bands were filtered as delta band (0.5 ~ 4 Hz), theta band (4 ~ 8 Hz), alpha band (8 ~ 13 Hz) and beta band (13 ~ 25Hz). The EEG data was entered into Brainwave soft (Brainwave 0.9.151.7.2 developed by C.S.; freely available at http://home.Kpn.nl/stam7883/brainwave.html).
Functional Connectivity
The functional connectivities of 64-channels EEG time series were analyzed by the phase lag index (PLI), a method of the asymmetry of the distribution of phase differences (Δϕ(tk) = 1,2,3,…,N) between two time-series (Stam et al., 2007). The PLI measures phase synchronization based on the asymmetry of distribution of instantaneous phase differences between two signals, which is determined by using the analytical signal based on the Hilbert transformation (Stam et al., 2007). Synchronization between two time series is calculated by the consistency of the phase lag of the non-zero to another signal (van Dellen et al., 2009). The instantaneous phase difference of each sampling point is analyzed by using the concept of signal analysis and Hilbert transform. The PLI have is high reliable of to measure asymmetric distribution measure functional interaction of these signal, which may be unaffected cannot make clear by volume conduction about a nonzero phase lag (van Dellen et al., 2014).
The PLI collected time series between two signals is used formula follow:
In here, || means absolute value, 〈〉 is use mean value and the PLI value range of 0 (no phase locking) to 1(complete synchronization) (Pischon et al., 2010).
Minimum Spanning Tree
The 64 × 64 weight adjacency matrix was obtained by PLI method. Then, MST method was used to analyze the weight matrix. The MST is not form a loop which is connected min edge M = N−1 and weight edges in the tree. In addition, using MST come into being a tree is will unique if weight graph is unique (Mare, 2008). To construct brain network defined node and edge is importance, so we defined 64 channels as nodes. The edges were detected by using MST method connection. The whole experiment processes were described in Figure 1.
Figure 1. An overview of the steps. First, artifact-free epochs are filtered for each frequency band (A). Secondly, the functional connectivity matrix based on Phase lag index (PLI) is calculated for each frequency band and epoch (B). Finally, using MST construct global networks and the tree view shows the hierarchical structure of the graph (C).
The following global metrics were calculated:
i) Degree (s): the tree node degree.
ii) Leaf fraction: number of tree nodes with exactly one connection to any other node.
iii) Diameter (d): the largest distance between any two nodes.
iv) Eccentricity: the shortest path distance between node i and any other tree node.
v) Betweenness centrality (BC): a network hub metric (Betzel et al., 2011) by where gjk is the shortest path between two tree nodes and gjk(i) is the number of those nodes that pass through node i.
vi) Hierarchy (TH): captures the ratio between a small diameter on one hand and overloading of hub nodes (Stam, 2014) by 2: .
vii) kappa (κ): captures the broadness of the degree distribution: (Stam and van Straaten, 2012; Stam et al., 2014).
Statistical Analysis
Two independent samples t-tests were used to examine the difference in the mean of two unrelated samples on the same variable by using SPSS 20.0 software. Pearson's correlation coefficients were calculated for performance of the Minimum spanning tree and smoking behaviors.
Results
Minimum Spanning Tree Network Analysis
Compared with nonsmokers, young smokers demonstrated decreased leaf fraction, kappa value, degree value and increased diameter, eccentricity value in alpha band (Table 2 and Figure 2). In detail, leaf fraction represented the whole tree distribution condition, which was decreased (t = −2.57, p = 0.009) in young smoker group. Diameter represented the small network short ways condition, which was increased (t = 2.387, p = 0.022) in young smoker. Eccentricity represented a tree leaf fraction and central correlation, which was increased (t = 2.371, p = 0.023) in young smoker. The kappa value responded to a tree high-degree and was decreased (t = −2.373, p = 0.023) in young smoker compare to nonsmoker. All changes between groups were shown in Figure 3. No significant changes in the delta, theta and beta band were found between the two groups.
Figure 2. MST matrices and MST graph in scalp view and tree view for the alpha band for nonsmoker and smoker. For illustrative the MST method was difference brain networks in smoker compare with nonsmoker.
Figure 3. Group averages for leaf fraction, eccentricity, diameter and kappa measures of the MST. Open bars refer to nonsmoker and filled bars to young smoker (*p < 0.05).
Correlation between Brain and Smoking Behaviors
Pearson's correlation coefficients were calculated for alpha band (Eccentricity, Diameter, Leaf fraction, kappa) and smoking behaviors (pack-years, FTND, onset time (age), CPD). Significant correlation was found between leaf fraction and onset times (Figure 4, r = 0.574, p = 0.008) in smokers. No other significant correlations were found.
Figure 4. Significant correlation was found between leaf fraction and onset times in young smoker (r = 0.574, p = 0.008).
Discussion
Complex brain networks may play a key role in the integrated processing of the brain function, which may be influenced by pathological processing. In present study, we analyzed EEG signals during the resting state to examine the topological characteristics of brain networks in young smokers by using PLI and MST methods. The major problem of EEG was volume conduction, which may cause overestimation of functional connectivity. The PLI measured the consistency of phase lags between EEG time series and provide more accurate results of functional connectivity (Stam et al., 2007; van Dellen et al., 2015; Yuan et al., 2016). Moreover, the MST method assessed the strongest connections in the network without forming loops, which solved the problem of arbitrary thresholding of connectivity values and preserved important information about network organization (Tewarie et al., 2015). Compared with nonsmokers, young smokers showed decreased leaf fraction, kappa value and increased diameter, eccentricity value in alpha band.
Leaf fraction indicates the extent of central organization in global network. High leaf fraction suggests that the communication of network is largely dependent on hub nodes. Kappa value indicates the broadness of the degree distribution (degree divergence). Diameter is the largest distance between any two nodes of the tree, which is related to the efficiency of global brain network organization. Information may be efficiently processed between remote brain regions with low diameter. In addition, eccentricity indicates longest distance between a reference node and any other node, which reflects the efficient information processing from the least central node (van Dellen et al., 2015; Fraga González et al., 2016). In present study, decreased leaf fraction, kappa value were found in young smoker, which were similar with previous EEG studies brain neural oscillations properties in smokers during the resting state (Rass et al., 2015). Decreased alpha band power in smokers might be associated with long-term consequences of nicotine use (Rass et al., 2015; Wilbanks et al., 2016). Despite of the significant advantages, the PLI analyses failed to reveal differences of the connectivity strength between nodes of the brain regions. The MST method avoided potential bias as thresholding value (Stam, 2014), which constructed of brain connectivity network (van Diessen et al., 2014). The MST analysis had been successfully applied to EEG data from several diseases (Delbeuck et al., 2003; van Dellen et al., 2014; Fraga González et al., 2016). In the current report, we found that smokers were associated with altered global brain networks. We also revealed the abnormal connectivity strength in the alpha band in young smokers by using the MST analyses. The optimal tree was configured by both short distances and no overload of the central tree nodes (Stam, 2014). In our study, the MST analysis presented us a disordered network of smokers as marked by decreases in leaf number indicating low network integration compare with nonsmokers. These results suggested a more path-like configuration in young smoker and a more star-like topology in nonsmoker in the tree topologies. The MST leaf fraction was lower for trees derived from regular networks and increased as these networks became more random and increased distances (Tewarie et al., 2015). Increased diameter value emerge to less communication between nodes of the network in smokers, which suggested the brain networks became more disorganized and decreased global efficiency in smokers. The eccentricity and kappa descript the configuration features of tree. In present study, increased eccentricity may reveal increased longest distance between a reference node and any other node in young male smokers, which reflects decreased efficient information processing from the least central node. Decreased kappa value was related to decreased broadness of the degree distribution, which may suggest slow synchronization of the tree in young male smokers (Stam and van Straaten, 2012).
Our results indicated less integrated network configuration in young smokers, which confirmed the decrease of global efficiency in alpha band in young smokers compared with nonsmokers. These findings may be associated with cholinergic deficits and cognitive impairment in young smokers. Positive correlation was found between leaf fraction and onset times in young smoker, which may suggest initial smoking was associated with the extent of central organization in global network. Young male smokers were associated with changes in brain network topology, including further decreased leaf number and increased eccentricity in the brain.
Limitation
The study is limited in its cross-sectional design, and it is not possible to infer the causality of the relationships between alpha band and severity of smoking. The relationships might be bi-directional and related to some latent variables. In addition, our participants were all male. Previous studies revealed clear gender effect in smoking effects (Kühn et al., 2012; Franklin et al., 2014). Another limitation of the study is the relatively small sample size, so the results might not be generalize able to wider populations. Further longitudinal studies with larger sample size would be required to improve the understanding of the EEG findings in the current study.
Conclusion
The abnormal topological characteristics of functional brain networks in alpha band were found in young smokers during the resting state. The findings may provide novel insights into the neural mechanism in smoking. We hoped our study could improve the understanding of smoking.
Author Contributions
SS, DY, KY, and XL conceived and designed the experiments; KY, SS, JC, YC, XZ, and YG performed the experiments; DY, KY, JC, YC, and TX analyzed the data; SS, DY, and KY wrote the article; YL and YB provided critical revision of the manuscript for important intellectual content. All authors critically reviewed content and approved final version for publication.
Funding
This paper is supported by the National Natural Science Foundation of China under Grant Nos. 81571753, 81571751, 61771266, 81701780, 81401488, 61179019, 81401478, 81271644, 81271546, 81271549, 81470816, 81471737, the Fundamental Research Funds for the Central Universities under the Grant Nos. JB151204, JB121405, the program for Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region NJYT-17-B11, the Natural Science Foundation of Inner Mongolia under Grant No. 2014BS0610, 2017MS(LH)0814, the program of Science and Technology in Universities of Inner Mongolia Autonomous Region NJZY17262, the Innovation Fund Project of Inner Mongolia University of Science and Technology Nos. 2015QNGG03, 2014QDL002, General Financial Grant the China Post- doctoral Science Foundation under Grant No. 2014M552416.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
Betzel, R. F., Erickson, M. A., Abell, M., O'donnell, B. F., Hetrick, W. P., and Sporns, O. (2011). Synchronization dynamics and evidence for a repertoire of network states in resting EEG. Front. Comput. Neurosci. 6:74. doi: 10.3389/fncom.2012.00074
Boersma, M., Smit, D. J., Boomsma, D. I., De Geus, E. J., Delemarre-van de Waal, H. A., and Stam, C. J. (2013). Growing trees in child brains: graph theoretical analysis of electroencephalography-derived minimum spanning tree in 5- and 7-year-old children reflects brain maturation. Brain Connect. 3, 50–60. doi: 10.1089/brain.2012.0106
Brody, A. L., Mandelkern, M. A., Olmstead, R. E., Jou, J., Tiongson, E., Allen, V., et al. (2007). Neural substrates of resisting craving during cigarette cue exposure. Biol. Psychiatry 62, 642–651. doi: 10.1016/j.biopsych.2006.10.026
Delbeuck, X., Linden, M. V. D., and Collette, F. (2003). Alzheimer' disease as a disconnection syndrome? Neuropsychol. Rev. 13, 79–92. doi: 10.1023/A:1023832305702
Fagerstrom, K.-O., and Schneider, N. G. (1989). Measuring nicotine dependence: a review of the Fagerstrom tolerance questionnaire. J. Behav. Med. 12, 159–182. doi: 10.1007/BF00846549
Fraga González, G., Van der Molen, M. J., Žarić, G., Bonte, M., Tijms, J., Blomert, L., et al. (2016). Graph analysis of EEG resting state functional networks in dyslexic readers. Clin. Neurophysiol. 127, 3165–3175. doi: 10.1016/j.clinph.2016.06.023
Franklin, T. R., Wetherill, R. R., Jagannathan, K., Johnson, B., Mumma, J., Hager, N., et al. (2014). The effects of chronic cigarette smoking on gray matter volume: influence of sex. PLoS ONE 9:e104102. doi: 10.1371/journal.pone.0104102
Gulley, J. M., and Juraska, J. M. (2013). The effects of abused drugs on adolescent development of corticolimbic circuitry and behavior. Neuroscience 249, 3–20. doi: 10.1016/j.neuroscience.2013.05.026
Kühn, S., Romanowski, A., Schilling, C., Mobascher, A., Warbrick, T., Winterer, G., et al. (2012). Brain grey matter deficits in smokers: focus on the cerebellum. Brain Struct. Funct. 217, 517–522. doi: 10.1007/s00429-011-0346-5
Lee, U., Kim, S., and Jung, K. Y. (2006). Classification of epilepsy types through global network analysis of scalp electroencephalograms. Phys. Rev. E 73, 89–132. doi: 10.1103/PhysRevE.73.041920
Leistedt, S. J., Coumans, N., Dumont, M., Lanquart, J. P., Stam, C. J., and Linkowski, P. (2009). Altered sleep brain functional connectivity in acutely depressed patients. Hum. Brain Mapp. 30, 2207–2219. doi: 10.1002/hbm.20662
Lydon, D. M., Wilson, S. J., Child, A., and Geier, C. F. (2014). Adolescent brain maturation and smoking: what we know and where we're headed. Neurosci. Biobehav. Rev. 45, 323–342. doi: 10.1016/j.neubiorev.2014.07.003
Mare, M. (2008). The saga of minimum spanning trees. Comput. Sci. Rev. 2, 165–221. doi: 10.1016/j.cosrev.2008.10.002
World Health Organization (2013). WHO Report on the Global Tobacco Epidemic. Enforcing Bans on Tobacco Advertising, Promotion and Sponsorship. Geneva: World Health Organisation.
Pischon, T., Boeing, H., Hoffmann, K., Doll, R., Peto, R., Boreham, J., et al. (2010). Effects of obesity and smoking on US life expectancy. N. Engl. J. Med. 2010, 855–857. doi: 10.1056/NEJMc1000079
Rass, O., Ahn, W. Y., and O'donnell, B. F. (2015). Resting-state EEG, impulsiveness, and personality in daily and nondaily smokers. Clin. Neurophysiol. 127, 409–418. doi: 10.1016/j.clinph.2015.05.007
Seeley, W. W., Crawford, R. K., Zhou, J., Miller, B. L., and Greicius, M. D. (2009). Neurodegenerative diseases target large-scale human brain networks. Neuron 62, 42–52. doi: 10.1016/j.neuron.2009.03.024
Stam, C. J. (2014). Modern network science of neurological disorders. Nat. Rev. Neurosci. 15, 683–695. doi: 10.1038/nrn3801
Stam, C. J., Nolte, G., and Daffertshofer, A. (2007). Phase lag index: assessment of functional connectivity from multi channel EEG and MEG with diminished bias from common sources. Hum. Brain Mapp. 28, 1178–1193. doi: 10.1002/hbm.20346
Stam, C. J., Tewarie, P., Van Dellen, E., Straaten, E. C. W. V., Hillebrand, A., and Van Mieghem, P. (2014). The trees and the forest: characterization of complex brain networks with minimum spanning trees. Int. J. Psychophysiol.92, 129–138. doi: 10.1016/j.ijpsycho.2014.04.001
Stam, C. J., and van Straaten, E. C. (2012). The organization of physiological brain networks. Clin. Neurophysiol. 123, 1067–1087. doi: 10.1016/j.clinph.2012.01.011
Tewarie, P., van Dellen, E., Hillebrand, A., and Stam, C. J. (2015). The minimum spanning tree: an unbiased method for brain network analysis. Neuroimage 104, 177–188. doi: 10.1016/j.neuroimage.2014.10.015
van Diessen, E., Otte, W. M., Braun, K. P., Stam, C. J., and Jansen, F. E. (2014). Does sleep deprivation alter functional EEG networks in children with focal epilepsy? Front. Syst. Neurosci. 8:67. doi: 10.3389/fnsys.2014.00067
van Dellen, E., de Waal, H., van der Flier, W. M., Lemstra, A. W., Slooter, A. J., Smits, L. L., et al. (2015). Loss of EEG network efficiency is related to cognitive impairment in dementia with lewy bodies. Mov. Disord. 30, 1785–1793. doi: 10.1002/mds.26309
van Dellen, E., Douw, L., Baayen, J. C., Heimans, J. J., Ponten, S. C., Vandertop, W. P., et al. (2009). Long-term effects of temporal lobe epilepsy on local neural networks: a graph theoretical analysis of corticography recordings. PLoS ONE 4:e8081. doi: 10.1371/journal.pone.0008081
van Dellen, E., Douw, L., Hillebrand, A., Hamer, P. C. D. W., Baayen, J. C., Heimans, J. J., et al. (2014). Epilepsy surgery outcome and functional network alterations in longitudinal MEG: a minimum spanning tree analysis. Neuroimage 86, 354–363. doi: 10.1016/j.neuroimage.2013.10.010
Wilbanks, H. E., Van Mohr, M., Potenza, M. N., Mayes, L. C., and Rutherford, H. J. (2016). Tobacco smoking and the resting maternal brain: a preliminary study of frontal EEG. Yale J. Biol. Med. 89, 115–122.
Keywords: smoking, electroencephalogram (EEG), resting state, minimum spanning tree (MST), leaf fraction
Citation: Su S, Yu D, Cheng J, Chen Y, Zhang X, Guan Y, Li Y, Bi Y, Xue T, Lu X and Yuan K (2017) Decreased Global Network Efficiency in Young Male Smoker: An EEG Study during the Resting State. Front. Psychol. 8:1605. doi: 10.3389/fpsyg.2017.01605
Received: 15 May 2017; Accepted: 01 September 2017;
Published: 15 September 2017.
Edited by:
Qinghua He, Southwest University, ChinaReviewed by:
Jiaojian Wang, University of Electronic Science and Technology of China, ChinaLinling Li, Shenzhen University, China
Copyright © 2017 Su, Yu, Cheng, Chen, Zhang, Guan, Li, Bi, Xue, Lu and Yuan. 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) or licensor 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: Dahua Yu, fmydh@imust.edu.cn
Kai Yuan, kyuan@xidian.edu.cn
†These authors have contributed equally to this work.