Abstract

Fractal patterns have been shown to change in resting- and task-state blood oxygen level-dependent signals in bipolar disorder patients. However, fractal characteristics of brain blood oxygen level-dependent signals when responding to external emotional stimuli in pediatric bipolar disorder remain unclear. Blood oxygen level-dependent signals of 20 PBD-I patients and 17 age- and sex-matched healthy controls were extracted while performing an emotional Go-Nogo task. Neural responses relevant to the task and Hurst exponent of the blood oxygen level-dependent signals were assessed. Correlations between clinical indices and Hurst exponent were estimated. Significantly increased activations were found in regions covering the frontal lobe, parietal lobe, temporal lobe, insula, and subcortical nuclei in PBD-I patients compared to healthy controls in contrast of emotional versus neutral distractors. PBD-I patients exhibited higher Hurst exponent in regions that involved in action control, such as superior frontal gyrus, inferior frontal gyrus, inferior temporal gyrus, and insula, with Hurst exponent of frontal orbital gyrus correlated with onset age. The present study exhibited overactivation, increased self-similarity and decreased complexity in cortical regions during emotional Go-Nogo task in patients relative to healthy controls, which provides evidence of an altered emotional modulation of cognitive control in pediatric bipolar disorder patients. Hurst exponent may be a fractal biomarker of neural activity in pediatric bipolar disorder.

Introduction

Bipolar disorder (BD) is a generalized neuropsychiatric disorder in which individuals exhibit inattention and decreased motivational behavior or goal-orientation in multiple functional areas (Yurgelun-Todd and Ross 2006; Zhang et al. 2011). BD can be divided into mania (BD-I) and hypomania (BD-II) based on whether there is a full-blown manic episode. BD-I is characterized by periodic episodes of mania and depression, usually accompanied by periods of relative relief of symptoms or euphoria (Robinson et al. 2008). It is believed that people with BD exhibit symptoms of emotional dysregulation, further leading to deficits in emotional recognition and assessment (Poletti et al. 2018).

Understanding how emotion is regulated is important in BD. It is believed that emotional dysregulation causes cognitive impairment through interaction between cognitive and emotional brain networks (Pavuluri et al. 2007; Wessa et al. 2007). Compared with controls, euthymic pediatric bipolar disorder (PBD) exhibited increased activation in higher-order visual perception areas and decreased activation in visual areas of the occipital cortex in a face processing task. This suggests that the patient has an inability to regulate emotional responses to the face and struggles to control sensory input (Pavuluri et al. 2007). In addition, in terms of emotion recognition and assessment, facial emotion is often used to study emotional processing. Assessing facial emotion recognition ability is helpful to understand the emotional regulation of BD patients (Pavuluri et al. 2007; Delvecchio et al. 2013). It was demonstrated that BD-I patients have difficulty in recognizing fearful and angry faces, and BD patients tend to judge extreme emotions as moderate or mild (Rocca et al. 2009). PBD patients showed greater amygdala activity and decreased ventrolateral prefrontal cortex (VLPFC) activity as participants passively viewed happy and angry faces (Pavuluri et al. 2007). In the recognition of negative emotions, the right inferior frontal gyrus (IFG) activation was reduced in BD patients (Hafeman et al. 2014). Researches using emotional Go-Nogo tasks have demonstrated that in euthymic BD, Nogo minus Go contrast revealed greater activation in IFG, insula, and orbitofrontal cortex (Hummer et al. 2013), comparison of emotional versus non-emotional distractors revealed greater activity in prefrontal cortex, cingulate cortex, parietal cortex and insula, and comparison between emotional and neutral distractors revealed an increased neural response in a fronto-striatal circuit (Wessa et al. 2007). In short, emotional regulation is a complex process, which is a “buffering” system that can respond flexibly to constantly changing conditions (Ortiz and Alda 2018). The system of BD become more organized over a longer time scale and less flexible, which different from that of healthy subjects (Bayani et al. 2017). The process of emotional variability might be better described using fractal dynamics.

Many physiological signals demonstrate fractal fluctuations on a wide time scale from seconds to hours, which is typical of self-similarity time patterns or statistical characteristics. Evidence suggested that the brain and its functional activities express fractal behavior on a limited scale. Fractal fluctuations are very complex, neither too random nor too regular. It is believed that fractal regulation may reflect integrity and adaptability of a physiological system in terms of maintaining internal stability meanwhile responding to external alterations, while the degradation of fractal regulation signifies the loss of functional flexibility and increased vulnerability to destructive events (Goldberger et al. 2002; Knapen et al. 2021). It has been shown that blood oxygen level-dependent (BOLD) signal of the brain exhibits fractal-like dynamics. That means the self-similarity in the time course of the BOLD signal (Gentili et al. 2015; Long et al. 2022) and it may be illustrated by a widely used indicator, Hurst exponent (HE), an index in the range of 0 to 1. According to the HE values, time series can be divided into three categories. (i) HE equals or close to 0.5 suggests a white noise series, indicating a process of with less memory and low predictability. (ii) HE between 0 and 0.5 indicates an anti-correlation time series with a dynamic of reversal in time. (iii) HE between 0.5 and 1 indicates a positively correlated series with long memory and it is more predictable (Gentili et al. 2015). By present changes in the degree of correlation in the signal, the HE value of the BOLD signal specifies if and how internal and external stressors affect brain function (Campbell and Weber 2022). HE has been used to investigate randomness and self-similarity of BOLD time series in healthy and diseased subjects and some significant results were reported. It was found that HE of resting-state fMRI signal was greater in gray matter than in white matter regions (Wink et al. 2008). The average HE value of the gray matter of whole brain became larger along with the increase of age. The most significant increase located in the parietal lobe and frontal lobe and some brain regions showed a downward trend of HE, mainly in emotion regulated brain regions including the temporal lobe, limbic lobe, and insula (Dong et al. 2018). Compared to healthy controls, higher HE of resting-state fMRI signals in medial and lateral temporal lobes, insula, dorsal cingulate/medial premotor cortex was found in Alzheimer’s disease subjects (Maxim et al. 2005). And higher HE at mean whole brain and regional levels was demonstrated in schizophrenia people while performing the Cyberball social exclusion task (Sokunbi et al. 2014). Furthermore, structural magnetic resonance imaging (MRI) research indicated increases in fractal dimensions (in patients relative to controls) in left lateral orbitofrontal cortex and right precuneus, and decreases in right caudal middle frontal, entorhinal cortex, and right pars orbitalis, and left fusiform and posterior cingulate cortices in BD (Nenadic et al. 2017). Considering the inherent relationship between brain structure and brain function, we speculate that the brain activity of BD, such as the neural response to external stimuli, will exhibit abnormal fractal characteristics in brain regions implicated in emotion regulation and cognitive control. There are few reports on this issue and further research is needed.

To test the hypotheses, a functional MRI (fMRI) Go-Nogo paradigm comprising both emotional and non-emotional task conditions was conducted in the present study. Fractal index of the BOLD signals of brain regions showing regional activations was analyzed, and relation between the fractal index and clinical indices was evaluated.

Materials and methods

Subjects

The data of the present study included 20 PBD-I patients from the Child and Adolescent Psychiatry Department of the Second Xiangya Hospital of Central South University in Changsha, Hunan Province. Seventeen age- and sex-matched healthy controls were randomly recruited from local secondary schools. The data for this study were collected from January 2012 to July 2014.

Inclusion criteria for PBD patients included the following: (i) meeting the criteria for BD of the Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition (DSM-IV); (ii) right-handed; (iii) aged 10–18 years; (iv) keeping body still during MRI scan; (v) Han nationality. Exclusion criteria for all subjects included the following: (i) IQ < 80; (ii) MRI contraindications such as claustrophobia, stents, and metal implants; (iii) having a neurological disease, such as injury to the head; (iv) other mental disorders, such as schizophrenia, personality disorders, or autism; (v) alcohol or drug addiction; and (vi) having received electroconvulsive therapy.

The present study was approved by the ethic committee of the Second Xiangya Hospital of Central South University. Written informed consent from all subjects and their guardians were obtained.

Clinical evaluations

Two experienced psychiatrists diagnosed all PBD patients by implementing the Emotional Disorders and Schizophrenia Scale for School-Age Children-Current and Lifetime Versions (Kaufman et al. 1997). The subjects and at least one guardian underwent a diagnostic semi-structured interview. Demographic information and clinical data were collected by questionnaires. The intelligence of all participants was assessed using the Wechsler Reduced Scale of Intelligence (WASI). Before the MRI scan, the Mood and Senses Questionnaire (MFQ; Chang et al. 2005) and Young Mania Rating Scale (YMRS; Young et al. 1978) were used to assess the severity of depression and mania symptoms respectively. YMRS is an 11-item clinician rating scale, and the score correlates with symptoms of mania. The higher the score, the more severe the degree of mania. MFQ is a 33-item measure to assess emotional symptoms in subjects. The more severe the degree of depression, the higher the score of the scale. Demographic data such as age of onset, duration of illness, number of episodes and psychotic symptoms were recorded. One patient fulfilled the criteria for a comorbid attention deficit hyperactivity disorder (ADHD) and generalized anxiety disorder (GAD). Two patients have comorbidity ADHD and GAD respectively. Three patients were medication-free, and 17 patients were taking psychotropic medication. Among those who were on regular medication, 1 received antidepressant, 1 received atypical neuroleptics, and 15 received a combined pharmacological treatment (e.g. lithium, valproate, and atypical antipsychotics).

Go-Nogo task paradigm

The task flow of fMRI in emotional Go-Nogo paradigm is shown in Fig. 1. Participants were asked to identify happy, sad and neutral facial expressions which from pictures of facial affect (Ekman and Friesen 1976), and perform keystroke tasks, including pressing a button when the target face was present and not responding when the distractor was present. The task consisted of eight blocks, six of which were emotional blocks (Go) with materials contained happy, sad and neutral faces, and the other two were control blocks (Nogo) with materials contained male, and female faces with neutral emotions. Each block lasts for 36 s and contains 16 targets and 8 distractors, which constantly appear through a pseudo-random sequence. Each stimulation lasted for 500 ms, and the stimulation interval was 1,000 ms. The subjects were asked to read the instructions (4,000 ms) before the start of each block and subsequently either respond or not respond to the stimulus. E-prime software was used to record the response time and accuracy of each subject.

Task flow of fMRI in emotional Go/Nogo paradigm. Each stimulation lasted for 500 ms, and the stimulation interval was 1,000 ms. instructions lasts for 4,000 ms before the start of each block. Each block lasts for 36 s and contains 16 targets and 8 distractors, which constantly appear through a pseudo-random sequence.
Fig 1

Task flow of fMRI in emotional Go/Nogo paradigm. Each stimulation lasted for 500 ms, and the stimulation interval was 1,000 ms. instructions lasts for 4,000 ms before the start of each block. Each block lasts for 36 s and contains 16 targets and 8 distractors, which constantly appear through a pseudo-random sequence.

MRI data acquisition and preprocessing

MRI images were collected on a 3.0 T Siemens Trio scanner (Siemens, Munich, Germany). During the MRI scanning, foam pads were placed on either side of the subject’s head to restrict head movement. Cotton plugs were used to reduce noise and protect the subject’s hearing. Structural images were obtained by using 3D magnetization and rapid acquisition gradient echo (3D MPRAGE) protocol. Detailed parameters covering the whole brain include repetition time (TR) =2,300 ms, echo time (TE) =2.03 ms, field of view (FOV) =256 × 256 mm2, flip angle =9o, matrix size =256 × 256, slice thickness = 1 mm, and gap =0 mm. Functional images were performed using T2-weighted single echo gradient imaging sequences with scan parameters as: TR = 2,000 ms, TE = 30 ms, FOV = 240 mm × 240 mm, flip angle = 90o, matrix size =64 × 64, slice thickness = 4 mm, and gap = 0.4 mm. The functional scanning lasted for 620 s, yielding a total of 310 volumes.

The fMRI signals were preprocessed in the following steps by SPM12. The first three volumes of each subject’s fMRI image were discarded in consideration of magnetization equilibration effects. The images were then sectioned for slice timing, realignment, covariate regression (Friston 24 parameter model, white matter signal, cerebrospinal fluid signal, and global signal) and spatial normalization. Finally, fMRI images were spatially smoothed using a 6 mm Gauss kernel FWHM.

BOLD changes associated with Go-Nogo task

SPM12 was used to analyze fMRI data by two level general linear model. At the subject level, the images of each subject were contrasted in three different ways. First, contrasts were made to observe task states versus resting states, which enables them to evaluate the influence of emotional valence in the task. Next, emotional versus non-emotional states contrasts were created to investigate the influence of emotional factors on the task performance. At last, emotional distractor versus neutral distractor were compared in emotional tasks to analyze the effect of emotional titer on response suppression. These outputs were input into second-level analyses to study neural correlates of tasks in each group (one-sample t-test) and determine regions that were significantly different between the PBD-I patients and the HCs (two-sample t-tests) (α < 0.05 with false discovery rate (FDR) corrected).

The brain atlas (Fan et al. 2016) that partitions the cerebral cortex into 246 ROIs, including 210 cortical sub-regions and 36 subcortical sub-regions, was used to identify brain regions where BOLD signals were significantly associated with tasks in the PBD-I and HCs groups. After overlapping the group activation map and the atlas, the fMRI signal of a region of interest (ROI) was extracted by averaging the intensity of all voxel signals in the ROI at each time point and used for further analyses.

HE calculation

Dispersion analysis, an effective method to detect the temporal complexity of a time series with requiring smaller datasets and exhibiting less bias (Bassingthwaighte and Raymond 1995), was utilized to derive the HE of the BOLD data series in the present study. It calculates descriptive statistical measures as a function of scale with local measure is the mean of time series of a window and the partition-based measure is the standard deviation of local means (Eke et al. 2002).The time series is divided into nonoverlapping windows of size m, and the standard deviation within the window SD(m) depends on the window size, m, in a power–law relationship:

The following practical formula is obtained by taking the logarithm

where m is the element size used to calculate SD, and m0 is the arbitrarily chosen reference size that is the number of points aggregated together in a time series (Bassingthwaighte and Raymond 1995). The HE may be read off as the slope of the line on the log–log plot. In the present study, m = 2 and m0 = 1.

Statistical analysis

Demographic data and task performance scores were analyzed using Chi-square test and independent sample t-test. A second-level random effects analysis was conducted to create task related T-statistics maps of the fMRI data, followed by a one-sample T-test and a two-sample T-test in SPM12. Brain regions where the BOLD signal showed significant correlation with the task and differences between PBD and HC groups were set ROIs. Averaged BOLD signals of the ROIs were extracted and the HE of the time series of the ROIs were calculated. Permutation testing was performed to test the null hypothesis that the HE showed no significant difference between the two groups other than by chance. The values were randomly relabeled 5,000 times and a distribution of the test statistic was obtained. The permutated P-values were obtained by the occupied null models with less than 5% or greater than 95% and corrected by FDR with P < 0.05. Correlations between clinical indices and values of HE were estimated with age and gender regressed. Due to the small sample of each group in the present study, effect size (Cohen’s d) was calculated to discover between-group difference of HE. A P-value < 0.05 was considered significant for all statistical analyses.

Results

Demographics and clinical data

As shown in Table 1, no significant difference was found in gender (P = 0.254), age (P = 0.064), educational years (P = 0.232), IQ (P = 0.703), and MFQ scores (P = 0.742) between the two groups (PBD-I, HC). There were significant differences between the two groups in scores of YMRS (P < 0.001). No significant difference was found in error of omission, error of response and reaction time between the two groups.

Table 1

Demographic characteristics and clinical variables.

CharacteristicsPBD-I (n = 20)HCs (n = 17)Statistic (T2)P-value
Gender (male/female)12/87/101.3030.254
Age (years)15.25 ± 1.6514.24 ± 1.561.9100.064
Education (years)8.15 ± 1.507.59 ± 1.281.2160.232
IQ105.60 ± 10.29106.71 ± 6.33−0.3850.703
YMRS scores25.43 ± 14.087.88 ± 3.305.010<0.001
MFQ scores6.60 ± 3.526.24 ± 3.110.3310.742
Onset age (years)13.90 ± 1.74
Illness duration (months)16.85 ± 13.41
Number of episodes2.60 ± 1.31
Psychotic symptoms (yes/no)6/14
Task performance
Emotional targets (Go)5.16 ± 1.734.26 ± 2.281.350.18
Neutral targets (Go)2.89 ± 1.873.27 ± 1.94−0.610.54
Control targets (Go)1.88 ± 1.051.64 ± 1.040.660.51
False response errors (%)
Emotional distractors (Nogo)3.18 ± 1.012.93 ± 1.510.590.55
Neutral distractors (Nogo)2.40 ± 1.322.14 ± 0.910.690.49
Control distractors (Nogo)2.26 ± 0.782.46 ± 1.33−0.560.57
Reaction time (ms)
Emotional targets (Go)590.09 ± 8.05604.63 ± 62.10−0.590.53
Neutral targets (Go)591.78 ± 78.88592.27 ± 69.02−0.020.98
Control targets (Go)565.57 ± 66.19555.72 ± 65.860.450.6
CharacteristicsPBD-I (n = 20)HCs (n = 17)Statistic (T2)P-value
Gender (male/female)12/87/101.3030.254
Age (years)15.25 ± 1.6514.24 ± 1.561.9100.064
Education (years)8.15 ± 1.507.59 ± 1.281.2160.232
IQ105.60 ± 10.29106.71 ± 6.33−0.3850.703
YMRS scores25.43 ± 14.087.88 ± 3.305.010<0.001
MFQ scores6.60 ± 3.526.24 ± 3.110.3310.742
Onset age (years)13.90 ± 1.74
Illness duration (months)16.85 ± 13.41
Number of episodes2.60 ± 1.31
Psychotic symptoms (yes/no)6/14
Task performance
Emotional targets (Go)5.16 ± 1.734.26 ± 2.281.350.18
Neutral targets (Go)2.89 ± 1.873.27 ± 1.94−0.610.54
Control targets (Go)1.88 ± 1.051.64 ± 1.040.660.51
False response errors (%)
Emotional distractors (Nogo)3.18 ± 1.012.93 ± 1.510.590.55
Neutral distractors (Nogo)2.40 ± 1.322.14 ± 0.910.690.49
Control distractors (Nogo)2.26 ± 0.782.46 ± 1.33−0.560.57
Reaction time (ms)
Emotional targets (Go)590.09 ± 8.05604.63 ± 62.10−0.590.53
Neutral targets (Go)591.78 ± 78.88592.27 ± 69.02−0.020.98
Control targets (Go)565.57 ± 66.19555.72 ± 65.860.450.6

Values are presented as mean ± standard deviation. Except for the Chi-squared test for gender, independent sample t-test was used for comparison of other variables. PBD-I, pediatric bipolar disorder type I; HCs: healthy controls; IQ: intelligence quotient; YMRS: Young Manic Rating Scale; MFQ: Mood and Feelings Questionnaire.

Table 1

Demographic characteristics and clinical variables.

CharacteristicsPBD-I (n = 20)HCs (n = 17)Statistic (T2)P-value
Gender (male/female)12/87/101.3030.254
Age (years)15.25 ± 1.6514.24 ± 1.561.9100.064
Education (years)8.15 ± 1.507.59 ± 1.281.2160.232
IQ105.60 ± 10.29106.71 ± 6.33−0.3850.703
YMRS scores25.43 ± 14.087.88 ± 3.305.010<0.001
MFQ scores6.60 ± 3.526.24 ± 3.110.3310.742
Onset age (years)13.90 ± 1.74
Illness duration (months)16.85 ± 13.41
Number of episodes2.60 ± 1.31
Psychotic symptoms (yes/no)6/14
Task performance
Emotional targets (Go)5.16 ± 1.734.26 ± 2.281.350.18
Neutral targets (Go)2.89 ± 1.873.27 ± 1.94−0.610.54
Control targets (Go)1.88 ± 1.051.64 ± 1.040.660.51
False response errors (%)
Emotional distractors (Nogo)3.18 ± 1.012.93 ± 1.510.590.55
Neutral distractors (Nogo)2.40 ± 1.322.14 ± 0.910.690.49
Control distractors (Nogo)2.26 ± 0.782.46 ± 1.33−0.560.57
Reaction time (ms)
Emotional targets (Go)590.09 ± 8.05604.63 ± 62.10−0.590.53
Neutral targets (Go)591.78 ± 78.88592.27 ± 69.02−0.020.98
Control targets (Go)565.57 ± 66.19555.72 ± 65.860.450.6
CharacteristicsPBD-I (n = 20)HCs (n = 17)Statistic (T2)P-value
Gender (male/female)12/87/101.3030.254
Age (years)15.25 ± 1.6514.24 ± 1.561.9100.064
Education (years)8.15 ± 1.507.59 ± 1.281.2160.232
IQ105.60 ± 10.29106.71 ± 6.33−0.3850.703
YMRS scores25.43 ± 14.087.88 ± 3.305.010<0.001
MFQ scores6.60 ± 3.526.24 ± 3.110.3310.742
Onset age (years)13.90 ± 1.74
Illness duration (months)16.85 ± 13.41
Number of episodes2.60 ± 1.31
Psychotic symptoms (yes/no)6/14
Task performance
Emotional targets (Go)5.16 ± 1.734.26 ± 2.281.350.18
Neutral targets (Go)2.89 ± 1.873.27 ± 1.94−0.610.54
Control targets (Go)1.88 ± 1.051.64 ± 1.040.660.51
False response errors (%)
Emotional distractors (Nogo)3.18 ± 1.012.93 ± 1.510.590.55
Neutral distractors (Nogo)2.40 ± 1.322.14 ± 0.910.690.49
Control distractors (Nogo)2.26 ± 0.782.46 ± 1.33−0.560.57
Reaction time (ms)
Emotional targets (Go)590.09 ± 8.05604.63 ± 62.10−0.590.53
Neutral targets (Go)591.78 ± 78.88592.27 ± 69.02−0.020.98
Control targets (Go)565.57 ± 66.19555.72 ± 65.860.450.6

Values are presented as mean ± standard deviation. Except for the Chi-squared test for gender, independent sample t-test was used for comparison of other variables. PBD-I, pediatric bipolar disorder type I; HCs: healthy controls; IQ: intelligence quotient; YMRS: Young Manic Rating Scale; MFQ: Mood and Feelings Questionnaire.

BOLD changes associated with Go-Nogo task

No significant regional brain activation difference was found in the contrast of all Go-Nogo conditions versus the rest and contrast emotional and non-emotional revealed between the PBD-I patients and HCs.

In the contrast of emotional distractor versus neutral distractor during Go-Nogo task, clusters that covering frontal lobe, parietal lobe, temporal lobe, limbic lobe, insula, and subcortical nucleus demonstrated increased brain regional activations in PBD1 patients than that of the HCs (FDR-corrected, Fig. 2). The brain atlas that partitions the cerebral cortex into 246 ROIs was used to identify the activated brain regions. After overlapping the group activation map and the atlas, 13 brain regions were survived including A10m (superior frontal gyrus, SFG), A45r (IFG), A12/47 l (frontal orbital gyrus, FOG), A4hf (precentral gyrus, PrG), A4ul (PrG), A1/2/3ll (paracentral lobule, PCL), A20iv (inferior temporal gyrus, ITG), TL (parahippocampal gyrus, PhG), A39c (inferior parietal lobule, IPL), A39rd (IPL), dIa (insula, INS), mAmyg (amygdala, Amyg), and cHipp (hippocampus, Hipp). Details of the 13 regions were described in Table 2. The fMRI signals the 13 brain regions were extracted for further HE calculation to evaluate the fractal dynamic characteristics of the signal.

Increased activations of PBD1 patients than HCs in contrast of emotional distractor versus neutral distractor during Go Nogo task (FDR-corrected). Active clusters were found covering frontal lobe, parietal lobe, temporal lobe, limbic lobe, insula, and subcortical nucleus. The brain atlas that partitions the cerebral cortex into 246 ROIs was used to identify the activated brain regions. After overlapping the group activation map and the atlas, 13 brain regions were survived including A10m (SFG), A45r (IFG), A12/47 l (FOG), A4hf (PrG), A4ul (PrG), A1/2/3ll (PCL), A20iv (ITG), TL (PhG), A39c (IPL),A39rd (IPL), dIa (insula), mAmyg (amygdala), and cHipp (hippocampus). The fMRI signals the brain regions were extracted for further HE calculation.
Fig 2

Increased activations of PBD1 patients than HCs in contrast of emotional distractor versus neutral distractor during Go Nogo task (FDR-corrected). Active clusters were found covering frontal lobe, parietal lobe, temporal lobe, limbic lobe, insula, and subcortical nucleus. The brain atlas that partitions the cerebral cortex into 246 ROIs was used to identify the activated brain regions. After overlapping the group activation map and the atlas, 13 brain regions were survived including A10m (SFG), A45r (IFG), A12/47 l (FOG), A4hf (PrG), A4ul (PrG), A1/2/3ll (PCL), A20iv (ITG), TL (PhG), A39c (IPL),A39rd (IPL), dIa (insula), mAmyg (amygdala), and cHipp (hippocampus). The fMRI signals the brain regions were extracted for further HE calculation.

Table 2

Thirteen brain regions showing increased activations in PBD-I patients after overlapping the group activation map and Human Brainnetome Atlas.

LobeGyrusCytoarchitectonic namePeak T valueCluster size
Frontal lobeSuperior frontal gyrusL.A10m2.66784
Inferior frontal gyrusL.A45r1.48243
Orbital gyrusR. A12/47 l2.27264
Precentral gyrusL.A4hf3.994121
L.A4ul3.303108
Paracentral lobuleL. A1/2/3ll2.58818
Temporal lobeInferior temporal gyrusL.A20iv2.97137
Parahippocampal gyrusL.TL2.39731
Parietal lobeInferior parietal lobuleL.A39c3.01126
R. A39rd2.674128
Insular lobeInsular gyrusR. dIa1.78960
Subcortical nucleiAmygdalaR. mAmyg0.5576
HippocampusL. cHipp3.79732
LobeGyrusCytoarchitectonic namePeak T valueCluster size
Frontal lobeSuperior frontal gyrusL.A10m2.66784
Inferior frontal gyrusL.A45r1.48243
Orbital gyrusR. A12/47 l2.27264
Precentral gyrusL.A4hf3.994121
L.A4ul3.303108
Paracentral lobuleL. A1/2/3ll2.58818
Temporal lobeInferior temporal gyrusL.A20iv2.97137
Parahippocampal gyrusL.TL2.39731
Parietal lobeInferior parietal lobuleL.A39c3.01126
R. A39rd2.674128
Insular lobeInsular gyrusR. dIa1.78960
Subcortical nucleiAmygdalaR. mAmyg0.5576
HippocampusL. cHipp3.79732

PBD-I: pediatric bipolar disorder type 1. L = left; R = right.

Table 2

Thirteen brain regions showing increased activations in PBD-I patients after overlapping the group activation map and Human Brainnetome Atlas.

LobeGyrusCytoarchitectonic namePeak T valueCluster size
Frontal lobeSuperior frontal gyrusL.A10m2.66784
Inferior frontal gyrusL.A45r1.48243
Orbital gyrusR. A12/47 l2.27264
Precentral gyrusL.A4hf3.994121
L.A4ul3.303108
Paracentral lobuleL. A1/2/3ll2.58818
Temporal lobeInferior temporal gyrusL.A20iv2.97137
Parahippocampal gyrusL.TL2.39731
Parietal lobeInferior parietal lobuleL.A39c3.01126
R. A39rd2.674128
Insular lobeInsular gyrusR. dIa1.78960
Subcortical nucleiAmygdalaR. mAmyg0.5576
HippocampusL. cHipp3.79732
LobeGyrusCytoarchitectonic namePeak T valueCluster size
Frontal lobeSuperior frontal gyrusL.A10m2.66784
Inferior frontal gyrusL.A45r1.48243
Orbital gyrusR. A12/47 l2.27264
Precentral gyrusL.A4hf3.994121
L.A4ul3.303108
Paracentral lobuleL. A1/2/3ll2.58818
Temporal lobeInferior temporal gyrusL.A20iv2.97137
Parahippocampal gyrusL.TL2.39731
Parietal lobeInferior parietal lobuleL.A39c3.01126
R. A39rd2.674128
Insular lobeInsular gyrusR. dIa1.78960
Subcortical nucleiAmygdalaR. mAmyg0.5576
HippocampusL. cHipp3.79732

PBD-I: pediatric bipolar disorder type 1. L = left; R = right.

HE results

As shown in Table 3, among the 13 brain regions that showed hyperactivation in the contrast of emotional versus neutral distractor during the Go-Nogo task, 8 brain regions demonstrated significantly higher HE values of fMRI signals in the PBD-I patients compared with that of HCs(FDR-corrected). The regions included left A10m_SFG, left A45r_IFG, right A12/47 l_ORG, right A4hf_PrG, left A20iv_ITG, left A39c_IPL, right A39rd_IPL, and right dIa_INS.

Table 3

PBD-I subjects exhibited significantly higher HE in eight brain regions than HCs.

Brain regionsBAPBD-I (n = 20)HCs (n = 17)p_FDRCohen’s d
L.A10m_SFG100.889 ± 0.0690.827 ± 0.0920.0250.762
L.A45r_IFG450.885 ± 0.0750.811 ± 0.1010.0250.831
R. A12/47 l_ORG470.889 ± 0.0860.813 ± 0.0840.0410.894
L.A4hf_ PrG60.889 ± 0.0700.829 ± 0.0810.0210.792
L.A4ul_ PrG60.822 ± 0.0940.805 ± 0.0960.3200.178
L. A1/2/3ll_ PCL0.844 ± 0.0850.811 ± 0.0880.1450.381
L.A20iv_ITG200.902 ± 0.0700.831 ± 0.0740.0180.985
L.TL_PhG300.876 ± 0.0990.830 ± 0.0730.0890.561
L.A39c_IPL390.896 ± 0.0610.821 ± 0.0880.0180.990
R. A39rd_IPL390.892 ± 0.0600.830 ± 0.0990.0250.757
R. dIa_INS480.884 ± 0.0690.830 ± 0.0750.0270.749
R. mAmyg_Amyg0.874 ± 0.0820.837 ± 0.0720.1010.479
L. cHipp_Hipp0.831 ± 0.0930.818 ± 0.1170.3630.123
Brain regionsBAPBD-I (n = 20)HCs (n = 17)p_FDRCohen’s d
L.A10m_SFG100.889 ± 0.0690.827 ± 0.0920.0250.762
L.A45r_IFG450.885 ± 0.0750.811 ± 0.1010.0250.831
R. A12/47 l_ORG470.889 ± 0.0860.813 ± 0.0840.0410.894
L.A4hf_ PrG60.889 ± 0.0700.829 ± 0.0810.0210.792
L.A4ul_ PrG60.822 ± 0.0940.805 ± 0.0960.3200.178
L. A1/2/3ll_ PCL0.844 ± 0.0850.811 ± 0.0880.1450.381
L.A20iv_ITG200.902 ± 0.0700.831 ± 0.0740.0180.985
L.TL_PhG300.876 ± 0.0990.830 ± 0.0730.0890.561
L.A39c_IPL390.896 ± 0.0610.821 ± 0.0880.0180.990
R. A39rd_IPL390.892 ± 0.0600.830 ± 0.0990.0250.757
R. dIa_INS480.884 ± 0.0690.830 ± 0.0750.0270.749
R. mAmyg_Amyg0.874 ± 0.0820.837 ± 0.0720.1010.479
L. cHipp_Hipp0.831 ± 0.0930.818 ± 0.1170.3630.123

Data are shown in mean ± standard deviation. L = left; R = right; SFG: superior frontal gyrus; IFG: inferior frontal gyrus; ORG: orbital gyrus; PrG: precentral gyrus; PCL: paracentral lobule; ITG: inferior temporal gyrus; PhG: parahippocampal gyrus; IPL: inferior parietal lobule; INS: insular gyrus; Amyg: amygdala; Hipp: hippocampus; BA: Brodmann area; p_FDR: FDR-corrected P-values

Table 3

PBD-I subjects exhibited significantly higher HE in eight brain regions than HCs.

Brain regionsBAPBD-I (n = 20)HCs (n = 17)p_FDRCohen’s d
L.A10m_SFG100.889 ± 0.0690.827 ± 0.0920.0250.762
L.A45r_IFG450.885 ± 0.0750.811 ± 0.1010.0250.831
R. A12/47 l_ORG470.889 ± 0.0860.813 ± 0.0840.0410.894
L.A4hf_ PrG60.889 ± 0.0700.829 ± 0.0810.0210.792
L.A4ul_ PrG60.822 ± 0.0940.805 ± 0.0960.3200.178
L. A1/2/3ll_ PCL0.844 ± 0.0850.811 ± 0.0880.1450.381
L.A20iv_ITG200.902 ± 0.0700.831 ± 0.0740.0180.985
L.TL_PhG300.876 ± 0.0990.830 ± 0.0730.0890.561
L.A39c_IPL390.896 ± 0.0610.821 ± 0.0880.0180.990
R. A39rd_IPL390.892 ± 0.0600.830 ± 0.0990.0250.757
R. dIa_INS480.884 ± 0.0690.830 ± 0.0750.0270.749
R. mAmyg_Amyg0.874 ± 0.0820.837 ± 0.0720.1010.479
L. cHipp_Hipp0.831 ± 0.0930.818 ± 0.1170.3630.123
Brain regionsBAPBD-I (n = 20)HCs (n = 17)p_FDRCohen’s d
L.A10m_SFG100.889 ± 0.0690.827 ± 0.0920.0250.762
L.A45r_IFG450.885 ± 0.0750.811 ± 0.1010.0250.831
R. A12/47 l_ORG470.889 ± 0.0860.813 ± 0.0840.0410.894
L.A4hf_ PrG60.889 ± 0.0700.829 ± 0.0810.0210.792
L.A4ul_ PrG60.822 ± 0.0940.805 ± 0.0960.3200.178
L. A1/2/3ll_ PCL0.844 ± 0.0850.811 ± 0.0880.1450.381
L.A20iv_ITG200.902 ± 0.0700.831 ± 0.0740.0180.985
L.TL_PhG300.876 ± 0.0990.830 ± 0.0730.0890.561
L.A39c_IPL390.896 ± 0.0610.821 ± 0.0880.0180.990
R. A39rd_IPL390.892 ± 0.0600.830 ± 0.0990.0250.757
R. dIa_INS480.884 ± 0.0690.830 ± 0.0750.0270.749
R. mAmyg_Amyg0.874 ± 0.0820.837 ± 0.0720.1010.479
L. cHipp_Hipp0.831 ± 0.0930.818 ± 0.1170.3630.123

Data are shown in mean ± standard deviation. L = left; R = right; SFG: superior frontal gyrus; IFG: inferior frontal gyrus; ORG: orbital gyrus; PrG: precentral gyrus; PCL: paracentral lobule; ITG: inferior temporal gyrus; PhG: parahippocampal gyrus; IPL: inferior parietal lobule; INS: insular gyrus; Amyg: amygdala; Hipp: hippocampus; BA: Brodmann area; p_FDR: FDR-corrected P-values

Correlation analysis

With age and gender regressed, significant correlation was found between HE of ORG and onset age in PBD-I subjects (Fig. 3).

Significant correlation was found between HE of frontal ORG(A12/47 l) and onset age in PBD-I subjects. Asterisk (*) represents the adjusted values with age and gender regressed. ORG: Orbital gyrus; HE: Hurst exponent; R: Right.
Fig 3

Significant correlation was found between HE of frontal ORG(A12/47 l) and onset age in PBD-I subjects. Asterisk (*) represents the adjusted values with age and gender regressed. ORG: Orbital gyrus; HE: Hurst exponent; R: Right.

Discussion

In the present study, an emotional stimulus in an emotional Go-Nogo task was conducted to evaluate emotional dysregulation process in PBD by investigating responses activation and alteration of fractal dynamic of BOLD signals. Increased neural responses and increased self-similarity were detected in PBD patients in predominantly cortical-limbic regions during the contrast of emotional distractor versus neutral distractor. No significant regional brain activation difference between PBD and HC subjects was observed in the contrast of Go-Nogo conditions versus resting state, and emontional versus non-emotional conditions. In addition, PBD adolescents exhibited no behavioral performance scarcity on the emotional Go-Nogo task.

It was reported that the response inhibition process in BD involves brain regions including the IFG, insula, dorsolateral prefrontal cortex, and parietal lobe to promote stimulus–response associations (Wessa et al. 2007). Consistent with this, the present study observed increased activation in a total of 13 brain regions including the above brain regions when comparing emotional and neutral distractors in PBD-I patients compared to healthy controls. The brain regions covered frontal lobe (SFG, IFG, ORG, PrG, PCL), temporal lobe (ITG, PhG), parietal lobe (IPL), insula, amygdala and hippocampus. The IFG extending into the right anterior insula was involved in refraining from action initiation (Cieslik et al. 2015; Fan et al. 2016). IPL is a component of frontoparietal network which is believed to guide attention to relevant stimuli, make trade-offs based on context, and maintain relevant information, so working memory and continuous attention are required (Seeley et al. 2007). Consistent with our results, increased brain activation to emotional versus neutral distractors was reported in orbitofrontal cortices and insula in euthymic bipolar patients. The increased orbitofrontal activation to happy distractors relative to neutral distractors was also found in manic BD patients (Elliott et al. 2000). The amygdala plays a key role in emotional information processing, and there is preliminary evidence of increased activation during mixed mood episodes in frontal and amygdala regions in BD (Chan et al. 2023). Resting glucose metabolism of euthymic old BD patients was found to be increased in the amygdala and PhG (Brooks 3rd et al. 2009). It was reported that manic patients showed elevated neuronal responses associated with emotional versus neutral targets in IFG (BA 47) and ITG (BA 20). On the other hand, most of the 13 brain regions discovered in the present study were key regions of ventrolateral circuit in corticlolimbic model related to BD (Brooks III and Vizueta 2014).The ventrolateral circuit consists of VLPFC, ventromedial striatum, globus pallidus, thalamus, anterior temporal cortex (BA20 and 38), and amygdala. It is believed to participate in responding to external emotional stimuli and regulating emotional expression. The overactivation of these brain regions discovered in the present study may cause damage to the ventrolateral circuit, thereby affecting PBD patients’ response to external stimuli. In a word, the overactivation of these regions may reflect a “compensatory” mechanism in which PBD patients need more efforts to engage the brain regions responsible for action inhibition and to achieve the same level of behavioral performance as healthy subjects.

Studies have verified that BOLD signals of the brain exhibit scale-free or fractal dynamics (Wei et al. 2013; Gentili et al. 2015). Fractal like dynamics represents the self-similarity of time series of BOLD signals. HE, an index between 0 and 1, is used to describe this self-similarity and a measure of brain activity complexity. Emotional regulation may be seen as a “buffering” system that allows for flexible responses to constantly changing conditions (Ortiz and Alda 2018). From a dynamic perspective, the characteristics that enable the system to adapt to these random changes are considered “complex” (Ortiz et al. 2021).The results of the present study indicated that, when comparing the emotional distractor versus neutral distractor during Go-Nogo task, the HE values of 8 among 13 overactivated brain regions were significantly higher in the patient group than in the normal group. The brain regions covering areas of frontal lobe (L.A10m_SFG, L.A45r_IFG, R. A12/47 l_ORG, L. A4hf_ PrG), temporal lobe (L.A20iv_ITG), parietal lobe (L.A39c_IPL, R. A39rd_IPL), and insula (R. dIa_INS)(Table 3). The enhanced HE means that the process during the Go-Nogo task have increased self-similarity and reduced complexity of BOLD signal in PBD-I patients. Furthermore, our correlation analysis results indicate that the decrease in complexity (increase in HE value) in ORG was negatively correlated with the age of onset of PBD. Consistent with our findings, reduced cortical complexity (gray matter fractal dimension) had been reported in frontal and temporal cortex that are heavily affected by BD subjects, which seem to reflect some neurodevelopmental antecedents (Squarcina et al. 2015; Nenadic et al. 2017). Furthermore, study of emotion and energy time series using paper-based visual analog scales exhibited a positive, short-term autocorrelation for both mood and energy series BD and HC groups with a more chaotic pattern in the time series for mood in HC than in BD participants (Ortiz et al. 2021), higher regularity of mood in patients with BD (Ortiz et al. 2015) and their unaffected first-degree relatives (Ortiz et al. 2019). Additionally, resting-state fMRI studies exhibited that brain dysfunction is associated with higher HE (lower fractal dimension) in individuals with autism spectrum disorder (Dona et al. 2017a) and mild traumatic brain injury (Dona et al. 2017b), with neurophysiological explanations as deficiency of adaptability and loss of functional components that accompany greater HE values. It has been suggested that healthy systems exhibited a more chaotic pattern, whereas damaged systems exhibit a more rigid pattern, lacking flexibility or resilience to adapt to the constantly changing demands of the environment (Ortiz et al. 2021). Brain regions shown increased HE in PBD subjects in the present study including dorsolateral prefrontal cortex, IFG, posterior parietal cortex, and insula are key regions of executive control network which involved in action initiation. The increase in BOLD signal HE in these brain regions indicates that the system responsible for responding to inhibition is not flexible enough for PBD-I patients to cope with the constantly changing needs of the environment when performing emotional Go-Nogo tasks. Thus, the BD may be believed as the result of a loss of complexity, rather than the loss of regularity (Ortiz et al. 2021).

In the present study, no significant difference was found in task performance between PBD-I patients and HCs. Similar results have been reported in euthymic and mix mood state BD subjects (Wessa et al. 2007; Morsel et al. 2017). It may suggest that BD patients have different neural mechanisms for achieving similar levels of performance as the HCs (Wessa et al. 2007). The compensatory over activation of key brain regions mentioned above, the increase of self-similarity and the decrease of complexity of neuronal activity process may contribute to the performance. On the other hand, this may also be related to the relatively simple nature of the Go Nogo task, while complex tasks, such as Stroop task, require more cognitive processes to be mobilized, which may lead to difficulty in response inhibition. The Go-Nogo task may directly enter the response inhibition process by load less other cognitive process. The advantage of utilizing simple task is that it displays changes of brain activity without being affected by differences in task performance, making an allowance for clear identification of neural abnormalities (Morsel et al. 2017).

Several limitations should be considered in the current study. Firstly, the sample size of each group in the study is relatively small. Studies involving more subjects are necessary to make more convincing conclusions in our future work. Secondly, medication exposure may influence the interpreting results of the present study. Though it has been reported that no significant effect of drug use in behavioral performance and neural activity between BD patients and HCs (Delvecchio et al. 2012; Jiang et al. 2020), more data from drug naïve BD patients will be extracted in our further research. Thirdly, the temporal sampling frequency (0.5 Hz) used to acquire the BOLD signal was relatively low, which may reduce the accuracy of fractal analysis. Higher sampling frequency may help to capture true dynamics of the neural activity process by means of fractal analysis. It is commanded that future research may use new MRI signal acquisition techniques to obtain BOLD signals and conduct fractal analysis and other nonlinear methods to mine more information.

Conclusions

The present study exhibited an overactivation, increased self-similarity and decreased complexity of BOLD signal in brain structures mediating response inhibition in PBD. The system that in charge of action inhibition is not flexible enough and requires greater effort to mobilize relevant brain regions to cope with cognitive tasks in PBD-I patients. HE may be a fractal biomarker of neural activity in PBD patients. Understanding the alteration in the fractal characteristics of the neural circuit related to response inhibition may provide useful clues for alleviating impulsive behaviors in BD, such as suicidal behavior and substance abuse.

Acknowledgments

We thank all the participants for their cooperation in the present study. All authors contributed to the article and approved the submitted version.

Author contributions

Yi-Bing Guo (Formal analysis, Methodology, Software, Validation, Visualization, Writing—original draft), Qian Xiao (Data curation, Funding acquisition, Resources), Qing Jiao (Conceptualization, Funding acquisition, Investigation, Project administration, Writing—review and editing), Weifang Cao (Funding acquisition, Writing—review and editing), Zhou Wu,Xiao-Tong Zhang (Investigation, Methodology, Software), Dong Cui (Supervision, Visualization), Guanghui Yu (Project administration, Supervision), Ruhai Dou (Investigation, Supervision), Linyan Su (Data curation, Resources), Guangming Lu (Data curation, Resources).

Funding

The present study was supported by the Funds for the National Natural Science Foundation of China (81371531 to Q.J.; 81901730 to W.-F.C.; 82201702 to Q.X.).

Conflict of interest statement: None declared.

Data availability

The data underlying this paper will be made available upon reasonable request the corresponding author.

References

Bassingthwaighte
JB
,
Raymond
GM
.
Evaluation of the dispersional analysis method for fractal time series
.
Ann Biomed Eng
.
1995
:
23
(
4
):
491
505
.

Bayani
A
,
Hadaeghi
F
,
Jafari
S
,
Murray
G
.
Critical slowing down as an early warning of transitions in episodes of bipolar disorder: a simulation study based on a computational model of circadian activity rhythms
.
Chronobiol Int
.
2017
:
34
(
2
):
235
245
.

Brooks
JO
III,
Vizueta
N
.
Diagnostic and clinical implications of functional neuroimaging in bipolar disorder
.
J Psychiatr Res
.
2014
:
57
:
12
25
.

Brooks
JO
III,
Hoblyn
JC
,
Woodard
SA
,
Rosen
AC
,
Ketter
TA
.
Corticolimbic metabolic dysregulation in euthymic older adults with bipolar disorder
.
J Psychiatr Res
.
2009
:
43
(
5
):
497
502
.

Campbell
OL
,
Weber
AM
.
Monofractal analysis of functional magnetic resonance imaging: an introductory review
.
Hum Brain Mapp
.
2022
:
43
(
8
):
2693
2706
.

Chan
CC
,
Alter
S
,
Hazlett
EA
,
Shafritz
KM
,
Yehuda
R
,
Goodman
M
,
Haznedar
MM
,
Szeszko
PR
.
Neural correlates of impulsivity in bipolar disorder: a systematic review and clinical implications
.
Neurosci Biobehav Rev
.
2023
:
147
:
105109
.

Chang
K
,
Barnea-Goraly
N
,
Karchemskiy
A
,
Simeonova
DI
,
Barnes
P
,
Ketter
T
,
Reiss
AL
.
Cortical magnetic resonance imaging findings in familial pediatric bipolar disorder
.
Biol Psychiatry
.
2005
:
58
(
3
):
197
203
.

Cieslik
EC
,
Mueller
VI
,
Eickhoff
CR
,
Langner
R
,
Eickhoff
SB
.
Three key regions for supervisory attentional control: evidence from neuroimaging meta-analyses
.
Neurosci Biobehav Rev
.
2015
:
48
:
22
34
.

Delvecchio
G
,
Fossati
P
,
Boyer
P
,
Brambilla
P
,
Falkai
P
,
Gruber
O
,
Hietala
J
,
Lawrie
SM
,
Martinot
JL
,
McIntosh
AM
, et al.
Common and distinct neural correlates of emotional processing in bipolar disorder and major depressive disorder: a voxel-based meta-analysis of functional magnetic resonance imaging studies
.
Eur Neuropsychopharmacol
.
2012
:
22
(
2
):
100
113
.

Delvecchio
G
,
Sugranyes
G
,
Frangou
S
.
Evidence of diagnostic specificity in the neural correlates of facial affect processing in bipolar disorder and schizophrenia: a meta-analysis of functional imaging studies
.
Psychol Med
.
2013
:
43
(
3
):
553
569
.

Dona
O
,
Hall
GB
,
Noseworthy
MD
.
Temporal fractal analysis of the rs-BOLD signal identifies brain abnormalities in autism spectrum disorder
.
PLoS One
.
2017a
:
12
(
12
):
e0190081
.

Dona
O
,
Noseworthy
MD
,
DeMatteo
C
,
Connolly
JF
.
Fractal analysis of brain blood oxygenation level dependent (BOLD) signals from children with mild traumatic brain injury (mTBI)
.
PLoS One
.
2017b
:
12
(
1
):
e0169647
.

Dong
J
,
Jing
B
,
Ma
X
,
Liu
H
,
Mo
X
,
Li
H
.
Hurst exponent analysis of resting-state fMRI signal complexity across the adult lifespan
.
Front Neurosci
.
2018
:
12
:
34
.

Eke
A
,
Herman
P
,
Kocsis
L
,
Kozak
LR
.
Fractal characterization of complexity in temporal physiological signals
.
Physiol Meas
.
2002
:
23
(
1
):
R1
R38
.

Ekman
P
,
Friesen
WV
.
Pictures of facial affect
.
Palo Alto, CA
:
Consulting Psychologists Press
;
1976
.

Elliott
R
,
Rubinsztein
JS
,
Sahakian
BJ
,
Dolan
RJ
.
Selective attention to emotional stimuli in a verbal go/no-go task: an fMRI study
.
Neuroreport
.
2000
:
11
(
8
):
1739
1744
.

Fan
L
,
Li
H
,
Zhuo
J
,
Zhang
Y
,
Wang
J
,
Chen
L
,
Yang
Z
,
Chu
C
,
Xie
S
,
Laird
AR
, et al.
The human Brainnetome atlas: a new brain atlas based on connectional architecture
.
Cereb Cortex
.
2016
:
26
(
8
):
3508
3526
.

Gentili
C
,
Vanello
N
,
Cristea
I
,
David
D
,
Ricciardi
E
,
Pietrini
P
.
Proneness to social anxiety modulates neural complexity in the absence of exposure: a resting state fMRI study using Hurst exponent
.
Psychiatry Res
.
2015
:
232
(
2
):
135
144
.

Goldberger
AL
,
Amaral
LA
,
Hausdorff
JM
,
Ivanov
P
,
Peng
CK
,
Stanley
HE
.
Fractal dynamics in physiology: alterations with disease and aging
.
Proc Natl Acad Sci U S A
.
2002
:
99
(
Suppl 1
):
2466
2472
.

Hafeman
DM
,
Bebko
G
,
Bertocci
MA
,
Fournier
JC
,
Bonar
L
,
Perlman
SB
,
Travis
M
,
Gill
MK
,
Diwadkar
VA
,
Sunshine
JL
, et al.
Abnormal deactivation of the inferior frontal gyrus during implicit emotion processing in youth with bipolar disorder: attenuated by medication
.
J Psychiatr Res
.
2014
:
58
:
129
136
.

Hummer
TA
,
Hulvershorn
LA
,
Karne
HS
,
Gunn
AD
,
Wang
Y
,
Anand
A
.
Emotional response inhibition in bipolar disorder: a functional magnetic resonance imaging study of trait- and state-related abnormalities
.
Biol Psychiatry
.
2013
:
73
(
2
):
136
143
.

Jiang
X
,
Fu
S
,
Yin
Z
,
Kang
J
,
Wang
X
,
Zhou
Y
,
Wei
S
,
Wu
F
,
Kong
L
,
Wang
F
, et al.
Common and distinct neural activities in frontoparietal network in first-episode bipolar disorder and major depressive disorder: preliminary findings from a follow-up resting state fMRI study
.
J Affect Disord
.
2020
:
260
:
653
659
.

Kaufman
J
,
Birmaher
B
,
Brent
D
,
Rao
U
,
Flynn
C
,
Moreci
P
,
Williamson
D
,
Ryan
N
.
Schedule for Affective Disorders and Schizophrenia for School-Age Children-Present and Lifetime Version (K-SADS-PL): initial reliability and validity data
.
J Am Acad Child Adolesc Psychiatry
.
1997
:
36
(
7
):
980
988
.

Knapen
SE
,
Li
P
,
Riemersma-van der Lek
RF
,
Verkooijen
S
,
Boks
MPM
,
Schoevers
RA
,
Scheer
F
,
Hu
K
.
Fractal biomarker of activity in patients with bipolar disorder
.
Psychol Med
.
2021
:
51
(
9
):
1562
1569
.

Long
Z
,
Li
J
,
Liao
H
,
Deng
L
,
Du
Y
,
Fan
J
,
Li
X
,
Miao
J
,
Qiu
S
,
Long
C
, et al.
A multi-modal and multi-atlas integrated framework for identification of mild cognitive impairment
.
Brain Sci
.
2022
:
12
(
6
):751.

Maxim
V
,
Sendur
L
,
Fadili
J
,
Suckling
J
,
Gould
R
,
Howard
R
,
Bullmore
E
.
Fractional Gaussian noise, functional MRI and Alzheimer's disease
.
NeuroImage
.
2005
:
25
(
1
):
141
158
.

Morsel
AM
,
Dhar
M
,
Hulstijn
W
,
Temmerman
A
,
Morrens
M
,
Sabbe
B
.
Inhibitory control in euthymic bipolar disorder: event related potentials during a Go/NoGo task
.
Clin Neurophysiol
.
2017
:
128
(
4
):
520
528
.

Nenadic
I
,
Yotter
RA
,
Dietzek
M
,
Langbein
K
,
Sauer
H
,
Gaser
C
.
Cortical complexity in bipolar disorder applying a spherical harmonics approach
.
Psychiatry Res Neuroimaging
.
2017
:
263
:
44
47
.

Ortiz
A
,
Alda
M
.
The perils of being too stable: mood regulation in bipolar disorder
.
J Psychiatry Neurosci
.
2018
:
43
(
6
):
363
365
.

Ortiz
A
,
Bradler
K
,
Garnham
J
,
Slaney
C
,
Alda
M
.
Nonlinear dynamics of mood regulation in bipolar disorder
.
Bipolar Disord
.
2015
:
17
(
2
):
139
149
.

Ortiz
A
,
Bradler
K
,
Garnham
J
,
Slaney
C
,
MacLean
S
,
Alda
M
.
Corrigendum to nonlinear dynamics of mood regulation in unaffected first-degree relatives of bipolar disorder patients [Journal of Affective Disorders 243 (2019) 274-279]
.
J Affect Disord
.
2019
:
245
:
16
.

Ortiz
A
,
Bradler
K
,
Mowete
M
,
MacLean
S
,
Garnham
J
,
Slaney
C
,
Mulsant
BH
,
Alda
M
.
The futility of long-term predictions in bipolar disorder: mood fluctuations are the result of deterministic chaotic processes
.
Int J Bipolar Disord
.
2021
:
9
(
1
):
30
.

Pavuluri
MN
,
O'Connor
MM
,
Harral
E
,
Sweeney
JA
.
Affective neural circuitry during facial emotion processing in pediatric bipolar disorder
.
Biol Psychiatry
.
2007
:
62
(
2
):
158
167
.

Poletti
S
,
Riberto
M
,
Vai
B
,
Ghiglino
D
,
Lorenzi
C
,
Vitali
A
,
Brioschi
S
,
Locatelli
C
,
Serretti
A
,
Colombo
C
, et al.
A glutamate transporter EAAT1 gene variant influences amygdala functional connectivity in bipolar disorder
.
J Mol Neurosci
.
2018
:
65
(
4
):
536
545
.

Robinson
JL
,
Monkul
ES
,
Tordesillas-Gutiérrez
D
,
Franklin
C
,
Bearden
CE
,
Fox
PT
,
Glahn
DC
.
Fronto-limbic circuitry in euthymic bipolar disorder: evidence for prefrontal hyperactivation
.
Psychiatry Res
.
2008
:
164
(
2
):
106
113
.

Rocca
CC
,
Heuvel
E
,
Caetano
SC
,
Lafer
B
.
Facial emotion recognition in bipolar disorder: a critical review
.
Braz J Psychiatry
.
2009
:
31
(
2
):
171
180
.

Seeley
WW
,
Menon
V
,
Schatzberg
AF
,
Keller
J
,
Glover
GH
,
Kenna
H
,
Reiss
AL
,
Greicius
MD
.
Dissociable intrinsic connectivity networks for salience processing and executive control
.
J Neurosci
.
2007
:
27
(
9
):
2349
2356
.

Sokunbi
MO
,
Gradin
VB
,
Waiter
GD
,
Cameron
GG
,
Ahearn
TS
,
Murray
AD
,
Steele
DJ
,
Staff RT
.
Nonlinear complexity analysis of brain FMRI signals in schizophrenia
.
PLoS One
.
2014
:
9
(
5
):
e95146
.

Squarcina
L
,
De Luca
A
,
Bellani
M
,
Brambilla
P
,
Turkheimer
FE
,
Bertoldo
A
.
Fractal analysis of MRI data for the characterization of patients with schizophrenia and bipolar disorder
.
Phys Med Biol
.
2015
:
60
(
4
):
1697
1716
.

Wei
M
,
Qin
J
,
Yan
R
,
Li
H
,
Yao
Z
,
Lu
Q
.
Identifying major depressive disorder using Hurst exponent of resting-state brain networks
.
Psychiatry Res
.
2013
:
214
(
3
):
306
312
.

Wessa
M
,
Houenou
J
,
Paillere-Martinot
ML
,
Berthoz
S
,
Artiges
E
,
Leboyer
M
,
Martinot
JL
.
Fronto-striatal overactivation in euthymic bipolar patients during an emotional Go/Nogo task
.
Am J Psychiatry
.
2007
:
164
(
4
):
638
646
.

Wink
AM
,
Bullmore
E
,
Barnes
A
,
Bernard
F
,
Suckling
J
.
Monofractal and multifractal dynamics of low frequency endogenous brain oscillations in functional MRI
.
Hum Brain Mapp
.
2008
:
29
(
7
):
791
801
.

Young
RC
,
Biggs
JT
,
Ziegler
VE
,
Meyer
DA
.
A rating scale for mania: reliability, validity and sensitivity
.
Br J Psychiatry
.
1978
:
133
(
5
):
429
435
.

Yurgelun-Todd
DA
,
Ross
AJ
.
Functional magnetic resonance imaging studies in bipolar disorder
.
CNS Spectr
.
2006
:
11
(
4
):
287
297
.

Zhang
L
,
Li
CT
,
Su
TP
,
Hu
XZ
,
Lanius
RA
,
Webster
MJ
,
Chung
MY
,
Chen
YS
,
Bai
YM
,
Barker
JL
, et al.
P11 expression and PET in bipolar disorders
.
J Psychiatr Res
.
2011
:
45
(
11
):
1426
1431
.

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (https://creativecommons.org/licenses/by-nc/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial re-use, please contact [email protected]