Generative discriminative models for multivariate inference and statistical mapping in medical imaging

This paper presents a general framework for obtaining interpretable multivariate discriminative models that allow efficient statistical inference for neuroimage analysis. The framework, termed generative discriminative machine (GDM), augments discriminative models with a generative regularization term. We demonstrate that the proposed formulation can be optimized in closed form and in dual space, allowing efficient computation for high dimensional neuroimaging datasets. Furthermore, we provide an analytic estimation of the null distribution of the model parameters, which enables efficient statistical inference and p-value computation without the need for permutation testing. We compared the proposed method with both purely generative and discriminative learning methods in two large structural magnetic resonance imaging (sMRI) datasets of Alzheimer's disease (AD) (n=415) and Schizophrenia (n=853). Using the AD dataset, we demonstrated the ability of GDM to robustly handle confounding variations. Using Schizophrenia dataset, we demonstrated the ability of GDM to handle multi-site studies. Taken together, the results underline the potential of the proposed approach for neuroimaging analyses.Comment: To appear in MICCAI 2018 proceeding

Advancing Statistical Inference For Population Studies In Neuroimaging Using Machine Learning

Modern neuroimaging techniques allow us to investigate the brain in vivo and in high resolution, providing us with high dimensional information regarding the structure and the function of the brain in health and disease. Statistical analysis techniques transform this rich imaging information into accessible and interpretable knowledge that can be used for investigative as well as diagnostic and prognostic purposes. A prevalent area of research in neuroimaging is group comparison, i.e., the comparison of the imaging data of two groups (e.g. patients vs. healthy controls or people who respond to treatment vs. people who don\u27t) to identify discriminative imaging patterns that characterize different conditions. In recent years, the neuroimaging community has adopted techniques from mathematics, statistics, and machine learning to introduce novel methodologies targeting the improvement of our understanding of various neuropsychiatric and neurodegenerative disorders. However, existing statistical methods are limited by their reliance on ad-hoc assumptions regarding the homogeneity of disease effect, spatial properties of the underlying signal and the covariate structure of data, which imposes certain constraints about the sampling of datasets. 1. First, the overarching assumption behind most analytical tools, which are commonly used in neuroimaging studies, is that there is a single disease effect that differentiates the patients from controls. In reality, however, the disease effect may be heterogeneously expressed across the patient population. As a consequence, when searching for a single imaging pattern that characterizes the difference between healthy controls and patients, we may only get a partial or incomplete picture of the disease effect. 2. Second, and importantly, most analyses assume a uniform shape and size of disease effect. As a consequence, a common step in most neuroimaging analyses it to apply uniform smoothing of the data to aggregate regional information to each voxel to improve the signal to noise ratio. However, the shape and size of the disease patterns may not be uniformly represented across the brain. 3. Lastly, in practical scenarios, imaging datasets commonly include variations due to multiple covariates, which often have effects that overlap with the searched disease effects. To minimize the covariate effects, studies are carefully designed by appropriately matching the populations under observation. The difficulty of this task is further exacerbated by the advent of big data analyses that often entail the aggregation of large datasets collected across many clinical sites. The goal of this thesis is to address each of the aforementioned assumptions and limitations by introducing robust mathematical formulations, which are founded on multivariate machine learning techniques that integrate discriminative and generative approaches. Specifically, 1. First, we introduce an algorithm termed HYDRA which stands for heterogeneity through discriminative analysis. This method parses the heterogeneity in neuroimaging studies by simultaneously performing clustering and classification by use of piecewise linear decision boundaries. 2. Second, we propose to perform regionally linear multivariate discriminative statistical mapping (MIDAS) toward finding the optimal level of variable smoothing across the brain anatomy and tease out group differences in neuroimaging datasets. This method makes use of overlapping regional discriminative filters to approximate a matched filter that best delineates the underlying disease effect. 3. Lastly, we develop a method termed generative discriminative machines (GDM) toward reducing the effect of confounds in biased samples. The proposed method solves for a discriminative model that can also optimally generate the data when taking into account the covariate structure. We extensively validated the performance of the developed frameworks in the presence of diverse types of simulated scenarios. Furthermore, we applied our methods on a large number of clinical datasets that included structural and functional neuroimaging data as well as genetic data. Specifically, HYDRA was used for identifying distinct subtypes of Alzheimer\u27s Disease. MIDAS was applied for identifying the optimally discriminative patterns that differentiated between truth-telling and lying functional tasks. GDM was applied on a multi-site prediction setting with severely confounded samples. Our promising results demonstrate the potential of our methods to advance neuroimaging analysis beyond the set of assumptions that limit its capacity and improve statistical power

Grey-matter texture abnormalities and reduced hippocampal volume are distinguishing features of schizophrenia

Neurodevelopmental processes are widely believed to underlie schizophrenia. Analysis of brain texture from conventional magnetic resonance imaging (MRI) can detect disturbance in brain cytoarchitecture. We tested the hypothesis that patients with schizophrenia manifest quantitative differences in brain texture that, alongside discrete volumetric changes, may serve as an endophenotypic biomarker. Texture analysis (TA) of grey matter distribution and voxel-based morphometry (VBM) of regional brain volumes were applied to MRI scans of 27 patients with schizophrenia and 24 controls. Texture parameters (uniformity and entropy) were also used as covariates in VBM analyses to test for correspondence with regional brain volume. Linear discriminant analysis tested if texture and volumetric data predicted diagnostic group membership (schizophrenia or control). We found that uniformity and entropy of grey matter differed significantly between individuals with schizophrenia and controls at the fine spatial scale (filter width below 2 mm). Within the schizophrenia group, these texture parameters correlated with volumes of the left hippocampus, right amygdala and cerebellum. The best predictor of diagnostic group membership was the combination of fine texture heterogeneity and left hippocampal size. This study highlights the presence of distributed grey-matter abnormalities in schizophrenia, and their relation to focal structural abnormality of the hippocampus. The conjunction of these features has potential as a neuroimaging endophenotype of schizophrenia

Structural and electrophysiological determinants of atrial cardiomyopathy identify remodeling discrepancies between paroxysmal and persistent atrial fibrillation

Background: Progressive atrial fibrotic remodeling has been reported to be associated with atrial cardiomyopathy (ACM) and the transition from paroxysmal to persistent atrial fibrillation (AF). We sought to identify the anatomical/structural and electrophysiological factors involved in atrial remodeling that promote AF persistency. Methods: Consecutive patients with paroxysmal (n = 134) or persistent (n = 136) AF who presented for their first AF ablation procedure were included. Patients underwent left atrial (LA) high-definition mapping (1,835 ± 421 sites/map) during sinus rhythm (SR) and were randomized to training and validation sets for model development and evaluation. A total of 62 parameters from both electro-anatomical mapping and non-invasive baseline data were extracted encompassing four main categories: (1) LA size, (2) extent of low-voltage-substrate (LVS), (3) LA voltages and (4) bi-atrial conduction time as identified by the duration of amplified P-wave (APWD) in a digital 12-lead-ECG. Least absolute shrinkage and selection operator (LASSO) and logistic regression were performed to identify the factors that are most relevant to AF persistency in each category alone and all categories combined. The performance of the developed models for diagnosis of AF persistency was validated regarding discrimination, calibration and clinical usefulness. In addition, HATCH score and C2HEST score were also evaluated for their performance in identification of AF persistency. Results: In training and validation sets, APWD (threshold 151 ms), LA volume (LAV, threshold 94 mL), bipolar LVS area < 1.0 mV (threshold 4.55 cm$^2$) and LA global mean voltage (GMV, threshold 1.66 mV) were identified as best determinants for AF persistency in the respective category. Moreover, APWD (AUC 0.851 and 0.801) and LA volume (AUC 0.788 and 0.741) achieved better discrimination between AF types than LVS extent (AUC 0.783 and 0.682) and GMV (AUC 0.751 and 0.707). The integrated model (combining APWD and LAV) yielded the best discrimination performance between AF types (AUC 0.876 in training set and 0.830 in validation set). In contrast, HATCH score and C2HEST score only achieved AUC < 0.60 in identifying individuals with persistent AF in current study. Conclusion: Among 62 electro-anatomical parameters, we identified APWD, LA volume, LVS extent, and mean LA voltage as the four determinant electrophysiological and structural factors that are most relevant for AF persistency. Notably, the combination of APWD with LA volume enabled discrimination between paroxysmal and persistent AF with high accuracy, emphasizing their importance as underlying substrate of persistent AF

Visual Feature Attribution using Wasserstein GANs

Attributing the pixels of an input image to a certain category is an important and well-studied problem in computer vision, with applications ranging from weakly supervised localisation to understanding hidden effects in the data. In recent years, approaches based on interpreting a previously trained neural network classifier have become the de facto state-of-the-art and are commonly used on medical as well as natural image datasets. In this paper, we discuss a limitation of these approaches which may lead to only a subset of the category specific features being detected. To address this problem we develop a novel feature attribution technique based on Wasserstein Generative Adversarial Networks (WGAN), which does not suffer from this limitation. We show that our proposed method performs substantially better than the state-of-the-art for visual attribution on a synthetic dataset and on real 3D neuroimaging data from patients with mild cognitive impairment (MCI) and Alzheimer's disease (AD). For AD patients the method produces compellingly realistic disease effect maps which are very close to the observed effects.Comment: Accepted to CVPR 201

Interaction between Neuroanatomical and Psychological Changes after Mindfulness-Based Training

Several cross-sectional studies have documented neuroanatomical changes in individuals with a long history of meditation, while a few evidences are available about the interaction between neuroanatomical and psychological changes even during brief exposure to meditation. Here we analyzed several morphometric indexes at both cortical and subcortical brain level, as well as multiple psychological dimensions, before and after a brief -8 weeks- Mindfulness Based Stress Reduction (MBSR) training program, in a group of 23 meditation naïve-subjects compared to age-gender matched subjects. We found a significant cortical thickness increase in the right insula and the somatosensory cortex of MBSR trainees, coupled with a significant reduction of several psychological indices related to worry, state anxiety, depression and alexithymia. Most importantly, an interesting correlation between the increase in right insula thickness and the decrease in alexithymia levels during the MBSR training were observed. Moreover, a multivariate pattern classification approach allowed to identify a cluster of regions more responsive to MBSR training across subjects. Taken together, these findings documented the significant impact of a brief MBSR training on brain structures, as well as stressing the idea of MBSR as a valuable tool for alexithymia modulation, also originally providing a plausible neurobiological evidence of a major role of right insula into mediating the observed psychological changes

Learning Interpretable Features of Graphs and Time Series Data

Graphs and time series are two of the most ubiquitous representations of data of modern time. Representation learning of real-world graphs and time-series data is a key component for the downstream supervised and unsupervised machine learning tasks such as classification, clustering, and visualization. Because of the inherent high dimensionality, representation learning, i.e., low dimensional vector-based embedding of graphs and time-series data is very challenging. Learning interpretable features incorporates transparency of the feature roles, and facilitates downstream analytics tasks in addition to maximizing the performance of the downstream machine learning models. In this thesis, we leveraged tensor (multidimensional array) decomposition for generating interpretable and low dimensional feature space of graphs and time-series data found from three domains: social networks, neuroscience, and heliophysics. We present the theoretical models and empirical results on node embedding of social networks, biomarker embedding on fMRI-based brain networks, and prediction and visualization of multivariate time-series-based flaring and non-flaring solar events

FUNCTIONAL NETWORK CONNECTIVITY IN HUMAN BRAIN AND ITS APPLICATIONS IN AUTOMATIC DIAGNOSIS OF BRAIN DISORDERS

The human brain is one of the most complex systems known to the mankind. Over the past 3500 years, mankind has constantly investigated this remarkable system in order to understand its structure and function. Emerging of neuroimaging techniques such as functional magnetic resonance imaging (fMRI) have opened a non-invasive in-vivo window into brain function. Moreover, fMRI has made it possible to study brain disorders such as schizophrenia from a different angle unknown to researchers before. Human brain function can be divided into two categories: functional segregation and integration. It is well-understood that each region in the brain is specialized in certain cognitive or motor tasks. The information processed in these specialized regions in different temporal and spatial scales must be integrated in order to form a unified cognition or behavior. One way to assess functional integration is by measuring functional connectivity (FC) among specialized regions in the brain. Recently, there is growing interest in studying the FC among brain functional networks. This type of connectivity, which can be considered as a higher level of FC, is termed functional network connectivity (FNC) and measures the statistical dependencies among brain functional networks. Each functional network may consist of multiple remote brain regions. Four studies related to FNC are presented in this work. First FNC is compared during the resting-state and auditory oddball task (AOD). Most previous FNC studies have been focused on either resting-state or task-based data but have not directly compared these two. Secondly we propose an automatic diagnosis framework based on resting-state FNC features for mental disorders such as schizophrenia. Then, we investigate the proper preprocessing for fMRI time-series in order to conduct FNC studies. Specifically the impact of autocorrelated time-series on FNC will be comprehensively assessed in theory, simulation and real fMRI data. At the end, the notion of autoconnectivity as a new perspective on human brain functionality will be proposed. It will be shown that autoconnectivity is cognitive-state and mental-state dependent and we discuss how this source of information, previously believed to originate from physical and physiological noise, can be used to discriminate schizophrenia patients with high accuracy

Classification of First-Episode Schizophrenia Patients and Healthy Subjects by Automated MRI Measures of Regional Brain Volume and Cortical Thickness

BACKGROUND: Although structural magnetic resonance imaging (MRI) studies have repeatedly demonstrated regional brain structural abnormalities in patients with schizophrenia, relatively few MRI-based studies have attempted to distinguish between patients with first-episode schizophrenia and healthy controls. METHOD: Three-dimensional MR images were acquired from 52 (29 males, 23 females) first-episode schizophrenia patients and 40 (22 males, 18 females) healthy subjects. Multiple brain measures (regional brain volume and cortical thickness) were calculated by a fully automated procedure and were used for group comparison and classification by linear discriminant function analysis. RESULTS: Schizophrenia patients showed gray matter volume reductions and cortical thinning in various brain regions predominantly in prefrontal and temporal cortices compared with controls. The classifiers obtained from 66 subjects of the first group successfully assigned 26 subjects of the second group with accuracy above 80%. CONCLUSION: Our results showed that combinations of automated brain measures successfully differentiated first-episode schizophrenia patients from healthy controls. Such neuroimaging approaches may provide objective biological information adjunct to clinical diagnosis of early schizophrenia