Many therapeutic interventions in psychiatry can be viewed as attempts to influence the brain's large-scale, dynamic network state transitions. Building on connectome-based graph analysis and control theory, Network Control Theory is emerging as a powerful tool to quantify network controllability-i.e., the influence of one brain region over others regarding dynamic network state transitions. If and how network controllability is related to mental health remains elusive. Here, from Diffusion Tensor Imaging data, we inferred structural connectivity and inferred calculated network controllability parameters to investigate their association with genetic and familial risk in patients diagnosed with major depressive disorder (MDD, n = 692) and healthy controls (n = 820). First, we establish that controllability measures differ between healthy controls and MDD patients while not varying with current symptom severity or remission status. Second, we show that controllability in MDD patients is associated with polygenic scores for MDD and psychiatric cross-disorder risk. Finally, we provide evidence that controllability varies with familial risk of MDD and bipolar disorder as well as with body mass index. In summary, we show that network controllability is related to genetic, individual, and familial risk in MDD patients. We discuss how these insights into individual variation of network controllability may inform mechanistic models of treatment response prediction and personalized intervention-design in mental health.

Multiplex families with a high prevalence of a psychiatric disorder are often examined to identify rare genetic variants with large effect sizes. In the present study, we analysed whether the risk for bipolar disorder (BD) in BD multiplex families is influenced by common genetic variants. Furthermore, we investigated whether this risk is conferred mainly by BD-specific risk variants or by variants also associated with the susceptibility to schizophrenia or major depression. In total, 395 individuals from 33 Andalusian BD multiplex families (166 BD, 78 major depressive disorder, 151 unaffected) as well as 438 subjects from an independent, BD case/control cohort (161 unrelated BD, 277 unrelated controls) were analysed. Polygenic risk scores (PRS) for BD, schizophrenia (SCZ), and major depression were calculated and compared between the cohorts. Both the familial BD cases and unaffected family members had higher PRS for all three psychiatric disorders than the independent controls, with BD and SCZ being significant after correction for multiple testing, suggesting a high baseline risk for several psychiatric disorders in the families. Moreover, familial BD cases showed significantly higher BD PRS than unaffected family members and unrelated BD cases. A plausible hypothesis is that, in multiplex families with a general increase in risk for psychiatric disease, BD development is attributable to a high burden of common variants that confer a specific risk for BD. The present analyses demonstrated that common genetic risk variants for psychiatric disorders are likely to contribute to the high incidence of affective psychiatric disorders in the multiplex families. However, the PRS explained only part of the observed phenotypic variance, and rare variants might have also contributed to disease development.

Common variation in the gene encoding the neuron-specific RNA splicing factor RNA Binding Fox-1 Homolog 1 (RBFOX1) has been identified as a risk factor for several psychiatric conditions, and rare genetic variants have been found causal for autism spectrum disorder (ASD). Here, we explored the genetic landscape of RBFOX1 more deeply, integrating evidence from existing and new human studies as well as studies in Rbfox1 knockout mice. Mining existing data from large-scale studies of human common genetic variants, we confirmed gene-based and genome-wide association of RBFOX1 with risk tolerance, major depressive disorder and schizophrenia. Data on six mental disorders revealed copy number losses and gains to be more frequent in ASD cases than in controls. Consistently, RBFOX1 expression appeared decreased in post-mortem frontal and temporal cortices of individuals with ASD and prefrontal cortex of individuals with schizophrenia. Brain-functional MRI studies demonstrated that carriers of a common RBFOX1 variant, rs6500744, displayed increased neural reactivity to emotional stimuli, reduced prefrontal processing during cognitive control, and enhanced fear expression after fear conditioning, going along with increased avoidance behaviour. Investigating Rbfox1 neuron-specific knockout mice allowed us to further specify the role of this gene in behaviour. The model was characterised by pronounced hyperactivity, stereotyped behaviour, impairments in fear acquisition and extinction, reduced social interest, and lack of aggression; it provides excellent construct and face validity as an animal model of ASD. In conclusion, convergent translational evidence shows that common variants in RBFOX1 are associated with a broad spectrum of psychiatric traits and disorders, while rare genetic variation seems to expose to early-onset neurodevelopmental psychiatric disorders with and without developmental delay like ASD, in particular. Studying the pleiotropic nature of RBFOX1 can profoundly enhance our understanding of mental disorder vulnerability.

Major Depressive Disorder (MDD) often is associated with significant cognitive dysfunction. We conducted a meta-analysis of genome-wide interaction of MDD and cognitive function using data from four large European cohorts in a total of 3510 MDD cases and 6057 controls. In addition, we conducted analyses using polygenic risk scores (PRS) based on data from the Psychiatric Genomics Consortium (PGC) on the traits of MDD, Bipolar disorder (BD), Schizophrenia (SCZ), and mood instability (MIN). Functional exploration contained gene expression analyses and Ingenuity Pathway Analysis (IPA®). We identified a set of significantly interacting single nucleotide polymorphisms (SNPs) between MDD and the genome-wide association study (GWAS) of cognitive domains of executive function, processing speed, and global cognition. Several of these SNPs are located in genes expressed in brain, with important roles such as neuronal development (REST), oligodendrocyte maturation (TNFRSF21), and myelination (ARFGEF1). IPA® identified a set of core genes from our dataset that mapped to a wide range of canonical pathways and biological functions (MPO, FOXO1, PDE3A, TSLP, NLRP9, ADAMTS5, ROBO1, REST). Furthermore, IPA® identified upstream regulator molecules and causal networks impacting on the expression of dataset genes, providing a genetic basis for further clinical exploration (vitamin D receptor, beta-estradiol, tadalafil). PRS of MIN and meta-PRS of MDD, MIN and SCZ were significantly associated with all cognitive domains. Our results suggest several genes involved in physiological processes for the development and maintenance of cognition in MDD, as well as potential novel therapeutic agents that could be explored in patients with MDD associated cognitive dysfunction.

Predicting antidepressant response has been a clinical challenge for mood disorder. Although several genome-wide association studies have suggested a number of genetic variants to be associated with antidepressant response, the sample sizes are small and the results are difficult to replicate. Previous animal studies have shown that knockout of the serotonin receptor 7 gene (HTR7) resulted in an antidepressant-like phenotype, suggesting it was important to antidepressant action. In this report, in the first stage, we used a cost-effective pooled-sequencing strategy to sequence the entire HTR7 gene and its regulatory regions to investigate the association of common variants in HTR7 and clinical response to four selective serotonin reuptake inhibitors (SSRIs: citalopram, paroxetine, fluoxetine and sertraline) in a retrospective cohort mainly consisting of subjects with bipolar disorder (n = 359). We found 80 single-nucleotide polymorphisms (SNPs) with false discovery rate < 0.05 associated with response to paroxetine. Among the significant SNPs, rs7905446 (T/G), which is located at the promoter region, also showed nominal significance (P < 0.05) in fluoxetine group. GG/TG genotypes for rs7905446 and female gender were associated with better response to two SSRIs (paroxetine and fluoxetine). In the second stage, we replicated this association in two independent prospective samples of SSRI-treated patients with major depressive disorder: the MARS (n = 253, P = 0.0169) and GENDEP studies (n = 432, P = 0.008). The GG/TG genotypes were consistently associated with response in all three samples. Functional study of rs7905446 showed greater activity of the G allele in regulating expression of HTR7. The G allele displayed higher luciferase activity in two neuronal-related cell lines, and estrogen treatment decreased the activity of only the G allele. Electrophoretic mobility shift assay suggested that the G allele interacted with CCAAT/enhancer-binding protein beta transcription factor (TF), while the T allele did not show any interaction with any TFs. Our results provided novel pharmacogenomic evidence to support the role of HTR7 in association with antidepressant response.