Skip to main content

Disrupted brain network dynamics and cognitive functions in methamphetamine use disorder: insights from EEG microstates



Dysfunction in brain network dynamics has been found to correlate with many psychiatric disorders. However, there is limited research regarding resting electroencephalogram (EEG) brain network and its association with cognitive process for patients with methamphetamine use disorder (MUD). This study aimed at using EEG microstate analysis to determine whether brain network dynamics in patients with MUD differ from those of healthy controls (HC).


A total of 55 MUD patients and 27 matched healthy controls were included for analysis. The resting brain activity was recorded by 64-channel electroencephalography. EEG microstate parameters and intracerebral current sources of each EEG microstate were compared between the two groups. Generalized linear regression model was used to explore the correlation between significant microstates with drug history and cognitive functions.


MUD patients showed lower mean durations of the microstate classes A and B, and a higher global explained variance of the microstate class C. Besides, MUD patients presented with different current density power in microstates A, B, and C relative to the HC. The generalized linear model showed that MA use frequency is negatively correlated with the MMD of class A. Further, the generalized linear model showed that MA use frequency, scores of Two-back task, and the error rate of MA word are correlated with the MMD and GEV of class B, respectively.


Intracranial current source densities of resting EEG microstates are disrupted in MUD patients, hence causing temporal changes in microstate topographies, which are correlated with attention bias and history of drug use.

Peer Review reports


Drug dependence is a severe medical and public health challenge globally. Methamphetamine (MA), one of the amphetamine-type stimulants, induces toxic effects on the brain, such as oxidative stress, mitochondrial toxicity, excitotoxicity and neuroinflammation [1]. Overall, these effects disrupt the structure and functioning of the brain. Previous studies reported that abuse of MA is associated with extensive cognitive impairment, encompassing the executive function, working memory, problem-solving, and attention deficit [2, 3]. It has been shown that these impairments are correlated to brain imaging changes in MA-dependent patients [2, 4]. Some of the brain regions damaged by MA include the dorsal prefrontal cortex, medial prefrontal cortex, anterior cingulate cortex, striatum, and insula. Notably, magnetic resonance imaging (MRI) can effectively diagnose and evaluate cognitive impairment in MA dependent patients. Besides, the typical characteristics of drug use history are also closely associated with cognitive impairment and structural changes of the brain of MA dependent patients [5,6,7].

Information processing by the brain is associated with the temporal and spatial dynamics of the whole brain resting network [8, 9]. Coordination and integration of multiple brain regions is an essential basis of understanding cognitive function [10, 11]. Studies have indicated that the spatial-temporal characteristics of brain networks are beneficial in defining and understanding the elements of cognitive impairment [12,13,14]. However, most previous studies on cognitive function in MA-dependent patients have focused on the behavioural changes and alterations in the spatial structure of the brain. Meanwhile, only a few have explored the temporal characteristics of cognitive processing [7, 15, 16]. Therefore, there is a dearth of information on cognitive impairment in patients with substance use disorder, such as methamphetamine use disorder (MUD).

In recent years, EEG microstate methods are increasingly being used to study the spatial and temporal resting state dynamics of brain networks [17]. Resting state EEG displays as periods of potential topographies and shows quasi-stable configurations called functional microstates [18]. Functionally, EEG microstates represent the momentary local states and may indicate spontaneous fluctuations in activity in large scale brain networks [19]. Also, EEG microstates reflect basic brain activity in spontaneous and event-related EEG researches [20, 21]. Most studies have clustered EEG microstate topographies into four standard microstate maps, and these categories represent significant brain configuration networks [17]. The various microstates are associated with different types of cognitive processes [22,23,24,25,26], suggesting that microstates might effectively analyze cognitive deficits resulting from a mental disorder and addiction. Britz and colleagues used simultaneous EEG-fMRI techniques to explore the correlation between microstates and distinct fMRI resting-state networks [22]. Previous research showed that the four microstates marked A, B, C, and D correspond to resting-state networks that had previously been identified to be associated with phonological processing, the optical network, the saliency network, and attention, respectively [22]. Besides, these four microstate classes are correlated with the various aspects of the default mode network [27], which play an important cognitive role in psychiatric diseases [28]. Notably, studies focusing on the association between resting-state EEG and cognition impairment in MUD patients mainly use the frequency and entropy analysis methods. The results of these studies suggested an important relationship between EEG signals and cognitive function [29,30,31,32]. The micro-state analysis method is beneficial in extracting the temporal information of EEG signals.

The findings of previous studies on several psychiatric and neurological disorders suggest that aberrant EEG microstates reflect abnormal information processing capabilities during the cognitive process [17, 21, 33,34,35]. However, the temporal dynamics of the brain networks in drug users have not been analyzed using the resting EEG microstate method. The present study, therefore, used this method to study the association between disrupted brain dynamic networks and cognitive function of MA users and the potential effect of drug use history. Through this, it is possible to provide novel strategies for understanding the dynamics of cognitive function disruption in MA-dependent patients and other drug users.

Therefore, this study aimed at investigating the resting state EEG microstates in MA-dependent patients. Since structural alterations have been identified using other brain imaging methods (i.e. MRI and PET), we hypothesized that relative to controls, a disrupted resting state EEG microstate may develop in MA-dependent patients. Besides, these disruptions could be associated with cognitive impairment and the characteristics of drug use history.



This study is a part of a larger project investigating the objective marker system for diagnosis of MUD and relapse warning among MA users in China. We recruited 57 participants with MUD from two compulsory rehabilitation centers in Shanghai and Yunnan provinces, China. All patients in the rehabilitation centers from June 2017 to August 2018 were invited to participate in the study. The recruitment criteria included: (1) minimum age of 18 years; (2) met the requirements of the fifth edition of the Diagnostic and Statistical Manual of Mental Disorders (DSM-5) for MUD; (3) used MA in the three months before recruitment in the study (4) normal vision and hearing, or within the normal range after correction (5) right-handed. The exclusion criteria were (1) diagnosis of other substance use disorders or other axes I psychiatric disorders; (2) history of head injury, neurological disorders, or loss of consciousness; (3) patients on any antidepressant or/and antipsychotic medication.

An additional 30 healthy controls (HC) were recruited from the local community through advertisement and word-of-mouth. The recruitment was conducted between June 2017 and November 2018, after which HC were matched with the MA groups by gender and age. Regarding the choice of matching variables, previous studies have shown that age significantly affects resting EEG signals and cognitive function [36]. Besides, previous studies have suggested that EEG signals and cognitive functions differ between males and females [37], so both age and gender were used as matching variables in this study. The enrollment of HC was done simultaneously with that of patients, so the matching process was part of recruitment. For HC, the exclusion criteria included (1) current substance abuse or dependence; (2) history of substance dependence other than nicotine; (3) a current or previous mental illness (or/and treatment); (4) history of head injury, neurological disorders, or loss of consciousness; (5) without normal vision or hearing, or out of range after correction; (6) left-handed. The study protocol was implemented in compliance with the principles of the Declaration of Helsinki and approved by the Ethical Committee of the Shanghai Mental Health Center. Informed consent was signed by all subjects before participation in this research.


Socio-demographic variables

After obtaining informed consent, trained research assistants at the Shanghai Mental Health Center used a Chinese version of the Addiction Severity Index (ASI-C) to collect baseline information from the subjects. The information was collected via face to face interviews in a private environment to avoid disturbance from others. Patients in the baseline study were interviewed after detoxification. The variables used included basic demographic characteristics (age, gender), history of MA use (accumulated years of MA use, frequency, and quantities of daily use), and alcohol use history. The psychometric properties of ASI-C have been validated in previous studies [38].

CogState battery

CogState Battery is a computerized cognitive tool designed to assess cognitive function, and consists of Two Back Task (TWOB, working memory), International Shopping List Task (ISL, verbal learning and memory), the Groton Maze Learning Task (GML, problem solving/ error monitoring), Social Emotional Cognition Task (SEC, social cognition), and Continuous Paired Association Learning Task (CPAL, spatial working memory). Detailed procedures and rules of the above paradigms have been provided in past research [39]. In these five tasks, the ISL score is the total correct response, the TWOB and SEC scores are the rates of accurate answers, while the GML and CPAL sores are the total error numbers. Good reliability and validity of the CogState Battery (Chinese version) were verified previously [39].

Stroop task

Attention bias was assessed by the MA Addiction Stroop Task [40]. Two target words ((“ICE” means “MA” and “SKATING” is the nickname for how Chinese drug users refer to taking MA)) and two neutral words ((“MESSY GRASS” and “RESOLVE”)) were involved in this paradigm.

The subjects were asked to press the buttons according to the color of the words presented (red, green, yellow, and blue), and ignore the meaning of the words. Each of the four words was presented eight times and was shown on the screen for 3000 ms each time. The word appeared randomly during the test, and the same word was set not to appear three times consecutively. This task included 64 MA-related word trials and 64 neutral word trials. Two outcomes were measured for this task: (1) Reaction time, which was obtained from button presses of the dominant hand on four color-marked buttons; (2) Errors, calculated as a ratio of the number of wrong color pressed by the subjects to the total number of trials recorded. More detailed information about this paradigm can be found in previous articles published by colleagues in our laboratory [40]. All assessments were done following the standardized instructions handbook by a trained researcher. The fixation cross and words were displayed on black background 75 cm from the eyes.

EEG recording and processing

The scalp EEG data were recorded from a high-density 64-channel sintered Ag/AgCl electrode cap (BrainCap; Asiacut, Germany). Two vertical electrooculography (EOG) electrodes were positioned above and below the right eye, and the other two horizontal EOG electrodes were placed on the outer canthi of both eyes. The reference electrode was located at the tip of the nose, while the ground electrode was placed on the forehead. The EEG was sampled continuously at a sampling frequency of 1000 Hz, and all impedances were kept below 5 kΩ. EEG was recorded in the eyes-closed resting state for 5–6 min.

EEG data were bandpass filtered at 0.1–70 Hz and eye movement artifacts were removed using independent component analysis. Manual screening for other artifacts was then conducted. In accordance with the steps detailed in previous research [18, 19], data were subsequently segmented into artifact-free 2 s epochs before microstate analysis, digitally band passed at 2–20 Hz and down-sampled to 250 Hz. The first 20 artifact-free epochs of each study participant were after re-referenced to the average reference and used for the analysis of resting EEG microstates. EEG data were processed using EEGLAB 13.0b [41] in MATLAB2013b (The MathWorks, Inc., Natick, US).

Microstate analysis

EEG epochs were transformed into global scalp sequences of momentary potential map distributions to analyze EEG microstates. We used the standard strategies for analysis of microstates described in previous researches [42,43,44]. The fitted microstate class topographies were computed using a modified version of the k-means spatial cluster analysis based on previous studies [19, 45]. Original topographies were clustered into four microstate classes according to the findings of a large normative study in healthy populations [19]. We then performed a two-steps k-means clustering analysis. At the individual level (first-step), global field power peaks of each participant were submitted for k-means clustering analysis, and four microstate classes were obtained per individual, which is consistent with previous studies [46,47,48]. In the second step, also known as the group-level, the dominant topographies across all individuals and groups were clustered to obtain four mean microstate classes for each group. Fitting of microstates was then conducted to calculate the related parameters, including mean microstate duration (MMD), the ratio of total time covered (RTT), global explained variance (GEV), and occurrence of each microstate in each participant. The MMD indicates the mean times of each microstate class, RTT indicates the mean ratio of total time covered by each of the four microstates, while GEV quantifies the extent to which each of the four microstates describes the data. Microstate analysis was performed using the free academic software Cartool ( site/cartoolcommunity/) [49].

Statistical analysis

The Student’s T-test or the Chi-square test was used to test for differences in baseline characteristics between patients with MUD and HC. Data with non-homogenous variance and/ or not normally distributed were analyzed using non-parametric tests. Two-way Analysis of Variance, performing with microstate class as within-subject factors and the subject group as between factor, was used to compare each microstate-related measurement. Two-tailed Student’s T-test was then used to compare significant results from each group. All microstate related parameters of each microstate were tested four times. The corrected T-test alpha value was set to 0.0125 based on the Bonferroni correction method, and the generalized linear regression model was used to explore the correlation between significant microstate related parameters and other factors (i.e. drug history, CogState Battery and Addiction Stroop Task).

We then investigated the intracerebral current sources of each microstate classes which best differentiated the groups in microstate-related parameters using standardized low-resolution brain electromagnetic tomography (sLORETA) [50]. The intracerebral current sources were calculated from the individual microstate templates of each class for all MA-dependent patients and HC. Electrode coordinates were transferred to 6239 cubic voxels corresponding to cortical grey matter and hippocampi in the sLORETA inverse solution space. Differences between the two groups were assessed by voxel-by-voxel unpaired sample log of F-ratio-tests for each microstate class. The statistical non-parametric mapping randomization test was used to correct critical probability threshold values for multiple comparisons. A total of 5000 permutations were used to determine the significance of each randomization test. Source analysis was performed using the sLORETA software (


Description of member characteristics

Five individuals (including 2 MA-dependent users and 3 HC) who had more than five bad channels were excluded in the subsequent analysis. To assess the sample size, we performed a power analysis using the microstates MMD difference between the two groups. The analysis results revealed that our sample size has a power of 0.82, suggesting reliability of the subsequent results. Table 1 summarizes the demographic characteristics and performance of each cognitive task (Cogstate Battery and Stroop Task) for both MA and HC samples. There were no significant differences in gender (t = 0.00, p = 0.96) and mean age (t = 0.53, p = 0.59) between MA-dependent patients and HC.

Table 1 Baseline characteristics of participants

Resting state EEG microstates

The topography classes observed in this study were similar to four classes reported in the previous literatures, and were labeled microstate A, B, C, and D accordingly (Fig. 1) [18, 19, 51]. The four microstate topographies included one with a right anterior-left posterior orientation (microstate A), one with a left anterior- right posterior orientation (microstate B), one with a central anterior-posterior orientation (microstate C), and one with a fronto-central extreme location (microstate D).

Fig. 1
figure 1

Topographical microstate maps. The figure shows spatial configuration of the four microstate classes (ad). Upper row: microstate configurations of healthy controls. Lower rows: microstate configurations of patients with methamphetamine use disorder. Different colors represent the different polarities. The maps are presented as seen from above angle

The microstate parameters of the four microstate classes are summarized in Table 2. No significant group effects were identified in the four microstate parameters, but a significant interaction between the groups and microstate classes were observed for GEV (F = 3.99, p = 0.01), MMD (F = 2.67, p = 0.05), and RTT (F = 3.31, p = 0.02). Subsequent analysis indicated significantly lower MMD (t = − 2.77, corrected p = 0.03) in microstate class A of MA patients relative to HC. For microstate B, we found a significantly lower MMD (t = − 3.30, corrected p = 0.01) and GEV (t = − 3.19, corrected p = 0.01) in MA patients relative to HC. In contrast, microstate class C presented with significantly higher GEV (t = 2.61, corrected p = 0.04) in MA patients compared to HC.

Table 2 Microstate parameter results summary

EEG microstate source differences

Analysis of current density power by sLORETA showed significant differences in microstates A, B, and C between MA patients and HC. For microstate A, suprathreshold cortical voxels (Threshold: Log-F = 0.70, p = 0.05) displayed increased activation of current density power in the inferior parietal lobule, precentral gyrus, middle frontal gyrus, middle temporal gyrus, superior temporal gyrus, superior parietal lobule, supramarginal gyrus, and insula (Table S1 and Fig. S1 in the Supplement). For microstate B, suprathreshold cortical voxels (Threshold: Log-F = 0.76, p = 0.05) displayed increased activation in the regions of the inferior parietal lobule, middle occipital gyrus, inferior frontal gyrus, cuneus, precuneus, superior occipital gyrus, supramarginal gyrus, and angular gyrus among others (Table S2 and Fig. S2 in the Supplement). For microstate C, increased activations (Threshold: Log-F = 0.80, p = 0.05) were observed within the inferior frontal gyrus, inferior temporal gyrus, cingulate gyrus, posterior cingulate, and the parahippocampal gyrus among others, as shown in Table S3 and Fig. S3.

Relationships between microstate parameters and other factors

A generalized linear model of class A indicated a negative correlation between MA use frequency (β = − 0.18, 95%CI = -0.74 to − 0.06, p = 0.02) and MMD. In Class B, a generalized linear model also indicated that MA use frequency (β = − 0.57, 95%CI = -0.96 to − 0.19, p = 0.00), scores of Two-back task (β = 0.45, 95%CI = 0.12 to 0.77, p = 0.01), an error rate of MA word (β = − 5.35, 95%CI = -8.52 to − 2.17, p = 0.00) were correlated with MMD. In class B also, we found that the MA use frequency (β = − 0.11, 95%CI = -0.17 to − 0.06, p = 0.00), scores of Two-back task (β = 0.06, 95%CI = 0.01 to 0.11, p = 0.01), and the error rate of MA word (β = − 0.58, 95%CI = -1.05 to − 0.10, p = 0.02) were correlated with GEV. None of other variables, including drug history, CogState Battery, and Addiction Stroop Task were found correlated with the MMD of class A, MMD of class B, or GEV of class B.

Also, no correlation was observed between the GEV of class C and drug history, CogState Battery, and Addiction Stroop Task.


To the best of our knowledge, this is the first study to investigate the resting state EEG microstate abnormality in drug abuse patients. Also, the present study is the first to explore the correlation between the alterations in brain resting network and cognitive functions (i.e. CogState Battery and Addiction Stroop Task), along with drug use history (i.e. accumulated years of MA use, frequency of MA use, daily dosage, and severity of MA use disorder). Here, we found that the temporal microstate profiles of MA-dependent patients are different from those of HC, and these microstate deviations could be explained by changes in intracranial current source densities. Moreover, our findings suggest that MA use frequency is associated with changes in microstates A and B, while variations in Addiction Stroop Task are correlated with changes in microstate B.


Microstates reflect the dynamic synchronization of neuronal functional networks. Relative to HC, we found that MA-dependent patients have altered microstates A, B, and C, but microstate D is comparable between the two groups. Although results on brain network studies of MA users are difficult to compare directly due to the different techniques and analytical methods adopted, the findings of several previous studies are consistent with our results. As Ahmadlou and his colleagues reported, chronic MA abusers exhibited disrupted Small-World brain network and global network deficit at the gamma band [52]. Functional Magnetic Resonance Imaging also revealed abnormal functional connectivity, which correlates with psychotic symptoms of MA patients [53]. Therefore, our results on brain networks partially concur with those of previous studies.

The various microstate parameters provide information about the underlying neural generators [21]. For instance, MMD is indicative of the stability of a microstate and its underlying neural assemblies [54, 55]. Relative to HC, we found that the MMD of MA-dependent patients is significantly shorter for microstates A and B. In contrast, microstates C and D were not significantly different between the two groups. This result suggests that MA-dependent patients have underlying instability in microstates A and B, but are more stable in microstates C and D. Similarly, patients with schizophrenia have reduced MMD in microstates A and B [56, 57], suggesting an electrophysiological similarity in the mechanism of schizophrenia and MA-induced psychosis [58,59,60]. GEV reflects the time covered in each microstate class and its underlying neural generators [21]. Relative to HC, we found that the GEV of patients with MA use disorder exhibited higher in microstate C, but lower for microstate B. Notably, patients with obsessive-compulsive disorder also have abnormalities in these two microstate classes [61]; that is, the occurrence of these microstates was significantly different from those of HC. Previous studies have also suggested similarities in the manifestation of visual and salience network networks in patients with obsessive-compulsive disorder and addiction [62,63,64,65,66,67]. Therefore, abnormalities in these brain networks may partly explain the alteration in the parameters in microstates B and C in these two populations [22].

Spatial configuration of microstate classes and the relationship between parameters

The findings of researches that combined resting state EEG and fMRI suggested that EEG microstate classes mirror various large-scale resting-state networks, and thus reflected the specific functional state of the mental process [22, 68]. For instance, Microstate A, which is closely associated with blood oxygen level-dependent (BOLD) activation of bilateral superior temporal and parietal cortex, was suggested to reflect the sensorimotor or auditory network [22, 69]. We conducted the sLORETA analysis to interpret the changes in microstate class A in MA-dependent patients and observed hyper-activation of the superior temporal and parietal area in these patients, relative to HC. Besides, we found that the microstate MMD of MA dependent patients was shorter than that of HC. Interestingly, we found that the MMD of microstate A was negatively correlated with the MA use frequency. Previous studies have reported that MA abuse impairs sensorimotor cortical plasticity and causes deficits in behavioural performance [70, 71]. The frequency of MA use was found to be associated with the volume of cortical gray matter [72]. Daumann and colleague suggested that lifetime use of MA is negatively correlated with individual gray matter [73]. Ruan and colleague also indicated that accumulated years of MA use are negatively correlated with the thickness of gray matter in the right superior temporal cortex [5]. These findings indicate that the aberrance of parameters of microstate A may reliably indicate the severity of MA use.

As one of the specific predictors of relapse, attention bias was assessed by the MA Addiction Stroop Task. We found that attention bias is correlated with aberrant temporal functioning of microstate B. The microstate B was found associated with negative BOLD activation in the striate and extrastriate cortex [22]. Similar to our findings, a previous study reported that the MMD of microstate B is significantly lower in patients with mild spastic diplegia [74]. Besides, patients with moderate spastic diplegia display various visual categorization impairments and attention deficits [75, 76]. Collectively, microstate class B may play an essential role in visual perception and attention. Therefore, the risk of attention bias in MA-dependent patients can be assessed by measurement of microstate B. Besides, the present study also demonstrated that MA use frequency is correlated with changes in the GEV and MMD of microstate B. Previous studies have shown that drug use history is associated with deficits in brain structure and cognitive impairment in MA-dependent patients [5, 16]. sLORETA analysis of microstate B further suggested that regions related to visual processing (including middle occipital gyrus, cuneus, superior occipital gyrus, and the precuneus) are hyperactive in patients with MA use disorder. Therefore, it is reasonable to conclude that the microstate B may be an explanation of cue-reactivity in MA-dependent patients. Interestingly, we also identified a positive correlation between working memory performance and the GEV of microstate B. It has been suggested that the parahippocampal gyrus and fusiform gyrus (specific brain area of visual network) transmit information relevant to the visual cue. Notably, their spatial activation patterns during selective maintenance of the cue-related picture type match those during the processing of visual working memory [77]. Therefore, abnormalities in microstate B also partially explain the severity of working memory deficits in patients with MA use disorder. Nevertheless, aberrant levels of microstate parameters do not indicate variations in the function of the corresponding brain region [57, 78]. Whether the behavioral performance is regulated by the balance between different EEG microstates and whether the performance is related to a single microstate parameter within a reasonable range are subjects of future research.

Limitations and strengths

There are some limitations to this study. Although the study participants were recruited from two regions, our sample size was relatively small. Therefore, future studies should use a larger sample size to increase the statistical power and the robustness of the research, and hence generate more reliable results for clinical applications. Second, although we screened all the possible patients in the two rehabilitation centers, our study population may not represent the characteristics of the MA patients at all various treatment states and different abstinence periods. Another shortcoming is that we did not simultaneously record MRI data, which provides a more direct picture of the activities in the brain regions. Therefore, we could not determine the correlation between EEG microstates and the activity of distinct brain areas. In addition, due to the cross-sectional nature of our study, the conclusion that the aberrance observed in the microstates resulted from MA dependence should be made with more caution. As such, a longitudinal study should be done to verify our findings. The strengths of this study are as follows. Firstly, the gender and age of our patient group and the HC were sufficiently matched, and therefore no significant differences between the groups. Secondly, the low-cost and noninvasive EEG recording used to detect the high time resolution network dynamics of the brain is reliable; hence the results are applicable in clinical practice.


Our study reveals that the microstate topographies were aberrant in patients with MUD. In general, the present study is the first to explain changes in the microstates of drug dependence patients and its correlation with relevant clinical features. Therefore, we provide preliminary evidence for exploring whether EEG microstates could be significant indicators of cognitive disruption in substance-dependent patients.

Availability of data and materials

The datasets used and/ or analyzed during the current study are available from the corresponding author on reasonable request and with consent from the Institutional Review Board (IRB).



Chinese version of the Addiction Severity Index


Blood oxygen level-dependent


Continuous Paired Association Learning Task


The fifth edition of the Diagnostic and Statistical Manual of Mental Disorders






Global explained variance


Groton Maze Learning Task


Health controls


International Shopping List Task




Mean microstate duration


Magnetic resonance imaging


Methamphetamine use disorder


Positron emission tomography


Ratio of total time covered


Standard deviation


Social Emotional Cognition Task


Standardized low-resolution brain electromagnetic tomography


Two Back Task


  1. Yang X, Wang Y, Li Q, Zhong Y, Chen L, Du Y, He J, Liao L, Xiong K, Yi CX, et al. The Main molecular mechanisms underlying methamphetamine- induced neurotoxicity and implications for pharmacological treatment. Front Mol Neurosci. 2018;11:186.

    PubMed  PubMed Central  Google Scholar 

  2. Sabrini S, Wang GY, Lin JC, Ian JK, Curley LE. Methamphetamine use and cognitive function: a systematic review of neuroimaging research. Drug Alcohol Depend. 2019;194:75–87.

    CAS  PubMed  Google Scholar 

  3. Dean AC, Groman SM, Morales AM, London ED. An evaluation of the evidence that methamphetamine abuse causes cognitive decline in humans. Neuropsychopharmacol. 2013;38(2):259–74.

    CAS  Google Scholar 

  4. Bernheim A, See RE, Reichel CM. Chronic methamphetamine self-administration disrupts cortical control of cognition. Neurosci Biobehav Rev. 2016;69:36–48.

    CAS  PubMed  PubMed Central  Google Scholar 

  5. Ruan X, Zhong N, Yang Z, Fan X, Zhuang W, Du J, Jiang H, Zhao M. Gray matter volume showed dynamic alterations in methamphetamine users at 6 and 12months abstinence: a longitudinal voxel-based morphometry study. Prog Neuro-Psychopharmacol Biol Psychiatry. 2018;81:350–5.

    CAS  Google Scholar 

  6. Zhang Z, Li P, Fan L, Tao F, Li Y, Jiang B, Liu J, Kou Z. Change of brain structure imaging of long-term withdrawal of methamphetamine-dependent patients. Zhong nan da xue xue bao Yi xue ban. 2018;43(1):41–8.

    PubMed  Google Scholar 

  7. Salo R, Nordahl TE, Buonocore MH, Natsuaki Y, Waters C, Moore CD, Galloway GP, Leamon MH. Cognitive control and white matter callosal microstructure in methamphetamine-dependent subjects: a diffusion tensor imaging study. Biol Psychiatry. 2009;65(2):122–8.

    CAS  PubMed  Google Scholar 

  8. Foster BL, He BJ, Honey CJ, Jerbi K, Maier A, Saalmann YB. Spontaneous neural dynamics and multi-scale network organization. Front Syst Neurosci. 2016;10:7.

    PubMed  PubMed Central  Google Scholar 

  9. Cabral J, Kringelbach ML, Deco G. Exploring the network dynamics underlying brain activity during rest. Prog Neurobiol. 2014;114:102–31.

    PubMed  Google Scholar 

  10. Goldberg Hermo X, Lemos Giraldez S, Fananas Saura L. A systematic review of the complex organization of human cognitive domains and their heritability. Psicothema. 2014;26(1):1–9.

    PubMed  Google Scholar 

  11. Palva S, Palva JM. Discovering oscillatory interaction networks with M/EEG: challenges and breakthroughs. Trends Cogn Sci. 2012;16(4):219–30.

    PubMed  Google Scholar 

  12. Jie B, Liu M, Shen D. Integration of temporal and spatial properties of dynamic connectivity networks for automatic diagnosis of brain disease. Med Image Anal. 2018;47:81–94.

    PubMed  PubMed Central  Google Scholar 

  13. Zaytseva Y, Fajnerova I, Dvoracek B, Bourama E, Stamou I, Sulcova K, Motyl J, Horacek J, Rodriguez M, Spaniel F. Theoretical modeling of cognitive dysfunction in schizophrenia by means of errors and corresponding brain networks. Front Psychol. 2018;9:1027.

    PubMed  PubMed Central  Google Scholar 

  14. Yuan J, Li X, Zhang J, Luo L, Dong Q, Lv J, Zhao Y, Jiang X, Zhang S, Zhang W, et al. Spatio-temporal modeling of connectome-scale brain network interactions via time-evolving graphs. NeuroImage. 2018;180(Pt B):350–69.

    PubMed  Google Scholar 

  15. Proebstl L, Kamp F, Koller G, Soyka M. Cognitive deficits in methamphetamine users: how strong is the evidence? Pharmacopsychiatry. 2018;51(6):243–50.

    CAS  PubMed  Google Scholar 

  16. Hart CL, Marvin CB, Silver R, Smith EE. Is cognitive functioning impaired in methamphetamine users? A critical review. Neuropsychopharmacol. 2012;37(3):586–608.

    CAS  Google Scholar 

  17. Michel CM, Koenig T. EEG microstates as a tool for studying the temporal dynamics of whole-brain neuronal networks: A review. NeuroImage. 2018;180(Pt B):577–93.

    PubMed  Google Scholar 

  18. Lehmann D, Faber PL, Galderisi S, Herrmann WM, Kinoshita T, Koukkou M, Mucci A, Pascual-Marqui RD, Saito N, Wackermann J, et al. EEG microstate duration and syntax in acute, medication-naive, first-episode schizophrenia: a multi-center study. Psychiatry Res. 2005;138(2):141–56.

    PubMed  Google Scholar 

  19. Koenig T, Prichep L, Lehmann D, Sosa PV, Braeker E, Kleinlogel H, Isenhart R, John ER. Millisecond by millisecond, year by year: normative EEG microstates and developmental stages. NeuroImage. 2002;16(1):41–8.

    PubMed  Google Scholar 

  20. Brandeis D, Lehmann D, Michel CM, Mingrone W. Mapping event-related brain potential microstates to sentence endings. Brain Topogr. 1995;8(2):145–59.

    CAS  PubMed  Google Scholar 

  21. Khanna A, Pascual-Leone A, Michel CM, Farzan F. Microstates in resting-state EEG: current status and future directions. Neurosci Biobehav Rev. 2015;49:105–13.

    PubMed  Google Scholar 

  22. Britz J, Van De Ville D, Michel CM. BOLD correlates of EEG topography reveal rapid resting-state network dynamics. NeuroImage. 2010;52(4):1162–70.

    PubMed  Google Scholar 

  23. Lehmann D, Strik WK, Henggeler B, Koenig T, Koukkou M. Brain electric microstates and momentary conscious mind states as building blocks of spontaneous thinking: I. visual imagery and abstract thoughts. Int J Psychophysiol. 1998;29(1):1–11.

    CAS  PubMed  Google Scholar 

  24. Britz J, Landis T, Michel CM. Right parietal brain activity precedes perceptual alternation of bistable stimuli. Cerebral cortex (New York, NY : 1991). 2009;19(1):55–65.

    Google Scholar 

  25. Mohr C, Michel CM, Lantz G, Ortigue S, Viaud-Delmon I, Landis T. Brain state-dependent functional hemispheric specialization in men but not in women. Cerebral cortex (New York, NY : 1991). 2005;15(9):1451–8.

    Google Scholar 

  26. Muller TJ, Koenig T, Wackermann J, Kalus P, Fallgatter A, Strik W, Lehmann D. Subsecond changes of global brain state in illusory multistable motion perception. J Neural Transmission (Vienna, Austria : 1996). 2005;112(4):565–76.

    Google Scholar 

  27. Pascual-Marqui RD, Lehmann D, Faber P, Milz P, Kochi K, Yoshimura M, Nishida K, Isotani T, Kinoshita T. The resting microstate networks (RMN): cortical distributions, dynamics, and frequency specific information flow. ArXiv. 2014;1411:1949.

    Google Scholar 

  28. Mevel K, Grassiot B, Chetelat G, Defer G, Desgranges B. Eustache F: [the default mode network: cognitive role and pathological disturbances]. Rev Neurol. 2010;166(11):859–72.

    CAS  PubMed  Google Scholar 

  29. Kalechstein AD, De la Garza R 2nd, Newton TF, Green MF, Cook IA, Leuchter AF. Quantitative EEG abnormalities are associated with memory impairment in recently abstinent methamphetamine-dependent individuals. J Neuropsychiatry Clin Neurosci. 2009;21(3):254–8.

    PubMed  PubMed Central  Google Scholar 

  30. Yun K, Park HK, Kwon DH, Kim YT, Cho SN, Cho HJ, Peterson BS, Jeong J. Decreased cortical complexity in methamphetamine abusers. Psychiatry Res. 2012;201(3):226–32.

    CAS  PubMed  Google Scholar 

  31. Newton TF, Kalechstein AD, Hardy DJ, Cook IA, Nestor L, Ling W, Leuchter AF. Association between quantitative EEG and neurocognition in methamphetamine-dependent volunteers. Clin Neurophysiol. 2004;115(1):194–8.

    CAS  PubMed  Google Scholar 

  32. Newton TF, Cook IA, Kalechstein AD, Duran S, Monroy F, Ling W, Leuchter AF. Quantitative EEG abnormalities in recently abstinent methamphetamine dependent individuals. Clin Neurophysiol. 2003;114(3):410–5.

    PubMed  Google Scholar 

  33. Musaeus CS, Nielsen MS, Hogh P. Microstates as disease and progression markers in patients with mild cognitive impairment. Front Neurosci. 2019;13:563.

    PubMed  PubMed Central  Google Scholar 

  34. Muthukrishnan SP, Ahuja N, Mehta N, Sharma R. Functional brain microstate predicts the outcome in a visuospatial working memory task. Behav Brain Res. 2016;314:134–42.

    PubMed  Google Scholar 

  35. Schumacher J, Peraza LR, Firbank M, Thomas AJ, Kaiser M, Gallagher P, O'Brien JT, Blamire AM, Taylor JP. Dysfunctional brain dynamics and their origin in Lewy body dementia. Brain. 2019;142(6):1767–82.

    PubMed  PubMed Central  Google Scholar 

  36. Barry RJ, De Blasio FM. EEG differences between eyes-closed and eyes-open resting remain in healthy ageing. Biol Psychol. 2017;129:293–304.

    PubMed  Google Scholar 

  37. Tomescu MI, Rihs TA, Rochas V, Hardmeier M, Britz J, Allali G, Fuhr P, Eliez S, Michel CM. From swing to cane: sex differences of EEG resting-state temporal patterns during maturation and aging. Dev Cogn Neurosci. 2018;31:58–66.

    CAS  PubMed  PubMed Central  Google Scholar 

  38. Zhao M, Li X, Hao W, Wang Z, Zhang M, Xu D. A preliminary study of the reliability and validity of the addiction severity index. Chin Med Res J. 2004;4:679–80.

    Google Scholar 

  39. Zhong N, Jiang H, Du J, Zhao Y, Sun H, Xu D, Li C, Zhuang W, Li X, Hashimoto K, et al. The cognitive impairments and psychological wellbeing of methamphetamine dependent patients compared with health controls. Prog Neuro-Psychopharmacol Biol Psychiatry. 2016;69:31–7.

    CAS  Google Scholar 

  40. Haifeng J, Wenxu Z, Hong C, Chuanwei L, Jiang D, Haiming S, Zhikang C, Din X, Jijun W, Min Z. P300 event-related potential in abstinent methamphetamine-dependent patients. Physiol Behav. 2015;149:142–8.

    PubMed  Google Scholar 

  41. Arnaud D, Scott M. EEGLAB: an open source toolbox for analysis of single-trial EEG dynamics including independent component analysis. J Neurosci Methods. 2004;134(1):9–21.

    Google Scholar 

  42. Strik WK, Lehmann D. Data-determined window size and space-oriented segmentation of spontaneous EEG map series. Electroencephalogr Clin Neurophysiol. 1993;87(4):169–74.

    CAS  PubMed  Google Scholar 

  43. Lehmann D, Ozaki H, Pal I. EEG alpha map series: brain micro-states by space-oriented adaptive segmentation. Electroencephalogr Clin Neurophysiol. 1987;67(3):271–88.

    CAS  PubMed  Google Scholar 

  44. Wackermann J, Lehmann D, Michel CM, Strik WK. Adaptive segmentation of spontaneous EEG map series into spatially defined microstates. Int J Psychophysiol. 1993;14(3):269–83.

    CAS  PubMed  Google Scholar 

  45. Pascual-Marqui RD, Michel CM, Lehmann D. Segmentation of brain electrical activity into microstates: model estimation and validation. IEEE Trans Biomed Eng. 1995;42(7):658–65.

    CAS  PubMed  Google Scholar 

  46. Drissi NM, Szakacs A, Witt ST, Wretman A, Ulander M, Stahlbrandt H, Darin N, Hallbook T, Landtblom AM, Engstrom M. Altered brain microstate dynamics in adolescents with narcolepsy. Front Hum Neurosci. 2016;10:369.

    PubMed  PubMed Central  Google Scholar 

  47. Kikuchi M, Koenig T, Munesue T, Hanaoka A, Strik W, Dierks T, Koshino Y, Minabe Y. EEG microstate analysis in drug-naive patients with panic disorder. PLoS One. 2011;6(7):e22912.

    CAS  PubMed  PubMed Central  Google Scholar 

  48. Corradini PL, Persinger MA. Spectral power, source localization and microstates to quantify chronic deficits from 'mild' closed head injury: correlation with classic neuropsychological tests. Brain Inj. 2014;28(10):1317–27.

    PubMed  Google Scholar 

  49. Brunet D, Murray MM, Michel CM. Spatiotemporal analysis of multichannel EEG: CARTOOL. Computational Intelligence Neuroscience. 2011;2011:813870.

    PubMed  Google Scholar 

  50. Pascual-Marqui RD. Standardized low-resolution brain electromagnetic tomography (sLORETA): technical details. Methods findings Exp Clin Pharmacol. 2002;24(Suppl D):5–12.

    Google Scholar 

  51. Andreou C, Faber PL, Leicht G, Schoettle D, Polomac N, Hanganu-Opatz IL, Lehmann D, Mulert C. Resting-state connectivity in the prodromal phase of schizophrenia: insights from EEG microstates. Schizophr Res. 2014;152(2–3):513–20.

    PubMed  Google Scholar 

  52. Ahmadlou M, Ahmadi K, Rezazade M, Azad-Marzabadi E. Global organization of functional brain connectivity in methamphetamine abusers. Clin Neurophysiol. 2013;124(6):1122–31.

    PubMed  Google Scholar 

  53. Ipser JC, Uhlmann A, Taylor P, Harvey BH, Wilson D, Stein DJ. Distinct intrinsic functional brain network abnormalities in methamphetamine-dependent patients with and without a history of psychosis. Addict Biol. 2018;23(1):347–58.

    CAS  PubMed  Google Scholar 

  54. Koenig T, Studer D, Hubl D, Melie L, Strik WK. Brain connectivity at different time-scales measured with EEG. Philos Trans R Soc Lond Ser B Biol Sci. 2005;360(1457):1015–23.

    CAS  Google Scholar 

  55. Stevens A, Lutzenberger W, Bartels DM, Strik W, Lindner K. Increased duration and altered topography of EEG microstates during cognitive tasks in chronic schizophrenia. Psychiatry Res. 1997;66(1):45–57.

    CAS  PubMed  Google Scholar 

  56. Irisawa S, Isotani T, Yagyu T, Morita S, Nishida K, Yamada K, Yoshimura M, Okugawa G, Nobuhara K, Kinoshita T. Increased omega complexity and decreased microstate duration in nonmedicated schizophrenic patients. Neuropsychobiol. 2006;54(2):134–9.

    Google Scholar 

  57. Nishida K, Morishima Y, Yoshimura M, Isotani T, Irisawa S, Jann K, Dierks T, Strik W, Kinoshita T, Koenig T. EEG microstates associated with salience and frontoparietal networks in frontotemporal dementia, schizophrenia and Alzheimer's disease. Clin Neurophysiol. 2013;124(6):1106–14.

    PubMed  Google Scholar 

  58. Li J, Zhu L, Guan F, Yan Z, Liu D, Han W, Chen T. Relationship between schizophrenia and changes in the expression of the long non-coding RNAs Meg3, Miat, Neat1 and Neat2. J Psychiatr Res. 2018;106:22–30.

    PubMed  Google Scholar 

  59. Dore G, Sweeting M. Drug-induced psychosis associated with crystalline methamphetamine. Australasian Psychiatry. 2006;14(1):86–9.

    PubMed  Google Scholar 

  60. Srisurapanont M, Arunpongpaisal S, Wada K, Marsden J, Ali R, Kongsakon R. Comparisons of methamphetamine psychotic and schizophrenic symptoms: a differential item functioning analysis. Prog Neuro-Psychopharmacol Biol Psychiatry. 2011;35(4):959–64.

    CAS  Google Scholar 

  61. Yoshimura M, Pascual-Marqui RD, Nishida K, Kitaura Y, Mii H, Saito Y, Ikeda S, Katsura K, Ueda S, Minami S, et al. Hyperactivation of the frontal control network revealed by symptom provocation in obsessive-compulsive disorder using EEG microstate and sLORETA analyses. Neuropsychobiology. 2019;77(4):176-85.

  62. Salo R, Buonocore MH, Leamon M, Natsuaki Y, Waters C, Moore CD, Galloway GP, Nordahl TE. Extended findings of brain metabolite normalization in MA-dependent subjects across sustained abstinence: a proton MRS study. Drug Alcohol Depend. 2011;113(2–3):133–8.

    PubMed  Google Scholar 

  63. Stern ER, Muratore AF, Taylor SF, Abelson JL, Hof PR, Goodman WK. Switching between internally and externally focused attention in obsessive-compulsive disorder: abnormal visual cortex activation and connectivity. Psychiatry research Neuroimaging. 2017;265:87–97.

    PubMed  Google Scholar 

  64. Zilverstand A, Huang AS, Alia-Klein N, Goldstein RZ. Neuroimaging impaired response inhibition and salience attribution in human drug addiction: a systematic review. Neuron. 2018;98(5):886–903.

    CAS  PubMed  PubMed Central  Google Scholar 

  65. Hazari N, Narayanaswamy JC, Venkatasubramanian G. Neuroimaging findings in obsessive-compulsive disorder: A narrative review to elucidate neurobiological underpinnings. Indian J Psychiatry. 2019;61(Suppl 1):S9–s29.

    PubMed  PubMed Central  Google Scholar 

  66. Liao J, Li T, Dong W, Wang J, Tian J, Liu J, Quan W, Yan J. Reduced prefrontal-temporal cortical activation during verbal fluency task in obsessive-compulsive disorder: a multi-channel near-infrared spectroscopy study. J Psychiatr Res. 2019;109:33–40.

    PubMed  Google Scholar 

  67. Nagase S. Effect of methamphetamine on the auditory system. Nihon Jibiinkoka Gakkai kaiho. 1993;96(2):203–11.

    CAS  PubMed  Google Scholar 

  68. Van de Ville D, Britz J, Michel CM. EEG microstate sequences in healthy humans at rest reveal scale-free dynamics. Proc Natl Acad Sci U S A. 2010;107(42):18179–84.

    PubMed  Google Scholar 

  69. Yuan H, Zotev V, Phillips R, Drevets WC, Bodurka J. Spatiotemporal dynamics of the brain at rest--exploring EEG microstates as electrophysiological signatures of BOLD resting state networks. NeuroImage. 2012;60(4):2062–72.

    PubMed  Google Scholar 

  70. Huang X, Chen YY, Shen Y, Cao X, Li A, Liu Q, Li Z, Zhang LB, Dai W, Tan T, et al. Methamphetamine abuse impairs motor cortical plasticity and function. Mol Psychiatry. 2017;22(9):1274–81.

    CAS  PubMed  PubMed Central  Google Scholar 

  71. Hadamitzky M, Markou A, Kuczenski R. Extended access to methamphetamine self-administration affects sensorimotor gating in rats. Behav Brain Res. 2011;217(2):386–90.

    CAS  PubMed  Google Scholar 

  72. Nakama H, Chang L, Fein G, Shimotsu R, Jiang CS, Ernst T. Methamphetamine users show greater than normal age-related cortical gray matter loss. Addiction (Abingdon, England). 2011;106(8):1474–83.

    Google Scholar 

  73. Daumann J, Koester P, Becker B, Wagner D, Imperati D, Gouzoulis-Mayfrank E, Tittgemeyer M. Medial prefrontal gray matter volume reductions in users of amphetamine-type stimulants revealed by combined tract-based spatial statistics and voxel-based morphometry. NeuroImage. 2011;54(2):794–801.

    CAS  PubMed  Google Scholar 

  74. Gao F, Jia H, Wu X, Yu D, Feng Y. Altered resting-state EEG microstate parameters and enhanced spatial complexity in male adolescent patients with mild spastic Diplegia. Brain Topogr. 2017;30(2):233–44.

    PubMed  Google Scholar 

  75. Wingert JR, Sinclair RJ, Dixit S, Damiano DL, Burton H. Somatosensory-evoked cortical activity in spastic diplegic cerebral palsy. Hum Brain Mapp. 2010;31(11):1772–85.

    PubMed  PubMed Central  Google Scholar 

  76. Himmelmann K, Uvebrant P. Function and neuroimaging in cerebral palsy: a population-based study. Dev Med Child Neurol. 2011;53(6):516–21.

    PubMed  Google Scholar 

  77. Han X, Berg AC, Oh H, Samaras D, Leung HC. Multi-voxel pattern analysis of selective representation of visual working memory in ventral temporal and occipital regions. NeuroImage. 2013;73:8–15.

    PubMed  PubMed Central  Google Scholar 

  78. Kikuchi M, Koenig T, Wada Y, Higashima M, Koshino Y, Strik W, Dierks T. Native EEG and treatment effects in neuroleptic-naive schizophrenic patients: time and frequency domain approaches. Schizophr Res. 2007;97(1–3):163–72.

    PubMed  Google Scholar 

Download references


The authors would like to thank all the patients and the healthy volunteers who took part in the sample collection and test.


The authors are supported by the National Key R&D Program of China (2017YFC1310400), National Nature Science Foundation (81771436, 81501148, U1502228, 81801319, 81601164), Shanghai Municipal Health and Family Planning Commission (2017ZZ02021), Municipal Human Resources Development Program for Outstanding Young Talents in Medical and Health Sciences in Shanghai (2017YQ013), Shanghai Key Laboratory of Psychotic Disorders (13DZ2260500), Program of Shanghai Academic Research Leader (17XD1403300), Clinical Research Center of Shanghai Jiao Tong University School of Medicine (DLY201818), Shanghai Science and Technology Committee (18411961200), Shanghai Municipal Science and Technology Major Project (2018SHZDZX05), and Shanghai Clinical Research Center for Mental Health (19MC1911100). The funders had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations



TC, HS, MZ, and HJ conceived and designed the research. TC, HS, HT, XL and YM achieved the data acquisition. TC and NZ performed the data analyses. CD, CZ, WS and JZ assisted with experiments. JB, DX, TL and QZ provided technical support. TC and HS drafted the manuscript. MZ and HJ revised the manuscript. All authors provided critical revision of the manuscript and approved the final version.

Corresponding author

Correspondence to Haifeng Jiang.

Ethics declarations

Ethics approval and consent to participate

The study protocol was implemented in conformity to the principles of the Declaration of Helsinki and approved by the Ethic Committee of the Shanghai Mental Health Center. Informed consent was signed by all subjects prior to participation in this research.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

Brain regions that showed increased activation of microstate A at rest in patients with methamphetamine use disorder. Table S2. Brain regions that showed increased activation of microstate B at rest in patients with methamphetamine use disorder. Table S3. Brain regions that showed increased activation of microstate C at rest in patients with methamphetamine use disorder. Fig. S1. sLORETA differences in the cortical distribution of electrical activity sources of microstate class A between patients with methamphetamine use disorder and health controls. Significant differences are displayed (yellow = patients>controls). Fig. S2. sLORETA differences in the cortical distribution of electrical activity sources of microstate class B between patients with methamphetamine use disorder and health controls. Significant differences are displayed (yellow = patients>controls). Fig. S3. sLORETA differences in the cortical distribution of electrical activity sources of microstate class C between patients with methamphetamine use disorder and health controls. Significant differences are displayed (yellow = patients>controls).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, T., Su, H., Zhong, N. et al. Disrupted brain network dynamics and cognitive functions in methamphetamine use disorder: insights from EEG microstates. BMC Psychiatry 20, 334 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: