Continuous evaluation of denoising strategies in resting-state fMRI connectivity using fMRIPrep and Nilearn.

Reducing contributions from non-neuronal sources is a crucial step in functional magnetic resonance imaging (fMRI) connectivity analyses. Many viable strategies for denoising fMRI are used in the literature, and practitioners rely on denoising benchmarks for guidance in the selection of an appropriate choice for their study. However, fMRI denoising software is an ever-evolving field, and the benchmarks can quickly become obsolete as the techniques or implementations change. In this work, we present a denoising benchmark featuring a range of denoising strategies, datasets and evaluation metrics for connectivity analyses, based on the popular fMRIprep software. The benchmark prototypes an implementation of a reproducible framework, where the provided Jupyter Book enables readers to reproduce or modify the figures on the Neurolibre reproducible preprint server (https://neurolibre.org/). We demonstrate how such a reproducible benchmark can be used for continuous evaluation of research software, by comparing two versions of the fMRIprep. Most of the benchmark results were consistent with prior literature. Scrubbing, a technique which excludes time points with excessive motion, combined with global signal regression, is generally effective at noise removal. Scrubbing was generally effective, but is incompatible with statistical analyses requiring the continuous sampling of brain signal, for which a simpler strategy, using motion parameters, average activity in select brain compartments, and global signal regression, is preferred. Importantly, we found that certain denoising strategies behave inconsistently across datasets and/or versions of fMRIPrep, or had a different behavior than in previously published benchmarks. This work will hopefully provide useful guidelines for the fMRIprep users community, and highlight the importance of continuous evaluation of research methods.

Open design of a reproducible videogame controller for MRI and MEG.

Videogames are emerging as a promising experimental paradigm in neuroimaging. Acquiring gameplay in a scanner remains challenging due to the lack of a scanner-compatible videogame controller that provides a similar experience to standard, commercial devices. In this paper, we introduce a videogame controller designed for use in the functional magnetic resonance imaging as well as magnetoencephalography. The controller is made exclusively of 3D-printed and commercially available parts. We evaluated the quality of our controller by comparing it to a non-MRI compatible controller that was kept outside the scanner. The comparison of response latencies showed reliable button press accuracies of adequate precision. Comparison of the subjects' motion during fMRI recordings of various tasks showed that the use of our controller did not increase the amount of motion produced compared to a regular MR compatible button press box. Motion levels during an ecological videogame task were of moderate amplitude. In addition, we found that the controller only had marginal effect on temporal SNR in fMRI, as well as on covariance between sensors in MEG, as expected due to the use of non-magnetic building materials. Finally, the reproducibility of the controller was demonstrated by having team members who were not involved in the design build a reproduction using only the documentation. This new videogame controller opens new avenues for ecological tasks in fMRI, including challenging videogames and more generally tasks with complex responses. The detailed controller documentation and build instructions are released under an Open Source Hardware license to increase accessibility, and reproducibility and enable the neuroimaging research community to improve or modify the controller for future experiments.

Brain decoding of the Human Connectome Project tasks in a dense individual fMRI dataset.

Brain decoding aims to infer cognitive states from patterns of brain activity. Substantial inter-individual variations in functional brain organization challenge accurate decoding performed at the group level. In this paper, we tested whether accurate brain decoding models can be trained entirely at the individual level. We trained several classifiers on a dense individual functional magnetic resonance imaging (fMRI) dataset for which six participants completed the entire Human Connectome Project (HCP) task battery >13 times over ten separate fMRI sessions. We evaluated nine decoding methods, from Support Vector Machines (SVM) and Multi-Layer Perceptron (MLP) to Graph Convolutional Neural Networks (GCN). All decoders were trained to classify single fMRI volumes into 21 experimental conditions simultaneously, using ∼7 h of fMRI data per participant. The best prediction accuracies were achieved with GCN and MLP models, whose performance (57-67 % accuracy) approached state-of-the-art accuracy (76 %) with models trained at the group level on >1 K hours of data from the original HCP sample. Our SVM model also performed very well (54-62 % accuracy). Feature importance maps derived from MLP -our best-performing model- revealed informative features in regions relevant to particular cognitive domains, notably in the motor cortex. We also observed that inter-subject classification achieved substantially lower accuracy than subject-specific models, indicating that our decoders learned individual-specific features. This work demonstrates that densely-sampled neuroimaging datasets can be used to train accurate brain decoding models at the individual level. We expect this work to become a useful benchmark for techniques that improve model generalization across multiple subjects and acquisition conditions.

The Canadian Open Neuroscience Platform-An open science framework for the neuroscience community.

The Canadian Open Neuroscience Platform (CONP) takes a multifaceted approach to enabling open neuroscience, aiming to make research, data, and tools accessible to everyone, with the ultimate objective of accelerating discovery. Its core infrastructure is the CONP Portal, a repository with a decentralized design, where datasets and analysis tools across disparate platforms can be browsed, searched, accessed, and shared in accordance with FAIR principles. Another key piece of CONP infrastructure is NeuroLibre, a preprint server capable of creating and hosting executable and fully reproducible scientific publications that embed text, figures, and code. As part of its holistic approach, the CONP has also constructed frameworks and guidance for ethics and data governance, provided support and developed resources to help train the next generation of neuroscientists, and has fostered and grown an engaged community through outreach and communications. In this manuscript, we provide a high-level overview of this multipronged platform and its vision of lowering the barriers to the practice of open neuroscience and yielding the associated benefits for both individual researchers and the wider community.

Continuous Evaluation of Denoising Strategies in Resting-State fMRI Connectivity Using fMRIPrep and Nilearn.

Reducing contributions from non-neuronal sources is a crucial step in functional magnetic resonance imaging (fMRI) connectivity analyses. Many viable strategies for denoising fMRI are used in the literature, and practitioners rely on denoising benchmarks for guidance in the selection of an appropriate choice for their study. However, fMRI denoising software is an ever-evolving field, and the benchmarks can quickly become obsolete as the techniques or implementations change. In this work, we present a denoising benchmark featuring a range of denoising strategies, datasets and evaluation metrics for connectivity analyses, based on the popular fMRIprep software. The benchmark is implemented in a fully reproducible framework, where the provided research objects enable readers to reproduce or modify core computations, as well as the figures of the article using the Jupyter Book project and the Neurolibre reproducible preprint server (https://neurolibre.org/). We demonstrate how such a reproducible benchmark can be used for continuous evaluation of research software, by comparing two versions of the fMRIprep software package. The majority of benchmark results were consistent with prior literature. Scrubbing, a technique which excludes time points with excessive motion, combined with global signal regression, is generally effective at noise removal. Scrubbing however disrupts the continuous sampling of brain images and is incompatible with some statistical analyses, e.g. auto-regressive modeling. In this case, a simple strategy using motion parameters, average activity in select brain compartments, and global signal regression should be preferred. Importantly, we found that certain denoising strategies behave inconsistently across datasets and/or versions of fMRIPrep, or had a different behavior than in previously published benchmarks. This work will hopefully provide useful guidelines for the fMRIprep users community, and highlight the importance of continuous evaluation of research methods. Our reproducible benchmark infrastructure will facilitate such continuous evaluation in the future, and may also be applied broadly to different tools or even research fields.

Genetic Heterogeneity Shapes Brain Connectivity in Psychiatry.

BACKGROUND: Polygenicity and genetic heterogeneity pose great challenges for studying psychiatric conditions. Genetically informed approaches have been implemented in neuroimaging studies to address this issue. However, the effects on functional connectivity of rare and common genetic risks for psychiatric disorders are largely unknown. Our objectives were to estimate and compare the effect sizes on brain connectivity of psychiatric genomic risk factors with various levels of complexity: oligogenic copy number variants (CNVs), multigenic CNVs, and polygenic risk scores (PRSs) as well as idiopathic psychiatric conditions and traits. METHODS: Resting-state functional magnetic resonance imaging data were processed using the same pipeline across 9 datasets. Twenty-nine connectome-wide association studies were performed to characterize the effects of 15 CNVs (1003 carriers), 7 PRSs, 4 idiopathic psychiatric conditions (1022 individuals with autism, schizophrenia, bipolar conditions, or attention-deficit/hyperactivity disorder), and 2 traits (31,424 unaffected control subjects). RESULTS: Effect sizes on connectivity were largest for psychiatric CNVs (estimates: 0.2-0.65 z score), followed by psychiatric conditions (0.15-0.42), neuroticism and fluid intelligence (0.02-0.03), and PRSs (0.01-0.02). Effect sizes of CNVs on connectivity were correlated to their effects on cognition and risk for disease (r = 0.9, p = 5.93 × 10-6). However, effect sizes of CNVs adjusted for the number of genes significantly decreased from small oligogenic to large multigenic CNVs (r = -0.88, p = 8.78 × 10-6). PRSs had disproportionately low effect sizes on connectivity compared with CNVs conferring similar risk for disease. CONCLUSIONS: Heterogeneity and polygenicity affect our ability to detect brain connectivity alterations underlying psychiatric manifestations.

Brain functional connectivity mirrors genetic pleiotropy in psychiatric conditions.

Pleiotropy occurs when a genetic variant influences more than one trait. This is a key property of the genomic architecture of psychiatric disorders and has been observed for rare and common genomic variants. It is reasonable to hypothesize that the microscale genetic overlap (pleiotropy) across psychiatric conditions and cognitive traits may lead to similar overlaps at the macroscale brain level such as large-scale brain functional networks. We took advantage of brain connectivity, measured by resting-state functional MRI to measure the effects of pleiotropy on large-scale brain networks, a putative step from genes to behaviour. We processed nine resting-state functional MRI datasets including 32 726 individuals and computed connectome-wide profiles of seven neuropsychiatric copy-number-variants, five polygenic scores, neuroticism and fluid intelligence as well as four idiopathic psychiatric conditions. Nine out of 19 pairs of conditions and traits showed significant functional connectivity correlations (rFunctional connectivity), which could be explained by previously published levels of genomic (rGenetic) and transcriptomic (rTranscriptomic) correlations with moderate to high concordance: rGenetic-rFunctional connectivity = 0.71 [0.40-0.87] and rTranscriptomic-rFunctional connectivity = 0.83 [0.52; 0.94]. Extending this analysis to functional connectivity profiles associated with rare and common genetic risk showed that 30 out of 136 pairs of connectivity profiles were correlated above chance. These similarities between genetic risks and psychiatric disorders at the connectivity level were mainly driven by the overconnectivity of the thalamus and the somatomotor networks. Our findings suggest a substantial genetic component for shared connectivity profiles across conditions and traits, opening avenues to delineate general mechanisms-amenable to intervention-across psychiatric conditions and genetic risks.

Functional connectivity subtypes associate robustly with ASD diagnosis.

Our understanding of the changes in functional brain organization in autism is hampered by the extensive heterogeneity that characterizes this neurodevelopmental disorder. Data driven clustering offers a straightforward way to decompose autism heterogeneity into subtypes of connectivity and promises an unbiased framework to investigate behavioral symptoms and causative genetic factors. Yet, the robustness and generalizability of functional connectivity subtypes is unknown. Here, we show that a simple hierarchical cluster analysis can robustly relate a given individual and brain network to a connectivity subtype, but that continuous assignments are more robust than discrete ones. We also found that functional connectivity subtypes are moderately associated with the clinical diagnosis of autism, and these associations generalize to independent replication data. We explored systematically 18 different brain networks as we expected them to associate with different behavioral profiles as well as different key regions. Contrary to this prediction, autism functional connectivity subtypes converged on a common topography across different networks, consistent with a compression of the primary gradient of functional brain organization, as previously reported in the literature. Our results support the use of data driven clustering as a reliable data dimensionality reduction technique, where any given dimension only associates moderately with clinical manifestations.

Predicting future cognitive decline from non-brain and multimodal brain imaging data in healthy and pathological aging.

Previous literature has focused on predicting a diagnostic label from structural brain imaging. Since subtle changes in the brain precede a cognitive decline in healthy and pathological aging, our study predicts future decline as a continuous trajectory instead. Here, we tested whether baseline multimodal neuroimaging data improve the prediction of future cognitive decline in healthy and pathological aging. Nonbrain data (demographics, clinical, and neuropsychological scores), structural MRI, and functional connectivity data from OASIS-3 (N = 662; age = 46-96 years) were entered into cross-validated multitarget random forest models to predict future cognitive decline (measured by CDR and MMSE), on average 5.8 years into the future. The analysis was preregistered, and all analysis code is publicly available. Combining non-brain with structural data improved the continuous prediction of future cognitive decline (best test-set performance: R2 = 0.42). Cognitive performance, daily functioning, and subcortical volume drove the performance of our model. Including functional connectivity did not improve predictive accuracy. In the future, the prognosis of age-related cognitive decline may enable earlier and more effective individualized cognitive, pharmacological, and behavioral interventions.

Deep learning models of cognitive processes constrained by human brain connectomes.

Decoding cognitive processes from recordings of brain activity has been an active topic in neuroscience research for decades. Traditional decoding studies focused on pattern classification in specific regions of interest and averaging brain activity over many trials. Recently, brain decoding with graph neural networks has been shown to scale at fine temporal resolution and on the full brain, achieving state-of-the-art performance on the human connectome project benchmark. The reason behind this success is likely the strong inductive connectome prior that enables the integration of distributed patterns of brain activity. Yet, the nature of such inductive bias is still poorly understood. In this work, we investigate the impact of the inclusion of multiple path lengths (through high-order graph convolution), the homogeneity of brain parcels (graph nodes), and the type of interactions (graph edges). We evaluate the decoding models on a large population of 1200 participants, under 21 different experimental conditions, acquired from the Human Connectome Project database. Our findings reveal that the optimal choice for large-scale cognitive decoding is to propagate neural dynamics within empirical functional connectomes and integrate brain dynamics using high-order graph convolutions. In this setting, the model exhibits high decoding accuracy and robustness against adversarial attacks on the graph architecture, including randomization in functional connectomes and lesions in targeted brain regions and networks. The trained model relies on biologically meaningful features for the prediction of cognitive states and generates task-specific graph representations resembling task-evoked activation maps. These results demonstrate that a full-brain integrative model is critical for the large-scale brain decoding. Our study establishes principles of how to effectively leverage human connectome constraints in deep graph neural networks, providing new avenues to study the neural substrates of human cognition at scale.

Data-Driven Analyses of Longitudinal Hippocampal Imaging Trajectories: Discrimination and Biomarker Prediction of Change Classes.

BACKGROUND: Hippocampal atrophy is a well-known biomarker of neurodegeneration, such as that observed in Alzheimer's disease (AD). Although distributions of hippocampal volume trajectories for asymptomatic individuals often reveal substantial heterogeneity, it is unclear whether interpretable trajectory classes can be objectively detected and used for prediction analyses. OBJECTIVE: To detect and predict hippocampal trajectory classes in a computationally competitive context using established AD-related risk factors/biomarkers. METHODS: We used biomarker/risk factor and longitudinal MRI data in asymptomatic adults from the AD Neuroimaging Initiative (n = 351; Mean = 75 years; 48.7% female). First, we applied latent class growth analyses to left (LHC) and right (RHC) hippocampal trajectory distributions to identify distinct classes. Second, using random forest analyses, we tested 38 multi-modal biomarkers/risk factors for their relative importance in discriminating the lower (potentially elevated atrophy risk) from the higher (potentially reduced risk) class. RESULTS: For both LHC and RHC trajectory distribution analyses, we observed three distinct trajectory classes. Three biomarkers/risk factors predicted membership in LHC and RHC lower classes: male sex, higher education, and lower plasma Aß1-42. Four additional factors selectively predicted membership in the lower LHC class: lower plasma tau and Aß1-40, higher depressive symptomology, and lower body mass index. CONCLUSION: Data-driven analyses of LHC and RHC trajectories detected three classes underlying the heterogeneous distributions. Machine learning analyses determined three common and four unique biomarkers/risk factors discriminating the higher and lower LHC/RHC classes. Our sequential analytic approach produced evidence that the dynamics of preclinical hippocampal trajectories can be predicted by AD-related biomarkers/risk factors from multiple modalities.

Open science datasets from PREVENT-AD, a longitudinal cohort of pre-symptomatic Alzheimer's disease.

To move Alzheimer Disease (AD) research forward it is essential to collect data from large cohorts, but also make such data available to the global research community. We describe the creation of an open science dataset from the PREVENT-AD (PResymptomatic EValuation of Experimental or Novel Treatments for AD) cohort, composed of cognitively unimpaired older individuals with a parental or multiple-sibling history of AD. From 2011 to 2017, 386 participants were enrolled (mean age 63 years old ± 5) for sustained investigation among whom 349 have retrospectively agreed to share their data openly. Repositories are findable through the unified interface of the Canadian Open Neuroscience Platform and contain up to five years of longitudinal imaging data, cerebral fluid biochemistry, neurosensory capacities, cognitive, genetic, and medical information. Imaging data can be accessed openly at https://openpreventad.loris.ca while most of the other information, sensitive by nature, is accessible by qualified researchers at https://registeredpreventad.loris.ca. In addition to being a living resource for continued data acquisition, PREVENT-AD offers opportunities to facilitate understanding of AD pathogenesis.

Unraveling reproducible dynamic states of individual brain functional parcellation.

Data-driven parcellations are widely used for exploring the functional organization of the brain, and also for reducing the high dimensionality of fMRI data. Despite the flurry of methods proposed in the literature, functional brain parcellations are not highly reproducible at the level of individual subjects, even with very long acquisitions. Some brain areas are also more difficult to parcellate than others, with association heteromodal cortices being the most challenging. An important limitation of classical parcellations is that they are static, that is, they neglect dynamic reconfigurations of brain networks. In this paper, we proposed a new method to identify dynamic states of parcellations, which we hypothesized would improve reproducibility over static parcellation approaches. For a series of seed voxels in the brain, we applied a cluster analysis to regroup short (3 min) time windows into "states" with highly similar seed parcels. We split individual time series of the Midnight scan club sample into two independent sets of 2.5 hr (test and retest). We found that average within-state parcellations, called stability maps, were highly reproducible (over 0.9 test-retest spatial correlation in many instances) and subject specific (fingerprinting accuracy over 70% on average) between test and retest. Consistent with our hypothesis, seeds in heteromodal cortices (posterior and anterior cingulate) showed a richer repertoire of states than unimodal (visual) cortex. Taken together, our results indicate that static functional parcellations are incorrectly averaging well-defined and distinct dynamic states of brain parcellations. This work calls to revisit previous methods based on static parcellations, which includes the majority of published network analyses of fMRI data. Our method may, thus, impact how researchers model the rich interactions between brain networks in health and disease.

Dissecting autism and schizophrenia through neuroimaging genomics.

Neuroimaging genomic studies of autism spectrum disorder and schizophrenia have mainly adopted a 'top-down' approach, beginning with the behavioural diagnosis, and moving down to intermediate brain phenotypes and underlying genetic factors. Advances in imaging and genomics have been successfully applied to increasingly large case-control studies. As opposed to diagnostic-first approaches, the bottom-up strategy begins at the level of molecular factors enabling the study of mechanisms related to biological risk, irrespective of diagnoses or clinical manifestations. The latter strategy has emerged from questions raised by top-down studies: why are mutations and brain phenotypes over-represented in individuals with a psychiatric diagnosis? Are they related to core symptoms of the disease or to comorbidities? Why are mutations and brain phenotypes associated with several psychiatric diagnoses? Do they impact a single dimension contributing to all diagnoses? In this review, we aimed at summarizing imaging genomic findings in autism and schizophrenia as well as neuropsychiatric variants associated with these conditions. Top-down studies of autism and schizophrenia identified patterns of neuroimaging alterations with small effect-sizes and an extreme polygenic architecture. Genomic variants and neuroimaging patterns are shared across diagnostic categories suggesting pleiotropic mechanisms at the molecular and brain network levels. Although the field is gaining traction; characterizing increasingly reproducible results, it is unlikely that top-down approaches alone will be able to disentangle mechanisms involved in autism or schizophrenia. In stark contrast with top-down approaches, bottom-up studies showed that the effect-sizes of high-risk neuropsychiatric mutations are equally large for neuroimaging and behavioural traits. Low specificity has been perplexing with studies showing that broad classes of genomic variants affect a similar range of behavioural and cognitive dimensions, which may be consistent with the highly polygenic architecture of psychiatric conditions. The surprisingly discordant effect sizes observed between genetic and diagnostic first approaches underscore the necessity to decompose the heterogeneity hindering case-control studies in idiopathic conditions. We propose a systematic investigation across a broad spectrum of neuropsychiatric variants to identify putative latent dimensions underlying idiopathic conditions. Gene expression data on temporal, spatial and cell type organization in the brain have also considerable potential for parsing the mechanisms contributing to these dimensions' phenotypes. While large neuroimaging genomic datasets are now available in unselected populations, there is an urgent need for data on individuals with a range of psychiatric symptoms and high-risk genomic variants. Such efforts together with more standardized methods will improve mechanistically informed predictive modelling for diagnosis and clinical outcomes.

Cerebral functional networks during sleep in young and older individuals.

Even though sleep modification is a hallmark of the aging process, age-related changes in functional connectivity using functional Magnetic Resonance Imaging (fMRI) during sleep, remain unknown. Here, we combined electroencephalography and fMRI to examine functional connectivity differences between wakefulness and light sleep stages (N1 and N2 stages) in 16 young (23.1 ± 3.3y; 7 women), and 14 older individuals (59.6 ± 5.7y; 8 women). Results revealed extended, distributed (inter-between) and local (intra-within) decreases in network connectivity during sleep both in young and older individuals. However, compared to the young participants, older individuals showed lower decreases in connectivity or even increases in connectivity between thalamus/basal ganglia and several cerebral regions as well as between frontal regions of various networks. These findings reflect a reduced ability of the older brain to disconnect during sleep that may impede optimal disengagement for loss of responsiveness, enhanced lighter and fragmented sleep, and contribute to age effects on sleep-dependent brain plasticity.

Functional annotation of human cognitive states using deep graph convolution.

A key goal in neuroscience is to understand brain mechanisms of cognitive functions. An emerging approach is "brain decoding", which consists of inferring a set of experimental conditions performed by a participant, using pattern classification of brain activity. Few works so far have attempted to train a brain decoding model that would generalize across many different cognitive tasks drawn from multiple cognitive domains. To tackle this problem, we proposed a multidomain brain decoder that automatically learns the spatiotemporal dynamics of brain response within a short time window using a deep learning approach. We evaluated the decoding model on a large population of 1200 participants, under 21 different experimental conditions spanning six different cognitive domains, acquired from the Human Connectome Project task-fMRI database. Using a 10s window of fMRI response, the 21 cognitive states were identified with a test accuracy of 90% (chance level 4.8%). Performance remained good when using a 6s window (82%). It was even feasible to decode cognitive states from a single fMRI volume (720ms), with the performance following the shape of the hemodynamic response. Moreover, a saliency map analysis demonstrated that the high decoding performance was driven by the response of biologically meaningful brain regions. Together, we provide an automated tool to annotate human brain activity with fine temporal resolution and fine cognitive granularity. Our model shows potential applications as a reference model for domain adaptation, possibly making contributions in a variety of domains, including neurological and psychiatric disorders.

Functional Connectivity of Successful Picture-Naming: Age-Specific Organization and the Effect of Engaging in Stimulating Activities.

Aging is a lifelong process that starts at birth. Throughout the course of their life, individuals are exposed to various levels of stimulating activities. A higher level of engagement in such activities is suspected to protect against the normal course of cognitive aging or the cognitive manifestations of age-related brain diseases. However, the exact mechanism underlying such protective action remains unclear. The concept of the neurocognitive reserve was introduced to refer to the hypothesis that engagement in stimulating activities shapes brain structure and function, thus indirectly allowing for better preserved cognitive abilities. Although it is known that word production is among the best-preserved cognitive abilities in aging, the underlying neurofunctional mechanisms that allow this relative preservation are still unknown, and it is still unclear how engagement in stimulating activities affects these processes. The objective of this study is to describe the brain functional connectivity patterns associated with picture-naming abilities in younger and older adults with varying levels of engagement in stimulating activities, as a proxy for neurocognitive reserve. A mediation analysis was applied to determine whether the association between reserve proxies and naming accuracy is dependent on task FC. Results show that naming accuracy depends on the posterior cingulate cortex (PCC) functional decoupling in both younger and older adults but through different pathways. While high-performing older adults rely on the asynchronization of this area from motor speech regions' activity, the best-performing younger adults rely on the functional decoupling with language-related regions. Mediation analysis reveals that the PCC decoupling mediates the relationship between the level of engagement in stimulating activities and naming accuracy in younger adults, but not in older adults. These findings suggest that reserve-related mechanisms may be more critical for naming in early adult life, while older adults' neurofunctional organization may benefit more from a lifetime of acquired knowledge.