Decomposition of brain diffusion imaging data uncovers latent schizophrenias with distinct patterns of white matter anisotropy

Fractional anisotropy (FA) analysis of diffusion tensor-images (DTI) has yielded inconsistent abnormalities in schizophrenia (SZ). Inconsistencies may arise from averaging heterogeneous groups of patients. Here we investigate whether SZ is a heterogeneous group of disorders distinguished by distinct patterns of FA reductions. We developed a Generalized Factorization Method (GFM) to identify biclusters (i.e., subsets of subjects associated with a subset of particular characteristics, such as low FA in specific regions). GFM appropriately assembles a collection of unsupervised techniques with Non-negative Matrix Factorization to generate biclusters, rather than averaging across all subjects and all their characteristics. DTI tract-based spatial statistics images, which output is the locally maximal FA projected onto the group white matter skeleton, were analyzed in 47 SZ and 36 healthy subjects, identifying 8 biclusters. The mean FA of the voxels of each bicluster was significantly different from those of other SZ subjects or 36 healthy controls. The eight biclusters were organized into four more general patterns of low FA in specific regions: 1) genu of corpus callosum (GCC), 2) fornix (FX) + external capsule (EC), 3) splenium of CC (SCC) + retrolenticular limb (RLIC) + posterior limb (PLIC) of the internal capsule, and 4) anterior limb of the internal capsule. These patterns were significantly associated with particular clinical features: Pattern 1 (GCC) with bizarre behavior, pattern 2 (FX + EC) with prominent delusions, and pattern 3 (SCC + RLIC + PLIC) with negative symptoms including disorganized speech. The uncovered patterns suggest that SZ is a heterogeneous group of disorders that can be distinguished by different patterns of FA reductions associated with distinct clinical features.