Multivariate consistency of resting-state fMRI connectivity maps acquired on a single individual over 2.5 years, 13 sites and 3 vendors

Longitudinal stability of individual brain plasticity patterns in blindness

The primary visual cortex (V1) in blindness is engaged in a wide spectrum of tasks and sensory modalities, including audition, touch, language, and memory. This widespread involvement raises questions regarding the constancy of its role and whether it might exhibit flexibility in its function over time, connecting to diverse network functions specific to task demands. This would suggest that reorganized V1 assumes a role like multiple-demand system regions. Alternatively, varying patterns of plasticity in blind V1 may be attributed to individual factors, with different blind individuals recruiting V1 preferentially for different functions. In support of this, we recently showed that V1 functional connectivity (FC) varies greatly across blind individuals. But do these represent stable individual patterns of plasticity, or are they driven more by instantaneous changes, like a multiple-demand system now inhabiting V1? Here, we tested whether individual FC patterns from the V1 of blind individuals are stable over time. We show that over two years, FC from the V1 is unique and highly stable in a small sample of repeatedly sampled congenitally blind individuals. Further, using multivoxel pattern analysis, we demonstrate that the unique reorganization patterns of these individuals allow decoding of participant identity. Together with recent evidence for substantial individual differences in V1 connectivity, this indicates that there may be a consistent role for V1 in blindness, which may differ for each individual. Further, it suggests that the variability in visual reorganization in blindness across individuals could be used to seek stable neuromarkers for sight rehabilitation and assistive approaches.

DeepComBat: A statistically motivated, hyperparameter-robust, deep learning approach to harmonization of neuroimaging data

Neuroimaging data acquired using multiple scanners or protocols are increasingly available. However, such data exhibit technical artifacts across batches which introduce confounding and decrease reproducibility. This is especially true when multi-batch data are analyzed using complex downstream models which are more likely to pick up on and implicitly incorporate batch-related information. Previously proposed image harmonization methods have sought to remove these batch effects; however, batch effects remain detectable in the data after applying these methods. We present DeepComBat, a deep learning harmonization method based on a conditional variational autoencoder and the ComBat method. DeepComBat combines the strengths of statistical and deep learning methods in order to account for the multivariate relationships between features while simultaneously relaxing strong assumptions made by previous deep learning harmonization methods. As a result, DeepComBat can perform multivariate harmonization while preserving data structure and avoiding the introduction of synthetic artifacts. We apply this method to cortical thickness measurements from a cognitive-aging cohort and show DeepComBat qualitatively and quantitatively outperforms existing methods in removing batch effects while preserving biological heterogeneity. Additionally, DeepComBat provides a new perspective for statistically motivated deep learning harmonization methods.

Myths and Methodologies: Understanding the health impact of head down bedrest for the benefit of older adults and astronauts. Study protocol of the Canadian Bedrest Study

Weightlessness during spaceflight can harm various bodily systems, including bone density, muscle mass, strength and cognitive functions. Exercise appears to somewhat counteract these effects. A terrestrial model for this is head-down bedrest (HDBR), simulating gravity loss. This mirrors challenges faced by older adults in extended bedrest and space environments. The first Canadian study, backed by the Canadian Space Agency, Canadian Institutes of Health Research, and Canadian Frailty Network, aims to explore these issues. The study seeks to: (1) scrutinize the impact of 14-day HDBR on physiological, psychological and neurocognitive systems, and (2) assess the benefits of exercise during HDBR. Eight teams developed distinct protocols, harmonized in three videoconferences, at the McGill University Health Center. Over 26 days, 23 participants aged 55-65 underwent baseline measurements, 14 days of -6° HDBR, and 7 days of recovery. Half did prescribed exercise thrice daily combining resistance and endurance exercise for a total duration of 1 h. Assessments included demographics, cardiorespiratory fitness, bone health, body composition, quality of life, mental health, cognition, muscle health and biomarkers. This study has yielded some published outcomes, with more forthcoming. Findings will enrich our comprehension of HDBR effects, guiding future strategies for astronaut well-being and aiding bedrest-bound older adults. By outlining evidence-based interventions, this research supports both space travellers and those enduring prolonged bedrest.

Replication of a neuroimaging biomarker for striatal dysfunction in psychosis

To bring biomarkers closer to clinical application, they should be generalizable, reliable, and maintain performance within the constraints of routine clinical conditions. The functional striatal abnormalities (FSA), is among the most advanced neuroimaging biomarkers in schizophrenia, trained to discriminate diagnosis, with post-hoc analyses indicating prognostic properties. Here, we attempt to replicate its diagnostic capabilities measured by the area under the curve (AUC) in receiver operator characteristic curves discriminating individuals with psychosis (n = 101) from healthy controls (n = 51) in the Human Connectome Project for Early Psychosis. We also measured the test-retest (run 1 vs 2) and phase encoding direction (i.e., AP vs PA) reliability with intraclass correlation coefficients (ICC). Additionally, we measured effects of scan length on classification accuracy (i.e., AUCs) and reliability (i.e., ICCs). Finally, we tested the prognostic capability of the FSA by the correlation between baseline scores and symptom improvement over 12 weeks of antipsychotic treatment in a separate cohort (n = 97). Similar analyses were conducted for the Yeo networks intrinsic connectivity as a reference. The FSA had good/excellent diagnostic discrimination (AUC = 75.4%, 95% CI = 67.0-83.3%; in non-affective psychosis AUC = 80.5%, 95% CI = 72.1-88.0%, and in affective psychosis AUC = 58.7%, 95% CI = 44.2-72.0%). Test-retest reliability ranged between ICC = 0.48 (95% CI = 0.35-0.59) and ICC = 0.22 (95% CI = 0.06-0.36), which was comparable to that of networks intrinsic connectivity. Phase encoding direction reliability for the FSA was ICC = 0.51 (95% CI = 0.42-0.59), generally lower than for networks intrinsic connectivity. By increasing scan length from 2 to 10 min, diagnostic classification of the FSA increased from AUC = 71.7% (95% CI = 63.1-80.3%) to 75.4% (95% CI = 67.0-83.3%) and phase encoding direction reliability from ICC = 0.29 (95% CI = 0.14-0.43) to ICC = 0.51 (95% CI = 0.42-0.59). FSA scores did not correlate with symptom improvement. These results reassure that the FSA is a generalizable diagnostic - but not prognostic - biomarker. Given the replicable results of the FSA as a diagnostic biomarker trained on case-control datasets, next the development of prognostic biomarkers should be on treatment-response data.

Longitudinal stability of individual brain plasticity patterns in blindness

The primary visual cortex (V1) in individuals born blind is engaged in a wide spectrum of tasks and sensory modalities, including audition, touch, language, and memory. This widespread involvement raises questions regarding the constancy of its role and whether it might exhibit flexibility in its function over time, connecting to diverse network functions in response to task-specific demands. This would suggest that reorganized V1 takes on a role similar to cognitive multiple-demand system regions. Alternatively, it is possible that the varying patterns of plasticity observed in the blind V1 can be attributed to individual factors, whereby different blind individuals recruit V1 for different functions, highlighting the immense idiosyncrasy of plasticity. In support of this second account, we have recently shown that V1 functional connectivity varies greatly across blind individuals. But do these represent stable individual patterns of plasticity or merely instantaneous changes, for a multiple-demand system now inhabiting V1? Here we tested if individual connectivity patterns from the visual cortex of blind individuals are stable over time. We show that over two years, fMRI functional connectivity from the primary visual cortex is unique and highly stable in a small sample of repeatedly sampled congenitally blind individuals. Further, using multivoxel pattern analysis, we demonstrate that the unique reorganization patterns of these individuals allow decoding of participant identity. Together with recent evidence for substantial individual differences in visual cortex connectivity, this indicates there may be a consistent role for the visual cortex in blindness, which may differ for each individual. Further, it suggests that the variability in visual reorganization in blindness across individuals could be used to seek stable neuromarkers for sight rehabilitation and assistive approaches.

Structural Fingerprinting of the Frontal Aslant Tract: Predicting Cognitive Control Capacity and Obsessive-Compulsive Symptoms

White matter of the human brain is influenced by common genetic variations and shaped by neural activity-dependent experiences. Variations in microstructure of cerebral white matter across individuals and even across fiber tracts might underlie differences in cognitive capacity and vulnerabilities to mental disorders. The frontoparietal and cingulo-opercular networks of the brain constitute the central system supporting cognitive functions, and functional connectivity of these networks has been used to distinguish individuals known as "functional fingerprinting." The frontal aslant tract (FAT) that passes through the two networks has been implicated in executive functions. However, whether FAT can be used as a "structural fingerprint" to distinguish individuals and predict an individual's cognitive function and dysfunction is unknown. Here we investigated the fingerprinting property of FAT microstructural profiles using three independent diffusion MRI datasets with repeated scans on human participants including both females and males. We found that diffusion and geometric profiles of FAT can be used to distinguish individuals with a high accuracy. Next, we demonstrated that fractional anisotropy in different FAT segments predicted distinct cognitive functions, including working memory, inhibitory control, and relational reasoning. Finally, we assessed the contribution of altered FAT microstructural profiles to cognitive dysfunction in unmedicated patients with obsessive-compulsive disorders. We found that the altered microstructure in FAT was associated with the severity of obsessive-compulsive symptoms. Collectively, our findings suggest that the microstructural profiles of FAT can identify individuals with a high accuracy and may serve as an imaging marker for predicting an individual's cognitive capacity and disease severity.SIGNIFICANCE STATEMENT The frontoparietal network and cingulo-opercular network of the brain constitute a dual-network architecture for human cognitive functions, and functional connectivity of these two networks can be used as a "functional fingerprint" to distinguish individuals. However, the structural underpinnings of these networks subserving individual heterogeneities in their functional connectivity and cognitive ability remain unknown. We show here that the frontal aslant tract (FAT) that passes through the two networks distinguishes individuals with a high accuracy. Further, we demonstrate that the diffusion profiles of FAT predict distinct cognitive functions in healthy subjects and are associated with the clinical symptoms in patients with obsessive-compulsive disorders. Our findings suggest that the FAT may serve as a unique structural fingerprint underlying individual cognitive capability.

Basal forebrain activity predicts functional degeneration in the entorhinal cortex in Alzheimer's disease

Recent models of Alzheimer's disease suggest the nucleus basalis of Meynert (NbM) as an early origin of structural degeneration followed by the entorhinal cortex (EC). However, the functional properties of NbM and EC regarding amyloid-β and hyperphosphorylated tau remain unclear. We analysed resting-state functional fMRI data with CSF assays from the Alzheimer's Disease Neuroimaging Initiative (n = 71) at baseline and 2 years later. At baseline, local activity, as quantified by fractional amplitude of low-frequency fluctuations, differentiated between normal and abnormal CSF groups in the NbM but not EC. Further, NbM activity linearly decreased as a function of CSF ratio, resembling the disease status. Finally, NbM activity predicted the annual percentage signal change in EC, but not the reverse, independent from CSF ratio. Our findings give novel insights into the pathogenesis of Alzheimer's disease by showing that local activity in NbM is affected by proteinopathology and predicts functional degeneration within the EC.

Confounding Effects on the Performance of Machine Learning Analysis of Static Functional Connectivity Computed from rs-fMRI Multi-site Data

Resting-state functional magnetic resonance imaging (rs-fMRI) is a non-invasive imaging technique widely used in neuroscience to understand the functional connectivity of the human brain. While rs-fMRI multi-site data can help to understand the inner working of the brain, the data acquisition and processing of this data has many challenges. One of the challenges is the variability of the data associated with different acquisitions sites, and different MRI machines vendors. Other factors such as population heterogeneity among different sites, with variables such as age and gender of the subjects, must also be considered. Given that most of the machine-learning models are developed using these rs-fMRI multi-site data sets, the intrinsic confounding effects can adversely affect the generalizability and reliability of these computational methods, as well as the imposition of upper limits on the classification scores. This work aims to identify the phenotypic and imaging variables producing the confounding effects, as well as to control these effects. Our goal is to maximize the classification scores obtained from the machine learning analysis of the Autism Brain Imaging Data Exchange (ABIDE) rs-fMRI multi-site data. To achieve this goal, we propose novel methods of stratification to produce homogeneous sub-samples of the 17 ABIDE sites, as well as the generation of new features from the static functional connectivity values, using multiple linear regression models, ComBat harmonization models, and normalization methods. The main results obtained with our statistical models and methods are an accuracy of 76.4%, sensitivity of 82.9%, and specificity of 77.0%, which are 8.8%, 20.5%, and 7.5% above the baseline classification scores obtained from the machine learning analysis of the static functional connectivity computed from the ABIDE rs-fMRI multi-site data.

Replication of a neuroimaging biomarker for striatal dysfunction in psychosis

To bring biomarkers closer to clinical application, they should be generalizable, reliable, and maintain performance within the constraints of routine clinical conditions. The functional striatal abnormalities (FSA), is among the most advanced neuroimaging biomarkers in schizophrenia, trained to discriminate diagnosis, with post-hoc analyses indicating prognostic properties. Here, we attempt to replicate its diagnostic capabilities measured by the area under the curve (AUC) in receiver operator characteristic curves discriminating individuals with psychosis (n=101) from healthy controls (n=51) in the Human Connectome Project for Early Psychosis. We also measured the test-retest (run 1 vs 2) and phase encoding direction (i.e., AP vs PA) reliability with intraclass correlation coefficients (ICC). Additionally, we measured effects of scan length on classification accuracy (i.e., AUCs) and reliability (i.e., ICCs). Finally, we tested the prognostic capability of the FSA by the correlation between baseline scores and symptom improvement over 12 weeks of antipsychotic treatment in a separate cohort (n=97). Similar analyses were conducted for the Yeo networks intrinsic connectivity as a reference. The FSA had good/excellent diagnostic discrimination (AUC=75.4%, 95%CI=67.0%-83.3%; in non-affective psychosis AUC=80.5%, 95%CI=72.1-88.0%, and in affective psychosis AUC=58.7%, 95%CI=44.2-72.0%). Test-retest reliability ranged between ICC=0.48 (95%CI=0.35-0.59) and ICC=0.22 (95%CI=0.06-0.36), which was comparable to that of networks intrinsic connectivity. Phase encoding direction reliability for the FSA was ICC=0.51 (95%CI=0.42-0.59), generally lower than for networks intrinsic connectivity. By increasing scan length from 2 to 10 minutes, diagnostic classification of the FSA increased from AUC=71.7% (95%CI=63.1%-80.3%) to 75.4% (95%CI=67.0%-83.3%) and phase encoding direction reliability from ICC=0.29 (95%CI=0.14-0.43) to ICC=0.51 (95%CI=0.42-0.59). FSA scores did not correlate with symptom improvement. These results reassure that the FSA is a generalizable diagnostic - but not prognostic - biomarker. Given the replicable results of the FSA as a diagnostic biomarker trained on case-control datasets, next the development of prognostic biomarkers should be on treatment-response data.

Replication of a neuroimaging biomarker for striatal dysfunction in psychosis

To bring biomarkers closer to clinical application, they should be generalizable, reliable, and maintain performance within the constraints of routine clinical conditions. The functional striatal abnormalities (FSA), is among the most advanced neuroimaging biomarkers in schizophrenia, trained to discriminate diagnosis, with post-hoc analyses indicating prognostic properties. Here, we attempt to replicate its diagnostic capabilities measured by the area under the curve (AUC) in receiver operator characteristic curves discriminating individuals with psychosis (n=101) from healthy controls (n=51) in the Human Connectome Project for Early Psychosis. We also measured the test-retest (run 1 vs 2) and phase encoding direction (i.e., AP vs PA) reliability with intraclass correlation coefficients (ICC). Additionally, we measured effects of scan length on classification accuracy (i.e., AUCs) and reliability (i.e., ICCs). Finally, we tested the prognostic capability of the FSA by the correlation between baseline scores and symptom improvement over 12 weeks of antipsychotic treatment in a separate cohort (n=97). Similar analyses were conducted for the Yeo networks intrinsic connectivity as a reference. The FSA had good/excellent diagnostic discrimination (AUC=75.4%, 95%CI=67.0%-83.3%; in non-affective psychosis AUC=80.5%, 95%CI=72.1-88.0%, and in affective psychosis AUC=58.7%, 95%CI=44.2-72.0%). Test-retest reliability ranged between ICC=0.48 (95%CI=0.35-0.59) and ICC=0.22 (95%CI=0.06-0.36), which was comparable to that of networks intrinsic connectivity. Phase encoding direction reliability for the FSA was ICC=0.51 (95%CI=0.42-0.59), generally lower than for networks intrinsic connectivity. By increasing scan length from 2 to 10 minutes, diagnostic classification of the FSA increased from AUC=71.7% (95%CI=63.1%-80.3%) to 75.4% (95%CI=67.0%-83.3%) and phase encoding direction reliability from ICC=0.29 (95%CI=0.14-0.43) to ICC=0.51 (95%CI=0.42-0.59). FSA scores did not correlate with symptom improvement. These results reassure that the FSA is a generalizable diagnostic - but not prognostic - biomarker. Given the replicable results of the FSA as a diagnostic biomarker trained on case-control datasets, next the development of prognostic biomarkers should be on treatment-response data.

Hormonal contraceptive exposure relates to changes in resting state functional connectivity of anterior cingulate cortex and amygdala

Introduction

Hormonal contraceptives (HCs), nowadays one of the most used contraceptive methods, downregulate endogenous ovarian hormones, which have multiple plastic effects in the adult brain. HCs usually contain a synthetic estrogen, ethinyl-estradiol, and a synthetic progestin, which can be classified as androgenic or anti-androgenic, depending on their interaction with androgen receptors. Both the anterior cingulate cortex (ACC) and the amygdala express steroid receptors and have shown differential functionality depending on the hormonal status of the participant and the use of HC. In this work, we investigated for the first time the relationship between ACC and amygdala resting state functional connectivity (rs-FC) and HC use duration, while controlling for progestin androgenicity.

Methods

A total of 231 healthy young women participated in five different magnetic resonance imaging studies and were included in the final analysis. The relation between HC use duration and (i) gray matter volume, (ii) fractional amplitude of low-frequency fluctuations, and (iii) seed-based connectivity during resting state in the amygdalae and ACC was investigated in this large sample of women.

Results

In general, rs-FC of the amygdalae with frontal areas, and between the ACC and temporoparietal areas, decreased the longer the HC exposure and independently of the progestin's androgenicity. The type of HC's progestin did show a differential effect in the gray matter volume of left ACC and the connectivity between bilateral ACC and the right inferior frontal gyrus.

Image harmonization: A review of statistical and deep learning methods for removing batch effects and evaluation metrics for effective harmonization

Magnetic resonance imaging and computed tomography from multiple batches (e.g. sites, scanners, datasets, etc.) are increasingly used alongside complex downstream analyses to obtain new insights into the human brain. However, significant confounding due to batch-related technical variation, called batch effects, is present in this data; direct application of downstream analyses to the data may lead to biased results. Image harmonization methods seek to remove these batch effects and enable increased generalizability and reproducibility of downstream results. In this review, we describe and categorize current approaches in statistical and deep learning harmonization methods. We also describe current evaluation metrics used to assess harmonization methods and provide a standardized framework to evaluate newly-proposed methods for effective harmonization and preservation of biological information. Finally, we provide recommendations to end-users to advocate for more effective use of current methods and to methodologists to direct future efforts and accelerate development of the field.

DeepComBat: A Statistically Motivated, Hyperparameter-Robust, Deep Learning Approach to Harmonization of Neuroimaging Data

Neuroimaging data from multiple batches (i.e. acquisition sites, scanner manufacturer, datasets, etc.) are increasingly necessary to gain new insights into the human brain. However, multi-batch data, as well as extracted radiomic features, exhibit pronounced technical artifacts across batches. These batch effects introduce confounding into the data and can obscure biological effects of interest, decreasing the generalizability and reproducibility of findings. This is especially true when multi-batch data is used alongside complex downstream analysis models, such as machine learning methods. Image harmonization methods seeking to remove these batch effects are important for mitigating these issues; however, significant multivariate batch effects remain in the data following harmonization by current state-of-the-art statistical and deep learning methods. We present DeepCombat, a deep learning harmonization method based on a conditional variational autoencoder architecture and the ComBat harmonization model. DeepCombat learns and removes subject-level batch effects by accounting for the multivariate relationships between features. Additionally, DeepComBat relaxes a number of strong assumptions commonly made by previous deep learning harmonization methods and is empirically robust across a wide range of hyperparameter choices. We apply this method to neuroimaging data from a large cognitive-aging cohort and find that DeepCombat outperforms existing methods, as assessed by a battery of machine learning methods, in removing scanner effects from cortical thickness measurements while preserving biological heterogeneity. Additionally, DeepComBat provides a new perspective for statistically-motivated deep learning harmonization methods.

Basal forebrain activity predicts functional degeneration in the entorhinal cortex and decreases with Alzheimer's Disease progression

BACKGROUND AND OBJECTIVES

Recent models of Alzheimer's Disease (AD) suggest the nucleus basalis of Meynert (NbM) as the origin of structural degeneration followed by the entorhinal cortex (EC). However, the functional properties of NbM and EC regarding amyloid-β and hyperphosphorylated tau remain unclear.

METHODS

We analyzed resting-state (rs)fMRI data with CSF assays from the Alzheimer's Disease Neuroimaging Initiative (ADNI, n=71) at baseline and two years later.

RESULTS

At baseline, local activity, as quantified by fractional amplitude of low-frequency fluctuations (fALFF), differentiated between normal and abnormal CSF groups in the NbM but not EC. Further, NbM activity linearly decreased as a function of CSF ratio, resembling the disease status. Finally, NbM activity predicted the annual percentage signal change in EC, but not the reverse, independent from CSF ratio.

DISCUSSION

Our findings give novel insights into the pathogenesis of AD by showing that local activity in NbM is affected by proteinopathology and predicts functional degeneration within the EC.

Functional magnetic resonance imaging of headache: Issues, best-practices, and new directions, a narrative review

OBJECTIVE

To ensure readers are informed consumers of functional magnetic resonance imaging (fMRI) research in headache, to outline ongoing challenges in this area of research, and to describe potential considerations when asked to collaborate on fMRI research in headache, as well as to suggest future directions for improvement in the field.

BACKGROUND

Functional MRI has played a key role in understanding headache pathophysiology, and mapping networks involved with headache-related brain activity have the potential to identify intervention targets. Some investigators have also begun to explore its use for diagnosis.

METHODS/RESULTS

The manuscript is a narrative review of the current best practices in fMRI in headache research, including guidelines on transparency and reproducibility. It also contains an outline of the fundamentals of MRI theory, task-related study design, resting-state functional connectivity, relevant statistics and power analysis, image preprocessing, and other considerations essential to the field.

CONCLUSION

Best practices to increase reproducibility include methods transparency, eliminating error, using a priori hypotheses and power calculations, using standardized instruments and diagnostic criteria, and developing large-scale, publicly available datasets.

Resting-state functional MRI in multicenter studies on multiple sclerosis: a report on raw data quality and functional connectivity features from the Italian Neuroimaging Network Initiative

The Italian Neuroimaging Network Initiative (INNI) is an expanding repository of brain MRI data from multiple sclerosis (MS) patients recruited at four Italian MRI research sites. We describe the raw data quality of resting-state functional MRI (RS-fMRI) time-series in INNI and the inter-site variability in functional connectivity (FC) features after unified automated data preprocessing. MRI datasets from 489 MS patients and 246 healthy control (HC) subjects were retrieved from the INNI database. Raw data quality metrics included temporal signal-to-noise ratio (tSNR), spatial smoothness (FWHM), framewise displacement (FD), and differential variation in signals (DVARS). Automated preprocessing integrated white-matter lesion segmentation (SAMSEG) into a standard fMRI pipeline (fMRIPrep). FC features were calculated on pre-processed data and harmonized between sites (Combat) prior to assessing general MS-related alterations. Across centers (both groups), median tSNR and FWHM ranged from 47 to 84 and from 2.0 to 2.5, and median FD and DVARS ranged from 0.08 to 0.24 and from 1.06 to 1.22. After preprocessing, only global FC-related features were significantly correlated with FD or DVARS. Across large-scale networks, age/sex/FD-adjusted and harmonized FC features exhibited both inter-site and site-specific inter-group effects. Significant general reductions were obtained for somatomotor and limbic networks in MS patients (vs. HC). The implemented procedures provide technical information on raw data quality and outcome of fully automated preprocessing that might serve as reference in future RS-fMRI studies within INNI. The unified pipeline introduced little bias across sites and appears suitable for multisite FC analyses on harmonized network estimates.

Comparing the test-retest reliability of resting-state functional magnetic resonance imaging metrics across single band and multiband acquisitions in the context of healthy aging

The identification of meaningful functional magnetic resonance imaging (fMRI) biomarkers requires measures that reliably capture brain performance across different subjects and over multiple scanning sessions. Recent developments in fMRI acquisition, such as the introduction of multiband (MB) protocols and in-plane acceleration, allow for increased scanning speed and improved temporal resolution. However, they may also lead to reduced temporal signal to noise ratio and increased signal leakage between simultaneously excited slices. These methods have been adopted in several scanning modalities including diffusion weighted imaging and fMRI. To our knowledge, no study has formally compared the reliability of the same resting-state fMRI (rs-fMRI) metrics (amplitude of low-frequency fluctuations; seed-to-voxel and region of interest [ROI]-to-ROI connectivity) across conventional single-band fMRI and different MB acquisitions, with and without in-plane acceleration, across three sessions. In this study, 24 healthy older adults were scanned over three visits, on weeks 0, 1, and 4, and, on each occasion, underwent a conventional single band rs-fMRI scan and three different rs-fMRI scans with MB factors 4 and 6, with and without in-plane acceleration. Across all three rs-fMRI metrics, the reliability scores were highest with MB factor 4 with no in-plane acceleration for cortical areas and with conventional single band for subcortical areas. Recommendations for future research studies are discussed.

Person-specific and precision neuroimaging: Current methods and future directions

Most neuroimaging studies of brain function analyze data in normalized space to identify regions of common activation across participants. These studies treat interindividual differences in brain organization as noise, but this approach can obscure important information about the brain's functional architecture. Recently, a number of studies have adopted a person-specific approach that aims to characterize these individual differences and explore their reliability and implications for behavior. A subset of these studies has taken a precision imaging approach that collects multiple hours of data from each participant to map brain function on a finer scale. In this review, we provide a broad overview of how person-specific and precision imaging techniques have used resting-state measures to examine individual differences in the brain's organization and their impact on behavior, followed by how task-based activity continues to add detail to these discoveries. We argue that person-specific and precision approaches demonstrate substantial promise in uncovering new details of the brain's functional organization and its relationship to behavior in many areas of cognitive neuroscience. We also discuss some current limitations in this new field and some new directions it may take.

MAGNIMS recommendations for harmonization of MRI data in MS multicenter studies

•

Sharing data from cooperative studies is essential to develop new biomarkers in MS.

•

Differences in MRI acquisition, analysis, storage represent a substantial constraint.

•

We review the state of the art and developments in the harmonization of MRI.

•

We provide recommendations to harmonize large MRI datasets in the MS field.

There is an increasing need of sharing harmonized data from large, cooperative studies as this is essential to develop new diagnostic and prognostic biomarkers. In the field of multiple sclerosis (MS), the issue has become of paramount importance due to the need to translate into the clinical setting some of the most recent MRI achievements. However, differences in MRI acquisition parameters, image analysis and data storage across sites, with their potential bias, represent a substantial constraint. This review focuses on the state of the art, recent technical advances, and desirable future developments of the harmonization of acquisition, analysis and storage of large-scale multicentre MRI data of MS cohorts. Huge efforts are currently being made to achieve all the requirements needed to provide harmonized MRI datasets in the MS field, as proper management of large imaging datasets is one of our greatest opportunities and challenges in the coming years. Recommendations based on these achievements will be provided here. Despite the advances that have been made, the complexity of these tasks requires further research by specialized academical centres, with dedicated technical and human resources. Such collective efforts involving different professional figures are of crucial importance to offer to MS patients a personalised management while minimizing consumption of resources.

Quantitative MRI Harmonization to Maximize Clinical Impact: The RIN-Neuroimaging Network

Neuroimaging studies often lack reproducibility, one of the cardinal features of the scientific method. Multisite collaboration initiatives increase sample size and limit methodological flexibility, therefore providing the foundation for increased statistical power and generalizable results. However, multisite collaborative initiatives are inherently limited by hardware, software, and pulse and sequence design heterogeneities of both clinical and preclinical MRI scanners and the lack of benchmark for acquisition protocols, data analysis, and data sharing. We present the overarching vision that yielded to the constitution of RIN-Neuroimaging Network, a national consortium dedicated to identifying disease and subject-specific in-vivo neuroimaging biomarkers of diverse neurological and neuropsychiatric conditions. This ambitious goal needs efforts toward increasing the diagnostic and prognostic power of advanced MRI data. To this aim, 23 Italian Scientific Institutes of Hospitalization and Care (IRCCS), with technological and clinical specialization in the neurological and neuroimaging field, have gathered together. Each IRCCS is equipped with high- or ultra-high field MRI scanners (i.e., ≥3T) for clinical or preclinical research or has established expertise in MRI data analysis and infrastructure. The actions of this Network were defined across several work packages (WP). A clinical work package (WP1) defined the guidelines for a minimum standard clinical qualitative MRI assessment for the main neurological diseases. Two neuroimaging technical work packages (WP2 and WP3, for clinical and preclinical scanners) established Standard Operative Procedures for quality controls on phantoms as well as advanced harmonized quantitative MRI protocols for studying the brain of healthy human participants and wild type mice. Under FAIR principles, a web-based e-infrastructure to store and share data across sites was also implemented (WP4). Finally, the RIN translated all these efforts into a large-scale multimodal data collection in patients and animal models with dementia (i.e., case study). The RIN-Neuroimaging Network can maximize the impact of public investments in research and clinical practice acquiring data across institutes and pathologies with high-quality and highly-consistent acquisition protocols, optimizing the analysis pipeline and data sharing procedures.

ENIGMA brain injury: Framework, challenges, and opportunities

Traumatic brain injury (TBI) is a major cause of disability worldwide, but the heterogeneous nature of TBI with respect to injury severity and health comorbidities make patient outcome difficult to predict. Injury severity accounts for only some of this variance, and a wide range of preinjury, injury-related, and postinjury factors may influence outcome, such as sex, socioeconomic status, injury mechanism, and social support. Neuroimaging research in this area has generally been limited by insufficient sample sizes. Additionally, development of reliable biomarkers of mild TBI or repeated subconcussive impacts has been slow, likely due, in part, to subtle effects of injury and the aforementioned variability. The ENIGMA Consortium has established a framework for global collaboration that has resulted in the largest-ever neuroimaging studies of multiple psychiatric and neurological disorders. Here we describe the organization, recent progress, and future goals of the Brain Injury working group.

Improving Between-Group Effect Size for Multi-Site Functional Connectivity Data via Site-Wise De-Meaning

Background: Multi-site functional MRI (fMRI) databases are becoming increasingly prevalent in the study of neurodevelopmental and psychiatric disorders. However, multi-site databases are known to introduce site effects that may confound neurobiological and measures such as functional connectivity (FC). Although studies have been conducted to mitigate site effects, these methods often result in reduced effect size in FC comparisons between controls and patients.

Methods: We present a site-wise de-meaning (SWD) strategy in multi-site FC analysis and compare its performance with two common site-effect mitigation methods, i.e., generalized linear model (GLM) and Combining Batches (ComBat) Harmonization. For SWD, after FC was calculated and Fisher z-transformed, the site-wise FC mean was removed from each subject before group-level statistical analysis. The above methods were tested on two multi-site psychiatric consortiums [Autism Brain Imaging Data Exchange (ABIDE) and Bipolar and Schizophrenia Network on Intermediate Phenotypes (B-SNIP)]. Preservation of consistent FC alterations in patients were evaluated for each method through the effect sizes (Hedge’s g) of patients vs. controls.

Results: For the B-SNIP dataset, SWD improved the effect size between schizophrenic and control subjects by 4.5–7.9%, while GLM and ComBat decreased the effect size by 22.5–42.6%. For the ABIDE dataset, SWD improved the effect size between autistic and control subjects by 2.9–5.3%, while GLM and ComBat decreased the effect size by up to 11.4%.

Conclusion: Compared to the original data and commonly used methods, the SWD method demonstrated superior performance in preserving the effect size in FC features associated with disorders.

Analysis of Functional Brain Network in MDD Based on Improved Empirical Mode Decomposition With Resting State EEG Data

At present, most brain functional studies are based on traditional frequency bands to explore the abnormal functional connections and topological organization of patients with depression. However, they ignore the characteristic relationship of electroencephalogram (EEG) signals in the time domain. Therefore, this paper proposes a network decomposition model based on Improved Empirical Mode Decomposition (EMD), it is suitable for time-frequency analysis of brain functional network. On the one hand, it solves the problem of mode mixing on original EMD method, especially on high-density EEG data. On the other hand, by building brain function networks on different intrinsic mode function (IMF), we can perform time-frequency analysis of brain function connections. It provides a new insight for brain function connectivity analysis of major depressive disorder (MDD). Experimental results found that the IMFs waveform decomposed by Improved EMD was more stable and the difference between IMFs was obvious, it indicated that the mode mixing can be effectively solved. Besides, the analysis of the brain network, we found that the changes in MDD functional connectivity on different IMFs, it may be related to the pathological changes for MDD. More statistical results on three network metrics proved that there were significant differences between MDD and normal controls (NC) group. In addition, the aberrant brain network structure of MDDs was also confirmed in the hubs characteristic. These findings may provide potential biomarkers for the clinical diagnosis of MDD patients.

A guide to the measurement and interpretation of fMRI test-retest reliability

The test-retest reliability of functional neuroimaging data has recently been a topic of much discussion. Despite early conflicting reports, converging reports now suggest that test-retest reliability is poor for standard univariate measures-namely, voxel- and region-level task-based activation and edge-level functional connectivity. To better understand the implications of these recent studies requires understanding the nuances of test-retest reliability as commonly measured by the intraclass correlation coefficient (ICC). Here we provide a guide to the measurement and interpretation of test-retest reliability in functional neuroimaging and review major findings in the literature. We highlight the importance of making choices that improve reliability so long as they do not diminish validity, pointing to the potential of multivariate approaches that improve both. Finally, we discuss the implications of recent reports of low test-retest reliability in the context of ongoing work in the field.

Effects of prenatal opioid exposure on functional networks in infancy

Prenatal opioid exposure has been linked to altered neurodevelopment and visual problems such as strabismus and nystagmus. The neural substrate underlying these alterations is unclear. Resting-state functional connectivity MRI (rsfMRI) is an advanced and well-established technique to evaluate brain networks. Few studies have examined the effects of prenatal opioid exposure on resting-state network connectivity in infancy. In this pilot study, we characterized network connectivity in opioid-exposed infants (n = 19) and controls (n = 20) between 4–8 weeks of age using both a whole-brain connectomic approach and a seed-based approach. Prenatal opioid exposure was associated with differences in distribution of betweenness centrality and connection length, with positive connections unique to each group significantly longer than common connections. The unique connections in the opioid-exposed group were more often inter-network connections while unique connections in controls and connections common to both groups were more often intra-network. The opioid-exposed group had smaller network volumes particularly in the primary visual network, but similar network strength as controls. Network topologies as determined by dice similarity index were different between groups, particularly in visual and executive control networks. These results may provide insight into the neural basis for the developmental and visual problems associated with prenatal opioid exposure.

Resample aggregating improves the generalizability of connectome predictive modeling

It is a longstanding goal of neuroimaging to produce reliable, generalizable models of brain behavior relationships. More recently, data driven predictive models have become popular. However, overfitting is a common problem with statistical models, which impedes model generalization. Cross validation (CV) is often used to estimate expected model performance within sample. Yet, the best way to generate brain behavior models, and apply them out-of-sample, on an unseen dataset, is unclear. As a solution, this study proposes an ensemble learning method, in this case resample aggregating, encompassing both model parameter estimation and feature selection. Here we investigate the use of resampled aggregated models when used to estimate fluid intelligence (fIQ) from fMRI based functional connectivity (FC) data. We take advantage of two large openly available datasets, the Human Connectome Project (HCP), and the Philadelphia Neurodevelopmental Cohort (PNC). We generate aggregated and non-aggregated models of fIQ in the HCP, using the Connectome Prediction Modelling (CPM) framework. Over various test-train splits, these models are evaluated in sample, on left-out HCP data, and out-of-sample, on PNC data. We find that a resample aggregated model performs best both within- and out-of-sample. We also find that feature selection can vary substantially within-sample. More robust feature selection methods, as detailed here, are needed to improve cross sample performance of CPM based brain behavior models.

The relationship between EEG and fMRI connectomes is reproducible across simultaneous EEG-fMRI studies from 1.5T to 7T

Both electroencephalography (EEG) and functional Magnetic Resonance Imaging (fMRI) are non-invasive methods that show complementary aspects of human brain activity. Despite measuring different proxies of brain activity, both the measured blood-oxygenation (fMRI) and neurophysiological recordings (EEG) are indirectly coupled. The electrophysiological and BOLD signal can map the underlying functional connectivity structure at the whole brain scale at different timescales. Previous work demonstrated a moderate but significant correlation between resting-state functional connectivity of both modalities, however there is a wide range of technical setups to measure simultaneous EEG-fMRI and the reliability of those measures between different setups remains unknown. This is true notably with respect to different magnetic field strengths (low and high field) and different spatial sampling of EEG (medium to high-density electrode coverage). Here, we investigated the reproducibility of the bimodal EEG-fMRI functional connectome in the most comprehensive resting-state simultaneous EEG-fMRI dataset compiled to date including a total of 72 subjects from four different imaging centers. Data was acquired from 1.5T, 3T and 7T scanners with simultaneously recorded EEG using 64 or 256 electrodes. We demonstrate that the whole-brain monomodal connectivity reproducibly correlates across different datasets and that a moderate crossmodal correlation between EEG and fMRI connectivity of r ≈ 0.3 can be reproducibly extracted in low- and high-field scanners. The crossmodal correlation was strongest in the EEG-β frequency band but exists across all frequency bands. Both homotopic and within intrinsic connectivity network (ICN) connections contributed the most to the crossmodal relationship. This study confirms, using a considerably diverse range of recording setups, that simultaneous EEG-fMRI offers a consistent estimate of multimodal functional connectomes in healthy subjects that are dominantly linked through a functional core of ICNs across spanning across the different timescales measured by EEG and fMRI. This opens new avenues for estimating the dynamics of brain function and provides a better understanding of interactions between EEG and fMRI measures. This observed level of reproducibility also defines a baseline for the study of alterations of this coupling in pathological conditions and their role as potential clinical markers.

Vascular Contributions to Neurodegeneration: Protocol of the COMPASS-ND Study

OBJECTIVE

To describe the neuroimaging and other methods for assessing vascular contributions to neurodegeneration in the Comprehensive Assessment of Neurodegeneration and Dementia (COMPASS-ND) study, a Canadian multi-center, prospective longitudinal cohort study, including reliability and feasibility in the first 200 participants.

METHODS

COMPASS-ND includes persons with Alzheimer's disease (AD; n = 150), Parkinson's disease (PD) and Lewy body dementias (LBDs) (200), mixed dementia (200), mild cognitive impairment (MCI; 400), subcortical ischemic vascular MCI (V-MCI; 200), subjective cognitive impairment (SCI; 300), and cognitively intact elderly controls (660). Magnetic resonance imaging (MRI) was acquired according to the validated Canadian Dementia Imaging Protocol and visually reviewed by either of two experienced readers blinded to clinical characteristics. Other relevant assessments include history of vascular disease and risk factors, blood pressure, height and weight, cholesterol, glucose, and hemoglobin A1c.

RESULTS

Analyzable data were obtained in 197/200 of whom 18 of whom were clinically diagnosed with V-MCI or mixed dementia. The overall prevalence of infarcts was 24.9%, microbleeds was 24.6%, and high white matter hyperintensity (WMH) was 31.0%. MRI evidence of a potential vascular contribution to neurodegeneration was seen in 12.9%-40.0% of participants clinically diagnosed with another condition such as AD. Inter-rater reliability was good to excellent.

CONCLUSION

COMPASS-ND will be a useful platform to study vascular brain injury and its association with risk factors, biomarkers, and cognitive and functional decline across multiple age-related neurodegenerative diseases. Initial findings show that MRI-defined vascular brain injury is common in all cognitive syndromes and is under-recognized clinically.

Brain network dynamics fingerprints are resilient to data heterogeneity

CONTEXT

Large multi-site neuroimaging datasets have significantly advanced our quest to understand brain-behavior relationships and to develop biomarkers of psychiatric and neurodegenerative disorders. Yet, such data collections come at a cost, as the inevitable differences across samples may lead to biased or erroneous conclusions.

OBJECTIVE

We aim to validate the estimation of individual brain network dynamics fingerprints and appraise sources of variability in large resting-state functional magnetic resonance imaging (rs-fMRI) datasets by providing a novel point of view based on data-driven dynamical models.

APPROACH

Previous work has investigated this critical issue in terms of effects on static measures, such as functional connectivity and brain parcellations. Here, we utilize dynamical models (Hidden Markov models - HMM) to examine how diverse scanning factors in multi-site fMRI recordings affect our ability to infer the brain's spatiotemporal wandering between large-scale networks of activity. Specifically, we leverage a stable HMM trained on the Human Connectome Project (homogeneous) dataset, which we then apply to an heterogeneous dataset of traveling subjects scanned under a multitude of conditions.

MAIN RESULTS

Building upon this premise, we first replicate previous work on the emergence of non-random sequences of brain states. We next highlight how these time-varying brain activity patterns are robust subject-specific fingerprints. Finally, we suggest these fingerprints may be used to assess which scanning factors induce high variability in the data.

SIGNIFICANCE

These results demonstrate that we can i) use large scale dataset to train models that can be then used to interrogate subject-specific data, ii) recover the unique trajectories of brain activity changes in each individual, but also iii) urge caution as our ability to infer such patterns is affected by how, where and when we do so.

A dataset of long-term consistency values of resting-state fMRI connectivity maps in a single individual derived at multiple sites and vendors using the Canadian Dementia Imaging Protocol

The impact of multisite acquisition on resting-state functional MRI (rsfMRI) connectivity has recently gained attention. We provide consistency values (Pearson's correlation) between rsfMRI connectivity maps of an adult volunteer (Csub) scanned 25 times over 3.5 years at 13 sites using the Canadian Dementia Imaging Protocol (CDIP, www.cdip-pcid.ca). This dataset was generated as part of the following article: Multivariate consistency of resting-state fMRI connectivity maps acquired on a single individual over 2.5 years, 13 sites and 3 vendors [1]. Acquired on three 3T scanner vendors (GE, Siemens and Philips), the Csub dataset is part of an ongoing effort to monitor the quality and comparability of MRI data collected across the Canadian Consortium on Neurodegeneration in Aging (CCNA) imaging network. The participant was scanned 25 times in the above-mentioned article: multiple times at six sites over a period of 2.5 years, and once at the remaining seven sites. Since then the participant was scanned an additional 45 times, allowing us to extend the dataset to 70 rsfMRI scans over a period of >4 years.

In addition, we provide intra- and inter-subject consistency values of rsfMRI connectivity maps derived from 26 adult participants belonging to the publicly released Hangzhou Normal University dataset (HNU1). All HNU1 participants underwent 10 rsfMRI scans over one month on a single 3T scanner (GE).

Connectivity maps of seven canonical networks were generated for each scan in the two datasets (Csub and HNU1). All consistency values, along with the scripts used to preprocess the rsfMRI data and generate connectivity maps and pairwise consistency values, have been made available on two public repositories, Github and Zenodo. We have also made available four Jupyter notebooks that use the provided consistency values to (a) generate interactive graphical summaries – 1 notebook, (b) perform statistical analyses - 2 notebooks, and (c) perform data-driven cluster analysis for the recovery of subject identity (i.e. rsfMRI fingerprinting) – 1 notebook. In addition, we provide two interactive dashboards that allow visualization of individual connectivity maps from the two datasets. Finally, we also provide minimally preprocessed rsfMRI data in Brain Imaging Data Standard (BIDS) format on all 70 scans in the extended dataset.

Structural and functional multi-platform MRI series of a single human volunteer over more than fifteen years

We present MRI data from a single human volunteer consisting in over 599 multi-contrast MR images (T1-weighted, T2-weighted, proton density, fluid-attenuated inversion recovery, T2* gradient-echo, diffusion, susceptibility-weighted, arterial-spin labelled, and resting state BOLD functional connectivity imaging) acquired in over 73 sessions on 36 different scanners (13 models, three manufacturers) over the course of 15+ years (cf. Data records). Data included planned data collection acquired within the Consortium pour l'identification précoce de la maladie Alzheimer - Québec (CIMA-Q) and Canadian Consortium on Neurodegeneration in Aging (CCNA) studies, as well as opportunistic data collection from various protocols. These multiple within- and between-centre scans over a substantial time course of a single, cognitively healthy volunteer can be useful to answer a number of methodological questions of interest to the community.