1
|
Achard S, Coeurjolly JF, de Micheaux PL, Lbath H, Richiardi J. Inter-regional correlation estimators for functional magnetic resonance imaging. Neuroimage 2023; 282:120388. [PMID: 37805021 DOI: 10.1016/j.neuroimage.2023.120388] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2023] [Revised: 09/05/2023] [Accepted: 09/22/2023] [Indexed: 10/09/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Sophie Achard
- Univ. Grenoble Alpes, CNRS, Inria, Grenoble-INP, LJK, 38000 Grenoble, France.
| | | | - Pierre Lafaye de Micheaux
- AMIS, Université Paul Valéry Montpellier 3, France; PreMeDICaL - Precision Medicine by Data Integration and Causal Learning, Inria Sophia Antipolis, France; Desbrest Institute of Epidemiology and Public Health, Univ Montpellier, INSERM, Montpellier, France
| | - Hanâ Lbath
- Univ. Grenoble Alpes, CNRS, Inria, Grenoble-INP, LJK, 38000 Grenoble, France
| | - Jonas Richiardi
- Department of Radiology, Lausanne University Hospital and University of Lausanne, Switzerland
| |
Collapse
|
2
|
Zhou Y, Müller HG, Zhu C, Chen Y, Wang JL, O'Muircheartaigh J, Bruchhage M, Deoni S, Bruchhage M, Carnell S, Deoni S, D’Sa V, Huentelman M, Klepac-Ceraj V, LeBourgeois M, Müller HG, O’Muircheartaigh J, Wang JL. Network evolution of regional brain volumes in young children reflects neurocognitive scores and mother's education. Sci Rep 2023; 13:2984. [PMID: 36804963 PMCID: PMC9941570 DOI: 10.1038/s41598-023-29797-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2022] [Accepted: 02/10/2023] [Indexed: 02/22/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Yidong Zhou
- Department of Statistics, University of California, Davis, Davis, CA, 95616, USA.
| | - Hans-Georg Müller
- Department of Statistics, University of California, Davis, Davis, CA, 95616, USA
| | - Changbo Zhu
- Department of Applied and Computational Mathematics and Statistics, University of Notre Dame, Notre Dame, IN, 46556, USA
| | - Yaqing Chen
- Department of Statistics, Rutgers University, New Brunswick, NJ, 08901, USA
| | - Jane-Ling Wang
- Department of Statistics, University of California, Davis, Davis, CA, 95616, USA
| | - Jonathan O'Muircheartaigh
- Centre for the Developing Brain, School of Biomedical Engineering and Imaging Sciences, King's College London, London, UK.,Department of Forensic and Neurodevelopmental Sciences, Institute of Psychiatry, Psychology and Neuroscience, King's College London, London, UK.,MRC Centre for Neurodevelopmental Disorders, King's College London, London, UK
| | - Muriel Bruchhage
- Department of Pediatrics, Warren Alpert Medical School at Brown University, Providence, USA.,Department of Diagnostic Imaging, Rhode Island Hospital, Providence, USA.,Institute of Social Sciences, Stavanger University, Stavanger, 4021, Norway
| | - Sean Deoni
- Maternal, Newborn, and Child Health Discovery and Tools, Bill and Melinda Gates Foundation, Seattle, WA, USA
| | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
3
|
Chen Y, Müller HG. Uniform convergence of local Fréchet regression with applications to locating extrema and time warping for metric space valued trajectories. Ann Stat 2022. [DOI: 10.1214/21-aos2163] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
Affiliation(s)
- Yaqing Chen
- Department of Statistics, University of California, Davis
| | | |
Collapse
|
4
|
Xiang J, Fan C, Wei J, Li Y, Wang B, Niu Y, Yang L, Lv J, Cui X. The Task Pre-Configuration Is Associated With Cognitive Performance Evidence From the Brain Synchrony. Front Comput Neurosci 2022; 16:883660. [PMID: 35603133 PMCID: PMC9120823 DOI: 10.3389/fncom.2022.883660] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2022] [Accepted: 04/05/2022] [Indexed: 11/13/2022] Open
Abstract
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.
Collapse
|
5
|
Exploration of Multiparameter Hematoma 3D Image Analysis for Predicting Outcome After Intracerebral Hemorrhage. Neurocrit Care 2021; 32:539-549. [PMID: 31359310 DOI: 10.1007/s12028-019-00783-8] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
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.
Collapse
|
6
|
Tibrewala R, Pedoia V, Bucknor M, Majumdar S. Principal Component Analysis of Simultaneous PET-MRI Reveals Patterns of Bone-Cartilage Interactions in Osteoarthritis. J Magn Reson Imaging 2020; 52:1462-1474. [PMID: 32207870 PMCID: PMC11090497 DOI: 10.1002/jmri.27146] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2019] [Revised: 03/02/2020] [Accepted: 03/03/2020] [Indexed: 01/04/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Radhika Tibrewala
- Department of Radiology and Biomedical Imaging, University of California, San Francisco, California, USA
| | - Valentina Pedoia
- Department of Radiology and Biomedical Imaging, University of California, San Francisco, California, USA
| | - Matthew Bucknor
- Department of Radiology and Biomedical Imaging, University of California, San Francisco, California, USA
| | - Sharmila Majumdar
- Department of Radiology and Biomedical Imaging, University of California, San Francisco, California, USA
| |
Collapse
|
7
|
Linear and Nonlinear EEG-Based Functional Networks in Anxiety Disorders. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2020; 1191:35-59. [PMID: 32002921 DOI: 10.1007/978-981-32-9705-0_3] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
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.
Collapse
|
8
|
Affiliation(s)
- Alexander Petersen
- Department of Statistics and Applied Probability, University of California, Santa Barbara, California 93106, U.S.A
| | - Hans-Georg Müller
- Department of Statistics, University of California, One Shields Avenue, Davis, California 95616, U.S.A
| |
Collapse
|