Functional Connectivity and Complexity Analyses of Resting-State fMRI in Pre-Adolescents with ADHD

Attention deficit hyperactivity disorder (ADHD) has been characterized by impairments among distributed functional brain networks, e.g., the frontoparietal network (FPN), default mode network (DMN), and reward and motivation-related circuits (RMN). In the current study, we evaluated the complexity and functional connectivity (FC) of resting state fMRI (rsfMRI) in pre-adolescents with ADHD for pathology-relevant networks. We leveraged data from the Adolescent Brain and Cognitive Development (ABCD) Study. The final study sample included 63 children with ADHD and 92 healthy control children matched on age, sex, and pubertal development status. For selected regions in relevant networks, ANCOVA compared multiscale entropy (MSE) and FC between the groups. Finally, differences in the association between MSE and FC were evaluated. We found significantly reduced MSE along with increased FC within the FPN of pre-adolescents with ADHD compared to matched healthy controls. Significant partial correlations between MSE and FC emerged in fewer regions in the participants with ADHD than in the controls. The observation of reduced entropy is consistent with existing literature using rsfMRI and other neuroimaging modalities. The current findings of complexity and FC in ADHD support hypotheses of altered function of inhibitory control networks in ADHD.


Introduction
Attention deficit hyperactivity disorder (ADHD) is characterized by a persistent pattern of inattention and/or hyperactivity and impulsivity that causes impairment in functioning and development.ADHD occurs worldwide and is estimated to affect 5% of children and approximately 2.5% of adults (Polanczyk, de Lima, Horta, Biederman, & Rohde, 2007).
In ADHD, altered functional connectivity has been reported in the default mode network (DMN), frontoparietal network (FPN), and reward and motivation-related circuits (RMN; Castellanos & Aoki, 2016;Posner, Park, & Wang, 2014).Importantly, FCN findings in ADHD are inconsistent across different studies and age groups.Three recent meta-analyses summarized the rsfMRI functional connectivity findings of ADHD vs. healthy controls based on seed-based analyses (Cortese, Aoki, Itahashi, Castellanos, & Eickhoff, 2021;Gao et al., 2019;Sutcubasi et al., 2020).Gao et al. (2019) found ADHD was associated with hyperconnectivity between the seeds in DMN and the left superior temporal gyrus and right supramarginal gyrus extending to the right angular gyrus.ADHD also showed hypoconnectivity between the seeds in DMN to the left middle frontal gyrus and subcallosal cingulate cortex.With respect to the seeds in FPN, Gao et al. (2019) reported hyperconnectivity to the left orbital frontal cortex, caudate, putamen, and insula and hypoconnectivity to the precentral gyri among those with ADHD.For the seeds in RMN, hyperconnectivity to the dorsolateral prefrontal cortex was observed in ADHD.Sutcubasi et al. (2020) reported somewhat similar results in DMN and FPN but no agreement for RMN.In terms of the seeds in DMN, they found reduced connectivity to the dorsal posterior cingulate cortex but elevated connectivity to the dorsomedial prefrontal cortex in ADHD, while for the seeds in FPN, ADHD also showed elevated connectivity to the anterior prefrontal cortex.Cortese et al. (2021) did find no significant spatial convergence of ADHD-related hyperconnectivity or hypoconnectivity across studies.However, these inconclusive patterns of hyper and hypoconnectivity in ADHD also persist when focusing on children and adolescents only (Gao et al., 2019;Sutcubasi et al., 2020) but suggesting a more widespread connectivity alteration as compared to adult ADHD groups that could potentially be attributed to neurodevelopmental changes or delays in functional connectivity from childhood to adulthood (Fair et al., 2009).
Alongside functional connectivity analyses that assess the signal coherence between nodes, nonlinear analyses of neural signals from fMRI that characterize the signal complexity within a node have been proposed as measures for information processing capacity of brain areas and networks (McDonough & Nashiro, 2014;Smith, Yan, & Wang, 2014; D. J. J. Wang et al., 2018;Z. Wang, Li, Childress, & Detre, 2014), or indices of pathological brain function (Fernández, Gómez, Hornero, & López-Ibor, 2013;Sun et al., 2020;Takahashi, 2013).Sample entropy (SampEn; Richman & Moorman, 2000) has attracted considerable attention in complexity analysis because of its simplicity and reduced dependency on time series length compared to other forms of entropy such as approximation entropy (Pincus, 1991).SampEn measures the randomness and predictability of a stochastic process and generally increases with greater complexity.Multiscale entropy (MSE) was subsequently introduced to differentiate complex processes from random fluctuations more accurately, by calculating SampEn of a signal at multiple coarse-grained time scales (Costa, Goldberger, & Peng, 2002).Because MSE is evaluated across different time scales, it is also capable to identify frequency-dependent neuropathophysiological processes in different brain regions.Two studies investigated the complexity of fMRI data in adults with ADHD (Guan et al., 2023;Sokunbi et al., 2013) and both reported reduced complexity in the resting state in the frontal cortical areas, but no study as of yet has investigated fMRI complexity differences in children or pre-adolescents when ADHD is initially diagnosed.Additionally, no study has assessed the relation between functional connectivity and MSE (McDonough & Nashiro, 2014; D. J. J. Wang et al., 2018) in preadolescents with ADHD.
In the current study, we aimed to address this critical gap in the extant literature and evaluated the functional connectivity and complexity of rsfMRI in pre-adolescents with ADHD for the three main networks commonly associated with ADHD: FPN, DMN, and RMN.We combined MSE complexity analysis with functional connectivity analysis to elucidate the functional alterations present in ADHD on network nodes, at node-to-node edges, and on nodal complexity and node-to-node connectivity interaction levels.According to the literature, we hypothesized that compared to the healthy controls, pre-adolescents with ADHD had lower MSE in the FPN and altered functional connectivity in the FPN, DMN, and RMN.Due to the lack of previous studies, we did not postulate a hypothesis in terms of the difference between ADHD and control groups on the association between MSE and functional connectivity.
Head motion was monitored during scan acquisition using real-time procedures to adjust scanning procedures as necessary.
. For the current study, we used rsfMRI run1 and run2.Spatial and temporal image preprocessing was performed in the CONN toolbox (Conn: fMRI functional connectivity toolbox).Default parameter setup was used including motion-realignment, coregistration and normalization to MNI152 template space and smoothing with a 6mm FWHM Gaussian kernel.In the denoising process, white matter, CSF, and motion realignment parameters were included for temporal preprocessing and denoising of rsfMRI data.While motion-scrubbing was included in the preprocessing for functional connectivity analysis, this step was excluded for complexity analysis.

Complexity Analysis
Voxel-wise multiscale entropy (MSE) was computed using the in-house developed LOFT Complexity Toolbox (github.com/kayjann/complexity)for each preprocessed and denoised rsfMRI sequence for each participant.SampEn is defined as the natural logarithm of the conditional probability that a pattern length of m points will repeat itself, excluding self-matches, for m + 1 points within a tolerance of r in a time series of length N (Richman & Moorman, 2000).Multiple "coarse-grained" time series of lengths N/1, N/2, …, and N/a were formed by averaging consecutive data points of increasing length.SampEn for each coarse-grained time series is calculated and MSE is the average of SampEn over all the temporal scales.
In the current MSE calculation, we set the pattern length m = 2, the sensitivity threshold r = 0.3, and the number of temporal scales a = 15.MSE was computed for each voxel for each rsfMRI time series for each participant.The MSE map was then obtained by averaging MSE over the two rsfMRI runs for each voxel for each participant.Finally, mean MSE within each region-of-interest (ROI) of the three networks was calculated for each participant: ROIs were .selected in the anterior cingulate gyrus for both FPN and DMN, middle frontal gyrus and superior frontal gyrus for FPN, supramarginal gyrus, medial frontal cortex, and posterior cingulate gyrus for DMN, and frontal orbital cortex, caudate, and putamen for RMN based on the Harvard-Oxford atlas (Desikan et al., 2006).
A factorial 2 (Group: ADHD, Control) by 2 (Sex: Male, Female) analysis of covariance (ANCOVA) on the mean MSE in each ROI was conducted with PDS (self-reported) and scansite as covariates.We were only interested in the main effects of Group here and in this study's other ANCOVA analyses.To reduce skewness and kurtosis (Bliss, Greenwood, & White, 1956), PDS scores were Rankit-transformed and then z-scored.These values were used as the continuous covariates in all the ANCOVAs in the current study.A Benjamini-Hochberg correction (Benjamini & Hochberg, 1995) was performed for each network individually.The same ANCOVA was repeated for SampEn at each temporal scale.
In order to evaluate the association between complexity and ADHD behavior symptoms, the partial correlation between MSE and CBCL DSM5 ADHD scores controlling for sex, PDS zscore, and scan-site was computed for the whole sample.We also computed the partial correlation between SampEn at each temporal scale and the CBCL DSM5 ADHD scores controlling for the same covariates for the whole sample.Note that the CBCL DSM5 ADHD scores were Rankit-transformed and then z-scored in this partial correlation analysis.

Functional Connectivity Analysis
After data preprocessing and denoising, the first-level analyses on functional connectivity were conducted for each participant across both runs using the CONN toolbox default setup.Time series were averaged over voxels in each ROI, correlated between each pair of ROIs using Pearson's correlation analysis, and then Fisher's Z-transformed.The same ANCOVA including PDS z-score and scan-site as covariates was then performed for functional connectivity (Fisher's Z-transformed) using a seed-to-seed analysis.A Benjamini-Hochberg correction (Benjamini & Hochberg, 1995) was performed on the resultant main effects of Group for each network.
In order to evaluate the association between functional connectivity and ADHD behavior symptoms, the partial correlation between functional connectivity and CBCL DSM5 ADHD scores controlling for sex, PDS z-score, and scan-site was computed for the whole sample.Note again that the CBCL DSM5 ADHD scores were Rankit-transformed and then z-scored in this partial correlation analysis.

Analysis of the Association between Complexity and Functional Connectivity
To investigate the association between complexity and functional connectivity, we computed the partial correlation (i.e., Fisher's z score) between functional connectivity (between each seed ROI and each other ROI) and MSE in each seed ROI for each group by controlling for sex, PDS z-score, and scan-site.A group comparison (ADHD vs. controls) was conducted on the partial correlations for each (seed-MSE, seed-to-seed connectivity) linkage.The same analysis was repeated for the partial correlation between SampEn at each temporal scale and functional connectivity.

Follow-Up Analyses
In order to investigate the potential medication effects, a factorial 2 (Medication: medicated ADHD, non-medicated ADHD) by 2 (Sex: Male, Female) ANCOVA on the mean MSE in each ROI was conducted with PDS z-score and scan-site as covariates.The same ANCOVA was repeated for SampEn at each temporal scale and for functional connectivity at each seed-to-seed edge.

Complexity Analysis
The average MSE maps over the two rsfMRI runs for the ADHD and control groups are displayed in Figures 1(a) and (b).Qualitatively, it can be seen that complexity in the grey matter was higher than in the white matter.Moreover, the ADHD group appeared to have lower MSE than the controls, specifically in the frontal cortex.The ANCOVA revealed significant main effects of Group in multiple ROIs, i.  2).
Figure S1 in the supplementary material displays the t value of ADHD vs. control for mean MSE and SampEn at each temporal scale in each ROI which showed a significant main effect of Group (p < 0.05, uncorrected).It reveals that at scale 1 (frequency = 1.25Hz), every ROI showed a significant main effect of Group except the left anterior supramarginal gyrus.As .
the scale number increased, some significant main effects of Group disappeared although all the ROIs in FPN remained till scale = 9 (frequency = 0.111Hz).At scale = 15 (frequency = 0.083Hz), there was no significant difference between the ADHD group and the control group in terms of SampEn (p > 0.05, uncorrected).
We further tested whether there was an association between MSE and the CBCL DSM5 ADHD score but the partial correlation analysis for the whole sample taking into account sex, PDS z-score, and scan-site did not reach statistical significance in any ROI (p > 0.05).We also computed the same partial correlation between SampEn and the CBCL DSM5 ADHD score at each temporal scale and none was significant (p > 0.05).

Functional Connectivity Analysis
Figures 2(a

Analysis of the Association between Complexity and Functional Connectivity
Figures 3(a  anterior supramarginal gyrus-to-medial frontal cortex).ADHD showed significantly higher partial correlations relative to the controls within RMN (p < 0.05, uncorrected): The linkage of (right putamen, right putamen-to-right caudate) was significant at each temporal scale; Some significant within-RMN linkages emerged at the higher scale number (scale > 10), i.e., (left putamen, left putamen-to-right caudate), (right caudate, right caudate-to-right putamen), and (right caudate, right caudate-to-left putamen).ADHD showed a significantly higher partial correlation relative to the controls within FPN only at scale = 1 in the linkage of (right superior frontal gyrus, right superior frontal gyrus-to-anterior cingulate gyrus) (p < 0.05, uncorrected).In terms of cross-network linkages, we observed ADHD had a lower partial correlation between SampEn and functional connectivity than the controls in the linkages of (frontal medial cortex, frontal medial cortex-to-right putamen) and (frontal medial cortex, frontal medial cortex-to-left putamen) at almost every temporal scale up to scale = 13 (p < 0.05, uncorrected).We also observed ADHD had a higher partial correlation than the controls in the linkage of (right putamen, right putamen-to-right middle frontal gyrus) at scale = 11 & 12 (p < 0.05, uncorrected).

Follow-Up Analyses
In order to investigate potential medication-related effects, a factorial 2 (Medication: medicated ADHD, non-medicated ADHD) by 2 (Sex: Male, Female) ANCOVA on the mean MSE in each ROI was conducted with PDS z-score and scan-site as covariates.Please refer to Table S3 in the supplementary material for the list of prescriptions.The ADHD participants who

Discussion
In this study, we report two main observations: (1) reduced rsfMRI complexity in areas of the FPN in pre-adolescents with ADHD as compared to the matched controls, and (2) increased functional connectivity between anterior portions of the FPN in the ADHD group relative to the control group.We did not observe a significant group difference in terms of the association between complexity and functional connectivity in any interested linkage.
The FPN has been hypothesized to play a role in instantiating and flexibly modulating cognitive control (Marek & Dosenbach, 2018).Cognitive control is commonly compromised in many forms of psychopathology, many of which emerge during childhood and adolescence, while the control function is refined throughout development.Emerging literature has suggested an important role of the FPN in the pathophysiology of ADHD.Previous structural MRI studies consistently reported ADHD was associated with abnormal morphometry (Albajara Sáenz, Villemonteix, & Massat, 2019;Frodl & Skokauskas, 2012) and developmental trajectories (M.Shaw et al., 2012;P. Shaw et al., 2007) in the prefrontal and cingulate structures.Both qualitative review (Rubia, 2011) and meta-analyses (Cortese et al., 2012;Hart, Radua, Nakao, Mataix-Cols, & Rubia, 2013) on task-fMRI reported hypoactivation in the regions of FPN in ADHD across different tasks.rsfMRI meta-analyses also suggested functional alternation of FPN in ADHD (Cortese et al., 2021;Gao et al., 2019;Sutcubasi et al., 2020).Specifically, Gao et al.
(2019) reported children/adolescents with ADHD demonstrate hyperconnectivity between the seeds in the FPN and the dorsal anterior cingulate cortex, a finding that our results corroborate.
In addition to altered functional connectivity in the FPN, the current study revealed significantly different complexity in pre-adolescents with ADHD relative to the controls thus extending existing literature and providing additional evidence for the hypothesis of altered cognitive control in ADHD.
The observation of significant reductions in MSE in rsfMRI time series in preadolescents with ADHD was consistent with the findings in adults with ADHD using the same modality (Guan et al., 2023;Sokunbi et al., 2013).Guan et al. (2023) observed significant .
decreases of SampEn in multiple ROIs in every investigated brain network under every atlas in subjects with ADHD compared to controls.They also observed reduced SampEn of the gray matter at every temporal scale under every atlas in ADHD.Sokunbi et al. (2013) reported significant reductions of SampEn in adults with ADHD relative to the controls in the superior frontal and anterior cingulate gyrus which was replicated for pre-adolescents with ADHD in the current study (Figure 1(C)).Our results aligned with existing literature in resting-state using other neuroimaging modalities such as EEG and MEG in which complexity was found lower in the participants with ADHD than in the healthy controls (Angulo-Ruiz, Muñoz, Rodríguez-Martínez, Cabello-Navarro, & Gómez, 2022;Fernández et al., 2009;Gu, Liu, & Woltering, 2022;Gómez et al., 2013;Hu et al., 2021;Rezaeezadeh, Shamekhi, & Shamsi, 2020).
While multiple previous studies support the conceptualization of ADHD as a dysconnectivity syndrome, findings relating to functional connectivity under rsfMRI have been highly heterogeneous in terms of the affected networks and network subcomponents, as well as the direction of effects.The current study observed enhanced connectivity within FPN, i.e., the right superior frontal gyrus and left middle frontal gyrus using the anterior cingulate gyrus as the seed in pre-adolescents with ADHD (Figure 2).This result was consistent with previous literature that reported enhanced connectivity between the anterior cingulate gyrus and middle frontal gyrus in children/adolescents with ADHD compared to healthy controls (Kumar, Arya, & Agarwal, 2021;Lin, Kessler, Tseng, & Gau, 2021;Liu et al., 2023) except for one study reporting decreased connectivity (Posner et al., 2013).Interestingly this increased functional connectivity is predominantly between anterior nodes of the network and one might speculate that this represents a developmental stage of the network of increased short-range connectivity before long-range distributed network connections are established (Fair et al., 2009).
To date, several studies have investigated the differences between psychiatric and nonpsychiatric cohorts in terms of complexity or/and functional connectivity using functional MRI time series.However, while an association between MSE and functional connectivity in healthy cohorts has been previously reported (McDonough & Nashiro, 2014; D. J. J. Wang et al., 2018), as far as we have known, the current study is the first effort to conduct a comparison on the association between complexity and functional connectivity between healthy and diseased groups.While we found significant positive relationships between MSE and functional connectivity in the FPN and negative associations in the RMN in the healthy control group, such relationships were almost entirely absent in the group with ADHD (Figure 3).However, the direct statistical comparison did not reveal a significant group difference in the relation between MSE and functional connectivity between participants with ADHD and the controls.Future research in ADHD or other neuropsychiatric disorders following this approach might reveal disease and network-specific alterations in the association between nodal complexity and nodeto-node connectivity, two complementary but interrelated network properties.
Several caveats should be considered with respect to the current results.First, many ADHD participants were prescribed ADHD medications.In order to evaluate the effect of various comorbidities such as oppositional defiant disorder/conduct disorder that might confound the results.The current results removed this concern and further confirmed the finding of reduced complexity.Moreover, the lack of a unified understanding of functional connectivity in ADHD across literature may be due in part to different rates of comorbid conditions in those studies.The current finding of enhanced functional connectivity in ADHD relative to the controls in the FPN might potentially contribute to clarifying the mystery and provide a greater understanding of this field.Third, due to the limitation of networkwise Benjamini-Hochberg correction, the group difference in terms of cross-network functional connectivity and crossnetwork partial correlation between MSE and functional connectivity was absent.However, we showed significant cross-network group differences in Tables S1 and S2 and Figures S2 and S3 in the supplementary material.The reader should keep in mind that those effects were based on uncorrected statistics.
In conclusion, our analyses support the feasibility and utility of temporal signal entropy in investigating the changes in complexity of rsfMRI signals in pre-adolescents with ADHD when compared to healthy controls.We found reduced complexity in the anterior cingulate gyrus, bilateral middle frontal gyrus, and bilateral superior frontal gyrus in the FPN and greater functional connectivity in the right superior frontal gyrus and left middle frontal gyrus using the anterior cingulate gyrus as the seed in pre-adolescents with ADHD.We suggest that combining analyses of complexity and functional connectivity is a useful and easily obtainable measure to reveal changes in ADHD brain dynamics.

Figure
FigureS3in the supplementary material displays the significant partial correlations consumed at least one prescription were categorized as medicated ADHD.The result indicated a main effect of Medication in the left anterior supramarginal gyrus (DMN; F = 4.563, p < 0.05, ߟ ଶ = 0.186).However, it did not survive the networkwise Benjamini-Hochberg correction (false discovery rate = 0.05).Therefore, we concluded there was no significant difference in MSE between the medicated ADHD participants and the non-medicated ADHD participants in any ROI.The same ANCOVA was replicated for SampEn at each temporal scale.Main effects were found in multiple regions at multiple temporal scales (p < 0.05, uncorrected).More details can be accessed from the supplementary material.However, none of them survived the Benjamini-Hochberg correction (false discovery rate = 0.05).The same ANCOVA was replicated for functional connectivity at each seed-to-seed edge.It revealed significant main effects of Medication at two edges, i.e., left middle frontal gyrus-toright frontal orbital cortex (FPN-to-RMN; F = 4.509, p < 0.05, ߟ ଶ = 0.184) and medial frontal gyrus-to-posterior cingulate gyrus (DMN-to-DMN; F = 4.574, p < 0.05, ߟ ଶ = 0.186).However, none of them survived the Benjamini-Hochberg correction (false discovery rate = 0.05).

Figure 1 .
Figure 1.(a) MSE map for the ADHD group.(b) MSE map for the control group.(c) t value map

Figure 2 .
Figure 2. Functional connectivity maps (Fisher's z) for (a) the ADHD group and (b) the control

Figure 3 .
Figure 3. Maps of partial correlations (Fisher's z) between MSE and functional connectivity for

Table 1 .
An overview of sample characteristics delineated by group.

Table 1 :
Pubertal development status was measured by the self-reported pubertal development scale

Table 2 .
Brain regions displaying significant main effects of Group, obtained from the ANCOVAs for MSE and functional connectivity.
a According to the Harvard-Oxford Atlas