A pilot multivariate parallel ICA study to investigate differential linkage between neural networks and genetic profiles in schizophrenia

Understanding genetic influences on both healthy and disordered brain function is a major focus in psychiatric neuroimaging. We utilized task-related imaging findings from an fMRI auditory oddball task known to be robustly associated with abnormal activation in schizophrenia, to investigate genomic factors derived from multiple single nucleotide polymorphisms (SNPs) from genes previously shown to be associated with schizophrenia. Our major aim was to investigate the relationship of these genomic factors to normal/abnormal brain functionality between controls and schizophrenia patients. We studied a Caucasian-only sample of 35 healthy controls and 31 schizophrenia patients. All subjects performed an auditory oddball task, which consists of detecting an infrequent sound within a series of frequent sounds. Each subject was characterized on 24 different SNP markers spanning multiple risk genes previously associated with schizophrenia. We used a recently developed technique named parallel independent component analysis (para-ICA) to analyze this multimodal data set (Liu et al., 2008). The method aims to identify simultaneously independent components of each modality (functional imaging, genetics) and the relationships between them. We detected three fMRI components significantly correlated with two distinct gene components. The fMRI components, along with their significant genetic profile (dominant SNP) correlations were as follows: (1) Inferior frontal-anterior/posterior cingulate-thalamus-caudate with SNPs from Brain derived neurotropic factor (BDNF) and dopamine transporter (DAT) [r=-0.51; p<0.0001], (2) superior/middle temporal gyrus-cingulate-premotor with SLC6A4_PR and SLC6A4_PR_AG (serotonin transporter promoter; 5HTTLPR) [r=0.27; p=0.03], and (3) default mode-fronto-temporal gyrus with Brain derived neurotropic factor and dopamine transporter (BDNF, DAT) [r=-0.25; p=0.04]. Functional components comprised task-relevant regions (including PFC, ACC, STG and MTG) frequently identified as abnormal in schizophrenia. Further, gene-fMRI combinations 1 (Z=1.75; p=0.03), 2 (Z=1.84; p=0.03) and 3 (Z=1.67; p=0.04) listed above showed significant differences between controls and patients, based on their correlated loading coefficients. We demonstrate a framework to identify interactions between "clusters" of brain function and of genetic information. Our results reveal the effect/influence of specific interactions, (perhaps epistastatic in nature), between schizophrenia risk genes on imaging endophenotypes representing attention/working memory and goal directed related brain function, thus establishing a useful methodology to probe multivariate genotype-phenotype relationships.