Dongdong Lin

Correspondence between fMRI and SNP data by group sparse canonical correlation analysis

Both genetic variants and brain region abnormalities are recognized as important factors for complex diseases (e.g., schizophrenia). In this paper, we investigated the correspondence between single nucleotide polymorphism (SNP) and brain activity measured by functional magnetic resonance imaging (fMRI) to understand how genetic variation influences the brain activity. A group sparse canonical correlation analysis method (group sparse CCA) was developed to explore the correlation between these two datasets which are high dimensional-the number of SNPs/voxels is far greater than the number of samples. Different from the existing sparse CCA methods (sCCA), our approach can exploit structural information in the correlation analysis by introducing group constraints. A simulation study demonstrates that it outperforms the existing sCCA. We applied this method to the real data analysis and identified two pairs of significant canonical variates with average correlations of 0.4527 and 0.4292 respectively, which were used to identify genes and voxels associated with schizophrenia. The selected genes are mostly from 5 schizophrenia (SZ)-related signalling pathways. The brain mappings of the selected voxles also indicate the abnormal brain regions susceptible to schizophrenia. A gene and brain region of interest (ROI) correlation analysis was further performed to confirm the significant correlations between genes and ROIs.

Interaction among subsystems within default mode network diminished in schizophrenia patients: A dynamic connectivity approach

Default mode network (DMN) has been reported altered in schizophrenia (SZ) using static connectivity analysis. However, the studies on dynamic characteristics of DMN in SZ are still limited. In this work, we compare dynamic connectivity within DMN between 82 healthy controls (HC) and 82 SZ patients using resting-state fMRI. Firstly, dynamic DMN was computed using a sliding time window method for each subject. Then, the overall connectivity strengths were compared between two groups. Furthermore, we estimated functional connectivity states using K-means clustering, and then investigated group differences with respect to the connectivity strengths in states, the dwell time in each state, and the transition times between states. Finally, graph metrics of time-varying connectivity patterns and connectivity states were assessed. Results suggest that measured by the overall connectivity, HC showed stronger inter-subsystem interaction than patients. Compared to HC, patients spent more time in the states with nodes sparsely connected. For each state, SZ patients presented relatively weaker connectivity strengths mainly in inter-subsystem. Patients also exhibited lower values in averaged node strength, clustering coefficient, global efficiency, and local efficiency than HC. In summary, our findings indicate that SZ show impaired interaction among DMN subsystems, with a reduced central role for posterior cingulate cortex (PCC) and anterior medial prefrontal cortex (aMPFC) hubs as well as weaker interaction between dorsal medial prefrontal cortex (dMPFC) subsystem and medial temporal lobe (MTL) subsystem. For SZ, decreased integration of DMN may be associated with impaired ability in making self-other distinctions and coordinating present mental states with episodic decisions about future.

Group sparse canonical correlation analysis for genomic data integration

BackgroundThe emergence of high-throughput genomic datasets from different sources and platforms (e.g., gene expression, single nucleotide polymorphisms (SNP), and copy number variation (CNV)) has greatly enhanced our understandings of the interplay of these genomic factors as well as their influences on the complex diseases. It is challenging to explore the relationship between these different types of genomic data sets. In this paper, we focus on a multivariate statistical method, canonical correlation analysis (CCA) method for this problem. Conventional CCA method does not work effectively if the number of data samples is significantly less than that of biomarkers, which is a typical case for genomic data (e.g., SNPs). Sparse CCA (sCCA) methods were introduced to overcome such difficulty, mostly using penalizations with l-1 norm (CCA-l1) or the combination of l-1and l-2 norm (CCA-elastic net). However, they overlook the structural or group effect within genomic data in the analysis, which often exist and are important (e.g., SNPs spanning a gene interact and work together as a group).ResultsWe propose a new group sparse CCA method (CCA-sparse group) along with an effective numerical algorithm to study the mutual relationship between two different types of genomic data (i.e., SNP and gene expression). We then extend the model to a more general formulation that can include the existing sCCA models. We apply the model to feature/variable selection from two data sets and compare our group sparse CCA method with existing sCCA methods on both simulation and two real datasets (human gliomas data and NCI60 data). We use a graphical representation of the samples with a pair of canonical variates to demonstrate the discriminating characteristic of the selected features. Pathway analysis is further performed for biological interpretation of those features.ConclusionsThe CCA-sparse group method incorporates group effects of features into the correlation analysis while performs individual feature selection simultaneously. It outperforms the two sCCA methods (CCA-l1 and CCA-group) by identifying the correlated features with more true positives while controlling total discordance at a lower level on the simulated data, even if the group effect does not exist or there are irrelevant features grouped with true correlated features. Compared with our proposed CCA-group sparse models, CCA-l1 tends to select less true correlated features while CCA-group inclines to select more redundant features.

Sparse representation based biomarker selection for schizophrenia with integrated analysis of fMRI and SNPs.

Integrative analysis of multiple data types can take advantage of their complementary information and therefore may provide higher power to identify potential biomarkers that would be missed using individual data analysis. Due to different natures of diverse data modality, data integration is challenging. Here we address the data integration problem by developing a generalized sparse model (GSM) using weighting factors to integrate multi-modality data for biomarker selection. As an example, we applied the GSM model to a joint analysis of two types of schizophrenia data sets: 759,075 SNPs and 153,594 functional magnetic resonance imaging (fMRI) voxels in 208 subjects (92 cases/116 controls). To solve this small-sample-large-variable problem, we developed a novel sparse representation based variable selection (SRVS) algorithm, with the primary aim to identify biomarkers associated with schizophrenia. To validate the effectiveness of the selected variables, we performed multivariate classification followed by a ten-fold cross validation. We compared our proposed SRVS algorithm with an earlier sparse model based variable selection algorithm for integrated analysis. In addition, we compared with the traditional statistics method for uni-variant data analysis (Chi-squared test for SNP data and ANOVA for fMRI data). Results showed that our proposed SRVS method can identify novel biomarkers that show stronger capability in distinguishing schizophrenia patients from healthy controls. Moreover, better classification ratios were achieved using biomarkers from both types of data, suggesting the importance of integrative analysis.

Predicting individualized clinical measures by a generalized prediction framework and multimodal fusion of MRI data

ABSTRACT Neuroimaging techniques have greatly enhanced the understanding of neurodiversity (human brain variation across individuals) in both health and disease. The ultimate goal of using brain imaging biomarkers is to perform individualized predictions. Here we proposed a generalized framework that can predict explicit values of the targeted measures by taking advantage of joint information from multiple modalities. This framework also enables whole brain voxel‐wise searching by combining multivariate techniques such as ReliefF, clustering, correlation‐based feature selection and multiple regression models, which is more flexible and can achieve better prediction performance than alternative atlas‐based methods. For 50 healthy controls and 47 schizophrenia patients, three kinds of features derived from resting‐state fMRI (fALFF), sMRI (gray matter) and DTI (fractional anisotropy) were extracted and fed into a regression model, achieving high prediction for both cognitive scores (MCCB composite r = 0.7033, MCCB social cognition r = 0.7084) and symptomatic scores (positive and negative syndrome scale [PANSS] positive r = 0.7785, PANSS negative r = 0.7804). Moreover, the brain areas likely responsible for cognitive deficits of schizophrenia, including middle temporal gyrus, dorsolateral prefrontal cortex, striatum, cuneus and cerebellum, were located with different weights, as well as regions predicting PANSS symptoms, including thalamus, striatum and inferior parietal lobule, pinpointing the potential neuromarkers. Finally, compared to a single modality, multimodal combination achieves higher prediction accuracy and enables individualized prediction on multiple clinical measures. There is more work to be done, but the current results highlight the potential utility of multimodal brain imaging biomarkers to eventually inform clinical decision‐making. HighlightsA generalized framework to predict explicit values of clinical measures using MRI.Locating imaging biomarkers/neuromarkers by whole brain voxel‐wise searching.A prediction model incorporating multi‐stage feature selection and multimodal info.May potentially help clinician to make early intervention for mental disorders.

Sparse models for correlative and integrative analysis of imaging and genetic data

The development of advanced medical imaging technologies and high-throughput genomic measurements has enhanced our ability to understand their interplay as well as their relationship with human behavior by integrating these two types of datasets. However, the high dimensionality and heterogeneity of these datasets presents a challenge to conventional statistical methods; there is a high demand for the development of both correlative and integrative analysis approaches. Here, we review our recent work on developing sparse representation based approaches to address this challenge. We show how sparse models are applied to the correlation and integration of imaging and genetic data for biomarker identification. We present examples on how these approaches are used for the detection of risk genes and classification of complex diseases such as schizophrenia. Finally, we discuss future directions on the integration of multiple imaging and genomic datasets including their interactions such as epistasis.

Integrating fMRI and SNP data for biomarker identification for schizophrenia with a sparse representation based variable selection method

BackgroundIn recent years, both single-nucleotide polymorphism (SNP) array and functional magnetic resonance imaging (fMRI) have been widely used for the study of schizophrenia (SCZ). In addition, a few studies have been reported integrating both SNPs data and fMRI data for comprehensive analysis.MethodsIn this study, a novel sparse representation based variable selection (SRVS) method has been proposed and tested on a simulation data set to demonstrate its multi-resolution properties. Then the SRVS method was applied to an integrative analysis of two different SCZ data sets, a Single-nucleotide polymorphism (SNP) data set and a functional resonance imaging (fMRI) data set, including 92 cases and 116 controls. Biomarkers for the disease were identified and validated with a multivariate classification approach followed by a leave one out (LOO) cross-validation. Then we compared the results with that of a previously reported sparse representation based feature selection method.ResultsResults showed that biomarkers from our proposed SRVS method gave significantly higher classification accuracy in discriminating SCZ patients from healthy controls than that of the previous reported sparse representation method. Furthermore, using biomarkers from both data sets led to better classification accuracy than using single type of biomarkers, which suggests the advantage of integrative analysis of different types of data.ConclusionsThe proposed SRVS algorithm is effective in identifying significant biomarkers for complicated disease as SCZ. Integrating different types of data (e.g. SNP and fMRI data) may identify complementary biomarkers benefitting the diagnosis accuracy of the disease.

Dynamic functional connectivity impairments in early schizophrenia and clinical high-risk for psychosis

ABSTRACT Individuals at clinical high‐risk (CHR) for psychosis are characterized by attenuated psychotic symptoms. Only a minority of CHR individuals convert to full‐blown psychosis. Therefore, there is a strong interest in identifying neurobiological abnormalities underlying the psychosis risk syndrome. Dynamic functional connectivity (DFC) captures time‐varying connectivity over short time scales, and has the potential to reveal complex brain functional organization. Based on resting‐state functional magnetic resonance imaging (fMRI) data from 70 healthy controls (HCs), 53 CHR individuals, and 58 early illness schizophrenia (ESZ) patients, we applied a novel group information guided ICA (GIG‐ICA) to estimate inherent connectivity states from DFC, and then investigated group differences. We found that ESZ patients showed more aberrant connectivities and greater alterations than CHR individuals. Results also suggested that disease‐related connectivity states occurred in CHR and ESZ groups. Regarding the dominant state with the highest contribution to dynamic connectivity, ESZ patients exhibited greater impairments than CHR individuals primarily in the cerebellum, frontal cortex, thalamus and temporal cortex, while CHR and ESZ populations shared common aberrances mainly in the supplementary motor area, parahippocampal gyrus and postcentral cortex. CHR‐specific changes were also found in the connections between the superior frontal gyrus and calcarine cortex in the dominant state. Our findings suggest that CHR individuals generally show an intermediate functional connectivity pattern between HCs and SZ patients but also have unique connectivity alterations.

Identifying dynamic functional connectivity biomarkers using GIG‐ICA: Application to schizophrenia, schizoaffective disorder, and psychotic bipolar disorder

Functional magnetic resonance imaging (fMRI) studies have shown altered brain dynamic functional connectivity (DFC) in mental disorders. Here, we aim to explore DFC across a spectrum of symptomatically‐related disorders including bipolar disorder with psychosis (BPP), schizoaffective disorder (SAD), and schizophrenia (SZ). We introduce a group information guided independent component analysis procedure to estimate both group‐level and subject‐specific connectivity states from DFC. Using resting‐state fMRI data of 238 healthy controls (HCs), 140 BPP, 132 SAD, and 113 SZ patients, we identified measures differentiating groups from the whole‐brain DFC and traditional static functional connectivity (SFC), separately. Results show that DFC provided more informative measures than SFC. Diagnosis‐related connectivity states were evident using DFC analysis. For the dominant state consistent across groups, we found 22 instances of hypoconnectivity (with decreasing trends from HC to BPP to SAD to SZ) mainly involving post‐central, frontal, and cerebellar cortices as well as 34 examples of hyperconnectivity (with increasing trends HC through SZ) primarily involving thalamus and temporal cortices. Hypoconnectivities/hyperconnectivities also showed negative/positive correlations, respectively, with clinical symptom scores. Specifically, hypoconnectivities linking postcentral and frontal gyri were significantly negatively correlated with the PANSS positive/negative scores. For frontal connectivities, BPP resembled HC while SAD and SZ were more similar. Three connectivities involving the left cerebellar crus differentiated SZ from other groups and one connection linking frontal and fusiform cortices showed a SAD‐unique change. In summary, our method is promising for assessing DFC and may yield imaging biomarkers for quantifying the dimension of psychosis. Hum Brain Mapp 38:2683–2708, 2017.

Bio marker identification for diagnosis of schizophrenia with integrated analysis of fMRI and SNPs

It is important to identify significant biomarkers such as SNPs for medical diagnosis and treatment. However, the size of a biological sample is usually far less than the number of measurements, which makes the problem more challenging. To overcome this difficulty, we propose a sparse representation based variable selection (SRVS) approach. A simulated data set was first tested to demonstrate the advantages and properties of the proposed method. Then, we applied the algorithm to a joint analysis of 759075 SNPs and 153594 functional magnetic resonance imaging (fMRJ) voxels in 208 subjects (92 cases/116 controls) to identify significant biomarkers for schizophrenia (SZ). When compared with previous studies, our proposed method located 20 genes out of the top 45 SZ genes that are publicly reported We also detected some interesting functional brain regions from the fMRI study. In addition, a leave one out (LOO) cross-validation was performed and the results were compared with that of a previously reported method, which showed that our method gave significantly higher classification accuracy. In addition, the identification accuracy with integrative analysis is much better than that of using single type of data, suggesting that integrative analysis may lead to better diagnostic accuracy by combining complementary SNP and fMRI data.