Abstract
The microbiome-gut-brain axis (MGBA) plays a critical role in schizophrenia (SZ). However, the underlying mechanisms of the interactions among the gut microbiome, brain networks, and symptom severity in SZ patients remain largely unknown. Fecal samples, structural and functional magnetic resonance imaging (MRI) data, and Positive and Negative Syndrome Scale (PANSS) scores were collected from 38 SZ patients and 38 normal controls, respectively. The data of 16S rRNA gene sequencing were used to analyze the abundance of gut microbiome and the analysis of human brain networks was applied to compute the nodal properties of 90 brain regions. A total of 1,691,280 mediation models were constructed based on 261 gut bacterial, 810 nodal properties, and 4 PANSS scores in SZ patients. A strong correlation between the gut microbiome and brain networks (r = 0.89, false discovery rate (FDR) -corrected p < 0.05) was identified. Importantly, the PANSS scores were linearly correlated with both the gut microbiome (r = 0.5, FDR-corrected p < 0.05) and brain networks (r = 0.59, FDR-corrected p < 0.05). The abundance of genus Sellimonas significantly affected the PANSS negative scores of SZ patients via the betweenness centrality of white matter networks in the inferior frontal gyrus and amygdala. Moreover, 19 significant mediation models demonstrated that the nodal properties of 7 brain regions, predominately from the systems of visual, language, and control of action, showed significant mediating effects on the PANSS scores with the gut microbiome as mediators. Together, our findings indicated the tripartite relationships among the gut microbiome, brain networks, and PANSS scores and suggested their potential role in the neuropathology of SZ.