Although it is being successfully implemented for exploration of the genome, discovery science has eluded the functional neuroimaging community. The core challenge remains the development of common paradigms for interrogating the myriad functional systems in the brain without the constraints of a priori hypotheses. Resting-state functional MRI (R-fMRI) constitutes a candidate approach capable of addressing this challenge. Imaging the brain during rest reveals large-amplitude spontaneous low-frequency (<0.1 Hz) fluctuations in the fMRI signal that are temporally correlated across functionally related areas. Referred to as functional connectivity, these correlations yield detailed maps of complex neural systems, collectively constituting an individual's "functional connectome." Reproducibility across datasets and individuals suggests the functional connectome has a common architecture, yet each individual's functional connectome exhibits unique features, with stable, meaningful interindividual differences in connectivity patterns and strengths. Comprehensive mapping of the functional connectome, and its subsequent exploitation to discern genetic influences and brain-behavior relationships, will require multicenter collaborative datasets. Here we initiate this endeavor by gathering R-fMRI data from 1,414 volunteers collected independently at 35 international centers. We demonstrate a universal architecture of positive and negative functional connections, as well as consistent loci of inter-individual variability. Age and sex emerged as significant determinants. These results demonstrate that independent R-fMRI datasets can be aggregated and shared. Highthroughput R-fMRI can provide quantitative phenotypes for molecular genetic studies and biomarkers of developmental and pathological processes in the brain. To initiate discovery science of brain function, the 1000 Functional Connectomes Project dataset is freely accessible at www.nitrc.org/projects/fcon_1000/.

Graph-based computational network analysis has proven a powerful tool to quantitatively characterize functional architectures of the brain. However, the test-retest (TRT) reliability of graph metrics of functional networks has not been systematically examined. Here, we investigated TRT reliability of topological metrics of functional brain networks derived from resting-state functional magnetic resonance imaging data. Specifically, we evaluated both short-term (<1 hour apart) and long-term (>5 months apart) TRT reliability for 12 global and 6 local nodal network metrics. We found that reliability of global network metrics was overall low, threshold-sensitive and dependent on several factors of scanning time interval (TI, long-term>short-term), network membership (NM, networks excluding negative correlations>networks including negative correlations) and network type (NT, binarized networks>weighted networks). The dependence was modulated by another factor of node definition (ND) strategy. The local nodal reliability exhibited large variability across nodal metrics and a spatially heterogeneous distribution. Nodal degree was the most reliable metric and varied the least across the factors above. Hub regions in association and limbic/paralimbic cortices showed moderate TRT reliability. Importantly, nodal reliability was robust to above-mentioned four factors. Simulation analysis revealed that global network metrics were extremely sensitive (but varying degrees) to noise in functional connectivity and weighted networks generated numerically more reliable results in compared with binarized networks. For nodal network metrics, they showed high resistance to noise in functional connectivity and no NT related differences were found in the resistance. These findings provide important implications on how to choose reliable analytical schemes and network metrics of interest.

Studies of resting-state fMRI have shown that blood oxygen level dependent (BOLD) signals giving rise to temporal correlation across voxels (or regions) are dominated by low-frequency fluctuations in the range of *0.01-0.1 Hz. These low-frequency fluctuations have been further divided into multiple distinct frequency bands (slow-5 and -4) based on earlier neurophysiological studies, though low sampling frequency of fMRI (*0.5 Hz) has substantially limited the exploration of other known frequency bands of neurophysiological origins (slow-3, -2, and -1). In this study, we used resting-state fMRI data acquired from 21 healthy subjects at a higher sampling frequency of 1.5 Hz to assess the presence of resting-state functional connectivity (RSFC) across multiple frequency bands: slow-5 to slow-1. The effect of different frequency bands on spatial extent and connectivity strength for known resting-state networks (RSNs) was also evaluated. RSNs were derived using independent component analysis and seed-based correlation. Commonly known RSNs, such as the default mode, the fronto-parietal, the dorsal attention, and the visual networks, were consistently observed at multiple frequency bands. Significant inter-hemispheric connectivity was observed between each seed and its contra lateral brain region across all frequency bands, though overall spatial extent of seed-based correlation maps decreased in slow-2 and slow-1 frequency bands. These results suggest that functional integration between brain regions at rest occurs over multiple frequency bands and RSFC is a multiband phenomenon. These results also suggest that further investigation of BOLD signal in multiple frequency bands for related cognitive processes should be undertaken.

There is a growing interest in studies of human brain networks using resting-state functional magnetic resonance imaging (fMRI). However, it is unclear whether and how brain networks measured during the resting-state exhibit comparable properties to brain networks during task performance. In the present study, we investigated meta-analytic coactivation patterns among brain regions based upon published neuroimaging studies, and compared the coactivation network configurations with those in the resting-state network. The strength of resting-state functional connectivity between two regions were strongly correlated with the coactivation strength. However, the coactivation network showed greater global efficiency, smaller mean clustering coefficient, and lower modularity compared with the resting-state network, which suggest a more efficient global information transmission and between system integrations during task performing. Hub shifts were also observed within the thalamus and the left inferior temporal cortex. The thalamus and the left inferior temporal cortex exhibited higher and lower degrees, respectively in the coactivation network compared with the resting-state network. These results shed light regarding the reconfiguration of the brain networks between task and resting-state conditions, and highlight the role of the thalamus in change of network configurations in task vs. rest.

Various studies have indicated that the thalamus is involved in controlling both cortico-cortical information flow and cortical communication with the rest of the brain. Detailed anatomy and functional connectivity patterns of the thalamocortical system are essential to understanding the cortical organization and pathophysiology of a wide range of thalamus-related neurological and neuropsychiatric diseases. The current study used resting-state fMRI to investigate the topography of the human thalamocortical system from a functional perspective. The thalamus-related cortical networks were identified by performing independent component analysis on voxel-based thalamic functional connectivity maps across a large group of subjects. The resulting functional brain networks were very similar to well-established resting-state network maps. Using these brain network components in a spatial regression model with each thalamic voxel's functional connectivity map, we localized the thalamic subdivisions related to each brain network. For instance, the medial dorsal nucleus was shown to be associated with the default mode, the bilateral executive, the medial visual networks; and the pulvinar nucleus was involved in both the dorsal attention and the visual networks. These results revealed that a single nucleus may have functional connections with multiple cortical regions or even multiple functional networks, and may be potentially related to the function of mediation or modulation of multiple cortical networks. This observed organization of thalamocortical system provided a reference for studying the functions of thalamic sub-regions. The importance of intrinsic connectivity-based mapping of the thalamocortical relationship is discussed, as well as the applicability of the approach for future studies.