Cognition and behavior depend on synchronized intrinsic brain activity which is organized into functional networks across the brain. Research has investigated how anatomical connectivity both shapes and is shaped by these networks, but not how anatomical connectivity interacts with intra-areal molecular properties to drive functional connectivity. Here, we present a novel linear model to explain functional connectivity in the mouse brain by integrating systematically obtained measurements of axonal connectivity, gene expression, and resting state functional connectivity MRI. The model suggests that functional connectivity arises from synergies between anatomical links and inter-areal similarities in gene expression. By estimating these interactions, we identify anatomical modules in which correlated gene expression and anatomical connectivity cooperatively, versus distinctly, support functional connectivity. Along with providing evidence that not all genes equally contribute to functional connectivity, this research establishes new insights regarding the biological underpinnings of coordinated brain activity measured by BOLD fMRI.
IntroductionThe brain is organized into a network of synchronized activity that has a complex and reproducible topological structure (1, 2). Resting state functional connectivity (FC) MRI, a technique which measures inter-areal correlations in spontaneous brain activity, has been particularly useful for studying functional network organization in both health and disease. Local and global features of this functional network are carefully calibrated to support healthy cognition (3) and network dysfunction is seen in numerous neurodevelopmental (4, 5) and neurodegenerative diseases (6, 7). Therefore, identifying the substrates that shape functional network organization is critical in linking molecular (e.g. gene transcription) and behavioral (e.g. psychometric) markers of disease to brain function.Despite an abundance of prior work examining the correspondence of large-scale functional and anatomical connectivity, the precise substrates that shape functional network organization remain unknown. Modeling approaches to predict FC networks based on macro-or meso-scale anatomical connectivity networks commonly simulate mass neuronal activity by optimizing parameters that describe local population dynamics as well as the contribution of inter-areal connectivity (8-10). These approaches allow for detailed theoretical exploration regarding the relative contributions of local dynamics vs global coupling, but are limited by a lack of empirical data regarding true areal differences in function. Furthermore, analytic measures of anatomical communication appear to predict FC at comparable values (11, 12), suggesting an upper limit to the predictive validity of models based on anatomical connectivity alone.The integration of diverse data from different scales of investigation in such models may enhance our understanding of how functional networks are shaped. Although the idea is intuitive to m...