Quantifying Individual Brain Connectivity with Functional Principal Component Analysis for Networks

Inter-regional correlation estimators for functional magnetic resonance imaging

Functional magnetic resonance imaging (fMRI) functional connectivity between brain regions is often computed using parcellations defined by functional or structural atlases. Typically, some kind of voxel averaging is performed to obtain a single temporal correlation estimate per region pair. However, several estimators can be defined for this task, with various assumptions and degrees of robustness to local noise, global noise, and region size. In this paper, we systematically present and study the properties of 9 different functional connectivity estimators taking into account the spatial structure of fMRI data, based on a simple fMRI data spatial model. These include 3 existing estimators and 6 novel estimators. We demonstrate the empirical properties of the estimators using synthetic, animal, and human data, in terms of graph structure, repeatability and reproducibility, discriminability, dependence on region size, as well as local and global noise robustness. We prove analytically the link between regional intra-correlation and inter-region correlation, and show that the choice of estimator has a strong influence on inter-correlation values. Some estimators, including the commonly used correlation of averages (ca), are positively biased, and have more dependence to region size and intra-correlation than robust alternatives, resulting in spatially-dependent bias. We define the new local correlation of averages estimator with better theoretical guarantees, lower bias, significantly lower dependence on region size (Spearman correlation 0.40 vs 0.55, paired t-test T=27.2, p=1.1e-47), at negligible cost to discriminative power, compared to the ca estimator. The difference in connectivity pattern between the estimators is not distributed uniformly throughout the brain, but rather shows a clear ventral-dorsal gradient, suggesting that region size and intra-correlation plays an important role in shaping functional networks defined using the ca estimator, and leading to non-trivial differences in their connectivity structure. We provide an open source R package and equivalent Python implementation to facilitate the use of the new estimators, together with preprocessed rat time-series.

Network evolution of regional brain volumes in young children reflects neurocognitive scores and mother's education

The maturation of regional brain volumes from birth to preadolescence is a critical developmental process that underlies emerging brain structural connectivity and function. Regulated by genes and environment, the coordinated growth of different brain regions plays an important role in cognitive development. Current knowledge about structural network evolution is limited, partly due to the sparse and irregular nature of most longitudinal neuroimaging data. In particular, it is unknown how factors such as mother's education or sex of the child impact the structural network evolution. To address this issue, we propose a method to construct evolving structural networks and study how the evolving connections among brain regions as reflected at the network level are related to maternal education and biological sex of the child and also how they are associated with cognitive development. Our methodology is based on applying local Fréchet regression to longitudinal neuroimaging data acquired from the RESONANCE cohort, a cohort of healthy children (245 females and 309 males) ranging in age from 9 weeks to 10 years. Our findings reveal that sustained highly coordinated volume growth across brain regions is associated with lower maternal education and lower cognitive development. This suggests that higher neurocognitive performance levels in children are associated with increased variability of regional growth patterns as children age.

The Task Pre-Configuration Is Associated With Cognitive Performance Evidence From the Brain Synchrony

Although many resting state and task state characteristics have been studied, it is still unclear how the brain network switches from the resting state during tasks. The current theory shows that the brain is a complex dynamic system and synchrony is defined to measure brain activity. The study compared the changes of synchrony between the resting state and different task states in healthy young participants (N = 954). It also examined the ability to switch from the resting state to the task-general architecture of synchrony. We found that the synchrony increased significantly during the tasks. And the results showed that the brain has a task-general architecture of synchrony during different tasks. The main feature of task-based reasoning is that the increase in synchrony of high-order cognitive networks is significant, while the increase in synchrony of sensorimotor networks is relatively low. In addition, the high synchrony of high-order cognitive networks in the resting state can promote task switching effectively and the pre-configured participants have better cognitive performance, which shows that spontaneous brain activity and cognitive ability are closely related. These results revealed changes in the brain network configuration for switching between the resting state and task state, highlighting the consistent changes in the brain network between different tasks. Also, there was an important relationship between the switching ability and the cognitive performance.

Exploration of Multiparameter Hematoma 3D Image Analysis for Predicting Outcome After Intracerebral Hemorrhage

BACKGROUND

Rapid diagnosis and proper management of intracerebral hemorrhage (ICH) play a crucial role in the outcome. Prediction of the outcome with a high degree of accuracy based on admission data including imaging information can potentially influence clinical decision-making practice.

METHODS

We conducted a retrospective multicenter study of consecutive ICH patients admitted between 2012-2017. Medical history, admission data, and initial head computed tomography (CT) scan were collected. CT scans were semiautomatically segmented for hematoma volume, hematoma density histograms, and sphericity index (SI). Discharge unfavorable outcomes were defined as death or severe disability (modified Rankin Scores 4-6). We compared (1) hematoma volume alone; (2) multiparameter imaging data including hematoma volume, location, density heterogeneity, SI, and midline shift; and (3) multiparameter imaging data with clinical information available on admission for ICH outcome prediction. Multivariate analysis and predictive modeling were used to determine the significance of hematoma characteristics on the outcome.

RESULTS

We included 430 subjects in this analysis. Models using automated hematoma segmentation showed incremental predictive accuracies for in-hospital mortality using hematoma volume only: area under the curve (AUC): 0.85 [0.76-0.93], multiparameter imaging data (hematoma volume, location, CT density, SI, and midline shift): AUC: 0.91 [0.86-0.97], and multiparameter imaging data plus clinical information on admission (Glasgow Coma Scale (GCS) score and age): AUC: 0.94 [0.89-0.99]. Similarly, severe disability predictive accuracy varied from AUC: 0.84 [0.76-0.93] for volume-only model to AUC: 0.88 [0.80-0.95] for imaging data models and AUC: 0.92 [0.86-0.98] for imaging plus clinical predictors.

CONCLUSIONS

Multiparameter models combining imaging and admission clinical data show high accuracy for predicting discharge unfavorable outcome after ICH.

Principal Component Analysis of Simultaneous PET-MRI Reveals Patterns of Bone-Cartilage Interactions in Osteoarthritis

BACKGROUND

Bone-cartilage interactions have been implicated in causing osteoarthritis (OA).

PURPOSE

To use [18 F]-NaF PET-MRI to 1) develop automatic image processing code in MatLab to create a model of bone-cartilage interactions and 2) find associations of bone-cartilage interactions with known manifestations of OA.

STUDY TYPE

Prospective study aimed to evaluate a data analysis method.

POPULATION

Twenty-nine patients with knee pain or joint stiffness.

FIELD STRENGTH/SEQUENCE

3T MRI (GE), 3D CUBE FSE, 3D combined T1 ρ/T2 MAPSS, [18F]-sodium fluoride, SIGNA TOF (OSEM).

ASSESSMENT

Correlation between MRI (cartilage) and PET (bone) quantitative parameters, bone-cartilage interactions model described by modes of variation as derived by principal component analysis (PCA), WORMS scoring on cartilage lesions, bone marrow abnormalities, subchondral cysts.

STATISTICAL TESTS

Linear regression, Pearson correlation.

RESULTS

Mode 1 was a positive predictor of the bone abnormality score (P = 0.0003, P = 0.001, P = 0.0007) and the cartilage lesion score (P = 0.03, P = 0.01, P = 0.02) in the femur, tibia, and patella, respectively. For the cartilage lesion scores, mode 5 was the most important positive predictor in the femur (P = 3.9E-06), and mode 2 were predictors, significant negative predictor in the tibia (P = 0.007). In the patella, mode 1 was a significant positive predictor of the bone abnormality score (P = 0.0007).

DATA CONCLUSION

By successfully building an automatic code to create a bone-cartilage interface, we were able to observe dynamic relationships between biochemical changes in the cartilage accompanied with bone remodeling, extended to the whole knee joint instead of simple colocalized observations, shedding light on the interactions that occur between bone and cartilage in OA. Evidence Level: 3 Technical Efficacy: Stage 3 J. Magn. Reson. Imaging 2020;52:1462-1474.

Linear and Nonlinear EEG-Based Functional Networks in Anxiety Disorders

Electrocortical network dynamics are integral to brain function. Linear and nonlinear connectivity applications enrich neurophysiological investigations into anxiety disorders. Discrete EEG-based connectivity networks are unfolding with some homogeneity for anxiety disorder subtypes. Attenuated delta/theta/beta connectivity networks, pertaining to anterior-posterior nodes, characterize panic disorder. Nonlinear measures suggest reduced connectivity of ACC as an executive neuro-regulator in germane "fear circuitry networks" might be more central than considered. Enhanced network complexity and theta network efficiency at rest define generalized anxiety disorder, with similar tonic hyperexcitability apparent in social anxiety disorder further extending to task-related/state functioning. Dysregulated alpha connectivity and integration of mPFC-ACC/mPFC-PCC relays implicated with attentional flexibility and choice execution/congruence neurocircuitry are observed in trait anxiety. Conversely, state anxiety appears to recruit converging delta and beta connectivity networks as panic, suggesting trait and state anxiety are modulated by discrete neurobiological mechanisms. Furthermore, EEG connectivity dynamics distinguish anxiety from depression, despite prevalent clinical comorbidity. Rethinking mechanisms implicated in the etiology, maintenance, and treatment of anxiety from the perspective of EEG network science across micro- and macroscales serves to shed light and move the field forward.