An Improved Fusion Paired Group Lasso Structured Sparse Canonical Correlation Analysis Based on Brain Imaging Genetics to Identify Biomarkers of Alzheimer’s Disease

Brain imaging genetics can demonstrate the complicated relationship between genetic factors and the structure or function of the humankind brain. Therefore, it has become an important research topic and attracted more and more attention from scholars. The structured sparse canonical correlation analysis (SCCA) model has been widely used to identify the association between brain image data and genetic data in imaging genetics. To investigate the intricate genetic basis of cerebrum imaging phenotypes, a great deal of other standard SCCA methods combining different interested structed have now appeared. For example, some models use group lasso penalty, and some use the fused lasso or the graph/network guided fused lasso for feature selection. However, prior knowledge may not be completely available and the group lasso methods have limited capabilities in practical applications. The graph/network guided approaches can use sample correlation to define constraints, thereby overcoming this problem. Unfortunately, this also has certain limitations. The graph/network conducted methods are susceptible to the sign of the sample correlation of the data, which will affect the stability of the model. To improve the efficiency and stability of SCCA, a sparse canonical correlation analysis model with GraphNet regularization (FGLGNSCCA) is proposed in this manuscript. Based on the FGLSCCA model, the GraphNet regularization penalty is imposed in our study and an optimization algorithm is presented to optimize the model. The structural Magnetic Resonance Imaging (sMRI) and gene expression data are used in this study to find the genotype and characteristics of brain regions associated with Alzheimer’s disease (AD). Experiment results shown that the new FGLGNSCCA model proposed in this manuscript is superior or equivalent to traditional methods in both artificially synthesized neuroimaging genetics data or actual neuroimaging genetics data. It can select essential features more powerfully compared with other multivariate methods and identify significant canonical correlation coefficients as well as captures more significant typical weight patterns which demonstrated its excellent ability in finding biologically important imaging genetic relations.

A multivariate to multivariate approach for voxel-wise genome-wide association analysis

The joint analysis of imaging-genetics data facilitates the systematic investigation of genetic effects on brain structures and functions with spatial specificity. We focus on voxel-wise genome-wide association analysis, which may involve trillions of single nucleotide polymorphism (SNP)-voxel pairs. We attempt to identify underlying organized association patterns of SNP-voxel pairs and understand the polygenic and pleiotropic networks on brain imaging traits. We propose a bi-clique graph structure (i.e., a set of SNPs highly correlated with a cluster of voxels) for the systematic association pattern. Next, we develop computational strategies to detect latent SNP-voxel bi-cliques and inference model for statistical testing. We further provide theoretical results to guarantee the accuracy of our computational algorithms and statistical inference. We validate our method by extensive simulation studies and then apply it to the whole genome genetic and voxel-level white matter integrity data collected from 1052 participants of the human connectome project (HCP). The results demonstrate multiple genetic loci influencing white matter integrity measures on splenium and genu of the corpus callosum.

A new Mendelian Randomization method to estimate causal effects of multivariable brain imaging exposures

The advent of simultaneously collected imaging-genetics data in large study cohorts provides an unprecedented opportunity to assess the causal effect of brain imaging traits on externally measured experimental results (e.g., cognitive tests) by treating genetic variants as instrumental variables. However, classic Mendelian Randomization methods are limited when handling high-throughput imaging traits as exposures to identify causal effects. We propose a new Mendelian Randomization framework to jointly select instrumental variables and imaging exposures, and then estimate the causal effect of multivariable imaging data on the outcome. We validate the proposed method with extensive data analyses and compare it with existing methods. We further apply our method to evaluate the causal effect of white matter microstructure integrity on cognitive function. The findings suggest that our method achieved better performance regarding sensitivity, bias, and false discovery rate compared to individually assessing the causal effect of a single exposure and jointly assessing the causal effect of multiple exposures without dimension reduction. Our application results indicated that WM measures across different tracts have a joint causal effect that significantly impacts the cognitive function among the participants from the UK Biobank.

Separating Clinical and Subclinical Depression by Big Data Informed Structural Vulnerability Index and Its impact on Cognition: ENIGMA Dot Product

Big Data neuroimaging collaborations including Enhancing Neuro Imaging Genetics through Meta-Analysis (ENIGMA) integrated worldwide data to identify regional brain deficits in major depressive disorder (MDD). We evaluated the sensitivity of translating ENIGMA-defined MDD deficit patterns to the individual level. We treated ENIGMA MDD deficit patterns as a vector to gauge the similarity between individual and MDD patterns by calculating ENIGMA dot product (EDP). We analyzed the sensitivity and specificity of EDP in separating subjects with (1) subclinical depressive symptoms without a diagnosis of MDD, (2) single episode MDD, (3) recurrent MDD, and (4) controls free of neuropsychiatric disorders. We compared EDP to the Quantile Regression Index (QRI; a linear alternative to the brain age metric) and the global gray matter thickness and subcortical volumes and fractional anisotropy (FA) of water diffusion. We performed this analysis in a large epidemiological sample of UK Biobank (UKBB) participants (N=17,053/19,265 M/F). Group-average increases in depressive symptoms from controls to recurrent MDD was mirrored by EDP (r2=0.85), followed by FA (r2=0.81) and QRI (r2=0.56). Subjects with MDD showed worse performance on cognitive tests than controls with deficits observed for 3 out of 9 cognitive tests administered by the UKBB. We calculated correlations of EDP and other brain indices with measures of cognitive performance in controls. The correlation pattern between EDP and cognition in controls was similar (r2=0.75) to the pattern of cognitive differences in MDD. This suggests that the elevation in EDP, even in controls, is associated with cognitive performance - specifically in the MDD-affected domains. That specificity was missing for QRI, FA or other brain imaging indices. In summary, translating anatomically informed meta-analytic indices of similarity using a linear vector approach led to better sensitivity to depressive symptoms and cognitive patterns than whole-brain imaging measurements or an index of accelerated aging.