Skip to main content

Proteomic signatures of schizophrenia-sourced iPSC-derived neural cells and brain organoids are similar to patients' postmortem brains



Schizophrenia is a complex and severe neuropsychiatric disorder, with a wide range of debilitating symptoms. Several aspects of its multifactorial complexity are still unknown, and some are accepted to be an early developmental deficiency with a more specifically neurodevelopmental origin. Understanding the timepoints of disturbances during neural cell differentiation processes could lead to an insight into the development of the disorder. In this context, human brain organoids and neural cells differentiated from patient-derived induced pluripotent stem cells are of great interest as a model to study the developmental origins of the disease.


Here we evaluated the differential expression of proteins of schizophrenia patient-derived neural progenitors (NPCs), early neurons, and brain organoids in comparison to healthy individuals. Using bottom-up shotgun proteomics with a label-free approach for quantitative analysis, we found multiple dysregulated proteins since NPCs, modified, and disrupted the 21DIV neuronal differentiation, and cerebral organoids. Our experimental methods have shown impairments in pathways never before found in patient-derived induced pluripotent stem cells studies, such as spliceosomes and amino acid metabolism; but also, those such as axonal guidance and synaptogenesis, in line with postmortem tissue studies of schizophrenia patients.


In conclusion, here we provide comprehensive, large-scale, protein-level data of different neural cell models that may uncover early events in brain development, underlying several of the mechanisms within the origins of schizophrenia.


Schizophrenia is one of the most prevalent neuropsychiatric disorders, though poorly understood regarding its molecular mechanisms. In general, initial symptoms are present at the beginning of adulthood; however, molecular roots of the disease have been linked to neurodevelopmental dysfunctions [1]. Several symptoms overlap with other disorders, such as major depression and bipolar disorders, making an early diagnosis and treatment strategies more difficult [2]. Both genetics and the environmental course on neurodevelopment have been associated with the disease onset. Thus, continued discoveries and characterization of new factors that contribute to disease onset and progression are key to understanding the complexity of the disorder [3].

Advances over the past few years are a result of employing a broader spectrum of tools to research postmortem brains and animal models, from brain imaging to cell-based studies, including induced pluripotent stem cells (iPSCs) [4, 5]. Furthermore, human iPSCs (hiPSCs) overcome the impracticality and poor accessibility of human brain cell types [6] and are providing the possibility to challenge and question neural cell fates, prior to or during the onset of a disease in question. Considerable progress has been made regarding the neural differentiation of human pluripotent stem cells into mature neurons and cerebral organoids [7]. Human neural progenitor cells (hNPC) form a useful cell system for high-throughput screening due to their homogeneity, along with low complexity and limited differentiation potential. In comparison, brain or cerebral organoids are complex, three-dimensional (3D) culture systems composed of multiple cell types that self-organize into various brain regions similar to those in vivo, including the cerebral cortex, ventral forebrain, midbrain-hindbrain boundary, and hippocampus [8,9,10]. The combination of different cell types in a complex 3D configuration can better simulate brain biology and function, allowing cerebral organoids to reproduce the function and architecture of the brain, especially regarding development and neuronal plasticity.

Through the combination of hiPSCs with neural organoid differentiation and label-free shotgun proteomics, we can expand our knowledge of expression dynamics, providing a deeper understanding of biological processes along the course of development, and simulating in vivo conditions [5, 11]. Therefore, here we evaluate schizophrenia using neural progenitor cells (NPCs), immature neurons, and brain organoids derived from schizophrenia patients and compared them to control subject-derived neural cells using shotgun label-free proteomics. Cells derived directly from schizophrenia patients offered an integrated view of protein expression during neurodevelopment, depicting compromised biochemical pathways. In addition, it established an in vitro schizophrenia platform with the potential to modify and manipulate compromised pathways using a plethora of chemical compounds. This preliminary global proteome analysis during cell development may contribute to the search for new therapies with a focus on the personalized treatment of the disease.

Materials and methods

Pluripotent stem cell culture

Human induced pluripotent stem cells (hiPSCs) were reprogrammed from three subjects diagnosed with the schizophrenia spectrum, available at Coriell [Cat# GM23760, RRID:CVCL_T804 (male, 26y, paranoid schizophrenia); Cat# GM23761B, RRID:CVCL_T807 (female, 27y, schizoaffective disorder); and Cat# GM23762B, RRID:CVCL_T808 (male, 23y, early onset schizophrenia)] [12]. Three control hiPSC cell lines were used, a Coriell Cat# GM23279A, RRID:CVCL_F178 (female, 36y, apparently healthy); and two reprogrammed at the D’Or Institute for Research and Education, CF1 (male, 37y, apparently healthy) and CF2 (male, 31y, apparently healthy) [13]. Organoids were additionally derived from human embryonic stem cell lines BR1 (Laboratory for Embryonic Stem Cell Research (LaNCE), University of São Paulo) [14] and H9 (RRID:CVCL_9773). Further information regarding those cell lines can be found in Additional file 1: Table S1. iPSCs and ESCs were cultured in mTeSR1 (Stemcell Technologies) or E8 medium (Thermo Scientific), on a Matrigel (BD Biosciences)-coated surface. Passaging of colonies was performed manually at 70% confluence and maintained at 37 °C in humidified air with 5% CO2. The establishment of hiPSCs and the derivation of NPCs, neurons, and organoids were carried out following international standards and with the approval of the local research ethics council (CAAE: 32385314.9.0000.5249).

Human neural progenitor cells and neuronal differentiation

To induce the differentiation of hiPSCs into neural progenitors, we used a previously described protocol [15]. In brief, the hiPSC cultures from six individuals (three control and three schizophrenia patients), at 70% confluence, were induced into a neural lineage in a defined adherent culture by retinoic acid and basic fibroblast growth factor (bFGF). After 18 days, neural tube-like structures can be collected and plated on dishes coated with 10 μg/mL of poly-L-ornithine and 2.5 μg/mL of laminin (Thermo Fisher Scientific), in N2B27 medium (DMEM-F12 with 1 × N2 and 1 × B27 supplements, 1% penicillin/streptomycin) supplemented with 25 ng/mL bFGF and 20 ng/mL EGF (Thermo Fisher Scientific). Human neural progenitor cells (hNPCs) migrate from the neural tube-like structures and were thereupon tested for the expression of neural markers before expansion. Cell expansion was done in N2B27 medium supplemented with EGF and bFGF, replaced every other day. hNPCs were expanded for no more than 5 passages. Characterization of this culture has been previously published [16].

For neuronal differentiation, hNPCs of each lineage were plated on poly-L-ornithine- and laminin-coated dishes. To trigger differentiation, 24 h after plating bFGF and EGF were removed from the N2B27 medium. Neuronal cultures were differentiated for 21 days in vitro (DIV). The whole medium was replaced every 5 days. Both hNPCs and neuronal cells were incubated at 37ºC and 5% CO2. After the period described, cells were collected for proteomics or immunocytochemistry.

Differentiation into brain organoids

For differentiation of schizophrenia and control human iPSCs and ESCs, cells were cultured in mTeSR1 medium (Stemcell Technologies) on Matrigel (BD Biosciences)-coated cell culture dishes. Differentiation of pluripotent stem cells into cerebral organoids was based on a previously described protocol [17, 18]. In summary, hiPSCs were passaged into single cells with Accutase (Merck Millipore) and inoculated in a spinner flask containing mTeSR1 medium supplemented with 10 μM Y-27632 (Rho-associated protein kinases inhibitor, iRock) (Merck Millipore) under constant rotation (40 rpm). The following day, the medium was replaced to initiate embryoid body formation. One week later, neural induction medium [DMEM/F12 1:1 supplemented with N2 (1x) supplement, 2 mM glutamax, 1% MEM-NEAA (Thermo Scientific), and heparin (1 μg/mL, Sigma)] was added. After 4 days, cellular aggregates were covered in Matrigel and cultured in differentiation medium [DMEM/F12:Neurobasal (1:1), supplemented with N2 (0.5x) and B27 minus vitamin A (1x) supplements, 2 mM glutamax, 0.5% MEM-NEAA, 0.2 μM 2-mercaptoethanol (Thermo Scientific), and 2.5 μg/mL insulin (Sigma)] for another 4 additional days. Subsequently, the medium was replaced with a neuronal differentiation medium, which comprises the same formulation, except for the use of 1 × B27 containing vitamin A (Thermo Scientific). This final medium was then changed every week for the differentiation of cerebral organoids, grown for 45 days, and collected for analyses. At 45-days, organoids present a well-organized structure, with both radial glia and a neuronal network [11, 17, 19], which is of interest in modeling initial developmental cues in schizophrenia.


Cells were grown on coverslips for 5 days (for NPCs) or directly differentiated into neurons (for 21 DIV), washed with PBS to remove the medium, and fixed in 4% paraformaldehyde (PFA). Cerebral organoids were collected after 45 days of differentiation and immediately fixed in 4% PFA, followed by incubation in sucrose solutions over an increasing gradient (10, 20, and 30%) in phosphate-buffered saline (PBS). Subsequently, the organoids were embedded in an optimal cutting temperature compound (OCT) and frozen in liquid nitrogen. The organoids were sectioned with a cryostat (Leica) into 20 μm-thick sections. After being fixed, cells and organoids were washed with PBS, permeabilized in a 0.3% Triton-X solution, and blocked in a 3% bovine serum albumin (BSA) solution before immunolabeling. Immunofluorescence was performed using the primary antibodies: anti-MAP2 (Sigma-Aldrich Cat# M1406, RRID:AB_477171), anti-class III β-tubulin (Millipore Cat# MAB1637, RRID:AB_2210524), anti-nestin (Millipore Cat# MAB5326, RRID:AB_2251134), anti-PAX6 (Santa Cruz Biotechnology Cat# sc-11357, RRID:AB_2159706), anti-SOX1 (Millipore Cat# AB15766, RRID:AB_870981), and anti-SOX2 (Millipore Cat# AB5603, RRID:AB_2286686). Secondary antibodies used were AlexaFluor 488 goat anti-rabbit (Molecular Probes Cat# A-11008, RRID:AB_143165), AlexaFluor 488 goat anti-mouse (Thermo Fisher Scientific Cat# A-11001, RRID:AB_2534069), and AlexaFluor 594 goat anti-mouse (Molecular Probes Cat# A-11032, RRID:AB_2534091). DAPI was used for nucleus staining. Images were acquired using a Leica TCS SP8 confocal microscope.

Sample preparation and processing

From each schizophrenia patient or control cell line, three independent NPCs (5 DIV) and neurons (21 DIV) samples, following protocols described earlier, were harvested, and pelleted in PBS. Similarly, five to six cerebral organoids were pooled from the spinner flask and pelleted in PBS to provide population variability within each patient or cell line. At least two batches of spinners were polled per cell line. Cell lines were then independently sampled and analyzed. Pellets of cells or organoids were removed from PBS and homogenized in lysis buffer [7 M urea, 2 M thiourea, 1% CHAPS, 70 mM DTT, and EDTA-free complete protease inhibitor cocktail (Roche)]. Sample lysates were kept on ice for approximately 20 min and centrifuged at 10,000 × g for 10 min at 4 °C; supernatants were collected, and protein content was quantified using a Qubit® 3.0 Fluorometer (Thermo Fisher Scientific). Each sample (100 µg) was subjected to a short SDS-PAGE run and a sequence of reduction, alkylation, and overnight trypsin digestion (1:50 w/w trypsin:total protein) at 37 °C. Peptides collected from this digestion were dried in a SpeedVac (Thermo Fisher Scientific) and stored at −80 °C until quantitative and qualitative shotgun mass spectrometry analyses.

Liquid chromatography-mass spectrometry

Peptides were separated by a two-dimensional nanoAcquity UPLC M-Class System (Waters Corporation, Milford, MA) liquid chromatographer coupled to a Synapt G2-Si mass spectrometer (Waters Corporation). Organoid samples were separated by two-dimension chromatography, for the first-dimension reverse-phase chromatography, peptides (5 µg) were loaded onto an M-Class BEH C18 Column (130 Å, 5 µm, 300 µm X 50 mm, Waters Corporation). Peptide fractionation was performed by increasing steps of acetonitrile concentration (13%, 18%, and 50% ACN). Peptide loads were directed to a second-dimensional separation on an HSS T3 Column (100 Å, 1.8 µm, 75 µm X 150 mm, Waters Corporation, Milford, MA), with a binary gradient of 7% to 40% ACN (v/v) over 54 min at a flow rate of 0.4 µL/min. The fractionation for NPC and 21 DIV neuronal samples was performed in one dimension. Peptide loads were directed to separation on an HSS T3 Column (100 Å, 1.8 µm, 75 µm × 150 mm, Waters Corporation, Milford, MA), with a binary gradient of 7% to 40% ACN (v/v) over 54 min at a flow rate of 0.4 µL/min. In all cases, peptides entered the mass spectrometer using nano-electrospray ionization in positive ion mode, nanoESI ( +). The MS analyses were performed using data-independent acquisition (DIA) enhanced with ion mobility separation (HDMSE). [Glu1]-Fibrinopeptide B human was used as the lock mass compound, which was sampled every 30 s. Each cell line was considered a biological sample, and at least duplicates (for organoids) or triplicates (for NPCs and 21-DIV neurons) were run. The LC–MS/MS method used was based on a previously described protocol [20].

Database search and quantitation

Raw data were aligned and processed in Progenesis® QI for proteomics version 3.0 (Waters). Protein identification and quantification were performed using the ion accounting algorithm with default parameters and searching against the Homo sapiens database—revised (Uniprot, version 2017/10). For protein identification, the following parameters were set: up to two missed cleavages for trypsin digestion; variable modification by oxidation (M) and fixed modification by carbamidomethyl (C), and False Discovery Rate (FDR) less than 1%. The quantitative analysis was carried out on the log2-values of the intensities after normalization by a software-calculated global scaling factor, and only proteins present in at least two out of three biological samples of both control and schizophrenia samples were considered as discovered. All proteins with an ANOVA p-value < 0.05, between the schizophrenia and control samples, were considered differentially expressed and selected for further analysis (Additional file 2: Table S2).

In silico analysis

Gene ontology and pathway enrichment were analyzed with DAVID (the Database for Annotation, Visualization and Integrated Discovery) [21, 22] using whole-genome background; and Metascape [23], using default settings (whole proteome) and the following databases: KEGG Pathway, GO Biological Processes, Reactome and CORUM. Protein networks and canonical pathways associated with differentially expressed proteins were identified using Ingenuity Pathway Analysis software (IPA, Ingenuity Systems, Qiagen, Redwood, CA, USA; The significance of biological functions was calculated using Fisher's exact test. Multiple correlation hypotheses were calculated with the Benjamini-Hochberg (B-H) approach using a 1% FDR threshold; the significance levels of the IPA tests were expressed as p-values. Differentially regulated proteins in schizophrenia postmortem brain tissue were selected from the literature by searching the Pubmed database ( using the following keywords: proteomics AND schizophrenia AND postmortem brain. Studies included were those with a list containing differentially regulated proteins, from several brain areas, such as the corpus callosum, anterior temporal lobe, anterior cingulate cortex, dorsolateral prefrontal cortex, mediodorsal thalamus, insular cortex, frontal cortex, Wernicke’s area, cerebellum, posterior cingulate cortex, caudate nucleus, and the hippocampus. The whole set of proteomic datasets was then compiled and is detailed in Additional file 3: Table S3. The compiled list of regulated proteins in postmortem brain areas of schizophrenia patients was then compared to the regulated proteomic data of PSC-derived schizophrenia versus control cerebral organoids, neurons, and NPCs. Comparative analyses were then performed in the Python programming language (v. 3.7.3). The overlap between proteins identified in our data and previous postmortem brain studies was visualized with a Chord diagram [24]. The over-representation analysis of brain regions and KEGG pathways was carried out using ClusterProfiler [25] in the R programming environment (v. 4.0).


Neurodevelopmental aspects of schizophrenia were investigated with the proteomic profiles of three models of neural differentiation: neural progenitors, immature neurons, and brain organoids. hiPSC from three schizophrenia patients and three paired controls were differentiated into neural cells and submitted to proteomic analyses, in a primary study to undercover dysregulated proteins, as presented in the workflow (Fig. 1). Each of the hiPSC samples from patients and controls (Fig. 1A) was differentiated into neural progenitor cells (NPCs) (Fig. 1B), young neurons with 21 days in vitro (21 DIV) (Fig. 1C), and brain organoids (cultured for 45 days) (Fig. 1D). Whole-cell proteomic profiles were generated using label-free quantitative proteomics (Fig. 1E) to uncover some of the molecular mechanisms of schizophrenia during neurodevelopment.

Fig. 1
figure 1

Schematic organization of the experimental design and analysis workflow. A Human iPSCs from schizophrenia patients and controls were differentiated into neural cell types: progenitors, neurons, and organoids. Representative photomicrographs of B neural progenitor cells (NPC), C neurons at 21 days in vitro, and D 45-day cerebral organoids. E Proteomics workflow processing. Label-free sample preparation (protein extraction and peptide digestion) followed by 2D-UPLC fractionation and on-line detection using HDMSE high-resolution MS/MS acquisition. Peptides and proteins were identified and quantified before functional annotation and other in silico analyses. F Differentially regulated proteins in cerebral organoids of schizophrenia patients were compared to available proteomics data from post-mortem brains of schizophrenia patients. Scale bars shown are B 400 µm, C 200 µm and D 1000 µm

In addition, the proteomic data from organoids, 21 DIV neurons, and NPCs data (schizophrenia compared to controls, p < 0.05, as in Additional file 2: Table S2) were later further compared to a compiled list of publicly available proteomics datasets (Additional file 3: Table S3) from postmortem brain tissue of schizophrenia patients (Fig. 1F).

Establishing in vitro models of schizophrenia

We first established cultures with progenitors and neurons differentiated from hiPSC cells of schizophrenia patients and controls, using a methodology we have previously characterized to provide mainly forebrain neural progenitors [26]. While NPCs were kept under a proliferative medium, neuronal differentiation was performed by removing FGF and EGF from the growth medium, switching the proliferation program to a neuronal one. At the end of each culture period, NPCs in both control and schizophrenia-derived cells stained positive for neural progenitor markers such as nestin and SOX2 (SRY-Box 2) (Fig. 2A). Differentiation of progenitors into young neurons was also induced in both schizophrenia and control cells. After 21 days of differentiation in vitro (21 DIV), morphological changes were observed in both control and schizophrenia cells, showing small cell bodies and a large number of cells with neurite elongation (Fig. 2A). Neurons at 21 DIV continued to stain for nestin and SOX2, with few cells staining for PAX6 (paired box 6), while several stained for β3-tubulin (TUBB3) (Fig. 2A), and fewer showed matured MAP2 staining (microtubule-associated protein 2) (Additional file 4: Figure S1A). Proteomic analysis also indicated the expression of nestin and TUBB3; in progenitors, the transcription factor PAX6 was observed, though it was not detected in neurons (Additional file 2: Table S2). Gene expression analyses (Additional file 4: Figure S1B) have shown that 21 DIV neurons express neuronal associated cytoskeleton, such as MAP2 and TUBB3, and synapses-related mRNA expression such as PSD-95 (postsynaptic density 95) and increased synaptophysin (SYP) in comparison to progenitors. In addition, as the differentiation protocol was not directed to particular brain regions, neuronal differentiation is heterogeneous. We observed the expression of genes related to the GABAergic systems, as shown by glutamic acid decarboxylase (GAD65), as well as genes involved in the differentiation of midbrain dopaminergic neurons, such as NURR1 (nuclear receptor-related 1 protein, or NR4A2) and TH (tyrosine hydroxylase), to be increasing in 21 DIV neurons in comparison to progenitors (Additional file 4: Figure S1B). Instead, the glutamatergic system, such as glutaminase (GLS), glutamate transporter (VGLUT1), and NMDA receptor (NMDAR), are present yet more prominent in progenitors. Glutamate is responsible for regulating neurogenesis, and its activity in progenitors influences neurodevelopment, secreting other neurotrophic factors. Our results point to a neuronal diversity toward increasing dopaminergic and GABAergic developmental cues in the differentiation protocol. In the human fetal cortex, NMDAR regulation in progenitors might directly increase the number of TH-positive neurons, regulating cell proliferation and the production of dopaminergic cells [27, 28].

Fig. 2
figure 2

Morphological characteristics between controls and schizophrenia patients are similar in progenitors (NPCs), neurons (21 DIV), and cerebral organoids. A Immunocytochemistry characterization of NPCs and 21 DIV neurons. Representative micrographs of control NPCs and neurons in upper panels (CTR), and schizophrenia NPCs and neurons in the lower panels (SCZ). Showing SOX2 and nestin; PAX6 and TUBB3. Scale bars = 100 µm. B Immunocytochemistry of controls and SCZ-derived cerebral organoids showing ventricle-like morphology and SOX2 and MAP2 staining, as well as TUBB3. Scale bars = 100 µm

To evaluate a more complex model of cellular organization in vitro, we generated cerebral organoids at 45 days, a point at which a neuronal network has already been initially formed [11, 17]. The whole-cerebral organoid protocol used has no cues for inducing a specific region, and the self-organized cytoarchitecture in both controls and schizophrenia-derived organoids was similar. The organoids from control and schizophrenia showed no major morphological differences, and 45-day-old cerebral organoids developed one to several putative ventricles, with young neurons identified by TUBB3-positive cells found at the subventricular zone in both control- and schizophrenia-derived organoids (Fig. 2B). Control organoids showed organized ventricle zones with SOX2-expressing progenitor cells and MAP2-expressing neuronal cells at the cortical layer, while schizophrenia patients had a larger area covered by SOX2 progenitors, also showing MAP2 distribution throughout the organoid (Fig. 2C). In summary, these morphological aspects and protein markers confirm the models as suitable for comparison of schizophrenia and control differences during neurodevelopment.

The proteomics of schizophrenia-derived neural cells

Proteomic analyses were performed at three different developmental models of hiPSC-derived neural cells: neural progenitor cells (NPCs), neurons at 21 days in vitro (21 DIV), and cerebral organoids at 45 days (Additional file 2: Table S2). The NPC dataset yielded a total of 1949 quantified proteins, of which 364 proteins were found to be dysregulated (p < 0.05) between schizophrenia-patients NPCs and controls. Of those proteins, 84% (306) were downregulated (Additional file 2: Table S2). In the 21 DIV neuronal model, 1833 proteins were quantified, 264 of which were considered dysregulated (p < 0.05) between schizophrenia and controls. Of these proteins, 70% (185) were downregulated. A slightly larger number of proteins was quantified from the 45-day cerebral organoids, which yielded 2177 quantified proteins. Between schizophrenia and control groups, 535 proteins were considered dysregulated (p < 0.05), 59% (317) of which were downregulated.

When comparing dysregulated proteins of schizophrenia versus control between the three cellular models, we found that some dysregulated proteins were common among NPC, neurons, and cerebral organoids (Fig. 3A, represented by purple lines), which points to several disrupted pathways at all developmental models (Fig. 3A, represented by blue lines). Twelve dysregulated proteins were commonly found in all models comparing schizophrenia with controls, most of them downregulated in the three models, which included i.e., 14–3-3 protein eta (YWHAH) and spectrin beta chain (SPTBN2); and others, such as poly [ADP-ribose] polymerase 1 (PARP1), were shown upregulated in organoids (Additional file 4: Figure S2).

Fig. 3
figure 3

Differential regulation of proteins in schizophrenia patients versus control cells from human PSC-derived organoids, neurons, and NPCs. A Chord diagram overviewing individual proteins (purple lines) and pathways (blue lines) shared among comparisons. The numbers of proteins (total and overlapped) are detailed. B Functional network tree representation of enriched pathways, processes, and GO terms. Terms were grouped by similarity metric (kappa scores) and the most statistically significant term within a cluster was set to represent the cluster, according to Metascape analyses. C Canonical pathway analyses were performed using the Ingenuity pathway analyses dataset. The blue plot represents pathway enrichment (Log2 adjusted p-value), and the purple/green plot represents z-scores prediction of activation/inhibition of those pathways. Purple represents activated scores, and green represent inhibited scores

We have selected high abundant proteins with varied unique peptides from the dysregulated list (YWHAH, MAP4, THY1, and CTNNB1), independently of being up- or downregulated, for validation by a multiple reaction monitoring method. Both analyses (shotgun and MRM) are different, and protein quantification can be distinct due to specific peptides chosen for intensity quantification, and post-translational modifications content (which were not included in these searches). Therefore, we have validated the presence of the protein by MRM qualitatively, yet not all were as well quantitatively equal (Additional file 4: Table S4).

Combining these observations of dysregulated proteins in schizophrenia (versus control), constituent pathways can be observed in more detail (Fig. 3B). Organoids and neurons have deregulation in mRNA processing, including splicing. This is followed by disruptions in protein synthesis and folding at all cell models analyzed, which are presumed to have downstream consequences on axon guidance, exocytosis, cell–cell adhesion, and cytoskeleton organization, disrupting hemostasis, and cell maintenance.

Ephrin B and ephrin receptor signaling canonical pathways were dysregulated in schizophrenia-derived cerebral organoids, NPCs, and neurons in comparison with controls. Our data indicate a major downregulation of proteins involved in Ephrin signaling (-Log10 p-value– organoids: 6.53; NPCs: 2.60; neurons: 1.24) and further prediction analysis using Ingenuity Pathway Analysis offered some possible outcomes (Fig. 3C). Another family of proteins found enriched in all analyses was 14–3-3 signaling (-Log10 p-value– organoids: 4.96; NPCs: 2.97; neurons: 2.44), indicating inhibition of this pathway. Axon guidance (-Log10 p-value – organoids: 7.68; NPCs: 2.63; neurons: 1.25) and synaptogenesis (-Log10 p-value – organoids: 3.74; NPCs: 3.77; neurons: 3.54) were also found dysregulated in all schizophrenia neural cells (Fig. 5). All these pathways have been previously associated with schizophrenia [29], and they perform essential roles in brain development.

Individual analyses of organoids, NPCs, and neurons derived from patients with schizophrenia are represented in Fig. 4 (A, C, and E). Furthermore, the protein–protein interaction (PPI) analyses performed by Metascape show that schizophrenia-derived cerebral organoids presented dysregulated proteins that are involved in mRNA splicing and RNA metabolism (Fig. 4B). Schizophrenia-derived NPCs presented proteins related to translation and oxidative phosphorylation (Fig. 4D). Lastly, proteins seen in schizophrenia-derived neurons (21 DIV) are associated with amino acid metabolic processes and cytosolic tRNA aminoacylation (Fig. 4F).

Fig. 4
figure 4

Pathway regulation and protein–protein interaction (PPI) visualization in organoids, neurons, and NPCs of schizophrenia patients versus controls regulated proteins. KEGG pathways found dysregulated by the differentially expressed proteins (DEP) from schizophrenia patients versus controls in A organoids, C NPCs, and E 21 DIV neurons are individually shown. Bars represent the number of DEPs found in a particular pathway and the color represents the enrichment adjusted p-value. The protein–protein interaction in B organoids connects proteins found regulated in mRNA splicing and metabolism of RNA. PPI of proteins regulated in D NPCs are related to translation, axon guidance, and oxidative phosphorylation. And PPI of proteins regulated in F neurons are part of cellular amino acid metabolic process and cytosolic tRNA aminoacylation. Circles in red indicate upregulation of proteins and green indicate downregulation of proteins. Circle size reflects the interaction degree of each protein. Analyses were performed using Metascape resources

Similarities among hiPSC-derived neural cells, brain organoids, and schizophrenia postmortem brains

Several of the results obtained with hiPSC-derived neural cells and organoids were similar to those previously described in studies with postmortem brain tissue from schizophrenia patients. Observing such similarities, we questioned whether organoids derived from patient cells could be used as an in vitro model to study aspects of the pathophysiology of schizophrenia. Thus, we compiled a list of proteomics data containing several of the available studies in postmortem brain proteomics from the literature considering only the dysregulated proteins (in comparison to healthy individuals) of several brain regions (Additional file 3: Table S3). Then, we have compared that dysregulated list to the hiPSC-derived neural cells and organoids data from schizophrenia patients versus healthy controls of this work (p < 0.05, Additional file 2: Table S2). The signaling pathways and proteins that were found in common between hiPSC-derived neural cells/organoids and brain tissue are shown in Fig. 5. A total of 126 proteins were found commonly regulated when comparing organoids and postmortem brain tissue list, while 50 were found commonly regulated between NPCs and brain tissue, and 56 between neurons and brain tissue (Fig. 5A).

Fig. 5
figure 5

Proteomic from iPSC-derived neural cells and organoids from schizophrenia patients present similar patterns when compared with post-mortem brain tissue of patients with schizophrenia. A Chord-diagram (upper) overviewing proteins (purple lines) shared among comparisons; upset plot (lower) indicating the number of proteins shared for each comparison. Vertical bar sizes (intersection size) represent the number of proteins in common among the set groups, where the purple dots (and connecting line) are pointing to where this intersection (or overlap) occurs; and the horizontal bars represent the number of proteins differentially regulated (SCZ versus CTR) in each dataset (set size). B Canonical pathway enrichment analysis against the KEGG database (upper), showing pathways regulated by protein present in each dataset, followed by an upset plot indicating the connections of those pathways among organoids, 21 DIV neurons, NPCs, and the postmortem brain. Highlighted terms (bold red) represent novel pathways found in this study, not covered in postmortem brain tissues. Color shades (purple to blue) are representative of the p-value, while circle sizes represent the number of proteins within respective canonical pathway. Vertical bars represent the number of pathways regulated in each dataset (set size); horizontal bars (intersection size) represent the number of pathways in common among the set group; green dots (and connecting lines) are pointing to where this intersection occurs. C Comparative analysis of total organoids dataset with proteins from different regions of the postmortem brain. Bars represent the number of proteins, while they were colored according to their adjusted p-value of the analyses overlap. DLPFC dorsolateral prefrontal cortex, WA Wernicke’s area, ACC nterior cingulate cortex

There is a high intersection of proteins between organoids and postmortem brain analyses of patients with schizophrenia (Fig. 5A), which is followed by a similar outcome in terms of canonical pathways (Fig. 5B), with pathways such as focal adhesion, protein processing, and pathways on neurodegeneration, corresponding to some of the several representative pathways modified. ErbB signaling pathway, RNA degradation and protein export, are for instance canonical pathways not commonly overrepresented in the postmortem dataset. The intersection shows higher similarity between organoids and postmortem brain (in comparison to NPC and neurons), thus, we followed with a comparative analysis of brain regions. Using proteins identified in the organoids, we compared different postmortem brain regions' specific expression, as to identify which of those regions would have greater similarity to the organoid data. Comparing those proteins on the dataset and the organoids derived, we found the dorsolateral prefrontal cortex (DLPFC), followed by Wernicke’s area (WA) and the anterior cingulate cortex (ACC) (Fig. 5C), as those with the stronger correlation, and therefore similarity to the organoids.


Due to the difficulty in obtaining brain tissue from fetuses, along with ethical considerations, animal models were almost exclusively used to study and understand neurodevelopment. hiPSC cultures, especially 3D models such as organoids, hold great promise for the study of neurodevelopment and complex, multifactorial diseases since these cells preserve the genetic background of patients [30, 31]. Considering the cell microenvironment necessary for development, it is of particular interest to study patient-derived 3D brain organoids in the disease context. And it is important to state that those are experimental cell models, which although not representing the actual human brain, can provide valuable data to understand the early human neurodevelopmental paths disrupted in disease.

The neurodevelopmental hypothesis of schizophrenia is largely accepted and is based on structural, cognitive, and motor abnormalities that are identified in childhood and coincide with the development of the disease as adults (reviewed in [32, 33]). Furthermore, environmental factors such as hypoxia, viral infections, perinatal injury, and maternal drug abuse have been considered major risk factors for the development of this disease [34,35,36,37,38,39,40]. It is for this reason that this study investigated the proteins expressed in brain organoids (45 days), neural progenitor cells, and immature neurons (21 DIV) to understand the neurodevelopmental machinery at work in the cells of patients with schizophrenia.

Early dysregulation of axonal guidance and synaptogenesis signaling in schizophrenia at the protein level

Several of the main dysregulations observed here are related to the development of neurons and their connections, which will drive the formation of synapses in the brain in development. The ephrin receptor signaling pathway (Fig. 3C) is one of those found dysregulated in all the models of neural cells studied. Previous studies using genetic approaches, both in vivo and in vitro, have shown the role of this pathway in important processes of neural development such as axon guidance and dendritic spine formation [41,42,43,44,45]. This canonical pathway was also found to be dysregulated in the postmortem brains of patients with schizophrenia [46].

The decrease of constituent proteins of the ephrin pathway impairs the formation of dendritic spines [41, 47] and postmortem tissue studies show that patients with schizophrenia have a reduced number of dendritic spines [48]. Both in organoids and NPCs, the ephrin receptor signaling pathway was predicted to be inhibited (IPA canonical pathway analysis, Fig. 3C). In contrast, this pathway was predicted to be activated in neurons, suggesting that regulatory mechanisms are being triggered in this cell type to reverse the dysregulation in synaptic formation observed in neural progenitor cells (NPCs). The ephrin pathway is known to regulate NMDA receptor activity [47], which has already been widely associated with the pathophysiology of schizophrenia, and glutamatergic dysfunction is another hypothesis for the etiology of the disease [49]. NMDA hypofunction in neurodevelopment causes changes in the number of dendritic spines and the development of the neural circuit [50]. Recent GWAS studies [51] have pointed to several genes associated with neuronal function, which implicates modifications in cell differentiation, the transmission of synapsis, and synaptic organization being altered in patients with schizophrenia. Those associations were mostly found in genes expressed in both excitatory and inhibitory neurons. And, although pre- and post-synaptic proteins are not shown altered in our models, possibly due to the models being young and still immature, dysregulations, mostly downregulations occur in disease, and several proteins participate in pathways leading to future synapsis impairment.

Another signaling pathway found dysregulated in all three analyses was 14–3-3 (Fig. 2C). These proteins are abundant in the mammalian brain and comprise 1% of the cerebral proteome [52], playing roles in cell cycle regulation, apoptosis, differentiation, and migration, as well as more complex processes such as neurotransmission, neuroplasticity, and synaptogenesis (reviewed in [53]). Proteins of the 14–3-3 pathway have also been specifically associated with the pathophysiology of schizophrenia [46, 54]. Functional knockout of 14–3-3 family proteins lead to deficits in neuronal migration and structural organization of cells of the cortex and hippocampus [55, 56], as well as behavioral changes and cognitive deficits [57, 58]. The protein 14–3-3 eta (YWHAH), in particular, was found downregulated in NPCs, neurons, and organoids. It is a protein located in a chromosomal region (22q12), which is considered a high-risk region for developing schizophrenia; and RNA levels of this protein have been found downregulated in postmortem cerebellum [59,60,61]. In addition, this protein has been found altered in four distinct postmortem brain regions of patients with schizophrenia [62,63,64].

The actin cytoskeleton signaling pathway, which was found to be dysregulated in NPCs, neurons, and organoids, and predicted to be inhibited, plays a role in the movement of growth cones, maintenance of neural polarization, formation of dendritic spines, and formation/stabilization of the synapse, functions that are known to be affected in patients with schizophrenia (reviewed in [65]). The Rho GTPase family regulates actin in neuronal morphology [66]. In NPCs and organoids, RhoA signaling and Rho family GTPases pathway were predicted to be inhibited; and both transforming protein RhoA (RHOA) and actin-related protein 2 (ACRT2) were found downregulated in the schizophrenia-patients’ cells. However, in neurons, we have found those canonical pathways predicted to be more active, with Rho guanine nucleotide exchange factor 4 (ARHGEF4) and Rho GTPase-activating protein 45 (ARHGAP45) upregulated in these cells. As the Rho GTPase family plays a role in the regulation of both axon guidance and synaptogenesis (Additional file 4: Figure S3 and 4, respectively), the activation of these proteins in neurons is probably a cellular response to restore actin signaling and consequently proper synaptic formation, which is impaired in NPCs derivate of patients with schizophrenia [67]. Synaptic impairment was also observed in a transcriptome study of brain organoids [68], and similarly, a decrease of synaptic proteins and neural connection was observed in an in vivo proteomic study of neural progenitor cells from schizophrenia patients [4]. Overall, these results indicate a possible dysregulation in the formation and maturation of synapses during development, especially regarding axon guidance and synaptogenesis and it is in line with previous hiPSC and postmortem tissue studies.

Organoids and immature neurons exhibit differences in main biochemical pathways and protein–protein interactions when compared with previous hiPSC studies

Spliceosome-associated proteins were enriched in patient-derived organoids. The splicing process is coordinated by five multi-megadalton ribonucleoproteins (snRNPs) forming a complex. Malfunctions in the splicing complex can lead to several diseases, many observed in the nervous system [69], mostly due to neurotransmission dysregulation, especially regarding the GABAergic, dopaminergic, and glutamatergic systems [70, 71]. Moreover, genes related to brain development and maturation demonstrate abnormal regulation of alternative splicing in schizophrenia [72].

Some proteins of the snRNP complex, the core component of the splicing operation, such as heterogeneous nuclear ribonucleoproteins L (HNRNPL), F (HNRNPF), and A2/B1 (HNRNPA2), along with heat shock 70 kDa protein 4L (HSPA4L), were found upregulated in the schizophrenia-derived organoids. These proteins have been previously found in schizophrenia studies involving postmortem brains and in vitro techniques [73,74,75,76]. No previous study, observing hiPSCs derived cells from schizophrenia patients had indicated dysregulation in this pathway, placing focus on this target for future hiPSC studies.

PPI interaction analyses revealed a main dysregulation in oxidative phosphorylation (OXPHOS) of patient-derived neural progenitor cells (NPCs). Defects in this pathway during neurodevelopment have been associated with several psychiatric diseases [77]. OXPHOS is especially important in neural development, a process that requires high energy demand, though it also supports signaling processes, calcium homeostasis, and reactive oxygen species (ROS) production (reviewed in [78]). Proteins of the OXPHOS complex have been found dysregulated in four postmortem brain regions of patients [79]. Previous studies have pointed to differential gene expression [68] related to mitochondria function and morphological changes [4] in mitochondria of hiPSC-derived NPCs from schizophrenia patients, which complemented by protein dysregulation, suggest functional impairment in this organelle in the disorder. Lastly, dysregulation of energy metabolism in schizophrenia has been observed in neuroimaging and postmortem brain studies [79, 80], including reduced OXPHOS activity in the prefrontal cortex of schizophrenia patients evaluated by integrated proteomics, metabolomics, and genomics [81].

On the other hand, the metabolism regulation of young neurons in schizophrenia patients, as found in the PPI network, falls under cellular amino acid metabolic processes and cytosolic tRNA aminoacylation. Several studies have shown that differences in the level of amino acids are associated with schizophrenia patients (reviewed in [82]). And this is the first time we are able to point those dysregulations on hiPSC-derived neurons in schizophrenia studies, indicating new targets for future in vitro neuronal developmental studies using iPSC. Amino acids, such as glycine and serine, play an important role in the neurotransmission system. In the context of NMDA receptor activation, both glycine and D-serine have the same function; however, it has been hypothesized that in neurodevelopment, D-serine may have greater importance than glycine, since it is found in abundance during this period, as well as being co-located with NMDA receptors [50]. Defects in glycine binding to its receptor have been linked to cognitive and negative symptoms of schizophrenia [83], while a decrease in D-serine has been observed in different brain regions, the blood, and the cerebral spinal fluid of schizophrenia patients (reviewed in [84]). Knockout mice for serine racemase are also used as an established model of schizophrenia because they satisfactorily mimic NMDA receptor hypofunction [85].

The prefrontal cortex stands out in schizophrenia development when comparing brain organoids and regions of postmortem tissue

The dorsolateral prefrontal cortex (DLPFC) is located in the prefrontal cortex and is involved in working memory and attention [86], both of which are known to be altered in patients with schizophrenia [87,88,89,90]. There is an excess of synapse pruning in the DLPFC that can reach up to 30% more than what occurs in a healthy individual (reviewed in [88]). This loss of synaptic connections in patients with schizophrenia usually occurs in late adolescence or early adulthood, a period marking the end of neurodevelopment and the point at which the onset of the disease usually occurs. In addition, this increase in synaptic pruning is related to the previous appearance of prodromal symptoms; that is, this dysregulation begins before the noticeable onset of the disease and is linked to synaptic formation [88].

Gulsuner et al. [91] mapped de novo mutations in fetal DLPFC, taking into consideration those who were later diagnosed with schizophrenia; genes for proteins that play a role in axon guidance, synaptic transmission, and transcriptional regulation were all highlighted. These results are consistent with those obtained in this study, reinforcing the role of the synaptic system, especially in the prefrontal cortex, in brain development and the emergence of schizophrenia from a pathophysiological point of view.

The results presented here have shown that several essential pathways for neurodevelopment, such as axon guidance and synaptogenesis, have constituent proteins with reduced abundance in schizophrenia. Similar neurodevelopmental aspects are observed in the three patients, and despite their distinct genetic and symptomatic backgrounds, differential proteins were consistently observed. The strength of our study using neural-derived models of hiPSCs from patients with schizophrenia is a pronounced 89% of dysregulated proteins found to be downregulated, which we hypothesize to be a delay or slowing in brain development, thus reinforcing the hypothesis of schizophrenia as a neurodevelopmental disease. Our findings, as a developmental model of the disease, reinforce several of the aspects studied in schizophrenia, as shown by the similarities to the dysregulation found in postmortem brain studies, and could be further confirmed with a larger dataset, as to specific focus on the pathways presented here. Those similarities indicate that distinct genetic aspects converge to similar and common dysregulation in schizophrenia, which we can already observe in initial differentiation, leading to disruptions in neurodevelopment common to the disease. In addition to causing problems during neurodevelopment, these deregulations persist in the adult individual and thus contribute to the multifactorial nature of the pathophysiology of schizophrenia. Much of the data found in these analyses agree with previous studies using hiPSC cells from patients with schizophrenia, both confirming the data from this study and reinforcing the findings of the previous experiments. Nonetheless, new overlaps were also found in amino acid metabolism and spliceosomes, indicating new targets of study, at both the level of individual proteins as well as the level of biological pathways to which these proteins belong. Despite having a reduced number of subjects whose cells were hiPSCs sources, these data confirm the viability of postmortem tissue in proteomic studies, and further support patient-derived cell cultures as a valid technique to study schizophrenia.

Data availability

The proteomic datasets generated for this study can be found in the PRIDE proteomics data repository ( with the accession numbers PXD026381 and PXD026593.


  1. Schmidt MJ, Mirnics K. Neurodevelopment, GABA system dysfunction, and schizophrenia. Neuropsychopharmacology. 2015;40(1):190–206.

    Article  PubMed  Google Scholar 

  2. Gonzalez-Pinto A, Gutierrez M, Mosquera F, Ballesteros J, Lopez P, Ezcurra J, et al. First episode in bipolar disorder: misdiagnosis and psychotic symptoms. J Affect Disord. 1998;50(1):41–4.

    Article  CAS  PubMed  Google Scholar 

  3. Kohane IS, Masys DR, Altman RB. The incidentalome: a threat to genomic medicine. JAMA. 2006;296(2):212–5.

    Article  CAS  PubMed  Google Scholar 

  4. Brennand K, Savas JN, Kim Y, Tran N, Simone A, Hashimoto-Torii K, et al. Phenotypic differences in hiPSC NPCs derived from patients with schizophrenia. Mol Psychiatry. 2014;20(3):361–8.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Brennand KJ. Personalized medicine in a dish: the growing possibility of neuropsychiatric disease drug discovery tailored to patient genetic variants using stem cells. Stem Cell Investig. 2017;1:91–1.

    Article  Google Scholar 

  6. Pedrosa E, Sandler V, Shah A, Carroll R, Chang C, Rockowitz S, et al. Development of patient-specific neurons in schizophrenia using induced pluripotent stem cells. J Neurogen. 2011;25(3):88–103.

    Article  CAS  Google Scholar 

  7. Kelava I, Lancaster MA. Stem cell models of human brain development. Cell Stem Cell. 2016;18(6):736–48.

    Article  CAS  PubMed  Google Scholar 

  8. Lancaster MA, Knoblich JA. Generation of cerebral organoids from human pluripotent stem cells. Nat Protoc. 2014;9(10):2329–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Qian X, Nguyen HN, Song MM, Hadiono C, Ogden SC, Hammack C, et al. Brain-region-specific organoids using mini-bioreactors for modeling ZIKV exposure. Cell. 2016;165(5):1238–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Camp JG, Badsha F, Florio M, Kanton S, Gerber T, Wilsch-Bräuninger M, et al. Human cerebral organoids recapitulate gene expression programs of fetal neocortex development. Proc Natl Acad Sci USA. 2015;112(51):15672–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Nascimento JM, Saia-Cereda VM, Sartore RC, da Costa RM, Schitine CS, Freitas HR, et al. Human cerebral organoids and fetal brain tissue share proteomic similarities. Front Cell Dev Biol. 2019;7:303.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Brennand KJ, Simone A, Jou J, Gelboin-Burkhart C, Tran N, Sangar S, et al. Modelling schizophrenia using human induced pluripotent stem cells. Nature. 2011;473(7346):221–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Sochacki J, Devalle S, Reis M, de Moraes Maciel R, da Silveira Paulsen B, Brentani H, et al. Stem Cell Res. 2016;17(1):97–101.

    Article  CAS  PubMed  Google Scholar 

  14. Fraga AM, Sukoyan M, Rajan P, Braga D, Iaconelli A Jr, Franco JG, et al. Establishment of a Brazilian line of human embryonic stem cells in defined medium: implications for cell therapy in an ethnically diverse population. Cell Transpl. 2011;20(3):431–40.

    Article  Google Scholar 

  15. Baharvand H, Mehrjardi N-Z, Hatami M, Kiani S, Rao M, Haghighi M-M. Neural differentiation from human embryonic stem cells in a defined adherent culture condition. Int J Dev Biol. 2007;51(5):371–8.

  16. Dakic V, Nascimento JM, Sartore RC, de Moraes MR, Araujo DB, Ribeiro S, et al. Short term changes in the proteome of human cerebral organoids induced by 5-MeO-DMT. Sci Rep. 2017;7(1):12863.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Sartore RC, Cardoso SC, Lages YVM, Paraguassu JM, Stelling MP, da Costa RFM, et al. Trace elements during primordial plexiform network formation in human cerebral organoids. PeerJ. 2017;5: e2927.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Lancaster MA, Renner M, Martin C-A, Wenzel D, Bicknell LS, Hurles ME, et al. Cerebral organoids model human brain development and microcephaly. Nature. 2013;501(7467):373–9.

    Article  CAS  PubMed  Google Scholar 

  19. Dezonne RS, Sartore RC, Nascimento JM, Saia-Cereda VM, Romão LF, Alves-Leon SV, et al. Derivation of functional human astrocytes from cerebral organoids. Sci Rep. 2017;7:45091.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Cassoli JS, Brandao-Teles C, Santana AG, Souza GHMF, Martins-de-Souza D. Ion mobility-enhanced data-independent acquisitions enable a deep proteomic landscape of oligodendrocytes. Proteomics. 2017;17(21):1700209.

    Article  Google Scholar 

  21. Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13.

    Article  PubMed  Google Scholar 

  22. Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.

    Article  CAS  PubMed  Google Scholar 

  23. Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Comms. 2019;10(1):1523.

    Article  Google Scholar 

  24. Krzywinski M, Schein J, Birol İ, Connors J, Gascoyne R, Horsman D, et al. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19(9):1639–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Dakic V, de Moraes MR, Drummond H, Nascimento JM, Trindade P, Rehen SK. Harmine stimulates proliferation of human neural progenitors. PeerJ. 2016;4: e2727.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Suzuki M, Nelson AD, Eickstaedt JB, Wallace K, Wright LS, Svendsen CN. Glutamate enhances proliferation and neurogenesis in human neural progenitor cell cultures derived from the fetal cortex. Eur J Neurosci. 2006;24(3):645–53.

    Article  PubMed  Google Scholar 

  28. Wegner F, Kraft R, Busse K, Schaarschmidt G, Härtig W, Schwarz SC, et al. Glutamate receptor properties of human mesencephalic neural progenitor cells: NMDA enhances dopaminergic neurogenesis in vitro. J Neurochem. 2009;111(1):204–16.

    Article  CAS  PubMed  Google Scholar 

  29. Wang Z, Li P, Wu T, Zhu S, Deng L, Cui G. Axon guidance pathway genes are associated with schizophrenia risk. Exp Ther Med. 2018; 16: 4519-4526.

  30. Quadrato G, Brown J, Arlotta P. The promises and challenges of human brain organoids as models of neuropsychiatric disease. Nat Med. 2016, 22, 1220–1228.

  31. Paşca SP. The rise of three-dimensional human brain cultures. Nature. 2018;553(7689):437–45.

    Article  PubMed  Google Scholar 

  32. Fatemi SH, Folsom TD. The neurodevelopmental hypothesis of schizophrenia. Revisited Schizophrenia Bull. 2009;35(3):528–48.

    Article  Google Scholar 

  33. Owen MJ, O’Donovan MC, Thapar A, Craddock N. Neurodevelopmental hypothesis of schizophrenia. Br J Psychiatry. 2011;198(3):173–5.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Abdolmaleky HM, Smith CL, Faraone SV, Shafa R, Stone W, Glatt SJ, et al. Methylomics in psychiatry: modulation of gene-environment interactions may be through DNA methylation. Am J Med Genet B Neuropsychiatr Genet. 2004;127B(1):51–9.

    Article  PubMed  Google Scholar 

  35. Ellman LM, Deicken RF, Vinogradov S, Kremen WS, Poole JH, Kern DM, et al. Structural brain alterations in schizophrenia following fetal exposure to the inflammatory cytokine interleukin-8. Schizophr Res; 2010;121(1–3):46–54.

  36. Cheslack-Postava K, Brown AS. Prenatal infection and schizophrenia: A decade of further progress. Schizophrenia Res; 2022; 247:7–15.

  37. Singh SM, Murphy B, O’Reilly RL. Involvement of gene-diet/drug interaction in DNA methylation and its contribution to complex diseases: from cancer to schizophrenia. Clin Genet. 2003;64(6):451–60.

    Article  CAS  PubMed  Google Scholar 

  38. Sørensen HJ, Mortensen EL, Schiffman J, Reinisch JM, Maeda J, Mednick SA. Early developmental milestones and risk of schizophrenia: a 45-year follow-up of the Copenhagen Perinatal Cohort. Schizophr Res. 2010;118(1–3):41–7.

    Article  PubMed  PubMed Central  Google Scholar 

  39. St Clair D, Xu M, Wang P, Yu Y, Fang Y, Zhang F, et al. Rates of adult schizophrenia following prenatal exposure to the Chinese famine of 1959–1961. JAMA. 2005;294(5):557–62.

    Article  CAS  PubMed  Google Scholar 

  40. Brown AS, Derkits EJ. Prenatal infection and schizophrenia: a review of epidemiologic and translational studies. Am J Psychiatry. 2010;167(3):261–80.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Henkemeyer M, Itkis OS, Ngo M, Hickmott PW, Ethell IM. Multiple EphB receptor tyrosine kinases shape dendritic spines in the hippocampus. J Cell Biol. 2003;163(6):1313–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Klein R. Bidirectional modulation of synaptic functions by Eph/ephrin signaling. Nat Neurosci. 2009;12(1):15–20.

    Article  CAS  PubMed  Google Scholar 

  43. Takeuchi S, Katoh H, Negishi M. Eph/ephrin reverse signalling induces axonal retraction through RhoA/ROCK pathway. J Biochem. 2015;158(3):245–52.

    Article  CAS  PubMed  Google Scholar 

  44. Lai K-O, Ip NY. Synapse development and plasticity: roles of ephrin/Eph receptor signaling. Curr Opin Neurobiol; 2009;19(3):275–83.

  45. Kao T-J, Kania A. Ephrin-mediated cis-attenuation of Eph receptor signaling is essential for spinal motor axon guidance. Neuron. 2011;71(1):76–91.

    Article  CAS  PubMed  Google Scholar 

  46. Saia-Cereda VM, Cassoli JS, Martins-de-Souza D, Nascimento JM. Psychiatric disorders biochemical pathways unraveled by human brain proteomics. Eur Arch Psychiatry Clin Neurosci. 2016;1:1–15.

    Google Scholar 

  47. Kayser MS, McClelland AC, Hughes EG, Dalva MB. Intracellular and trans-synaptic regulation of glutamatergic synaptogenesis by EphB receptors. J Neurosci. 2006;26(47):12152–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Moyer CE, Shelton MA, Sweet RA. Dendritic spine alterations in schizophrenia. Neurosci Lett. 2015;601:46–53.

    Article  CAS  PubMed  Google Scholar 

  49. Javitt DC, Zukin SR. Recent advances in the phencyclidine model of schizophrenia. Am J Psychiatry. 1991;148(10):1301–8.

    Article  CAS  PubMed  Google Scholar 

  50. van Horn MR, Sild M, Ruthazer ES. D-serine as a gliotransmitter and its roles in brain development and disease. Front Cell Neurosci. 2013;1:7–39.

    Google Scholar 

  51. Trubetskoy V, Pardiñas AF, Qi T, Panagiotaropoulou G, Awasthi S, Bigdeli TB, et al. Mapping genomic loci implicates genes and synaptic biology in schizophrenia. Nature. 2022;604(7906):502–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Cornell B, Toyo-Oka K. 14-3-3 proteins in brain development: neurogenesis, neuronal migration and neuromorphogenesis. Front Mol Neurosci. 2017;10:318.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Zhang J, Zhou Y. 14–3–3 Proteins in Glutamatergic Synapses. Neural Plast; 2018;2018:8407609.

  54. Föcking M, Dicker P, English JA, Schubert KO, Dunn MJ, Cotter DR. Common proteomic changes in the hippocampus in schizophrenia and bipolar disorder and particular evidence for involvement of cornu ammonis regions 2 and 3. Arch Gen Psychiatry. 2011;68(5):477–88.

    Article  PubMed  Google Scholar 

  55. Xu Y, Yue W, Yao-Shugart Y, Li S, Cai L, Li Q, et al. Exploring transcription factors-microRNAs co-regulation networks in schizophrenia. Schizophr Bull; 2016;42(4):1037–45.

  56. Toyooka K, Asama K, Watanabe Y, Muratake T, Takahashi M, Someya T, et al. Decreased levels of brain-derived neurotrophic factor in serum of chronic schizophrenic patients. Psychiatry Res. 2002;110(3):249–57.

    Article  CAS  PubMed  Google Scholar 

  57. Stark KL, Xu B, Bagchi A, Lai W-S, Liu H, Hsu R, et al. Altered brain microRNA biogenesis contributes to phenotypic deficits in a 22q11-deletion mouse model. Nat Genet. 2008;40(6):751–60.

    Article  CAS  PubMed  Google Scholar 

  58. Clapcote SJ, Lipina TV, Millar JK, Mackie S, Christie S, Ogawa F, et al. Behavioral phenotypes of Disc1 missense mutations in mice. Neuron. 2007;54(3):387–402.

    Article  CAS  PubMed  Google Scholar 

  59. Vawter MP, Barrett T, Cheadle C, Sokolov BP, et al. Application of cDNA microarrays to examine gene expression differences in schizophrenia. Brain Res Bull; 2001; 55(5): 641–50.

  60. Bell R, Munro J, Russ C, Powell JF, Bruinvels A, Kerwin RW, et al. Systematic screening of the 14–3–3 eta (eta) chain gene for polymorphic variants and case-control analysis in schizophrenia. Am J Med Genet; 2000;96(6):736–43.

  61. Muratake T, Hayashi S, Ichikawa T, Kumanishi T, Ichimura Y, Kuwano R, et al. Structural organization and chromosomal assignment of the human 14-3-3 eta chain gene (YWHAH). Genomics. 1996;36(1):63–9.

    Article  CAS  PubMed  Google Scholar 

  62. Saia-Cereda VM, Cassoli JS, Schmitt A, Falkai P, Nascimento JM, Martins-de-Souza D. Proteomics of the corpus callosum unravel pivotal players in the dysfunction of cell signaling, structure, and myelination in schizophrenia brains. Eur Arch Psychiatry Clin Neurosci. 2015;265(7):601–12.

    Article  PubMed  Google Scholar 

  63. Martins-de-Souza D, Gattaz WF, Schmitt A, Maccarrone G, Hunyadi-Gulyás E, Eberlin MN, et al. Proteomic analysis of dorsolateral prefrontal cortex indicates the involvement of cytoskeleton, oligodendrocyte, energy metabolism and new potential markers in schizophrenia. J Psychiatr Res. 2009;43(11):978–86.

    Article  PubMed  Google Scholar 

  64. Middleton FA, Peng L, Lewis DA, Levitt P, Mirnics K. Altered expression of 14–3–3 genes in the prefrontal cortex of subjects with schizophrenia. Neuropsychopharmacology. 2005;30(5):974–83.

    Article  CAS  PubMed  Google Scholar 

  65. Konietzny A, Bär J, Mikhaylova M. Dendritic actin cytoskeleton: structure, functions, and regulations. Front Cell Neurosci; 2017;11:147.

  66. Kuhn TB, Meberg PJ, Brown MD, Bernstein BW, Minamide LS, Jensen JR, et al. Regulating actin dynamics in neuronal growth cones by ADF/cofilin and rho family GTPases. J Neurobiol. 2000;44(2):126–44.

    Article  CAS  PubMed  Google Scholar 

  67. Irie F, Yamaguchi Y. EphB receptors regulate dendritic spine development via intersectin, Cdc42 and N-WASP. Nat Neurosci. 2002;5(11):1117–8.

    Article  CAS  PubMed  Google Scholar 

  68. Kathuria A, Lopez-Lengowski K, Jagtap SS, McPhie D, Perlis RH, Cohen BM, et al. Transcriptomic landscape and functional characterization of induced pluripotent stem cell-derived cerebral organoids in schizophrenia. JAMA Psychiat. 2020;77(7):745–54.

    Article  Google Scholar 

  69. Ward AJ, Cooper TA. The pathobiology of splicing. J Pathol. 2010;220(2):152–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Huntsman MM, Tran BV, Potkin SG, Bunney WEJ, Jones EG. Altered ratios of alternatively spliced long and short gamma2 subunit mRNAs of the gamma-amino butyrate type A receptor in prefrontal cortex of schizophrenics. Proc Natl Acad Sci USA. 1998;95(25):15066–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Park E, Iaccarino C, Lee J, Kwon I, Baik SM, Kim M, et al. Regulatory roles of heterogeneous nuclear ribonucleoprotein M and Nova-1 protein in alternative splicing of dopamine D2 receptor pre-mRNA. J Biol Chem; 2011;286(28):25301–8.

  72. Morikawa T, Manabe T. Aberrant regulation of alternative pre-mRNA splicing in schizophrenia. Neurochem Int. 2010;57(7):691–704.

    Article  CAS  PubMed  Google Scholar 

  73. Cassoli JS, Iwata K, Steiner J, Guest PC, Turck CW, Nascimento JM, et al. Effect of MK-801 and clozapine on the proteome of cultured human oligodendrocytes. Front Cell Neurosci. 2016;10:1.

    Article  Google Scholar 

  74. Saia-Cereda VM, Santana AG, Schmitt A, Falkai P, Martins-de-Souza D. The Nuclear proteome of white and gray matter from schizophrenia postmortem brains. Mol Neuropsychiatry; 2017;3(1):37–52.

  75. Zhao L, Mandler MD, Yi H, Feng Y. Quaking I controls a unique cytoplasmic pathway that regulates alternative splicing of myelin-associated glycoprotein. Proc Natl Acad Sci USA. 2010;107(44):19061–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Iwata K, Matsuzaki H, Manabe T, Mori N. Altering the expression balance of hnRNP C1 and C2 changes the expression of myelination-related genes. Psychiatry Res; 2011 Dec;190(2–3):364–6.

  77. Falk MJ. Neurodevelopmental manifestations of mitochondrial disease. J Dev Behav Pediatr; 2010;31(7):610–21.

  78. Bergman O, Ben-Shachar D. Mitochondrial Oxidative Phosphorylation System (OXPHOS) Deficits in Schizophrenia: Possible Interactions with Cellular Processes. Canadian journal of psychiatry Revue canadienne de psychiatrie. Can J Psychiatry; 2016;61(8):457–69.

  79. Maurer I, Zierz S, Möller H. Evidence for a mitochondrial oxidative phosphorylation defect in brains from patients with schizophrenia. Schizophr Res; 2001;48(1):125–36.

  80. Rice MW, Smith KL, Roberts RC, Perez-Costas E, Melendez-Ferro M. Assessment of cytochrome C oxidase dysfunction in the substantia nigra/ventral tegmental area in schizophrenia. PLoS ONE. 2014;9(6): e100054.

    Article  PubMed  PubMed Central  Google Scholar 

  81. Prabakaran S, Swatton JE, Ryan MM, Huffaker SJ, Huang JT-J, Griffin JL, et al. Mitochondrial dysfunction in schizophrenia: evidence for compromised brain metabolism and oxidative stress. Mol Psychiatry. 2004;9(7):684–97–643.

  82. Saleem S, Shaukat F, Gul A, Arooj M, Malik A. Potential role of amino acids in pathogenesis of schizophrenia.. Int J Health Sci (Qassim); 2017;11(3):63–8.

  83. Labrie V, Lipina T, Roder JC. Mice with reduced NMDA receptor glycine affinity model some of the negative and cognitive symptoms of schizophrenia. Psychopharmacology. 2008;200(2):217–30.

    Article  CAS  PubMed  Google Scholar 

  84. MacKay M-AB, Kravtsenyuk M, Thomas R, Mitchell ND, Dursun SM, Baker GB. D-Serine: potential therapeutic agent and/or biomarker in schizophrenia and depression?. Front Psychiatry; 2019;10:25.

  85. Balu DT, Li Y, Puhl MD, Benneyworth MA, Basu AC, Takagi S, et al. Multiple risk pathways for schizophrenia converge in serine racemase knockout mice, a mouse model of NMDA receptor hypofunction. Proc Natl Acad Sci USA. 2013;110(26):E2400–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  86. Kane MJ, Engle RW. The role of prefrontal cortex in working-memory capacity, executive attention, and general fluid intelligence: an individual-differences perspective. Psychon Bull Rev; 2002;9(4):637–71.

  87. Gold JM, Robinson B, Leonard CJ, Hahn B, Chen S, McMahon RP, et al. Selective attention, working memory, and executive function as potential independent sources of cognitive dysfunction in schizophrenia. Schizophr Bull; 2018;44(6):1227–34.

  88. Bennett MR. Synapse formation and regression in the cortex during adolescence and in schizophrenia. Med J Aust; 2009;190(S4):S14–6.

  89. Wible CG, Anderson J, Shenton ME, Kricun A, Hirayasu Y, Tanaka S, et al. Prefrontal cortex, negative symptoms, and schizophrenia: an MRI study. Psychiatry Res; 2001;108(2):65–78.

  90. Abbruzzese M, Bellodi L, Ferri S, Scarone S. Frontal lobe dysfunction in schizophrenia and obsessive-compulsive disorder: a neuropsychological study. Brain Cogn; 1995;27(2):202–12.

  91. Gulsuner S, Walsh T, Watts AC, Lee MK, Thornton AM, Casadei S, et al. Spatial and temporal mapping of de novo mutations in schizophrenia to a fetal prefrontal cortical network. Cell. 2013;154(3):518–29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references


The authors thank the excellent technical support of Ludmila SS Bastos, Ismael Carlos da Silva Gomes, Marcelo do Nascimento Costa, Gabriela Lopes Vitória, Jarek Sochacki, and Paulo Baldasso.


This work was supported by the São Paulo Research Foundation (FAPESP), grants 2014/21035-0, 2016/07332-7, 2017/25588-1, 2018/14666-4, 2018/01410-1, 2019/00098-7, 2019/05747-3 and 2019/25957-2. Other funds were provided by the Fund for Support to Teaching, Research and Outreach Activities (FAEPEX), National Council for Scientific and Technological Development (CNPq), and the Instituto Nacional de Ciência e Tecnologia de Neurociência Translacional (INCT-INNT), Foundation for Research Support in the State of Rio de Janeiro (FAPERJ), in addition to intramural grants from D’Or Institute for Research and Education.

Author information

Authors and Affiliations



JMN, DMS, and SKR conceptualized and designed the study; JMN performed iPSC-derived cell experiments. JMN, VMS-C, and GSZ performed experiments, collected, assembled, and analyzed data; GRO analyzed data; JMN and VMSC wrote the paper, prepared figures, and revised the final manuscript. VCC performed and analyzed MRM experiments. BJS reviewed and contributed to the discussion. SKR and DMS provided financial support and coordinated the study. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Juliana Minardi Nascimento, Stevens K. Rehen or Daniel Martins-de-Souza.

Ethics declarations

Ethics approval and consent to participate

This study was performed following international standards and with the approval of the D’Or Institute for Research and Education ethics council (CAAE: 32385314.9.0000.5249).

Consent for publication

Not applicable.

Competing interests

The authors have declared no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

Cell lines information

Additional file 2: Table S2.

Protein expression of iPSC-derived NPCs, neurons (21 DIV), and cerebral organoids from schizophrenia patients.

Additional file 3: Table S3.

Compilation of proteomic studies in different regions of the postmortem brain.

Additional file 4.

Additional tables and figures.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Nascimento, J.M., Saia-Cereda, V.M., Zuccoli, G.S. et al. Proteomic signatures of schizophrenia-sourced iPSC-derived neural cells and brain organoids are similar to patients' postmortem brains. Cell Biosci 12, 189 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Neuroproteomics
  • Schizophrenia
  • Brain organoids
  • Neural cells
  • Human pluripotent stem cells