The connectivity of the human frontal pole cortex, and a theory of its involvement in exploit versus explore.

The frontal pole is implicated in humans in whether to exploit resources versus explore alternatives. Effective connectivity, functional connectivity, and tractography were measured between six human frontal pole regions and for comparison 13 dorsolateral and dorsal prefrontal cortex regions, and the 360 cortical regions in the Human Connectome Project Multi-modal-parcellation atlas in 171 HCP participants. The frontal pole regions have effective connectivity with Dorsolateral Prefrontal Cortex regions, the Dorsal Prefrontal Cortex, both implicated in working memory; and with the orbitofrontal and anterior cingulate cortex reward/non-reward system. There is also connectivity with temporal lobe, inferior parietal, and posterior cingulate regions. Given this new connectivity evidence, and evidence from activations and damage, it is proposed that the frontal pole cortex contains autoassociation attractor networks that are normally stable in a short-term memory state, and maintain stability in the other prefrontal networks during stable exploitation of goals and strategies. However, if an input from the orbitofrontal or anterior cingulate cortex that expected reward, non-reward, or punishment is received, this destabilizes the frontal pole and thereby other prefrontal networks to enable exploration of competing alternative goals and strategies. The frontal pole connectivity with reward systems may be key in exploit versus explore.

Dynamical differential covariance recovers directional network structure in multiscale neural systems.

Investigating neural interactions is essential to understanding the neural basis of behavior. Many statistical methods have been used for analyzing neural activity, but estimating the direction of network interactions correctly and efficiently remains a difficult problem. Here, we derive dynamical differential covariance (DDC), a method based on dynamical network models that detects directional interactions with low bias and high noise tolerance under nonstationarity conditions. Moreover, DDC scales well with the number of recording sites and the computation required is comparable to that needed for covariance. DDC was validated and compared favorably with other methods on networks with false positive motifs and multiscale neural simulations where the ground-truth connectivity was known. When applied to recordings of resting-state functional magnetic resonance imaging (rs-fMRI), DDC consistently detected regional interactions with strong structural connectivity in over 1,000 individual subjects obtained by diffusion MRI (dMRI). DDC is a promising family of methods for estimating connectivity that can be generalized to a wide range of dynamical models and recording techniques and to other applications where system identification is needed.

Discovering the effective connectome of the brain with dynamic Bayesian DAG learning.

Understanding the complex mechanisms of the brain can be unraveled by extracting the Dynamic Effective Connectome (DEC). Recently, score-based Directed Acyclic Graph (DAG) discovery methods have shown significant improvements in extracting the causal structure and inferring effective connectivity. However, learning DEC through these methods still faces two main challenges: one with the fundamental impotence of high-dimensional dynamic DAG discovery methods and the other with the low quality of fMRI data. In this paper, we introduce Bayesian Dynamic DAG learning with M-matrices Acyclicity characterization (BDyMA) method to address the challenges in discovering DEC. The presented dynamic DAG enables us to discover direct feedback loop edges as well. Leveraging an unconstrained framework in the BDyMA method leads to more accurate results in detecting high-dimensional networks, achieving sparser outcomes, making it particularly suitable for extracting DEC. Additionally, the score function of the BDyMA method allows the incorporation of prior knowledge into the process of dynamic causal discovery which further enhances the accuracy of results. Comprehensive simulations on synthetic data and experiments on Human Connectome Project (HCP) data demonstrate that our method can handle both of the two main challenges, yielding more accurate and reliable DEC compared to state-of-the-art and traditional methods. Additionally, we investigate the trustworthiness of DTI data as prior knowledge for DEC discovery and show the improvements in DEC discovery when the DTI data is incorporated into the process.

Negative association between neurovascular coupling and cortical gray matter volume during the lifespan.

Recent studies have established the moment-to-moment turnover of the blood-oxygen-level-dependent signal (TBOLD) at resting state as a key measure of local cortical brain function. Here, we sought to extend that line of research by evaluating TBOLD in 70 cortical areas with respect to corresponding brain volume, age, and sex across the lifespan in 1,344 healthy participants including 633 from the Human Connectome Project (HCP)-Development cohort (294 males and 339 females, age range 8-21 yr) and 711 healthy participants from HCP-Aging cohort (316 males and 395 females, 36-90 yr old). In both groups, we found that 1) TBOLD increased with age, 2) volume decreased with age, and 3) TBOLD and volume were highly significantly negatively correlated, independent of age. The inverse association between TBOLD and volume was documented in nearly all 70 brain areas and for both sexes, with slightly stronger associations documented for males. The strong correspondence between TBOLD and volume across age and sex suggests a common influence such as chronic neuroinflammation contributing to reduced cortical volume and increased TBOLD across the lifespan.NEW & NOTEWORTHY We report a significant negative association between resting functional magnetic resonance imaging (fMRI) blood-oxygen-level-dependent (BOLD) signal turnover (TBOLD) and cortical gray matter volume across the lifespan, such that TBOLD increased whereas volume decreased. We attribute this association to a hypothesized chronic, low-grade neuroinflammation, probably induced by various neurotropic pathogens, including human herpes viruses known to be dormant in the brain in a latent state and reactivated by stress, fever, and various environmental exposures, such as ultraviolet light.

Discovering individual fingerprints in resting-state functional connectivity using deep neural networks.

Non-negligible idiosyncrasy due to interindividual differences is an ongoing issue in resting-state functional MRI (rfMRI) analysis. We show that a deep neural network (DNN) can be employed for individual identification by learning important features from the time-varying functional connectivity (FC) of rfMRI in the Human Connectome Project. We employed the trained DNN to identify individuals from an independent dataset acquired at our institution. The results revealed that the DNN could successfully identify 300 individuals with an error rate of 2.9% using 15 s time-window and 870 individuals with an error rate of 6.7%. A trained DNN with nonlinear hidden layers led to the proposal of the "fingerprint of FC" (fpFC) as representative edges of individual FC. The fpFCs for individuals exhibited commonly important and individual-specific edges across time-window lengths (from 5 min to 15 s). Furthermore, the utility of our model for another group of subjects was validated, supporting the feasibility of our technique in the context of transfer learning. In conclusion, our study offers an insight into the discovery of the intrinsic mode of the human brain using whole-brain resting-state FC and DNNs.

Structural-functional connectivity mapping of the insular cortex: a combined data-driven and meta-analytic topic mapping.

In this study, we examined structural and functional profiles of the insular cortex and mapped associations with well-described functional networks throughout the brain using diffusion tensor imaging (DTI) and resting-state functional connectivity (RSFC) data. We used a data-driven method to independently estimate the structural-functional connectivity of the insular cortex. Data were obtained from the Human Connectome Project comprising 108 adult participants. Overall, we observed moderate to high associations between the structural and functional mapping scores of 3 different insular subregions: the posterior insula (associated with the sensorimotor network: RSFC, DTI = 50% and 72%, respectively), dorsal anterior insula (associated with ventral attention: RSFC, DTI = 83% and 83%, respectively), and ventral anterior insula (associated with the frontoparietal: RSFC, DTI = 42% and 89%, respectively). Further analyses utilized meta-analytic decoding maps to demonstrate specific cognitive and affective as well as gene expression profiles of the 3 subregions reflecting the core properties of the insular cortex. In summary, given the central role of the insular in the human brain, our results revealing correspondence between DTI and RSFC mappings provide a complementary approach and insight for clinical researchers to identify dysfunctional brain organization in various neurological disorders associated with insular pathology.

Prefrontal and somatosensory-motor cortex effective connectivity in humans.

Effective connectivity, functional connectivity, and tractography were measured between 57 cortical frontal and somatosensory regions and the 360 cortical regions in the Human Connectome Project (HCP) multimodal parcellation atlas for 171 HCP participants. A ventral somatosensory stream connects from 3b and 3a via 1 and 2 and then via opercular and frontal opercular regions to the insula, which then connects to inferior parietal PF regions. This stream is implicated in "what"-related somatosensory processing of objects and of the body and in combining with visual inputs in PF. A dorsal "action" somatosensory stream connects from 3b and 3a via 1 and 2 to parietal area 5 and then 7. Inferior prefrontal regions have connectivity with the inferior temporal visual cortex and orbitofrontal cortex, are implicated in working memory for "what" processing streams, and provide connectivity to language systems, including 44, 45, 47l, TPOJ1, and superior temporal visual area. The dorsolateral prefrontal cortex regions that include area 46 have connectivity with parietal area 7 and somatosensory inferior parietal regions and are implicated in working memory for actions and planning. The dorsal prefrontal regions, including 8Ad and 8Av, have connectivity with visual regions of the inferior parietal cortex, including PGs and PGi, and are implicated in visual and auditory top-down attention.

Altered neural associations with cognitive and emotional functions in cannabis dependence.

Negative emotional state has been found to correlate with poor cognitive performance in cannabis-dependent (CD) individuals, but not healthy controls (HCs). To examine the neural substrates underlying such unusual emotion-cognition coupling, we analyzed the behavioral and resting state fMRI data from the Human Connectome Project and found opposite brain-behavior associations in the CD and HC groups: (i) although the cognitive performance was positively correlated with the within-network functional connectivity strength and segregation (i.e. clustering coefficient and local efficiency) of the cognitive network in HCs, these correlations were inversed in CDs; (ii) although the cognitive performance was positively correlated with the within-network Granger effective connectivity strength and integration (i.e. characteristic path length) of the cognitive network in CDs, such associations were not significant in HCs. In addition, we also found that the effective connectivity strength within cognition network mediated the behavioral coupling between emotional state and cognitive performance. These results indicate a disorganization of the cognition network in CDs, and may help improve our understanding of substance use disorder.

Sex differences in default mode network connectivity in healthy aging adults.

Women show an increased lifetime risk of Alzheimer's disease (AD) compared with men. Characteristic brain connectivity changes, particularly within the default mode network (DMN), have been associated with both symptomatic and preclinical AD, but the impact of sex on DMN function throughout aging is poorly understood. We investigated sex differences in DMN connectivity over the lifespan in 595 cognitively healthy participants from the Human Connectome Project-Aging cohort. We used the intrinsic connectivity distribution (a robust voxel-based metric of functional connectivity) and a seed connectivity approach to determine sex differences within the DMN and between the DMN and whole brain. Compared with men, women demonstrated higher connectivity with age in posterior DMN nodes and lower connectivity in the medial prefrontal cortex. Differences were most prominent in the decades surrounding menopause. Seed-based analysis revealed higher connectivity in women from the posterior cingulate to angular gyrus, which correlated with neuropsychological measures of declarative memory, and hippocampus. Taken together, we show significant sex differences in DMN subnetworks over the lifespan, including patterns in aging women that resemble changes previously seen in preclinical AD. These findings highlight the importance of considering sex in neuroimaging studies of aging and neurodegeneration.

Neuronal activity induces symmetry breaking in neurodegenerative disease spreading.

Dynamical systems on networks typically involve several dynamical processes evolving at different timescales. For instance, in Alzheimer's disease, the spread of toxic protein throughout the brain not only disrupts neuronal activity but is also influenced by neuronal activity itself, establishing a feedback loop between the fast neuronal activity and the slow protein spreading. Motivated by the case of Alzheimer's disease, we study the multiple-timescale dynamics of a heterodimer spreading process on an adaptive network of Kuramoto oscillators. Using a minimal two-node model, we establish that heterogeneous oscillatory activity facilitates toxic outbreaks and induces symmetry breaking in the spreading patterns. We then extend the model formulation to larger networks and perform numerical simulations of the slow-fast dynamics on common network motifs and on the brain connectome. The simulations corroborate the findings from the minimal model, underscoring the significance of multiple-timescale dynamics in the modeling of neurodegenerative diseases.

Variability of the Surface Area of the V1, V2, and V3 Maps in a Large Sample of Human Observers.

How variable is the functionally defined structure of early visual areas in human cortex and how much variability is shared between twins? Here we quantify individual differences in the best understood functionally defined regions of cortex: V1, V2, V3. The Human Connectome Project 7T Retinotopy Dataset includes retinotopic measurements from 181 subjects (109 female, 72 male), including many twins. We trained four "anatomists" to manually define V1-V3 using retinotopic features. These definitions were more accurate than automated anatomical templates and showed that surface areas for these maps varied more than threefold across individuals. This threefold variation was little changed when normalizing visual area size by the surface area of the entire cerebral cortex. In addition to varying in size, we find that visual areas vary in how they sample the visual field. Specifically, the cortical magnification function differed substantially among individuals, with the relative amount of cortex devoted to central vision varying by more than a factor of 2. To complement the variability analysis, we examined the similarity of visual area size and structure across twins. Whereas the twin sample sizes are too small to make precise heritability estimates (50 monozygotic pairs, 34 dizygotic pairs), they nonetheless reveal high correlations, consistent with strong effects of the combination of shared genes and environment on visual area size. Collectively, these results provide the most comprehensive account of individual variability in visual area structure to date, and provide a robust population benchmark against which new individuals and developmental and clinical populations can be compared.SIGNIFICANCE STATEMENT Areas V1, V2, and V3 are among the best studied functionally defined regions in human cortex. Using the largest retinotopy dataset to date, we characterized the variability of these regions across individuals and the similarity between twin pairs. We find that the size of visual areas varies dramatically (up to 3.5×) across healthy young adults, far more than the variability of the cerebral cortex size as a whole. Much of this variability appears to arise from inherited factors, as we find very high correlations in visual area size between monozygotic twin pairs, and lower but still substantial correlations between dizygotic twin pairs. These results provide the most comprehensive assessment of how functionally defined visual cortex varies across the population to date.

Multi-modal multi-resolution atlas of the human neonatal cerebral cortex based on microstructural similarity.

The neonatal period is a critical window for the development of the human brain and may hold implications for the long-term development of cognition and disorders. Multi-modal connectome studies have revealed many important findings underlying the adult brain but related studies were rare in the early human brain. One potential challenge is the lack of an appropriate and unbiased parcellation that combines structural and functional information in this population. Using 348 multi-modal MRI datasets from the developing human connectome project, we found that the information fused from the structural, diffusion, and functional MRI was relatively stable across MRI features and showed high reproducibility at the group level. Therefore, we generated automated multi-resolution parcellations (300 - 500 parcels) based on the similarity across multi-modal features using a gradient-based parcellation algorithm. In addition, to acquire a parcellation with high interpretability, we provided a manually delineated parcellation (210 parcels), which was approximately symmetric, and the adjacent areas around each boundary were statistically different in terms of the integrated similarity metric and at least one kind of original features. Overall, the present study provided multi-resolution and neonate-specific parcellations of the cerebral cortex based on multi-modal MRI properties, which may facilitate future studies of the human connectome in the early development period.

PPA: Principal parcellation analysis for brain connectomes and multiple traits.

Our understanding of the structure of the brain and its relationships with human traits is largely determined by how we represent the structural connectome. Standard practice divides the brain into regions of interest (ROIs) and represents the connectome as an adjacency matrix having cells measuring connectivity between pairs of ROIs. Statistical analyses are then heavily driven by the (largely arbitrary) choice of ROIs. In this article, we propose a human trait prediction framework utilizing a tractography-based representation of the brain connectome, which clusters fiber endpoints to define a data-driven white matter parcellation targeted to explain variation among individuals and predict human traits. This leads to Principal Parcellation Analysis (PPA), representing individual brain connectomes by compositional vectors building on a basis system of fiber bundles that captures the connectivity at the population level. PPA eliminates the need to choose atlases and ROIs a priori, and provides a simpler, vector-valued representation that facilitates easier statistical analysis compared to the complex graph structures encountered in classical connectome analyses. We illustrate the proposed approach through applications to data from the Human Connectome Project (HCP) and show that PPA connectomes improve power in predicting human traits over state-of-the-art methods based on classical connectomes, while dramatically improving parsimony and maintaining interpretability. Our PPA package is publicly available on GitHub, and can be implemented routinely for diffusion image data.

Identifying and reverting the adverse effects of white matter hyperintensities on cortical surface analyses.

The Human Connectome Project (HCP)-style surface-based brain MRI analysis is a powerful technique that allows precise mapping of the cerebral cortex. However, the strength of its surface-based analysis has not yet been tested in the older population that often presents with white matter hyperintensities (WMHs) on T2-weighted (T2w) MRI (hypointensities on T1w MRI). We investigated T1-weighted (T1w) and T2w structural MRI in 43 healthy middle-aged to old participants. Juxtacortical WMHs were often misclassified by the default HCP pipeline as parts of the gray matter in T1w MRI, leading to incorrect estimation of the cortical surfaces and cortical metrics. To revert the adverse effects of juxtacortical WMHs, we incorporated the Brain Intensity AbNormality Classification Algorithm into the HCP pipeline (proposed pipeline). Blinded radiologists performed stereological quality control (QC) and found a decrease in the estimation errors in the proposed pipeline. The superior performance of the proposed pipeline was confirmed using an originally-developed automated surface QC based on a large database. Here we showed the detrimental effects of juxtacortical WMHs for estimating cortical surfaces and related metrics and proposed a possible solution for this problem. The present knowledge and methodology should help researchers identify adequate cortical surface biomarkers for aging and age-related neuropsychiatric disorders.

Changes of gray matter volumes of subcortical regions across the lifespan: a Human Connectome Project study.

We assessed changes in gray matter volume (GMV) of nine subcortical regions (accumbens, amygdala, brainstem, caudate, cerebellar cortex, pallidum, putamen, thalamus, and ventral diencephalon) across the lifespan in a large sample of participants in the Human Connectome Project (n = 2,458, 5-90 yr old, 1,113 males and 1,345 females). 3T MRI data were acquired using a harmonized protocol and were processed in an identical way for all brains. GMVs of individual regions were adjusted for estimated total intracranial volume and regressed against age. We found highly statistically significant changes in GMV with age (P < 0.001) that were distinct among areas and mostly consistent between sexes, as follows. 1) The GMVs of accumbens, caudate, putamen, and cerebellum decreased with age in a linear fashion. The rate of decrease was steeper in males than in females for all regions. 2) The GMVs of the amygdala, pallidum, thalamus, ventral diencephalon, and brainstem changed with age in a quadratic fashion, i.e., increasing first and decreasing afterward. The estimated age at the peak (vertex) of the parabola was 51.8 yr for the brainstem and 28.0-37.9 yr for the other regions. The peak occurred earlier in males than in females, by an average of 8 yr, with the exception of the brainstem, where the age at the peak was very similar in both sexes. These results confirm previous findings and offer new insights into region-specific age-related changes in subcortical brain GMVs.NEW & NOTEWORTHY We report mixed effects of age on subcortical grey matter volume (GMV) during lifespan (n = 2458, 5-90 yr old, 1113 male, 1345 female). Striatal and cerebellar GMVs decreased linearly with age, more steeply in males. In contrast, GMVs of the amygdala, pallidum, thalamus, ventral diencephalon, and brainstem changed in a quadratic fashion, increasing first and decreasing afterward, with males peaking earlier than females in all regions but the brainstem where they peaked at nearly the same time.

Changes of cortical gray matter volume during development: a Human Connectome Project study.

We assessed changes in gray matter volume of 35 cerebrocortical regions in a large sample of participants in the Human Connectome Project-Development (n = 649, 6-21 yr old, 299 males and 350 females). The same protocol for MRI data acquisition and processing was used for all brains. Volumes of individual areas were adjusted for estimated total intracranial volume and linearly regressed against age. We found changes of volume with age that were distinct among areas and consistent between sexes, as follows: 1) the overall cortical volume decreased significantly with age; 2) the volumes of 30/35 areas also decreased significantly with age; 3) the volumes of the hippocampal cortex (hippocampus, parahippocampal, and entorhinal) and that of pericalcarine cortex did not show significant age-related changes; and 4) the volume of the temporal pole increased significantly with age. The rates of volume reduction with age did not differ significantly between the two sexes, except for areas of the parietal lobe where males showed statistically significantly higher volume reduction with age than females. These results, obtained from a large sample of male and female participants, and acquired and processed in the same way, confirm previous findings, offer new insights into region-specific age-related changes in cortical brain volume, and are discussed in the context of the hypothesis that reduction in cortical volume may be partly due to a background, low-grade chronic neuroinflammation inflicted by common viruses residing latently in the brain, notably viruses of the human herpes family.NEW & NOTEWORTHY We report mixed effects of age on cortical gray matter volume during development in a large sample of 649 participants studied in an identical manner (6-21 yr old, 299 males, 350 females). Volumes of 30/35 cortical areas decreased with age, temporal pole increased, and pericalcarine and hippocampal cortex (hippocampus, parahippocampal, and entorhinal) did not change. These findings were very similar in both sexes and provide a solid base for assessing region-specific cortical changes during development.

Differential reduction of gray matter volume with age in 35 cortical areas in men (more) and women (less).

It is known that brain volume decreases with age. Here, we assessed the rate of this decrease in gray matter volume of 35 cortical regions in a large sample of healthy participants (n = 712, age range 36-90 yr) of the Human Connectome Project-Aging. We evaluated the difference in this rate between men (n = 316) and women (n = 396) and found that the volumes of cortical areas decreased by an average of 5.25%/decade, with the highest rate of decrease observed in the rostral anterior cingulate cortex (7.28%/decade). The rate of decrease was higher in men than in women in general and in 30/35 (85.7%) areas in particular, involving most prominently the cingulate lobe. These findings could serve as a normative reference for clinical conditions that manifest with abnormal brain atrophy.NEW & NOTEWORTHY This study showed an overall decrease of cortical gray matter with age but with different rates of volume reduction in different areas, with smaller decrease rates in women than in men. The highest volume reduction rate was observed for the rostral anterior cingulate cortex, an area linked to depression. These findings could serve as a normative reference for clinical conditions that manifest with abnormal brain atrophy.

Associations between cortisol stress responses and limbic volume and thickness in young adults: An exploratory study.

The investigation of the relationship between neural measures of limbic structures and hypothalamic pituitary adrenal axis responses to acute stress exposure in healthy young adults has so far focused in particular on task-based and resting state functional connectivity studies. Thus, the present study examined the association between limbic volume and thickness measures and acute cortisol responses to the psychosocial stress paradigm ScanSTRESS. Using Permutation Analysis of Linear Models controlling for sex, age and total brain volume, the associations between (sex-specific) cortisol increases and human connectome project style anatomical variables of limbic structures (i.e. volume and thickness) were investigated in 66 healthy and young (18-33 years) subjects (35 men, 31 women taking oral contraceptives). In addition, exploratory (sex-specific) bivariate correlations between cortisol increases and structural measures were conducted. The present data provide interesting new insights into the involvement of striato-limbic structures in psychosocial stress processing, suggesting that acute cortisol stress responses are also associated with mere structural measures of the human brain. Thus, our preliminary findings suggest that not only situation- and context-dependent reactions of the limbic system (i.e. blood oxygenation level-dependent reactions) are related to acute (sex-specific) cortisol stress responses but also basal and somewhat more constant structural measures. Our study hereby paves the way for further analyses in this context and highlights the relevance of the topic.

High-accuracy machine learning techniques for functional connectome fingerprinting and cognitive state decoding.

The human brain is a complex network comprised of functionally and anatomically interconnected brain regions. A growing number of studies have suggested that empirical estimates of brain networks may be useful for discovery of biomarkers of disease and cognitive state. A prerequisite for realizing this aim, however, is that brain networks also serve as reliable markers of an individual. Here, using Human Connectome Project data, we build upon recent studies examining brain-based fingerprints of individual subjects and cognitive states based on cognitively demanding tasks that assess, for example, working memory, theory of mind, and motor function. Our approach achieves accuracy of up to 99% for both identification of the subject of an fMRI scan, and for classification of the cognitive state of a previously unseen subject in a scan. More broadly, we explore the accuracy and reliability of five different machine learning techniques on subject fingerprinting and cognitive state decoding objectives, using functional connectivity data from fMRI scans of a high number of subjects (865) across a number of cognitive states (8). These results represent an advance on existing techniques for functional connectivity-based brain fingerprinting and state decoding. Additionally, 16 different functional connectome (FC) matrix construction pipelines are compared in order to characterize the effects of different aspects of the production of FCs on the accuracy of subject and task classification, and to identify possible confounds.

Explainable fMRI-based brain decoding via spatial temporal-pyramid graph convolutional network.

Brain decoding, aiming to identify the brain states using neural activity, is important for cognitive neuroscience and neural engineering. However, existing machine learning methods for fMRI-based brain decoding either suffer from low classification performance or poor explainability. Here, we address this issue by proposing a biologically inspired architecture, Spatial Temporal-pyramid Graph Convolutional Network (STpGCN), to capture the spatial-temporal graph representation of functional brain activities. By designing multi-scale spatial-temporal pathways and bottom-up pathways that mimic the information process and temporal integration in the brain, STpGCN is capable of explicitly utilizing the multi-scale temporal dependency of brain activities via graph, thereby achieving high brain decoding performance. Additionally, we propose a sensitivity analysis method called BrainNetX to better explain the decoding results by automatically annotating task-related brain regions from the brain-network standpoint. We conduct extensive experiments on fMRI data under 23 cognitive tasks from Human Connectome Project (HCP) S1200. The results show that STpGCN significantly improves brain-decoding performance compared to competing baseline models; BrainNetX successfully annotates task-relevant brain regions. Post hoc analysis based on these regions further validates that the hierarchical structure in STpGCN significantly contributes to the explainability, robustness and generalization of the model. Our methods not only provide insights into information representation in the brain under multiple cognitive tasks but also indicate a bright future for fMRI-based brain decoding.