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.

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

Studies using resting-state functional magnetic resonance imaging (rsfMRI) are increasingly collecting data at multiple sites in order to speed up recruitment or increase sample size. The main objective of this study was to assess the long-term consistency of rsfMRI connectivity maps derived at multiple sites and vendors using the Canadian Dementia Imaging Protocol (CDIP, www.cdip-pcid.ca). Nine to 10 min of functional BOLD images were acquired from an adult cognitively healthy volunteer scanned repeatedly at 13 Canadian sites on three scanner makes (General Electric, Philips and Siemens) over the course of 2.5 years. The consistency (spatial Pearson's correlation) of rsfMRI connectivity maps for seven canonical networks ranged from 0.3 to 0.8, with a negligible effect of time, but significant site and vendor effects. We noted systematic differences in data quality (i.e. head motion, number of useable time frames, temporal signal-to-noise ratio) across vendors, which may also confound some of these results, and could not be disentangled in this sample. We also pooled the long-term longitudinal data with a single-site, short-term (1 month) data sample acquired on 26 subjects (10 scans per subject), called HNU1. Using randomly selected pairs of scans from each subject, we quantified the ability of a data-driven unsupervised cluster analysis to match two scans of the same subjects. In this "fingerprinting" experiment, we found that scans from the Canadian subject (Csub) could be matched with high accuracy intra-site (>95% for some networks), but that the accuracy decreased substantially for scans drawn from different sites and vendors, even falling outside of the range of accuracies observed in HNU1. Overall, our results demonstrate good multivariate stability of rsfMRI measures over several years, but substantial impact of scanning site and vendors. How detrimental these effects are will depend on the application, yet our results demonstrate that new methods for harmonizing multisite analysis represent an important area for future work.

Sex beyond the genitalia: The human brain mosaic.

Whereas a categorical difference in the genitals has always been acknowledged, the question of how far these categories extend into human biology is still not resolved. Documented sex/gender differences in the brain are often taken as support of a sexually dimorphic view of human brains ("female brain" or "male brain"). However, such a distinction would be possible only if sex/gender differences in brain features were highly dimorphic (i.e., little overlap between the forms of these features in males and females) and internally consistent (i.e., a brain has only "male" or only "female" features). Here, analysis of MRIs of more than 1,400 human brains from four datasets reveals extensive overlap between the distributions of females and males for all gray matter, white matter, and connections assessed. Moreover, analyses of internal consistency reveal that brains with features that are consistently at one end of the "maleness-femaleness" continuum are rare. Rather, most brains are comprised of unique "mosaics" of features, some more common in females compared with males, some more common in males compared with females, and some common in both females and males. Our findings are robust across sample, age, type of MRI, and method of analysis. These findings are corroborated by a similar analysis of personality traits, attitudes, interests, and behaviors of more than 5,500 individuals, which reveals that internal consistency is extremely rare. Our study demonstrates that, although there are sex/gender differences in the brain, human brains do not belong to one of two distinct categories: male brain/female brain.

NIDM-Terms: community-based terminology management for improved neuroimaging dataset descriptions and query.

The biomedical research community is motivated to share and reuse data from studies and projects by funding agencies and publishers. Effectively combining and reusing neuroimaging data from publicly available datasets, requires the capability to query across datasets in order to identify cohorts that match both neuroimaging and clinical/behavioral data criteria. Critical barriers to operationalizing such queries include, in part, the broad use of undefined study variables with limited or no annotations that make it difficult to understand the data available without significant interaction with the original authors. Using the Brain Imaging Data Structure (BIDS) to organize neuroimaging data has made querying across studies for specific image types possible at scale. However, in BIDS, beyond file naming and tightly controlled imaging directory structures, there are very few constraints on ancillary variable naming/meaning or experiment-specific metadata. In this work, we present NIDM-Terms, a set of user-friendly terminology management tools and associated software to better manage individual lab terminologies and help with annotating BIDS datasets. Using these tools to annotate BIDS data with a Neuroimaging Data Model (NIDM) semantic web representation, enables queries across datasets to identify cohorts with specific neuroimaging and clinical/behavioral measurements. This manuscript describes the overall informatics structures and demonstrates the use of tools to annotate BIDS datasets to perform integrated cross-cohort queries.

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.

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.

A Standardized Protocol for Efficient and Reliable Quality Control of Brain Registration in Functional MRI Studies.

Automatic alignment of brain anatomy in a standard space is a key step when processing magnetic resonance imaging for group analyses. Such brain registration is prone to failure, and the results are therefore typically reviewed visually to ensure quality. There is however no standard, validated protocol available to perform this visual quality control (QC). We propose here a standardized QC protocol for brain registration, with minimal training overhead and no required knowledge of brain anatomy. We validated the reliability of three-level QC ratings (OK, Maybe, Fail) across different raters. Nine experts each rated N = 100 validation images, and reached moderate to good agreement (kappa from 0.4 to 0.68, average of 0.54 ± 0.08), with the highest agreement for "Fail" images (Dice from 0.67 to 0.93, average of 0.8 ± 0.06). We then recruited volunteers through the Zooniverse crowdsourcing platform, and extracted a consensus panel rating for both the Zooniverse raters (N = 41) and the expert raters. The agreement between expert and Zooniverse panels was high (kappa = 0.76). Overall, our protocol achieved a good reliability when performing a two level assessment (Fail vs. OK/Maybe) by an individual rater, or aggregating multiple three-level ratings (OK, Maybe, Fail) from a panel of experts (3 minimum) or non-experts (15 minimum). Our brain registration QC protocol will help standardize QC practices across laboratories, improve the consistency of reporting of QC in publications, and will open the way for QC assessment of large datasets which could be used to train automated QC systems.

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.

Mutations associated with neuropsychiatric conditions delineate functional brain connectivity dimensions contributing to autism and schizophrenia.

16p11.2 and 22q11.2 Copy Number Variants (CNVs) confer high risk for Autism Spectrum Disorder (ASD), schizophrenia (SZ), and Attention-Deficit-Hyperactivity-Disorder (ADHD), but their impact on functional connectivity (FC) remains unclear. Here we report an analysis of resting-state FC using magnetic resonance imaging data from 101 CNV carriers, 755 individuals with idiopathic ASD, SZ, or ADHD and 1,072 controls. We characterize CNV FC-signatures and use them to identify dimensions contributing to complex idiopathic conditions. CNVs have large mirror effects on FC at the global and regional level. Thalamus, somatomotor, and posterior insula regions play a critical role in dysconnectivity shared across deletions, duplications, idiopathic ASD, SZ but not ADHD. Individuals with higher similarity to deletion FC-signatures exhibit worse cognitive and behavioral symptoms. Deletion similarities identified at the connectivity level could be related to the redundant associations observed genome-wide between gene expression spatial patterns and FC-signatures. Results may explain why many CNVs affect a similar range of neuropsychiatric symptoms.

Analysis of task-based functional MRI data preprocessed with fMRIPrep.

Functional magnetic resonance imaging (fMRI) is a standard tool to investigate the neural correlates of cognition. fMRI noninvasively measures brain activity, allowing identification of patterns evoked by tasks performed during scanning. Despite the long history of this technique, the idiosyncrasies of each dataset have led to the use of ad-hoc preprocessing protocols customized for nearly every different study. This approach is time consuming, error prone and unsuitable for combining datasets from many sources. Here we showcase fMRIPrep (http://fmriprep.org), a robust tool to prepare human fMRI data for statistical analysis. This software instrument addresses the reproducibility concerns of the established protocols for fMRI preprocessing. By leveraging the Brain Imaging Data Structure to standardize both the input datasets (MRI data as stored by the scanner) and the outputs (data ready for modeling and analysis), fMRIPrep is capable of preprocessing a diversity of datasets without manual intervention. In support of the growing popularity of fMRIPrep, this protocol describes how to integrate the tool in a task-based fMRI investigation workflow.

A highly predictive signature of cognition and brain atrophy for progression to Alzheimer's dementia.

BACKGROUND: Clinical trials in Alzheimer's disease need to enroll patients whose cognition will decline over time, if left untreated, in order to demonstrate the efficacy of an intervention. Machine learning models used to screen for patients at risk of progression to dementia should therefore favor specificity (detecting only progressors) over sensitivity (detecting all progressors), especially when the prevalence of progressors is low. Here, we explore whether such high-risk patients can be identified using cognitive assessments and structural neuroimaging by training machine learning tools in a high-specificity regime. RESULTS: A multimodal signature of Alzheimer's dementia was first extracted from the ADNI1 dataset. We then validated the predictive value of this signature on ADNI1 patients with mild cognitive impairment (N = 235). The signature was optimized to predict progression to dementia over 3 years with low sensitivity (55.1%) but high specificity (95.6%), resulting in only moderate accuracy (69.3%) but high positive predictive value (80.4%, adjusted for a "typical" 33% prevalence rate of true progressors). These results were replicated in ADNI2 (N = 235), with 87.8% adjusted positive predictive value (96.7% specificity, 47.3% sensitivity, 85.1% accuracy). CONCLUSIONS: We found that cognitive measures alone could identify high-risk individuals, with structural measurements providing a slight improvement. The signature had comparable receiver operating characteristics to standard machine learning tools, yet a marked improvement in positive predictive value was achieved over the literature by selecting a high-specificity operating point. The multimodal signature can be readily applied for the enrichment of clinical trials.

A correspondence between individual differences in the brain's intrinsic functional architecture and the content and form of self-generated thoughts.

Although neural activity often reflects the processing of external inputs, intrinsic fluctuations in activity have been observed throughout the brain. These may relate to patterns of self-generated thought that can occur while not performing goal-driven tasks. To understand the relationship between self-generated mental activity and intrinsic neural fluctuations, we developed the New York Cognition Questionnaire (NYC-Q) to assess the content and form of an individual's experiences during the acquisition of resting-state fMRI data. The data were collected as a part of the Nathan Kline Rockland Enhanced sample. We decomposed NYC-Q scores using exploratory factor analysis and found that self-reported thoughts clustered into distinct dimensions of content (future related, past related, positive, negative, and social) and form (words, images, and specificity). We used these components to perform an individual difference analysis exploring how differences in the types of self-generated thoughts relate to whole brain measures of intrinsic brain activity (fractional amplitude of low frequency fluctuations, regional homogeneity, and degree centrality). We found patterns of self-generated thoughts related to changes that were distributed across a wide range of cortical areas. For example, individuals who reported greater imagery exhibited greater low frequency fluctuations in a region of perigenual cingulate cortex, a region that is known to participate in the so-called default-mode network. We also found certain forms of thought were associated with other areas, such as primary visual cortex, the insula, and the cerebellum. For example, individuals who reported greater future thought exhibited less homogeneous neural fluctuations in a region of lateral occipital cortex, a result that is consistent with the claim that particular types of self-generated thought depend on processes that are decoupled from sensory processes. These data provide evidence that self-generated thought is a heterogeneous category of experience and that studying its content can be helpful in understanding brain dynamics.