Shiwan Tao1
Yamin Zhang2,3
Qiang Wang1
Chunxia Qiao1
Wei Deng2,3
Sugai Liang2,3
Jinxue Wei1
Wei Wei2,3
Hua Yu2,3Xiaojing Li2,3
Mingli Li1
Wanjun Guo2,3Xiaohong Ma1
Liansheng Zhao1
Tao Li2,3,4*- 1Mental Health Center and Psychiatric Laboratory, West China Hospital of Sichuan University, Chengdu, Sichuan, China
- 2Department of Neurobiology, Affiliated Mental Health Center & Hangzhou Seventh People’s Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang, China
- 3NHC and CAMS Key Laboratory of Medical Neurobiology, MOE Frontier Science Center for Brain Science and Brain-Machine Integration, School of Brain Science and Brain Medicine, Zhejiang University, Hangzhou, Zhejiang, China
- 4Guangdong-Hong Kong-Macao Greater Bay Area Center for Brain Science and Brain-Inspired Intelligence, Guangzhou, China
Emerging evidence has demonstrated overlapping biological abnormalities underlying schizophrenia (SCZ), bipolar disorder (BP), and major depressive disorder (MDD); these overlapping abnormalities help explain the high heterogeneity and the similarity of patients within and among diagnostic categories. This study aimed to identify transdiagnostic subtypes of these psychiatric disorders based on lipidomics abnormalities. We performed discriminant analysis to identify lipids that classified patients (N = 349, 112 with SCZ, 132 with BP, and 105 with MDD) and healthy controls (N = 198). Ten lipids that mainly regulate energy metabolism, inflammation, oxidative stress, and fatty acylation of proteins were identified. We found two subtypes (named Cluster 1 and Cluster 2 subtypes) across patients with SCZ, BP, and MDD by consensus clustering analysis based on the above 10 lipids. The distribution of clinical diagnosis, functional impairment measured by Global Assessment of Functioning (GAF) scales, and brain white matter abnormalities measured by fractional anisotropy (FA) and radial diffusivity (RD) differed in the two subtypes. Patients within the Cluster 2 subtype were mainly SCZ and BP patients and featured significantly elevated RD along the genu of corpus callosum (GCC) region and lower GAF scores than patients within the Cluster 1 subtype. The SCZ and BP patients within the Cluster 2 subtype shared similar biological patterns; that is, these patients had comparable brain white matter abnormalities and functional impairment, which is consistent with previous studies. Our findings indicate that peripheral lipid abnormalities might help identify homogeneous transdiagnostic subtypes across psychiatric disorders.
1 Introduction
Schizophrenia (SCZ), bipolar disorder (BP), and major depressive disorder (MDD) are three common psychiatric disorders with a heavy disease burden (Vigo et al., 2016). In clinical practice, it is an issue that the boundaries among the various diagnoses are not clearly distinct from each other. Patients with different diagnoses usually present baffling similarities to each other, such as the “with psychotic features” item of BP and MDD diagnoses and the emotional dysfunctions in SCZ. This is partly because rather than diagnostic objective criteria or biological markers, the current psychiatric diagnosis nosology relies on descriptive information elicited from self-report history and clinical observation (Scadding, 1996; Craddock and Mynors-Wallis, 2014; Heckers and Kendler, 2020). This hampers the diagnostic accuracy of these psychiatric disorders.
Increasing evidence has demonstrated that there are overlapping biological characteristics across SCZ, BP, and MDD, such as genetic risk factors (Cross-Disorder Group of the Psychiatric Genomics, 2013; Cross-Disorder Group of the Psychiatric Genomics et al., 2013; Ruderfer et al., 2014; Cross-Disorder Group of the Psychiatric Genomics Consortium. Electronic address and Cross-Disorder Group of the Psychiatric Genomics, 2019; Andlauer et al., 2021), brain structure and functional abnormalities (Meda et al., 2014; Godwin et al., 2018; Kelly et al., 2018; Li et al., 2018; Favre et al., 2019; van Velzen et al., 2020) and cognitive impairment (Czobor et al., 2007; Reichenberg et al., 2009; Millan et al., 2012; Barch and Sheffield, 2014; Reilly and Sweeney, 2014; Tamminga et al., 2014). All of these findings imply that the diagnostic classes are not distinct entities, and the descriptive diagnosis nosology has fundamental flaws. Therefore, it is necessary to identify natural biological homogeneous subtypes across different psychiatric disorders.
Lipid metabolites are downstream biochemical end products that are more close to phenotypes than genomics and proteomics. As an essential part of systems biology, lipidomics could comprehensively illuminate the lipid metabolic profile of individuals and identify changes related to phenotype (Patti et al., 2012; Zhao et al., 2014). Plasma lipid alterations, therefore, are sensitive and specific to several observed risk factors for psychiatric disorders, including genetic variations, brain white matter (WM) structural abnormalities, and oxidative stress and inflammation. For example, the ABCD1 gene mutation caused very long-chain fatty acid accumulation in brain WM, which led to psychiatric symptoms (Kitchin et al., 1987; Kemp et al., 2016). Plasma lipids, such as triglyceride, were also reported to be associated with brain WM microstructural changes and axonal degeneration (Iriondo et al., 2021). Notably, derived from peripheral essential omega-6 and omega-3 polyunsaturated fatty acids, lipid-derived mediators serve as pro/anti-inflammatory mediators regulating brain inflammation (Laye et al., 2018). Brain tissues are susceptible to oxidative stress due to their high oxygen consumption and unsaturated fatty acid enrichment, which have been reported to be associated with SCZ, BP, and MDD (Salim, 2017; Cobley et al., 2018).
Lipidomics has recently been developed as a powerful tool to investigate the natural characteristics of SCZ, BP, and MDD. The peripheral lipidomics profile alterations of these psychiatric disorders have been pervasively characterized using ultrahigh-performance liquid chromatography-tandem mass spectrometry (UHPLC–MS/MS) technology and have served as promising biomarkers for early diagnosis and clinical outcome prediction (Oresic et al., 2011; Brunkhorst-Kanaan et al., 2019; Zhou et al., 2019; Bot et al., 2020; Hussain et al., 2020; Zhuo et al., 2020; Dickens et al., 2021). Therefore, taking advantage of lipidomics analysis may help identify biologically homogeneous subtypes across these psychiatric disorders. In this study, first, we investigated the peripheral lipidomics profile abnormalities between psychiatric patients (with SCZ, BP, and MDD) and healthy controls (HCs) by discriminant analysis, and identify the most contributory lipids for classification. Then, based on these identified lipids, we further investigated the potential subtypes across SCZ, BP, and MDD by consensus clustering analysis. To comprehensively profile the differences in these potential biological subtypes, we further described and compared the brain WM microstructure and clinical features of these subtypes.
2 Methods
2.1 Participants
All participants, who were right-handed Han Chinese and aged 16–55 years old, were interviewed by at least two trained psychiatrists using the Structured Clinical Interview of the Diagnostic and Statistical Manual of Mental Disorders, 4th Edition, Text Revision (DSM-IV-TR)—Patient Version (SCID-P). A total of 547 participants (112 patients with SCZ, 132 with BP, 105 with MDD, and 198 HCs) were recruited from West China Hospital of Sichuan University between 2014 and 2019. The inclusion criteria for patients were as follows: 1) fulfilment of one of the DSM-IV-TR criteria for SCZ, BP, or MDD; 2) Han Chinese; 3) right-handed; 4) education achievement of more than 6 years; and 5) scores on Wechsler’s intelligence test equal to or higher than 70. The exclusion criteria for patients were as follows: 1) comorbidity with other DSM-IV-TR axis I or axis II disorders (such as alcohol and substance abuse); 2) presence of organic brain diseases, neurological diseases or somatic diseases undergoing drug treatment (such as diabetes); 3) any history of head trauma; 4) any physical therapies, such as electroconvulsive therapy, undergone within the past 6 months before magnetic resonance imaging (MRI) scan; 5) any contraindication to perform MRI scan; 6) pregnant or breastfeeding; and 7) Wechsler’s intelligence test scores less than 70. In this study, all SCZ patients were first-episode and drug-naïve. There were 77 MDD patients and 63 BP patients who were drug-naïve, and 28 MDD and 69 BP patients who were not drug-naïve but had at least a two-week wash-out period.
HCs were enrolled via online and local advertisements. They were screened for any mental disorder by the SCID—Non-Patient Version (SCID-NP). The exclusion criteria for HCs were similar to those for patients. Moreover, HCs with first-degree relatives with DSM-IV-TR axis I or II disorders were excluded.
2.2 Ethical principles
This study abided by the guidelines of the Declaration of Helsinki and was approved by the Institutional Ethics Committee of West China Hospital, Sichuan University. After the study procedure had been fully explained, written informed consent was obtained from all participants and their guardians if participants were less than 18 years old.
2.3 Clinical assessment
We used the Global Assessment of Functioning (GAF) scale to evaluate functional impairment in all patients. The Positive and Negative Syndrome Scale (PANSS), the Young Mania Rating Scale (YMRS), the Hamilton Anxiety Scale (HAMA), and the Hamilton Depression Scale (HAMD) were used to assess symptom severity in patients as appropriate. Clinical features, including onset age, total duration of illness period (TDP), current duration of illness period (CDP), duration of untreated period (DUP), current episode state, BP I or II subtype for BP, and the number of episodes for MDD, were also documented.
2.4 Lipidomics data acquisition and preprocessing
Peripheral blood was collected in EDTA tubes from all participants on the same day they were enrolled in this study. Lipid extraction, UHPLC‒MS/MS analysis, and lipid qualitative and quantitative identification are described in the Supplementary Methods. A total of 7212 lipid features in the positive polarity model and 4,898 lipid features in the negative polarity model were obtained. The lipidomic data were preprocessed by the “statTarget” (version 1.22.0) (Luan et al., 2018) and “MetaboAnalystR” (version 3.0.3) (Pang et al., 2020) packages in R software (version 4.1.0). We performed preprocessing steps as follows: 1) drift signal correction using the quality control-based robust locally estimated scatterplot smoothing (LOESS) signal correction (QC-RLSC) algorithm (Dunn et al., 2011); 2) a quality assurance procedure to remove metabolic features with relative standard deviation (RSD) >20%, which was calculated for all QC samples (Dunn et al., 2011); 3) log2 transformation and ComBat batch effect correction (Johnson et al., 2007); and 4) interquartile range (IQR) data filtering. The quality control results of the lipidomics data are described and depicted in Supplementary Figure S1. After the preprocessing steps, a total of 1,164 lipids remained for discriminant analysis.
2.5 DTI data acquisition and preprocessing
Brain WM microstructural abnormalities were measured by fractional anisotropy (FA) and radial diffusivity (RD) using diffusion tensor imaging (DTI) scans. The FA indicates the underlying characteristics of white matter microstructure, such as the directionality of axonal fibres, diameter, and density (Basser, 1995; Basser and Pierpaoli, 1996). RD is considered an indicator of myelin sheath thickness, reflecting myelin damage (Song et al., 2002; Song et al., 2005). Altered FA or RD in some regions indicated the brain white matter microstructures abnormalities here. The DTI scan parameters are described in the Supplementary Methods. Raw images were processed by MRIcroN (http://www.mricro.com), DTIPrep, and FMRIB Software Library (FSL) (version 5.0.8). The imaging format was converted by MRIcroN, and then the imaging quality was checked by DTIPrep (translation <2 mm, rotation <0.5 mm). Individual images that met the quality control criteria were included for subsequent procedures (6 samples were excluded after checking the imaging quality). The preprocessing steps included motion and eddy current correction, gradient direction reorientation, and brain mask estimation to remove the nonbrain spaces. After calculating diffusion tensor metrics, normalization and linear/nonlinear registration were also performed to allow comparison across participants. Brain regions of interest (ROIs) were defined by the JHU-ICBM-DTI-81 WM labels atlas (n = 48). The z scores of the mean FA and RD in each ROI were calculated for further statistical analysis.
2.6 Statistical analysis
2.6.1 Demographic characteristics
Demographic characteristics including age, sex, educational attainment (years), and body mass index (BMI) of different groups (psychiatric patients and HCs) were compared using the independent test or chi-square test. BMI was calculated as weight divided by height squared (kg/m2). All the analyses above were calculated in R software.
2.6.2 Discriminant analyses for patients and HCs based on lipidomics data
Preprocessed lipids were further analysed using the “mixOmics” (version 6.16.3) package in R software (Rohart et al., 2017). The data were centred on zero mean and unit variance (auto scaling). Principal component analysis (PCA) was used to check the homogeneity of the samples and determine whether QC samples were tightly clustered together. After removing the outliers, we developed a sparse partial least square-discriminant analysis (sPLS-DA) model, a supervised machine learning analysis, to identify the lipids that contributed most to the classification of psychiatric and HC groups. Parameter tuning processes were performed using the tune function to determine the optimal parameters. The performance of the tuning sPLS-DA model obtained with a balanced error rate (BER) was estimated with 7-fold cross validation and repeated 1,000 times. The optimal parameters, including the number of components and variables, were selected when the tuning model had a low classification error rate. The performance of the optimal sPLS-DA model was estimated by using the perf function, with 7-fold cross validation repeated 1,000 times. Evaluated indexes included BER and overall classification error rate (prediction distances were calculated by max, centroids, and Mahalanobis distance) and areas under the receiver operating curve (AUCs). We also performed univariate analysis and two-sample Wilcoxon rank-sum tests to complement the multivariate analysis, followed by false discovery rate (FDR) adjustment. The most important variables (lipids) for differentiating the psychiatric and HC groups satisfied the following cut-off criteria: 1) AUCs of the sPLS-DA model >0.8, 2) variable importance in projection (VIP) scores >1, 3) occurrence frequency of the lipids >0.8 after performing 1,000 times of cross validation, and 4) p value <0.05 after FDR adjustment. A total of 10 significantly altered lipids met all the above criteria to differentiate between patients and HCs.
2.6.3 Identifying lipid-based subtypes utilizing consensus cluster
We developed an unsupervised cluster model to investigate the potential subtypes within the group of psychiatric patients using the data of 10 identified lipids. Consensus partitioning was performed and summarized by the “cola” package (Gu et al., 2021). Features for consensus partitioning were calculated by four methods: standard deviation (SD), median absolute deviation (MAD), coefficient of variation (CV), and ability to correlate to other rows (ATC). Partitioning methods included hierarchical clustering (hclust), k-means clustering (kmeans), partitioning around medoids (pam), and spherical k-means clustering (skmeans). The partitioning step was repeatedly executed 50 times for each partitioning method. The mean silhouette score and concordance were calculated to evaluate the cluster models and select the optimal number (k) of subtypes. The SD-skmeans model generated an optimal k of 2.
2.6.4 The differences between lipid-based subtypes across multiple-level data
We compared the differences between the two subtypes in terms of clinical features (including the global functional impairment measured by GAF scale scores; symptoms severity measured by PANSS scores, YMRS scores, HAMA scores, and HAMD scores in patients as appropriate; onset age; TDP; CDP; and DUP), and brain WM microstructural alterations (measured by FA and RD) in R software. A two-sample t test was performed to compare the difference in GAF scores and ROI-based FA and RD data between the subtypes (followed by FDR adjustment). The 48 WM regional FA and RD effect sizes of subtypes (Cohen’s d) were also calculated.
3 Results
3.1 Demographic characteristics
We removed 7 individuals (3 patients and 4 HCs) after the lipidomic data quality control process (Supplementary Figure S1). The demographic characteristics of the remaining 346 psychiatric patients and 194 HCs are described in Table 1. There were no significant differences between patients and HCs in terms of age, sex, or BMI. The mean educational attainment years of participants in the HC group were higher than those of participants in the psychiatric group (p < 0.001).
TABLE 1. Comparison of demographic characteristics between the psychiatric patient and healthy control groups.
3.2 Discriminant analyses for patients and HCs
3.2.1 Choosing optimal parameters from the tuning model
The performance of tuning the sPLS-DA model is displayed inFigure 1A. The balanced classification error rates were decreased when more components were added to the model. In the tuning model, the first two components (composed of 2 lipid features selected from the first component and 20 lipid features selected from the second component) were sufficient to achieve good performance (error rate = 0.046 ± 0.005, 7-fold cross validation repeated 1,000 times).
FIGURE 1. The sPLS-DA model for differentiating psychiatric patients and healthy controls using lipidomic data. (A) Balanced error rates (BERs) decreased when more components were added to the tuning sPLS-DA model. Here, the first 2 components (light blue line) were sufficient to achieve good performance (error rate = 0.046 ± 0.005, 1,000×7-fold cross-validation), and the optimal features of each component are indicated as a diamond plot. (B) The sPLS-DA sample plot with ellipse circles indicating the 95% confidence interval. The first two components of the sPLS-DA model differentiating the psychiatric patient group (orange triangle) from the HC group (blue circle). (C) Classification performance per component (overall and BER) for three prediction distances using repeated cross-validation (1,000×7-fold). All the classification error rates were lower than 0.06. (D) The ROC curve of the sPLS-DA model, and the AUC = 0.986. sPLS-DA, sparse partial least squares discriminant analysis; ROC, receiver operating curve; AUC, area under the receiver operating curve; HC, healthy control group; Psychosis, psychiatric patient group.
3.2.2 Identifying contributing lipids for the classification of psychiatric patients and HCs
Figure 1B displays the sPLS-DA sample plot. The first two components accurately distinguished psychiatric patients from HCs. Figure 1C displays the BER and overall error rates of the two components for three prediction distances (7-fold cross validation, repeated 1,000 times). All classification error rates were less than 0.06 (details in Supplementary Table S1). The receiver operating curve (ROC), as an additional measure that helped reflect the performance of the sPLS-DA model, is depicted in Figure 1D, and the AUC = 0.986. The low classification error rates and high AUC indicate that the previous tuning process led to a final sPLS-DA model that achieved good performance. According to the cut-off criteria mentioned in the methods, 10 lipids were selected (Table 2).
TABLE 2. Identified differential lipids for classifying psychiatric patients and healthy controls.
3.3 Consensus cluster analysis within the group of psychiatric patients
Consensus clustering was performed among the psychiatric patients. The skmeans model generated stable partitions compared to other methods, especially when combined with SD (details in Supplementary Table S2). The confident samples with silhouette scores >0.5 (N = 319) were classified into two stable subtypes named the Cluster 1 and Cluster 2 subtypes (mean silhouette = 0.8; concordance = 0.9). The consensus heatmap (Figure 2A) provides a visual representation of how consistent two samples were in the same subtype. The PCA plot (Figure 2B) also confirmed that the two subtypes were separate from each other.
FIGURE 2. The consensus clustering analysis identified two subgroups within psychiatric patients. (A) The consensus matrix heatmap plot visualized the stability of the two subtypes. The “p1” and “p2” labels refer to the probability of the sample being classified into “class 1” and “class 2,” respectively, after clustering was repeated 50 times. “Prob” refers to the calculated probability of the sample being classified into the corresponding subgroup. The “consensus” legend refers to how consistent two samples were in the same subgroup. (B) The PCA plot confirmed that there were two subgroups of psychiatric patients. The confident samples (silhouette score >0.5) are classified into two subgroups obviously separated from each other, and the ambiguous samples (silhouette score <0.5) are indicated by crosses on the plot.
3.4 The differences in lipid-based subtypes across multiple-level data
3.4.1 Demographic characteristics of the two subtypes
There were no significant differences between the two lipid-based subtypes in terms of demographic characteristics (age, sex, educational attainment, and BMI) (Table 3).
TABLE 3. Comparison of demographic characteristics and functional impairment assessment between the lipid-based subgroups.
3.4.2 The differences in clinical features of the two subtypes
Cluster 1 included 179 patients (52.65%), and Cluster 2 included 140 patients (41.18%). The clinical diagnosis distribution varied between the two subtypes (x2 = 65.81, p < 0.001) (Table 3). Cluster 1 consisted of 33 (18%) patients with SCZ, 59 (33%) patients with BP and 87 (49%) patients with MDD, and Cluster 2 consisted of 65 (46%) patients with SCZ, 64 (46%) patients with BP and 11 (8%) patients with MDD (Figure 3A). A higher proportion of patients with MDD was in Cluster 1 (89%) than in Cluster 2 (11%). In contrast, more patients with SCZ were allocated to Cluster 2 (66%) than to Cluster 1 (34%). Patients with BP were uniformly distributed in Cluster 1 (52%) and Cluster 2 (48%) (Figure 3B). In regard to clinical features, Cluster 2 patients (50.17 ± 13.83) showed significantly lower GAF scores than Cluster 1 patients (54.96 ± 13.08) (t = 2.94, p = 0.0036) (Table 3). The clinical features of schizophrenia did not show any difference between the two clusters (Supplementary Table S3); the distribution of bipolar I and bipolar II disorder differed in the two clusters (x2 = 4.87, p = 0.027), and a higher proportion of bipolar I patients occurred in Cluster 2 (61%) than in Cluster 1 (39%) (Supplementary Table S4). In addition, the HAMA scores of MDD patients with a Cluster 2 (11.36 ± 5.73) subtype were significantly lower than those with a Cluster 1 subtype (16.07 ± 5.61) (t = 2.56, p = 0.024) (Supplementary Table S5). Table 4 provides a general summary schema to summarize the comparison results of the clinical features of SCZ, BP and MDD allocated to the two subtypes.
FIGURE 3. Distribution of clinical diagnoses in the two subgroups. (A) Cluster 1 included 179 psychiatric patients, consisting of 33 (18%) patients with SCZ, 59 (33%) patients with BP and 87 (49%) patients with MDD, and Cluster 2 included 140 psychiatric patients. (B) A higher proportion of MDD patients was present in Cluster 1 (89%) than in Cluster 2 (11%). There were more patients with SCZ in Cluster 2 (66%) than in Cluster 1 (34%). Patients with BP were uniformly distributed in Cluster 1 (52%) and Cluster 2 (48%). SCZ, schizophrenia; BP, bipolar disorder; MDD, depressive disorder. *p < 0.05; **p < 0.01; ***p < 0.001.
TABLE 4. Comparison of clinical features of SCZ, BP and MDD patients between the lipid-based subgroups.
3.4.3 Brain white matter alterations between the two subtypes
After FDR adjustment, patients in the Cluster 2 group showed significantly increased RD (1.169 ± 0.768) compared to those in the Cluster 1 group (0.857 ± 0.771) (Cohen’s d = 0.405; t = −3.591; p.adj = 0.018), mainly along the genu of corpus callosum (GCC) (Figure 4; Supplementary Table S6). Patients within the Cluster 2 subtype showed trends of decreased FA along the fornix (including the column and body of the fornix) (p = 0.018, Cohen’s d = 0.266) and right posterior thalamic radiation (p = 0.032, Cohen’s d = 0.243) and increased FA mainly along the left hippocampus region (p = 0.044, Cohen’s d = 0.226) compared to patients within the Cluster 1 subtype, although significance did not survive FDR adjustment (Supplementary Table S7). Supplementary Table S8 provides the association of identified lipids and brain WM alterations in psychiatric patients.
FIGURE 4. Radial diffusivity (RD) differences between the two subgroups for 48 white matter brain regions that represent the major fasciculi. The colour bar (red–yellow) indicates the mean effect size of the group (Cohen’s d). The genu of the corpus callosum (green arrow) showed significantly increased RD in patients in the Cluster 2 subgroup (Cohen’s d = 0.405; p.adj = 0.018).
4 Discussion
In this study, discriminant analysis identified 10 disease-specific lipids that contribute to the classification of psychiatric patients (including patients with SCZ, BP, and MDD) and HCs. We further found two lipid-based subtypes (named the Cluster 1 and Cluster 2 subtypes) within the psychiatric patients utilizing cluster analysis. The two subtypes differed in clinical features and brain WM abnormalities. The clinical diagnosis distribution significantly differed in the two subtypes: patients with BP were uniformly distributed in the two subtypes, but a higher proportion of patients with MDD (89%) was noted in Cluster 1, and a higher proportion of patients with SCZ (66%) was noted in Cluster 2. Patients in Cluster 2 showed significantly lower GAF scores than those in Cluster 1. Moreover, the patients within Cluster 2 showed significantly increased RD in the GCC, decreased FA trends in the fornix and posterior thalamic radiation, and increased FA trend in the hippocampus.
Patients within the Cluster 2 subtype mainly consisted of those with SCZ and BP (a total of 92%). Previous studies have indicated that SCZ and BP are characterized by similar biological patterns, such as high genetic correlation (Cross-Disorder Group of the Psychiatric Genomics et al., 2013) and comparable WM abnormalities. A large-scale meta-analysis has reported that patients with SCZ/BP (but not in MDD) shared limbic system (such as the fornix) abnormalities (Koshiyama et al., 2020), and posterior thalamic contraction (Mamah et al., 2016). In this study, the interesting constituent ratio of clinical diagnosis, and the decreased FA trend along the fornix and posterior thalamic region in patients within the Cluster 2 subtype were supported and consistent with previous findings. We also observed that bipolar I and bipolar II disorder distributions differed between the two subtypes. This finding is consistent with previous studies that found biological heterogeneity between bipolar I and bipolar II disorder (Charney et al., 2017; Huang et al., 2022). HAMA scores of MDD patients differed in the two subtypes. However, considering that only 11 MDD patients were allocated to Cluster 2, the small sample size may not satisfy the statistical power. Other clinical features of SCZ, BP and MDD showed no differences between the two lipid-based subtypes, which support that there are mismatch boundaries between biological subtypes and clinical diagnosis based on descriptive data. The differential findings of the two subtypes support the high similarity among SCZ and BP patients.
In addition to the lipid-based biological pattern, Cluster 2 patients also presented differential WM abnormalities measured by RD mainly along the GCC region and lower GAF scores. Brain WM abnormalities in the corpus callosum have been widely and consistently reported across several psychiatric disorders by meta-analyses, especially in the GCC of patients with SCZ (Kelly et al., 2018; Koshiyama et al., 2018; Favre et al., 2019; van Velzen et al., 2020). Lower GAF scores indicate severer psychological, social and occupational functioning impairment. The GCC is the bend of the anterior corpus callosum; thus, it facilitates prefrontal interhemispheric connectivity and relates to social competence, planning and memory performance, etc. (Paul et al., 2007). We speculate that the greater functioning impairment of patients within the Cluster 2 subtype are potentially the consequences of structural abnormalities in the GCC. As RD is a specific index reflecting the demyelination or morphology abnormalities of fibre tracts (Song et al., 2002; Song et al., 2005), the significantly elevated RD in the GCC region of Cluster 2 patients might reflect greater brain WM lesions here. In summary, these findings indicate that the lipid-based subtypes across psychiatric disorders also showed differential multiple-level biological characteristics.
The 10 identified lipids from the sPLS-DA model relate to several abnormal physiological processes, including inflammation and oxidative stress, brain structural or functional abnormality regulation, and metabolic deterioration. 10-Nitro-9Z,12Z-octadecadienoic acid (nitrolinoleic acid, LNO2) is rich in human plasma and red cell membranes. It acts as a lipid-derived mediator in activating antioxidant signalling pathways (Kalyanaraman, 2004; Koutoulogenis and Kokotos, 2021). LNO2 also exhibits robust cell signalling activities as an anti-inflammatory (Coles et al., 2002; Schopfer et al., 2005; Wright et al., 2006; Koutoulogenis and Kokotos, 2021). In this study, decreased plasma LNO2 might indicate the vulnerable anti-inflammatory status of psychiatric patients. In addition, there are several other identified lipids associated with inflammation and oxidative stress. Of note, change in hexadecanedioic acid level was reported to be related to inflammatory status, and it contributed to the classification of SCZ and HCs in previous studies (Cui et al., 2020; Qian et al., 2021). The diacylglyceryl-N,N,N-trimethylhomoserine (DGTS)16:0/18:1 level is considered a biomarker reflecting low oxidative stability among wheat varieties (Wei et al., 2021). A meta-analysis has revealed that first-episode psychiatric patients exhibited a proinflammatory and vulnerable antioxidant status (Fraguas et al., 2019). Recently, a proposed hypothesis illuminated immune/inflammatory-mediated alteration of brain WM in the limbic system as the main pathophysiological mechanism of psychiatric disorders (Magioncalda and Martino, 2022). In this study, DGTS 16:0/18:1 and hexadecandioic acid positively related to RD of GCC and FA of left hippocampus. These above-altered lipids might indicate inflammation and oxidative stress imbalance in psychiatric patients.
Although there is insufficient evidence, previous studies have implicated that the altered 12-tridecynoic acid and 4-amino-3-methylbutanoic acid levels might be related to the regulation of brain structural abnormalities. The Wnt signalling pathways are important in modulating synapse growth and synaptic plasticity in humans, and altered Wnt signalling was documented in patients with SCZ and BP (Tabares-Seisdedos and Rubenstein, 2009; Hoseth et al., 2018). 12-Tridecynoic acid is one of the lipids that participates in the fatty acylation/deacylation of Wnt proteins (Gao and Hannoush, 2014; Torres et al., 2019), which are necessary for their biofunction (Willert et al., 2003; Rios-Esteves et al., 2014). However, there is no direct evidence linking the changes in fatty acid levels to Wnt protein activation. In this study, the effect of increased 12-tridecynoic acid levels on the brain structure is unknown, and further studies could perhaps investigate the relationship between them by evaluating Wnt signalling pathways. 4-Amino-3-methylbutanoic acid is a 3-substituted γ-aminobutyric acid (GABA) analogue with greater affinity for GABA receptors in the human brain (Nicholson et al., 1979). Moreover, it could raise GABA levels by increasing L-glutamic acid decarboxylase (GAD) activity in the mouse brain and produce an anticonvulsant effect (Silverman et al., 1991; Taylor et al., 1992). However, there is also no evidence linking decreased 4-amino-3-methylbutanoic acid concentrations with the function of central GABAergic neurons, which requires more research. These peripheral lipid alterations may provide clues and broaden our understanding of the mechanisms underlying brain structure abnormalities, which is one of the main pathogenic mechanisms of psychiatric disorders.
Unhealthy dietary and behaviour patterns have recently been noted as risk factors for the metabolic deterioration of patients with SCZ, BP, and MDD (Beyer and Payne, 2016; Vancampfort et al., 2017). A previous randomized crossover trial reported that the high fiber consumption dietary intervention decreased plasma 9,12-octadecadienal level. It is considered to be involved in mediating the positive effect of a healthy diet on maintaining satiety and preventing obesity (Lankinen et al., 2011). In this study, increased 9,12-octadecadienal level in psychiatric patients might reflect the unhealthy dietary pattern (such as low consumption of fiber and fruit) of these patients (Dipasquale et al., 2013). Caprylic acid is important in regulating food intake behaviour by esterifying ghrelin, which is a key peptide hormone with orexigenic biofunction (Kojima et al., 1999; Kojima and Kangawa, 2002; Delporte, 2013). Previous randomized controlled trial studies have reported that ingestion of caprylic acid helps stimulate food intake behaviour, and has been used to treat anorexia nervosa (Kawai et al., 2017) and cachectic patients (Ashitani et al., 2009). Disordered eating behaviours are common among SCZ, BP (such as binge eating, food cravings, and night eating), and MDD (emotional and external eating) patients, which were occurred in the initial onset and cannot all be attributed to the side effects of drug treatment (Paans et al., 2018; Stogios et al., 2020; Sankaranarayanan et al., 2021). In this study, the increased caprylic acid level might provide clues about the disordered eating behaviour among psychiatric patients. Apart from attention to dietary patterns and disordered eating behaviour, an unhealthy behaviour pattern is another important risk factor for metabolic deterioration. Of note, 20-oxo-22,23,24,25,26,27-hexanorvitamin D3 is only synthesized by skin tissue in humans through ultraviolet B (UVB) induced physicochemical processes (Slominski et al., 2012). Since UVB is essential in the synthetic process, it is reasonable to infer that the significantly decreased 20-oxo-22,23,24,25,26,27-hexanorvitamin D3 levels in psychiatric patients may be attributed to lower sunlight exposure, which is associated with unhealthy behavioural patterns (such as sedentary behaviour). A previous meta-analysis reported that patients with SCZ, BP and MDD have significant sedentary behaviour (average 476 min per day) during waking hours and low activity (38.4 min per day) (Vancampfort et al., 2017). Above all, the alterations in 9,12-octadecadienal, caprylic acid, and 20-oxo-22,23,24,25,26,27-hexanorvitamin D3 levels might reflect the unhealthy dietary and behavioural pattern of the mechanism that underlying metabolic deterioration of psychiatric patients.
There are some limitations to this study. First, as antipsychotic drugs affect lipid metabolism, we tried our best to recruit drug-naïve patients. In this study, all the recruited schizophrenia patients were first-episode and drug-naïve, however, drug-naïve bipolar disorder patients were hard to recruit due to the diagnosis delay and high misdiagnosis rate (Culpepper, 2014; Fritz et al., 2017). We recruited 63 (47.8%) drug-naïve BP patients and set at least 2 weeks wash-out period for other BP patients (the median current duration of illness period was 2 months). When it comes to MDD patients, the main treatment strategies are SSRI/SNRI drugs, even so, we recruited 77 (77.3%) drug-naïve patients and also set at least 2 weeks wash-out period for other MDD patients. Then, because sample collection at a single center with lower variability may restrict the generalization of these findings. We will conduct independent sample validation in the future study to make these findings more robust and convincing. Moreover, the biological functions of some identified lipids are attractive, such as LNO2 and caprylic acid. Although previous studies have evidenced their biological function in psychiatric patients, further studies could better elucidate the effects of these lipids on psychiatric diseases. For example, adding inflammatory factors examination, and the questionnaire about the dietary and behaviour patterns.
5 Conclusion
In conclusion, our findings suggested that peripheral blood lipidomic profile alterations could help identify homogeneous transdiagnostic subtypes across psychiatric disorders consisting of SCZ, BP and MDD. One of the subtypes that mainly consisted of patients with SCZ and BP represented more severe brain WM abnormalities and functional impairments. It is suggested that lipid-based subtypes might help identify patients with differential biological characterizations.
Data availability statement
The datasets presented in this article are not readily available because According to the Regulations of the People’s Republic of China on the Administration of Human Genetic Resources, which came into effect on 1 July 2019, those who provide or open to use human genetic resources information to foreign organizations, individuals and institutions established or actually controlled by them shall submit the information to the administrative department of science and technology under The State Council for the record and submit backup information. The plasma sample analyzed in this article are also within the limits of this regulation. Requests to access the datasets should be directed to litaozjusc@zju.edu.cn.
Ethics statement
The studies involving human participants were reviewed and approved by the Institutional Ethics Committee of West China Hospital, Sichuan University. Written informed consent to participate in this study was provided by all participants. For participants who were less than 18 years old, their legal guardian/next of kin also provided written informed consent.
Author contributions
ST, YZ, and QW designed the analysis strategies. CQ, WD, HY, XL, ML, WG, and XM contributed to clinical data collection and assessment. JW and LZ contributed to lipidomics data acquisition. SL and WW contributed to guiding the imaging data analysis. ST performed the statistical analysis. ST and YZ interpreted the results and wrote the manuscript. TL overviewed the whole experiment and approved the final manuscript.
Funding
This work was partly funded by the National Nature Science Foundation of China Key Project (TL and PS, 81920108018); the Key R&D Program of Zhejiang (TL 2022C03096); the Special Foundation for Brain Research from Science and Technology Program of Guangdong (2018B030334001); the Project for Hangzhou Medical Disciplines of Excellence & Key Project for Hangzhou Medical Disciplines; and the Special Foundation for Brain Research from Science and Technology Program of Guangdong (2018B030334001).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2022.969575/full#supplementary-material
References
Andlauer, T. F. M., Guzman-Parra, J., Streit, F., Strohmaier, J., Gonzalez, M. J., Gil Flores, S., et al. (2021). Bipolar multiplex families have an increased burden of common risk variants for psychiatric disorders. Mol. Psychiatry 26 (4), 1286–1298. doi:10.1038/s41380-019-0558-2
Ashitani, J., Matsumoto, N., and Nakazato, M. (2009). Effect of octanoic acid-rich formula on plasma ghrelin levels in cachectic patients with chronic respiratory disease. Nutr. J. 8, 25. doi:10.1186/1475-2891-8-25
Barch, D. M., and Sheffield, J. M. (2014). Cognitive impairments in psychotic disorders: common mechanisms and measurement. World Psychiatry 13 (3), 224–232. doi:10.1002/wps.20145
Basser, P. J. (1995). Inferring microstructural features and the physiological state of tissues from diffusion-weighted images. NMR Biomed. 8 (7-8), 333–344. doi:10.1002/nbm.1940080707
Basser, P. J., and Pierpaoli, C. (1996). Microstructural and physiological features of tissues elucidated by quantitative-diffusion-tensor MRI. J. Magn. Reson. B 111 (3), 209–219. doi:10.1006/jmrb.1996.0086
Beyer, J. L., and Payne, M. E. (2016). Nutrition and bipolar depression. Psychiatr. Clin. North Am. 39 (1), 75–86. doi:10.1016/j.psc.2015.10.003
Bot, M., Milaneschi, Y., Al-Shehri, T., Amin, N., Garmaeva, S., Onderwater, G. L. J., et al. (2020). Metabolomics profile in depression: A pooled analysis of 230 metabolic markers in 5283 cases with depression and 10, 145 controls. Biol. Psychiatry 87 (5), 409–418. doi:10.1016/j.biopsych.2019.08.016
Brunkhorst-Kanaan, N., Klatt-Schreiner, K., Hackel, J., Schroter, K., Trautmann, S., Hahnefeld, L., et al. (2019). Targeted lipidomics reveal derangement of ceramides in major depression and bipolar disorder. Metabolism. 95, 65–76. doi:10.1016/j.metabol.2019.04.002
Charney, A. W., Ruderfer, D. M., Stahl, E. A., Moran, J. L., Chambert, K., Belliveau, R. A., et al. (2017). Evidence for genetic heterogeneity between clinical subtypes of bipolar disorder. Transl. Psychiatry 7 (1), e993. doi:10.1038/tp.2016.242
Cobley, J. N., Fiorello, M. L., and Bailey, D. M. (2018). 13 reasons why the brain is susceptible to oxidative stress. Redox Biol. 15, 490–503. doi:10.1016/j.redox.2018.01.008
Coles, B., Bloodsworth, A., Eiserich, J. P., Coffey, M. J., McLoughlin, R. M., Giddings, J. C., et al. (2002). Nitrolinoleate inhibits platelet activation by attenuating calcium mobilization and inducing phosphorylation of vasodilator-stimulated phosphoprotein through elevation of cAMP. J. Biol. Chem. 277 (8), 5832–5840. doi:10.1074/jbc.M105209200
Craddock, N., and Mynors-Wallis, L. (2014). Psychiatric diagnosis: impersonal, imperfect and important. Br. J. Psychiatry 204 (2), 93–95. doi:10.1192/bjp.bp.113.133090
Cross-Disorder Group of the Psychiatric Genomics, C. (2013). Identification of risk loci with shared effects on five major psychiatric disorders: a genome-wide analysis. Lancet (London, Engl. 381 (9875), 1371–1379. doi:10.1016/S0140-6736(12)62129-1
Cross-Disorder Group of the Psychiatric Genomics, C., Lee, S. H., Ripke, S., Neale, B. M., Faraone, S. V., Purcell, S. M., et al. (2013). Genetic relationship between five psychiatric disorders estimated from genome-wide SNPs. Nat. Genet. 45 (9), 984–994. doi:10.1038/ng.2711
Cross-Disorder Group of the Psychiatric Genomics Consortium. Electronic address., and Cross-Disorder Group of the Psychiatric Genomics (2019). Genomic relationships, novel loci, and pleiotropic mechanisms across eight psychiatric disorders. Cell 179 (7), 1469–1482. e1411. doi:10.1016/j.cell.2019.11.020
Cui, G., Qing, Y., Hu, X., Wang, P., Sun, L., Yang, X., et al. (2020). Serum metabolomic profiling based on fourier transform-ion cyclotron resonance-mass spectrometry: do the dysfunctions of metabolic pathways reveal a universal risk of oxidative stress in schizophrenia? Antioxid. Redox Signal. 33 (10), 679–688. doi:10.1089/ars.2020.8141
Culpepper, L. (2014). Misdiagnosis of bipolar depression in primary care practices. J. Clin. Psychiatry 75 (3), e05. doi:10.4088/JCP.13019tx1c
Czobor, P., Jaeger, J., Berns, S. M., Gonzalez, C., and Loftus, S. (2007). Neuropsychological symptom dimensions in bipolar disorder and schizophrenia. Bipolar Disord. 9 (1-2), 71–92. doi:10.1111/j.1399-5618.2007.00428.x
Delporte, C. (2013). Structure and physiological actions of ghrelin. Sci. (Cairo) 2013, 518909. doi:10.1155/2013/518909
Dickens, A. M., Sen, P., Kempton, M. J., Barrantes-Vidal, N., Iyegbe, C., Nordentoft, M., et al. (2021). Dysregulated lipid metabolism precedes onset of psychosis. Biol. Psychiatry 89 (3), 288–297. doi:10.1016/j.biopsych.2020.07.012
Dipasquale, S., Pariante, C. M., Dazzan, P., Aguglia, E., McGuire, P., and Mondelli, V. (2013). The dietary pattern of patients with schizophrenia: a systematic review. J. Psychiatr. Res. 47 (2), 197–207. doi:10.1016/j.jpsychires.2012.10.005
Dunn, W. B., Broadhurst, D., Begley, P., Zelena, E., Francis-McIntyre, S., Anderson, N., et al. (2011). Procedures for large-scale metabolic profiling of serum and plasma using gas chromatography and liquid chromatography coupled to mass spectrometry. Nat. Protoc. 6 (7), 1060–1083. doi:10.1038/nprot.2011.335
Favre, P., Pauling, M., Stout, J., Hozer, F., Sarrazin, S., Abe, C., et al. (2019). Widespread white matter microstructural abnormalities in bipolar disorder: evidence from mega- and meta-analyses across 3033 individuals. Neuropsychopharmacology 44 (13), 2285–2293. doi:10.1038/s41386-019-0485-6
Fraguas, D., Diaz-Caneja, C. M., Ayora, M., Hernandez-Alvarez, F., Rodriguez-Quiroga, A., Recio, S., et al. (2019). Oxidative stress and inflammation in first-episode psychosis: a systematic review and meta-analysis. Schizophr. Bull. 45 (4), 742–751. doi:10.1093/schbul/sby125
Fritz, K., Russell, A. M. T., Allwang, C., Kuiper, S., Lampe, L., and Malhi, G. S. (2017). Is a delay in the diagnosis of bipolar disorder inevitable? Bipolar Disord. 19 (5), 396–400. doi:10.1111/bdi.12499
Gao, X., and Hannoush, R. N. (2014). Single-cell imaging of Wnt palmitoylation by the acyltransferase porcupine. Nat. Chem. Biol. 10 (1), 61–68. doi:10.1038/nchembio.1392
Godwin, D., Alpert, K. I., Wang, L., and Mamah, D. (2018). Regional cortical thinning in young adults with schizophrenia but not psychotic or non-psychotic bipolar I disorder. Int. J. Bipolar Disord. 6 (1), 16. doi:10.1186/s40345-018-0124-x
Gu, Z., Schlesner, M., and Hubschmann, D. (2021). cola: an R/Bioconductor package for consensus partitioning through a general framework. Nucleic Acids Res. 49 (3), e15. doi:10.1093/nar/gkaa1146
Heckers, S., and Kendler, K. S. (2020). The evolution of Kraepelin's nosological principles. World Psychiatry 19 (3), 381–388. doi:10.1002/wps.20774
Hoseth, E. Z., Krull, F., Dieset, I., Morch, R. H., Hope, S., Gardsjord, E. S., et al. (2018). Exploring the Wnt signaling pathway in schizophrenia and bipolar disorder. Transl. Psychiatry 8 (1), 55. doi:10.1038/s41398-018-0102-1
Huang, Y., Liu, Y., Wu, Y., Tang, Y., Zhang, M., Liu, S., et al. (2022). Patterns of convergence and divergence between bipolar disorder type I and type II: Evidence from integrative genomic analyses. Front. Cell Dev. Biol. 10, 956265. doi:10.3389/fcell.2022.956265
Hussain, G., Anwar, H., Rasul, A., Imran, A., Qasim, M., Zafar, S., et al. (2020). Lipids as biomarkers of brain disorders. Crit. Rev. Food Sci. Nutr. 60 (3), 351–374. doi:10.1080/10408398.2018.1529653
Iriondo, A., Garcia-Sebastian, M., Arrospide, A., Arriba, M., Aurtenetxe, S., Barandiaran, M., et al. (2021). Plasma lipids are associated with white matter microstructural changes and axonal degeneration. Brain Imaging Behav. 15 (2), 1043–1057. doi:10.1007/s11682-020-00311-9
Johnson, W. E., Li, C., and Rabinovic, A. (2007). Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics 8 (1), 118–127. doi:10.1093/biostatistics/kxj037
Kalyanaraman, B. (2004). Nitrated lipids: a class of cell-signaling molecules. Proc. Natl. Acad. Sci. U. S. A. 101 (32), 11527–11528. doi:10.1073/pnas.0404309101
Kawai, K., Nakashima, M., Kojima, M., Yamashita, S., Takakura, S., Shimizu, M., et al. (2017). Ghrelin activation and neuropeptide Y elevation in response to medium chain triglyceride administration in anorexia nervosa patients. Clin. Nutr. ESPEN 17, 100–104. doi:10.1016/j.clnesp.2016.10.001
Kelly, S., Jahanshad, N., Zalesky, A., Kochunov, P., Agartz, I., Alloza, C., et al. (2018). Widespread white matter microstructural differences in schizophrenia across 4322 individuals: Results from the ENIGMA schizophrenia DTI working group. Mol. Psychiatry 23 (5), 1261–1269. doi:10.1038/mp.2017.170
Kemp, S., Huffnagel, I. C., Linthorst, G. E., Wanders, R. J., and Engelen, M. (2016). Adrenoleukodystrophy - neuroendocrine pathogenesis and redefinition of natural history. Nat. Rev. Endocrinol. 12 (10), 606–615. doi:10.1038/nrendo.2016.90
Kitchin, W., Cohen-Cole, S. A., and Mickel, S. F. (1987). Adrenoleukodystrophy: frequency of presentation as a psychiatric disorder. Biol. Psychiatry 22 (11), 1375–1387. doi:10.1016/0006-3223(87)90072-2
Kojima, M., Hosoda, H., Date, Y., Nakazato, M., Matsuo, H., and Kangawa, K. (1999). Ghrelin is a growth-hormone-releasing acylated peptide from stomach. Nature 402 (6762), 656–660. doi:10.1038/45230
Kojima, M., and Kangawa, K. (2002). Ghrelin, an orexigenic signaling molecule from the gastrointestinal tract. Curr. Opin. Pharmacol. 2 (6), 665–668. doi:10.1016/s1471-4892(02)00220-5
Koshiyama, D., Fukunaga, M., Okada, N., Morita, K., Nemoto, K., Usui, K., et al. (2020). White matter microstructural alterations across four major psychiatric disorders: mega-analysis study in 2937 individuals. Mol. Psychiatry 25 (4), 883–895. doi:10.1038/s41380-019-0553-7
Koshiyama, D., Fukunaga, M., Okada, N., Morita, K., Nemoto, K., Yamashita, F., et al. (2018). Role of frontal white matter and corpus callosum on social function in schizophrenia. Schizophr. Res. 202, 180–187. doi:10.1016/j.schres.2018.07.009
Koutoulogenis, G. S., and Kokotos, G. (2021). Nitro fatty acids (NO2-FAs): an emerging class of bioactive fatty acids. Molecules 26 (24), 7536. doi:10.3390/molecules26247536
Lankinen, M., Schwab, U., Seppanen-Laakso, T., Mattila, I., Juntunen, K., Mykkanen, H., et al. (2011). Metabolomic analysis of plasma metabolites that may mediate effects of rye bread on satiety and weight maintenance in postmenopausal women. J. Nutr. 141 (1), 31–36. doi:10.3945/jn.110.131656
Laye, S., Nadjar, A., Joffre, C., and Bazinet, R. P. (2018). Anti-inflammatory effects of omega-3 fatty acids in the brain: physiological mechanisms and relevance to pharmacology. Pharmacol. Rev. 70 (1), 12–38. doi:10.1124/pr.117.014092
Li, J., Tang, Y., Womer, F., Fan, G., Zhou, Q., Sun, W., et al. (2018). Two patterns of anterior insular cortex functional connectivity in bipolar disorder and schizophrenia. World J. Biol. Psychiatry 19 (3), S115-S123–S123. doi:10.1080/15622975.2016.1274051
Luan, H., Ji, F., Chen, Y., and Cai, Z. (2018). statTarget: a streamlined tool for signal drift correction and interpretations of quantitative mass spectrometry-based omics data. Anal. Chim. Acta 1036, 66–72. doi:10.1016/j.aca.2018.08.002
Magioncalda, P., and Martino, M. (2022). A unified model of the pathophysiology of bipolar disorder. Mol. Psychiatry 27 (1), 202–211. doi:10.1038/s41380-021-01091-4
Mamah, D., Alpert, K. I., Barch, D. M., Csernansky, J. G., and Wang, L. (2016). Subcortical neuromorphometry in schizophrenia spectrum and bipolar disorders. Neuroimage. Clin. 11, 276–286. doi:10.1016/j.nicl.2016.02.011
Meda, S. A., Ruano, G., Windemuth, A., O'Neil, K., Berwise, C., Dunn, S. M., et al. (2014). Multivariate analysis reveals genetic associations of the resting default mode network in psychotic bipolar disorder and schizophrenia. Proc. Natl. Acad. Sci. U. S. A. 111 (19), E2066–E2075. doi:10.1073/pnas.1313093111
Millan, M. J., Agid, Y., Brune, M., Bullmore, E. T., Carter, C. S., Clayton, N. S., et al. (2012). Cognitive dysfunction in psychiatric disorders: characteristics, causes and the quest for improved therapy. Nat. Rev. Drug Discov. 11 (2), 141–168. doi:10.1038/nrd3628
Nicholson, S. H., Suckling, C. J., and Iversen, L. L. (1979). GABA analogues: conformational analysis of effects on [3H]GABA binding to postsynaptic receptors in human cerebellum. J. Neurochem. 32 (1), 249–252. doi:10.1111/j.1471-4159.1979.tb04539.x
Oresic, M., Tang, J., Seppanen-Laakso, T., Mattila, I., Saarni, S. E., Saarni, S. I., et al. (2011). Metabolome in schizophrenia and other psychotic disorders: a general population-based study. Genome Med. 3 (3), 19. doi:10.1186/gm233
Paans, N. P. G., Bot, M., van Strien, T., Brouwer, I. A., Visser, M., and Penninx, B. (2018). Eating styles in major depressive disorder: Results from a large-scale study. J. Psychiatr. Res. 97, 38–46. doi:10.1016/j.jpsychires.2017.11.003
Pang, Z., Chong, J., Li, S., and Xia, J. (2020). MetaboAnalystR 3.0: toward an optimized workflow for global metabolomics. Metabolites 10 (5), E186. doi:10.3390/metabo10050186
Patti, G. J., Yanes, O., and Siuzdak, G. (2012). Innovation: Metabolomics: the apogee of the omics trilogy. Nat. Rev. Mol. Cell Biol. 13 (4), 263–269. doi:10.1038/nrm3314
Paul, L. K., Brown, W. S., Adolphs, R., Tyszka, J. M., Richards, L. J., Mukherjee, P., et al. (2007). Agenesis of the corpus callosum: genetic, developmental and functional aspects of connectivity. Nat. Rev. Neurosci. 8 (4), 287–299. doi:10.1038/nrn2107
Qian, G., Xu, L., Qin, J., Huang, H., Zhu, L., Tang, Y., et al. (2021). Leukocyte proteomics coupled with serum metabolomics identifies novel biomarkers and abnormal amino acid metabolism in Kawasaki disease. J. Proteomics 239, 104183. doi:10.1016/j.jprot.2021.104183
Reichenberg, A., Harvey, P. D., Bowie, C. R., Mojtabai, R., Rabinowitz, J., Heaton, R. K., et al. (2009). Neuropsychological function and dysfunction in schizophrenia and psychotic affective disorders. Schizophr. Bull. 35 (5), 1022–1029. doi:10.1093/schbul/sbn044
Reilly, J. L., and Sweeney, J. A. (2014). Generalized and specific neurocognitive deficits in psychotic disorders: utility for evaluating pharmacological treatment effects and as intermediate phenotypes for gene discovery. Schizophr. Bull. 40 (3), 516–522. doi:10.1093/schbul/sbu013
Rios-Esteves, J., Haugen, B., and Resh, M. D. (2014). Identification of key residues and regions important for porcupine-mediated wnt acylation. J. Biol. Chem. 289 (24), 17009–17019. doi:10.1074/jbc.M114.561209
Rohart, F., Gautier, B., Singh, A., and Le Cao, K. A. (2017). mixOmics: An R package for 'omics feature selection and multiple data integration. PLoS Comput. Biol. 13 (11), e1005752. doi:10.1371/journal.pcbi.1005752
Ruderfer, D. M., Fanous, A. H., Ripke, S., McQuillin, A., Amdur, R. L., Schizophrenia Working Group of the Psychiatric Genomics, C., et al. (2014). Schizophrenia working group of the psychiatric genomics polygenic dissection of diagnosis and clinical dimensions of bipolar disorder and schizophrenia. Mol. Psychiatry 19 (9), 1017–1024. doi:10.1038/mp.2013.138
Salim, S. (2017). Oxidative stress and the central nervous system. J. Pharmacol. Exp. Ther. 360 (1), 201–205. doi:10.1124/jpet.116.237503
Sankaranarayanan, A., Johnson, K., Mammen, S. J., Wilding, H. E., Vasani, D., Murali, V., et al. (2021). Disordered eating among people with schizophrenia spectrum disorders: a systematic review. Nutrients 13 (11), 3820. doi:10.3390/nu13113820
Scadding, J. G. (1996). Essentialism and nominalism in medicine: logic of diagnosis in disease terminology. Lancet 348 (9027), 594–596. doi:10.1016/s0140-6736(96)02049-1
Schopfer, F. J., Lin, Y., Baker, P. R., Cui, T., Garcia-Barrio, M., Zhang, J., et al. (2005). Nitrolinoleic acid: an endogenous peroxisome proliferator-activated receptor gamma ligand. Proc. Natl. Acad. Sci. U. S. A. 102 (7), 2340–2345. doi:10.1073/pnas.0408384102
Silverman, R. B., Andruszkiewicz, R., Nanavati, S. M., Taylor, C. P., and Vartanian, M. G. (1991). 3-Alkyl-4-aminobutyric acids: the first class of anticonvulsant agents that activates L-glutamic acid decarboxylase. J. Med. Chem. 34 (7), 2295–2298. doi:10.1021/jm00111a053
Slominski, A. T., Kim, T. K., Chen, J., Nguyen, M. N., Li, W., Yates, C. R., et al. (2012). Cytochrome P450scc-dependent metabolism of 7-dehydrocholesterol in placenta and epidermal keratinocytes. Int. J. Biochem. Cell Biol. 44 (11), 2003–2018. doi:10.1016/j.biocel.2012.07.027
Song, S. K., Sun, S. W., Ramsbottom, M. J., Chang, C., Russell, J., and Cross, A. H. (2002). Dysmyelination revealed through MRI as increased radial (but unchanged axial) diffusion of water. Neuroimage 17 (3), 1429–1436. doi:10.1006/nimg.2002.1267
Song, S. K., Yoshino, J., Le, T. Q., Lin, S. J., Sun, S. W., Cross, A. H., et al. (2005). Demyelination increases radial diffusivity in corpus callosum of mouse brain. Neuroimage 26 (1), 132–140. doi:10.1016/j.neuroimage.2005.01.028
Stogios, N., Smith, E., Asgariroozbehani, R., Hamel, L., Gdanski, A., Selby, P., et al. (2020). Exploring patterns of disturbed eating in psychosis: a scoping review. Nutrients 12 (12), E3883. doi:10.3390/nu12123883
Tabares-Seisdedos, R., and Rubenstein, J. L. (2009). Chromosome 8p as a potential hub for developmental neuropsychiatric disorders: Implications for schizophrenia, autism and cancer. Mol. Psychiatry 14 (6), 563–589. doi:10.1038/mp.2009.2
Tamminga, C. A., Pearlson, G., Keshavan, M., Sweeney, J., Clementz, B., and Thaker, G. (2014). Bipolar and schizophrenia network for intermediate phenotypes: outcomes across the psychosis continuum. Schizophr. Bull. 40 (2), S131–S137. doi:10.1093/schbul/sbt179
Taylor, C. P., Vartanian, M. G., Andruszkiewicz, R., and Silverman, R. B. (1992). 3-alkyl GABA and 3-alkylglutamic acid analogues: two new classes of anticonvulsant agents. Epilepsy Res. 11 (2), 103–110. doi:10.1016/0920-1211(92)90044-t
Torres, V. I., Godoy, J. A., and Inestrosa, N. C. (2019). Modulating Wnt signaling at the root: porcupine and Wnt acylation. Pharmacol. Ther. 198, 34–45. doi:10.1016/j.pharmthera.2019.02.009
van Velzen, L. S., Kelly, S., Isaev, D., Aleman, A., Aftanas, L. I., Bauer, J., et al. (2020). White matter disturbances in major depressive disorder: a coordinated analysis across 20 international cohorts in the ENIGMA MDD working group. Mol. Psychiatry 25 (7), 1511–1525. doi:10.1038/s41380-019-0477-2
Vancampfort, D., Firth, J., Schuch, F. B., Rosenbaum, S., Mugisha, J., Hallgren, M., et al. (2017). Sedentary behavior and physical activity levels in people with schizophrenia, bipolar disorder and major depressive disorder: a global systematic review and meta-analysis. World Psychiatry 16 (3), 308–315. doi:10.1002/wps.20458
Vigo, D., Thornicroft, G., and Atun, R. (2016). Estimating the true global burden of mental illness. Lancet. Psychiatry 3 (2), 171–178. doi:10.1016/S2215-0366(15)00505-2
Wei, C. Y., Zhu, D., Ryback, B., and Nyström, L. (2021). Evaluation of the oxidative stability of wholegrain flours using nontargeted ultraperformance liquid chromatography coupled with electrospray ionization quadrupole time-of-flight (UPLC-ESI-QTOF) metabolomics. ACS Food Sci. Technol. 1 (7), 1151–1159. doi:10.1021/acsfoodscitech.0c00161
Willert, K., Brown, J. D., Danenberg, E., Duncan, A. W., Weissman, I. L., Reya, T., et al. (2003). Wnt proteins are lipid-modified and can act as stem cell growth factors. Nature 423 (6938), 448–452. doi:10.1038/nature01611
Wright, M. M., Schopfer, F. J., Baker, P. R., Vidyasagar, V., Powell, P., Chumley, P., et al. (2006). Fatty acid transduction of nitric oxide signaling: nitrolinoleic acid potently activates endothelial heme oxygenase 1 expression. Proc. Natl. Acad. Sci. U. S. A. 103 (11), 4299–4304. doi:10.1073/pnas.0506541103
Zhao, Y. Y., Wu, S. P., Liu, S., Zhang, Y., and Lin, R. C. (2014). Ultra-performance liquid chromatography-mass spectrometry as a sensitive and powerful technology in lipidomic applications. Chem. Biol. Interact. 220, 181–192. doi:10.1016/j.cbi.2014.06.029
Zhou, X., Liu, L., Lan, X., Cohen, D., Zhang, Y., Ravindran, A. V., et al. (2019). Polyunsaturated fatty acids metabolism, purine metabolism and inosine as potential independent diagnostic biomarkers for major depressive disorder in children and adolescents. Mol. Psychiatry 24 (10), 1478–1488. doi:10.1038/s41380-018-0047-z
Keywords: schizophrenia, bipolar disorder, major depressive disorder, lipidomics profile, brain whiter matter
Citation: Tao S, Zhang Y, Wang Q, Qiao C, Deng W, Liang S, Wei J, Wei W, Yu H, Li X, Li M, Guo W, Ma X, Zhao L and Li T (2022) Identifying transdiagnostic biological subtypes across schizophrenia, bipolar disorder, and major depressive disorder based on lipidomics profiles. Front. Cell Dev. Biol. 10:969575. doi: 10.3389/fcell.2022.969575
Received: 15 June 2022; Accepted: 01 August 2022;
Published: 05 September 2022.
Edited by:
Weihua Yue, Peking University Sixth Hospital, ChinaReviewed by:
Paola Magioncalda, Taipei Medical University, TaiwanJiansong Zhou, Second Xiangya Hospital, Central South University, China
Hao Yan, Peking University, China
Copyright © 2022 Tao, Zhang, Wang, Qiao, Deng, Liang, Wei, Wei, Yu, Li, Li, Guo, Ma, Zhao and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tao Li, litaozjusc@zju.edu.cn