Cold, clumpy accretion onto an active supermassive black hole.

Supermassive black holes in galaxy centres can grow by the accretion of gas, liberating energy that might regulate star formation on galaxy-wide scales. The nature of the gaseous fuel reservoirs that power black hole growth is nevertheless largely unconstrained by observations, and is instead routinely simplified as a smooth, spherical inflow of very hot gas. Recent theory and simulations instead predict that accretion can be dominated by a stochastic, clumpy distribution of very cold molecular clouds--a departure from the 'hot mode' accretion model--although unambiguous observational support for this prediction remains elusive. Here we report observations that reveal a cold, clumpy accretion flow towards a supermassive black hole fuel reservoir in the nucleus of the Abell 2597 Brightest Cluster Galaxy (BCG), a nearby (redshift z = 0.0821) giant elliptical galaxy surrounded by a dense halo of hot plasma. Under the right conditions, thermal instabilities produce a rain of cold clouds that fall towards the galaxy's centre, sustaining star formation amid a kiloparsec-scale molecular nebula that is found at its core. The observations show that these cold clouds also fuel black hole accretion, revealing 'shadows' cast by the molecular clouds as they move inward at about 300 kilometres per second towards the active supermassive black hole, which serves as a bright backlight. Corroborating evidence from prior observations of warmer atomic gas at extremely high spatial resolution, along with simple arguments based on geometry and probability, indicate that these clouds are within the innermost hundred parsecs of the black hole, and falling closer towards it.

Test-retest reliability of dynamic functional connectivity in resting state fMRI.

While static functional connectivity (sFC) of resting state fMRI (rfMRI) measures the average functional connectivity (FC) over the entire rfMRI scan, dynamic FC (dFC) captures the temporal variations of FC at shorter time windows. Although numerous studies have implemented dFC analyses, only a few studies have investigated the reliability of dFC and this limits the biological interpretation of dFC. Here, we used a large cohort (N = 820) of subjects and four rfMRI scans from the Human Connectome Project to systematically explore the relationship between sFC, dFC and their test-retest reliabilities through intra-class correlation (ICC). dFC ICC was explored through the sliding window approach with three dFC statistics (standard deviation, ALFF, and excursion). Excursion demonstrated the highest dFC ICC and the highest age prediction accuracy. dFC ICC was generally higher at window sizes less than 40 s. sFC and dFC were negatively correlated. Compared to sFC, dFC was less reliable. While sFC and sFC ICC were positively correlated, dFC and dFC ICC were negatively correlated, indicating that FC that was more dynamic was less reliable. Intra-network FCs in the frontal-parietal, default mode, sensorimotor and visual networks demonstrated high sFC and low dFC. Moreover, ICCs of both sFC and dFC in these regions were higher. The above results were consistent across two brain atlases and independent component analysis-based networks, multiple window sizes and all three dFC statistics. In summary, dFC is less reliable than sFC and additional experiments are required to better understand the neurophysiological relevance of dFC.

Functional connectivity predicts gender: Evidence for gender differences in resting brain connectivity.

Prevalence of certain forms of psychopathology, such as autism and depression, differs between genders and understanding gender differences of the neurotypical brain may provide insights into risk and protective factors. In recent research, resting state functional magnetic resonance imaging (rfMRI) is widely used to map the inherent functional networks of the brain. Although previous studies have reported gender differences in rfMRI, the robustness of gender differences is not well characterized. In this study, we use a large data set to test whether rfMRI functional connectivity (FC) can be used to predict gender and identify FC features that are most predictive of gender. We utilized rfMRI data from 820 healthy controls from the Human Connectome Project. By applying a predefined functional template and partial least squares regression modeling, we achieved a gender prediction accuracy of 87% when multi-run rfMRI was used. Permutation tests confirmed that gender prediction was reliable ( p<.001). Effects of motion, age, handedness, blood pressure, weight, and brain volume on gender prediction are discussed. Further, we found that FC features within the default mode (DMN), fronto-parietal and sensorimotor networks contributed most to gender prediction. In the DMN, right fusiform gyrus and right ventromedial prefrontal cortex were important contributors. The above regions have been previously implicated in aspects of social functioning and this suggests potential gender differences in social cognition mediated by the DMN. Our findings demonstrate that gender can be reliably predicted using rfMRI data and highlight the importance of controlling for gender in brain imaging studies.

Leveraging Clinical Imaging Archives for Radiomics: Reliability of Automated Methods for Brain Volume Measurement.

Purpose To validate the use of thick-section clinically acquired magnetic resonance (MR) imaging data for estimating total brain volume (TBV), gray matter (GM) volume (GMV), and white matter (WM) volume (WMV) by using three widely used automated toolboxes: SPM ( www.fil.ion.ucl.ac.uk/spm/ ), FreeSurfer ( surfer.nmr.mgh.harvard.edu ), and FSL (FMRIB software library; Oxford Centre for Functional MR Imaging of the Brain, Oxford, England, https://fsl.fmrib.ox.ac.uk/fsl ). Materials and Methods MR images from a clinical archive were used and data were deidentified. The three methods were applied to estimate brain volumes from thin-section research-quality brain MR images and routine thick-section clinical MR images acquired from the same 38 patients (age range, 1-71 years; mean age, 22 years; 11 women). By using these automated methods, TBV, GMV, and WMV were estimated. Thin- versus thick-section volume comparisons were made for each method by using intraclass correlation coefficients (ICCs). Results SPM exhibited excellent ICCs (0.97, 0.85, and 0.83 for TBV, GMV, and WMV, respectively). FSL exhibited ICCs of 0.69, 0.51, and 0.60 for TBV, GMV, and WMV, respectively, but they were lower than with SPM. FreeSurfer exhibited excellent ICC of 0.63 only for TBV. Application of SPM's voxel-based morphometry on the modulated images of thin-section images and interpolated thick-section images showed fair to excellent ICCs (0.37-0.98) for the majority of brain regions (88.47% [306924 of 346916 voxels] of WM and 80.35% [377 282 of 469 502 voxels] of GM). Conclusion Thick-section clinical-quality MR images can be reliably used for computing quantitative brain metrics such as TBV, GMV, and WMV by using SPM. © RSNA, 2017 Online supplemental material is available for this article.

Wavelet-based fMRI analysis: 3-D denoising, signal separation, and validation metrics.

We present a novel integrated wavelet-domain based framework (w-ICA) for 3-D denoising functional magnetic resonance imaging (fMRI) data followed by source separation analysis using independent component analysis (ICA) in the wavelet domain. We propose the idea of a 3-D wavelet-based multi-directional denoising scheme where each volume in a 4-D fMRI data set is sub-sampled using the axial, sagittal and coronal geometries to obtain three different slice-by-slice representations of the same data. The filtered intensity value of an arbitrary voxel is computed as an expected value of the denoised wavelet coefficients corresponding to the three viewing geometries for each sub-band. This results in a robust set of denoised wavelet coefficients for each voxel. Given the de-correlated nature of these denoised wavelet coefficients, it is possible to obtain more accurate source estimates using ICA in the wavelet domain. The contributions of this work can be realized as two modules: First, in the analysis module we combine a new 3-D wavelet denoising approach with signal separation properties of ICA in the wavelet domain. This step helps obtain an activation component that corresponds closely to the true underlying signal, which is maximally independent with respect to other components. Second, we propose and describe two novel shape metrics for post-ICA comparisons between activation regions obtained through different frameworks. We verified our method using simulated as well as real fMRI data and compared our results against the conventional scheme (Gaussian smoothing+spatial ICA: s-ICA). The results show significant improvements based on two important features: (1) preservation of shape of the activation region (shape metrics) and (2) receiver operating characteristic curves. It was observed that the proposed framework was able to preserve the actual activation shape in a consistent manner even for very high noise levels in addition to significant reduction in false positive voxels.

WITHDRAWN: Erratum to "Does function follow form?: Methods to fuse structural and functional brain images show decreased linkage in schizophrenia" [NeuroImage 49 (2010) 2626-2637].

This article has been withdrawn at the request of the author(s) and/or editor. The Publisher apologizes for any inconvenience this may cause. The full Elsevier Policy on Article Withdrawal can be found at http://www.elsevier.com/locate/withdrawalpolicy.

Does function follow form?: methods to fuse structural and functional brain images show decreased linkage in schizophrenia.

When both structural magnetic resonance imaging (sMRI) and functional MRI (fMRI) data are collected they are typically analyzed separately and the joint information is not examined. Techniques that examine joint information can help to find hidden traits in complex disorders such as schizophrenia. The brain is vastly interconnected, and local brain morphology may influence functional activity at distant regions. In this paper we introduce three methods to identify inter-correlations among sMRI and fMRI voxels within the whole brain. We apply these methods to examine sMRI gray matter data and fMRI data derived from an auditory sensorimotor task from a large study of schizophrenia. In Method 1 the sMRI-fMRI cross-correlation matrix is reduced to a histogram and results show that healthy controls (HC) have stronger correlations than do patients with schizophrenia (SZ). In Method 2 the spatial information of sMRI-fMRI correlations is retained. Structural regions in the cerebellum and frontal regions show more positive and more negative correlations, respectively, with functional regions in HC than in SZ. In Method 3 significant sMRI-fMRI inter-regional links are detected, with regions in the cerebellum showing more significant positive correlations with functional regions in HC relative to SZ. Results from all three methods indicate that the linkage between gray matter and functional activation is stronger in HC than SZ. The methods introduced can be easily extended to comprehensively correlate large data sets.

A method to fuse fMRI tasks through spatial correlations: applied to schizophrenia.

Single task analysis methods of functional MRI brain data, though useful, are not able to evaluate the joint information between tasks. Data fusion of multiple tasks that probe different cognitive processes provides knowledge of the joint information and may be important in order to better understand complex disorders such as schizophrenia. In this article, we introduce a simple but effective technique to fuse two tasks by computing the histogram of correlations for all possible combinations of whole brain voxels. The approach was applied to data derived from healthy controls and patients with schizophrenia from four different tasks, auditory oddball (target), auditory oddball (novel), Sternberg working memory, and sensorimotor. It was found that in four out of six task combinations patients' intertask correlations were more positively correlated than controls', in one combination the controls showed more positive correlations and in another there was no significant difference. The robustness of this result was checked with several testing techniques. The four task combinations for which patients had more positive correlation occurred at different scanning sessions and the task combination that showed the opposite result occurred within the same scanning session. Brain regions that showed high intertask correlations were found for both groups and regions that correlated differently between the two groups were identified. The approach introduced finds interesting results and new differential features that cannot be achieved through traditional methods.

Divide and Conquer: Sub-Grouping of ASD Improves ASD Detection Based on Brain Morphometry.

Low success (<60%) in autism spectrum disorder (ASD) classification using brain morphometry from the large multi-site ABIDE dataset and inconsistent findings on brain morphometric abnormalities in ASD can be attributed to the ASD heterogeneity. In this study, we show that ASD brain morphometry is highly heterogeneous, and demonstrate that the heterogeneity can be mitigated and classification improved if autism severity (AS), verbal IQ (VIQ) and age are used with morphometric features. Morphometric features from structural MRIs (sMRIs) of 734 males (ASD: 361, controls: 373) of ABIDE were derived using FreeSurfer. Applying the Random Forest classifier, an AUC of 0.61 was achieved. Adding VIQ and age to morphometric features, AUC improved to 0.68. Sub-grouping the subjects by AS, VIQ and age improved the classification with the highest AUC of 0.8 in the moderate-AS sub-group (AS = 7-8). Matching subjects on age and/or VIQ in each sub-group further improved the classification with the highest AUC of 0.92 in the low AS sub-group (AS = 4-5). AUC decreased with AS and VIQ, and was the lowest in the mid-age sub-group (13-18 years). The important features were mainly from the frontal, temporal, ventricular, right hippocampal and left amygdala regions. However, they highly varied with AS, VIQ and age. The curvature and folding index features from frontal, temporal, lingual and insular regions were dominant in younger subjects suggesting their importance for early detection. When the experiments were repeated using the Gradient Boosting classifier similar results were obtained. Our findings suggest that identifying brain biomarkers in sub-groups of ASD can yield more robust and insightful results than searching across the whole spectrum. Further, it may allow identification of sub-group specific brain biomarkers that are optimized for early detection and monitoring, increasing the utility of sMRI as an important tool for early detection of ASD.

Approaches to Capture Variance Differences in Rest fMRI Networks in the Spatial Geometric Features: Application to Schizophrenia.

Identification of functionally connected regions while at rest has been at the forefront of research focusing on understanding interactions between different brain regions. Studies have utilized a variety of approaches including seed based as well as data-driven approaches to identifying such networks. Most such techniques involve differentiating groups based on group mean measures. There has been little work focused on differences in spatial characteristics of resting fMRI data. We present a method to identify between group differences in the variability in the cluster characteristics of network regions within components estimated via independent vector analysis (IVA). IVA is a blind source separation approach shown to perform well in capturing individual subject variability within a group model. We evaluate performance of the approach using simulations and then apply to a relatively large schizophrenia data set (82 schizophrenia patients and 89 healthy controls). We postulate, that group differences in the intra-network distributional characteristics of resting state network voxel intensities might indirectly capture important distinctions between the brain function of healthy and clinical populations. Results demonstrate that specific areas of the brain, superior, and middle temporal gyrus that are involved in language and recognition of emotions, show greater component level variance in amplitude weights for schizophrenia patients than healthy controls. Statistically significant correlation between component level spatial variance and component volume was observed in 19 of the 27 non-artifactual components implying an evident relationship between the two parameters. Additionally, the greater spread in the distance of the cluster peak of a component from the centroid in schizophrenia patients compared to healthy controls was observed for seven components. These results indicate that there is hidden potential in exploring variance and possibly higher-order measures in resting state networks to better understand diseases such as schizophrenia. It furthers comprehension of how spatial characteristics can highlight previously unexplored differences between populations such as schizophrenia patients and healthy controls.

Sex and Age Effects of Functional Connectivity in Early Adulthood.

Functional connectivity (FC) in resting-state functional magnetic resonance imaging (rs-fMRI) is widely used to find coactivating regions in the human brain. Despite its widespread use, the effects of sex and age on resting FC are not well characterized, especially during early adulthood. Here we apply regression and graph theoretical analyses to explore the effects of sex and age on FC between the 116 AAL atlas parcellations (a total of 6670 FC measures). rs-fMRI data of 494 healthy subjects (203 males and 291 females; age range: 22-36 years) from the Human Connectome Project were analyzed. We report the following findings. (1) Males exhibited greater FC than females in 1352 FC measures (1025 survived Bonferroni correction; [Formula: see text]). In 641 FC measures, females exhibited greater FC than males but none survived Bonferroni correction. Significant FC differences were mainly present in frontal, parietal, and temporal lobes. Although the average FC values for males and females were significantly different, FC values of males and females exhibited large overlap. (2) Age effects were present only in 29 FC measures and all significant age effects showed higher FC in younger subjects. Age and sex differences of FC remained significant after controlling for cognitive measures. (3) Although sex [Formula: see text] age interaction did not survive multiple comparison correction, FC in females exhibited a faster cross-sectional decline with age. (4) Male brains were more locally clustered in all lobes but the cerebellum; female brains had a higher clustering coefficient at the whole-brain level. Our results indicate that although both male and female brains show small-world network characteristics, male brains were more segregated and female brains were more integrated. Findings of this study further our understanding of FC in early adulthood and provide evidence to support that age and sex should be controlled for in FC studies of young adults.

Spatial Variance in Resting fMRI Networks of Schizophrenia Patients: An Independent Vector Analysis.

Spatial variability in resting functional MRI (fMRI) brain networks has not been well studied in schizophrenia, a disease known for both neurodevelopmental and widespread anatomic changes. Motivated by abundant evidence of neuroanatomical variability from previous studies of schizophrenia, we draw upon a relatively new approach called independent vector analysis (IVA) to assess this variability in resting fMRI networks. IVA is a blind-source separation algorithm, which segregates fMRI data into temporally coherent but spatially independent networks and has been shown to be especially good at capturing spatial variability among subjects in the extracted networks. We introduce several new ways to quantify differences in variability of IVA-derived networks between schizophrenia patients (SZs = 82) and healthy controls (HCs = 89). Voxelwise amplitude analyses showed significant group differences in the spatial maps of auditory cortex, the basal ganglia, the sensorimotor network, and visual cortex. Tests for differences (HC-SZ) in the spatial variability maps suggest, that at rest, SZs exhibit more activity within externally focused sensory and integrative network and less activity in the default mode network thought to be related to internal reflection. Additionally, tests for difference of variance between groups further emphasize that SZs exhibit greater network variability. These results, consistent with our prediction of increased spatial variability within SZs, enhance our understanding of the disease and suggest that it is not just the amplitude of connectivity that is different in schizophrenia, but also the consistency in spatial connectivity patterns across subjects.

Inter-Method Discrepancies in Brain Volume Estimation May Drive Inconsistent Findings in Autism.

Previous studies applying automatic preprocessing methods on Structural Magnetic Resonance Imaging (sMRI) report inconsistent neuroanatomical abnormalities in Autism Spectrum Disorder (ASD). In this study we investigate inter-method differences as a possible cause behind these inconsistent findings. In particular, we focus on the estimation of the following brain volumes: gray matter (GM), white matter (WM), cerebrospinal fluid (CSF), and total intra cranial volume (TIV). T1-weighted sMRIs of 417 ASD subjects and 459 typically developing controls (TDC) from the ABIDE dataset were estimated using three popular preprocessing methods: SPM, FSL, and FreeSurfer (FS). Brain volumes estimated by the three methods were correlated but had significant inter-method differences; except TIVSPM vs. TIVFS, all inter-method differences were significant. ASD vs. TDC group differences in all brain volume estimates were dependent on the method used. SPM showed that TIV, GM, and CSF volumes of ASD were larger than TDC with statistical significance, whereas FS and FSL did not show significant differences in any of the volumes; in some cases, the direction of the differences were opposite to SPM. When methods were compared with each other, they showed differential biases for autism, and several biases were larger than ASD vs. TDC differences of the respective methods. After manual inspection, we found inter-method segmentation mismatches in the cerebellum, sub-cortical structures, and inter-sulcal CSF. In addition, to validate automated TIV estimates we performed manual segmentation on a subset of subjects. Results indicate that SPM estimates are closest to manual segmentation, followed by FS while FSL estimates were significantly lower. In summary, we show that ASD vs. TDC brain volume differences are method dependent and that these inter-method discrepancies can contribute to inconsistent neuroimaging findings in general. We suggest cross-validation across methods and emphasize the need to develop better methods to increase the robustness of neuroimaging findings.

The predictive power of structural MRI in Autism diagnosis.

Diagnosis of Autism Spectrum Disorder (ASD) using structural magnetic resonance imaging (sMRI) of the brain has been a topic of significant research interest. Previous studies using small datasets with well-matched Typically Developing Controls (TDC) report high classification accuracies (80-96%) but studies using the large heterogeneous ABIDE dataset report accuracies less than 60%. In this study we investigate the predictive power of sMRI in ASD using 373 ASD and 361 TDC male subjects from the ABIDE. Brain morphometric features were derived using FreeSurfer and classification was performed using three different techniques: Random Forest (RF), Support Vector Machine (SVM) and Gradient Boosting Machine (GBM). Although high classification accuracies were possible in individual sites (with a maximum of 97% in Caltech), the highest classification accuracy across all sites was only 60% (sensitivity = 57%, specificity = 64%). However, the accuracy across all sites improved to 67% when IQ and age information were added to morphometric features. Across all three classifiers, volume and surface area had more discriminative power. In general, important features for classification were present in the frontal and temporal regions and these regions have been implicated in ASD. This study also explores the effect of demographics and behavioral measures on the predictive power of sMRI. Results show that classification accuracy increases with autism severity and that ASD detection with sMRI is easier before the age of 10 years.

Enhanced disease characterization through multi network functional normalization in fMRI.

Conventionally, structural topology is used for spatial normalization during the pre-processing of fMRI. The co-existence of multiple intrinsic networks which can be detected in the resting brain are well-studied. Also, these networks exhibit temporal and spatial modulation during cognitive task vs. rest which shows the existence of common spatial excitation patterns between these identified networks. Previous work (Khullar et al., 2011) has shown that structural and functional data may not have direct one-to-one correspondence and functional activation patterns in a well-defined structural region can vary across subjects even for a well-defined functional task. The results of this study and the existence of the neural activity patterns in multiple networks motivates us to investigate multiple resting-state networks as a single fusion template for functional normalization for multi groups of subjects. We extend the previous approach (Khullar et al., 2011) by co-registering multi group of subjects (healthy control and schizophrenia patients) and by utilizing multiple resting-state networks (instead of just one) as a single fusion template for functional normalization. In this paper we describe the initial steps toward using multiple resting-state networks as a single fusion template for functional normalization. A simple wavelet-based image fusion approach is presented in order to evaluate the feasibility of combining multiple functional networks. Our results showed improvements in both the significance of group statistics (healthy control and schizophrenia patients) and the spatial extent of activation when a multiple resting-state network applied as a single fusion template for functional normalization after the conventional structural normalization. Also, our results provided evidence that the improvement in significance of group statistics lead to better accuracy results for classification of healthy controls and schizophrenia patients.

ICA-fNORM: Spatial Normalization of fMRI Data Using Intrinsic Group-ICA Networks.

A common pre-processing challenge associated with group level fMRI analysis is spatial registration of multiple subjects to a standard space. Spatial normalization, using a reference image such as the Montreal Neurological Institute brain template, is the most common technique currently in use to achieve spatial congruence across multiple subjects. This method corrects for global shape differences preserving regional asymmetries, but does not account for functional differences. We propose a novel approach to co-register task-based fMRI data using resting state group-ICA networks. We posit that these intrinsic networks (INs) can provide to the spatial normalization process with important information about how each individual's brain is organized functionally. The algorithm is initiated by the extraction of single subject representations of INs using group level independent component analysis (ICA) on resting state fMRI data. In this proof-of-concept work two of the robust, commonly identified, networks are chosen as functional templates. As an estimation step, the relevant INs are utilized to derive a set of normalization parameters for each subject. Finally, the normalization parameters are applied individually to a different set of fMRI data acquired while the subjects performed an auditory oddball task. These normalization parameters, although derived using rest data, generalize successfully to data obtained with a cognitive paradigm for each subject. The improvement in results is verified using two widely applied fMRI analysis methods: the general linear model and ICA. Resulting activation patterns from each analysis method show significant improvements in terms of detection sensitivity and statistical significance at the group level. The results presented in this article provide initial evidence to show that common functional domains from the resting state brain may be used to improve the group statistics of task-fMRI data.

Correlations of diffusion tensor imaging values and symptom scores in patients with schizophrenia.

Abnormalities in white matter (WM) brain regions are attributed as a possible biomarker for schizophrenia (SZ). Diffusion tensor imaging (DTI) is used to capture WM tracts. Psychometric tests that evaluate the severity of symptoms of SZ are clinically used in the diagnosis process. In this study we investigate the correlates of scalar DTI measures, such as fractional anisotropy, mean diffusivity, axial diffusivity, and radial diffusivity with behavioral test scores. The correlations were found by different schemes: mean correlation with WM atlas regions and multiple regression of DTI values with test scores. The corpus callosum, superior longitudinal fasciculus right and inferior longitudinal fasciculus left were found to be having high correlations with test scores.

A method to classify schizophrenia using inter-task spatial correlations of functional brain images.

The clinical heterogeneity of schizophrenia (scz) and the overlap of self reported and observed symptoms with other mental disorders makes its diagnosis a difficult task. At present no laboratory-based or image-based diagnostic tool for scz exists and such tools are desired to support existing methods for more precise diagnosis. Functional magnetic resonance imaging (fMRI) is currently employed to identify and correlate cognitive processes related to scz and its symptoms. Fusion of multiple fMRI tasks that probe different cognitive processes may help to better understand hidden networks of this complex disorder. In this paper we utilize three different fMRI tasks and introduce an approach to classify subjects based on inter-task spatial correlations of brain activation. The technique was applied to groups of patients and controls and its validity was checked with the leave-one-out method. We show that the classification rate increases when information from multiple tasks are combined.