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: Schizophrenia (SCZ) is a severe and chronic psychiatric disorder with premature age-related physiological changes. However, numerous previous studies examined the epigenetic age acceleration in SCZ patients and yielded inconclusive results. In this study, we propose to explore the epigenetic age acceleration in drug-naive first-episode SCZ (FSCZ) patients and investigate whether epigenetic age acceleration is associated with antipsychotic treatment, psychotic symptoms, cognition, and subcortical volumes. Methods: We assessed the epigenetic age in 38 drug-naive FSCZ patients and 38 healthy controls by using three independent clocks, including Horvath, Hannum and Levine algorithms. The epigenetic age measurements in SCZ patients were repeated after receiving 8 weeks risperidone monotherapy. Results: Our findings showed significantly positive correlations between epigenetic ages assessed by three clocks and chronological age in both FSCZ patients and healthy controls. Compared with healthy controls, drug-naive FSCZ patients have a significant epigenetic age deceleration in Horvath clock (p = 0.01), but not in Hannum clock (p = 0.07) and Levine clock (p = 0.43). The epigenetic ages of Hannum clock (p = 0.002) and Levine clock (p = 0.01) were significantly accelerated in SCZ patients after 8-week risperidone treatment. However, no significant associations between epigenetic age acceleration and psychotic symptoms, cognitive function, as well as subcortical volumes were observed in FSCZ patients. Conclusion: These results demonstrate that distinct epigenetic clocks are sensitive to different aspects of aging process. Further investigations with comprehensive epigenetic clock analyses and large samples are required to confirm our findings.

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: To improve the treatment outcomes of patients with schizophrenia, research efforts have focused on identifying brain-based markers of treatment response. Personal characteristics regarding disease-related behaviors likely stem from inter-individual variability in the organization of brain functional systems. This study aimed to track dimension-specific changes in psychotic symptoms following risperidone treatment using individual-level functional connectivity (FC). Methods: A reliable cortical parcellation approach that accounts for individual heterogeneity in cortical functional anatomy was used to localize functional regions in a longitudinal cohort, consisting of 42 drug-naive first-episodes schizophrenia (FES) patients at baseline and after 8 weeks of risperidone treatment. FC was calculated in individually specified brain regions and used to predict the baseline severity and improvement of positive and negative symptoms in FES. Results: Distinct sets of individual-specific FC were separately associated with the positive and negative symptom burden at baseline, which could be used to track the corresponding symptom resolution in FES patients following risperidone treatment. Between-network connections of the fronto-parietal network (FPN) contributed the most to predicting the positive symptom domain. A combination of between-network connections of the default mode network, FPN, and within-network connections of the FPN contributed markedly to the prediction model of negative symptom. Conclusion: This novel study, which accounts for individual brain variation, take a step toward establishing individual-specific theranostic biomarkers 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.