Hypothesis testing of matrix graph model with application to brain connectivity analysis.

The multiplicity problem has become increasingly important in genetic studies as the capacity for high-throughput genotyping has increased. The control of False Discovery Rate (FDR) (Benjamini and Hochberg. [1995] J. R. Stat. Soc. Ser. B 57:289-300) has been adopted to address the problems of false positive control and low power inherent in high-volume genome-wide linkage and association studies. In many genetic studies, there is often a natural stratification of the m hypotheses to be tested. Given the FDR framework and the presence of such stratification, we investigate the performance of a stratified false discovery control approach (i.e. control or estimate FDR separately for each stratum) and compare it to the aggregated method (i.e. consider all hypotheses in a single stratum). Under the fixed rejection region framework (i.e. reject all hypotheses with unadjusted p-values less than a pre-specified level and then estimate FDR), we demonstrate that the aggregated FDR is a weighted average of the stratum-specific FDRs. Under the fixed FDR framework (i.e. reject as many hypotheses as possible and meanwhile control FDR at a pre-specified level), we specify a condition necessary for the expected total number of true positives under the stratified FDR method to be equal to or greater than that obtained from the aggregated FDR method. Application to a recent Genome-Wide Association (GWA) study by Maraganore et al. ([2005] Am. J. Hum. Genet. 77:685-693) illustrates the potential advantages of control or estimation of FDR by stratum. Our analyses also show that controlling FDR at a low rate, e.g. 5% or 10%, may not be feasible for some GWA studies.

Frontiers Events is a rapidly growing calendar management system dedicated to the scheduling of academic events. This includes announcements and invitations, participant listings and search functionality, abstract handling and publication, related events and post-event exchanges. Whether an organizer or participant, make your event a Frontiers Event!

Multivariate statistics are often available as well as necessary in hypothesis tests. We study how to use such statistics to control not only false discovery rate (FDR) but also positive FDR (pFDR) with good power. We show that FDR can be controlled through nested regions of multivari- ate p-values of test statistics. If the distributions of the test statistics are known, then the regions can be constructed explicitly to achieve FDR con- trol with maximum power among procedures satisfying certain conditions. On the other hand, our focus is where the distributions are only partially known. Under certain conditions, a type of nested regions are proposed and shown to attain (p)FDR control with asymptotically maximum power as the pFDR control level approaches its attainable limit. The procedure based on the nested regions is compared with those based on other nested regions that are easier to construct as well as those based on more straightforward combinations of the test statistics. AMS 2000 subject classifications: Primary 62G10, 62H15; secondary 62G20.

Parcellation of the human cortex has important implications in neuroscience. Parcellation is often a crucial requirement before meaningful regional analysis can occur. The human cortex can be parcellated into distinct regions based on structural features, such as gyri and sulci. Brain network patterns in a given region with respect to its neighbors, known as connectional fingerprints, can be used to parcellate the cortex. Distinct imaging modalities might provide complementary information for brain parcellation. Here, we established functional connectivity with time series data from functional MRI (fMRI) combined with a correlation map of cortical thickness obtained from T1-weighted MRI. We aimed to extend the previous study, which parcellated the medial frontal cortex (MFC) using functional connectivity, and to test the value of additional information regarding cortical thickness. Two types of network information were used to parcellate the MFC into two sub-regions with spectral and Ward’s clustering approaches. The MFC region was defined using manual delineation based on in-house data (n = 12). Parcellation was applied to independent large-scale data obtained from the Human Connectome Project (HCP, n = 248). Agreement between parcellation using fMRI- and thickness-driven connectivity yielded dice coefficient overlaps of 0.74 (Ward’s clustering) and 0.54 (spectral clustering). We also explored whole brain connectivity using the MFC sub-regions as seed regions based on these two types of information. The results of whole brain connectivity analyses were also consistent for both types of information. We observed that an inter-regional correlation map derived from cortical thickness strongly reflected the underlying functional connectivity of MFC region.

Diffusion imaging and brain connectivity analyses can monitor white matter deterioration, revealing how neural pathways break down in aging and Alzheimer's disease (AD). Here we tested how AD disrupts the 'rich club' effect - a network property found in the normal brain - where high-degree nodes in the connectivity network are more heavily interconnected with each other than expected by chance. We analyzed 3-Tesla whole-brain diffusionweighted images (DWI) from 66 subjects (22 AD/44 normal elderly). We performed whole-brain tractography based on the orientation distribution functions. Connectivity matrices were compiled, representing the proportion of detected fibers interconnecting 68 cortical regions. As expected, AD patients had a lower nodal degree (average number of connections) in cortical regions implicated in the disease. Unexpectedly, the normalized rich club coefficient was higher in AD. AD disrupts cortical networks by removing connections; when these networks are thresholded, organizational properties are disrupted leading to additional new biomarkers of AD.

Major depressive disorder (MDD) is one of the most common mental health disorders that can affect sleep, mood, appetite, and behavior of people. Multimodal neuroimaging data, such as functional and structural magnetic resonance imaging (MRI) scans, have been widely used in computer-aided detection of MDD. However, previous studies usually treat these two modalities separately, without considering their potentially complementary information. Even though a few studies propose integrating these two modalities, they usually suffer from significant inter-modality data heterogeneity. In this paper, we propose an adaptive multimodal neuroimage integration (AMNI) framework for automated MDD detection based on functional and structural MRIs. The AMNI framework consists of four major components: (1) a graph convolutional network to learn feature representations of functional connectivity networks derived from functional MRIs, (2) a convolutional neural network to learn features of T1-weighted structural MRIs, (3) a feature adaptation module to alleviate inter-modality difference, and (4) a feature fusion module to integrate feature representations extracted from two modalities for classification. To the best of our knowledge, this is among the first attempts to adaptively integrate functional and structural MRIs for neuroimaging-based MDD analysis by explicitly alleviating inter-modality heterogeneity. Extensive evaluations are performed on 533 subjects with resting-state functional MRI and T1-weighted MRI, with results suggesting the efficacy of the proposed method.

Executive functions are supported by the dynamic coordination of neural activity over large-scale networks. The properties of large-scale brain coordination processes, however, remain unclear. Using tools combining MEG and intracranial EEG with brain connectivity analyses, we provide evidence that visuomotor behaviors, a hallmark of executive functions, are mediated by the interplay of multiple and spatially overlapping subnetworks. These subnetworks span visuomotor-related areas, the cortico-cortical and cortico-subcortical interactions of which evolve rapidly and reconfigure over timescales relevant for behavior. Visuomotor-related functional connectivity dynamics are characterized by sparse connections, nonstationarity, switching dynamics, and areal flexibility. We suggest that these properties represent key aspects of large-scale functional networks and cognitive architectures.

Differential network analysis investigates how the network of connected genes changes from one condition to another and has become a prevalent tool to provide a deeper and more comprehensive understanding of the molecular etiology of complex diseases. Based on the asymptotically normal estimation of large Gaussian graphical model (GGM) in the high-dimensional setting, we developed a computationally efficient test for differential network analysis through testing the equality of two precision matrices, which summarize the conditional dependence network structures of the genes. Additionally, we applied a multiple testing procedure to infer the differential network structure with false discovery rate (FDR) control. Through extensive simulation studies with different combinations of parameters including sample size, number of vertices, level of heterogeneity and graph structure, we demonstrated that our method performed much better than the current available methods in terms of accuracy and computational time. In real data analysis on lung adenocarcinoma, we revealed a differential network with 3503 nodes and 2550 edges, which consisted of 50 clusters with an FDR threshold at 0.05. Many of the top gene pairs in the differential network have been reported relevant to human cancers. Our method represents a powerful tool of network analysis for high-dimensional biological data.

The central nervous system, specifically the brain, is implicated in the development of exertional fatigue under a hot environment. Diverse neuroimaging techniques have been used to visualize the brain activity during or after exercise. Notably, the use of magnetic resonance imaging (MRI) has become prevalent due to its excellent spatial resolution and versatility. This review evaluates the significance and limitations of various brain MRI techniques in exercise studies-brain volumetric analysis, functional MRI, functional connectivity MRI, and arterial spin labeling. The review aims to provide a summary on the neural basis of exertional fatigue and proposes future directions for brain MRI studies. A systematic literature search was performed where a total of thirty-seven brain MRI studies associated with exercise, fatigue, or related physiological factors were reviewed. The findings suggest that with moderate dehydration, there is a decrease in total brain volume accompanied with expansion of ventricular volume. With exercise fatigue, there is increased activation of sensorimotor and cognitive brain areas, increased thalamo-insular activation and decreased interhemispheric connectivity in motor cortex. Under passive hyperthermia, there are regional changes in cerebral perfusion, a reduction in local connectivity in functional brain networks and an impairment to executive function. Current literature suggests that the brain structure and function are influenced by exercise, fatigue, and related physiological perturbations. However, there is still a dearth of knowledge and it is hoped that through understanding of MRI advantages and limitations, future studies will shed light on the central origin of exertional fatigue in the heat.

BackgroundFor gene expression or gene association studies with a large number of hypotheses the number of measurements per marker in a conventional single-stage design is often low due to limited resources. Two-stage designs have been proposed where in a first stage promising hypotheses are identified and further investigated in the second stage with larger sample sizes. For two types of two-stage designs proposed in the literature we derive multiple testing procedures controlling the False Discovery Rate (FDR) demonstrating FDR control by simulations: designs where a fixed number of top-ranked hypotheses are selected and designs where the selection in the interim analysis is based on an FDR threshold. In contrast to earlier approaches which use only the second-stage data in the hypothesis tests (pilot approach), the proposed testing procedures are based on the pooled data from both stages (integrated approach).ResultsFor both selection rules the multiple testing procedures control the FDR in the considered simulation scenarios. This holds for the case of independent observations across hypotheses as well as for certain correlation structures. Additionally, we show that in scenarios with small effect sizes the testing procedures based on the pooled data from both stages can give a considerable improvement in power compared to tests based on the second-stage data only.ConclusionThe proposed hypothesis tests provide a tool for FDR control for the considered two-stage designs. Comparing the integrated approaches for both selection rules with the corresponding pilot approaches showed an advantage of the integrated approach in many simulation scenarios.

Functional MRI resting state and connectivity studies of brain focus on neural fluctuations at low frequencies which share power with physiological fluctuations originating from lung and heart. Due to the lack of automated software to process physiological signals collected at high magnetic fields, a gap exists in the processing pathway between the acquisition of physiological data and its use in fMRI software for both physiological noise correction and functional analyses of brain activation and connectivity. To fill this gap, we developed an open source, physiological signal processing program, called PhysioNoise, in the python language. We tested its automated processing algorithms and dynamic signal visualization on resting monkey cardiac and respiratory waveforms. PhysioNoise consistently identifies physiological fluctuations for fMRI noise correction and also generates covariates for subsequent analyses of brain activation and connectivity.

False discovery rate (FDR) controlling procedures provide important statistical guarantees for replicability in signal identification based on multiple hypotheses testing. In many fields of study, FDR controling procedures are used in high-dimensional (HD) analyses to discover features that are truly associated with the outcome. In some recent applications, data on the same set of candidate features are independently collected in multiple different studies. For example, gene expression data are collected at different facilities and with different cohorts, to identify the genetic biomarkers of multiple types of cancers. These studies provide us with opportunities to identify signals by considering information from different sources (with potential heterogeneity) jointly. This paper is about how to provide FDR control guarantees for the tests of union null hypotheses of conditional independence. We present a knockoff-based variable selection method (Simultaneous knockoffs) to identify mutual signals from multiple independent datasets, providing exact FDR control guarantees under finite sample settings. This method can work with very general model settings and test statistics. We demonstrate the performance of this method with extensive numerical studies and two real-dataexamples.

Brain Connectivity: A Journal of Clinical Neurology, Neuroscience, and Neuroimaging Advancing the Field of Neurology.

Inferring brain connectivity network and quantifying the significance of interactions between brain regions are of paramount importance in neuroscience. Although there have recently emerged some tests for graph inference based on independent samples, there is no readily available solution to test the change of brain network for paired and correlated samples. In this article, we develop a paired test of matrix graphs to infer brain connectivity network when the groups of samples are correlated. The proposed test statistic is both bias corrected and variance corrected, and achieves a small estimation error rate. The subsequent multiple testing procedure built on this test statistic is guaranteed to asymptotically control the false discovery rate at the pre-specified level. Both the methodology and theory of the new test are considerably different from the two independent samples framework, owing to the strong correlations of measurements on the same subjects before and after the stimulus activity. We illustrate the efficacy of our proposal through simulations and an analysis of an Alzheimer's Disease Neuroimaging Initiative dataset.

Chronic pain is a common and severely disabling disease whose treatment is often unsatisfactory. Insights into the brain mechanisms of chronic pain promise to advance the understanding of the underlying pathophysiology and might help to develop disease markers and novel treatments. Here, we systematically exploited the potential of electroencephalography to determine abnormalities of brain function during the resting state in chronic pain. To this end, we performed state-of-the-art analyses of oscillatory brain activity, brain connectivity, and brain networks in 101 patients of either sex suffering from chronic pain. The results show that global and local measures of brain activity did not differ between chronic pain patients and a healthy control group. However, we observed significantly increased connectivity at theta (4-8 Hz) and gamma (>60 Hz) frequencies in frontal brain areas as well as global network reorganization at gamma frequencies in chronic pain patients. Furthermore, a machine learning algorithm could differentiate between patients and healthy controls with an above-chance accuracy of 57%, mostly based on frontal connectivity. These results suggest that increased theta and gamma synchrony in frontal brain areas are involved in the pathophysiology of chronic pain. Although substantial challenges concerning the reproducibility of the findings and the accuracy, specificity, and validity of potential electroencephalography-based disease markers remain to be overcome, our study indicates that abnormal frontal synchrony at theta and gamma frequencies might be promising targets for noninvasive brain stimulation and/or neurofeedback approaches.