Background: Volumetric alterations of subcortical structures as predictors of antipsychotic treatment response have been previously corroborated, but less is known about whether their morphological covariance relates to treatment outcome and is driven by gene expression and epigenetic modifications.Methods: Subcortical volumetric covariance was analyzed by using baseline T1-weighted magnetic resonance imaging (MRI) in 38 healthy controls and 38 drug-naïve first-episode schizophrenia patients. Patients were treated with 8-week risperidone monotherapy and divided into responder and non-responder groups according to the Remission in Schizophrenia Working Group (RSWG). We utilized partial least squares (PLS) regression to examine the spatial associations between gene expression of subcortical structures from a publicly available transcriptomic dataset and between-group variances of structural covariance. The peripheral DNA methylation (DNAm) status of a gene of interest (GOI), overlapping between genes detected in the PLS and 108 schizophrenia candidate gene loci previously reported, was examined in parallel with MRI scanning.Results: In the psychotic symptom dimension, non-responders had a higher baseline structural covariance in the putamen–hippocampus–pallidum–accumbens pathway compared with responders. For disorganized symptoms, significant differences in baseline structural covariant connections were found in the putamen–hippocampus–pallidum–thalamus circuit between the two subgroups. The imaging variances related to psychotic symptom response were spatially related to the expression of genes enriched in neurobiological processes and dopaminergic pathways. The DNAm of GOI demonstrated significant associations with patients’ improvement of psychotic symptoms.Conclusion: Baseline subcortical structural covariance and peripheral DNAm may relate to antipsychotic treatment response. Phenotypic variations in subcortical connectome related to psychotic symptom response may be transcriptomically and epigenetically underlaid. This study defines a roadmap for future studies investigating multimodal imaging epigenetic biomarkers for treatment response in schizophrenia.

Antipsychotic treatment has been conceived to alter brain connectivity, but it is unclear how the changes of network phenotypes relate to the underlying transcriptomics. Given DNA methylation (DNAm) may alter transcriptional levels, we further integrated an imaging-transcriptomic-epigenetic analysis to explore multi-omics treatment response biomarkers. Forty-two treatment-naive first-episode schizophrenia patients were scanned by TI weighted (T1W) imaging and DTI before and after 8-week risperidone monotherapy, and their peripheral blood genomic DNAm values were examined in parallel with MRI scanning. Morphometric similarity network (MSN) quantified with DTI and T1W data were used as a marker of treatment-related alterations in interareal cortical connectivity. We utilized partial least squares (PLS) to examine spatial associations between treatment-related MSN variations and cortical transcriptomic data obtained from the Allen Human Brain Atlas. Longitudinal MSN alterations were related to treatment response on cognitive function and general psychopathology symptoms, while DNAm values of 59 PLS1 genes were on negative and positive symptoms. Virtual-histology transcriptomic analysis linked the MSN alterations with the neurobiological, cellular and metabolic pathways or processes, and assigned MSN-related genes to multiple cell types, specifying neurons and glial cells as contributing most to the transcriptomic associations of longitudinal changes in MSN. We firstly reveal how brain-wide transcriptional levels and cell classes capture molecularly validated cortical connectivity alterations after antipsychotic treatment. Our findings represent a vital step towards the exploration of treatment response biomarkers on the basis of multiple omics rather than a single omics type as a strategy for advancing precise care.

Background: Mounting evidence from diffusion tensor imaging (DTI) and epigenetic studies, respectively, confirmed the abnormal alterations of brain white matter integrity and DNA methylation (DNAm) in schizophrenia. However, few studies have been carried out in the same sample to simultaneously explore the WM pathology relating to clinical behaviors, as well as the DNA methylation basis underlying the WM deficits.Methods: We performed DTI scans in 42 treatment-naïve first-episode schizophrenia patients and 38 healthy controls. Voxel-based method of fractional anisotropy (FA) derived from DTI was used to assess WM integrity. Participants' peripheral blood genomic DNAm status, quantified by using Infinium® Human Methylation 450K BeadChip, was examined in parallel with DTI scanning. Participants completed Digit Span test and Trail Making test, as well as Positive and Negative Syndrome Scale measurement. We acquired genes that are differentially expressed in the brain regions with abnormal FA values according to the Allen anatomically comprehensive atlas, obtained DNAm levels of the corresponding genes, and then performed Z-test to compare the differential epigenetic-imaging associations (DEIAs) between the two groups.Results: Significant decreases of FA values in the patient group were in the right middle temporal lobe WM, right cuneus WM, right anterior cingulate WM, and right inferior parietal lobe WM, while the significant increases were in the bilateral middle cingulate WM (Ps < 0.01, GRF correction). Abnormal FA values were correlated with patients' clinical symptoms and cognitive impairments. In the DEIAs, patients showed abnormal couple patterns between altered FA and DNAm components, for which the enriched biological processes and pathways could be largely grouped into three biological procedures: the neurocognition, immune, and nervous system.Conclusion: Schizophrenia may not cause widespread neuropathological changes, but subtle alterations affecting local cingulum WM, which may play a critical role in positive symptoms and cognitive impairments. This imaging-epigenetics study revealed for the first time that DNAm of genes enriched in neuronal, immunologic, and cognitive processes may serve as the basis in the effect of WM deficits on clinical behaviors in schizophrenia.

Background: Evidence from functional magnetic resonance imaging (fMRI) studies of schizophrenia suggests that interindividual variation in the stationary striatal functional circuit may be correlated with antipsychotic treatment response. However, little is known about the role of the dynamic striatum-related network in predicting patients’ clinical improvement. The spontaneous coactivation pattern (CAP) technique has recently been found to be important for elucidating the non-stationary nature of functional brain networks.Methods: Forty-two drug-naive first-episode schizophrenia patients underwent fMRI and T1W imaging before and after 8 weeks of risperidone monotherapy. The striatum was divided into 3 subregions, including the putamen, pallidum, and caudate. Spontaneous CAPs and CAP states were utilized to measure the dynamic characteristics of brain networks. We used DPARSF and Dynamic Brain Connectome software to analyze each subregion-related CAP and CAP state for each group and then compared the between-group differences in the neural network biomarkers. We used Pearson’s correlation analysis to determine the associations between the neuroimaging measurements with between-group differences and the improvement in patients’ psychopathological symptoms.Results: In the putamen-related CAPs, patients showed significantly increased intensity in the bilateral thalamus, bilateral supplementary motor areas, bilateral medial, and paracingulate gyrus, left paracentral lobule, left medial superior frontal gyrus, and left anterior cingulate gyrus compared with healthy controls. After treatment, thalamic signals in the putamen-related CAP 1 showed a significant increase, while the signals of the medial and paracingulate gyrus in the putamen-related CAP 3 revealed a significant decrease. The increase in thalamic signal intensity in the putamen-related CAP 1 was significantly and positively correlated with the percentage reduction in PANSS_P.Conclusion: This study is the first to combine striatal CAPs and fMRI to explore treatment response-related biomarkers in the early phase of schizophrenia. Our findings suggest that dynamic changes in CAP states in the putamen-thalamus circuit may be potential biomarkers for predicting patients’ variation in the short-term treatment response of positive symptoms.

Background: Moyamoya angiopathy (MMA) can be potentially missed in the initial magnetic resonance (MR) examination without MR angiography (MRA). The aim of this study was to develop an optimal deep learning model based on T2-weighted imaging (T2WI) for MMA detection. Methods: This retrospective multicenter study included MMA patients, control group patients with normal MRA and patients with cerebrovascular disease except MMA from seven hospitals (site 1 to site 7). Five models, namely shallow convolutional neural network (SCNN), LeNet-5 Convolutional Neural Network (LeNet), Visual Geometry Group Network (VGG), Residual Neural Network (ResNet) and Dense Convolutional Network (DenseNet), were used for training and validation. The model training and internal validation were performed with data from sites 1–4. Data from sites 5–7 were used for independent external validation, and the optimal model was selected according to the results of accuracy. Chi-squared test was used to further verify the influence of different MR manufacturers, field strength, age at the MR examination and MRA score on the optimal model. Results: A total of 1,038 MMA patients, 1,211 normal MRA and 271 patients with cerebrovascular disease except MMA were included. DenseNet showed the highest accuracy (0.859, 95% CI: 0.833, 0.884) in the independent external validation, which was not significantly different from that of VGG (0.834, 95% CI: 0.807, 0.861) and ResNet (0.855, 95% CI: 0.829, 0.880) but was significantly higher than that of SCNN (0.631, 95% CI: 0.595, 0.665; P<0.001) and LeNet (0.563, 95% CI: 0.527, 0.599; P=0.001). The accuracy of 1.5 T data was higher than 3.0 T data (χ²=6.559, P=0.01). The accuracy of MMA with MRA who scored more than 5 was higher than that scoring ≤5 (≤5 vs. 6–10: χ²=10.734, P=0.001; ≤5 vs. ≥11: χ²=10.369, P=0.001). Conclusions: DenseNet based on T2WI can be used to screen MMA, outperforming SCNN, LeNet, VGG and ResNet. The MRA score of MMA affected the DenseNet accuracy.