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.

Brain morphological variability between whites and African Americans: the importance of racial identity in brain imaging research

In a segregated society, marked by a historical background of inequalities, there is a consistent under-representation of ethnic and racial minorities in biomedical research, causing disparities in understanding genetic and acquired diseases as well as in the effectiveness of clinical treatments affecting different groups. The repeated inclusion of small and non-representative samples of the population in neuroimaging research has led to generalization bias in the morphological characterization of the human brain. A few brain morphometric studies between Whites and African Americans have reported differences in orbitofrontal volumetry and insula cortical thickness. Nevertheless, these studies are mostly conducted in small samples and populations with cognitive impairment. For this reason, this study aimed to identify brain morphological variability due to racial identity in representative samples. We hypothesized that, in neurotypical young adults, there are differences in brain morphometry between participants with distinct racial identities. We analyzed the Human Connectome Project (HCP) database to test this hypothesis. Brain volumetry, cortical thickness, and cortical surface area measures of participants identified as Whites (n = 338) or African Americans (n = 56) were analyzed. Non-parametrical permutation analysis of covariance between these racial identity groups adjusting for age, sex, education, and economic income was implemented. Results indicated volumetric differences in choroid plexus, supratentorial, white matter, and subcortical brain structures. Moreover, differences in cortical thickness and surface area in frontal, parietal, temporal, and occipital brain regions were identified between groups. In this regard, the inclusion of sub-representative minorities in neuroimaging research, such as African American persons, is fundamental for the comprehension of human brain morphometric diversity and to design personalized clinical brain treatments for this population.

The effect of processing pipelines, input images and age on automatic cortical morphology estimates

BACKGROUND AND OBJECTIVE

Magnetic resonance imaging of the brain allows to enrich the study of the relationship between cortical morphology, healthy ageing, diseases and cognition. Since manual segmentation of the cerebral cortex is time consuming and subjective, many software packages have been developed. FreeSurfer (FS) and Advanced Normalization Tools (ANTs) are the most used and allow as inputs a T1-weighted (T1w) image or its combination with a T2-weighted (T2w) image. In this study we evaluated the impact of different software and input images on cortical estimates. Additionally, we investigated whether the variation of the results depending on software and inputs is also influenced by age.

METHODS

For 240 healthy subjects, cortical thickness was computed with ANTs and FreeSurfer. Estimates were derived using both the T1w image and adding the T2w image. Significant effects due to software, input images and age range were investigated with ANOVA statistical analysis. Moreover, the accuracy of the cortical thickness estimates was assessed based on their age-prediction precision.

RESULTS

Using FreeSurfer and ANTs with T1w or T1w-T2w images resulted in significant differences in the cortical thickness estimates. These differences change with the age range of the subjects. Regardless of the images used, the more recent FS version tested exhibited the best performances in terms of age prediction.

CONCLUSIONS

Our study points out the importance of i) consistently processing data using the same tool; ii) considering the software, input images and the age range of the subjects when comparing multiple studies.

Investigating Habenula Functional Connectivity and Reward-Related Activity in Obesity Using Human Connectome Project Data

Introduction: The habenula, a brain region involved in aversion, might negatively modulate caloric intake. Functional magnetic resonance imaging (fMRI) studies reported associations between weight loss and habenula functional connectivity. However, whether habenula resting-state functional connectivity (rsFC) and reward-related activity are altered in obesity is yet unknown. Methods: Using data from the Human Connectome Project, we included 300 subjects with various body mass indexes (BMIs) and a healthy long-term blood glucose (hemoglobin A1c [HbA1c]). In addition, we investigated a potential BMI × HbA1c interaction in a separate cohort including subjects with prediabetes (n = 72). Habenula rsFC was assessed using a region of interest (ROI)-to-ROI analysis. Furthermore, a separate analysis using gambling task fMRI data focused on reward-related habenula activity. Results: We did not find an association between BMI and habenula rsFC for any of the ROIs. For the exploratory analysis of the BMI × HbA1c effect, a significant interaction effect was found for the habenula-ventral tegmental area (VTA) connection, but this did not survive multiple comparisons correction. Monetary punishment compared with reward activated the bilateral habenula in the BMI sample, but this activity was not associated with BMI. Discussion: In conclusion, we did not find evidence for an association between BMI and habenula rsFC or reward-related activity. However, there might be an interaction between BMI and HbA1c for the habenula-VTA rsFC, suggestive of a role of the habenula in glucose regulation. Future studies should focus on metabolic parameters in their experimental design to confirm our findings and explore the precise role of the habenula in metabolism.

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.

Variability of visual field maps in human early extrastriate cortex challenges the canonical model of organization of V2 and V3

Visual field maps in human early extrastriate areas (V2 and V3) are traditionally thought to form mirror-image representations which surround the primary visual cortex (V1). According to this scheme, V2 and V3 form nearly symmetrical halves with respect to the calcarine sulcus, with the dorsal halves representing lower contralateral quadrants, and the ventral halves representing upper contralateral quadrants. This arrangement is considered to be consistent across individuals, and thus predictable with reasonable accuracy using templates. However, data that deviate from this expected pattern have been observed, but mainly treated as artifactual. Here, we systematically investigate individual variability in the visual field maps of human early visual cortex using the 7T Human Connectome Project (HCP) retinotopy dataset. Our results demonstrate substantial and principled inter-individual variability. Visual field representation in the dorsal portions of V2 and V3 was more variable than in their ventral counterparts, including substantial departures from the expected mirror-symmetrical patterns. In addition, left hemisphere retinotopic maps were more variable than those in the right hemisphere. Surprisingly, only one-third of individuals had maps that conformed to the expected pattern in the left hemisphere. Visual field sign analysis further revealed that in many individuals the area conventionally identified as dorsal V3 shows a discontinuity in the mirror-image representation of the retina, associated with a Y-shaped lower vertical representation. Our findings challenge the current view that inter-individual variability in early extrastriate cortex is negligible, and that the dorsal portions of V2 and V3 are roughly mirror images of their ventral counterparts.

A unified global tractography framework for automatic visual pathway reconstruction

The human visual pathway starts from the retina, passes through the retinogeniculate visual pathway, the optic radiation, and finally connects to the primary visual cortex. Diffusion MRI tractography is the only technology that can noninvasively reconstruct the visual pathway. However, complete and accurate visual pathway reconstruction is challenging because of the skull base environment and complex fiber geometries. Specifically, the optic nerve within the complex skull base environment can cause abnormal diffusion signals. The crossing and fanning fibers at the optic chiasm, and a sharp turn of Meyer's loop at the optic radiation, contribute to complex fiber geometries of the visual pathway. A fiber trajectory distribution (FTD) function-based tractography method of our previous work and several high sensitivity tractography methods can reveal these complex fiber geometries, but are accompanied by false-positive fibers. Thus, the related studies of the visual pathway mostly applied the expert region of interest selection strategy. However, interobserver variability is an issue in reconstructing an accurate visual pathway. In this paper, we propose a unified global tractography framework to automatically reconstruct the visual pathway. We first extend the FTD function to a high-order streamline differential equation for global trajectory estimation. At the global level, the tractography process is simplified as the estimation of global trajectory distribution coefficients by minimizing the cost between trajectory distribution and the selected directions under the prior guidance by introducing the tractography template as anatomic priors. Furthermore, we use a deep learning-based method and tractography template prior information to automatically generate the mask for tractography. The experimental results demonstrate that our proposed method can successfully reconstruct the visual pathway with high accuracy.

Tensorial independent component analysis reveals social and reward networks associated with major depressive disorder

Major depressive disorder (MDD) has been associated with changes in functional brain connectivity. Yet, typical analyses of functional connectivity, such as spatial independent components analysis (ICA) for resting-state data, often ignore sources of between-subject variability, which may be crucial for identifying functional connectivity patterns associated with MDD. Typically, methods like spatial ICA will identify a single component to represent a network like the default mode network (DMN), even if groups within the data show differential DMN coactivation. To address this gap, this project applies a tensorial extension of ICA (tensorial ICA)-which explicitly incorporates between-subject variability-to identify functionally connected networks using functional MRI data from the Human Connectome Project (HCP). Data from the HCP included individuals with a diagnosis of MDD, a family history of MDD, and healthy controls performing a gambling and social cognition task. Based on evidence associating MDD with blunted neural activation to rewards and social stimuli, we predicted that tensorial ICA would identify networks associated with reduced spatiotemporal coherence and blunted social and reward-based network activity in MDD. Across both tasks, tensorial ICA identified three networks showing decreased coherence in MDD. All three networks included ventromedial prefrontal cortex, striatum, and cerebellum and showed different activation across the conditions of their respective tasks. However, MDD was only associated with differences in task-based activation in one network from the social task. Additionally, these results suggest that tensorial ICA could be a valuable tool for understanding clinical differences in relation to network activation and connectivity.

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.

A method for reconstruction of interpretable brain networks from transient synchronization in resting-state BOLD fluctuations

Resting-state (rs) fMRI has been widely used to examine brain-wide large-scale spatiotemporal architectures, known as resting-state networks (RSNs). Recent studies have focused on the temporally evolving characteristics of RSNs, but it is unclear what temporal characteristics are reflected in the networks. To address this issue, we devised a novel method for voxel-based visualization of spatiotemporal characteristics of rs-fMRI with a time scale of tens of seconds. We first extracted clusters of dominant activity-patterns using a region-of-interest approach and then used these temporal patterns of the clusters to obtain voxel-based activation patterns related to the clusters. We found that activation patterns related to the clusters temporally evolved with a characteristic temporal structure and showed mutual temporal alternations over minutes. The voxel-based representation allowed the decoding of activation patterns of the clusters in rs-fMRI using a meta-analysis of functional activations. The activation patterns of the clusters were correlated with behavioral measures. Taken together, our analysis highlights a novel approach to examine brain activity dynamics during rest.

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.

Interpretive JIVE: Connections with CCA and an application to brain connectivity

Joint and Individual Variation Explained (JIVE) is a model that decomposes multiple datasets obtained on the same subjects into shared structure, structure unique to each dataset, and noise. JIVE is an important tool for multimodal data integration in neuroimaging. The two most common algorithms are R.JIVE, an iterative approach, and AJIVE, which uses principal angle analysis. The joint structure in JIVE is defined by shared subspaces, but interpreting these subspaces can be challenging. In this paper, we reinterpret AJIVE as a canonical correlation analysis of principal component scores. This reformulation, which we call CJIVE, (1) provides an intuitive view of AJIVE; (2) uses a permutation test for the number of joint components; (3) can be used to predict subject scores for out-of-sample observations; and (4) is computationally fast. We conduct simulation studies that show CJIVE and AJIVE are accurate when the total signal ranks are correctly specified but, generally inaccurate when the total ranks are too large. CJIVE and AJIVE can still extract joint signal even when the joint signal variance is relatively small. JIVE methods are applied to integrate functional connectivity (resting-state fMRI) and structural connectivity (diffusion MRI) from the Human Connectome Project. Surprisingly, the edges with largest loadings in the joint component in functional connectivity do not coincide with the same edges in the structural connectivity, indicating more complex patterns than assumed in spatial priors. Using these loadings, we accurately predict joint subject scores in new participants. We also find joint scores are associated with fluid intelligence, highlighting the potential for JIVE to reveal important shared structure.

Predictive models for social functioning in healthy young adults: A machine learning study integrating neuroanatomical, cognitive, and behavioral data

Poor social functioning is an emerging public health problem associated with physical and mental health consequences. Developing prognostic tools is critical to identify individuals at risk for poor social functioning and guide interventions. We aimed to inform prediction models of social functioning by evaluating models relying on bio-behavioral data using machine learning. With data from the Human Connectome Project Healthy Young Adult sample (age 22-35, N = 1,101), we built Support Vector Regression models to estimate social functioning from variable sets of brain morphology to behavior with increasing complexity: 1) brain-only model, 2) brain-cognition model, 3) cognition-behavioral model, and 4) combined brain-cognition-behavioral model. Predictive accuracy of each model was assessed and the importance of individual variables for model performance was determined. The combined and cognition-behavioral models significantly predicted social functioning, whereas the brain-only and brain-cognition models did not. Negative affect, psychological wellbeing, extraversion, withdrawal, and cortical thickness of the rostral middle-frontal and superior-temporal regions were the most important predictors in the combined model. Results demonstrate that social functioning can be accurately predicted using machine learning methods. Behavioral markers may be more significant predictors of social functioning than brain measures for healthy young adults and may represent important leverage points for preventative intervention.

Inferring a Cognitive Architecture from Multitask Neuroimaging Data: A Data-Driven Test of the Common Model of Cognition Using Granger Causality

Cognitive architectures (i.e., theorized blueprints on the structure of the mind) can be used to make predictions about the effect of multiregion brain activity on the systems level. Recent work has connected one high-level cognitive architecture, known as the "Common Model of Cognition," to task-based functional MRI data with great success. That approach, however, was limited in that it was intrinsically top-down, and could thus only be compared with alternate architectures that the experimenter could contrive. In this paper, we propose a bottom-up method to infer a cognitive architecture directly from brain imaging data itself, overcoming this limitation. Specifically, Granger causality modeling was applied to the same task-based fMRI data to infer a network of causal connections between brain regions based on their functional connectivity. The resulting network shares many connections with those proposed by the Common Model of Cognition but also suggests important additions likely related to the role of episodic memory. This combined top-down and bottom-up modeling approach can be used to help formalize the computational instantiation of cognitive architectures and further refine a comprehensive theory of cognition.

A Multinomial Processing Tree Model of the 2-back Working Memory Task

The n -back task is a widely used behavioral task for measuring working memory and the ability to inhibit interfering information. We develop a novel model of the commonly used 2-back task using the cognitive psychometric framework provided by Multinomial Processing Trees. Our model involves three parameters: a memory parameter, corresponding to how well an individual encodes and updates sequence information about presented stimuli; a decision parameter corresponding to how well participants execute choices based on information stored in memory; and a base-rate parameter corresponding to bias for responding "yes" or "no". We test the parameter recovery properties of the model using existing 2-back experimental designs, and demonstrate the application of the model to two previous data sets: one from social psychology involving faces corresponding to different races (Stelter and Degner, British Journal of Psychology 109:777-798, 2018), and one from cognitive neuroscience involving more than 1000 participants from the Human Connectome Project (Van Essen et al., Neuroimage 80:62-79, 2013). We demonstrate that the model can be used to infer interpretable individual-level parameters. We develop a hierarchical extension of the model to test differences between stimulus conditions, comparing faces of different races, and comparing face to non-face stimuli. We also develop a multivariate regression extension to examine the relationship between the model parameters and individual performance on standardized cognitive measures including the List Sorting and Flanker tasks. We conclude by discussing how our model can be used to dissociate underlying cognitive processes such as encoding failures, inhibition failures, and binding failures.

Dynamical differential covariance recovers directional network structure in multiscale neural systems

We sense, move, and think by dynamical interactions between neurons. It is now possible to simultaneously record from many individual neurons and brain regions. Methods for analyzing these large-scale recordings are needed that can reveal how the patterns of activity give rise to behavior. We developed dynamical differential covariance (DDC), an efficient, intuitive, and robust way to analyze these recordings and validated it on simulations of model neural networks where the ground truth was known. It can estimate not only the presence of a connection but also which direction the information is flowing in a network between neurons or cortical areas. We applied DDC to recordings from functional magnetic resonance imaging in humans and confirmed predicted connectivity with direct measurements.

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.

Cerebello-basal Ganglia Networks and Cortical Network Global Efficiency

The cerebellum (CB) and basal ganglia (BG) each have topographically distinct functional subregions that are functionally and anatomically interconnected with cortical regions through discrete thalamic loops and with each other via disynaptic connections, with previous work detailing high levels of functional connectivity between these phylogenetically ancient regions. It was posited that this CB-BG network provides support for cortical systems processing, spanning cognitive, emotional, and motor domains, implying that subcortical network measures are strongly related to cortical network measures (Bostan & Strick, 2018); however, it is currently unknown how network measures within distinct CB-BG networks relate to cortical network measures. Here, 122 regions of interest comprising cognitive and motor CB-BG networks and 7 canonical cortical resting-state were used to investigate whether the integration (quantified using global efficiency, GE) of cognitive CB-BG network (CCBN) nodes and their segregation from motor CB-BG network (MCBN) nodes is related to cortical network GE and segregation in 233 non-related, right-handed participants (Human Connectome Project-1200). CCBN GE positively correlated with GE in the default mode, motor, and auditory networks and MCBN GE positively correlated with GE in all networks, except the default mode and emotional. MCBN segregation was related to motor network segregation. These findings highlight the CB-BG network's potential role in cortical networks associated with executive function, task switching, and verbal working memory. This work has implications for understanding cortical network organization and cortical-subcortical interactions in healthy adults and may help in determining biomarkers and deciphering subcortical differences seen in disease states.

A More Female-Characterized Resting-State Brain: Graph Similarity Analyses of Sex Influence on the Human Brain Intrinsic Functional Network

It remains unclear whether human species exhibits sexual dimorphism in brain activities, and how the dimorphisms associated with sex-characterized behaviors. Here, in a large dataset from Human Connectome Project, we investigated sex differences of resting-state network structure by using local and global network graph similarity analysis. The "typical male" and "typical female" resting-state networks were highly similar. However, we found significant inter-sex difference in all local brain networks compared with sex-label permutations. The global and many local network topologies showed significant higher intra-female similarity, while males' network topologies were more dissimilar to each other. Additionally, by using global graph similarity analysis, we found that female individuals whose brain network were more similar to the average pattern present lower social-related anger, lower social distress and better companionships, while similar effects were not detected for males. Our study confirms the existence of sex-related resting-state network topology. Female's intrinsic brain is closer to a typical pattern than male's, and they may more fulfill the "similarity breeds connection" principle in building social ties.

The positive-negative mode link between brain connectivity, demographics and behaviour: a pre-registered replication of Smith et al. (2015)

In mental health research, it has proven difficult to find measures of brain function that provide reliable indicators of mental health and well-being, including susceptibility to mental health disorders. Recently, a family of data-driven analyses have provided such reliable measures when applied to large, population-level datasets. In the current pre-registered replication study, we show that the canonical correlation analysis (CCA) methods previously developed using resting-state magnetic resonance imaging functional connectivity and subject measures (SMs) of cognition and behaviour from healthy adults are also effective in measuring well-being (a 'positive-negative axis') in an independent developmental dataset. Our replication was successful in two out of three of our pre-registered criteria, such that a primary CCA mode's weights displayed a significant positive relationship and explained a significant amount of variance in both functional connectivity and SMs. The only criterion that was not successful was that compared to other modes the magnitude of variance explained by the primary CCA mode was smaller than predicted, a result that could indicate a developmental trajectory of a primary mode. This replication establishes a signature neurotypical relationship between connectivity and phenotype, opening new avenues of research in neuroscience with clear clinical applications.

White matter association tracts underlying language and theory of mind: An investigation of 809 brains from the Human Connectome Project

Language and theory of mind (ToM) are the cognitive capacities that allow for the successful interpretation and expression of meaning. While functional MRI investigations are able to consistently localize language and ToM to specific cortical regions, diffusion MRI investigations point to an inconsistent and sometimes overlapping set of white matter tracts associated with these two cognitive domains. To further examine the white matter tracts that may underlie these domains, we use a two-tensor tractography method to investigate the white matter microstructure of 809 participants from the Human Connectome Project. 20 association white matter tracts (10 in each hemisphere) are uniquely identified by leveraging a neuroanatomist-curated automated white matter tract atlas. The fractional anisotropy (FA), mean diffusivity (MD), and number of streamlines (NoS) are measured for each white matter tract. Performance on neuropsychological assessments of semantic memory (NIH Toolbox Picture Vocabulary Test, TPVT) and emotion perception (Penn Emotion Recognition Test, PERT) are used to measure critical subcomponents of the language and ToM networks, respectively. Regression models are constructed to examine how structural measurements of left and right white matter tracts influence performance across these two assessments. We find that semantic memory performance is influenced by the number of streamlines of the left superior longitudinal fasciculus III (SLF-III), and emotion perception performance is influenced by the number of streamlines of the right SLF-III. Additionally, we find that performance on both semantic memory & emotion perception is influenced by the FA of the left arcuate fasciculus (AF). The results point to multiple, overlapping white matter tracts that underlie the cognitive domains of language and ToM. Results are discussed in terms of hemispheric dominance and concordance with prior investigations.

An extended Human Connectome Project multimodal parcellation atlas of the human cortex and subcortical areas

A modified and extended version, HCPex, is provided of the surface-based Human Connectome Project-MultiModal Parcellation atlas of human cortical areas (HCP-MMP v1.0, Glasser et al. 2016). The original atlas with 360 cortical areas has been modified in HCPex for ease of use with volumetric neuroimaging software, such as SPM, FSL, and MRIcroGL. HCPex is also an extended version of the original atlas in which 66 subcortical areas (33 in each hemisphere) have been added, including the amygdala, thalamus, putamen, caudate nucleus, nucleus accumbens, globus pallidus, mammillary bodies, septal nuclei and nucleus basalis. HCPex makes available the excellent parcellation of cortical areas in HCP-MMP v1.0 to users of volumetric software, such as SPM and FSL, as well as adding some subcortical regions, and providing labelled coronal views of the human brain.

Prediction of life satisfaction from resting-state functional connectome

BACKGROUND

Better life satisfaction (LS) is associated with better psychological and psychiatric outcomes. To the best of our knowledge, no studies have examined prediction models for LS.

METHODS

Using resting-state functional magnetic resonance imaging (R-fMRI) data from the Human Connectome Project (HCP) Young Adult S1200 dataset, we examined whether LS is predictable from intrinsic functional connectivity (iFC). All the HCP data were subdivided into either discovery (n = 100) or validation (n = 766) datasets. Using R-fMRI data in the discovery dataset, we computed a matrix of iFCs between brain regions. Ridge regression, in combination with principal component analysis and 10-fold cross-validation, was used to predict LS. Prediction performance was evaluated by comparing actual and predicted LS scores. The generalizability of the prediction model obtained from the discovery dataset was evaluated by applying this model to the validation dataset.

RESULTS

The model was able to successfully predict LS in the discovery dataset (r = 0.381, p < .001). The model was also able to successfully predict the degree of LS (r = 0.137, 5000-repetition permutation test p = .006) in the validation dataset, suggesting that our model is generalizable to the prediction of LS in young adults. iFCs stemming from visual, ventral attention, or limbic networks to other networks (such as the dorsal attention network and default mode network) were likely to contribute positively toward predicted LS scores. iFCs within ventral attention and limbic networks also positively contributed to predicting LS. On the other hand, iFCs stemming from the visual and cerebellar networks to other networks were likely to contribute negatively to the predicted LS scores.

CONCLUSION

The present findings suggest that LS is predictable from the iFCs. These results are an important step toward identifying the neural basis of life satisfaction.

Mixed-effects multilevel analysis followed by canonical correlation analysis is an effective fMRI tool for the investigation of idiosyncrasies

We report that regions-of-interest (ROIs) associated with idiosyncratic individual behavior can be identified from functional magnetic resonance imaging (fMRI) data using statistical approaches that explicitly model individual variability in neuronal activations, such as mixed-effects multilevel analysis (MEMA). We also show that the relationship between neuronal activation in fMRI and behavioral data can be modeled using canonical correlation analysis (CCA). A real-world dataset for the neuronal response to nicotine use was acquired using a custom-made MRI-compatible apparatus for the smoking of electronic cigarettes (e-cigarettes). Nineteen participants smoked e-cigarettes in an MRI scanner using the apparatus with two experimental conditions: e-cigarettes with nicotine (ECIG) and sham e-cigarettes without nicotine (SCIG) and subjective ratings were collected. The right insula was identified in the ECIG condition from the χ² -test of the MEMA but not from the t-test, and the corresponding activations were significantly associated with the similarity scores (r = -.52, p = .041, confidence interval [CI] = [-0.78, -0.17]) and the urge-to-smoke scores (r = .73, p <.001, CI = [0.52, 0.88]). From the contrast between the two conditions (i.e., ECIG > SCIG), the right orbitofrontal cortex was identified from the χ² -tests, and the corresponding neuronal activations showed a statistically meaningful association with similarity (r = -.58, p = .01, CI = [-0.84, -0.17]) and the urge to smoke (r = .34, p = .15, CI = [0.09, 0.56]). The validity of our analysis pipeline (i.e., MEMA followed by CCA) was further evaluated using the fMRI and behavioral data acquired from the working memory and gambling tasks available from the Human Connectome Project.

Delay discounting, cognitive ability, and personality: What matters?

Steep delay discounting is associated with problems such as addiction, obesity, and risky sexual behavior that are frequently described as reflecting impulsiveness and lack of self-control, but it may simply indicate poor cognitive functioning. The present investigation took advantage of the unique opportunity provided by the Human Connectome Project (N=1,206) to examine the relation between delay discounting and 11 cognitive tasks as well as the Big Five fundamental personality traits. With income level and education statistically controlled, discounting was correlated with only four of the 11 cognitive abilities evaluated, although the rs were all small (<.20). Importantly, the two discounting measures loaded on their own factor. Discounting was not correlated with Neuroticism or Conscientiousness, traits related to psychometric impulsiveness and self-control. These findings suggest that steep delay discounting is not simply an indicator of poor cognitive functioning or psychometric impulsiveness but an important individual difference characteristic in its own right.

Reduced Segregation Between Cognitive and Emotional Processes in Cannabis Dependence

Addiction is characterized by an erosion of cognitive control toward drug taking that is accentuated by negative emotional states. Here we tested the hypothesis that enhanced interference on cognitive control reflects a loss of segregation between cognition and emotion in addiction. We analyzed Human Connectome Project data from 1206 young adults, including 89 with cannabis dependence (CD). Two composite factors, one for cognition and one for emotion, were derived using principal component (PC) analyses. Component scores for these PCs were significantly associated in the CD group, such that negative emotionality correlated with poor cognition. However, the corresponding component scores were uncorrelated in matched controls and nondependent recreational cannabis users (n = 87). In CD, but not controls or recreational users, functional magnetic resonance imaging activations to emotional stimuli (angry/fearful faces > shapes) correlated with activations to cognitive demand (working memory; 2-back > 0-back). Canonical correlation analyses linked individual differences in cognitive and emotional component scores with brain activations. In CD, there was substantial overlap between cognitive and emotional brain-behavior associations, but in controls, associations were more restricted to the cognitive domain. These findings support our hypothesis of impaired segregation between cognitive and emotional processes in CD that might contribute to poor cognitive control under conditions of increased emotional demand.

Delineating the Decussating Dentato-rubro-thalamic Tract and Its Connections in Humans Using Diffusion Spectrum Imaging Techniques

The objective of this study was to identify the decussating dentato-rubro-thalamic tract (d-DRTT) and its afferent and efferent connections in healthy humans using diffusion spectrum imaging (DSI) techniques. In the present study, the trajectory and lateralization of the d-DRTT was explored using data from subjects in the Massachusetts General Hospital-Human Connectome Project adult diffusion dataset. The afferent and efferent networks that compose the cerebello-thalamo-cerebral pathways were also reconstructed. Correlation analysis was performed to identify interrelationships between subdivisions of the cerebello-dentato-rubro-thalamic and thalamo-cerebral connections. The d-DRTT was visualized bilaterally in 28 subjects. According to a normalized quantitative anisotropy and lateralization index evaluation, the left and right d-DRTT were relatively symmetric. Afferent regions were found mainly in the posterior cerebellum, especially the entire lobule VII (crus I, II and VIIb). Efferent fibers mainly are projected to the contralateral frontal cortex, including the motor and nonmotor regions. Correlations between cerebello-thalamic connections and thalamo-cerebral connections were positive, including the lobule VIIa (crus I and II) to the medial prefrontal cortex (MPFC) and the dorsolateral prefrontal cortex and lobules VI, VIIb, VIII, and IX, to the MPFC and motor and premotor areas. These results provide DSI-based tratographic evidence showing segregated and parallel cerebellar outputs to cerebral regions. The posterior cerebellum may play an important role in supporting and handling cognitive activities through d-DRTT. Future studies will allow for a more comprehensive understanding of cerebello-cerebral connections.

Boost in Test-Retest Reliability in Resting State fMRI with Predictive Modeling

Recent studies found low test-retest reliability in functional magnetic resonance imaging (fMRI), raising serious concerns among researchers, but these studies mostly focused on the reliability of individual fMRI features (e.g., individual connections in resting state connectivity maps). Meanwhile, neuroimaging researchers increasingly employ multivariate predictive models that aggregate information across a large number of features to predict outcomes of interest, but the test-retest reliability of predicted outcomes of these models has not previously been systematically studied. Here we apply 10 predictive modeling methods to resting state connectivity maps from the Human Connectome Project dataset to predict 61 outcome variables. Compared with mean reliability of individual resting state connections, we find mean reliability of the predicted outcomes of predictive models is substantially higher for all 10 modeling methods assessed. Moreover, improvement was consistently observed across all scanning and processing choices (i.e., scan lengths, censoring thresholds, volume- vs. surface-based processing). For the most reliable methods, the reliability of predicted outcomes was mostly, though not exclusively, in the "good" range (above 0.60). Finally, we identified three mechanisms that help to explain why predicted outcomes of predictive models have higher reliability than individual imaging features. We conclude that researchers can potentially achieve higher test-retest reliability by making greater use of predictive models.

The morphometry of left cuneus mediating the genetic regulation on working memory

Working memory is a basic human cognitive function. However, the genetic signatures and their biological pathway remain poorly understood. In the present study, we tried to clarify this issue by exploring the potential associations and pathways among genetic variants, brain morphometry and working memory performance. We first carried out association analyses between 2‐back accuracy and 212 image‐derived phenotypes from 1141 Human Connectome Project (HCP) subjects using a linear mixed model (LMM). We found a significantly positive correlation between the left cuneus volume and 2‐back accuracy (T = 3.615, p = 3.150e−4, Cohen's d = 0.226, corrected using family‐wise error [FWE] method). Based on the LMM‐based genome‐wide association study (GWAS) on the HCP dataset and UK Biobank 33 k GWAS summary statistics, we identified eight independent single nucleotide polymorphisms (SNPs) that were reliably associated with left cuneus volume in both UKB and HCP dataset. Within the eight SNPs, we found a negative correlation between the rs76119478 polymorphism and 2‐back accuracy accuracy (T = −2.045, p = .041, Cohen's d = −0.129). Finally, an LMM‐based mediation analysis elucidated a significant effect of left cuneus volume in mediating rs76119478 polymorphism on the 2‐back accuracy (indirect effect = −0.007, 95% BCa CI = [−0.045, −0.003]). These results were also replicated in a subgroup of Caucasians in the HCP population. Further fine mapping demonstrated that rs76119478 maps on intergene CTD‐2315A10.2 adjacent to protein‐encoding gene DAAM1, and is significantly associated with L3HYPDH mRNA expression. Our study suggested this new variant rs76119478 may regulate the working memory through exerting influence on the left cuneus volume.

A Connectivity-Based Psychometric Prediction Framework for Brain-Behavior Relationship Studies

The recent availability of population-based studies with neuroimaging and behavioral measurements opens promising perspectives to investigate the relationships between interindividual variability in brain regions' connectivity and behavioral phenotypes. However, the multivariate nature of connectivity-based prediction model severely limits the insight into brain-behavior patterns for neuroscience. To address this issue, we propose a connectivity-based psychometric prediction framework based on individual regions' connectivity profiles. We first illustrate two main applications: 1) single brain region's predictive power for a range of psychometric variables and 2) single psychometric variable's predictive power variation across brain region. We compare the patterns of brain-behavior provided by these approaches to the brain-behavior relationships from activation approaches. Then, capitalizing on the increased transparency of our approach, we demonstrate how the influence of various data processing and analyses can directly influence the patterns of brain-behavior relationships, as well as the unique insight into brain-behavior relationships offered by this approach.

A Whole-Cortex Probabilistic Diffusion Tractography Connectome

The WU-Minn Human Connectome Project (HCP) is a publicly-available dataset containing state-of-the-art structural magnetic resonance imaging (MRI), functional MRI (fMRI), and diffusion MRI (dMRI) for over a thousand healthy subjects. While the planned scope of the HCP included an anatomic connectome, resting-state fMRI (rs-fMRI) forms the bulk of the HCP's current connectomic output. We address this by presenting a full-cortex connectome derived from probabilistic diffusion tractography and organized into the HCP-MMP1.0 atlas. Probabilistic methods and large sample sizes are preferable for whole-connectome mapping as they increase the fidelity of traced low-probability connections. We find that overall, connection strengths are lognormally distributed and decay exponentially with tract length, that connectivity reasonably matches macaque histologic tracing in homologous areas, that contralateral homologs and left-lateralized language areas are hyperconnected, and that hierarchical similarity influences connectivity. We compare the dMRI connectome to existing rs-fMRI and cortico-cortico-evoked potential connectivity matrices and find that it is more similar to the latter. This work helps fulfill the promise of the HCP and will make possible comparisons between the underlying structural connectome and functional connectomes of various modalities, brain states, and clinical conditions.

New Insights Into the Anatomy, Connectivity and Clinical Implications of the Middle Longitudinal Fasciculus

The middle longitudinal fascicle (MdLF) is a long, associative white matter tract connecting the superior temporal gyrus (STG) with the parietal and occipital lobe. Previous studies show different cortical terminations, and a possible segmentation pattern of the tract. In this study, we performed a post-mortem white matter dissection of 12 human hemispheres and an in vivo deterministic fiber tracking of 24 subjects acquired from the Human Connectome Project to establish whether a constant organization of fibers exists among the MdLF subcomponents and to acquire anatomical information on each subcomponent. Moreover, two clinical cases of brain tumors impinged on MdLF territories are reported to further discuss the anatomical results in light of previously published data on the functional involvement of this bundle. The main finding is that the MdLF is consistently organized into two layers: an antero-ventral segment (aMdLF) connecting the anterior STG (including temporal pole and planum polare) and the extrastriate lateral occipital cortex, and a posterior-dorsal segment (pMdLF) connecting the posterior STG, anterior transverse temporal gyrus and planum temporale with the superior parietal lobule and lateral occipital cortex. The anatomical connectivity pattern and quantitative differences between the MdLF subcomponents along with the clinical cases reported in this paper support the role of MdLF in high-order functions related to acoustic information. We suggest that pMdLF may contribute to the learning process associated with verbal-auditory stimuli, especially on left side, while aMdLF may play a role in processing/retrieving auditory information already consolidated within the temporal lobe.

Structural connectivity of the human massa intermedia: A probabilistic tractography study

The role of massa intermedia (MI) is poorly understood in humans. Recent studies suggest its presence may play a role in normal human neurocognitive function while prior studies have shown the absence of MI correlated with psychiatric disorders. There is growing evidence that MI is likely a midline white matter conduit, responsible for interhemispheric connectivity, similar to other midline commissures. MI presence was identified in an unrelated sample using the Human Connectome Project database. MI structural connectivity maps were created and gray matter target regions were identified using probabilistic tractography of the whole brain. Probabilistic tractography revealed an extensive network of connections between MI and limbic, frontal and temporal lobes as well as insula and pericalcarine cortices. Women compared to men had stronger connectivity via their MI. The presented results support the role of MI as a midline commissure with strong connectivity to the amygdala, hippocampus, and entorhinal cortex.

A Whole-Cortex Probabilistic Diffusion Tractography Connectome

The WU-Minn Human Connectome Project (HCP) is a publicly-available dataset containing state-of-the-art structural magnetic resonance imaging (MRI), functional MRI (fMRI), and diffusion MRI (dMRI) for over a thousand healthy subjects. While the planned scope of the HCP included an anatomic connectome, resting-state fMRI (rs-fMRI) forms the bulk of the HCP's current connectomic output. We address this by presenting a full-cortex connectome derived from probabilistic diffusion tractography and organized into the HCP-MMP1.0 atlas. Probabilistic methods and large sample sizes are preferable for whole-connectome mapping as they increase the fidelity of traced low-probability connections. We find that overall, connection strengths are lognormally distributed and decay exponentially with tract length, that connectivity reasonably matches macaque histologic tracing in homologous areas, that contralateral homologs and left-lateralized language areas are hyperconnected, and that hierarchical similarity influences connectivity. We compare the dMRI connectome to existing rs-fMRI and cortico-cortico-evoked potential connectivity matrices and find that it is more similar to the latter. This work helps fulfill the promise of the HCP and will make possible comparisons between the underlying structural connectome and functional connectomes of various modalities, brain states, and clinical conditions.

fMRI volume classification using a 3D convolutional neural network robust to shifted and scaled neuronal activations

Deep-learning methods based on deep neural networks (DNNs) have recently been successfully utilized in the analysis of neuroimaging data. A convolutional neural network (CNN) is a type of DNN that employs a convolution kernel that covers a local area of the input sample and moves across the sample to provide a feature map for the subsequent layers. In our study, we hypothesized that a 3D-CNN model with down-sampling operations such as pooling and/or stride would have the ability to extract robust feature maps from the shifted and scaled neuronal activations in a single functional MRI (fMRI) volume for the classification of task information associated with that volume. Thus, the 3D-CNN model would be able to ameliorate the potential misalignment of neuronal activations and over-/under-activation in local brain regions caused by imperfections in spatial alignment algorithms, confounded by variability in blood-oxygenation-level-dependent (BOLD) responses across sessions and/or subjects. To this end, the fMRI volumes acquired from four sensorimotor tasks (left-hand clenching, right-hand clenching, auditory attention, and visual stimulation) were used as input for our 3D-CNN model to classify task information using a single fMRI volume. The classification performance of the 3D-CNN was systematically evaluated using fMRI volumes obtained from various minimal preprocessing scenarios applied to raw fMRI volumes that excluded spatial normalization to a template and those obtained from full preprocessing that included spatial normalization. Alternative classifier models such as the 1D fully connected DNN (1D-fcDNN) and support vector machine (SVM) were also used for comparison. The classification performance was also assessed for several k-fold cross-validation (CV) schemes, including leave-one-subject-out CV (LOOCV). Overall, the classification results of the 3D-CNN model were superior to that of the 1D-fcDNN and SVM models. When using the fully-processed fMRI volumes with LOOCV, the mean error rates (± the standard error of the mean) for the 3D-CNN, 1D-fcDNN, and SVM models were 2.1% (± 0.9), 3.1% (± 1.2), and 4.1% (± 1.5), respectively (p = 0.041 from a one-way ANOVA). The error rates for 3-fold CV were higher (2.4% ± 1.0, 4.2% ± 1.3, and 10.1% ± 2.0; p < 0.0003 from a one-way ANOVA). The mean error rates also increased considerably using the raw fMRI 3D volume data without preprocessing (26.2% for the 3D-CNN, 75.0% for the 1D-fcDNN, and 75.0% for the SVM). Furthermore, the ability of the pre-trained 3D-CNN model to handle shifted and scaled neuronal activations was demonstrated in an online scenario for five-class classification (i.e., four sensorimotor tasks and the resting state) using the real-time fMRI of three participants. The resulting classification accuracy was 78.5% (± 1.4), 26.7% (± 5.9), and 21.5% (± 3.1) for the 3D-CNN, 1D-fcDNN, and SVM models, respectively. The superior performance of the 3D-CNN compared to the 1D-fcDNN was verified by analyzing the resulting feature maps and convolution filters that handled the shifted and scaled neuronal activations and by utilizing an independent public dataset from the Human Connectome Project.

How Tasks Change Whole-Brain Functional Organization to Reveal Brain-Phenotype Relationships

Functional connectivity (FC) calculated from task fMRI data better reveals brain-phenotype relationships than rest-based FC, but how tasks have this effect is unknown. In over 700 individuals performing seven tasks, we use psychophysiological interaction (PPI) and predictive modeling analyses to demonstrate that task-induced changes in FC successfully predict phenotype, and these changes are not simply driven by task activation. Activation, however, is useful for prediction only if the in-scanner task is related to the predicted phenotype. To further characterize these predictive FC changes, we develop and apply an inter-subject PPI analysis. We find that moderate, but not high, task-induced consistency of the blood-oxygen-level-dependent (BOLD) signal across individuals is useful for prediction. Together, these findings demonstrate that in-scanner tasks have distributed, phenotypically relevant effects on brain functional organization, and they offer a framework to leverage both task activation and FC to reveal the neural bases of complex human traits, symptoms, and behaviors.

Scalable Bayesian matrix normal graphical models for brain functional networks

Recently, there has been an explosive growth in graphical modeling approaches for estimating brain functional networks. In a detailed study, we show that surprisingly, standard graphical modeling approaches for fMRI data may not yield accurate estimates of the brain network due to the inability to suitably account for temporal correlations. We propose a novel Bayesian matrix normal graphical model that jointly models the temporal covariance and the brain network under a separable structure for the covariance to obtain improved estimates. The approach is implemented via an efficient optimization algorithm that computes the maximum-a-posteriori network estimates having desirable theoretical properties and which is scalable to high dimensions. The proposed method leads to substantial gains in network estimation accuracy compared to standard brain network modeling approaches as illustrated via extensive simulations. We apply the method to resting state fMRI data from the Human Connectome Project involving a large number of time scans and brain regions, to study the relationships between fluid intelligence and functional connectivity, where it is not computationally feasible to apply existing matrix normal graphical models. Our proposed approach led to the detection of differences in connectivity between high and low fluid intelligence groups, whereas these differences were less pronounced or absent using the graphical lasso.

Decoding and mapping task states of the human brain via deep learning

Support vector machine (SVM)-based multivariate pattern analysis (MVPA) has delivered promising performance in decoding specific task states based on functional magnetic resonance imaging (fMRI) of the human brain. Conventionally, the SVM-MVPA requires careful feature selection/extraction according to expert knowledge. In this study, we propose a deep neural network (DNN) for directly decoding multiple brain task states from fMRI signals of the brain without any burden for feature handcrafts. We trained and tested the DNN classifier using task fMRI data from the Human Connectome Project's S1200 dataset (N = 1,034). In tests to verify its performance, the proposed classification method identified seven tasks with an average accuracy of 93.7%. We also showed the general applicability of the DNN for transfer learning to small datasets (N = 43), a situation encountered in typical neuroscience research. The proposed method achieved an average accuracy of 89.0 and 94.7% on a working memory task and a motor classification task, respectively, higher than the accuracy of 69.2 and 68.6% obtained by the SVM-MVPA. A network visualization analysis showed that the DNN automatically detected features from areas of the brain related to each task. Without incurring the burden of handcrafting the features, the proposed deep decoding method can classify brain task states highly accurately, and is a powerful tool for fMRI researchers.

Toward a Neural Model of the Openness-Psychoticism Dimension: Functional Connectivity in the Default and Frontoparietal Control Networks

Psychosis proneness has been linked to heightened Openness to Experience and to cognitive deficits. Openness and psychotic disorders are associated with the default and frontoparietal networks, and the latter network is also robustly associated with intelligence. We tested the hypothesis that functional connectivity of the default and frontoparietal networks is a neural correlate of the openness-psychoticism dimension. Participants in the Human Connectome Project (N = 1003) completed measures of psychoticism, openness, and intelligence. Resting state functional magnetic resonance imaging was used to identify intrinsic connectivity networks. Structural equation modeling revealed relations among personality, intelligence, and network coherence. Psychoticism, openness, and especially their shared variance were related positively to default network coherence and negatively to frontoparietal coherence. These associations remained after controlling for intelligence. Intelligence was positively related to frontoparietal coherence. Research suggests that psychoticism and openness are linked in part through their association with connectivity in networks involving experiential simulation and cognitive control. We propose a model of psychosis risk that highlights roles of the default and frontoparietal networks. Findings echo research on functional connectivity in psychosis patients, suggesting shared mechanisms across the personality-psychopathology continuum.

Partial Least Squares Regression Performs Well in MRI-Based Individualized Estimations

Estimation of individuals' cognitive, behavioral and demographic (CBD) variables based on MRI has attracted much research interest in the past decade, and effective machine learning techniques are of great importance for these estimations. Partial least squares regression (PLSR) is an attractive machine learning technique that can accommodate both single- and multi-label learning in a simple framework, while its potential for MRI-based estimations of CBD variables remains to be explored. In this study, we systemically investigated the performance of PLSR in MRI-based estimations of individuals' CBD variables, especially its performance in simultaneous estimation of multiple CBD variables (multi-label learning). We performed the study on the dataset included in the HCP S1200 release. Resting state functional connections (RSFCs) were used as features, and a total of 10 CBD variables (e.g., age, gender, grip strength, and picture vocabulary) were estimated. The results showed that PLSR performed well in both single- and multi-label learning. In fact, the present estimations were better than those reported in literatures, as indicated by stronger correlations between the estimated and actual CBD variables, as well as high gender classification accuracy (97.8% in this study). Moreover, the RSFCs that contributed to the estimations exhibited strong correlations with the CBD variable estimated, that is, PLSR algorithm automatically selected the RSFCs closely related to one CBD variable to establish predictive models for the variable. Besides, the estimation accuracies based on RSFCs among 100, 200, and 300 regions of interest (ROIs) were higher than those based on RSFCs among 15, 25, and 50 ROIs; the estimation accuracies based on RSFCs evaluated using partial correlation were higher than those based on RSFCs evaluated using full correlation. In addition to the aforementioned virtues, PLSR is efficient in model training and testing, and it is simple and easy to use. Therefore, PLSR can be a favorable choice for future MRI-based estimations of CBD variables.

Dynamic functioning of transient resting-state coactivation networks in the Human Connectome Project

Resting‐state analyses evaluating large‐scale brain networks have largely focused on static correlations in brain activity over extended time periods, however emerging approaches capture time‐varying or dynamic patterns of transient functional networks. In light of these new approaches, there is a need to classify common transient network states (TNS) in terms of their spatial and dynamic properties. To fill this gap, two independent resting state scans collected in 462 healthy adults from the Human Connectome Project were evaluated using coactivation pattern analysis to identify (eight) TNS that recurred across participants and over time. These TNS spatially overlapped with prototypical resting state networks, but also diverged in notable ways. In particular, analyses revealed three TNS that shared cortical midline overlap with the default mode network (DMN), but these “complex” DMN states also encompassed distinct regions that fall beyond the prototypical DMN, suggesting that the DMN defined using static methods may represent the average of distinct complex‐DMN states. Of note, dwell time was higher in “complex” DMN states, challenging the idea that the prototypical DMN, as a single unit, is the dominant resting‐state network as typically defined by static resting state methods. In comparing the two resting state scans, we also found high reliability in the spatial organization and dynamic activities of network states involving DMN or sensorimotor regions. Future work will determine whether these TNS defined by coactivation patterns are in other samples, and are linked to fundamental cognitive properties.

Associations Between Drinking and Cortical Thickness in Younger Adult Drinkers: Findings From the Human Connectome Project

BACKGROUND

Previous neuroimaging studies examining relations between alcohol misuse and cortical thickness have revealed that increased drinking quantity and alcohol-related problems are associated with thinner cortex. Although conflicting regional effects are often observed, associations are generally localized to frontal regions (e.g., dorsolateral prefrontal cortex [DLPFC], inferior frontal gyrus [IFG], and anterior cingulate cortex). Inconsistent findings may be attributed to methodological differences, modest sample sizes, and limited consideration of sex differences.

METHODS

This study examined neuroanatomical correlates of drinking quantity and heavy episodic drinking in a large sample of younger adults (N = 706; Mage  = 28.8; 51% female) using magnetic resonance imaging data from the Human Connectome Project. Exploratory analyses examined neuroanatomical correlates of executive function (flanker task) and working memory (list sorting).

RESULTS

Hierarchical linear regression models (controlling for age, sex, education, income, smoking, drug use, twin status, and intracranial volume) revealed significant inverse associations between drinks in past week and frequency of heavy drinking and cortical thickness in a majority of regions examined. The largest effect sizes were found for frontal regions (DLPFC, IFG, and the precentral gyrus). Follow-up regression models revealed that the left DLPFC was uniquely associated with both drinking variables. Sex differences were also observed, with significant effects largely specific to men.

CONCLUSIONS

This study adds to the understanding of brain correlates of alcohol use in a large, gender-balanced sample of younger adults. Although the cross-sectional methodology precludes causal inferences, these findings provide a foundation for rigorous hypothesis testing in future longitudinal investigations.

Resting-state fMRI effective connectivity between the bed nucleus of the stria terminalis and amygdala nuclei

The bed nucleus of the stria terminalis (BNST) and the laterobasal nucleus (LB), centromedial nucleus (CM), and superficial nucleus (SF) of the amygdala form an interconnected dynamical system, whose combined activity mediates a variety of behavioral and autonomic responses in reaction to homeostatic challenges. Although previous research provided deeper insight into the structural and functional connections between these nuclei, studies investigating their resting-state functional magnetic resonance imaging (fMRI) connectivity were solely based on undirected connectivity measures. Here, we used high-quality data of 391 subjects from the Human Connectome Project to estimate the effective connectivity (EC) between the BNST, the LB, CM, and SF through spectral dynamic causal modeling, the relation of the EC estimates with age and sex as well as their stability over time. Our results reveal a time-stable asymmetric EC structure with positive EC between all amygdala nuclei, which strongly inhibited the BNST while the BNST exerted positive influence onto all amygdala nuclei. Simulation of the impulse response of the estimated system showed that this EC structure shapes partially antagonistic (out of phase) activity flow between the BNST and amygdala nuclei. Moreover, the BNST-LB and BNST-CM EC parameters were less negative in males. In conclusion, our data points toward partially separated information processing between BNST and amygdala nuclei in the resting-state.

Visual Organization of the Default Network

The default network (DN) is a brain network with correlated activities spanning frontal, parietal, and temporal cortical lobes. The DN activates for high-level cognition tasks and deactivates when subjects are actively engaged in perceptual tasks. Despite numerous observations, the role of DN deactivation remains unclear. Using computational neuroimaging applied to a large dataset of the Human Connectome Project (HCP) and to two individual subjects scanned over many repeated runs, we demonstrate that the DN selectively deactivates as a function of the position of a visual stimulus. That is, we show that spatial vision is encoded within the DN by means of deactivation relative to baseline. Our results suggest that the DN functions as a set of high-level visual regions, opening up the possibility of using vision-science tools to understand its putative function in cognition and perception.

From Matrices to Knowledge: Using Semantic Networks to Annotate the Connectome

The connectome is regarded as the key to brain function in health and disease. Structural and functional neuroimaging enables us to measure brain connectivity in the living human brain. The field of connectomics describes the connectome as a mathematical graph with its connection strengths being represented by connectivity matrices. Graph theory algorithms are used to assess the integrity of the graph as a whole and to reveal brain network biomarkers for brain diseases; however, the faulty wiring of single connections or subnetworks as the structural correlate for neurological or mental diseases remains elusive. We describe a novel approach to represent the knowledge of human brain connectivity by a semantic network – a formalism frequently used in knowledge management to describe the semantic relations between objects. In our novel approach, objects are brain areas and connectivity is modeled as semantic relations among them. The semantic network turns the graph of the connectome into an explicit knowledge base about which brain areas are interconnected. Moreover, this approach can semantically enrich the measured connectivity of an individual subject by the semantic context from ontologies, brain atlases and molecular biological databases. Integrating all measurements and facts into one unified feature space enables cross-modal comparisons and analyses. We used a query mechanism for semantic networks to extract functional, structural and transcriptome networks. We found that in general higher structural and functional connectivity go along with a lower differential gene expression among connected brain areas; however, subcortical motor areas and limbic structures turned out to have a localized high differential gene expression while being strongly connected. In an additional explorative use case, we could show a localized high availability of fkbp5, gmeb1, and gmeb2 genes at a connection hub of temporo-limbic brain networks. Fkbp5 is known for having a role in stress-related psychiatric disorders, while gmeb1 and gmeb2 encode for modulator proteins of the glucocorticoid receptor, a key receptor in the hormonal stress system. Semantic networks tremendously ease working with multimodal neuroimaging and neurogenetics data and may reveal relevant coincidences between transcriptome and connectome networks.

Genetic influences on cortical myelination in the human brain

Cortical myelination, which is essential for interneuronal communication and neurodevelopment, has been reported to be under genetic control. However, the degree to which genes contribute to the variability of myelination, the pattern of genetic control, and how genes influence the organization of myelination are largely unknown. To answer these questions, the present study calculated heritability estimates for myelination of the cortical regions using the high quality structural magnetic resonance imaging (MRI) scans from the Human Connectome Project pedigree cohort (n = 873, 383/490 M/F, 22-36 years of age). Then, we used transcriptional profiles to evaluate the contribution of myelination-related genes (data from the Allen Human Brain Atlas) to explain interregional variations in cortical myelination. Our results showed that all the cortical areas were modestly to moderately influenced by genetic factors (h² = 29%-66%, all Ps < 0.05 after Bonferroni correction). The genetic control of cortical myelination showed bilateral symmetry and an anterior-to-posterior gradation. A bivariate model indicated that the regions are strongly genetically correlated with their homologs in the opposite cerebral hemisphere. A cross-modal analysis did not find a correlation between cortical myelination and the expression levels of myelination-related genes. This could have been due to the small number of samples with expression data in each cortical region. Overall, our findings suggest that cortical myelination is shaped by genetic factors and may be useful to bridge the underlying genetic variants and the cognitive functioning and related neuropsychiatric disorders.

The Relation Between White Matter Microstructure and Network Complexity: Implications for Processing Efficiency

Brain structure has been proposed to facilitate as well as constrain functional interactions within brain networks. Simulation models suggest that integrity of white matter (WM) microstructure should be positively related to the complexity of BOLD signal - a measure of network interactions. Using 121 young adults from the Human Connectome Project, we empirically tested whether greater WM integrity would be associated with greater complexity of the BOLD signal during rest via multiscale entropy. Multiscale entropy measures the lack of predictability within a given time series across varying time scales, thus being able to estimate fluctuating signal dynamics within brain networks. Using multivariate analysis techniques (Partial Least Squares), we found that greater WM integrity was associated with greater network complexity at fast time scales, but less network complexity at slower time scales. These findings implicate two separate pathways through which WM integrity affects brain function in the prefrontal cortex - an executive-prefrontal pathway and a perceptuo-occipital pathway. In two additional samples, the main patterns of WM and network complexity were replicated. These findings support simulation models of WM integrity and network complexity and provide new insights into brain structure-function relationships.

Behavioral and Health Correlates of Resting-State Metastability in the Human Connectome Project

Metastability is currently considered a fundamental property of the functional configuration of brain networks. The present study sought to generate a normative reference framework for the metastability of the major resting-state networks (RSNs) (resting-state metastability dataset) and discover their association with demographic, behavioral, physical and cognitive features (non-imaging dataset) from 818 participants of the Human Connectome Project. Using sparse canonical correlation analysis, we found that the metastability and non-imaging datasets showed significant but modest interdependency. Notable associations between the metastability variate and the non-imaging features were observed for higher-order cognitive ability and indicators of physical well-being. The intra-class correlation coefficient between the sibling pairs in the sample was very low which argues against a significant familial influence on RSN metastability.

High-resolution whole-brain diffusion MRI at 7T using radiofrequency parallel transmission

PURPOSE

Investigating the utility of RF parallel transmission (pTx) for Human Connectome Project (HCP)-style whole-brain diffusion MRI (dMRI) data at 7 Tesla (7T).

METHODS

Healthy subjects were scanned in pTx and single-transmit (1Tx) modes. Multiband (MB), single-spoke pTx pulses were designed to image sagittal slices. HCP-style dMRI data (i.e., 1.05-mm resolutions, MB2, b-values = 1000/2000 s/mm² , 286 images and 40-min scan) and data with higher accelerations (MB3 and MB4) were acquired with pTx.

RESULTS

pTx significantly improved flip-angle detected signal uniformity across the brain, yielding ∼19% increase in temporal SNR (tSNR) averaged over the brain relative to 1Tx. This allowed significantly enhanced estimation of multiple fiber orientations (with ∼21% decrease in dispersion) in HCP-style 7T dMRI datasets. Additionally, pTx pulses achieved substantially lower power deposition, permitting higher accelerations, enabling collection of the same data in 2/3 and 1/2 the scan time or of more data in the same scan time.

CONCLUSION

pTx provides a solution to two major limitations for slice-accelerated high-resolution whole-brain dMRI at 7T; it improves flip-angle uniformity, and enables higher slice acceleration relative to current state-of-the-art. As such, pTx provides significant advantages for rapid acquisition of high-quality, high-resolution truly whole-brain dMRI data.

Sensitivity analysis of human brain structural network construction

Network neuroscience leverages diffusion-weighted magnetic resonance imaging and tractography to quantify structural connectivity of the human brain. However, scientists and practitioners lack a clear understanding of the effects of varying tractography parameters on the constructed structural networks. With diffusion images from the Human Connectome Project (HCP), we characterize how structural networks are impacted by the spatial resolution of brain atlases, total number of tractography streamlines, and grey matter dilation with various graph metrics. We demonstrate how injudicious combinations of highly refined brain parcellations and low numbers of streamlines may inadvertently lead to disconnected network models with isolated nodes. Furthermore, we provide solutions to significantly reduce the likelihood of generating disconnected networks. In addition, for different tractography parameters, we investigate the distributions of values taken by various graph metrics across the population of HCP subjects. Analyzing the ranks of individual subjects within the graph metric distributions, we find that the ranks of individuals are affected differently by atlas scale changes. Our work serves as a guideline for researchers to optimize the selection of tractography parameters and illustrates how biological characteristics of the brain derived in network neuroscience studies can be affected by the choice of atlas parcellation schemes.

Diffusion tractography has been proven to be a promising noninvasive technique to study the network properties of the human brain. However, how various tractography and network construction parameters affect network properties has not been studied using a large cohort of high-quality data. We utilize data provided by the Human Connectome Project to characterize the changes to network properties induced by varying the brain parcellation atlas scales, the number of reconstructed tractography tracks, and the degree of grey matter dilation with graph metrics. We illustrate the importance of increasing the reconstructed track sampling rate when higher atlas scales are used. In addition to changing the raw values of graph metrics, we find that the ranks of individuals relative to the population metric distributions are altered. We further discuss how the dependency of graph metric ranks can affect the brain characteristics derived in group comparison studies using network neuroscience techniques.