Neurodevelopmental trajectories, polygenic risk, and lipometabolism in vulnerability and resilience to schizophrenia
BMC Psychiatry volume 23, Article number: 153 (2023)
Schizophrenia (SZ) arises from a complex interplay involving genetic and molecular factors. Early intervention of SZ hinges upon understanding its vulnerability and resiliency factors in study of SZ and genetic high risk for SZ (GHR).
Herein, using integrative and multimodal strategies, we first performed a longitudinal study of neural function as measured by amplitude of low frequency function (ALFF) in 21 SZ, 26 GHR, and 39 healthy controls to characterize neurodevelopmental trajectories of SZ and GHR. Then, we examined the relationship between polygenic risk score for SZ (SZ-PRS), lipid metabolism, and ALFF in 78 SZ, and 75 GHR in cross-sectional design to understand its genetic and molecular substrates.
Across time, SZ and GHR diverge in ALFF alterations of the left medial orbital frontal cortex (MOF). At baseline, both SZ and GHR had increased left MOF ALFF compared to HC (P < 0.05). At follow-up, increased ALFF persisted in SZ, yet normalized in GHR. Further, membrane genes and lipid species for cell membranes predicted left MOF ALFF in SZ; whereas in GHR, fatty acids best predicted and were negatively correlated (r = -0.302, P < 0.05) with left MOF.
Our findings implicate divergence in ALFF alteration in left MOF between SZ and GHR with disease progression, reflecting vulnerability and resiliency to SZ. They also indicate different influences of membrane genes and lipid metabolism on left MOF ALFF in SZ and GHR, which have important implications for understanding mechanisms underlying vulnerability and resiliency in SZ and contribute to translational efforts for early intervention.
Schizophrenia (SZ) manifests from a complex interplay of genes and molecules that shape brain structure and function across the lifespan. It appears to be a highly heritable neurodevelopmental disorder [1, 2] with progression of neural alterations in early development leading to psychotic onset usually following puberty. Understanding the vulnerability and resilency factors in SZ is critical for early intervention of this debilitating illness. Studies of SZ and nonpsychotic individuals at genetic high risk for SZ (GHR) and longitudinal designs are especially important for elucidating vulnerability and resiliency to SZ and preventive interventions . GHR individuals, which are defined as those with first-degree relatives with SZ, have almost a tenfold increased risk of developing SZ . Similarities in neural alterations between SZ and GHR have been observed in early development [5,6,7]. However, prior studies indicate that most GHR do not develop SZ in their lifetime and early neurodevelopmental similarities normalize with age in GHR but not in SZ [7,8,9,10]. These studies implicate resiliency factors in GHR and suggest divergent neurodevelopmental trajectories between SZ and GHR.
Longitudinal neuroimaging studies is of increasing importance in characterizing neural trajectory. However, there are significant challenges in performing longitudinal studies in human subjects including participant attrition over time, feasibility of multimodal data collection at every timepoint, and sufficient funding to sustain such study. Some of these challenges could be mitigated by cross-sectional studies that consist of large sample sizes and multimodal data (e.g., genetic and metabolic) collection allowing integrative analysis. Further, previous studies implicate both genetic and metabolic factors in SZ etiology and pathophysiology [11, 12]. Along with neuroimaging, genomic and metabolomic data may provide a more comprehensive perspective on the development of SZ and related vulnerability and resiliency factors.
Genome-wide association studies (GWAS) offer a powerful approach to understanding the genetic basis for brain alterations. While each individual disease-associated single nucleotide polymorphism (SNP) carries only a subtle increase in SZ risk (with odds ratios in the range of 1.1 to 1.2), the cumulative risk of the various SNPs amount to a measure called the polygenic risk score (PRS), capturing the polygenic nature of complex disorders like SZ . Previous studies have examined the effect of PRS for SZ (SZ-PRSs) on functional brain alterations of SZ [14,15,16,17]. Studies have shown that SZ-PRS could predict mnemonic hippocampal activity  and is associated with left dorsolateral prefrontal cortex (PFC) inefficiency of SZ . A family study also suggests that SZ-PRS affects early neurodevelopment and indicates an increased risk of developing the disorder . Therefore, imaging genetics using PRS could lead to better understanding of genetic determinants of neurodevelopmental deficits in SZ .
Metabolomics is also a promising technology to understand brain alterations at the molecular level. Lipids are especially important regulators of brain function and are increasingly implicated in neuropsychiatric disease . One study suggested that changes in lipometabolism presented early in the development of SZ . Levels of lipids were correlated with the severity of negative symptoms in SZ [22, 23]. Moreover, lipid abnormalities have been shown in both postmortem brain tissue and peripheral blood of SZ [24, 25]. In addition, an animal study found that lipid washout led to a change in brain structure . Altogether, these findings demonstrate the importance of lipometabolism in brain alterations associated with SZ.
SZ is a multifaceted disorder arising from a complex interplay of genetic and molecular factors. Thus, integrative and multimodal strategies are needed to comprehensively understand SZ development and pathophysiology. Amplitude of low frequency fluctuations (ALFF), an efficient index of local spontaneous neuronal activity at rest , is used to examine neural function. ALFF exhibits moderate to substantial test–retest reliability  ensuring a high upper bound for its validity. Evidence has suggested ALFF alterations both in SZ and GHR [29, 30], reflecting their altered neuronal function. In this study, we first performed a longitudinal neuroimaging study of ALFF in SZ and GHR to characterize neurodevelopmental trajectories of SZ and GHR. Using divergent ALFF alterations in SZ and GHR across time, we then examined the influence of polygenic risk and lipometabolism on neural function in SZ and GHR in cross-sectional designs to comprehensively understand genetic and molecular substrates of altered neural function.
Materials and methods
The study was approved by the Medical Science Research Ethics Committee of the First Affiliated Hospital of China Medical University (approval reference number 25–1). All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. All participants provided written informed consent by themselves or by their parents/guardians if they were under 18 years old after a complete description of the study. SZ and GHR participants were recruited from the inpatient and outpatient services at Shenyang Mental Health Center and the Department of Psychiatry at First Affiliated Hospital of China Medical University. Healthy controls (HC) participants were recruited from the local community by advertisement.
All components of the study were conducted at a single site and included both longitudinal and cross-sectional study cohorts, aged 13–45 years. All participants were evaluated by 2 trained psychiatrists to determine the presence or absence of Axis I psychiatric diagnoses using the Structured Clinical Interview for Diagnostic and Statistical Manual of Mental Disorders-IV-Text Revision (DSM-IV) Axis I Disorders (SCID) in those 18 years old and older and the Schedule for Affective Disorders and Schizophrenia for School-Age Children-present and Lifetime Version (K-SADS-PL) in those younger than 18 years. SZ participants met DSM-IV diagnostic criteria for SZ and not any other Axis I disorder. GHR participants were first-degree relatives of individuals with SZ and did not meet criteria for any DSM-IV Axis I disorder. HC participants did not have current or lifetime Axis I disorder or history of psychotic, mood, or other Axis I disorders in first-degree relatives as determined by detailed family history. Participants were excluded if any of the following were present: (1) the existence of substance/alcohol abuse or dependence or concomitant major medical disorder, (2) any magnetic resonance imaging (MRI) contraindications, and (3) history of head trauma with loss of consciousness for ≥ 5 min or any neurological disorder. Symptom severity was measured using the Brief Psychiatric Rating Scale (BPRS).
Longitudinal ALFF alterations
Participants were recruited and participated between April 2011 and April 2019. A total of 86 participants were included in this portion of the study, including 21 SZ, 26 GHR, and 39 HC. At baseline, the mean age was 25.93 (7.65) years, aged 13–45 years. There were no significant differences in age among the SZ, GHR, and HC groups; however significant differences were observed in sex (χ2 = 12.82, p = 0.002). All participants underwent clinical and resting state functional MRI (R-fMRI) assessment at baseline and follow-up at least 10 months after initial scan (ranged 11–67 months). The mean duration between baseline and follow-up scans was 26.16 (14.28) months. Diagnoses were confirmed again at follow-up (using SCID or K-SADS-PL) with no diagnostic change for any participants.
The influence of polygenic risk and lipometabolism on ALFF alterations
A total of 298 participants, consisting of 78 SZ, 75 GHR, and 145 HC, were included in this portion of the study. Baseline measures for 51 participants from the longitudinal study were included for this portion. Blood samples were obtained within 24 h of scanning. All participants underwent venipuncture between 10:00 a.m. and 3 p.m. Venous blood samples were centrifuged at 2,000 rpm for 10 min, and then stored at -80 °C for genotyping and lipidomic measurements.
MRI Acquisition and Data Processing
MRI data acquisition
MRI data were acquired using a GE signa HDX 3.0 T scanner (General Electric, Milwaukee, USA) with a standard 8-channel head coil at the First Affiliated Hospital of China Medical University, Shenyang, China. Functional images were collected with a gradient-echo planar imaging (EPI-GRE) sequence. The parameters were as follows: TR = 2000 ms, TE = 30 ms, flip angle = 90°, field of view = 240 × 240 mm2, matrix = 64 × 64. Thirty-five axial slices were collected with 3 mm thickness without gap. The scan lasted for 6 min and 40 s, resulting in 200 volumes. Participants were instructed to rest and relax with their eyes closed but remain awake during scanning.
Data preprocessing and ALFF calculation
Preprocessing of all functional images was performed using SPM12 (www.fil.ion.ucl.ac.uk/spm/) and DPARSF . The first 10 time points were discarded for magnetic field stabilization and allowing participants to adapt to the scanning environment. The subsequent preprocessing steps included slice time correction and head motion correction. Each participant’s motion was assessed by means of translation/rotation, and an exclusion criterion (translation > 3 mm, rotation > 3° in each direction) was set. Next, the corrected functional images were normalized to MNI space using the EPI template in SPM12, resampled to 3 mm isotropic voxels, and further smoothed via a Gaussian kernel with a 6 mm full-width at half-maximum. Linear detrending was also performed. Finally, voxel-wise ALFF (0.01–0.08 Hz) maps were calculated for each participant. ALFF (0.01–0.08 Hz) of the BOLD signal, which is considered to be physiologically meaningful and related to regional spontaneous neural activity , was used to identify regional cerebral function.
Polygenic risk scores (PRS)
Genotyping and quality control
Venous blood was collected from all the participants of cross-sectional study, and genomic DNA was extracted from whole blood according to standard procedures. DNA samples were genome-wide genotyped using the Illumina Global Screening Array-24 v1.0 BeadChip, which provides data for 642,824 fixed genetic variants, in addition to 53,411 customized variants. The quality control parameters used for the exclusion of SNPs were as follows: (1) SNPs with minor allele frequency (MAF) < 1%, (2) SNPs with call rate < 95%, or (3) deviation of a SNP from Hardy–Weinberg equilibrium p < 10–5. We then applied the following quality control criteria to exclude participants: (1) individuals with excessive missingness > 5%, (2) gender mismatch, and (3) an estimation of identity-by-descent > 0.90.
Imputation and calculation of PRS
This process has been previously described . Briefly, genotype imputation was performed by a commercial imputation engine named GenoImpute. We constructed PRS in our sample using the PRSice software (www.PRSice.info) based on common SNP risk effects derived from summary statistics from the international SZ GWAS results conducted by the Psychiatric Genomics Consortium (PGC-SZ) including 33,426 SZ cases and 32,541 controls. P-value-informed clumping was performed with a cutoff of r2 = 0.1 in a 250-kb window in order to account for linkage disequilibrium (LD) among SNPs. PRS were calculated using the following twelve PRSs at different p value thresholds (pT) (0.0001, 0.001, 0.01, 0.02, 0.03, 0.04, 0.05, 0.1, 0.2, 0.3, 0.4, and 0.5) for each study participant.
Bioinformatics enrichment analyses
All the SNPs in SZ-PRS under a certain p value threshold were extracted and mapped to the corresponding genes where they were located based on dbSNP database. A gene list was obtained and uploaded to the online tool DAVID Bioinformatics Resources v6.8 (https://david.ncifcrf.gov/,) [34, 35] for the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG)  pathway analyses. The functions of genes were annotated with three GO terms: biological process (BP), cellular component (CC) and molecular function (MF). Multiple testing corrections were performed with Benjamini method (significance level at 0.05).
Lipidomics profile acquirement.
Untargeted metabolomics analysis were conducted in this study. The measurement details can be found in Additional file 1, including sample preparation, metabolic profiling data acquisition, and data processing. Lipidomics profiles (lipids and lipid-related metabolites) were selected from untargeted metabolomics data to conduct subsequent statistical analysis.
Longitudinal ALFF alterations
Two-way analysis of variance with a repeated measures design was implemented in SPM with group (SZ, GHR or HC) as a between-subject factor, time (baseline and follow-up) as within-subject factor, and age, gender and inter-MRI interval as covariates. Clusters that had main effects of group, time, and especially group-by-time interaction on ALFF values were identified. Significance was set at voxel-level inference of P < 0.05 with Gaussian random field (GRF) correction for cluster-level inference of P < 0.05.
For regions with significant group-by-time interaction, ALFF values were extracted for additional post-hoc analyses. At baseline and follow-up, general linear model was performed among SZ, GHR and HC, with age and gender as covariates. Baseline to follow‐up comparisons were performed at SZ, GHR or HC by the paired t-test. Statistical significance was set to P < 0.05.
The influence of polygenic risk and lipometabolism on ALFF alterations
ANOVAs (analyses of variance) or chi-square tests were used to examine participants’ demographic characteristics, clinical characteristics, SZ-PRS and lipometabolism accordingly. Results were considered significant at P < 0.05.
For participants in this portion of the study, ALFF values were extracted from the regions with significant group by time interaction in the longitudinal study for regression and correlation analyses with SZ-PRS and lipid profiles.
SZ-PRS and ALFF alterations
In SZ and GHR, random forest was used to identify the SZ-PRS that could predict significant ALFF alterations identified in the longitudinal study. Random forest is an ensemble method composed of a number of decision trees for regression prediction. Five-fold cross-validation was performed to validate regression results. Subsequently, the prediction model’s performance was evaluated based on mean absolute error (MAE, the average absolute error between the predicted values and the real value) and Pearson correlation coefficient (the correlation between the predicted ALFF values and the actual values). For each SZ-PRS, we separately trained a prediction model in SZ and GHR to predict ALFF values in regions with significant group-by-time interaction in the longitudinal study. A total of 24 models were trained and evaluated.
To further explore the relationships between predictive SZ-PRSs and significant ALFF alterations, partial correlation analyses, controlling for age and sex, were performed in SZ or GHR. Results were considered statistically significant at P < 0.05.
Lipometabolism and ALFF alterations
Similar to above, regression analyses were performed to identify the lipids and their metabolites that strongly predicted significant ALFF alterations in SZ and GHR. The lipid profiling consists of 474 lipids or lipid-related metabolites that can be divided into 19 lipid species. For each lipid species, we separately trained a prediction model in SZ and GHR to predict ALFF values in regions with significant group-by-time interaction. A total of 38 models were trained and evaluated. Further details can be found in Additional file 1.
To further explore the relationships between predictive lipid species and significant ALFF alterations, partial correlation analyses were performed in SZ or GHR, with age, gender and BMI as covariates. Results were considered statistically significant at P < 0.05.
We included two parts in this study. First, we performed a longitudinal neuroimaging study of ALFF in SZ and GHR to identify the key regions with differentiating neural features (Fig. 1A). After that, using divergent ALFF alterations in SZ and GHR across time, we then examined the influence of polygenic risk and lipometabolism on neural function and their relationship with neural function in SZ and GHR in cross-sectional designs (Fig. 1B).
Demographic and Clinical Data
Demographic and clinical information for individuals in the longitudinal study sample is presented in Table 1. There were no significant differences in age among the SZ, GHR, and HC groups at baseline. Significant differences were observed in age at follow-up (F = 3.9, P < 0.05, Table 1). Significant differences were observed in gender (χ2 = 12.82, P < 0.05, Table 1). Significant differences were observed in BPRS scores among the three groups at baseline and follow-up (Table 1).
Demographic and clinical information for individuals in the cross-sectional study sample is presented in Table 2. There were significant differences in age among the SZ, GHR, and HC groups (F = 11.18, P < 0.001). Significant differences were observed in gender (χ2 = 13.57, P = 0.001, Table 2). Significant differences were observed in BPRS scores among the three groups (Table 2).
Longitudinal ALFF alterations
Group-by-time interaction effect on ALFF was identified in left medial orbital frontal cortex (MOF) (Fig. 2A). Post-hoc analyses revealed that at baseline both SZ and GHR groups had increased ALFF in left MOF compared to HC. At follow-up, compared to HC, increased ALFF in left MOF persisted in SZ but not in GHR (Fig. 2B). Notably, compared to HC, GHR did not significantly differ in ALFF in left MOF at follow-up, but normalized across time (Fig. 2B). In SZ, ALFF in left MOF did not significantly change from baseline to follow-up (P > 0.05). In GHR, ALFF in left MOF significantly decreased from baseline to follow-up, while ALFF in left MOF significantly increased in HC.
The influence of polygenic risk and lipometabolism on ALFF alterations
SZ-PRS and ALFF in left MOF
In SZ, the prediction model for SZ-PRS at pT 0.0001 (MAE = 0.153, r = 0.25, P < 0.05), pT 0.2 (MAE = 0.152, r = 0.25, P < 0.05) and pT 0.3 (MAE = 0.145, r = 0.24, P < 0.05) performed relatively better. In GHR, the prediction model of SZ-PRS at pT 0.001 (MAE = 0.121, r = 0.31, P < 0.05) had relatively better performance. Prediction models for other SZ-PRS had poor performance (P > 0.05). SZ-PRS that had better prediction performance were identified as predictive SZ-PRS for SZ and GHR, respectively. Prediction performance for different SZ-PRS are summarized in Additional file 1 (Table S1).
Significant positive correlation was found between left MOF ALFF and predictive SZ-PRS at pT 0.0001 (r = 0.349, P < 0.005) in SZ. There was no significant correlation between left MOF ALFF and predictive SZ-PRS at pT 0.001 in GHR (P > 0.05).
Bioinformatics enrichment analyses were performed for predictive SZ-PRS at pT 0.0001, which was positively correlated with left MOF in SZ as described above. SZ-PRS genes were significantly enriched in five GO terms of cellular component (Fig. 3A), involved in postsynaptic density, postsynaptic membrane, cell junction, dendritic spine and plasma membrane and were mainly located within neuronal structures, particularly cell membranes. Analyses also indicated the significance of two enriched KEGG pathways (Fig. 3B) including AMP-activated protein kinase (AMPK) signaling pathway and axon guidance.
Lipometabolism and ALFF in left MOF
In SZ and GHR, all models had good prediction performance (P < 0.05). The five lipid species with the least MAE were identified as predictive lipid species for SZ and GHR, respectively. The predictive lipid species for SZ were coenzyme, phosphatidylinositol, lysophosphatidylcholine, sphingomyelin, and ceramide (Table 3). For GHR, these predictive lipid species were fatty acids, ceramide, sphingosines, phosphatidylserine and diglyceride (Table 3). Prediction performance of other lipid species can be found in Additional file 1 (Table S2).
Partial correlation between ALFF values in left MOF and mean lipid levels were performed in SZ and GHR for their respective predictive lipid species. In SZ, there was no significant correlation between ALFF in left MOF and mean levels of these five species (P > 0.05). In GHR, left MOF ALFF and mean fatty acid level had significant negative correlation (r = -0.302, P < 0.05) (Fig. 4).
In this study, we used a guided and multimodal approach to examine the longitudinal alterations in neural function as measured by ALFF in GHR and SZ and then the influence of polygenic risk and lipometabolism on divergent ALFF alterations in these two groups to understand vulnterability and resiliency to SZ. We first performed a longitudinal study of ALFF alterations in SZ and GHR to determine key regions with differentiating neural features between the two groups. Differentiating features would reflect disease vulnerability that significantly contribute to disease state and resiliency to SZ. Our longitudinal study found that left MOF appears to be a key region in differentiating SZ and nonpsychotic GHR. We then examined the effects of polygenic risk and lipometabolism on left MOF ALFF, as well as their relationship to ALFF in left MOF. In SZ, SZ-PRS at pT 0.0001 predicted left MOF ALFF and had significant positive correlation with it, and bioinformatics enrichment analyses revealed that genes involved in neuron structure and mainly relate to cell membranes may drive this correlation. Correspondingly, lipodomics showed that lipid species involved in cell membranes best predicted left MOF in SZ. In GHR, fatty acids best predicted ALFF in left MOF and had significant negative correlation with left MOF ALFF. These lipid differences may reflect vulnerability and resiliency to SZ. Interestingly, the longitudinal study found normalization of left MOF ALFF in GHR. Taken together with fatty acid findings in GHR, fatty acids may mediate normalization of ALFF in left MOF. Further studies are needed for more definitive conclusions regarding the interplay between polygenic risk, lipometabolism, and left MOF ALFF.
The study of brain development encompasses evaluation of the structural, functional, and network-level changes that occur across the lifespan. Increasing evidence suggested that the brain achieve a maximum level of maturity around the age of 20–30 years, and then begin to decline [37,38,39,40,41,42]. The longitudinal findings herein support divergent neural trajectories in SZ and GHR during development and across the lifespan. Specifically, the left MOF appears to be a region of divergent ALFF trajectory in SZ and GHR. Left MOF, part of prefrontal cortex (PFC), is thought to contribute to higher cognitive functions . The shared ALFF increases in left MOF in SZ and GHR at baseline could represent shared genetic liability for SZ and may be a neural marker of disease vulnerability. Prefrontal deficits have been previously shown in both SZ and GHR at early ages [10, 44, 45]. Conversely, normalization of ALFF in left MOF at follow-up may reflect resiliency or compensatory mechanisms in GHR individuals who do not develop SZ. Altogether, our findings suggest that left MOF is a critical brain region in understanding the developmental processes leading to SZ.
Consistent with the known heritability of SZ, this study found genetic influences on neural measures in SZ. We found that ALFF in left MOF was positively correlated with predictive SZ-PRS at pT 0.0001 in SZ. This finding is consistent with previous studies in which higher SZ-PRS were associated with functional alterations of PFC [15, 46, 47]. Bioinformatics enrichment analyses indicated that genes involved in neuron structure, especially its cell membranes, likely influence ALFF in left MOF. Many studies report altered neuron structure in the PFC in SZ, including altered postsynaptic density , postsynaptic membrane , cell junction , dendritic spine density , and plasma membrane  as well as axon guidance [52, 53], in line with our findings. Lipodomic findings further supported the bioinformatics findings. Our analyses showed that phospholipid and ceramide, which are key components of cell membranes, best predicted left MOF ALFF in SZ. Collectively, our findings implicate the importance of cell membrane alterations in the etiology and pathophysiology of SZ.
However in GHR, fatty acids best predicted ALFF in left MOF, different from SZ. Further, fatty acids negatively correlated with left MOF ALFF in GHR. This correlation along with the observed normalization of ALFF in GHR in the longitudinal study suggest fatty acids may contribute to resiliency in GHR. Of note, the lipodomic analyses included baseline measures of a subset of participants in the longitudinal study. Fatty acids play a number of key roles in metabolism, including suppliers of energy and signaling molecules . The passage of fatty acids from the blood to the brain have been previously shown . Reduced levels of fatty acids in plasma and red blood cells have been observed in SZ [23, 56]. Disturbances in PFC fatty acid composition have also been reported in SZ [57, 58], suggesting disturbed PFC fatty acid concentrations as a pathological aspect of SZ. However, little is known about fatty acid alterations in GHR. Nevertheless, there is some evidence to support the role of fatty acids in GHR resilience. Fatty acids have the ability to regulate brain-derived neurotrophic factor (BDNF), which has been associated with the pathophysiology of SZ . Supplementation of fatty acids have been suggested to normalize levels of BDNF and reduce oxidative damage , and improve psychiatric symptoms . Altogether, fatty acids may be involved in resiliency in GHR and counteract or compensate for the effects of genetic vulnerability. Further studies are needed.
There are several limitations in this study. Most SZ participants were taking psychotropic medications at the time of the study, and thus there may be confounding effects as GHR participants were not taking psychotropic medications. The longitudinal sample was only of moderate size. Further studies are needed in a larger and unmedicated sample of SZ and GHR to confirm our results. Moreover, gender and age factors were not controlled in the recruitment of participants to ensure the sample size of longitudinal study, which may confound our findings. We therefore performed analysis to test for interactions between gender and group, and age and group. At baseline and follow-up, there was no significant interaction effect of group x gender and group x age on left MOF ALFF in the longitudinal study. There was also no significant interaction effect of group x gender and group x age on the lipids and on left MOF ALFF in the cross-sectional study. In addition, there was no significant linear correlation between ALFF in left MOF and lipid species that best predicted left MOF ALFF in SZ. There also was no significant correlation between left MOF ALFF and predictive SZ-PRSs in GHR. Possible explanations for these findings include complex relationships among polygenic risk, lipometabolism, and ALFF that cannot be simply accounted for by linear correlations. Therefore, further studies are needed to investigate the complex interplay between SZ-PRS, lipids, and ALFF in left MOF in SZ and GHR.
Altogether, using a guided, multimodal approach, our findings implicate left MOF as a key region with divergent neural trajectory between SZ and GHR. Genomic and lipodomic analyses implicate the importance of membrane genes and lipids in differential ALFF alterations and their progression in left MOF in SZ and GHR, reflecting potential indicators of vulnerability and resiliency to SZ. The findings herein suggest that specific lipids and left MOF function may be important targets for early intervention or prevention in GHR, although further studies are warranted for definitve conclusions. These findings may have important implications for understanding mechanisms underlying vulnerability and resiliency in SZ and contribute to translational efforts for early intervention and prevention in SZ.
Availability of data and materials
The datasets analyzed in the current study are available from the corresponding author on reasonable request.
Genetic high risk for SZ
Amplitude of low frequency function
Polygenic risk score for SZ
Medial orbital frontal cortex
Genome-wide association studies
Single nucleotide polymorphism
Diagnostic and Statistical Manual of Mental Disorders-IV-Text Revision
Statistical Manual of Mental Disorders
Schedule for Affective Disorders and Schizophrenia for School-Age Children-present and Lifetime Version
Magnetic resonance imaging
Brief Psychiatric Rating Scale
Minor allele frequency
Kyoto Encyclopedia of Genes and Genomes
Gaussian random field
Analyses of variance
Mean absolute error
Insel TR. Rethinking schizophrenia. Nature. 2010;468:187–93.
Rapoport JL, Giedd JN, Gogtay N. Neurodevelopmental model of schizophrenia: update 2012. Mol Psychiatry. 2012;17:1228–38.
Woods SW, Choi J, Mamah D. Most at-risk individuals will not develop a mental disorder: the limited predictive strength of risk factors. World Psychiatry. 2021;20:223–4.
Tsuang M. Schizophrenia: genes and environment. Biol Psychiat. 2000;47:210–20.
Paus T, Keshavan M, Giedd JN. Why do many psychiatric disorders emerge during adolescence? Nat Rev Neurosci. 2008;9:947–57.
Pantelis C, Yücel M, Bora E, Fornito A, Testa R, Brewer WJ, et al. Neurobiological markers of illness onset in psychosis and schizophrenia: The search for a moving target. Neuropsychol Rev. 2009;19:385–98.
Zalesky A, Pantelis C, Cropley V, Fornito A, Cocchi L, McAdams H, et al. Delayed Development of Brain Connectivity in Adolescents With Schizophrenia and Their Unaffected Siblings. JAMA Psychiat. 2015;72:900–8.
Ordonez AE, Luscher ZI, Gogtay N. Neuroimaging findings from childhood onset schizophrenia patients and their non-psychotic siblings. Schizophr Res. 2016;173:124–31.
Gogtay N, Hua X, Stidd R, Boyle CP, Lee S, Weisinger B, et al. Delayed White Matter Growth Trajectory in Young Nonpsychotic Siblings of Patients With Childhood-Onset Schizophrenia. Arch Gen Psychiatry. 2012;69:875–84.
Mattai AA, Weisinger B, Greenstein D, Stidd R, Clasen L, Miller R, et al. Normalization of cortical gray matter deficits in nonpsychotic siblings of patients with childhood-onset schizophrenia. J Am Acad Child Adolesc Psychiatry. 2011;50:697–704.
Landek-Salgado MA, Faust TE, Sawa A. Molecular substrates of schizophrenia: homeostatic signaling to connectivity. Mol Psychiatry. 2016;21:10–28.
Kaddurah-Daouk R, McEvoy J, Baillie RA, Lee D, Yao JK, Doraiswamy PM, et al. Metabolomic mapping of atypical antipsychotic effects in schizophrenia. Mol Psychiatry. 2007;12:934–45.
Schizophrenia Working Group of the Psychiatric Genomics C. Biological insights from 108 schizophrenia-associated genetic loci. Nature. 2014;511:421–7.
Lancaster TM, Ihssen N, Brindley LM, Tansey KE, Mantripragada K, O’Donovan MC, et al. Associations between polygenic risk for schizophrenia and brain function during probabilistic learning in healthy individuals. Hum Brain Mapp. 2016;37:491–500.
Walton E, Turner J, Gollub RL, Manoach DS, Yendiki A, Ho BC, et al. Cumulative genetic risk and prefrontal activity in patients with schizophrenia. Schizophr Bull. 2013;39:703–11.
Chen Q, Ursini G, Romer AL, Knodt AR, Mezeivtch K, Xiao E, et al. Schizophrenia polygenic risk score predicts mnemonic hippocampal activity. Brain. 2018;141:1218–28.
Dezhina Z, Ranlund S, Kyriakopoulos M, Williams SCR, Dima D. A systematic review of associations between functional MRI activity and polygenic risk for schizophrenia and bipolar disorder. Brain Imaging Behav. 2019;13:862–77.
Neilson E, Bois C, Clarke TK, Hall L, Johnstone EC, Owens DGC, et al. Polygenic risk for schizophrenia, transition and cortical gyrification: a high-risk study. Psychol Med. 2018;48:1532–9.
Dima D, Breen G. Polygenic risk scores in imaging genetics: usefulness and applications. J Psychopharmacol. 2015;29:867–71.
Zhao YY, Miao H, Cheng XL, Wei F. Lipidomics: Novel insight into the biochemical mechanism of lipid metabolism and dysregulation-associated disease. Chem Biol Interact. 2015;240:220–38.
McEvoy J, Baillie RA, Zhu H, Buckley P, Keshavan MS, Nasrallah HA, et al. Lipidomics reveals early metabolic changes in subjects with schizophrenia: effects of atypical antipsychotics. PLoS ONE. 2013;8:e68717.
Schneider M, Levant B, Reichel M, Gulbins E, Kornhuber J, Müller CP. Lipids in psychiatric disorders and preventive medicine. Neurosci Biobehav Rev. 2017;76:336–62.
Sethom M, Fares S, Bouaziz N, Melki W, Jemaa R, Feki M, et al. Polyunsaturated fatty acids deficits are associated with psychotic state and negative symptoms in patients with schizophrenia. Prostaglandins Leukot Essent Fatty Acids. 2010;83:131–6.
Tang Y, Kong L, Wu F, Womer F, Jiang W, Cao Y, et al. Decreased functional connectivity between the amygdala and the left ventral prefrontal cortex in treatment-naive patients with major depressive disorder: a resting-state functional magnetic resonance imaging study. Psychol Med. 2013;43:1921–7.
Orešič M, Seppänen-Laakso T, Sun D, Tang J, Therman S, Viehman R, et al. Phospholipids and insulin resistance in psychosis: a lipidomics study of twin pairs discordant for schizophrenia. Genome Med. 2012;4:1.
Leuze C, Aswendt M, Ferenczi E, Liu CW, Hsueh B, Goubran M, et al. The separate effects of lipids and proteins on brain MRI contrast revealed through tissue clearing. Neuroimage. 2017;156:412–22.
Zang YF, He Y, Zhu CZ, Cao QJ, Sui MQ, Liang M, et al. Altered baseline brain activity in children with ADHD revealed by resting-state functional MRI. Brain Dev. 2007;29:83–91.
Zuo XN, Xing XX. Test-retest reliabilities of resting-state FMRI measurements in human brain functional connectomics: a systems neuroscience perspective. Neurosci Biobehav Rev. 2014;45:100–18.
Hoptman MJ, Zuo XN, Butler PD, Javitt DC, D’Angelo D, Mauro CJ, et al. Amplitude of low-frequency oscillations in schizophrenia: a resting state fMRI study. Schizophr Res. 2010;117:13–20.
Tang Y, Chen K, Zhou Y, Liu J, Wang Y, Driesen N, et al. Neural activity changes in unaffected children of patients with schizophrenia: A resting-state fMRI study. Schizophr Res. 2015;168:360–5.
Yan CG, Zang YF. DPARSF: A MATLAB Toolbox for “Pipeline” Data Analysis of Resting-State fMRI. Front Syst Neurosci. 2010;4:13.
Cordes D, Haughton VM, Arfanakis K, Carew JD, Turski PA, Moritz CH, et al. Frequencies contributing to functional connectivity in the cerebral cortex in “resting-state” data. Am J Neuroradiol. 2001;22:1326–33.
Chang M, Womer FY, Gong X, Chen X, Tang L, Feng R, et al. Identifying and validating subtypes within major psychiatric disorders based on frontal-posterior functional imbalance via deep learning. Mol Psychiatry. 2021;26:2991–3002.
da Huang W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37:1–13.
da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.
Somerville LH. Searching for Signatures of Brain Maturity: What Are We Searching For? Neuron. 2016;92:1164–7.
Lebel C, Walker L, Leemans A, Phillips L, Beaulieu C. Microstructural maturation of the human brain from childhood to adulthood. Neuroimage. 2008;40:1044–55.
Sowell ER, Peterson BS, Thompson PM, Welcome SE, Henkenius AL, Toga AW. Mapping cortical change across the human life span. Nat Neurosci. 2003;6:309–15.
Walhovd KB, Fjell AM, Reinvang I, Lundervold A, Dale AM, Eilertsen DE, et al. Effects of age on volumes of cortex, white matter and subcortical structures. Neurobiol Aging. 2005;26:1261–70 (discussion 75-8).
Kochunov P, Williamson DE, Lancaster J, Fox P, Cornell J, Blangero J, et al. Fractional anisotropy of water diffusion in cerebral white matter across the lifespan. Neurobiol Aging. 2012;33:9–20.
Dosenbach NU, Nardos B, Cohen AL, Fair DA, Power JD, Church JA, et al. Prediction of individual brain maturity using fMRI. Science. 2010;329:1358–61.
du Boisgueheneuc F, Levy R, Volle E, Seassau M, Duffau H, Kinkingnehun S, et al. Functions of the left superior frontal gyrus in humans: a lesion study. Brain. 2006;129:3315–28.
Gogtay N, Greenstein D, Lenane M, Clasen L, Sharp W, Gochman P, et al. Cortical brain development in nonpsychotic siblings of patients with childhood-onset schizophrenia. Arch Gen Psychiatry. 2007;64:772–80.
Unschuld PG, Buchholz AS, Varvaris M, van Zijl PC, Ross CA, Pekar JJ, et al. Prefrontal brain network connectivity indicates degree of both schizophrenia risk and cognitive dysfunction. Schizophr Bull. 2014;40:653–64.
Kauppi K, Westlye LT, Tesli M, Bettella F, Brandt CL, Mattingsdal M, et al. Polygenic Risk for Schizophrenia Associated With Working Memory-related Prefrontal Brain Activation in Patients With Schizophrenia and Healthy Controls. Schizophr Bull. 2014;41:736–43.
Walton E, Geisler D, Lee PH, Hass J, Turner JA, Liu J, et al. Prefrontal inefficiency is associated with polygenic risk for schizophrenia. Schizophr Bull. 2014;40:1263–71.
Catts VS, Derminio DS, Hahn CG, Weickert CS. Postsynaptic density levels of the NMDA receptor NR1 subunit and PSD-95 protein in prefrontal cortex from people with schizophrenia. NPJ Schizophr. 2015;1:15037.
Lisowski P, Wieczorek M, Goscik J, Juszczak GR, Stankiewicz AM, Zwierzchowski L, et al. Effects of chronic stress on prefrontal cortex transcriptome in mice displaying different genetic backgrounds. J Mol Neurosci. 2013;50:33–57.
Glantz LA, Lewis DA. Decreased dendritic spine density on prefrontal cortical pyramidal neurons in schizophrenia. Arch Gen Psychiatry. 2000;57:65–73.
Bauer D, Gupta D, Harotunian V, Meador-Woodruff JH, McCullumsmith RE. Abnormal expression of glutamate transporter and transporter interacting molecules in prefrontal cortex in elderly patients with schizophrenia. Schizophr Res. 2008;104:108–20.
Gilabert-Juan J, Sáez AR, Lopez-Campos G, Sebastiá-Ortega N, González-Martínez R, Costa J, et al. Semaphorin and plexin gene expression is altered in the prefrontal cortex of schizophrenia patients with and without auditory hallucinations. Psychiatry Res. 2015;229:850–7.
Harris LW, Lockstone HE, Khaitovich P, Weickert CS, Webster MJ, Bahn S. Gene expression in the prefrontal cortex during adolescence: implications for the onset of schizophrenia. BMC Med Genomics. 2009;2:28.
de Carvalho C, Caramujo MJ. The Various Roles of Fatty Acids. Molecules. 2018;23(10):2583.
Mitchell RW, On NH, Del Bigio MR, Miller DW, Hatch GM. Fatty acid transport protein expression in human brain and potential role in fatty acid transport across human brain microvessel endothelial cells. J Neurochem. 2011;117:735–46.
Evans D, Parikh V, Khan M, Coussons C, Buckley PF, Mahadik S. Red blood cell membrane essential fatty acid metabolism in early psychotic patients following antipsychotic drug treatment. Prostaglandins Leukot Essent Fatty Acids. 2003;69:393–9.
Taha AY, Cheon Y, Ma K, Rapoport SI, Rao JS. Altered fatty acid concentrations in prefrontal cortex of schizophrenic patients. J Psychiatr Res. 2013;47:636–43.
McNamara RK, Jandacek R, Rider T, Tso P, Hahn CG, Richtand NM, et al. Abnormalities in the fatty acid composition of the postmortem orbitofrontal cortex of schizophrenic patients: gender differences and partial normalization with antipsychotic medications. Schizophr Res. 2007;91:37–50.
Rao JS, Ertley RN, Lee HJ, DeMar JC Jr, Arnold JT, Rapoport SI, et al. n-3 polyunsaturated fatty acid deprivation in rats decreases frontal cortex BDNF via a p38 MAPK-dependent mechanism. Mol Psychiatry. 2007;12:36–46.
Wu A, Ying Z, Gomez-Pinilla F. Dietary omega-3 fatty acids normalize BDNF levels, reduce oxidative damage, and counteract learning disability after traumatic brain injury in rats. J Neurotrauma. 2004;21:1457–67.
Messamore E, McNamara RK. Detection and treatment of omega-3 fatty acid deficiency in psychiatric practice: Rationale and implementation. Lipids Health Dis. 2016;15:25.
The authors thank the members of the patient and healthy participants for taking part in the study.
The authors were supported by research grants from the National Science Fund for Distinguished Young Scholars (81725005 to F.W.), National Natural Science Foundation of China (62176129 to X.Z.), National Key Research and Development Program of China (2018YFC1311600 and 2016YFC1306900 to Y.T., 2022YFC2405603 to X.Z.), Jiangsu Provincial Department of Human Resources and Social Security (2022ZB425 to J.D.), NSFC-Guangdong Joint Fund (U20A6005 to F.W.), Jiangsu Provincial Key Research and Development Program (BE2021617 to F.W.), Liaoning Revitalization Talents Program (XLYC1808036 to Y.T.), China Postdoctoral Science Foundation (2022M721681 to J.Z.)
Ethics approval and consent to participate
The study was approved by the Medical Science Research Ethics Committee of the First Affiliated Hospital of China Medical University (approval reference number 25–1). All participants provided written informed consent by themselves or by their parents/guardians if they were under 18 years old after a complete description of the study. All methods were carried out in accordance with relevant guidelines and regulations.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Main effects of group on ALFF among SZ, GHR and HC. Figure S2. Main effects of time on ALFF among SZ, GHR and HC between baseline and follow-up. Figure S3. Group-by-time interaction on ALFF among SZ, GHR and HC in cross-sectional study. Figure S4. Significant differences in SZ-PRS among groups. Table S1. The prediction results of SZ-PRSs for ALFF values of left MOF in SZ or GHR. Table S2. The prediction results of other lipid species for ALFF values of left MOF in SZ or GHR.
About this article
Cite this article
Duan, J., Gong, X., Womer, F.Y. et al. Neurodevelopmental trajectories, polygenic risk, and lipometabolism in vulnerability and resilience to schizophrenia. BMC Psychiatry 23, 153 (2023). https://doi.org/10.1186/s12888-023-04597-z